Search NASA⌕ Search

SEARCH · Search NASA

Results for “variable density”

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 1,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Evaluating the Influence of CAM5 Aerosol Configuration on Simulated Tropical Cyclones in the North Atlantic

This study examines the influence of prescribed and prognostic aerosol model configurations on the formation of tropical cyclones (TCs) in the North Atlantic Ocean in Community Atmosphere Model version 5 (CAM5). The impact of aerosol parameterization is examined by investigating storm track density, genesis density, potential intensity, and genesis potential index. This work shows that both CAM5 configurations simulate reduced storm frequency when compared to observations and that differences in TC climatology between the model configurations can be explained by differences in the large-scale environment. The analysis shows that simulation with the prognostic aerosol parameterization scheme reasonably captures the observed interannual variability in tropical cyclones and aerosols (i.e., dust) in the North Atlantic, while simulation with the prescribed configuration (climatology) is less favorable. The correlation between dust and TCs in observations (i.e., reanalysis and satellite datasets) is shown to be negative, and this relationship was also found for the prognostic aerosol configuration despite an overall decrease in the frequency of TCs. This indicates that, to accurately replicate certain aspects of TC interannual variability, the aerosol configuration within CAM5 needs to account for the appropriate dust variability.

Huff, J. Jacob A. (ORCID:0000000329143817)↗
DOE OSTI2016

Observations of the lensed quasar Q2237+0305 with Canaricam at GTC

We present new mid-IR observations of the quadruply lensed quasar Q2237+0305 taken with CanariCam on the Gran Telescopio Canarias. Mid-IR emission by hot dust, unlike the optical and near-IR emission from the accretion disk, is unaffected by the interstellar medium (extinction/scattering) or stellar microlensing. We compare these “true” ratios to the (stellar) microlensed flux ratios observed in the optical/near-IR to constrain the structure of the quasar accretion disk. We find a half-light radius of R{sub 1/2}=3.4{sub −2.1}{sup +5.3}√(〈M〉/0.3M{sub ⊙}) lt-day at λ{sub rest}=1736 A-ring and an exponent for the temperature profile R∝λ{sup p} of p=0.79±0.55, where p=4/3 for a standard thin-disk model. If we assume that the differences in the mid-IR flux ratios measured over the years are due to microlensing variability, we find a lower limit for the size of the mid-IR-emitting region of R{sub 1/2}≳200 √(〈M〉/0.3M{sub ⊙}) lt-day. We also test for the presence of substructure/satellites by comparing the observed mid-IR flux ratios with those predicted from smooth lens models. We can explain the differences if the surface density fraction in satellites near the lensed images is α=0.033{sub −0.019}{sup +0.046} for a singular isothermal ellipsoid plus external shear mass model or α=0.013{sub −0.008}{sup +0.019} for a mass model combining ellipsoidal NFW and de Vaucouleurs profiles in an external shear.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Typology of Cities for Systems and Modeling for Accelerated Research in Transportation (SMART) Mobility Consortium

Mobility patterns, technology adoption and associated energy outcomes vary tremendously across settlement types. This poster discusses how a highly geographically resolved exploration of the social, economic, techno-infrastructural and environmental domains is key to understanding observed variations in transportation technology adoption and associated mobility and energy outcomes. Current socioeconomic and mobility data sets at the census block group level are integrated in a hierarchical clustering approach to show how variations in mobility and energy outcomes are shaped by these domains - which may enhance progress on investments, or effectively inform planning and other decisions for state-wide goals. In a pilot analysis in New York, the clustering produces four settlement types to predict dependent variables of electric vehicle (EV) adoption rates, commute mode, vehicle fuel economy, and vehicles per household. This typology shows EV adoption rates among the core urban population, which is wealthy and highly-educated, are high - 3 EVs/1,000 vehicles versus 1 EV/1,000 vehicles among the urban working class. Commuting mode is closely linked with population and employment density - more than 90% of core urbanites use transit or active modes, compared with only 22% of suburbanites and 17% of rural residents. Household vehicle ownership also varies, with two vehicles per household in rural areas and only 0.5 in core urban settings.Important findings on differences among the rural, suburban, urban, and urban core settlement types suggest a need to explore how best to manage and anticipate very different types of services that may be supportive in achieving energy-efficient and affordable mobility systems state-wide.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2026

Dynamical Dark Energy Imprints in the Lyman-Alpha Forest

The nature of dark energy (DE) remains elusive, even though it constitutes the dominant energy-density component of the Universe and drives the late-time acceleration of cosmic expansion. By combining measurements of the expansion history from baryon acoustic oscillations, supernova surveys, and cosmic microwave background data, the Dark Energy Spectroscopic Instrument (DESI) Collaboration has inferred that the DE equation of state may evolve over time. The profound implications of a time-variable, ``dynamical" DE (DDE) that departs from a cosmological constant motivate the need for independent observational tests. In this work, we use cosmological hydrodynamical simulations of structure formation to investigate how DDE affects the properties of the Lyman-Alpha ``forest'' of absorption features produced by neutral hydrogen in the cosmic web. We find that DDE models consistent with the DESI constraints induce a spectral tilt in the forest transmitted flux power spectrum, imprinting a scale- and redshift-dependent signature relative to standard Lambda-CDM cosmologies. These models also yield higher intergalactic medium temperatures and reduced Lyman-Alpha opacity compared to Lambda-CDM. We discuss the observational implications of these trends as potential avenues for independent confirmation of DDE.

Garza, Diego [UC, Santa Cruz, Inst. Part. Phys.] (↗
DOE OSTI2026

Photophysical Properties and Phase Behavior of Ultrawide Photovoltaic Bandgap Cesium–Lead-Based Triple Halide Perovskites

Metal halide perovskite films in the top cell of triple-junction tandems require bandgaps around 2.0 eV to achieve current matching, assuming that the middle absorbing layer is the commonly used FAPbI 3 composition and the bottom cell has a bandgap around 1.1 eV. Unfortunately, mixed organic/inorganic metal halide perovskites that have the necessary Br content to reach a bandgap of 2.0 eV segregate into iodine-rich and bromine-rich phases under illumination, limiting their obtainable voltage. Previous reports have shown improved photostability using either Cs-based inorganic compositions or Cl incorporation on the X-site. Here, we investigate the inorganic triple halide compositional space CsPb­(I 1–x–y Br y Cl x ) 3 where bandgaps near 2.0 eV are expected based on the knowledge that CsPbI2Br has a bandgap of 1.90 eV. Incorporation of Cl occurs readily for x ≤ 0.07–0.10 within perovskites with a Br content of 0.3 ≤ y ≤ 0.42. When x >0.1, X-ray diffraction and photoluminescence (PL) measurements indicate that multiple compositional phases form. We hypothesize that the variable sizes of the three halide ions are not supported within the rigid Cs lattice, resulting in the formation of multiple compositional phases. The photoluminescence quantum yield of the single-phase compositional space–CsPb­(I 1–x–y Br y Cl x ) 3 where x ≤ 0.07was typically 0.001–0.004%, most likely as a result of a high defect density, including mobile iodine species. PL light-soaking measurements of many perovskite compositions with bandgaps in the range of 1.89–2.05 eV demonstrate that phase segregation occurs when initial bandgaps are above 1.95 eV regardless of halide content: indicating further iodide oxidation and corresponding migration under illumination. The conclusion is that further compositional or additive engineering is necessary for the development of inorganic triple halide compositions that accomplish the elusive goal of fabricating high-quality and photostable 2.0 eV films for use in multijunction tandems.

Electrical conductivity↗
DOE OSTI2020

Subsurface Nitrogen Dissociation Kinetics in Lithium Metal from Metadynamics

The dissociation of molecular nitrogen in lithium is of interest for several promising technologies, such as the catalytic synthesis of ammonia in ambient or mild conditions. In this work we simulate nitrogen dissociation in the lithium BCC (110) surface at ambient and elevated temperatures using density functional theory (DFT) metadynamics simulations. The rate constants at temperatures of 300, 400, and 500 K are calculated by statistical analysis of the reaction time distributions from the accelerated simulations. This approach finds and estimates rate constants for transition pathways out of the initial state; the required input is the stable initial state and a reasonable choice of collective variable. A single collective variable is used in this case: the N–N distance. The results are robust to changes in metadynamics parameters, and the reaction time distributions follow the expected exponential distribution. We show that the metadynamics-derived rate constants are in agreement with results from the conventional harmonic approximation approach using a climbing image nudged elastic band (NEB) transition state search. The reaction barriers from metadynamics and the NEB/harmonic approximation agree to within 0.02–0.04 eV at all temperatures studied. This paper demonstrates that the harmonic approximation provides an accurate description of the rate constants for nitrogen dissociation in lithium metal, even at temperatures near or above the melting point of lithium, lending credence to previous and future theoretical studies using this approximation. Moreover, this work demonstrates a step toward the automated exploration and discovery of reaction mechanisms and associated rate constants for elementary surface-catalyzed reactions using DFT-based metadynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Atomic-Scale Tracking of Topological Defect Motion and Incommensurate Charge Order Melting

Charge order pervades the phase diagrams of quantum materials where it competes with superconducting and magnetic phases, hosts electronic phase transitions and topological defects, and couples to the lattice generating intricate structural distortions. Incommensurate charge order is readily stabilized in manganese oxides, where it is associated with anomalous electronic and magnetic properties, but its nanoscale structural inhomogeneity complicates precise characterization and understanding of its relationship with competing phases. Leveraging atomic-resolution variable-temperature cryogenic scanning transmission electron microscopy, we characterize the thermal evolution of charge order as it transforms from its ground state in a model manganite system. We find that mobile networks of discommensurations and dislocations generate phase inhomogeneity and induce global incommensurability in an otherwise lattice-locked modulation. Driving the order to melt at high temperatures, the discommensuration density grows, and regions of order locally decouple from the lattice periodicity. Published by the American Physical Society 2025

Schnitzer, Noah (ORCID:0000000210421146)↗
DOE OSTI2025

A general kinetic framework for dislocation mobility derived from probabilistic cellular automaton simulations of the kink-pair mechanism

Dislocation mobility laws are essential components of dislocation-density-based crystal plasticity models. For dislocations governed by the kink-pair mechanism, however, existing formulations are typically limited to specific regimes due to the com plex interplay between stochastic kink-pair nucleation and lateral kink migration. In this work, we develop a general kinetic framework that expresses the average dislocation velocity as a function of mechanism-level variables: positive/negative kink pair nucleation rates, kink migration velocity, dislocation segment length, critical kink-pair width, and kink height. Probabilis tic cellular automaton simulations are used to capture the behavior of conceptual dislocation segments between the limiting conditions of migration outpacing nucleation on the one end and nucleation outpacing migration on the other. An elemen tary functional form that captures the system dynamics is suggested and fitted against the simulation results. This framework remains valid for arbitrary combinations of the six variables and is, therefore, compatible with any admissible constitutive re lations that describe their stress and temperature dependence. Comparisons with established approaches and experimental results confirm the robustness and physical consistency of the formulation, making it broadly applicable to material systems in which dislocation motion is governed by the kink-pair mechanism.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Strain-induced optical band gap variation of SnO 2 films

In this paper, thickness dependent strain relaxation effects are utilized to study the impact of crystal anisotropy on the optical band gap of epitaxial SnO 2 films grown by pulsed laser deposition on (0001)-oriented sapphire substrates. An X-ray diffraction analysis reveals that all films are under tensile biaxial in-plane strain and that strain relaxation occurs with increasing thickness. Variable angle spectroscopic ellipsometry shows that the optical band gap of the SnO 2 films continuously increases with increasing film thickness. This increase in the band gap is linearly related to the strain state of the films, which indicates that the main origin of the band gap change is strain relaxation. The experimental observation is in excellent agreement with results from density functional theory for biaxial in-plane strain. Our research demonstrates that strain is an effective way to tune the band gap of SnO 2 films and suggests that strain engineering is an appealing route to tailor the optical properties of oxide semiconductors.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Surface chemical reactivity of ultrathin Pd(111) films on Ru(0001): Importance of orbital symmetry in the application of the d -band model

The chemical bonding of adsorbate molecules on transition-metal surfaces is strongly influenced by the hybridization between the molecular orbitals and the metal d-band. The strength of this interaction is often correlated with the location of the metal d-band center relative to the Fermi level. Here, we exploit finite size effects in the electronic structure of ultrathin Pd(111) films grown on Ru(0001) to tune their reactivity by changing the film thickness one atom layer at a time, while keeping all other variables unchanged. Interestingly, while bulk Pd(111) is reactive toward oxygen, Pd(111) films below five monolayers are surprisingly inert. This observation is fully in line with the d-band model prediction when applied to the orbitals involved in the bonding. The shift of the d-band center with film thickness is primarily attributed to shifts in the partial density of states associated with the 4d xz and 4d yz orbitals. This study provides an in-depth look into the orbital specific contributions to the surface chemical reactivity, providing new insights that could be useful in surface catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Interfacial Charge Transfer Pathways in Photoelectrochemical H 2 Evolution by a Single-Component Molecular Catalyst on a Conductive Metal Oxide

Dye-sensitized photoelectrosynthesis cells traditionally combine a photosensitizing dye to harvest light and a catalyst to generate chemical fuels on a semiconductor. Here, in this work, a photoactive catalyst capable of both light absorption and fuel formation, [Cp*Ir(4,4′-Y 2 -bpy)Cl][Cl] (Cp* = pentamethylcyclopentadienyl, bpy = 2,2′-bipyridine, Y = CH 2 PO 3 H 2 ), is anchored to a mesoporous tin-doped indium oxide (ITO) electrode and facilitates photoelectrochemical H 2 evolution in water without the need for additional photosensitizers or sacrificial reductants. Spectroelectrochemistry indicates a single-site H 2 evolution mechanism involving charge injection to ITO, in contrast to the bimetallic mechanism observed in solution. Cyclic voltammetry and variable-potential chronoamperometry under illumination probe competing pathways via interfacial electron transfer between the Ir hydride excited state and the conductive ITO electrode. A Marcus theory framework provides reorganization energies for competing interfacial electron transfer pathways from potential-dependent quantum yield measurements. The light-driven H 2 evolution catalysis on conducting oxide proceeds with high Faradaic efficiency and current densities comparable to photoelectrodes utilizing p-type semiconductors. By uncovering the principal electron transfer pathways that govern photocatalytic efficiency, this study establishes design principles for hybrid molecular photoelectrocatalysts.

catalysts↗
DOE OSTI2021

Coordination and activation of nitrous oxide by iron zeolites

Iron-containing zeolites are heterogeneous catalysts that exhibit remarkable activity in the selective oxidation of inert hydrocarbons and catalytic decomposition of nitrous oxide (N 2 O). The reduction of N 2 O is critical to both these functions, but experimental data tracking the iron active sites during N 2 O binding and activation are limited. Here, the N 2 O-ligated Fe(ii) active site in iron-exchanged zeolite beta is isolated and characterized by variable-temperature Mössbauer, diffuse reflectance UV-vis-NIR and Fourier transform infrared spectroscopy. N 2 O binds through the terminal nitrogen atom with substantial backbonding from the Fe(ii) centre at low temperature. At higher temperatures, the Fe–N 2 O interaction is weakened, facilitating isomerization to the O-bound form, which is competent in O-atom transfer. Density functional theory calculations show the geometric and electronic structure requirements for N 2 O binding and activation. A geometric distortion imposed by the zeolite lattice plays an important role in activating N 2 O. As a result, this highlights a mechanism for structural control over function in Fe-zeolite catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Renormalized site density functional theory

Site density functional theory (SDFT) provides a rigorous framework for statistical mechanics analysis of inhomogeneous molecular liquids. The key defining feature of these systems is the presence of two very distinct interactions scales (intra- and inter-molecular), and as such proper description of both effects is critical to the accuracy of the calculations. Current SDFT applications utilize the same approximation scheme for both interaction motifs, which negatively impacts the results. Dual space methodology, used in this work, alleviates this issue by providing the flexibility of evaluating part of the interactions in traditional field representation. For molecular liquid this translates into retaining density representation for inter-molecular interactions but describing stiff intra-molecular remainder by more appropriate conventional field based methods. This opens the way to decouple analysis of the two interactions scales - the idea which is developed further in this work for the case of homogeneous reference approximation of inter-molecular interactions. We demonstrate that by defining new collective variables, the behavior of the original molecular liquid system at the inter-molecular level can be transformed to resemble that of an effective simple fluid mixture. The latter is linked to intra-molecular scale through renormalized interaction parameters. We illustrate this renormalization procedure for several types of diatomic liquids, showing that this approach cures many of the shortcomings of existing SDFT methods.

Chuev, Gennady N.↗
DOE OSTI2025

Metagenome-assembled genomes from topsoils collected during NEON campaign in East River, CO (06/14/2018-06/28/2018)

The Watershed Function Science Focus Area (WF SFA) at Lawrence Berkeley National Lab is working to build a mechanistic understanding of the distribution and dynamics of biogeochemical processes in mountainous watersheds and their response to perturbation. In June 2018, the NEON (National Ecological Observatory Network) Airborne Observatory Platform (AOP) performed a taskable airborne imaging campaign to collect visible to shortwave infrared (VSWIR) imaging spectroscopy and LiDAR data across 330 km2 in the Upper East River at Crested Butte, CO. We conducted a parallel ground sampling campaign to sample vegetation traits, as well as soil physical, chemical, and microbiological characteristics. We collected these samples from 438 sites across 12 locations spanning much of the elevation, topographic, and geologic variability across the study area. A subset of 250 samples were used for soil metagenomics which is presented here. In addition, at each site, vegetation samples were collected to measure species-specific leaf water content and leaf mass area, foliar elemental composition and foliar CN stable isotope ratios. Soil samples were collected to measure soil physical properties which include bulk density and soil texture analysis. A suite of soil chemical properties was measured from the samples collected at each site, including pH, organic matter, concentrations exchangeable cations, total elemental composition, and the concentrations of extractable N pools (e.g. total free amino acids, ammonium, nitrate, dissolved organic N, and total dissolved N). Additionally, we have measured soil microbial biomass CN stoichiometry. Here, we present 1982 metagenome-assembled genomes (MAGs) for the bacterial and archaeal community from topsoil collected from during NEON 2018 campaign. All metagenomes were sequenced at JGI (Joint Genome Institute) (GOLD Study ID: Gs0149986). Metagenomes were assembled using JGI Metagenome Workflow (10.1128/mSystems.00804-20). The dataset includes (1) zip files for 1982 MAG fasta files (neon_genomes1-5.tar.gz, split into 5 tarballs to keep tarballs under 0.5 GB), (2) neon_Gs0149986_samples_soilproperties_metagenomes.csv: the sample information together with the accession numbers for the underlying metagenomes and the associated soil physical and chemical measurements in NMDC (National Microbiome Data Collaborative) compliant format, (3) neon_Gs0149986.kml: location bounding box file for the sampled locations, (4) samples.csv: sample metadata file used to register Internationall Generic Sample Numbers (IGSNs), (5) flmd.csv: file level metadata file, and (6) dd.csv: data dictionary file. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

2018 NEON and 2025 CHESS Campaigns↗
DOE OSTI2023

A Modeling Study on Ammonia and Ammonia/Hydrogen Kinetics for Gas Turbine Engines

The use of ammonia as a fuel source in gas turbine engine power cycles represents an attractive means to decarbonize the energy sector due to higher energy density and achieving liquid state at far lower pressures compared to pure hydrogen. However, due to low flammability and a propensity for high NOx emissions, its use is not without challenge. Here, a number of 0D and 1D modeling tools were utilized to study the combustion characteristics of ammonia and ammonia/hydrogen mixtures, examining basic fundamental properties such as laminar flame speed, variability among existing chemical kinetic mechanisms, and considering the use of two-stage rich-lean combustion strategies to achieve low NOx emissions.

42 ENGINEERING↗
DOE OSTI2018

Robust Period Estimation Using Mutual Information for Multiband Light Curves in the Synoptic Survey Era

The Large Synoptic Survey Telescope (LSST) will produce an unprecedented amount of light curves using six optical bands. Robust and efficient methods that can aggregate data from multidimensional sparsely sampled time-series are needed. Here, we present a new method for light curve period estimation based on quadratic mutual information (QMI). The presented method does not assume a particular model for the light curve nor its underlying probability density and it is robust to non-Gaussian noise and outliers. By combining the QMI from several bands the true period can be estimated even when no single-band QMI yields the period. Period recovery performance as a function of average magnitude and sample size is measured using 30,000 synthetic multiband light curves of RR Lyrae and Cepheid variables generated by the LSST Operations and Catalog simulators. The results show that aggregating information from several bands is highly beneficial in LSST sparsely sampled time-series, obtaining an absolute increase in period recovery rate up to 50%. We also show that the QMI is more robust to noise and light curve length (sample size) than the multiband generalizations of the Lomb–Scargle and AoV periodograms, recovering the true period in 10%–30% more cases than its competitors. A python package with efficient Cython implementations of the QMI and other methods is provided.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Modeling the distribution of the endangered Jemez Mountains salamander (Plethodon neomexicanus) in relation to geology, topography, and climate

The Jemez Mountains salamander (Plethodon neomexicanus; hereafter JMS) is an endangered salamander restricted to the Jemez Mountains in north-central New Mexico, United States. This strictly terrestrial and lungless species requires moist surface conditions for activities such as mating and foraging. Threats to its current habitat include fire suppression and ensuing severe fires, changes in forest composition, habitat fragmentation, and climate change. Forest composition changes resulting from reduced fire frequency and increased tree density suggest that its current aboveground habitat does not mirror its historically successful habitat regime. However, because of its limited habitat area and underground behavior, we hypothesized that geology and topography might play a significant role in the current distribution of the salamander. We modeled the distribution of the JMS using a machine learning algorithm to assess how geology, topography, and climate variables influence its distribution. The best habitat suitability model indicates that geology type and maximum winter temperature (November to March) were most important in predicting the distribution of the salamander (23.5% and 50.3% permutation importance, respectively). Minimum winter temperature was also an important variable (21.4%), suggesting this also plays a role in salamander habitat. Our habitat suitability map reveals low uncertainty in model predictions, and we found slight discrepancies between the designated critical habitat and the most suitable areas for the JMS. Because geological features are important to its distribution, we recommend that geological and topographical data are considered, both during survey design and in the description of localities of JMS records once detected.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

A reduced order model to analytically infer atmospheric CO 2 concentration from stomatal and climate data

To address questions related to the acceleration or deceleration of the global hydrological cycle or links between the carbon and water cycles over land, reliable data for past climatic conditions based on proxies are required. In particular, the reconstruction of palaeoatmospheric CO 2 content (C a ) is needed to assist the separation of natural from anthropogenic C a variability and to explore phase relations between C a and air temperature T a time series. Both T a and C a are needed to fingerprint anthropogenic signatures in vapor pressure deficit, a major driver used to explain acceleration or deceleration phases in the global hydrological cycle. Current approaches to C a reconstruction rely on a robust inverse correlation between measured stomatal density in leaves (ν) of many plant taxa and C a . There are two methods that exploit this correlation: The first uses calibration curves obtained from extant species assumed to represent the fossil taxa, thereby restricting the suitable taxa to those existing today. The second is a hybrid eco-hydrological/physiological approach that determines C a with the aid of systems of equations based on quasi-instantaneous leaf-gas exchange theories and fossil stomatal data collected along with other measured leaf anatomical traits and parameters. In this contribution, a reduced order model (ROM) is proposed that derives C a from a single equation incorporating the aforementioned stomatal data, basic climate (e.g. temperature), estimated biochemical parameters of assimilation and isotope data. The usage of the ROM is then illustrated by applying it to isotopic and anatomical measurements from three extant species. The ROM derivation is based on a balance between the biochemical demand and atmospheric supply of CO 2 that leads to an explicit expression linking stomatal conductance to internal CO 2 concentration (C i ) and C a . The resulting expression of stomatal conductance from the carbon economy of the leaf is then equated to another expression derived from water vapor gas diffusion that includes anatomical traits. When combined with isotopic measurements for long-term C i /C a , C a can be analytically determined and is interpreted as the time-averaged C a that existed over the life-span of the leaf. Key advantages of the proposed ROM are: 1) the usage of isotopic data provides constraints on the reconstructed atmospheric CO 2 concentration from ν, 2) the analytical form of this approach permits direct links between parameter uncertainties and reconstructed C a , and 3) the time-scale mismatch between the application of instantaneous leaf-gas exchange expressions constrained with longer-term isotopic data is reconciled through averaging rules and sensitivity analysis. The latter point was rarely considered in prior reconstruction studies that combined models of leaf-gas exchange and isotopic data to reconstruct C a from ν. The proposed ROM is not without its limitations given the need to a priori assume a parameter related to the control on photosynthetic rate. Our work here further explores immanent constraints for the aforementioned photosynthetic parameter.

54 ENVIRONMENTAL SCIENCES↗