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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 505 records · Page 28

All sources NASA NTRS DOE OSTI

NASA NTRS2020

Mechanisms Linked to Recent Ozone Decreases in the Northern Hemisphere Lower Stratosphere

Previous studies have robustly identified a decrease since 1998 in lower stratospheric (LS) ozone in the Northern Hemisphere (NH). While this ozone decrease is qualitatively explained as resulting from changes in the large‐scale circulation, there is not yet a quantitative mechanistic explanation of these changes. Here we explore the drivers of recent ozone changes using two different configurations of the Goddard Earth Observing System (GEOS) general circulation model. The first configuration of GEOS includes a full chemistry module and is constrained with meteorological fields from the Modern‐Era Retrospective analysis for Research and Applications, Version 2 (MERRA‐2). This configuration (M2GMI) is used to analyze an idealized tracer that covaries closely with ozone on interannual and decadal timescales, revealing that recent ozone decreases in the NH subtropics are associated with a poleward expansion of upwelling in the NH LS, with reduced (enhanced) downwelling over northern subtropics (midlatitudes). The second configuration of GEOS is a free‐running version of the GEOS Chemistry‐Climate Model (CCM) that is used to perform a ten‐member ensemble of free‐running simulations. Comparisons of the two configurations reveal that, while the free‐running model can produce negative ozone changes in the NH LS, the magnitude of these changes is significantly weaker, relative to both M2GMI and MERRA‐2; moreover, these weaker ozone decreases are consistent with weaker simulated changes in the residual circulation. Finally, we examine the GEOS model results in the broader context of the hindcast simulations performed as part of Phase 1 of the Chemistry Climate Modeling Initiative. We show that the majority of the free‐running simulations considered here also exhibit weaker long‐term residual circulation changes, compared to reanalyses.

Ozone↗

NASA NTRS2020

On the Role of Trend and Variability of Hydroxyl Radical (OH) in the Global Methane Budget

Decadal trends and interannual variations in the hydroxyl radical (OH), while poorly constrained at present, are critical for understanding the observed evolution of atmospheric methane (CH4). Through analyzing the OH fields simulated by the model ensemble of the Chemistry-Climate Model Initiative (CCMI), we find (1) the negative OH anomalies during the El Niño years mainly corresponding to the enhanced carbon monoxide (CO) emissions from biomass burning and (2) a positive OH trend during 1980-2010 dominated by the elevated primary production and the reduced loss of OH due to decreasing CO after 2000. Both two-box model inversions and variational 4D inversions suggest that ignoring the negative anomaly of OH during the El Niño years leads to a large overestimation of the increase in global CH4 emissions by up to10±3Tg yr-1to match the observed CH4 increase over these years. Not accounting for the increasing OH trends given by the CCMI models leads to an underestimation of the CH4 emission increase by 23±9Tg yr-1from 1986 to 2010. The variational inversion estimated CH4 emissions show that the tropical regions contribute most to the uncertainties related to OH. This study highlights the significant impact of climate and chemical feedbacks related to OH on the top-down estimates of the global CH4 budget.

Hydroxyl radical (OH)↗

NASA NTRSDate not supplied

Reconstructing Subgrid Cloud Variability Guided by CloudSat/CALIPSO Observations

Predicting realistic cloud subgrid variability remains a challenge for Global Climate Models (GCMs) even though it can play an important role for proper representation of processes pertaining to cloud microphysics, precipitation, and radiation, but also for comparisons with satellite observations which are of much higher resolution than model grids. The diagnostic approach of subgrid cloud variability is commonly handled by subcolumn cloud generators. For a specific GCM, one ideally wants the subgrid variability used for comparisons with satellite observations to be created by the same generator and with the same rules as the one used for model integration. With this in mind, we have embarked in an effort to test and improve cloud subcolumn generators appropriate for GCMs. For this purpose, we use cloud (hydrometeor) products from active observations by the CloudSat radar (CPR) and the CALIPSO lidar (CALIOP). Cloud products from active sensors while suffering significant sampling and coverage drawbacks have the advantage of resolving both horizontal and vertical variability. The main question is: given a profile of cloud condensate mean and variance, can we create a subgrid cloud field that is statistically similar to the observed subgrid cloud field? By “statistically”, we suggest that we do not aspire to reproduce “well” each individual subgrid cloud field of a GCM-scale region, but that our generator performs well for a large ensemble of cases. We simulate radar, passive imager, and radiation flux fields from the observed 2D cloud fields and create one-point statistics; we use the profiles of cloud condensate mean and variance as input to the generator to create subgrid cloud fields; we compare the statistics; we adjust the rules of the generator to create the best possible agreement between the radar, imager and radiation field statistics. In this process, the active observations have actually a dual role: they provide the actual subgrid cloud field which can be reduced to a profile of mean and variance used by the subcolumn generator, but they also provide the rules needed by the generator, such as measures (e.g, decorrelation length) of the vertical overlap of cloud fraction and of the condensate horizontal variability. This presentation will show our progress using Cloudsat and CALIPSO products in this dual fashion.

cloud↗

NASA NTRSDate not supplied

Reconstructing Subgrid Cloud Variability Guided by CloudSat/CALIPSO Observations

Predicting realistic cloud subgrid variability remains a challenge for Global Climate Models (GCMs) even though it can play an important role for proper representation of processes pertaining to cloud microphysics, precipitation, and radiation, but also for comparisons with satellite observations which are of much higher resolution than model grids. The diagnostic approach of subgrid cloud variability is commonly handled by subcolumn cloud generators. For a specific GCM, one ideally wants the subgrid variability used for comparisons with satellite observations to be created by the same generator and with the same rules as the one used for model integration. With this in mind, we have embarked in an effort to test and improve cloud subcolumn generators appropriate for GCMs. For this purpose, we use cloud (hydrometeor) products from active observations by the CloudSat radar (CPR) and the CALIPSO lidar (CALIOP). Cloud products from active sensors while suffering signi􀂦cant sampling and coverage drawbacks have the advantage of resolving both horizontal and vertical variability. The main question is: given a profile of cloud condensate mean and variance, can we create a subgrid cloud field that is statistically similar to the observed subgrid cloud field? By “statistically”, we suggest that we do not aspire to reproduce “well” each individual subgrid cloud field of a GCM-scale region, but that our generator performs well for a large ensemble of cases. We simulate radar, passive imager, and radiation flux fields from the observed 2D cloud fields and create one-point statistics; we use the profiles of cloud condensate mean and variance as input to the generator to create subgrid cloud fields; we compare the statistics; we adjust the rules of the generator to create the best possible agreement between the radar, imager and radiation field statistics. In this process, the active observations have actually a dual role: they provide the actual subgrid cloud field which can be reduced to a profile of mean and variance used by the subcolumn generator, but they also provide the rules needed by the generator, such as measures (e.g, decorrelation length) of the vertical overlap of cloud fraction and of the condensate horizontal variability. This presentation will show our progress using Cloudsat and CALIPSO products in this dual fashion.

clouds↗

NASA NTRSDate not supplied

Extreme Precipitation in the Southern US Great Plains in the Spring of 2015: Mechanisms and Prediction

During May of 2015, the southern US Great Plains and adjacent Gulf Coast region experienced more than twice the long-term mean precipitation, making it the wettest May since 1895. We investigate the physical mechanisms associated with this event using a suite of large-ensemble regional replay AGCM simulations from the NASA-GEOS model. In these simulations, certain regions of the globe are constrained to closely follow observations while the remainder of the domain is free running, allowing for the isolation of the remote regions that were important for the event. Preliminary analysis provides evidence that the extreme southern US precipitation was linked in part to positive precipitation anomalies in the central and eastern tropical Pacific via a wave train, which ultimately caused anomalous moisture flux from the Gulf of Mexico. An analysis of Subseasonal Experiment (SubX) model output was conducted to explore the subseasonal prediction skill of the event. Several models are able the predict the presence of positive precipitation anomalies in or near the southern US at lead times exceeding 10 days, albeit with errors in the locations and magnitude of the heaviest precipitation anomalies. A more thorough investigation with version 2 of NASA’s GEOS-S2S model shows that the prediction skill stems from the model’s ability to reasonably predict the positive tropical Pacific precipitation anomalies and the initiation of the Rossby wave train that is believed to be linked to the event. The potential causes for limitations in the prediction skill of this event will be explored.

Great Plains↗

DOE OSTI2024

Structural mechanisms for binding and activation of a contact-quenched fluorophore by RhoBAST

The fluorescent light-up aptamer RhoBAST, which binds and activates the fluorophore–quencher conjugate tetramethylrhodamine-dinitroaniline with high affinity, super high brightness, remarkable photostability, and fast exchange kinetics, exhibits excellent performance in super-resolution RNA imaging. Here we determine the co-crystal structure of RhoBAST in complex with tetramethylrhodamine-dinitroaniline to elucidate the molecular basis for ligand binding and fluorescence activation. The structure exhibits an asymmetric “A”-like architecture for RhoBAST with a semi-open binding pocket harboring the xanthene of tetramethylrhodamine at the tip, while the dinitroaniline quencher stacks over the phenyl of tetramethylrhodamine instead of being fully released. Molecular dynamics simulations show highly heterogeneous conformational ensembles with the contact-but-unstacked fluorophore–quencher conformation for both free and bound tetramethylrhodamine-dinitroaniline being predominant. The simulations also show that, upon RNA binding, the fraction of xanthene-dinitroaniline stacked conformation significantly decreases in free tetramethylrhodamine-dinitroaniline. This highlights the importance of releasing dinitroaniline from xanthene tetramethylrhodamine to unquench the RhoBAST–tetramethylrhodamine-dinitroaniline complex. Using SAXS and ITC, we characterized the magnesium dependency of the folding and binding mode of RhoBAST in solution and indicated its strong structural robustness. The structures and binding modes of relevant fluorescent light-up aptamers are compared, providing mechanistic insights for rational design and optimization of this important fluorescent light-up aptamer-ligand system.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS2022

The Impact of Winds on AMOC in a Fully-Coupled Climate Model

Here we investigate the role of the atmospheric circulation in AMOC by comparing a fully-coupled large ensemble, a forced-ocean simulation, and new experiments using a fully-coupled global climate model where winds above the boundary layer are nudged towards reanalysis. When winds are nudged north of 45°N, agreement with RAPID array observations of AMOC at 26.5°N improves across several metrics. The phasing of interannual variability is well-captured due to the response of the local Ekman component in both wind-nudging and forced-ocean simulations, however the variance remains underestimated. The mean AMOC strength is substantially reduced relative to the fully-coupled model large ensemble, which is biased high, due to the impact of winds on surface buoyancy fluxes over the subpolar gyre. Nudging winds towards observations also reduces the 1979-2016 trend in AMOC, suggesting that improvement in the representation of the high-latitude atmosphere is important for projecting long-term AMOC changes.

AMOC↗

DOE OSTI2020

Reactive Molecular Dynamics Simulations and Quantum Chemistry Calculations To Investigate Soot-Relevant Reaction Pathways for Hexylamine Isomers

Sooting tendencies of a series of nitrogen-containing hydrocarbons (NHCs) have been recently characterized experimentally using the yield sooting index (YSI) methodology. This work aims to identify soot-relevant reaction pathways for three selected C6H15N amines, namely, dipropylamine (DPA), diisopropylamine (DIPA), and 3,3-dimethylbutylamine (DMBA) using ReaxFF molecular dynamics (MD) simulations and quantum mechanical (QM) calculations and to interpret the experimentally observed trends. ReaxFF MD simulations are performed to determine the important intermediate species and radicals involved in the fuel decomposition and soot formation processes. QM calculations are employed to extensively search for chemical reactions involving these species and radicals based on the ReaxFF MD results and also to quantitatively characterize the potential energy surfaces. Specifically, ReaxFF simulations are carried out in the NVT ensemble at 1400, 1600, and 1800 K, where soot has been identified to form in the YSI experiment. These simulations account for the interactions among test fuel molecules and pre-existing radicals and intermediate species generated from rich methane combustion, using a recently proposed simulation framework. ReaxFF simulations predict that the reactivity of the amines decrease in the order DIPA > DPA > DMBA, independent of temperature. Both QM calculations and ReaxFF simulations predict that C2H4, C3H6, and C4H8 are the main nonaromatic soot precursors formed during the decomposition of DPA, DIPA, and DMBA, respectively, and the associated reaction pathways are identified for each amine. Both theoretical methods predict that sooting tendency increases in the order DPA, DIPA, and DMBA, consistent with the experimentally measured trend in YSI. This work demonstrates that sooting tendencies and soot-relevant reaction pathways of fuels with unknown chemical kinetics can be identified efficiently through combined ReaxFF and QM simulations. Overall, predictions from ReaxFF simulations and QM calculations are consistent, in terms of fuel reactivity, major intermediates, and major nonaromatic soot precursors.

fuel decomposition↗

DOE OSTI2022

Dark Energy Survey Year 3 results: A 2.7% measurement of baryon acoustic oscillation distance scale at redshift 0.835

We present angular diameter measurements obtained by measuring the position of baryon acoustic oscillations (BAO) in an optimized sample of galaxies from the first three years of Dark Energy Survey data (DES Y3). The sample consists of 7 million galaxies distributed over a footprint of 4100 deg 2 with 0.6 < z photo < 1.1 and a typical redshift uncertainty of 0.03 ( 1 + z ) . The sample selection is the same as in the BAO measurement with the first year of DES data, but the analysis presented here uses three times the area, extends to higher redshift, and makes a number of improvements, including a fully analytical BAO template, the use of covariances from both theory and simulations, and an extensive preunblinding protocol. We used two different statistics; angular correlation function and power spectrum, and validate our pipeline with an ensemble of over 1500 realistic simulations. Both statistics yield compatible results. We combine the likelihoods derived from angular correlations and spherical harmonics to constrain the ratio of comoving angular diameter distance D M at the effective redshift of our sample to the sound horizon scale at the drag epoch. We obtain D M ( z eff = 0.835 ) / r d = 18.92 ± 0.51 , which is consistent with, but smaller than, the Planck prediction assuming flat Λ CDM , at the level of 2.3 σ . The analysis was performed blind and is robust to changes in a number of analysis choices. It represents the most precise BAO distance measurement from imaging data to date, and is competitive with the latest transverse ones from spectroscopic samples at z > 0.75 . When combined with DES 3 x 2 pt + SNIa , they lead to improvements in H 0 and Ω m constraints by ∼ 20 % .

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2015

Rising Temperatures Reduce Global Wheat Production

Crop models are essential tools for assessing the threat of climate change to local and global food production. Present models used to predict wheat grain yield are highly uncertain when simulating how crops respond to temperature. Here we systematically tested 30 different wheat crop models of the Agricultural Model Intercomparison and Improvement Project against field experiments in which growing season mean temperatures ranged from 15 degrees C to 32◦ degrees C, including experiments with artificial heating. Many models simulated yields well, but were less accurate at higher temperatures. The model ensemble median was consistently more accurate in simulating the crop temperature response than any single model, regardless of the input information used. Extrapolating the model ensemble temperature response indicates that warming is already slowing yield gains at a majority of wheat-growing locations. Global wheat production is estimated to fall by 6% for each degree C of further temperature increase and become more variable over space and time.

agriculture↗

NASA NTRS2020

Large-scale features and evaluation of the PMIP4-CMIP6 midHolocene simulations

The mid-Holocene (6000 years ago) is a standard time period for the evaluation of the simulated response of global climate models using palaeoclimate reconstructions. The latest mid-Holocene simulations are a palaeoclimate entry card for the Palaeoclimate Model Intercomparison Project (PMIP4) component of the current phase of the Coupled Model Intercomparison Project (CMIP6) – hereafter referred to as PMIP4-CMIP6. Here we provide an initial analysis and evaluation of the results of the experiment for the mid-Holocene. We show that state-of-the-art models produce climate changes that are broadly consistent with theory and observations, including increased summer warming of the Northern Hemisphere and associated shifts in tropical rainfall. Many features of the PMIP4-CMIP6 simulations were present in the previous generation (PMIP3-CMIP5) of simulations. The PMIP4-CMIP6 ensemble for the mid-Holocene has a global mean temperature change of −0.3 K, which is −0.2 K cooler than the PMIP3-CMIP5 simulations predominantly as a result of the prescription of realistic greenhouse gas concentrations in PMIP4-CMIP6. Biases in the magnitude and the sign of regional responses identified in PMIP3-CMIP5, such as the amplification of the northern African monsoon, precipitation changes over Europe, and simulated aridity in mid-Eurasia, are still present in the PMIP4-CMIP6 simulations. Despite these issues, PMIP4-CMIP6 and the mid-Holocene provide an opportunity both for quantitative evaluation and derivation of emergent constraints on the hydrological cycle, feedback strength, and potentially climate sensitivity.

mid-Holocene simulations↗

DOE OSTI2026

CONUS-wide Projected Flood Frequency and Uncertainty Estimates, Version 1.0

This dataset presents a large-ensemble of CONUS-wide projected flood frequency and uncertainty estimates across ~2.7 million NHDPlusV2 river reaches over the CONUS. The framework producing this dataset leverages a multi-model, uncertainty-aware modeling framework that allows evaluating shifts in flood frequences at the stream reach level across the CONUS. CONUS-wide ensemble streamflow projections generated from hydrologic simulations driven by downscaled and bias-corrected Coupled Model Intercomparison Project Phase 6 (CMIP6) outputs are used to derive these flood frequency and uncertainty estimates over the period 1980 - 2099. A spatially consistent regional L-moment algorithm is applied across clusters defined by the US Hydrologic Unit Code Subregions (HUC4s and HUC8s) and NHDPlusV2 stream orders to estimate flood frequencies. The dataset also includes at-site based flood estimates that allow for the comparison between local and regional approach-based estimates, assess projected changes, and characterize their uncertainties. For more reliable estimation of rare flood frequencies such as 500 and 1000-year return periods, super-ensemble based estimates are also included in the dataset. This dataset is derived to support the "Impact-Informed Dam Safety Risk Assessment for Securing Hydropower Assests" project for the US Department of Energy (DOE) Hydropower and Hydrokinetic Office (H2O). For further details, refer to Kao et al. (2022), Ghimire et al. (2023), Ghimire et al. (2025), and Hosking and Wallis (1997).

Ghimire, Ganesh [ORNL] (ORCID:0000000242843941)↗

DOE OSTI2022

Dynamical tuning of the chemical potential to achieve a target particle number in grand canonical Monte Carlo simulations

We present a method to facilitate Monte Carlo simulations in the grand canonical ensemble given a target mean particle number. The method imposes a fictitious dynamics on the chemical potential, to be run concurrently with the Monte Carlo sampling of the physical system. Corrections to the chemical potential are made according to time-averaged estimates of the mean and variance of the particle number, with the latter being proportional to thermodynamic compressibility. Here, we perform a variety of tests, and in all cases find rapid convergence of the chemical potential—inexactness of the tuning algorithm contributes only a minor part of the total measurement error for realistic simulations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Simulation of evapotranspiration and yield of maize: An inter-comparison among 41 maize models

Accurate simulation of crop water use (evapotranspiration, ET) can help crop growth models to assess the likely effects of climate change on future crop productivity, as well as being an aid for irrigation scheduling for today's growers. To determine how well maize (Zea mays L.) growth models can simulate ET, an initial inter-comparison study was conducted in 2019 under the umbrella of AgMIP (Agricultural Model Inter-Comparison and Improvement Project). Herein, we present results of a second inter-comparison study of 41 maize models that was conducted using more comprehensive datasets from two additional sites - Mead, Nebraska, USA and Bushland, Texas, USA. There were 20 treatment-years with varying irrigation levels over multiple seasons at both sites. ET was measured using eddy covariance at Mead and using large weighing lysimeters at Bushland. A wide range in ET rates was simulated among the models, yet several generally were able to simulate ET rates adequately. The ensemble median values were generally close to the observations, but a few of the models sometimes performed better than the median. Many of the models that did well at simulating ET for the Mead site did poorly for drier, windy days at the Bushland site, suggesting they need to improve how they handle humidity and wind. Additional variability came from the approaches used to simulate soil water evaporation. Fortunately, several models were identified that did well at simulating soil water evaporation, canopy transpiration, biomass accumulation, and grain yield. These models were older and have been widely used, which suggests that a larger number of users have tested these models over a wider range of conditions leading to their improvement. These revelations of the better approaches are leading to model improvements and more accurate simulations of ET.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2021

Simulated effects of sample size and grain neighborhood on the modeling of extreme value fatigue response

Assessing the size of representative volume elements (RVEs) for fatigue-related applications is challenging. A RVE relevant to random microstructure requires a volume of material that is sufficiently large to capture the grain/phase heterogeneity that captures all statistical moments of the distribution of the driving force for fatigue crack formation at “hot spot” grains. Consequently, the large size of a microstructure RVE required to study fatigue phenomena is largely computationally intractable and difficult to explore. A more realistic objective in this work is to systematically study, as a function of the size of a statistical sample of microstructure, trends towards convergence of the simulated distribution of driving force for fatigue crack formation. Our present work accordingly leverages the recently developed open-source PRISMS-Fatigue framework to examine the trends in convergence of extreme value distributions (EVD) of Fatigue Indicator Parameters (FIPs) in progressively larger polycrystalline microstructure realizations of FCC Al alloy 7075-T6 using crystal plasticity finite element method simulations. The results are compared to the traditional method in which ensembles of statistical volume elements (SVEs) are simulated to build up statistics intended to approximate those associated with a larger volume of material. The convergence of EVDs with increase of size of a SVE of microstructure is closely related to the extent of grain nearest neighbor (NN) interactions. Accordingly, the sensitivity of the local micromechanical response at hot spot grains is quantitatively investigated by systematically varying the orientations of NN grains. Results indicate that SVEs with cubic crystallographic texture tend towards convergence of the EVD of FIPs with tens of thousands of grains while the random and rolled textures require larger volumes. Simple relationships based on microstructure parameters (e.g., Schmid Factor, grain size, NN misorientation) do not completely correlate to fatigue hot spot grains. Finally, the sensitivity of the extreme value fatigue response at hot spot grains extends to the 3rd NN when a single neighborhood grain orientation is altered.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Advancing Organized Convection Representation in the Unified Model: Implementing and Enhancing Multiscale Coherent Structure Parameterization

To address the effect of stratiform latent heating on meso- to large-scale circulations, an enhanced implementation of the Multiscale Coherent Structure Parameterization (MCSP) is developed for the Met Office Unified Model. MCSP represents the top-heavy stratiform latent heating from under-resolved organized convection in general circulation models. We couple the MCSP with a mass-flux convection scheme (CoMorph-A) to improve storm lifecycle continuity. The improved MCSP trigger is specifically designed for mixed-phase deep convective cloud, combined with a background vertical wind shear, both known to be crucial for stratiform development. We also test a cloud top temperature dependent convective-stratiform heating partitioning, in contrast to the earlier fixed partitioning. Assessments from ensemble weather forecasts and decadal simulations demonstrate that MCSP directly reduces cloud deepening and precipitation areas by moderating mesoscale circulations. Indirectly, it amends tropical precipitation biases, notably correcting dry and wet biases over India and the Indian Ocean, respectively. Remarkably, the scheme outperforms a climate model ensemble by improving seasonal precipitation cycle predictions in these regions. The scheme also improves Madden-Julian Oscillation (MJO) spectra, achieving better alignment with observational and reanalysis data by intensifying the simulated MJO over the Indian Ocean during phases 4 to 5. However, the scheme increases precipitation overestimation over the Western Pacific. Shifting from fixed to temperature-dependent convective-stratiform partitioning reduces the Pacific precipitation overestimation and further improves the seasonal cycle in India. Spatially correlated biases highlight the necessity for advances beyond deterministic approaches to align MCSP with environmental conditions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Molecular simulation data for 'Data-guided Multi-Map variables for ensemble refinement of molecular movies'

These trajectories, scripts, and analysis performed on Summit underly the work published as 'Data-guided Multi-Map variables for ensemble refinement of molecular movies'. The trajectories include equilibrium and non-equilibrium sampling of ADK, CODH, and FLPP3, the scripts used to build the systems, and the scripts used to analyze the output. The directory structure is explained further in an internal README file.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

Investigating the linkage between simulated precipitation climatology and ENSO-related precipitation anomaly based on multi-model and perturbed parameter ensembles

It remains a challenge for climate models to correctly capture the relationship between precipitation and ENSO. This studyexamines the linkage between the simulated precipitation climatology and ENSO-related precipitation anomaly during borealwinter based on the multi-model ensemble from the Atmospheric Model Intercomparison Project Phase 5 (AMIP5) and perturbedparameter ensemble (PPE) with the Beijing Climate Center (BCC) atmospheric model. The AMIP5 models have large biases insimulating the tropical precipitation anomaly during El Niño, such as the shifts of the Inter-Tropical Convergence Zone (ITCZ)and South Pacific Convergence Zone (SPCZ). The inter-model differences show that the precipitation change in response to seasurface temperature (SST) change increases with enhanced precipitation climatology. The ENSO-related precipitation anomalycan also be related to the spatial distribution of the mean-state precipitation. The simulated ITCZ/SPCZ displacements aresignificantly correlated with the spatial precipitation-SST relationship in the mean state. Models with stronger mean-stateprecipitation-SST relationship also produce stronger SPCZ/ITCZ displacements. The mean-state tropical precipitation also hasstrong impacts on the ENSO-related precipitation anomalies over East Asia. In the BCC PPE, the connection between the mean-state precipitation and ENSO-related precipitation anomaly is overall consistent with that in AMIP5. Parameters associated withthe low-cloud and deep convection processes are the most influential ones for the precipitation simulations in BCC. Comparedwith the version in AMIP5, the new BCC model can better simulate the precipitation climatology and the relationship betweenENSO and precipitation over southern China. These results have important implications for model development and model-basedclimate predictions.

Yang, Ben↗