The role of the mean state on MJO simulation in CESM2 ensemble simulation
Explore the source record for details and available documents.
Daehyun Kang↗
SEARCH · Search NASA
Results for “ensemble simulation”
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.
The macroscopic behavior of cumulus ensembles simulated by a cumulus ensemble model
A number of simulation experiments are performed with the 2D UCLA cumulus ensemble model to study the macroscopic behavior of cumulus convection under a variety of different large-scale and underlying surface conditions. Specifically, the modulation of cumulus activity by the imposed large-scale processes and the eddy kinetic energy (EKE) budget are investigated in detail. In all simulations, cumulus convection is rather strongly modulated by large-scale advective processes in spite of the existence of some nonmodulated high-frequency fluctuations. The modulation exhibits some phase delays, however, when the basic wind shear is strong. This is presumably associated with the existence of mesoscale convective organization. The EKE budget analysis shows that the net eddy buoyancy generation rate is nearly zero for a wind range of cumulus ensembles.
Xu, Kuan-Man↗
Ensemble Simulation Techniques and Fast Randomized Algorithms
The major goals of the project were to develop and analyze new ensemble simulation techniques, including trajectory stratification and preconditioned MCMC techniques, as well as develop fast numerical linear algebra techniques closely related to ensemble simulation ideas. The trajectory stratification techniques involve simulating in parallel short trajectory fragments of a Markov process confined to a specific region of space‐time and then patching together the statistics gathered to assemble estimates of very general dynamical properties. We have also developed this approach for rare event simulation and extended the techniques to applications requiring a more general framework (such as electronic structure calculations). The preconditioned MCMC techniques involve simulating multiple Markov chains in parallel and then using information from the ensemble to speed the mixing of each individual chain. The fast randomized linear algebra methods are motivated by the diffusion Monte Carlo technique, but are applicable to finding the dominant eigenvalue of (almost) general matrices. For most non‐negative matrices, the schemes result in an error (compared to the power method) that is constant in the dimension of the problem. For more general matrices, we see a very clear sublinear cost trend in computational tests.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Bridging Hydrological Ensemble Simulation and Learning Using Deep Neural Operators
Ensemble-based simulation and learning (ESnL) has long been used in hydrology for parameter inference, but computational demands of process-based ESnL can be quite high. To address this issue, we propose a deep neural operator learning approach. Neural operators are generic machine learning algorithms that can learn functional mappings between infinite-dimensional spaces, providing a highly flexible tool for scientific machine learning. Our approach is built upon DeepONet, a specific deep neural operator, and is designed to address several common problems in hydrology, namely, model parameter estimation, prediction at ungaged locations, and uncertainty quantification. Here we demonstrate the effectiveness of our DeepONet-based workflow using an existing large model ensemble created for an eastern U.S. watershed that is instrumented with 10 streamflow gages. Results suggest DeepONet achieves high efficiency in learning an ML surrogate model from the model ensemble, with the modified Kling-Gupta Efficiency exceeding 0.9 on holdout test sets. Parameter inference, carried out using the trained DeepONet surrogate model and genetic algorithm, also yields robust results. Additionally, we formulate and train a separate DeepONet model for physics-informed, seq-to-seq streamflow forecasting, which further reduces biases in the pre-trained DeepONet surrogate model. While this study focuses primarily on a single watershed, our approach is general and may be extended to enable learning from model ensembles across multiple basins or models. Thus, this research represents a significant contribution to the application of hybrid machine learning in hydrology.
54 ENVIRONMENTAL SCIENCES↗
Comparing multi-model ensemble simulations with observations and decadal projections of upper atmospheric variations following the Hunga eruption
The Hunga Tonga–Hunga Ha'apai Model–Observation Comparison (HTHH–MOC) project aims to comprehensively investigate the evolution of volcanic water vapor and sulfur emissions and their subsequent atmospheric impacts and underlying response mechanisms using state-of-the-art global climate models. This study evaluates multi-model ensemble simulations participating in the HTHH–MOC free-run experiment with climate projections for 10 years (2022–2032). Model results are evaluated against satellite observations to assess their ability to reproduce the observed evolution of stratospheric water vapor, aerosols, temperature, and ozone from 2022 to 2024. The participating models accurately capture the observed distribution patterns and associated upper atmospheric responses, providing confidence for their future projections. Model simulations suggest that the Hunga eruption-induced stratospheric water vapor anomaly lasts 4–7 years, with a water vapor e-folding time of 31–43 months. This prolonged water vapor perturbation leads to significant stratospheric and mesospheric cooling, resulting in significant ozone loss in the upper stratosphere and lower mesosphere for 7–10 years. Comparisons between simulations with both SO 2 and H 2 O emissions and those with H 2 O-only emissions indicate that the pronounced dipole response with upper-stratospheric cooling and lower-stratospheric warming is driven by the combined effects of SO 2 and H 2 O injections. These results highlight the prolonged atmospheric impacts of the Hunga eruption and the potential critical role of stratospheric water vapor in modulating long-term atmospheric chemistry and dynamics.
Zhuo, Zhihong [Univ. of Quebec, Montreal, QC (Cana↗
A 12-year (1987-1998) Ensemble Simulation of the US Climate with a Variable Resolution Stretched Grid GCM
The variable-resolution stretched-grid (SG) GEOS (Goddard Earth Observing System) GCM has been used for limited ensemble integrations with a relatively coarse, 60 to 100 km, regional resolution over the U.S. The experiments have been run for the 12-year period, 1987-1998, that includes the recent ENSO cycles. Initial conditions 1-2 days apart are used for ensemble members. The goal of the experiments is analyzing the long-term SG-GCM ensemble integrations in terms of their potential in reducing the uncertainties of regional climate simulation while producing realistic mesoscales. The ensemble integration results are analyzed for both prognostic and diagnostic fields. A special attention is devoted to analyzing the variability of precipitation over the U.S. The internal variability of the SG-GCM has been assessed. The ensemble means appear to be closer to the verifying analyses than the individual ensemble members. The ensemble means capture realistic mesoscale patterns, especially those of induced by orography. Two ENSO cycles have been analyzed in terms their impact on the U.S. climate, especially on precipitation. The ability of the SG-GCM simulations to produce regional climate anomalies has been confirmed. However, the optimal size of the ensembles depending on fine regional resolution used, is still to be determined. The SG-GCM ensemble simulations are performed as a preparation or a preliminary stage for the international SGMIP (Stretched-Grid Model Intercomparison Project) that is under way with participation of the major centers and groups employing the SG-approach for regional climate modeling.
Fox-Rabinovitz, Michael S.↗
Diurnal evolution of non-precipitating marine stratocumuli in a large-eddy simulation ensemble
Abstract. We explore the cloud system evolution of non-precipitating marine stratocumuli with a focus on the impacts of the diurnal cycle and free-tropospheric (FT) humidity based on an ensemble of 244 large-eddy simulations generated by perturbing initial thermodynamic profiles and aerosol conditions. Cases are categorized based on their degree of decoupling and the cloud liquid water path (LWPc, based on model columns with cloud optical depths greater than one). A budget analysis method is proposed to analyze the evolution of cloud water in both coupled and decoupled boundary layers. More coupled clouds start with a relatively low LWPc and cloud fraction (fc) but experience the least decrease in LWPc and fc during the daytime. More decoupled clouds undergo greater daytime reduction in LWPc and fc, especially those with higher LWPc at sunrise because they suffer from faster weakening of net radiative cooling. During the nighttime, a positive correlation between FT humidity and the LWPc emerges, consistent with higher FT humidity reducing both radiative cooling and the humidity jump, both of which reduce entrainment and increase LWPc. The LWPc is more likely to decrease during the nighttime for a larger LWPc and greater inversion base height (zi), conditions under which entrainment dominates as turbulence develops. In the morning, the rate of the LWPc reduction depends on the LWPc at sunrise, zi, and the degree of decoupling, with distinct contributions from subsidence and radiation.
Chen, Yao-Sheng (ORCID:0000000208354132)↗
Ensemble Simulations on Leadership Computing Systems
Scientific productivity can be enhanced through workflow management tools, relieving large High Performance Computing (HPC) system users from the tedious tasks of scheduling and designing the complex computational execution of scientific applications. This paper presents a study on the usage of ensemble workflow tools to accelerate science using the Summit and Frontier supercomputing systems. The research aims to connect science domain simulations using Oak Ridge Leadership Computing Facility (OLCF) supercomputing platforms with ensemble workflow methods in order to accelerate HPC-enabled discovery and boost scientific impact. We present the coupling, porting and optimization of Radical-Cybertools on three applications: Chroma, NAMD and LAMMPS. The tools augment traditional HPC monolithic runs with a pilot scheduler. Lessons-learned are discussed for physics, biology and materials science applications. We discuss intrinsic limitations of coupling and porting ensemble workflow tools to applications that run on large HPC systems. The origins of technical challenges and their solutions developed during the implementation process are discussed. Data management strategies, OLCF’s policies for ensembles, and natively supported workflow tools are also summarized.
Georgiadou, Antigoni [ORNL] (ORCID:000000020977631↗
Ensemble Simulations of the 2012 July 12 Coronal Mass Ejection with the Constant-turn Flux Rope Model
Flux-rope-based magnetohydrodynamic modeling of coronal mass ejections (CMEs) is a promising tool for prediction of the CME arrival time and magnetic field at Earth. In this work, we introduce a constant-turn flux rope model and use it to simulate the 2012 July 12 16:48 CME in the inner heliosphere. We constrain the initial parameters of this CME using the graduated cylindrical shell (GCS) model and the reconnected flux in posteruption arcades. We correctly reproduce all the magnetic field components of the CME at Earth, with an arrival time error of approximately 1 hr. We further estimate the average subjective uncertainties in the GCS fittings by comparing the GCS parameters of 56 CMEs reported in multiple studies and catalogs. We determined that the GCS estimates of the CME latitude, longitude, tilt, and speed have average uncertainties of 5°.74, 11°.23, 24°.71, and 11.4%, respectively. Using these, we have created 77 ensemble members for the 2012 July 12 CME. We found that 55% of our ensemble members correctly reproduce the sign of the magnetic field components at Earth. We also determined that the uncertainties in GCS fitting can widen the CME arrival time prediction window to about 12 hr for the 2012 July 12 CME. On investigating the forecast accuracy introduced by the uncertainties in individual GCS parameters, we conclude that the half-angle and aspect ratio have little impact on the predicted magnetic field of the 2012 July 12 CME, whereas the uncertainties in longitude and tilt can introduce relatively large spread in the magnetic field predicted at Earth.
Talwinder Singh↗
Diurnal evolution of non-precipitating marine stratocumuli in a large-eddy simulation ensemble
Explore the source record for details and available documents.
54 ENVIRONMENTAL SCIENCES↗
The Impact of the Atlantic Cold Tongue on West African Monsoon Onset in Regional Model Simulations for 1998-2002
The Atlantic cold tongue (ACT) develops during spring and early summer near the Equator in the Eastern Atlantic Ocean and Gulf of Guinea. The hypothesis that the ACT accelerates the timing of West African monsoon (WAM) onset is tested by comparing two regional climate model (RM3) simulation ensembles. Observed sea surface temperatures (SST) that include the ACT are used to force a control ensemble. An idealized, warm SST perturbation is designed to represent lower boundary forcing without the ACT for the experiment ensemble. Summer simulations forced by observed SST and reanalysis boundary conditions for each of five consecutive years are compared to five parallel runs forced by SST with the warm perturbation. The article summarizes the sequence of events leading to the onset of the WAM in the Sahel region. The representation of WAM onset in RM3 simulations is examined and compared to Tropical Rainfall Measuring Mission (TRMM), Global Precipitation Climatology Project (GPCP) and reanalysis data. The study evaluates the sensitivity of WAM onset indicators to the presence of the ACT by analysing the differences between the two simulation ensembles. Results show that the timing of major rainfall events and therefore theWAM onset in the Sahel are not sensitive to the presence of the ACT. However, the warm SST perturbation does increase downstream rainfall rates over West Africa as a consequence of enhanced specific humidity and enhanced northward moisture flux in the lower troposphere.
Regional model simulation↗
Impacts of Atmosphere-Ocean Coupling on Southern Hemisphere Climate Change
Climate in the Southern Hemisphere (SH) has undergone significant changes in recent decades. These changes are closely linked to the shift of the Southern Annular Mode (SAM) towards its positive polarity, which is driven primarily by Antarctic ozone depletion. There is growing evidence that Antarctic ozone depletion has significant impacts on Southern Ocean circulation change. However, it is poorly understood whether and how ocean feedback might impact the SAM and climate change in the SH atmosphere. This outstanding science question is investigated using the Goddard Earth Observing System Coupled Atmosphere-Ocean-Chemistry Climate Model(GEOS-AOCCM).We perform ensemble simulations of the recent past (1960-2010) with and without the interactive ocean. For simulations without the interactive ocean, we use sea surface temperatures and sea ice concentrations produced by the interactive ocean simulations. The differences between these two ensemble simulations quantify the effects of atmosphere-ocean coupling. We will investigate the impacts of atmosphere-ocean coupling on stratospheric processes such as Antarctic ozone depletion and Antarctic polar vortex breakup. We will address whether ocean feedback affects Rossby wave generation in the troposphere and wave propagation into the stratosphere. Another focuson this study is to assess how ocean feedback might affect the tropospheric SAM response to Antarctic ozone depletion
Chemistry-Climate Modeling↗
Present and future aerosol impacts on Arctic climate change in the GISS-E2.1 Earth system model
The Arctic is warming 2 to 3 times faster than the global average, partly due to changes in short-lived climate forcers (SLCFs) including aerosols. In order to study the effects of atmospheric aerosols in this warming, recent past (1990–2014) and future (2015–2050) simulations have been carried out using the GISS-E2.1 Earth system model to study the aerosol burdens and their radiative and climate impacts over the Arctic (>60°N), using anthropogenic emissions from the Eclipse V6b and the Coupled Model Intercomparison Project Phase 6 (CMIP6) databases, while global annual mean greenhouse gas concentrations were prescribed and kept fixed in all simulations. Results showed that the simulations have underestimated observed surface aerosol levels, in particular black carbon (BC) and sulfate (SO2−4), by more than 50 %, with the smallest biases calculated for the atmosphere-only simulations, where winds are nudged to reanalysis data. CMIP6 simulations performed slightly better in reproducing the observed surface aerosol concentrations and climate parameters, compared to the Eclipse simulations. In addition, simulations where atmosphere and ocean are fully coupled had slightly smaller biases in aerosol levels compared to atmosphere-only simulations without nudging. Arctic BC, organic aerosol (OA), and SO2−4 burdens decrease significantly in all simulations by 10 %–60 % following the reductions of 7 %–78 % in emission projections, with the Eclipse ensemble showing larger reductions in Arctic aerosol burdens compared to the CMIP6 ensemble. For the 2030–2050 period, the Eclipse ensemble simulated a radiative forcing due to aerosol–radiation interactions (RFARI) of −0.39 ± 0.01 W/sq. m, which is −0.08 W/sq. m larger than the 1990–2010 mean forcing (−0.32 W/sq. m), of which −0.24 ± 0.01 W/sq. m was attributed to the anthropogenic aerosols. The CMIP6 ensemble simulated a RFARI of −0.35 to −0.40 W/sq. m for the same period, which is −0.01 to −0.06 W/sq. m larger than the 1990–2010 mean forcing of −0.35 W/sq. m. The scenarios with little to no mitigation (worst-case scenarios) led to very small changes in the RFARI, while scenarios with medium to large emission mitigations led to increases in the negative RFARI, mainly due to the decrease in the positive BC forcing and the decrease in the negative SO2−4 forcing. The anthropogenic aerosols accounted for −0.24 to −0.26 W/sq. m of the net RFARI in 2030–2050 period, in Eclipse and CMIP6 ensembles, respectively. Finally, all simulations showed an increase in the Arctic surface air temperatures throughout the simulation period. By 2050, surface air temperatures are projected to increase by 2.4 to 2.6 °C in the Eclipse ensemble and 1.9 to 2.6 °C in the CMIP6 ensemble, compared to the 1990–2010 mean. Overall, results show that even the scenarios with largest emission reductions leads to similar impact on the future Arctic surface air temperatures and sea-ice extent compared to scenarios with smaller emission reductions, implying reductions of greenhouse emissions are still necessary to mitigate climate change.
Aerosols↗
The Statistical Emulators of GGCMI Phase 2: Responses of Year-to-Year Variation of Crop Yield to CO 2 , Temperature, Water and Nitrogen Perturbations
Understanding the impact of climate change on year-to-year variation of crop yield is critical to global food stability and security. While crop model emulators are believed to be lightweight tools to replace the models, few emulators have been developed to capture such interannual variation of crop yield in response to climate variability. In this study, we developed a statistical emulator with a machine learning algorithm to reproduce the response of year-to-year variation of four crop yields to CO 2 (C), temperature (T), water (W), and nitrogen (N) perturbations defined in the Global Gridded Crop Model Intercomparison Project (GGCMI) phase 2. The emulators were able to explain more than 52 % of the variance of simulated yield and performed well in capturing the year-to-year variation of global average and gridded crop yield over current croplands in the baseline. With the changes in CO 2 –temperature–water–nitrogen (CTWN) perturbations, the emulators could reproduce the year-to-year variation of crop yield well over most current cropland. The variation of R and the mean absolute error was small under the single CTWN perturbations and dual-factor perturbations. These emulators thus provide statistical response surfaces of yield, including both its mean and interannual variability, to climate factors. They could facilitate spatiotemporal downscaling of crop model simulation, projecting the changes in crop yield variability in the future and serving as a lightweight tool for multi-model ensemble simulation. The emulators enhanced the flexibility of crop yield estimates and expanded the application of large-ensemble simulations of crop yield under climate change.
Climate change↗
A Greening Future Elevates Flash Drought Risk in Northern Mid‐to‐High Latitudes
Flash droughts have become a growing concern, as they can emerge rapidly and increase the risk of crop failure. Although past studies have investigated the meteorological drivers and future changes of flash drought, why flash drought is more frequent over humid and vegetated regions remains underexplored. This study delves further into the mechanism by which vegetation regulates flash drought and its future change using observations from multiple data sets and large ensemble simulations from three Earth system models. On an interannual timescale, both observations and simulations show robust increases in flash drought frequency and a higher flash-to-sub-seasonal drought ratio during spring or antecedent conditions with dense vegetation, supporting the important role of vegetation in flash drought occurrence, especially in the northern mid-to-high latitudes. In the latter regions, the large ensemble simulations show robust increases in flash drought (e.g., 67% and 46% increases in Eastern U.S. and North Asia in 2050–2100 relative to 1950–2000 under the high emission scenario), where the growing season is lengthening. Although greening might suggest reduced drought stress, it drives precipitation-soil moisture-evapotranspiration decoupling by increasing evapotranspiration partitioning to transpiration. As transpiration can access deep soil water through the plant root system, its increased portion can weaken the constraints of concurrent precipitation on evapotranspiration, thus accelerating soil moisture depletion under high evaporative demand, driving a slow-to-rapid drought transition. How vegetation regulates flash drought by regulating surface moisture budget is supported by observations and simulations. Although warming supports early planting, agriculture may increasingly be threatened by surging flash drought risk.
Drought↗
Data for Physiological Controls on Carbon Fluxes and Biomass Production in Miscanthus: Insights From a Process- Based Agroecosystem Model
Biomass crops serve as essential feedstocks for renewable energy and bioproducts and play a critical role in achieving lower emissions in the transportation sector. However, dedicated perennial biomass crops such as Miscanthus × giganteus (Miscanthus) remain underrepresented in process- based agroecosystem models, limiting robust evaluation of their economic and environmental performance. In this study, we developed a data- constrained representation of the sterile triploid Miscanthus (IL clone) within the process- based model ecosys, integrating global sensitivity analysis, ensemble simulation, and parameter calibration. Planting, harvesting, and fertilization practices consistent with field management were incorporated, and phenology was constrained using PhenoCam- derived Green Chromatic Coordinate (GCC) data. Using the Morris global sensitivity analysis method, we identified 11 key physiological parameters governing plant carbon, water, and nutrient relations, particularly processes associated with CO2 assimilation. We then conducted ensemble simulations by perturbing these parameters and calibrated the model against eddy covariance fluxes and field- measured biomass. Building on the calibrated operating state, parameter- response analyses show that different photosynthetic processes influence productivity in different ways. Protein allocation determines whether productivity increases toward a higher level, whereas electron transport capacity controls additional gains once protein allocation approaches saturation. These findings demonstrate that parameter importance depends on physiological context and on which photosynthetic processes remain limiting. Calibration and validation against observations show that ecosys can reliably reproduce carbon and water fluxes, as well as both above- and belowground biomass, with post- calibration GPP R2 improving from 0.67 to 0.95 during the calibration period and remaining high during validation (R2 = 0.95). These results provide a mechanistic foundation for regional simulations and sustainable bioenergy assessments.
Carbon↗
Physiological Controls on Carbon Fluxes and Biomass Production in Miscanthus: Insights From a Process‐Based Agroecosystem Model
Biomass crops serve as essential feedstocks for renewable energy and bioproducts and play a critical role in achieving lower emissions in the transportation sector. However, dedicated perennial biomass crops such as Miscanthus × giganteus (Miscanthus) remain underrepresented in process-based agroecosystem models, limiting robust evaluation of their economic and environmental performance. In this study, we developed a data-constrained representation of the sterile triploid Miscanthus (IL clone) within the process-based model ecosys, integrating global sensitivity analysis, ensemble simulation, and parameter calibration. Planting, harvesting, and fertilization practices consistent with field management were incorporated, and phenology was constrained using PhenoCam-derived Green Chromatic Coordinate (GCC) data. Using the Morris global sensitivity analysis method, we identified 11 key physiological parameters governing plant carbon, water, and nutrient relations, particularly processes associated with CO 2 assimilation. We then conducted ensemble simulations by perturbing these parameters and calibrated the model against eddy covariance fluxes and field-measured biomass. Building on the calibrated operating state, parameter-response analyses show that different photosynthetic processes influence productivity in different ways. Protein allocation determines whether productivity increases toward a higher level, whereas electron transport capacity controls additional gains once protein allocation approaches saturation. These findings demonstrate that parameter importance depends on physiological context and on which photosynthetic processes remain limiting. Calibration and validation against observations show that ecosys can reliably reproduce carbon and water fluxes, as well as both above- and belowground biomass, with post-calibration GPP R 2 improving from 0.67 to 0.95 during the calibration period and remaining high during validation (R 2 = 0.95). These results provide a mechanistic foundation for regional simulations and sustainable bioenergy assessments.
ecosys↗