Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ensemble simulations”

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 109 records · Page 6

WRF-ELM v1.0: a regional climate model to study land–atmosphere interactions over heterogeneous land use regions

Abstract. The Energy Exascale Earth System Model (E3SM) Land Model (ELM) is a state-of-the-art land surface model that simulates the intricate interactions between the terrestrial land surface and other components of the Earth system. Originating from the Community Land Model (CLM) version 4.5, ELM has been under active development, with added new features and functionality, including plant hydraulics, radiation–topography interaction, subsurface multiphase flow, and more explicit land use and management practices. This study integrates ELM v2.1 with the Weather Research and Forecasting (WRF; WRF-ELM) model through a modified Lightweight Infrastructure for Land Atmosphere Coupling (LILAC) framework, enabling affordable high-resolution regional modeling by leveraging ELM's innovative features alongside WRF's diverse atmospheric parameterization options. This framework includes a top-level driver for variable communication between WRF and ELM and Earth System Modeling Framework (ESMF) caps for the WRF atmospheric component and ELM workflow control, encompassing initialization, execution, and finalization. Importantly, this LILAC–ESMF framework demonstrates a more modular approach compared to previous coupling efforts between WRF and land surface models. It maintains the integrity of ELM's source code structure and facilitates the transfer of future developments in ELM to WRF-ELM. To test the ability of the coupled model to capture land–atmosphere interactions over regions with a variety of land uses and land covers, we conducted high-resolution (4 km) WRF-ELM ensemble simulations over the Great Lakes region (GLR) in the summer of 2018 and systematically compared the results against observations, reanalysis data, and WRF-CTSM (WRF coupled with the Community Terrestrial Systems Model). In general, the coupled WRF-ELM model has reasonably captured the spatial distribution of surface state variables and fluxes across the GLR, particularly over the natural vegetation areas. The evaluation results provide a baseline reference for further improvements in ELM in the regional application of high-resolution weather and climate predictions. Our work serves as an example to the model development community for expanding an advanced land surface model's capability to represent fully-coupled land–atmosphere interactions at fine spatial scales. The development and release of WRF-ELM marks a significant advancement for the ELM user community, providing opportunities for fine-scale regional representation, parameter calibration in coupled mode, and examination of new schemes with atmospheric feedback.

54 ENVIRONMENTAL SCIENCES↗

Seasonally Dependent Future Changes in the U.S. Midwest Hydroclimate and Extremes

This study investigates the responses of the hydroclimate and extremes in the U.S. Midwest to global warming, based on ensemble projections of phase 6 of the Coupled Model Intercomparison Project and the multimodel initial-condition large-ensemble simulations. The precipitation response features a seasonally dependent change with increased precipitation in April–May but reduced precipitation in July–August. In this work, the late-spring wetting is attributed to the enhanced low-level moisture-transporting southerlies, which are induced by regional sea level pressure anomalies linked to the poleward shift of the North American westerly jet (NAWJ). The late-summer drying is attributed to the weakened storm track, which is also linked to the poleward NAWJ shift. The seasonally dependent future changes of the Midwest precipitation are analogous to its climatological seasonal progression, which increases over late spring as the NAWJ approaches the Midwest and decreases over late summer as the NAWJ migrates away. In response to the mean precipitation changes, extremely wet late springs (April–May precipitation above the 99th percentile of the historical period) and extremely dry late summers (below the 1st percentile) will occur much more frequently, implying increased late-spring floods and late-summer droughts. Future warming in the Midwest is amplified in late summer due to the reduced precipitation. With amplified background warming and increased occurrence, future late-summer droughts will be more devastating. Our results highlight that, under a time-invariant poleward jet shift, opposite precipitation changes arise before and after the peak rainy month, leading to substantial increases in the subseasonal extremes. The severity of such climate impacts is obscured in projections of the rainy-season mean.

54 ENVIRONMENTAL SCIENCES↗

Land cover change-induced decline in terrestrial gross primary production over the conterminous United States from 2001 to 2016

As one of the most dynamic aspects of global environmental change, land cover change (LCC) has a profound impact on terrestrial carbon sequestration. However, LCC-induced carbon fluxes are still the most uncertain terms in global and regional carbon budgets. Ecosystem gross primary production (GPP) is the total carbon uptake by vegetation through photosynthesis, serving as a major control on ecosystem function and land carbon balance during and after the modification of the land surface. However, accurately capturing LCC-induced GPP changes requires both high-quality land cover data and controlling for variation driven by other environmental factors such as climate. In this study, we comprehensively examined the effects of LCC on annual GPP trends over the conterminous United States (CONUS) from 2001 to 2016 using the USGS National Land Cover Database, a remote sensing-driven ecosystem model, and the Google Earth Engine cloud computing platform. We designed a series of model experiments to identify LCC effects on GPP by controlling climate effects. During the study period, LCC exerted a strong negative effect on total GPP across the CONUS ([-2.2, -1.8] Tg C yr -2 ), while climate had smaller positive effects ([0.17, 0. 92] Tg C yr -2 ). The LCC-induced reduction of GPP was mainly caused by net forest loss ([-1.98, -1.39] Tg C yr -2 ) and urban expansion ([-2.03, -1.92] Tg C yr -2 ), but was partially offset by increases in crop area ([+0.66, +0.79] Tg C yr -2 ). Ensemble simulations from TRENDY did not capture the strong negative LCC influences on GPP, likely due to limitations of the adopted land use/cover data. Overall, our study provides a novel perspective on LCC-induced GPP changes, which could help to improve our understanding of ecosystem function changes and constrain the estimation of land carbon balance in the context of anthropogenic activity and climate change.

54 ENVIRONMENTAL SCIENCES↗

Entropy-based feature selection for capturing impacts in Earth system models with abrupt forcing

This paper presents the development of a new entropy-based feature selection method for identifying and quantifying impacts. Here, impacts are defined as statistically significant differences in spatio-temporal fields when comparing datasets with and without an external forcing in an Earth system model. Temporal feature selection is performed by first computing the cross-fuzzy entropy to quantify similarity of patterns between two datasets and then applying changepoint detection to identify regions of statistically constant entropy. The method is used to capture temperate north surface cooling from a 9-member simulation ensemble of the Mt. Pinatubo volcanic eruption, which injected 10 Tg of SO 2 into the stratosphere. The results estimate a mean difference decrease in near surface air temperature of -0.560 K with a 99% confidence interval between -0.864 K and -0.257 K between April and November of 1992, one year following the eruption. A sensitivity analysis with decreasing SO 2 injection revealed that the impact is statistically significant at 5 Tg but not at 3 Tg. Using identified features, a dependency graph model based on a 9-day lag had significantly fewer nodes than a graph based on monthly means. Furthermore, this demonstrates our method’s ability to perform dimension reduction while still uncovering source-to-impact pathways.

Changepoint detection↗

Large Scale Study of Ligand–Protein Relative Binding Free Energy Calculations: Actionable Predictions from Statistically Robust Protocols

The accurate and reliable prediction of protein–ligand binding affinities can play a central role in the drug discovery process as well as in personalized medicine. Of considerable importance during lead optimization are the alchemical free energy methods that furnish an estimation of relative binding free energies (RBFE) of similar molecules. Recent advances in these methods have increased their speed, accuracy, and precision. This is evident from the increasing number of retrospective as well as prospective studies employing them. However, such methods still have limited applicability in real-world scenarios due to a number of important yet unresolved issues. Here, we report the findings from a large data set comprising over 500 ligand transformations spanning over 300 ligands binding to a diverse set of 14 different protein targets which furnish statistically robust results on the accuracy, precision, and reproducibility of RBFE calculations. We use ensemble-based methods which are the only way to provide reliable uncertainty quantification given that the underlying molecular dynamics is chaotic. These are implemented using TIES (Thermodynamic Integration with Enhanced Sampling). Results achieve chemical accuracy in all cases. Ensemble simulations also furnish information on the statistical distributions of the free energy calculations which exhibit non-normal behavior. We find that the “enhanced sampling” method known as replica exchange with solute tempering degrades RBFE predictions. We also report definitively on numerous associated alchemical factors including the choice of ligand charge method, flexibility in ligand structure, and the size of the alchemical region including the number of atoms involved in transforming one ligand into another. Our findings provide a key set of recommendations that should be adopted for the reliable application of RBFE methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Inverse Thermodynamics: Designing Interactions for Targeted Phase Behavior

The traditional goal of inverse self-assembly is to design interactions that drive particles toward a desired target structure. However, achieving successful self-assembly also requires tuning the thermodynamic conditions under which the structure is stable. In this work, we extend the inverse design paradigm to explicitly address this challenge by developing a framework for inverse thermodynamics, i.e., the design of interaction potentials that realize specific thermodynamic behavior. As a step in this direction, using patchy particle mixtures as a model system, we demonstrate how precise control over both bonding topology and bond energetics enables the programming of targeted phase behavior. In particular, we establish design principles for azeotropic demixing and show how to create mixtures that exhibit azeotropy at any prescribed composition. Our predictions are validated through Gibbs-ensemble simulations [Panagiotopoulos, Mol. Phys. 1987, 61, 813−826]. These results highlight the necessity of coupling structural design with thermodynamic engineering, and provide a blueprint for controlling complex phase behavior in multicomponent systems.

Azeotropes↗

The Community Earth System Model version 2 (CESM2)

An overview of the Community Earth System Model version 2 (CESM2) is provided, including a discussion of the challenges encountered during its development and how they were addressed. In addition, an evaluation of a pair of CESM2 long pre-industrial control and historical ensemble simulations is presented. These simulations were performed using the nominal 1° horizontal resolution configuration of the coupled model with both the “low-top” (40 km, with limited chemistry) and “high-top” (140 km, with comprehensive chemistry) versions of the atmospheric component. CESM2 contains many substantial science and infrastructure improvements and new capabilities since its previous major release, CESM1, resulting in improved historical simulations in comparison to CESM1 and available observations. These include major reductions in low latitude precipitation and short-wave cloud forcing biases; better representation of the Madden-Julian Oscillation; better El Niño – Southern Oscillation-related teleconnections; and a global land carbon accumulation trend that agrees well with observationally-based estimates. Most tropospheric and surface features of the low- and high-top simulations are very similar to each other, so these improvements are present in both configurations. CESM2 has an equilibrium climate sensitivity of 5.1-5.3°C, larger than in CESM1, primarily due to a combination of relatively small changes to cloud microphysics and boundary layer parameters. In contrast, CESM2’s transient climate response of 1.9-2.0°C is comparable to that of CESM1. The model outputs from these and many other simulations are available to the research community, and theyrepresent CESM2’s contributions to the Coupled Model Intercomparison Project phase 6 (CMIP6).

Danabasoglu, Gokhan↗

Coupled Climate Responses to Recent Australian Wildfire and COVID-19 Emissions Anomalies Estimated in CESM2

In this work, multiple 50-member ensemble simulations with the Community Earth System Model version 2 are performed to estimate the coupled climate responses to the 2019–2020 Australian wildfires and COVID-19 pandemic policies. The climate response to the pandemic is found to be weak generally, with global-mean net top-of-atmosphere radiative anomalies of +0.23 ± 0.14 W m –2 driving a gradual global warming of 0.05 ± 0.04 K by the end of 2022. While regional anomalies are detectable in aerosol burdens and clear-sky radiation, few significant anomalies exist in other fields due to internal variability. In contrast, the simulated response to Australian wildfires is a strong and rapid cooling, peaking globally at –0.95 ± 0.15 W m –2 in late 2019 with a global cooling of 0.06 ± 0.04 K by mid-2020. Transport of fire aerosols throughout the Southern Hemisphere increases albedo and drives a strong interhemispheric radiative contrast, with simulated responses that are consistent generally with those to a Southern Hemisphere volcanic eruption.

54 ENVIRONMENTAL SCIENCES↗

Does Regional Hydroclimate Change Scale Linearly With Global Warming?

Many aspects of climate change scale linearly with global warming. However, nonlinear changes are possible, especially in the context of hydroclimate, and under emissions scenarios with stabilized global temperature, as aspired to by current climate targets. In CMIP5 and 6, a progressively larger land area shows nonlinear changes as a function of global warming when considering precipitation, evaporation, and soil moisture, with the latter showing nonlinearity over ~50% of global land. Here, using ensemble simulations with the Community Earth System Model 1, in which individual forcing factors are held constant, we illustrate how nonadditive responses to anthropogenic greenhouse gases and industrial and fire-related aerosols can yield complex soil moisture changes in certain regions. This complexity contributes to uncertainty in regional soil moisture projections and suggests that the timing of, as well as model response uncertainty to, future aerosol reductions will have significant impacts on regional hydroclimate change as global temperatures stabilize.

54 ENVIRONMENTAL SCIENCES↗

An Observationally Trained Markov Model for MJO Propagation

A Markovian stochastic model is developed for studying the propagation of the Madden-Julian Oscillation (MJO). This model represents the daily changes in real time multivariate MJO (RMM) indices as random functions of their current state and background conditions. The probability distribution function of the RMM changes is obtained using a machine learning algorithm trained to maximize MJO forecast skills using observed daily indices of RMM and different modes of variability. Skillful forecasts are obtained for lead times between 8 and 27 days. Large ensemble simulations by the stochastic model show that with monsoonal changes in the background state, MJO propagation across the Maritime Continent (MC) is most likely to be disrupted in boreal spring and summer when MJO events propagate from favorable conditions over the Indian Ocean to unfavorable ones over the MC, and predictability is higher during spring and summer when MJO activity is away from the MC region.

54 ENVIRONMENTAL SCIENCES↗

Different Strategies of Stratospheric Aerosol Injection Would Significantly Affect Climate Extreme Mitigation

Abstract Stratospheric aerosol injection (SAI) has been proposed as a potential supplement to mitigate some climate impacts of anthropogenic warming. Using Community Earth System Model ensemble simulation results, we analyze the response of temperature and precipitation extremes to two different SAI strategies: one injects SO 2 at the equator to stabilize global mean temperature and the other injects SO 2 at multiple locations to stabilize global mean temperature as well as the interhemispheric and equator‐to‐pole temperature gradients. Our analysis shows that in the late 21st century, compared with the present‐day climate, both equatorial and multi‐location injection lead to reduced hot extremes in the tropics, corresponding to overcooling of the mean climate state. In mid‐to‐high latitude regions, in comparison to the present‐day climate, substantial decreases in cold extremes are observed under both equatorial and multi‐location injection, corresponding to residual winter warming of the mean climate state. Both equatorial and multi‐location injection reduce precipitation extremes in the tropics below the present‐day level, associated with the decrease in mean precipitation. Overall, for most regions, temperature and precipitation extremes show reduced change in response to multi‐location injection than to equatorial injection, corresponding to reduced mean climate change for multi‐location injection. In comparison with equatorial injection, in response to multi‐location injection, most land regions experience fewer years with significant change in cold extremes from the present‐day level, and most tropical regions experience fewer years with significant change in hot extremes. The design of SAI strategies to mitigate anthropogenic climate extremes merits further study.

54 ENVIRONMENTAL SCIENCES↗

Contributions of Climate Change and ENSO Variability to Future Precipitation Extremes Over California

The El Niño/Southern Oscillation (ENSO) affects the occurrence frequency and intensity of extreme precipitation through modulations of regional heat and moisture fluxes. California experiences particularly strong ENSO influences and models project different to its extreme precipitation. It remains unclear how diverse projections of future precipitation extremes relate to inter-model differences for those changing signals. Here, we use “large ensemble” simulations with multiple climate models along with the Coupled Model Intercomparison Project Phase 5 to investigate the range of precipitation extreme changes over California and the influences from ENSO-related teleconnections. We found that precipitation amount increases are much larger during El Niño relative to La Niña years, mainly caused by the differences in frequency of extreme events during different phases. The ENSO-driven effect is even larger than the overall climate change signal for the most extreme events, implying uncertainties from inter-model differences in ENSO-related SST variability for extreme precipitation changes.

54 ENVIRONMENTAL SCIENCES↗

Distilling the Evolving Contributions of Anthropogenic Aerosols and Greenhouse Gases to Large‐Scale Low‐Frequency Surface Ocean Changes Over the Past Century

Abstract Anthropogenic aerosols (AER) and greenhouse gases (GHG)—the leading drivers of the forced historical change—produce different large‐scale climate response patterns, with correlations trending from negative to positive over the past century. To understand what caused the time‐evolving comparison between GHG and AER response patterns, we apply a low‐frequency component analysis to historical surface ocean changes from CESM1 single‐forcing large‐ensemble simulations. While GHG response is characterized by its first leading mode, AER response consists of two distinct modes. The first one, featuring long‐term global AER increase and global cooling, opposes GHG response patterns up to the mid‐twentieth century. The second one, featuring multidecadal variations in AER distributions and interhemispheric asymmetric surface ocean changes, appears to reinforce the GHG warming effect over recent decades. AER thus can have both competing and synergistic effects with GHG as their emissions change temporally and spatially.

Dong, Yue↗

Assessing Radiative Feedbacks and Their Contribution to the Arctic Amplification Measured by Various Metrics

Arctic amplification (AA), characterized by a more rapid surface air temperature (SAT) warming in the Arctic than the global average, is a major feature of global climate warming. Various metrics have been used to quantify AA based on SAT anomalies, trends, or variability, and they can yield quite different conclusions regarding the magnitude and temporal patterns of AA. This study examines and compares various AA metrics for their temporal consistency in the region north of 70°N from the early twentieth to the early 21st century using observational data and reanalysis products. We also quantify contributions of different radiative feedback mechanisms to AA based on short-term climate variability in reanalysis and model data using the Kernel-Gregory approach. Albedo and lapse rate feedbacks are positive and comparable, with albedo feedback being the leading contributor for all AA metrics. The net cloud feedback, which has large uncertainties, depends strongly on the data sets and AA metrics used. By quantifying the influence of internal variability on AA and related feedbacks based on global climate model ensemble simulations, we find that water vapor and cloud feedbacks are most heavily affected by internal variability.

54 ENVIRONMENTAL SCIENCES↗

Sensitivity of Anthropogenic Cloud Droplet Number Change to Preindustrial Emission Inventories and Physics Parameterizations

Understanding the preindustrial (PI) to present-day (PD) change in cloud droplet number concentration (ΔNd) is important to constrain the anthropogenic influence on clouds and aid future model projections. Perturbed parameter ensemble simulations of PI and PD conditions in the Energy Exascale Earth System Model (E3SMv3) reveal two locations with consistent negative ΔNd signal across parameter space–the Southeastern US and the United Kingdom–despite a trend of positive ΔNd globally from increased fossil fuel consumption. These negative ΔNd signals are driven by higher anthropogenic biofuel emissions in the PI. By varying the PI aerosol mass and/or number emissions across these data sets' uncertainty ranges, we can change the magnitude and sign of ΔNd and the effective radiative forcing from aerosol-cloud interactions (ERFaci) regionally and globally, with near-zero ΔNd leading to near-zero ERFaci. These results highlight the need to constrain PI emission uncertainty to better understand Earth system responses to aerosol-cloud interactions.

Brown, Hunter Y. [University of Wyoming, Laramie, ↗

Emergence of seasonal delay of tropical rainfall during 1979–2019

Fundamental to Earth’s climate, tropical rainfall marks a dominant annual cycle characterized by its amplitude and phase, representing the range between wet and dry seasons and their onset timing, respectively. Under global warming, previous studies revealed enhanced amplitude of rainfall annual cycle over ocean and delayed phase over land in the tropics, underpinned by solid theories. It is unclear however whether these changes have already emerged in the observational records. Here, we report a robust delay of 3.02 and 3.43 days over the northern tropical land and Sahel respectively during 1979-2018 based on five observation datasets. Furthermore, most of the delay is captured by multiple ensemble simulations driven by external forcings, with greenhouse gases (GHG) more important over the northern tropical land and anthropogenic aerosols (AER) more important over the Sahel. Both GHG and AER delay the annual cycle of regional precipitation similarly by producing a moister and more ocean-like atmosphere, thus increasing its lag in response to seasonal solar forcing. Our findings provide new evidence of the human footprints on Earth’s hydrological cycle at continental and regional scales. As GHG and AER continue to increase and decrease respectively, these seasonal delays will be further amplified in the future.

54 ENVIRONMENTAL SCIENCES↗

Enhanced risk of concurrent regional droughts with increased ENSO variability and warming

Spatially compounding extremes pose substantial threats to globally interconnected socio-economic systems. In this study we use multiple large ensemble simulations of the high-emissions scenario to show increased risk of compound droughts during the boreal summer over ten global regions. Relative to the late twentieth century, the probability of compound droughts increases by ~40% and ~60% by the middle and late twenty-first century, respectively, with a disproportionate increase in risk across North America and the Amazon. These changes contribute to an approximately ninefold increase in agricultural area and population exposure to severe compound droughts with continued fossil-fuel dependence. ENSO is the predominant large-scale driver of compound droughts with 68% of historical events occurring during El Niño or La Niña conditions. With ENSO teleconnections remaining largely stationary in the future, a ~22% increase in frequency of ENSO events combined with projected warming drives the elevated risk of compound droughts.

54 ENVIRONMENTAL SCIENCES↗

Emergent constraints on future methane emissions from global wetlands

Future methane (CH 4 ) emissions from natural wetlands are predicted to increase due to global warming, leading to positive feedback on climate change. However, the magnitude of this increase remains highly uncertain. Here we present novel ensemble simulations of seven state-of-the-art terrestrial biosphere models to estimate wetland CH 4 emissions (eCH 4 ) during the twenty-first century. Our estimates suggest that for every 1 °C increase in global land surface temperature, there is a 24 ± 10 Tg CH 4 yr −1 increase in eCH 4 . We also identify an emergent relationship between contemporary temperature dependence and projected eCH 4 . When constrained by 163 site-year eddy-covariance measurements of eCH 4 , we show that wetland emissions can increase by 50–60% by the 2090s relative to the 2010s under a high-warming scenario. The projected decadal increase in eCH 4 from the 2010–2019 baseline to the 2030s would very likely (90% probability) offset an amount equivalent in scale to 8–10% of anthropogenic eCH 4 at the 2020 level, comparable to the reductions committed under the Global Methane Pledge. However, the constraint is dominated by mid- and high-latitude observations, with limited tropical coverage, and uncertainties in projected wetland inundation contribute substantially to uncertainty in eCH 4 . Our findings reduce the uncertainty in projected wetland methane–climate feedback and highlight its potential impacts on methane mitigation efforts to slow global warming.

Zhang, Zhen [Chinese Academy of Sciences (CAS), Be↗