Search NASA⌕ Search

SEARCH · Search NASA

Results for “Ensemble Impacts”

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 361 records · Page 20

Computational investigation of hysteresis and phase equilibria of n-alkanes in a metal-organic framework with both micropores and mesopores

Abstract Adsorption hysteresis is a phenomenon related to phase transitions that can impact applications such as gas storage and separations in porous materials. Computational approaches can greatly facilitate the understanding of phase transitions and phase equilibria in porous materials. In this work, adsorption isotherms for methane, ethane, propane, and n-hexane were calculated from atomistic grand canonical Monte Carlo (GCMC) simulations in a metal-organic framework having both micropores and mesopores to better understand hysteresis and phase equilibria between connected pores of different size and the external bulk fluid. At low temperatures, the calculated isotherms exhibit sharp steps accompanied by hysteresis. As a complementary simulation method, canonical (NVT) ensemble simulations with Widom test particle insertions are demonstrated to provide additional information about these systems. The NVT+Widom simulations provide the full van der Waals loop associated with the sharp steps and hysteresis, including the locations of the spinodal points and points within the metastable and unstable regions that are inaccessible to GCMC simulations. The simulations provide molecular-level insight into pore filling and equilibria between high- and low-density states within individual pores. The effect of framework flexibility on adsorption hysteresis is also investigated for methane in IRMOF-1.

36 MATERIALS SCIENCE↗

Data-driven closure modeling for hypersonic turbulent flows

The Reynolds-averaged Navier–Stokes (RANS) equations remain a workhorse technology for simulating compressible fluid flows of practical interest. Due to model-form errors, however, RANS models can yield erroneous predictions that preclude their use on mission-critical problems. This report summarizes work performed from FY22-FY24 focused on improving RANS models for hypersonic flows using data-driven modeling and scientific machine learning. In this work we: 1. Investigate the current capabilities of RANS models in Sandia’s parallel aerodynamics and re-entry code (SPARC) for hypersonic flows with a focus on shock boundary layer interactions (SBLIs), 2. Assess several established corrections that exist in the literature aimed at improving predictions for SBLIs, 3. Develop improved models for the Reynolds stress tensor using tensor-basis neural networks, 4. Develop a neural-network-based variable turbulent Prandtl number model to reduce errors in wall heating in SBLIs. 5. Begin future investigations including employing the LIFE framework to improve wall heating predictions in SBLIs as well as the ensemble Kalman filter. We find that current RANS models in SPARC are deficient for complex SBLI flows. In particular, no current model jointly predicts wall heat flux, wall shear stress, and wall pressure with reasonable accuracy. Existing corrections help, but do not alleviate this issue altogether. The development of improved models for the Reynolds stress tensor via tensor-basis neural networks results in more predictive RANS models across a suite of low-speed and high-speed cases. For hypersonic boundary layers, the inclusion of the wall-normal Reynolds stress via TBNNs has an appreciable impact on the wall-normal momentum balance and wall quantities. However, we find that improvements to the Reynolds stress tensor do not address the over-prediction in wall heat flux in SBLIs. We find that a neural-network-based variable turbulent Prandtl number model systematically and substantially improves wall heating predictions for a range of SBLI cases.

97 MATHEMATICS AND COMPUTING↗

A multi‐variable framework for selecting WRF physics configurations at convection‐permitting scales: An Amazon wet‐season case study

Tropical convection over rainforests modulates atmospheric circulation and the energy and hydrological cycles across multiple scales. However, the scarcity of observations still limits our understanding of these processes. Although numerical models are utilized to investigate atmospheric physical processes, their performance depends on the choice of parameterizations and grid resolution. Here, in this work, we introduce and apply a multi‐variable, multi‐physics framework to evaluate and rank the Advanced Research Weather Research and Forecasting (WRF–ARW) configurations at 1‐km resolution over the central Amazon during the wet season. A 48‐member ensemble combines three land‐surface models (LSMs), four planetary boundary layer (PBL), and four microphysics (MP) schemes. Model performance for seven convection‐related near‐surface and boundary‐layer variables is assessed using Taylor diagrams and the Taylor skill score (TSS), analysis of variance (ANOVA)‐based sensitivity metrics, and non‐parametric rank tests. We then construct a combined, weighted TSS to identify configurations that are comparatively robust across variables. Results showed that most configurations reproduce near‐surface temperature, sensible and latent heat fluxes, and boundary‐layer height reasonably well, whereas humidity and rainfall remain challenging. LSM choice has the strongest impact on the surface fluxes and a secondary influence on near‐surface temperature and humidity, PBL schemes dominate boundary‐layer height, and MP schemes exert the largest control on rainfall. No single configuration is optimal for all variables, but the combination of Noah (LSM), Yonsei University (PBL), and Morrison (MP) emerges as the most robust configuration for this case, with the WRF single‐moment six‐class scheme (WSM6) providing a competitive, computationally cheaper MP alternative. The framework is general and can be applied to other regions, seasons, and convective regimes.

Amazon rainforest↗

Matilda v1.0: An R package for probabilistic climate projections using a reduced complexity climate model

A primary advantage to using reduced complexity climate models (RCMs) has been their ability to quickly conduct probabilistic climate projections, a key component of uncertainty quantification in many impact studies and multisector systems. Providing frameworks for such analyses has been a target of several RCMs used in studies of the future co-evolution of the human and Earth systems. In this paper, we present Matilda, an open-science R software package that facilitates probabilistic climate projection analysis, implemented here using the Hector simple climate model in a seamless and easily applied framework. The primary goal of Matilda is to provide the user with a turn-key method to build parameter sets from literature-based prior distributions, run Hector iteratively to produce perturbed parameter ensembles (PPEs), weight ensembles for realism against observed historical climate data, and compute probabilistic projections for different climate variables. This workflow gives the user the ability to explore viable parameter space and propagate uncertainty to model ensembles with just a few lines of code. The package provides significant freedom to select different scoring criteria and algorithms to weight ensemble members, as well as the flexibility to implement custom criteria. Additionally, the architecture of the package simplifies the process of building and analyzing PPEs without requiring significant programming expertise, to accommodate diverse use cases. We present a case study that provides illustrative results of a probabilistic analysis of mean global surface temperature as an example of the software application.

54 ENVIRONMENTAL SCIENCES↗

Assessing the representation of the Australian carbon cycle in global vegetation models

Australia plays an important role in the global terrestrial carbon cycle on inter-annual timescales. While the Australian continent is included in global assessments of the carbon cycle such as the global carbon budget, the performance of dynamic global vegetation models (DGVMs) over Australia has rarely been evaluated. We assessed simulations of net biome production (NBP) and the carbon stored in vegetation between 1901 to 2018 from 13 DGVMs (TRENDY v8 ensemble). We focused our analysis on Australia's short-term (inter-annual) and long-term (decadal to centennial) terrestrial carbon dynamics. The TRENDY models simulated differing magnitudes of NBP on inter-annual timescales, and these differences resulted in significant differences in long-term vegetation carbon accumulation (–4.7 to 9.5 Pg C). We compared the TRENDY ensemble to several satellite-derived datasets and showed that the spread in the models' simulated carbon storage resulted from varying changes in carbon residence time rather than differences in net carbon uptake. Differences in simulated long-term accumulated NBP between models were mostly due to model responses to land-use change. The DGVMs also simulated different sensitivities to atmospheric carbon dioxide (CO 2 ) concentration, although notably, the models with nutrient cycles did not simulate the smallest NBP response to CO 2 . Our results suggest that a change in the climate forcing did not have a large impact on the carbon cycle on long timescales. However, the inter-annual variability in precipitation drives the year-to-year variability in NBP. We analysed the impact of key modes of climate variability, including the El Niño–Southern Oscillation (ENSO) and the Indian Ocean Dipole (IOD), on NBP. While the DGVMs agreed on sign of the response of NBP to El Niño and La Niña and to positive and negative IOD events, the magnitude of inter-annual variability in NBP differed strongly between models. In addition, we find that differences in the timing of simulated phenology and fire dynamics are associated with differences in simulated or prescribed vegetation cover and process representation. We further find model disagreement in simulated vegetation carbon, phenology, and apparent carbon residence time, indicating that the models have different types and coverage of vegetation across Australia (whether prescribed or emergent). Our study highlights the need to evaluate parameter assumptions and the key processes that drive vegetation dynamics, such as phenology, mortality, and fire, in an Australian context to reduce uncertainty across models.

54 ENVIRONMENTAL SCIENCES↗

CMIP7 data request: impacts and adaptation priorities and opportunities

The Coupled Model Intercomparison Project Phase 7 (CMIP7) undertook an extensive process to gather community input and refine data requests related to impacts and adaptation applications of Earth System Model (ESM) outputs. The Impacts and Adaptation (I&A) Data Request Team worked with CMIP7 leadership to distribute an open solicitation across many communities that use climate model outputs requesting inputs for new and existing variables, the most applicable temporal characteristics, and groupings of variables that together allow for specific application opportunities. This input was then collated and translated into CMIP7 standard templates for inclusion in the broader data request, leading to 13 I&A data request opportunities, 60 variable groups and 539 unique variables sought by vulnerability, impacts, adaptation, and climate services user communities. Here, we describe these opportunities and variable groups, as well as new insights into how ESM groups can prioritize outputs that set off a chain of further analyses, ultimately informing decisions impacting society and natural systems. These include an emphasis on high-resolution outputs to allow further modeling of climate impacts at regional and local scales, improved representation of extreme weather events, enhanced accuracy of downscaling and bias-adjustment techniques, and support for more detailed assessments for decision-making in adaptation and mitigation strategies. There is also broad interest in more extensive provisioning of two-dimensional variables at the Earth's surface, prioritizing experiments that enhance our understanding of both the recent past and future scenarios, and providing outputs that allow further downscaling and bias adjustment. We emphasize that variable groups are the fundamental level at which to engage with the I&A data request, matching the scale of input and the way output provision enables specific I&A applications. Given resource constraints, we applaud CMIP7 efforts to foster strong engagement and communication between ESM groups and the I&A team to build consensus around prudent compromises in priority variables, temporal resolutions, simulation experiments, time subsets, and ensemble members.

Ruane, Alex C. [NASA Goddard Inst. for Space Studi↗

Estimators and Fusers for Fiber Delay Estimation Using Environmental Measurements

The properties of deployed network fiber are affected by environmental factors due to their exposure to the elements. Particularly for quantum networks, the resultant delay variations may have significant impacts due to the extreme sensitivity of synchronization, coincidence counting, and other critical operations. In this paper, the delays of 15 km aerial-inground fiber connections are measured, and effects due to temperature, humidity and wind speed are analyzed over multiple periods spanning four seasons of a year. Machine learning methods are first utilized to reveal surprisingly pronounced effects of humidity on the delay, in addition to the expected temperature and its seasonal variations. Estimator and fusion methods are developed to estimate the delay using temperature, humidity and wind speed measurements, by utilizing smooth Gaussian Process Regression (GPR) and nonsmooth Ensemble of Trees (EOT) methods. Measurements from winter and summer periods are temporally fused using twelve different methods, and eight methods provide estimates for the delay throughout the year with median test errors under 1.28%. The results reveal distinct temperature-humidity trends across the seasons, and the ability of estimator and temporal fusion methods to exploit them for estimating the delay. These results constitute a case study of machine learning analytical results, wherein generalization equations explain the performance of various estimator and fuser methods.

Rao, Nageswara [ORNL] (ORCID:0000000234085941)↗

Cross-Layered Distributed Data-Driven Framework for Enhanced Smart Grid Cyber-Physical Security

Smart Grid (SG) research and development has drawn much attention from academia, industry and government due to the great impact it will have on society, economics and the environment. Securing the SG is a considerably significant challenge due the increased dependency on communication networks to assist in physical process control, exposing them to various cyber-threats. In addition to attacks that change measurement values using False Data Injection (FDI) techniques, attacks on the communication network may disrupt the power system's real-time operation by intercepting messages, or by flooding the communication channels with unnecessary data. Addressing these attacks requires a cross-layer approach. In this paper a cross-layered strategy is presented, called Cross-Layer Ensemble CorrDet with Adaptive Statistics(CECD-AS), which integrates the detection of faulty SG measurement data as well as inconsistent network inter-arrival times and transmission delays for more reliable and accurate anomaly detection and attack interpretation. Numerical results show that CECD-AS can detect multiple False Data Injections, Denial of Service (DoS) and Man In The Middle (MITM) attacks with a high F1-score compared to current approaches that only use SG measurement data for detection such as the traditional physics-based State Estimation, Ensemble CorrDet with Adaptive Statistics strategy and other machine learning classification-based detection schemes.

cyber-physical security↗

An orthologous gene coevolution network provides insight into eukaryotic cellular and genomic structure and function

The evolutionary rates of functionally related genes often covary. We present a gene coevolution network inferred from examining nearly 3 million orthologous gene pairs from 332 budding yeast species spanning ~400 million years of evolution. Network modules provide insight into cellular and genomic structure and function. Examination of the phenotypic impact of network perturbation using deletion mutant data from the baker’s yeast Saccharomyces cerevisiae, which were obtained from previously published studies, suggests that fitness in diverse environments is affected by orthologous gene neighborhood and connectivity. Mapping the network onto the chromosomes of S. cerevisiae and Candida albicans revealed that coevolving orthologous genes are not physically clustered in either species; rather, they are often located on different chromosomes or far apart on the same chromosome. The coevolution network captures the hierarchy of cellular structure and function, provides a roadmap for genotype-to-phenotype discovery, and portrays the genome as a linked ensemble of genes.

59 BASIC BIOLOGICAL SCIENCES↗

HOD-dependent systematics for luminous red galaxies in the DESI 2024 BAO analysis

In this paper, we present the estimation of systematics related to the halo occupation distribution (HOD) modeling in the baryon acoustic oscillations (BAO) distance measurement of the Dark Energy Spectroscopic Instrument (DESI) 2024 analysis. This paper focuses on the study of HOD systematics for luminous red galaxies (LRG). We consider three different HOD models for LRGs, including the base 5-parameter vanilla model and two extensions to it, that we refer to as baseline and extended models. The baseline model is described by the 5 vanilla HOD parameters, an incompleteness factor and a velocity bias parameter, whereas the extended one also includes a galaxy assembly bias and a satellite profile parameter. We utilize the 25 dark matter simulations available in the AbacusSummit simulation suite at z=0.8 and generate mock catalogs for our different HOD models. To test the impact of the HOD modeling in the position of the BAO peak, we run BAO fits for all these sets of simulations and compare the best-fit BAO-scaling parameters α iso and α AP between every pair of HOD models. Furthermore, we do this for both Fourier and configuration spaces independently, using post-reconstruction measurements. We find a 3.3σ detection of HOD systematic for α AP in configuration space with an amplitude of 0.19%. For the other cases, we did not find a 3σ detection, and we decided to compute a conservative estimation of the systematic using the ensemble of shifts between all pairs of HOD models. By doing this, we quote a systematic with an amplitude of 0.07% in α iso for both Fourier and configuration spaces; and of 0.09% in α AP for Fourier space.

79 ASTRONOMY AND ASTROPHYSICS↗

ETHOS: An automated framework to generate multi-fidelity constitutive data tables and propagate uncertainties to hydrodynamic simulations

Accurate constitutive data, such as equations of state and plasma transport coefficients, are necessary for reliable hydrodynamic simulations of plasma systems such as fusion targets, planets, and stars. Here, we develop a framework for automatically generating transport-coefficient tables using a parameterized model that incorporates data from both high-fidelity sources (e.g., density functional theory calculations and reference experiments) and lower-fidelity sources (e.g., average-atom and analytic models). The framework incorporates uncertainties from these multi-fidelity sources, generating ensembles of optimally diverse tables that are suitable for uncertainty quantification of hydrodynamic simulations. We illustrate the utility of the framework with magnetohydrodynamic simulations of magnetically launched flyer plates, which are used to measure material properties in pulsed-power experiments. We explore how changes in the uncertainties assigned to the multi-fidelity data sources propagate to changes in simulation outputs and find that our simulations are most sensitive to uncertainties near the melting transition. The presented framework enables computationally efficient uncertainty quantification that readily incorporates new high-fidelity measurements or calculations and identifies plasma regimes where additional data will have high impact.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Noise-induced transitions of the Atlantic Meridional Overturning Circulation in CMIP5 models

By studying transition probabilities of the Atlantic Meridional Overturning Circulation (AMOC) in an ensemble of CMIP5 climate models, we revisit one of the stability indicators of the AMOC, i.e. the freshwater transport carried by the AMOC at the southern boundary of the Atlantic basin. A correction to this indicator, based on the transition probabilities, is suggested to measure whether an AMOC state is in a multiple equilibrium regime or not. As a consequence, the AMOC of all CMIP5 models considered is in a multiple equilibrium regime and hence, in principle, a collapsed AMOC state should exist in each of these models. The results further demonstrate the dependence of the Atlantic surface freshwater flux on the AMOC and the impact of extreme events in the AMOC on temperatures in the North Atlantic region.

54 ENVIRONMENTAL SCIENCES↗

Uncertainty-quantified phenomenological optical potentials for single-nucleon scattering

Optical-model potentials (OMPs) continue to play a key role in nuclear reaction calculations. However, the uncertainty of phenomenological OMPs in widespread use—inherent to any parametric model trained on data—has not been fully characterized, and its impact on downstream users of OMPs remains unclear. Here we assign well-calibrated uncertainties for two representative global OMPs, those of Koning-Delaroche and Chapel Hill '89, using Markov-chain Monte Carlo for parameter inference. By comparing the canonical versions of these OMPs against the experimental data originally used to constrain them, we show how a lack of outlier rejection and a systematic underestimation of experimental uncertainties contributes to bias of, and overconfidence in, best-fit parameter values. Our updated, uncertainty-quantified versions of these OMPs address these issues and yield complete covariance information for potential parameters. Scattering predictions generated from our ensembles show improved performance both against the original training corpora of experimental data and against a new “test” corpus comprising many of the experimental single-nucleon scattering data collected over the last twenty years. Finally, we apply our uncertainty-quantified OMPs to two case studies of application-relevant cross sections. We conclude that, for many common applications of OMPs, including OMP uncertainty should become standard practice. Furthermore, to facilitate their immediate use, digital versions of our updated OMPs and related tools for forward uncertainty propagation are included as Supplemental Material.

150 ≤ A ≤ 189↗

Evolution of the Internal Climate Modes under Future Warming

Abstract Climate modes play an important role in weather and climate variability over multiple spatial and temporal scales. This research assesses Earth system model (ESM) projections of the spatiotemporal characteristics of key internal climate modes (NAM, SAM, PNA, ENSO, PDO, and AMO) under high (SSP585) and low (SSP126) radiative forcing scenarios and contextualizes those projections using historical fidelity. Time series analyses are used to assess trends and mode phase characteristics are summarized for the historical period and for the end of the twenty-first century. Spatial patterns are compared to infer morphological changes. Shifts in the power spectra are used to examine changes in variability at subannual, interannual, and interdecadal scales. Changes in time-lagged correlations are used to capture the evolution of first-order interactions. While differences in historical skill are predominantly ESM dependent, changing mode characteristics in a warmer climate also exhibit variability between individual ensemble realizations. NAM, SAM, and ENSO tend to evolve toward increased prevalence of the positive phase up to 2100 across the multimodel ensemble while the PNA and PDO exhibit little trend but increasing phase intensity. AMO characteristics are shown to depend on the method used to remove the external signal. ESMs that show higher historical fidelity tend to show more modest changes in those modes under global nonstationarity. Changes in mode interactions are found to be highly ESM dependent but exhibit broadly similar behavior to historical relationships. These findings have implications for our understanding of internal variability and make clear that the choice of ESM, and even the ESM realization, matters for applications of climate projections. Significance Statement Internal modes of variability are important to understand due to their impact on local, regional, and global weather and climate patterns. Future climate changes will not only be affected by the variability arising from these modes, but the modes will themselves change in response to the changing climate. Spatial and temporal aspects of the modes are assessed from projections of future climate and related to how well they are captured in the historical climate. This yields some measure of confidence in the changes exhibited by the models. In most cases, when historically skillful models exhibit changes that are different from those produced by less skillful models, they tend to produce more modest changes. These results, as well as the variability between model outcomes, mean decisions on which ESM to use for projections of the future climate matter significantly.

Meteorology & Atmospheric Sciences↗

Quantifying the impact of temporal analysis of products reactor initial state uncertainties on kinetic parameters

The temporal analysis of products (TAP) reactor, a transient kinetic tool, provides users with information as the catalyst state evolves. However, the state of the art for TAP uncertainty quantification only considers the experimental noise present in the outlet flow signal. Additional sources of uncertainty, including, but not limited to, initial surface coverages, the catalyst zone location, the inert void fraction, and the gas pulse intensity and pulse delay, are not included. For this reason, a framework for quantifying all uncertainty sources present in TAP experiments is presented and applied to a carbon monoxide oxidation case study. Herein, two methods for quantifying these sources of uncertainty are introduced. The first utilizes initial state sensitivities to approximate the parameter variances, as well as to provide insights into the structural certainty of the model. The second generates parameter confidence distributions through an ensemble-based sampling algorithm. This initial state covariance matrix can ultimately be merged with the experimental noise covariance matrix, providing a unified description of the parameter uncertainties for a single TAP experiment.

36 MATERIALS SCIENCE↗

Simulation of Continental Shallow Cumulus Populations Using an Observation-Constrained Cloud-System Resolving Model

Continental shallow cumulus (ShCu) clouds observed on 30 August 2016 during the Holistic Interactions of Shallow Clouds, Aerosols, and Land-Ecosystems (HI-SCALE) field campaign are simulated by using an observation-constrained cloud-system resolving model. On this day, ShCu forms over Oklahoma and southern Kansas and some of these clouds transition to deeper, precipitating convection during the afternoon. We apply a four-dimensional ensemble-variational (4DEnVar) hybrid technique in the Community Gridpoint Statistical Interpolation (GSI) system to assimilate operational data sets and unique boundary layer measurements including a Raman lidar, radar wind profilers, radiosondes, and surface stations collected by the U.S. Department of Energy's (DOE) Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) atmospheric observatory into the Weather Research and Forecasting (WRF) model to ascertain how improved environmental conditions can influence forecasts of ShCu populations and the transition to deeper convection. Independent observations from aircraft, satellite, as well as ARM's remote sensors are used to evaluate model performance in different aspects. Several model experiments are conducted to identify the impact of data assimilation (DA) on the prediction of clouds evolution. The analyses indicate that ShCu populations are more accurately reproduced after DA in terms of cloud initiation time and cloud base height, which can be attributed to an improved representation of the ambient meteorological conditions and the convective boundary layer. Extending the assimilation to 18 UTC (local noon) also improved the simulation of shallow-to-deep transitions of convective clouds.

54 ENVIRONMENTAL SCIENCES↗

A High‐Resolution Tropical Mesoscale Convective System Reanalysis (TMeCSR)

Abstract Modern global reanalysis products have greatly accelerated meteorological research in synoptic‐to‐planetary‐scale phenomena. However, their use in studying tropical mesoscale convective systems (MCSs) and their regional‐to‐global impact has mostly been limited to supplying initial and boundary conditions for MCS‐resolving simulations and providing information about the large‐scale environments of MCSs. These limitations are due to difficulties in resolving tropical MCS dynamics in the relatively low‐resolution global models and that tropical MCSs often occur over poorly observed regions. In this work, a Tropical MCS‐resolving Reanalysis product (TMeCSR) was created over a region with frequent tropical MCSs. This region spans the tropical Indian Ocean, tropical continental Asia, Maritime Continent, and Western Pacific. TMeCSR is produced by assimilating all‐sky infrared radiances from geostationary satellites and other conventional observations into an MCS‐resolving regional model using the Ensemble Kalman Filter. The resulting observation‐constrained high‐resolution (9‐km grid spacing) data set is available hourly during the boreal summer (June‐August) of 2017, during which widespread severe flooding occurred. Comparisons of TMeCSR and European Center for Medium Range Weather Forecast Reanalysis version 5 (ERA5) against independent satellite retrievals indicate that TMeCSR's cloud and multiscale rain fields are better than those of ERA5. Furthermore, TMeCSR better captured the diurnal variability of rainfall and the statistical characteristics of MCSs. Forecasts initialized from TMeCSR also have more accurate rain and clouds than those initialized from ERA5. The TMeCSR and ERA5 forecasts have similar performances with respect to sounding and surface observations. These results indicate that TMeCSR is a promising MCS‐resolving data set for tropical MCS studies.

Chan, Man‐Yau↗

Systematic Underestimation of Canopy Conductance Sensitivity to Drought by Earth System Models

The response of vegetation canopy conductance (g c ) to changes in moisture availability ($γ^{m}_{gc}$) is a major source of uncertainty in climate projections. While vegetation typically reduces stomatal conductance during drought, accurately modeling how and to what degree stomata respond to changes in moisture availability at global scales is particularly challenging, because no global scale g c observations exist. Here, we leverage a collection of satellite, reanalysis and station-based near-surface air and surface temperature estimates, which are physically and statistically linked to $γ^{m}_{gc}$ due to the local cooling effect of g c through transpiration, to develop a novel emergent constraint of $γ^{m}_{gc}$ in an ensemble of Earth System Models (ESMs). We find that ESMs systematically underestimate $γ^{m}_{gc}$ by ~33%, particularly in grasslands, croplands, and savannas in semi-arid and bordering regions of the Central United States, Central Europe, Southeastern South America, Southern Africa, Eastern Australia, and parts of East Asia. We show that this underestimation occurs because ESMs inadequately reduce g c when soil moisture decreases. As g c controls carbon, water and energy fluxes, the misrepresentation of modeled $γ^{m}_{gc}$ contributes to biases in ESM projections of gross primary production, transpiration, and temperature during droughts. Our results suggest that the severity and duration of droughts may be misrepresented in ESMs due to the impact of sustained g c on both soil moisture dynamics and the biosphere-atmosphere feedbacks that affect local temperatures and regional weather patterns.

54 ENVIRONMENTAL SCIENCES↗