Search NASA⌕ Search

SEARCH · Search NASA

Results for “proxy”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Hydraulic characterization of a fault zone from fracture distribution

A quantitative assessment of how faults control the migration of geofluids is critical in many areas of geosciences. Here we integrated geological fieldwork, quantitative analysis of the fractures distribution and numerical modeling to build a geometrical representation of a fault zone and to characterize its hydraulic properties. Our target is a fault located in the Majella Mountain (Italy). We collected 21 scan lines across the fault profile in order to characterize its architecture. The numerical modeling of the fracture network of the damage zones and their hydraulic parameters was performed using both commercial (Move ® ) and open source software (dfnWorks and PFLOTRAN). Move ® was used to build a representative model of the fault zone using fracture spacing as a proxy, and to model the hydraulic parameters of the different fault domains. dfnWorks and PFLOTRAN were employed to infer the hydraulic parameters of the damage zones of the fault and then upscale these properties to an equivalent continuum domain, suitable for fluid flow simulations through the whole fault zone. Our findings show how even in a relatively small area it is possible to describe changes in terms of hydraulic properties of a fault zone and to build models capable to represent these variations.

58 GEOSCIENCES↗
DOE OSTI2023

Aperture size distribution, length, and preferential location of bed-parallel veins in shale

Bed-parallel, calcite-filled veins (BPVs) are common in shale formations, and although they have been widely described in other studies, little is known about their population aperture size distribution. To address this knowledge gap, we analyzed BPV sizes in outcrops and cores from the Vaca Muerta Formation, Neuquén Basin, Argentina; in two cores from the Marcellus Formation, Appalachian Basin, northeast Pennsylvania; and in one core from the Wolfcamp Shale, Delaware Basin, West Texas. Nine out of ten aperture size populations follow a negative exponential distribution, with one following a weak power law. Bed-parallel vein size distribution and intensity vary among formations and within the same shale. We define three groups of distributions: (1) Vaca Muerta outcrops, with the highest BPV intensity and the largest BPVs (cumulative frequency of 4.9 BPVs per meter [BPVs/m] for apertures 0.265 mm to 8.7 cm); (2) Vaca Muerta cores with a similar BPV intensity overall but with no apertures wider than 1.2 cm; and (3) Vaca Muerta, Wolfcamp, and Marcellus cores with the fewest BPVs (cumulative frequency up to 0.63 BPVs/m) and very few wider than 1 cm. Aperture and length in two outcrop data sets are weakly positively correlated and follow power laws with exponents of 0.44 and 0.49. Mechanical interfaces at boundaries between different lithologies exert a strong control on BPV location, with 65–75% of observed interfaces having BPVs along them. Only 25–30% of the BPVs occur at observed material interfaces, however, and unless subtle, unobserved mechanical layering is present, other factors must also control location. BPV intensity and organic richness (TOC) from Vaca Muerta well logs are correlated in some instances but not in others, indicating TOC is not always a good proxy for BPV location or intensity. Furthermore, these findings provide useful information for modeling of hydraulic fracture treatments where BPVs may influence development of the stimulated fracture network, for example by limiting height growth.

58 GEOSCIENCES↗
DOE OSTI2022

Identification of the optimal carbon fiber shape for cost-specific compressive performance

We report carbon fiber composites offer superior mechanical performance compared to nearly all other useful materials for the design of structures. However, for cost-driven industries, such as with the wind energy and vehicle industries, the cost of commercial carbon fiber materials is often prohibitive for their usage compared to alternatives. This paper develops an approach to optimize fiber geometries for use in carbon fiber reinforced polymers to increase the compressive strength per unit cost. Compressive strength is a composite property that depends on the fiber, matrix, and interface, and an exact analytic expression does not exist that can accurately represent these complicated relationships. The approach taken instead is to use a weighted summation between the fiber cross-sectional area moment of inertia and perimeter as a proxy for compressive strength, with different weightings explored within the paper. Analyses are performed to identify optimal fiber geometries that increase the cost-specific compressive strength based on various assumptions and desired fiber volume fraction. Robust optimal shapes are identified which outperform circular fibers due to increases in area moment of inertia and perimeter, as well as decreases in carbon fiber processing costs.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Recrystallization driven softening and heating rate dependencies of FeCrAl nuclear fuel cladding during accident transients

A refined understanding of FeCrAl cladding behavior during rapid transients is critical for its potential deployment in light-water reactors. Current assessments focus on transient burst testing metrics such as balloon geometry, burst temperature, and hoop stress, often used as proxies for simpler conventional tensile properties. However, directly correlating isothermal tensile and creep data with accident transient scenarios remains a challenge, although it is essential for high fidelity model development. Recent modeling based on tensile tests up to 800 °C, conducted with both immediate loading and a 10-minute soak, showed that immediate loading better predicts experimental burst temperatures, indicating a thermal softening effect. Building upon this observation, the current study connects transient performance, microstructural evolution, and high-temperature tensile properties by leveraging results from C26M claddings burst tests performed at heating rates of 1–50 °C/s and hoop stresses from 25 to 100 MPa. At 25 MPa, rupture temperatures varied by only 6 °C, but at 100 MPa, the difference reached 116 °C, with faster heating yielding higher burst temperatures. Microstructural analysis identified recrystallization as the primary cause of heating rate-dependent softening, eliminating prior cold-working. In-situ thermomechanical data linked ballooning onset to localized instabilities, similar to ultimate tensile strength behavior in conventional tensile tests. High heating rates correlated with immediate loading tensile data, while lower rates matched soaked data. Furthermore, by linking burst performance to microstructural evolution and tensile properties, this work provides a foundation for more accurate modeling of FeCrAl claddings and potentially other Fe-based materials under accident conditions.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Revisiting Bayesian constraints on the transport coefficients of QCD

Multistage models based on relativistic viscous hydrodynamics have proven successful in describing hadron measurements from relativistic nuclear collisions. These measurements are sensitive to the shear and the bulk viscosities of QCD and provide a unique opportunity to constrain these transport coefficients. Bayesian analyses can be used to obtain systematic constraints on the viscosities of QCD, through methodical model-to-data comparisons. In this manuscript, we discuss recent developments in Bayesian analyses of heavy ion collision data. Furthermore, we highlight the essential role of closure tests in validating a Bayesian analysis before comparison with measurements. We discuss the role of the emulator that is used as proxy for the multistage theoretical model. We use an ongoing Bayesian analysis of soft hadron measurements by the JETSCAPE Collaboration as context for the discussion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Plant wax and carbon isotope response to heat and drought in the conifer Juniperus monosperma

Plant waxes, including n-alkanes, are commonly used for a wide range of paleo-applications. Several common traits of n-alkanes that are used as paleo-proxies include chain length distribution and average chain length (ACL), as well as plant wax carbon and hydrogen isotopic compositions. The effect of climate on plant wax traits has been the subject of many studies, but a common challenge with modern calibrations is disentangling the effects of species (genetic), temperature, and precipitation from one another. Here, in this study, we explore the effect of temperature and drought, independently and combined, on plant wax composition of the species Juniper monosperma in a large ecosystem-scale field manipulation experiment. We find that n-alkane concentrations significantly increase with temperature, but other parameters (including ACL) are not affected. These results support physiological studies that identify n-alkanes as an important barrier to water loss within the plant cuticle. Combined with prior studies, it appears that changes in ACL within sediments are likely controlled by changes in species composition rather than directly by changes in climate. We find little variation in the carbon isotopic composition (δ 13 C) of n-alkanes across the treatments whereas bulk leaf δ 13 C values are higher in the heat and drought treatment. Because leaf δ 13 C values represent a weighted C assimilation signal, these values reflect differences in leaf gas exchange among treatments, whereas the n-alkanes are synthesized when water availability is higher and differences among treatments are not significant enough to influence their values. These results have important implications for using n-alkane traits, including ACL and δ 13 C values, for paleoenvironmental reconstructions.

Compound-specific isotope analysis↗
DOE OSTI2022

Ocean circulation and climate variability in the northern South China Sea during the Greek Minimum derived from coral Δ 14 C and Sr/Ca records

Ocean circulation transports heat, salt and nutrients, and has profound impacts on the marine environment and climate change. However, the seasonal to centennial variations of the paleocirculation of the South China Sea (SCS) and their driving forces are still unclear due to limited records. Here we reconstructed a high-resolution Δ 14 C record of a coral from Sanya to investigate the seasonal variations of the SCS ocean circulation around 2300 cal yr BP during a grand solar minimum (GSM) period. The seasonal Δ 14 C variability shows the influence of coastal upwelling caused by the East Asian Summer Monsoon (EASM) in summer and the effect of water intrusion from the western Pacific Ocean driven by the Kuroshio Intrusion (KI) in winter. Here we compiled the marine radiocarbon reservoir correction (ΔR) records since 2500 cal yr BP, which indicate a gradual decrease of upwelling significantly correlated with the EASM. Comparisons with other climate records suggest that both the EASM and KI may regulate the ocean circulation variability on centennial time scales. As a good index of the sea surface temperature (SST), our Sr/Ca record and spectral analysis results show a low temperature and a low frequency of El Niño-Southern Oscillation (ENSO) events around 2300 cal yr BP. Moreover, the composite coral Sr/Ca-SST and ENSO index records since 2500 cal yr BP show relatively low SST and weak ENSO during GSMs, supporting the modulation of tropical SST and ENSO by solar activity. Our study has provided high-resolution proxy data and revealed the driving forces of ocean circulation and climate change in the SCS at multiple time scales, which should be considered in further modeling work.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Optimization with the OpenACC-to-FPGA framework on the Arria 10 and Stratix 10 FPGAs

The reconfigurable computing paradigm with field programmable gate arrays (FPGAs) has received renewed interest in the high-performance computing field due to FPGAs’ unique combination of performance and energy efficiency. However, difficulties in programming and optimizing FPGAs have prevented them from being widely accepted as general-purpose computing devices. In accelerator-based heterogeneous computing, portability across diverse heterogeneous devices is also an important issue, but the unique architectural features in FPGAs make this difficult to achieve. To address these issues, a directive-based, high-level FPGA programming and optimization framework was previously developed. In this work, developed optimizations were combined holistically using the directive-based approach to show that each individual benchmark requires a unique set of optimizations to maximize performance. We perform this exploration on Intel Arria 10 and Stratix 10 FPGAs. We also explored the relationships between performance, resource usages, and compilation times, and investigated implications for performance portability. Finally, we present an initial evaluation of a real-world proxy application, LULESH.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Characterization of nuclear breakup as a function of hard-scattering kinematics using dijets measured by ATLAS in p+Pb collisions

This Letter analyzes the sensitivity of event geometry estimators to the initial-state kinematics of hard scattering in proton–lead collisions. This analysis uses dijets as a proxy for the parton–parton scattering configuration, correlating it with event geometry estimators, namely the energy deposited in the Zero-Degree Calorimeter and the transverse energy recorded in the Forward Calorimeter in the Pb-going direction. The analysis uses data recorded by the ATLAS detector at the Large Hadron Collider with a nucleon–nucleon center-of-mass energy of 8.16 TeV, corresponding to an integrated luminosity of 56 nb − 1 . The jets are measured within the pseudorapidity interval −2.8 < η < 4.5, where positive η values correspond to the direction of the proton beam. Results are presented as a function of the Bjorken-x of the parton originating from the proton, xp . Both event geometry estimators are found to be dependent on xp , with the energy deposited in the Zero-Degree Calorimeter about six times less sensitive to xp compared with the transverse energy deposited in the Forward Calorimeter.

Aad, G↗
DOE OSTI2020

Late Holocene history of Tokotan Lake (Kuril Archipelago, Russian Far East): The use of lacustrine records for paleoclimatic reconstructions from geologically dynamic settings

A high-resolution, multiproxy analysis of a core from Tokotan Lake, Urup Island, illustrates the complex paleoenvironmental signals that can be preserved in sedimentary deposits from geologically active areas, such as the Kuril Archipelago. Diatom and geochemical analyses of the Tokotan core suggest a shallowing and lessening of the lake area over the last ~3450 cal BP. Without other proxies, the record could be interpreted as showing a general decrease in precipitation over the Late Holocene. However, palynological data indicate no change in vegetation, thereby implying no significant climatic variation. Neither of these paleoclimatic interpretations agrees with Late Holocene reconstructions from the Kuril Archipelago or the Russian mainland. Paleomagnetic data imply variable sediment sources during the lake's history, indicating alterations within the Tokotan catchment (e.g., shifts in gradients or sediment loads in the inflowing streams). Such changes would be particularly important in the Shabalina River, which today has formed a partial dam across the lake's outlet. Here, the results of the Tokotan study indicate that other factors (e.g., volcanism, seismic activity, sea-level changes) rather than climate were more probable determinants of the Late Holocene lake record. Given the geologic dynamism of the region, it is perhaps not surprising that paleo-records from the Kuril Islands display an inconsistency in the numbers and/or duration of proposed Holocene climatic oscillations, underscoring the importance of multiproxy investigations in areas that are or have been geologically dynamic.

58 GEOSCIENCES↗
DOE OSTI2021

Climate reconstructions for the Last Glacial Maximum from a simple cirque glacier in Fiordland, New Zealand

Mountain glacier records offer important constraints on the timing and magnitude of climate variability during the last glacial cycle. Existing moraine chronologies from the central Southern Alps indicate maximum ice extent was achieved during marine isotope stages 3–4, followed by repeated advances of similar, but gradually declining extent during marine isotope stage 2, until onset of the glacial termination. Questions remain over the precise role of climate in driving these changes, as most existing moraine chronologies come from large, complex former valley glacier systems, where non-climatic influences such as changing bed topography and proglacial lake formation may have influenced glacier length changes. In this work, we address this problem via a new cosmogenic 10Be chronology and equilibrium line altitude reconstruction from a cirque glacier situated in Fiordland, New Zealand. Our chronology shows moraine deposition at 32 ± 11 ka, 18.7 ± 0.2 ka, 18.1 ± 0.1 ka, and c. 17.2 ± 0.3 ka. The simple geometry of the former glacier supports the role of climate in driving a net decline of regional ice volume during marine isotope stages 3–2. Close spacing and good preservation of the 19–17 ka moraines permits 2D glacier reconstruction which suggests the equilibrium line altitude remained depressed by c.1130 m (equivalent to 5.8 ± 0.6 °C colder than present) during this interval. Onset of warming after 17.2 ± 0.2 ka is consistent with climate proxy evidence for a sustained southward shift in the southern westerly winds, which may have promoted deglaciation via shifting ocean currents and promoting increases in atmospheric carbon dioxide.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Cosmogenic 10Be constraints on deglacial snowline rise in the Southern Alps, New Zealand

We report geochronological dating of glacial landforms, such as terminal and lateral moraines, is useful for determining the extent and timing of past glaciation and for reconstructing the magnitude and rate of past climate changes. In the Southern Alps of New Zealand, well-dated glacial geomorphological records constrain the last glacial cycle across much of the Waitaki River basin (e.g. Ōhau, Pukaki, Tekapo) but its southern sector such as the Ahuriri River valley remains comparatively unconstrained. Recently, there has been debate on the scale and rapidity of mountain glacier retreat during the last glacial termination, particularly the 20–17 ka period in New Zealand. Missing from this debate is well-constrained equilibrium-line altitude (ELA) and associated temperature reconstructions, particularly over the period around 17 ka, which can help us to develop a more complete picture of how past temperature changes drove glacier retreat. Here we report the first glacial chronology dataset from the Last Glacial Maximum (LGM) and subsequent deglaciation from the Ahuriri River valley, Southern Alps, New Zealand (44°23'54"S, 169°39'48"E) based on 38 beryllium-10 ( 10 Be) surface-exposure ages from terminal moraine systems and glaciated bedrock situated at the lower and middle sections of the valley. Our results show that the former Ahuriri Glacier reached its maximum extent at 19.8 ± 0.3 ka, which coincides with the global Last Glacial Maximum. By 16.7 ± 0.3 ka, the glacier had retreat ~18 km up-valley suggesting at least ~43% glacier-length loss relative to its full LGM extent. This deglaciation was accompanied by the formation of a shallow proglacial lake. Using the accumulation area ratio (AAR) method, we estimate that the ELA was lower than present by ~880 m (~1120 m a.s.l.) at 19.8 ± 0.3 ka, and ~770 m lower (~1230 m a.s.l.) at 16.7 ± 0.3 ka. Applying an estimate for temperature lapse rate, this ELA anomaly implies that local air temperature was 5 ± 1 °C colder than present (1981–2010) at 19.8 ± 0.3 ka, while it was 4.4 ± 0.9 °C colder at 16.7 ± 0.3 ka, assuming no change in precipitation. The substantial glacier retreat in response to a relatively small accompanying increases in ELA (110 m) and temperature (0.6 °C) may have been a result of the high glacier-length sensitivity of this glacier system due to its low gradient of former ice surface. Our low warming estimate differs markedly from other deglaciation studies, specifically from Rakaia River valley, which reports a much larger temperature increase at the onset of the last deglaciation. This precisely-dated moraine record along with reconstructed ELA as proxies for atmospheric conditions, provides new insight into post LGM glacier behaviour and climate conditions in New Zealand.

58 GEOSCIENCES↗
DOE OSTI2025

Paleoclimatic implications of glacial fluctuations in the Sierra Nevada del Cocuy, northern Andes, Colombia, during the Lateglacial and Holocene

The reconstruction of former mountain glaciers from geomorphic mapping and cosmogenic-nuclide surface-exposure dating provides a unique opportunity to infer patterns of past terrestrial climate variability. Tropical mountain glaciers are particularly valuable as there are comparatively few terrestrial climate proxies at equatorial latitudes relative to higher latitudes. As the single largest climate zone on Earth, the tropics play an outsized role in mediating global climate via the ocean-atmosphere transfer of latent heat and water vapor. Nonetheless, there remains a persistent gap in our understanding of how the tropics influenced – or were influenced by – the high-magnitude climate shifts of the Late Pleistocene, and whether this high-energy region simply responded to extratropical forcing or was itself a driver of global climatic change. To help address this knowledge gap, we analyzed geologic evidence for past glacial fluctuations in three adjacent valleys in the Sierra Nevada del Cocuy, the highest subrange of the Eastern Cordillera in the Colombian Andes, to provide a terrestrial record of atmospheric temperature during the latter part of Termination 1. Coupled with geomorphic mapping and paleo-snowline reconstructions, our beryllium-10 glacial chronology indicates that glaciers in the humid inner tropics underwent pronounced growth and gradual decay during the Antarctic Cold Reversal (14.5–12.8 ka) and Younger Dryas (12.8–11.7 ka) periods, respectively, following a trend that, according to directly dated moraine records from throughout both polar hemispheres, appears to have been global. While the specific mechanism(s) behind this large-scale behavior remains to be corroborated, we revisit the hypothesis that ocean-atmosphere heat transfer and water vapor flux are key drivers of abrupt Lateglacial temperature fluctuations. Subsequent to the Lateglacial, deglaciation of the Sierra Nevada del Cocuy accelerated during the Early Holocene, a pattern also observed in other tropical glacier records. More recently, the magnitude of snowline rise and glacier retreat over the last two centuries supports the view that modern tropospheric warming is anomalously strong at least relative to the last ∼16,000 years.

Andes↗
DOE OSTI2020

Downwind coning concept rotor for a 25 MW offshore wind turbine

The size of offshore wind turbines over the next decade is expected to continually increase due to reduced balance of station costs per MW and also the higher wind energy at increased altitudes that can lead to higher capacity factors. However, there are challenges that may limit the degree of upscaling which is possible. In this paper, a two-bladed downwind turbine system is upscaled from 13.2 MW to 25 MW, by redesigning aerodynamics, structures, and controls. In particular, three 25 MW rotors have been developed: V1 is the upscaled model, V2 is a partial redesigned model, and V3 is a fully redesigned model. Despite their radically large sizes, it is found that these 25 MW turbine rotors satisfy this limited set of design drivers at the rated condition and that larger blade lengths are possible with cone-wise load-alignment. In addition, flapwise morphing (varying the cone angle with a wind-speed schedule) is investigated in terms of minimizing mean and fluctuating root bending loads using steady inflow proxies for the maximum and damage equivalent load moments. The resulting series of 25 MW rotors, which are the largest ever designed, can be a useful baseline for additional development and assessment.

17 WIND ENERGY↗
DOE OSTI2022

Machine learning and shallow groundwater chemistry to identify geothermal prospects in the Great Basin, USA

This study discovers various geothermal prospects in the Great Basin, USA based on shallow groundwater chemical (geochemical) data. The geochemical data are expected to include hidden (latent) information that is a proxy for geothermal prospectivity. We processed the sparse geochemical data in the Great Basin at 14,341 locations including 18 attributes. Next, a non-negative matrix factorization with customized k-means clustering is applied to the geochemical data matrix that automatically finds three hidden geothermal signatures representing modestly, moderately, and highly confident geothermal prospects. The algorithm also evaluated the probability of occurrence of these types of resources through the studied region. There is a consistency between regional geothermal prospectivity as estimated by our ML methodology and the traditional play fairway analysis conducted over a portion of the study area. We also identify the dominant data attributes associated with each signature. Finally, our ML analyses allow us to reconstruct attributes from sparse into continuous over the study domain. The predicted continuous attributes can be used for future detailed geothermal explorations in the Great Basin.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Towards multi-fidelity deep learning of wind turbine wakes

We report engineering wake models that accurately predict wake in a computationally efficient manner are very important for tasks such as layout optimization and control of wind farms. In this paper, we explore an application of deep learning (DL) to learn the wake model from hierarchies of physics-based approaches ranging from analytical models to an approximate form of the Reynolds-averaged Navier-Stokes equations. We first illustrate the application of principal component analysis to obtain a lower-dimensional representation that allows a computationally tractable training and deployment of DL models. Then, the DL model is trained to learn the mapping from input parameter space to the principal components, which are then used to reconstruct the three-dimensional flow field. Additionally, we investigate a composite framework consisting of two neural networks to learn the correlation between low- and high-fidelity data with Gauss and curl models treated as proxies for low- and high-fidelity models, respectively. The prediction from both DL models matches well with the high-fidelity data with a maximum relative percentage error for the kinetic energy flux of <1%. This work opens up possibilities for data-efficient construction of surrogate models for wake prediction that can be used to study the influence of wind speed and yaw angles on wind farm power production.

17 WIND ENERGY↗
DOE OSTI2021

TROPOMI observations allow for robust exploration of the relationship between solar-induced chlorophyll fluorescence and terrestrial gross primary production

Solar-induced chlorophyll fluorescence (SIF) observed by satellites has advanced the monitoring of terrestrial photosynthesis regionally and globally. The relationship between SIF and gross primary production (GPP) at leaf, canopy, and ecosystem scales has received tremendous attention in recent years. It remains controversial whether the SIF-GPP relationship at the ecosystem scale is universal or dependent upon vegetation type. New SIF observations from the TROPOspheric Monitoring Instrument (TROPOMI) with unprecedented high spatial and temporal resolution provide a new opportunity to elucidate the SIF-GPP relationship. Here, we examine the SIF-GPP relationship for seven major vegetation types across the U.S. with TROPOMI SIF and in-situ GPP data for 83 eddy covariance flux sites. We find that TROPOMI SIF shows a strong and consistent relationship with tower-based GPP at both satellite footprint and grid-cell levels. The slope of the SIF-GPP relationship is similar among all the vegetation types except croplands, demonstrating a nearly universal (converging to ~13.5 g C m -2 d -1 /W m -2 μm -1 sr -1) rather than vegetation type-specific SIF-GPP relationship. This confirms that TROPOMI SIF can be used as a proxy for GPP across a wide variety of vegetation types and can also be used to quantify GPP by avoiding uncertainty associated with land cover maps. The C 4 crops have a much higher slope than the C 3 crops, and therefore croplands tend to have a higher slope than C 3 -dominated vegetation types (e.g., forests, shrublands, savannas). We also find that the TROPOMI SIF is well correlated with GPP under normal or wetter conditions, while their relationship becomes weaker under water stress. In conclusion, our TROPOMI-based study could improve our understanding of the SIF-GPP relationship at the ecosystem scale and advance the mapping of GPP globally with SIF observations from space.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Decoupling of greenness and gross primary productivity as aridity decreases

Ecosystem primary productivity is a key ecological process influencing many ecosystem services, including carbon storage. Thus, clarifying how primary productivity in terrestrial ecosystems responds to climatic variability can reveal key mechanisms that will drive future changes in the global carbon budget. Satellite products of canopy greenness are widely used as proxies for vegetation productivity to evaluate how ecosystems respond to climate variability. However, to what degree inter-annual variations in productivity are consistent with greenness and how this relationship varies spatially remains unclear. Here we investigated the strength of the coupling between inter-annual variations in leaf area index (LAI, a measure of greenness) and ecosystem gross primary productivity (GPP) derived from eddy covariance towers, i.e., the r2 of the LAI-GPP relationship. Overall, inter-annual GPP and LAI were highly coupled (i.e., high r2) in arid grasslands, but were fully decoupled in mesic evergreen broadleaf forests, indicating that this relationship varies strongly along aridity gradients. A possible mechanism of the spatial variation in the LAI-GPP relationship is that the tradeoff between ecosystem structure (LAI) and physiology (photosynthesis per unit leaf area) becomes stronger in more humid climates. Land models overestimated the r2 of LAI-GPP correlation for most ecosystem types and failed to capture the spatial pattern along aridity gradients. We conclude that relying on greenness products for evaluating inter-annual changes in vegetation productivity may bias assessments, especially in tropical rainforest ecosystems. Our findings may also reconcile observed disparities between responses in greenness and GPP during drought in Amazon forests.

54 ENVIRONMENTAL SCIENCES↗