Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Spatial dependence”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Novel Cell-Type-Specific Drought-Responsive Proteins in Root Tips of Field-Grown Perennial Switchgrass

The root-tip region of plants, including the root cap, forms the most basal terminal of the root and exhibits a high degree of cellular complexity in terms of morphology, cytological function, and interaction with environmental cues in the soil. Cells in this region follow a developmental trajectory, transitioning from stem cells to meristematic cells, and ultimately to fully differentiated cell types. However, our understanding of root-tip cell-type specific proteomic responses to abiotic stresses, such as drought, particularly under field conditions, remains limited. This study aimed to identify spatially resolved, cell type-specific proteomes in switchgrass (Panicum virgatum) root tips under drought stress. Root tips were collected from seven-year-old, field-grown switchgrass ‘Alamo’ plants excavated under both well-watered and long-term drought conditions. Cell type-specific proteins were identified using laser capture microdissection (LCM) coupled with nanoPOTS (Nanodroplet Processing in One Pot for Trace Samples) and nano-LC-MS proteomics analysis. Five distinct cell types were targeted: (1) cells in the quiescent center and stem cell niche (QuC), (2) protodermal epidermal cells (PEC) in the meristematic zone, (3) epidermal cells in the transition and elongation zones above the root cap (Epi), (4) peripheral root cap cells (PRC), forming 2–3 layers below the PEC and 1–2 layers above the root border cells, and (5) columella root cap cells (Col) comprising of the columella initials and a single underlying layer of cells undergoing active growth. Principal component analysis (PCA) revealed clear separation among the five targeted cell types, confirming distinct proteomic profiles. Proteins predominantly enriched in each cell type were linked to distinct cellular functions, with QuC cells showing involvement in chromosomal behavior, DNA replication, and mitosis—key processes for stem cell niche regulation. Drought stress resulted in alterations of proteostasis, as evidenced by significant decreases in ribosomal proteins and increases in protein synthesis inhibitors. Moreover, drought stress induced unique cell-type–specific proteins involved in phytohormone biosynthesis and signaling pathways, including auxin, cytokinin, and jasmonic acid. In particular, QuC cells were more highly enriched in proteins associated with DNA repair and mitotic processes. Metabolic pathways related to amino acids, carbohydrates, and lipids were differentially affected in a cell-type–dependent manner, whereas general stress-responsive proteins exhibited consistent changes across all five cell types. Overall, this study provides unique spatially resolved, cell-type-specific proteomic profiles in root tips, representing a significant advancement in our understanding of the cellular mechanisms underlying plant responses to drought stress in natural field conditions.

perennial grass↗
DOE OSTI2025

Plastoquinone redox status influences carboxysome integrity via a RpaA ‐ and reactive oxygen species‐dependent regulatory network

SUMMARY Carboxysomes are bacterial microcompartments that encapsulate Rubisco and are a core component of the cyanobacterial carbon concentration mechanism (CCM). While carboxysome number, size, and spatial organization vary in different environmental conditions (CO 2 , light availability, redox state, temperature, and light quality), the molecular mechanisms underlying this potentially adaptive process remain elusive. Herein, we observe that mutants of the circadian rhythm/metabolism factor, Regulator of Phycobilisome Association A (RpaA), exhibit a striking breakdown of carboxysomes under certain environmental conditions. We find that conditions leading to overreduction of the plastoquinone (PQ) pool (mixotrophic growth, high irradiance, or chemical inhibition of electron transfer from PQ to the cytochromeb 6 fcomplex) are accompanied by an elevated generation of reactive oxygen species (ROS) and correlate with the loss of carboxysome integrity. Carboxysome breakdown is reversed by environmental conditions or chemical inhibitors that prevent PQ overreduction and accompanying ROS generation. Taken together, our data support a novel link between the redox status of the PQ pool and carboxysome integrity. Our results have implications for the fundamental understanding of cyanobacterial energy‐balancing pathways and may indicate new research directions for understanding how the carboxysome is remodeled in response to changing environments.

Plant Sciences↗
DOE OSTI2021

Confronting lattice parton distributions with global QCD analysis

We present the first Monte Carlo based global QCD analysis of spin-averaged and spin-dependent parton distribution functions (PDFs) that includes nucleon isovector matrix elements in coordinate space from lattice QCD. We investigate the degree of universality of the extracted PDFs when the lattice and experimental data are treated under the same conditions within the Bayesian likelihood analysis. For the unpolarized sector, we find rather weak constraints from the current lattice data on the phenomenological PDFs, and difficulties in describing the lattice matrix elements at large spatial distances. In contrast, for the polarized PDFs we find good agreement between experiment and lattice data, with the latter providing significant constraints on the spin-dependent isovector quark and antiquark distributions.

Sato, Nobuo↗
DOE OSTI2020

Revisiting the first Sandia Fracture Challenge with transient deformation heating and strain localization considerations

In this study, we present a systematic experimental study and accompanying theoretical analysis of the dependence of ductile fracture on strain localization, strain hardening rates, deformation-induced thermal softening, and transient heat conduction under spatially uniform as well as spatially heterogeneous deformation. Spatially uniform cases are studied via standard dogbone shaped specimens subject to uniaxial tension. Spatially heterogeneous deformation is studied via the so-called Sandia Fracture Challenge (SFC) specimen, which is a standard compact tension specimen modified with three machined holes in front of a blunt notch. We utilize the same precipitation hardened martensitic stainless steel used in the first SFC experiments, i.e. 15-5 PH with an H1075 heat treatment. We also study 15-5 PH in Condition A and with the H900 heat treatment, each of which has a different hardening behavior and ductility. We find that the ductility does not correlate with the deformation to first crack initiation in the SFC specimen. Instead, hardening rates are better correlated. Moreover, a re-examination finds that an accurate calibration of the hardening rates is strongly correlated with accurate blind predictions of ductile fracture in the SFC specimens. Given that thermal softening can greatly affect hardening rates, we provided a thermographic analysis of deformation-induced heating and transient heat transfer in both the dogbone shaped samples and the SFC specimens. Transient heating in stainless steels is found to have first-order effects on the strain at the onset of necking, strain to fracture, as well as strain localization and ductile fracture in complex geometries even under so-called quasi-static loading rates.

42 ENGINEERING↗
DOE OSTI2025

Sylvester-preconditioned adaptive-rank implicit time integrators for advection-diffusion equations with variable coefficients

Here, we consider the adaptive-rank integration of multi-dimensional time-dependent advection-diffusion partial differential equations (PDEs) with variable coefficients. We employ a standard finite-difference method for spatial discretization coupled with high-order diagonally implicit Runge-Kutta temporal schemes. The discrete equation is a generalized Sylvester equation (GSE), which we solve with a projection-based adaptive-rank algorithm structured around two key strategies: (i) constructing dimension-wise subspaces using a novel atypical extended Krylov strategy, and (ii) efficiently solving the basis coefficient matrix with a preconditioned GMRES solver. The low-rank decomposition is performed in 2D using SVD and with high-order SVD (HOSVD) in 3D to represent the tensor in a compressed Tucker format. For d-dimensional problems (here, d = 2 or 3), the computational complexity and memory storage of the approach are found numerically to scale as and $\mathscr{O}(Nr^2) + \mathscr{O} (r^{d+1})$ and $\mathscr{O}(Nr) + \mathscr{O} (r^{d})$, respectively, with the one-dimensional resolution and the maximal rank during the Krylov iteration (which we find to be largely independent of on our numerical examples). We present numerical examples that illustrate the advertised properties of the algorithm.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

High-Resolution Measurement of Soil Thermal Properties and Moisture Content Using a Novel Heated Fiber Optics Approach

Hydrological parameters are scale dependent. Efficient monitoring techniques capable of measuring hydrological parameters, such as soil moisture content (θ), over a wide range of spatial scales are essential for understanding the complexity of water and energy movement across the landscape. Techniques to measure θ over spatial scales in the range from centimeters to thousands of meters, however, are sorely lacking. Recent improvements in the distributed temperature sensing (DTS) technology supported the development of novel techniques to fill that gap. However, improvements in the accuracy and applicability of DTS techniques are still needed. This study investigates the possibility of improving the accuracy of the fiber optics dual-probe heat-pulse (FO-DPHP) DTS technique by using a new design to maintain the spacing between the FO-DPHP probes and by introducing a novel data interpretation approach. The accuracy of the novel FO-DPHP design was tested at different θ in a sand column experiment. The FO-DPHP measurements obtained using traditional and novel data interpretation approaches were compared against independent measurements from several calibrated soil water content (EC5) sensors. Monte-Carlo analyses were also performed to assess the impact of DTS measurement errors on the accuracy achieved using the data interpretation approaches. The novel design and data interpretation approach allowed for accurate measurements of soil thermal properties and θ without the need to perform a hard-to-achieve soil-specific calibration. Measured θ had mean errors and standard deviations <0.03 and <0.01 m 3 m -3 , respectively, for moisture conditions ranging from dry to near saturation. The standard deviation in the measured heat capacity was <0.01 MJ m -3 K -1 .

42 ENGINEERING↗
DOE OSTI2022

Turbulent flow characteristics in an 84-pin rod bundle for typical and damaged spacer grids

Hexagonal rod bundles arranged in a tightly packed triangular lattice are extensively used for heat transfer and energy generation applications. Staggered spacer grids are used to maintain the structural integrity of gas-cooled fast reactor (GFR) fuel assemblies, while inducing localized turbulence in flow. Damage to these spacer grids results in a disruption of flow fields within these hexagonal fuel bundles. Experimental flow visualizations are critical to identify the differences in local flow properties that the structural damage may cause. This experimental research investigates the flow-field characteristics at a near-wall and center plane in a prototypical 84-pin GFR fuel assembly. Newly installed typical spacers and spacers subject to naturally occurring damage due to material degradation over prolonged experimentation were investigated. Velocity fields were acquired by utilizing the matched-index-of-refraction method to obtain time-resolved particle image velocimetry measurements for a Reynolds number of 12 000. Reynolds decomposition statistical results divulged differences in the time-averaged velocity, velocity fluctuations, flow anisotropy, and Reynolds stress distributions. Galilean decomposition demarcated the influence of spacer grid damage on the velocity fields. To extract turbulent structures and elucidate mechanisms of flow instabilities, proper orthogonal decomposition analysis was employed. Reduced order flow reconstructions enabled the application of vortex identification algorithms to determine the spatial and statistical characteristics of vortices generated. This research work provides unique experimental data on the spacer grid condition-dependent flow. The results offer a deeper understanding of fluid dynamics behavior to support GFR rod bundle design efforts and computational fluid dynamics model validation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A red-emitting carborhodamine for monitoring and measuring membrane potential

Biological membrane potentials, or voltages, are a central facet of cellular life. Optical methods to visualize cellular membrane voltages with fluorescent indicators are an attractive complement to traditional electrode-based approaches, since imaging methods can be high throughput, less invasive, and provide more spatial resolution than electrodes. Recently developed fluorescent indicators for voltage largely report changes in membrane voltage by monitoring voltage-dependent fluctuations in fluorescence intensity. However, it would be useful to be able to not only monitor changes but also measure values of membrane potentials. This study discloses a fluorescent indicator which can address both. We describe the synthesis of a sulfonated tetramethyl carborhodamine fluorophore. When this carborhodamine is conjugated with an electron-rich, methoxy (-OMe) containing phenylenevinylene molecular wire, the resulting molecule, CRhOMe, is a voltage-sensitive fluorophore with red/far-red fluorescence. Using CRhOMe, changes in cellular membrane potential can be read out using fluorescence intensity or lifetime. In fluorescence intensity mode, CRhOMe tracks fast-spiking neuronal action potentials (APs) with greater signal-to-noise than state-of-the-art BeRST 1 (another voltage-sensitive fluorophore). CRhOMe can also measure values of membrane potential. The fluorescence lifetime of CRhOMe follows a single exponential decay, substantially improving the quantification of membrane potential values using fluorescence lifetime imaging microscopy (FLIM). The combination of red-shifted excitation and emission, mono-exponential decay, and high voltage sensitivity enable fast FLIM recording of APs in cardiomyocytes. The ability to both monitor and measure membrane potentials with red light using CRhOMe makes it an important approach for studying biological voltages.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Spatial Profile and Beam Composition Measurements of the University of Wisconsin Inertial Electrostatic Confinement Laboratory Ion Guns

A custom designed and manufactured set of ion guns has been in use at the University of Wisconsin Inertial Electrostatic Confinement Laboratory for both beam fusion experiments and materials implantation experiments. For the first time, direct measurements have been made on the spatial profiles and the mass compositions of He and D ion beams produced by these guns. The results validate assumptions about the circular Gaussian spatial profiles for both He and D ion beams. Mass composition measurements of the He beam identified a pressure-dependent minimum impurity content of 15% N + . The D beam contained relative molecular ion fractions of 58% D 3 + , 32% D 2 + , and 10% D + with impurities of 15% to 20% D 2 O + . Finally, a new experimental platform, the Ion Beam and Source Analyzer was developed to perform these experiments on the ion guns used to irradiate candidate fusion materials.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Aberration-corrected STEM imaging of 2D materials: Artifacts and practical applications of threefold astigmatism

High-resolution scanning transmission electron microscopy (HR-STEM) with spherical aberration correction enables researchers to peer into two-dimensional (2D) materials and correlate the material properties with those of single atoms. The maximum intensity of corrected electron beam is confined in the area having sub-angstrom size. Meanwhile, the residual threefold astigmatism of the electron probe implies a triangular shape distribution of the intensity, whereas its tails overlap and thus interact with several atomic species simultaneously. The result is the resonant modulation of contrast that interferes the determination of phase transition of 2D materials. Here, we theoretically reveal and experimentally determine the origin of resonant modulation of contrast and its unintended impact on violating the power-law dependence of contrast on coordination modes between transition metal and chalcogenide atoms. The finding illuminates the correlation between atomic contrast, spatially inequivalent chalcogenide orientation, and residual threefold astigmatism on determining the atomic structure of emerging 2D materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Classical Benchmarks for Variational Quantum Eigensolver Simulations of the Hubbard Model

Simulating the Hubbard model is of great interest to a wide range of applications within condensed matter physics, however its solution on classical computers remains challenging in dimensions larger than one. The relative simplicity of this model, embodied by the sparseness of the Hamiltonian matrix, allows for its efficient implementation on quantum computers, and for its approximate solution using variational algorithms such as the variational quantum eigensolver. While these algorithms have been shown to reproduce the qualitative features of the Hubbard model, their quantitative accuracy in terms of producing true ground state energies and other properties, and the dependence of this accuracy on the system size and interaction strength, the choice of variational ansatz, and the degree of spatial inhomogeneity in the model, remains unknown. Here we present a rigorous classical benchmarking study, demonstrating the potential impact of these factors on the accuracy of the variational solution of the Hubbard model on quantum hardware, for systems with up to 32 qubits. We find that even when using the most accurate wavefunction ansätze for the Hubbard model, the error in its ground state energy and wavefunction plateaus for larger lattices, while stronger electronic correlations magnify this issue. Concurrently, spatially inhomogeneous parameters and the presence of off-site Coulomb interactions only have a small effect on the accuracy of the computed ground state energies. Our study highlights the capabilities and limitations of current approaches for solving the Hubbard model on quantum hardware, and we discuss potential future avenues of research.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Microstructural Topology as a Prescriptor for Quantum Coherence: Towards A Unified Framework for Decoherence in Superconducting Qubits

In superconducting quantum circuits, decoherence improvements are frequently obtained through process interventions that simultaneously modify surface chemistry, microstructural topology, and device geometry, leaving mechanistic attribution structurally underdetermined. Predictive materials engineering requires measurable structural statistics to be separated from geometry-dependent coupling coefficients into independently testable factors. We introduce the concept of classical and quantum microstructure. In that context, we formulate a channel-wise separable framework for decoherence in superconducting transmon qubits in which each loss channel is described by a reduced prescriptor. Here, a channel-specific microstructural state variable is determined independently of device geometry, and a geometry-dependent coupling functional is computable from field solutions without reference to surface chemistry. We derive this product form from a spatially resolved kernel representation and establish a perturbative separability criterion that defines the regime where independent variation of the variables is valid. The framework specifies five prescriptor classes for dominant loss pathways in transmon-class devices. Falsifiability is operationalized through a pre-committed 2x2 experimental protocol in which the variables must satisfy independent ratio checks within propagated uncertainty. A Minimum-Dataset Specification standardizes reporting for cross-laboratory inference. Part I establishes the conceptual and mathematical architecture; coordinated experimental validation is reserved for Part II.

Dravid, Vinayak P. [Northwestern U.]↗
DOE OSTI2022

Radio-Frequency Rapid Thermal Processing Enabling Spatial Phase Transformation and Nanocrystallization of Soft Magnetic Amorphous Alloys

Thermal processing of soft magnetic amorphous and nanocrystalline alloys is explored under the influence of radio-frequency induction-heating techniques. Direct induction-heating concepts based on longitudinal and transverse flux heating are examined and the details of electromagnetic fields interaction with metallic strips are discussed by analytical calculations as well as finite element analysis. Initial experimental results confirming spatial control of phase transformations and nanocrystallization within a single strip of Finemet Fe-based amorphous ribbons are reported. Here, the degree to which primary and secondary crystallization temperature are achieved depends on the spacing between the ribbon relative to the induction coil as well as the coil design and configuration. For transverse coil configurations, the local temperature and therefore microstructural evolution is different across the lateral dimension of processed ribbons, with reduced gap sizes producing enhanced peak temperatures and larger temperature distributions with greater spatial variation in microstructure. In addition, indirect susceptor-based induction heating under tension is performed and the impact of microstructure is demonstrated. Herein, potential for exploiting spatially optimized phase transformations is illustrated through electromagnetic field–assisted processing in a scalable manufacturing process with amorphous and nanocrystalline soft magnetic alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Generic two-phase coexistence in a type-2 Schloegl model for autocatalysis on a square lattice: Analysis via heterogeneous master equations

Schloegl's second model (also known as the quadratic contact process) on a square lattice involves spontaneous annihilation of particles at lattice sites at rate p, and their autocatalytic creation at unoccupied sites with n ≥ 2 occupied neighbors at rate k n . Here, kinetic Monte Carlo (KMC) simulation reveals that these models exhibit a nonequilibrium discontinuous phase transition with generic two-phase coexistence: the p value for equistability of coexisting populated and vacuum states, p eq (S), depends on the orientation or slope, S, of a planar interface separating those phases. The vacuum state displaces the populated state for p > eq (S), and the opposite applies for p < p eq (S) for 0 < S < ∞. The special “combinatorial” rate choice k n = n(n–1)/12 facilitates an appealing simplification of the exact master equations for the evolution of spatially heterogeneous states in the model, which aids analytic investigation of these equations via hierarchical truncation approximations. Truncation produces coupled sets of lattice differential equations which can describe orientation-dependent interface propagation and equistability. The pair approximation predicts that p eq (max) = p eq (S = 1) = 0.09645 and p eq (min) = p eq (S → ∞) = 0.08827, values deviating less than 15% from KMC predictions. In the pair approximation, a perfect vertical interface is stationary for all p < p eq (S = ∞) = 0.08907, a value exceeding p eq (S → ∞). One can regard an interface for large S → ∞ as a vertical interface decorated with isolated kinks. For p < p eq (S = ∞), the kink can move in either direction along this otherwise stationary interface depending upon p, but for p = p eq (min) the kink is also stationary.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Assessing Cloud and Precipitation Properties on Temporal and Spatial Scales Using LASSO Simulations over ENA

Low clouds and precipitation representation remain a major source of uncertainty in Earth System Models (ESMs), particularly due to challenges in representing their sub-grid variability and scale-dependent sampling. This study evaluates the performance of preliminary simulations from the Large-Eddy Simulation (LES) ARM Symbiotic Simulation and Observation (LASSO) project over the Eastern North Atlantic (ENA), with a focus on liquid water path (LWP), ice water path (IWP), cloud fraction (CF), and surface precipitation simulated across closed-cell, open-cell, and transitional cloud regimes. Using LES (100 m horizontal grid spacing) driven by ERA5 and MERRA-2 reanalyses, we assess the representativeness of ground-based point observations by analyzing their correspondence to model-resolved spatial and temporal means. Results suggest that observational sampling of at least 6 hours is required to achieve consistency with domain-scale averages, in particular for observations that exhibit pronounced sub-grid heterogeneity, such as precipitation. ERA5-forced simulations exhibit improved spatial coherence and agreement with domain-averaged quantities when compared to MERRA-2 runs, with performance discrepancies largest for convective cloud conditions due to differences in forcing fidelity and temporal resolution. These findings highlight the importance of regime-aware model evaluation strategies and potentially demonstrate how LES can inform observation-model comparison practices and the development of cloud and precipitation parameterizations in ESMs.

Liang, Jiakun [University of Hawai'i at Manoa] (OR↗
DOE OSTI2023

Spatially anisotropic 𝑆=1 square-lattice antiferromagnet with single-ion anisotropy realized in a Ni(II) pyrazine-𝑛,𝑛′-dioxide coordination polymer

The Ni(NCS) 2 (pyzdo) 2 coordination polymer is found to be an S = 1 spatially anisotropic square lattice with easy-axis single-ion anisotropy. This conclusion is based upon considering in concert the experimental probes x-ray diffraction, magnetic susceptibility, magnetic-field-dependent heat capacity, muon-spin relaxation, neutron diffraction, neutron spectroscopy, and pulsed-field magnetization. Long-range antiferromagnetic (AFM) order develops at T N = 18.5 K. Although the samples are polycrystalline, there is an observable spin-flop transition and saturation of the magnetization at ≈80 T. Linear spin-wave theory yields spatially anisotropic exchanges within an AFM square lattice, J x =0.235 meV, J y =2.014 meV, and an easy-axis single-ion anisotropy D=−1.622 meV (after renormalization). The anisotropy of the exchanges is supported by density functional theory.

2-dimensional systems↗
DOE OSTI2020

In situ X-ray spatial profiling reveals uneven compression of electrode assemblies and steep lateral gradients in lithium-ion coin cells

Coin cells are used extensively as test devices in battery research for evaluation of new materials and optimization of cycling protocols. Here, in situ X-ray diffraction profilometry is used to characterize spatial distribution of the active materials, lithiation, and phase distribution in electrodes of NCM523/graphite coin cells. The X-ray data indicate uneven areal compression of the electrode assembly in such cells, which we trace to a specific design feature that leads to elastic deformation of a metal spacer. Steep lithiation gradients observed in the electrodes imply radially-dependent resistivity, for which uneven compression of the separator is a likely cause. Electrochemical model calculations suggest that variable porosity of the polymer separator would account for the salient features of spatial profiles observed in these coin cells.

25 ENERGY STORAGE↗
DOE OSTI2020

Factors controlling the oxygen isotopic composition of lacustrine authigenic carbonates in Western China: implications for paleoclimate reconstructions

Abstract In the carbonate-water system, at equilibrium, the oxygen isotopic composition of carbonate is dependent not only on the temperature but also on the isotopic composition of host water in which the carbonate is formed. In this study, lake surface sediment and water samples were collected from 33 terminal lakes in Western China to evaluate controls on the oxygen isotopic composition of lacustrine authigenic carbonates (δ 18 O carb ) and its spatial distribution. Our results show that water oxygen isotopic composition (δ 18 O water ) rather than lake summer water temperature (T water ), is the main determinant of δ 18 O carb , irrespective of whether oxygen isotope equilibrium is achieved. There are significant linear correlations between δ 18 O carb and elevation, as well as that between δ 18 O carb and latitude for lakes located on the Tibetan Plateau. In Western China, the spatial distribution of δ 18 O carb is consistent with that of δ 18 O water , and is ultimately controlled by the isotopic composition of local precipitation (δ 18 O precipitation ) that depends on the source of water vapor. Therefore, changes in δ 18 O carb can be predominantly interpreted as variations of δ 18 O water , which in turn represent changes in δ 18 O precipitation for paleoclimate reconstructions in this region, and may be relevant to studies of other areas.

59 BASIC BIOLOGICAL SCIENCES↗