Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “lateral”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 289 records · Page 16

Time‐Lapse Image Classification Using a Diffractive Neural Network

Diffractive deep neural networks (D 2 NNs), comprised of spatially engineered passive surfaces, collectively process optical input information at the speed of light propagation through a thin diffractive volume, without any external computing power. Diffractive networks were demonstrated to achieve all‐optical object classification and perform universal linear transformations. Herein, a “time‐lapse” image classification scheme using a diffractive network is demonstrated for the first time, significantly advancing its classification accuracy and generalization performance on complex input objects by using the lateral movements of the input objects and/or the diffractive network, relative to each other. In a different context, such relative movements of the objects and/or the camera are routinely being used for image super‐resolution applications; inspired by their success, a time‐lapse diffractive network is designed to benefit from the complementary information content created by controlled or random lateral shifts. The design space and performance limits of time‐lapse diffractive networks are numerically explored, revealing a blind testing accuracy of 62.03% on the optical classification of objects from the CIFAR‐10 dataset. This constitutes the highest inference accuracy achieved so far using a single diffractive network on the CIFAR‐10 dataset. Time‐lapse diffractive networks will be broadly useful for the spatiotemporal analysis of input signals using all‐optical processors.

36 MATERIALS SCIENCE↗

Riverscape heterogeneity in estimated Chinook Salmon emergence phenology and implications for size and growth

Abstract Many salmonid‐bearing rivers exhibit thermal and hydrologic heterogeneity at multiple spatial and temporal scales, but how this translates into spatiotemporal patterns of fry emergence is poorly understood. Understanding this variability is important because emergence timing determines the biophysical conditions fish first experience (e.g., temperature, flow, and food supply), thereby influencing growth opportunities and survival during this critical life stage. We predicted spring Chinook Salmon ( Oncorhynchus tshawytscha ) emergence phenology across four northeastern Oregon subbasins over 5–9 years using empirical spawning and temperature data. We then related interannual emergence timing estimates to juvenile salmon size and growth rates at consistent sampling locations. There were clear longitudinal patterns of predicted emergence timing in each subbasin: The shape of these patterns was consistent among years, but not among subbasins. In two subbasins, emergence occurred progressively later with distance upstream, whereas in the other two subbasins emergence was earliest at upstream sites. Within each year, median emergence dates among sites within each subbasin ranged between 44 and 58 days. This spatial variation was comparable to interannual variation, with median emergence dates for a given location in each subbasin ranging between 47 and 74 days among years. Contrary to our expectations, juvenile salmon were not larger in years with earlier emergence, owing to slower estimated spring and summer growth rates compared to years with later emergence. Despite large interannual variation in estimated emergence dates, these results suggest that other factors (e.g., stream flow, temperature, and density‐dependence) were more important than growth duration in determining juvenile salmon growth rates and size among years. We demonstrated considerable spatial and interannual variation in emergence phenology within these subbasins. Understanding how this variation translates to spatiotemporal patterns of juvenile salmon habitat use, growth, and survival has important implications for guiding restoration efforts and understanding how climate change may impact these populations.

Kaylor, Matthew J.↗

Earthquake-induced impact of base-isolated buildings: theory, numerical modeling, and design solutions

Earthquake shaking more intense than that used to size the horizontal clearance between a base–isolated building and near–rigid perimeter moat wall will result in hard impact, producing high–frequency, high–amplitude acceleration response in the structure and supported equipment. Herein this paper provides a design solution for the damaging effects of hard impact by installing a compliant engineered element in the load path between the base–isolated building and the moat wall, resulting in soft impact and a much smaller acceleration response. The engineered element assumed herein is a commercial–off–the–shelf marine fender with mechanical properties determined by physical testing. The attachment of a flexible engineered element, with well–defined stiffness and damping, to a near–rigid moat wall, simplifies the numerical modeling of the building–moat wall system and eliminates the need to bound the lateral stiffness of the wall for impact calculations. The simple model of the engineered element can be implemented in commercial finite element codes. Theory is developed for two–sided impact of a single–degree–of–freedom oscillator. Analytical solutions are derived for the shifted first–mode frequency of the impacted oscillator and for its free–vibration response. The shifted first–mode frequency is a function of the composite lateral stiffness of the isolator–engineered element assembly and its earthquake–induced displacement. Local peaks in the spectral response of the impacted oscillator form at odd integer multiples of the shifted first–mode frequency. The analytical solutions can be used to verify, in part, the numerical model used for impact analysis.

42 ENGINEERING↗

River–aquifer interactions enhancing evapotranspiration in a semiarid riparian zone: A modelling study

The hydrologic flows across the river–aquifer interface play an important role in groundwater dynamics and biogeochemical reactions within the subsurface; however, little is known about the effects of river–aquifer interactions on land surface processes. In this study, we developed a fully coupled three-dimensional (3D) land surface and subsurface model at a high resolution (~1 km) that accounts for high-frequency hydrologic exchange flow conditions to investigate how river–aquifer interactions modulate surface water budgets in the Upper Columbia-Priest Rapids watershed, a typical semiarid watershed located in the northwestern United States where river stage fluctuates in response to reservoir releases changing. Our results show that the spatiotemporal dynamics of river–aquifer interactions are highly heterogeneous, driven mainly by river-stage fluctuations. Adding 6.64 × 10 6 m 3 year –1 of water over the watershed from the river to groundwater owing to the lateral flow, river–aquifer interactions led to an increase in soil evaporation and transpiration supplied by higher soil moisture content, particularly in deeper subsurface. In a hypothetic future scenarios where a 5-m rise in river stage was assumed, the hydrologic flow exchange rates were intensified, resulting in higher surface water over the entire watershed. Overall, lateral flow induced by river–aquifer exchanges leads to an increase in evapotranspiration of ~75% in the historical period and of ~83% in the hypothetical future scenario. Finally, our study demonstrates the potential of coupled model as an effective tool for understanding river–aquifer–land surface interactions, and indicates that river–aquifer interactions fundamentally alter the water balance of the riparian zone for the semiarid watershed and will likely become more frequent and intense in the future under the effects of climate change.

54 ENVIRONMENTAL SCIENCES↗

An algorithm for detecting and quantifying disturbance and recovery in high-frequency time series

Determining when a disturbance has occurred, its severity, and when the system recovered, is important to numerous questions in the aquatic sciences. This problem can be conceptualized as the timing and degree of perturbation from a typical state, and when the system returns to that typical state. We present an algorithm for detecting disturbance and recovery designed for high-frequency time series, e.g., data produced by automated sampling devices in instrumented buoys and flux towers. The algorithm quantifies differences in the empirical cumulative distribution functions of moving windows over reference and evaluation periods, and is sensitive to changes in the mean, variance, and higher statistical moments. Tests on simulated data show it accurately identifies disturbance and recovery. Three case studies illustrate the application of our algorithm in different empirical settings. A case study on dissolved oxygen in a Florida, USA estuary following a hurricane identified the disturbance and recovery 73 d later. A case study on air temperature and net ecosystem exchange in the Florida everglades identified cold snaps coinciding with periods of reduced carbon uptake. A case study on rotifer abundance following zebra mussel invasion in the Hudson River, NY showed rotifer collapse following invasion and recovery over a decade later. Methods such as ours can improve understanding response to disturbance and facilitate comparative and synthetic study of disturbance impacts across ecosystems.

54 ENVIRONMENTAL SCIENCES↗

On the Equilibrium Electrostatic Potential and Light-Induced Charge Redistribution in Halide Perovskite Structures

Lead halide perovskites are semiconductor materials which are employed as nonintentionally doped absorbers inserted between two selective carrier transport layers (SCTL), realizing a p-i-n or n-i-p heterojunction. In our study, we have developed and investigated a lateral device, based on methylammonium lead iodide (MAPbI3) in which the p-i-n heterojunction develops in the horizontal direction. Our research suggests that the effective doping level in the MAPbI3 film should be very low, below 1012 cm-3. Along the vertical direction, this doping level is not enough to screen the electric field of the buried heterojunction with the SCTL. The perovskite work function is therefore affected by the work function of the SCTL underneath. From drift-diffusion simulations, we show that intrinsic perovskite-SCTL structures develop mV range surface photovoltages (SPVs) under continuous illumination. However, perovskite-SCTL structures can develop SPVs of hundreds of mV, as confirmed by our measurements. We therefore analyzed the compatibility between low doping and low defect densities in the perovskite layer and such high SPV values using numerical modeling. It is shown that these high SPV values could originate from electronic processes due to large band offsets in the buried perovskite-SCTL heterojunctions, or at the SCTL-transparent conductive oxide (TCO) buried heterojunction. However, such electronic processes can hardly explain the long SPV persistence after switching off the illumination.

defects↗

Hexagonal Boron Nitride: Physical Properties, Hydride Vapor‐Phase Epitaxy Growth of Large‐Diameter Quasi‐Bulk Wafers and Applications

Hexagonal boron nitride ( h ‐BN), with its ultrawide bandgap and 2D structure, holds an immense promise for advanced semiconductor applications. Scaling bulk crystals to large‐diameter wafers, crucial for complex device fabrication, remains a challenge with high temperature, high pressure, and metal flux solution methods. To address this, recent efforts have focused on hydride vapor‐phase epitaxy (HVPE) for producing large diameter thick h ‐BN quasi‐bulk wafers, reaching hundreds of micrometers. These HVPE‐grown quasi‐bulk crystals exhibit excellent c ‐axis long‐range order. Notably, the in‐plane mobility‐lifetime products for both electrons and holes surpass 10 −4 cm 2 V −1 , which are two orders of magnitude greater than the out‐of‐plane (vertical) values, highlighting the potential for high‐performance devices leveraging superior lateral transport. Lateral detectors fabricated from 100 μm thick B‐10 enriched h‐ BN wafers have achieved a record 60% thermal neutron detection efficiency. Based on its physical properties, h‐ BN appears to be an outstanding material of choice for light‐triggered electronic power switches capable of supporting high ‐ voltage and high ‐ power operations. These recent advancements in large‐diameter h ‐BN quasi‐bulk crystal growth, enabled by HVPE, pave the way for applications spanning deep UV photonics, high‐power electronics, high‐efficiency neutron detection, and quantum information technologies, establishing h ‐BN as both a versatile active semiconductor and an ideal substrate.

Jiang, Hongxing [Department of Electrical and Comp↗

Understanding the Effects of Inhomogeneities at the Back Interface of CdTe‐Based Solar Cells Using 2D Modeling

One-dimensional modeling cannot capture lateral inhomogeneities in CdTe-based devices. Here, we use 2D modeling to investigate the role of varying energetics at the back interface. We consider improvements in the back interface layer (BIL) through either reducing back surface recombination velocity (BSRV) or decreasing the downward band bending near the back interface. We show that when the BSRV is reduced, but strong downward band bending remains, there is no change in the device performance until the BSRV of 90% of the back interface is improved by the BIL. On the other hand, any coverage with a BIL that improves band bending results in device improvements. We use band bending, back interface recombination current densities, and voltage dependent current flow through the device to understand these improvements. The modeling shows that lateral flow of carriers greatly affects device performance, which is not captured in parallel diode modeling, and demonstrates improved understanding with 2D modeling.

2D modelling↗

FAST.Farm Development and Validation of Structural Load Prediction Against Large Eddy Simulations

FAST.Farm is a mid-fidelity engineering tool developed by the National Renewable Energy Laboratory targeted at accurately and efficiently predicting wind turbine power production and structural loading in wind farm settings, including wake interactions between turbines. FAST.Farm is based on several principles of the dynamic wake meandering (DWM) model, but also addresses limitations of previous DWM implementations. Previous FAST.Farm studies have shown the similarities and differences between FAST.Farm and large eddy simulations for rigid turbine cases. The objective of this work is to quantify the ability of FAST.Farm to accurately predict turbine structural response in a small wind farm. This is done by comparing FAST.Farm structural response to SOWFA-OpenFAST results for three laterally-aligned turbines. The purpose of this study is to characterize the similarities and differences between FAST.Farm and a higher fidelity model for predicting turbine structural response in a wind farm for differing atmospheric inflows. Strong statistical agreement was found between FAST.Farm and SOWFA-OpenFAST structural response for the non-waked upstream turbine, and good agreement was found for the downstream turbines for most structural quantities. Higher differences were seen for downstream turbines with low ambient turbulence intensity or yawed turbines, suggesting areas for FAST.Farm wake dynamics modeling improvements. For all cases and turbines, small statistical differences were seen between blade deflections and bending moments, with larger differences for tower-top and tower-base bending moments. Overall, the results establish confidence for applying FAST.Farm to wind farm power and loads analyses and identify areas where further model validations and model improvements should be targeted.

49 EE - Wind and Water Power Program - Wind (EE-4W↗

Blade planform design optimization to enhance turbine wake control

Abstract This study considers optimizing the planform of wind turbine blades to ultimately enhance wind plant controls, namely, wake steering strategies. Adjoint‐enabled unsteady actuator line simulations are carried out to obtain gradients for optimization of several different performance objectives with respect to blade chord length at 10 locations along blade span. We demonstrate different blade design optimizations that can maximize time‐averaged lateral wake deflection, entrainment of kinetic energy, or total power of multiple turbines. Our optimized designs can produce a 4+% increase in wake deflection, a 4× increase in vertical kinetic energy entrainment, or a 3.6% increase in power when compared with the baseline case. While lateral wake deflection is only modestly sensitive to chord changes, we find that increasing the outboard chord length can dramatically increase kinetic energy entrainment, resulting in faster wake recovery and gains in net power. While this work develops only a few case studies emphasizing relative performance improvements and general trends, these results show the promise of a framework that combines mid‐fidelity computation with adjoint‐based optimization for control and design problems.

WIND ENERGY↗

Synchrotron-based diffraction-enhanced imaging and diffraction-enhanced imaging combined with CT X-ray imaging systems to image seeds at 30 keV

Utilized the upgraded Synchrotron-based non-destructive Diffraction-enhanced imaging and Diffraction-enhanced imaging coupled with CT X-ray imaging systems to image the chickpea seeds, to enhance the contrast in plant root architecture, visibility of fine structures of root architecture growth and some aspects of physiology at acceptable level. DEI-CT images were acquired with 30 keV synchrotron X-rays. A series of DEI-CT slices were assembled together, to form a 3D data set. DEI-CT images explored more structural information and morphology. Noticed detailed anatomical, physiological observations, and contrast mechanisms. Furthermore, with these systems, some of the complex plant traits, root morphology, growth of laterals and subsequent laterals can be visualized directly.

36 MATERIALS SCIENCE↗

Breeding latitude is associated with the timing of nesting and migration around the annual calendar among Purple Martin (Progne subis) populations

Timing is critical for species that migrate seasonally, as timing mismatches can lead to reduced reproductive success and population declines. We investigated the effect of breeding latitude on the migration timing of a long-distance migratory songbird (Purple Martin, Progne subis subis) by using light-level geolocators deployed across breeding sites between 26° and 53° latitude (n = 329). We used nesting data (1st egg dates; n=28,165) to explore latitudinal patterns in the timing of breeding. Using generalized linear mixed models (GLMMs) we found that breeding latitude was highly associated with timing around the calendar, explaining population-level timing across both spring and fall migration. Populations breeding at more northern latitudes migrated later in fall (mean days/degree latitude =2), in spring (mean days/degree latitude =3), and nested later (days/degree latitude = 3) than populations breeding at more southern latitudes. Generalized additive mixed models (GAMMs) revealed significant non-linearity between breeding latitude and spring and fall migration timing, and the timing of nesting. In spring, the most northern breeders (> 40 degrees N) were more synchronous in their timing, departing their nonbreeding sites and arriving earlier at their breeding sites relative to their breeding latitude. These most northern breeders also had similar nest timing and fall migration departure timing. Breeding latitude was also associated with seasonal partitioning of migratory and sedentary periods. The most northern breeders spent 10% more time on fall migration, 2% more at wintering grounds, 1% more on spring migration, and 13% less time at their breeding sites. Our results support the notion that selection for timing at breeding latitudes can have year-round effects on the timing of seasonal events in long-distance migratory songbirds. Future studies should examine whether these timing patterns associated with latitude limit phenotypic responses to climate change.

54 ENVIRONMENTAL SCIENCES↗

Effects of volitional emigration timing and smolt size on survival and age-at-return in a Pacific salmon hatchery population

Abstract Many Pacific salmon populations are returning from sea at younger ages and smaller sizes. Hatchery culture, management practices, and environmental factors influence juvenile release size and emigration timing, which in turn affect important demographic characteristics in returning adults. We analyzed data from approximately 345,000 tagged spring Chinook Salmon juveniles exiting Cle Elum Supplementation and Research Facility (Yakima River, Washington, USA) acclimation sites over thirteen brood years (2002–2014), evaluating smolt size, emigration timing, river flow, juvenile survival, and age-at-return. We observed a relationship between size and volitional exit timing of smolts from acclimation sites, with larger fish tending to emigrate earlier than smaller fish. Early emigration was also coincident with lower river flows near acclimation sites. Later emigration timing was correlated with an increase in apparent survival of juveniles to Bonneville Dam (500–530 km downstream of acclimation sites), but also with a lower rate of survival to return from sea. In general, for juveniles successfully emigrating downstream of Bonneville Dam, age-at-return increased with decreasing juvenile fish size and later emigration timing. Our results support a growing body of evidence that hatchery practices may result in larger smolts that tend to return at younger ages. Given the biological and economic consequences of younger age-at-maturation, methods to reverse this trend should be further explored and implemented.

Bosch, William J. (ORCID:0000000300447390)↗

Elucidating the Onset of Plasticity in Sliding Contacts Using Differential Computational Orientation Tomography

Abstract Depending on the mechanical and thermal energy introduced to a dry sliding interface, the near-surface regions of the mated bodies may undergo plastic deformation. In this work, we use large-scale molecular dynamics simulations to generate “differential computational orientation tomographs” (dCOT) and thus highlight changes to the microstructure near tribological FCC alloy surfaces, allowing us to detect subtle differences in lattice orientation and small distances in grain boundary migration. The analysis approach compares computationally generated orientation tomographs with their undeformed counterparts via a simple image analysis filter. We use our visualization method to discuss the acting microstructural mechanisms in a load- and time-resolved fashion, focusing on sliding conditions that lead to twinning, partial lattice rotation, and grain boundary-dominated processes. Extracting and laterally averaging the color saturation value of the generated tomographs allows us to produce quantitative time- and depth-resolved maps that give a good overview of the progress and severity of near-surface deformation. Corresponding maps of the lateral standard deviation in the color saturation show evidence of homogenization processes occurring in the tribologically loaded microstructure, frequently leading to the formation of a well-defined separation between deformed and undeformed regions. When integrated into a computational materials engineering framework, our approach could help optimize material design for tribological and other deformation problems. Graphic Abstract .

36 MATERIALS SCIENCE↗

The Surface-Topography Challenge: A Multi-Laboratory Benchmark Study to Advance the Characterization of Topography

Surface performance is critically influenced by topography in virtually all real-world applications. The current standard practice is to describe topography using one of a few industry-standard parameters. The most commonly reported number is Ra, the average absolute deviation of the height from the mean line (at some, not necessarily known or specified, lateral length scale). However, other parameters, particularly those that are scale-dependent, influence surface and interfacial properties; for example the local surface slope is critical for visual appearance, friction, and wear. The present Surface-Topography Challenge was launched to raise awareness for the need of a multi-scale description, but also to assess the reliability of different metrology techniques. In the resulting international collaborative effort, 153 scientists and engineers from 64 research groups and companies across 20 countries characterized statistically equivalent samples from two different surfaces: a “rough” and a “smooth” surface. The results of the 2088 measurements constitute the most comprehensive surface description ever compiled. We find wide disagreement across measurements and techniques when the lateral scale of the measurement is ignored. Consensus is established through scale-dependent parameters while removing data that violates an established resolution criterion and deviates from the majority measurements at each length scale. Our findings suggest best practices for characterizing and specifying topography. The public release of the accumulated data and presented analyses enables global reuse for further scientific investigation and benchmarking.

42 ENGINEERING↗

Biomass Recalcitrance in Willow Under Two Biological Conversion Paradigms: Enzymatic Hydrolysis and Anaerobic Digestion

Biomass recalcitrance, the inherent resistance of plants towards deconstruction, negatively affects the viability of biorefineries. This trait is not only dictated by the properties of the biomass but also by the conversion system used and its interactions with specific features of the biomass. Here, biomass recalcitrance to anaerobic digestion (AD) was assessed using a biomethanation potential (BMP) assay. Plant material (n = 94) was selected from a large population of natural Salix viminalis accessions, previously evaluated for biomass recalcitrance using hydrothermal pretreatment-enzymatic hydrolysis. Correlations between yields from the two biological conversion systems were evaluated, as well as the influence of biomass compositional features, analyzed by pyrolysis-molecular beam mass spectrometry (py-MBMS), and other biomass physical properties on conversion performance. BMP values averaged 198.0 Nml CH 4 /g biomass after 94 days, ranging from 28.6 to 245.9. S lignin and carbohydrate-derived spectral features were positively correlated with performance under both systems, whereas G lignin, p-coumaric acid, and ferulic acid-derived ions were negatively correlated with yields and rates. Most spectral features were more strongly correlated with enzymatic hydrolysis yields compared to methane production. For early-stage methane production and rate, recalcitrance factors were similar compared to enzymatic hydrolysis, with weaker correlations observed at later timepoints. The results suggest that although variation in methane potential was considerably lower than enzymatic hydrolysis yields, a reduced recalcitrance under this system will still be of importance to improve early conversion rates. Spectral features of low methane-producing samples indicate the presence of inhibitory substances, warranting further study.

59 BASIC BIOLOGICAL SCIENCES↗

Diagnostic assessment of reservoir response to fracturing: a case study from Hydraulic Fracturing Test Site (HFTS) in Midland Basin

Abstract This paper outlines a data collection and diagnostics case study involving multiple horizontal shale wells. We look at well production profiles using rate transient analysis, differences in near wellbore complexity, geologic variations within the area of interest, as well as compositional differences in the rocks based on cores obtained from within the stimulated reservoir. The Hydraulic Fracturing Test Site is a multi-well experiment involving stimulation of unconventional shale wells in the southeastern Midland portion of the Permian Basin. The targeted formations include both the upper as well as the middle Wolfcamp formations, also referred alternatively as Wolfcamp A and Wolfcamp B. Data integration and analysis shared in this paper help us understand the various geologic controls impacting well productivity, particularly the wide variance observed between the Wolfcamp A and Wolfcamp B formations. Rate transient analysis indicates similar system permeabilities for stimulated wells. However, we observe higher effective fracture half-lengths for upper Wolfcamp wells. Using observations from 3D seismic interpretations (such as pad scale faults) as well as petrophysical and image log data, we highlight the substantial differences in stimulation as we move along the well laterals from the heel toward the toe sections. These differences are further reconciled with observations from zones with high data density at the core locations through stimulated rock, as well as independent data such as microseismic emissions. At the test site, Wolfcamp A was found to be relatively quartz rich with significant heterogeneity whereas Wolfcamp B is richer in clay and organic content. This impacts the geomechanical characteristics of the rock mass with much higher natural fracture density in the shallower interval. Thus, the fracture growth is more uniform in the deeper interval and more heterogeneous with branching likely in upper interval. Increased complexity also leads to consistently better productivity from the wells in the shallower interval as demonstrated from RTA results. This case study is unique because it provides valuable insights from actual sampling of the stimulated zones in hydraulically fractured wells and helps understand impact of various factors that contribute toward variability in well production. The findings from this study provides insights into need for optimization of completion designs in the various Wolfcamp landing zones, such as optimization of cluster or fracture spacing in various Wolfcamp intervals. In addition, it provides a useful template for data collection and research direction in future field test sites of similar nature in unconventional reservoirs.

Energy & Fuels↗

Interfacial X-Ray Scattering from Small Surfaces: Adapting Mineral-Fluid Structure Methods for Microcrystalline Materials

Crystal truncation rod (CTR) X-ray diffraction is an invaluable tool for measuring mineral surface and adsorbate structures, and has been applied to several environmentally and geochemically important systems. Traditionally, the method has been restricted to single crystals with lateral dimensions >3 mm. Minerals that meet this size criterion represent a minute fraction of those that are relevant to interfacial geochemistry questions, however. Crystal screening, data collection, and CTR measurement methods have been developed for crystals of <0.3 mm in lateral size using the manganese oxide mineral chalcophanite (ZnMn 3 O 7 •3H 2 O) as a case study. Furthermore, this work demonstrates the feasibility of applying the CTR technique to previously inaccessible surfaces, opening up a large suite of candidate substrates for future study.

58 GEOSCIENCES↗