The role of the mean state on MJO simulation in CESM2 ensemble simulation
Explore the source record for details and available documents.
SEARCH · Search NASA
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.
Explore the source record for details and available documents.
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.
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.
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.
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.
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
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.
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.
In this study we examine the links between tropical heating, the Madden Julian Oscillation (MJO)/Intraseasonal Oscillation (ISO), and the Asian monsoon. We are particularly interested in isolating the nature of the poleward propagation of the ISO/MJO in the monsoon region. We examine both observations and idealized "MJO heating" experiments employing the NASA Seasonal-Interannual Prediction Project (NSIPP) atmospheric general circulation model (AGCM). In the idealized 10-member ensemble simulations, the model is forced by climatological SST and an idealized eastward propagating heating profile that is meant to mimic the canonical heating associated with the MJO in the Indian Ocean and western Pacific. In order to understand the impact of SST on the off equatorial convection (or Rossby-wave response), a second set of 10-member ensemble simulations is carried out with the climatological SSTs shifted in time by 6-months. The observational analysis highlights the strong link between the Indian summer monsoon and the tropical ISO/MJO activity and heating. This includes the well-known meridional propagation that affects the summer monsoons of both hemispheres. The AGCM experiments with the idealized eastward propagating MJO-like heating reproduce the observed meridional propagation including the observed seasonal differences. The impact of the SSTs are to enhance the magnitude of the propagation into the summer hemispheres. The results suggest that the winter/summer differences associated with the MJO/ISO are auxiliary features that depend on the MJO's environment (basic state and boundary conditions) and are not the result of fundamental differences in the MJO itself.
Meaningful climate predictions should be accompanied by the corresponding uncertainty range. Common methods for estimating the uncertainty range are based on the spread of ensemble predictions. However, a simulation ensemble is not necessarily a proper sample of the real distribution of the climate, and therefore, the ensemble spread cannot be interpreted as the actual uncertainty. We propose a new method that links between the ensemble spread and the uncertainty without relying on any assumptions regarding the distribution of the ensemble predictions. The method is tested using CMIP5 1981-2010 decadal predictions and is shown to outperform other common methods.
A statistical ensemble of large eddy simulations (LES) is run simultaneously for the same flow. The information provided by the different large scale velocity fields is used to propose an ensemble averaged version of the dynamic model. This produces local model parameters that only depend on the statistical properties of the flow. An important property of the ensemble averaged dynamic procedure is that it does not require any spatial averaging and can thus be used in fully inhomogeneous flows. Also, the ensemble of LES's provides statistics of the large scale velocity that can be used for building new models for the subgrid-scale stress tensor. The ensemble averaged dynamic procedure has been implemented with various models for three flows: decaying isotropic turbulence, forced isotropic turbulence, and the time developing plane wake. It is found that the results are almost independent of the number of LES's in the statistical ensemble provided that the ensemble contains at least 16 realizations.
Various studies have documented the effects of modern ]day irrigation on regional and global climate, but none, to date, have considered the time ]varying impact of steadily increasing irrigation rates on climate during the 20th century. We investigate the impacts of observed irrigation changes over this century with two ensemble simulations using an atmosphere general circulation model. Both ensembles are forced with transient climate forcings and observed sea surface temperatures from 1902 to 2000; one ensemble includes irrigation specified by a time ]varying data set of irrigation water withdrawals. Early in the century, irrigation is primarily localized over southern and eastern Asia, leading to significant cooling in boreal summer (June.August) over these regions. This cooling spreads and intensifies by century fs end, following the rapid expansion of irrigation over North America, Europe, and Asia. Irrigation also leads to boreal winter (December.February) warming over parts of North America and Asia in the latter part of the century, due to enhanced downward longwave fluxes from increased near ]surface humidity. Precipitation increases occur primarily downwind of the major irrigation areas, although precipitation in parts of India decreases due to a weaker summer monsoon. Irrigation begins to significantly reduce temperatures and temperature trends during boreal summer over the Northern Hemisphere midlatitudes and tropics beginning around 1950; significant increases in precipitation occur in these same latitude bands. These trends reveal the varying importance of irrigation ]climate interactions and suggest that future climate studies should account for irrigation, especially in regions with unsustainable irrigation resources.
The climate change impact and adaptation simulations from the Agricultural Model Intercomparison and Improvement Project (AgMIP) for wheat provide a unique dataset of multi-model ensemble simulations for 60 representative global locations covering all global wheat mega environments. The multi-model ensemble reported here has been thoroughly benchmarked against a large number of experimental data, including different locations, growing season temperatures, atmospheric CO2 concentration, heat stress scenarios, and their interactions. In this paper, we describe the main characteristics of this global simulation dataset. Detailed cultivar, crop management, and soil datasets were compiled for all locations to drive 32 wheat growth models. The dataset consists of 30-year simulated data including 25 output variables for nine climate scenarios, including Baseline (1980-2010) with 360 or 550 ppm CO2, Baseline +2oC or +4oC with 360 or 550 ppm CO2, a mid-century climate change scenario (RCP8.5, 571 ppm CO2), and 1.5°C (423 ppm CO2) and 2.0oC (487 ppm CO2) warming above the pre-industrial period (HAPPI). This global simulation dataset can be used as a benchmark from a well-tested multi-model ensemble in future analyses of global wheat. Also, resource use efficiency (e.g., for radiation, water, and nitrogen use) and uncertainty analyses under different climate scenarios can be explored at different scales. The DOI for the dataset is 10.5281/zenodo.4027033 (AgMIP-Wheat, 2020), and all the data are available on the data repository of Zenodo (http://doi.org/10.5281/zenodo.4027033). Two scientific publications have been published based on some of these data here.
Five annual climate cycles (1998-2002) are simulated for continental Africa and adjacent oceans by a regional atmospheric model (RM3). RM3 horizontal grid spacing is 0.44deg at 28 vertical levels. Each of 2 simulation ensembles is driven by lateral boundary conditions from each of 2 alternative reanalysis data sets. One simulation downs cales National Center for Environmental Prediction reanalysis 2 (NCPR2) and the other the European Centre for Medium Range Weather Forecasts Interim reanalysis (ERA-I). NCPR2 data are archived at 2.5deg grid spacing, while a recent version of ERA-I provides data at 0.75deg spacing. ERA-I-forced simulations are recomrp. ended by the Coordinated Regional Downscaling Experiment (CORDEX). Comparisons of the 2 sets of simulations with each other and with observational evidence assess the relative performance of each downscaling system. A third simulation also uses ERA-I forcing, but degraded to the same horizontal resolution as NCPR2. RM3-simulated pentad and monthly mean precipitation data are compared to Tropical Rainfall Measuring Mission (TRMM) data, gridded at 0.5deg, and RM3-simulated circulation is compared to both reanalyses. Results suggest that each downscaling system provides advantages and disadvantages relative to the other. The RM3/NCPR2 achieves a more realistic northward advance of summer monsoon rains over West Africa, but RM3/ERA-I creates the more realistic monsoon circulation. Both systems recreate some features of JulySeptember 1999 minus 2002 precipitation differences. Degrading the resolution of ERA-I driving data unrealistically slows the monsoon circulation and considerably diminishes summer rainfall rates over West Africa. The high resolution of ERA-I data, therefore, contributes to the quality of the downscaling, but NCPR2laterai boundary conditions nevertheless produce better simulations of some features.
The Earth possesses an internal magnetic field (geomagnetic field) generated by convection in the outer core (geodynamo). Previous efforts have been focused along two distinct paths: (1) numerical geodynamo modeling to understand the origin of the geomagnetic field, and the mechanisms of geomagnetic secular variations (SV); and (2) geomagnetic field modeling to map the spatial/temporal variations of the field from geomagnetic data, and to derive core properties, e.g. inversion of core flow near the core-mantle boundary (CMB). Geomagnetic data assimilation is a new approach emerged over the past 5 years: surface observations are assimilated with geodynamo models for better understanding of the core dynamical state, and accurately prediction of SV. In collaboration with several geomagnetic research groups, we have developed the first working geomagnetic data assimilation system, Modular, Scalable, Self-consistent, and Three-dimensional (MoSST) DAS, that includes the MoSST numerical dynamo model; 7000 years of geomagnetic field maps from several field models utilizing satellite and ground observatory data, historical magnetic records and archeo/paleo magnetic data; and an ensemble based optimal interpolation (01) assimilation algorithm. With this system, we have demonstrated clearly that the assimilated core dynamical state is substantially different from those of pure geodynamo simulations. Ensemble assimilation runs also show the convergence of the assimilated solutions inside the core, suggesting that the simulation state is pulled closer to the truth via data assimilation. The forecasts from this system are also very accurate: the 5-year forecast of the geomagnetic field agrees very well with the observations; and the 5-year secular variation forecast is more accurate than the IGRF SV forecast models in the past. Using geomagnetic records up to 2009, we have made an SV forecast for the period from 2010-2015, and is a candidate SV model for IGRF-11.
Studies have shown the correlation and connection of surface temperatures across the globe, ocean and land, related to Tropical SSTs especially El Nino. This climate variability greatly influences regional weather and hydroclimate extremes (e.g. drought and flood). In this paper, we evaluate the relationship of temperatures across the tropical oceans and continents in MERRA-2, and also in a newly developed MERRA-2 AMIP ensemble simulation (M2AMIP). M2AMIP uses the same model and spatial resolution as MERRA-2, producing the same output diagnostics over 10 ensemble members. Composite El Nino temperature data are compared with observations to evaluate the land/sea contrast, variations and phase relationship. The temperature variations are related to surface heat fluxes and the atmospheric temperatures and transport, to identify the processes that lead to the lagged redistribution of heat in the tropics and beyond. Discernable cloud, radiation and data assimilation changes accompany the onset of El Nino affecting continental regions through the progression to and following the peak values. While the model represents these variations in general, regional strengths and weaknesses can be identified.
A bi-stable mode of the Atlantic Meridional Overturning Circulation (AMOC) is found in a 10-member ensemble simulation of the SSP2-4.5 scenario using the NASA GISS-E2-1-G climate model. Local feedbacks in the subpolar North Atlantic region in conjunction with internal variability in sea-ice transport and melt play a critical role in causing the divergent behavior of the AMOC in the ensemble members. While other fully coupled models have demonstrated the important role of surface freshening in leading to AMOC shutdown, either through hosing experiments or increased precipitation and greenhouse gas warming at high latitudes, in the GISS simulations, there are no external freshwater perturbations. This is the first time that a CMIP-class model has shown such a bifurcation across an initial condition ensemble.
A 10-member ensemble simulation using the NASA GISS-E2-1-G climate model shows a clear bifurcation in the Atlantic Meridional Overturning Circulation (AMOC) strength under the SSP2-4.5 future scenario, leading to 8 stable AMOC-on and 2 stable AMOC-off climate states after 400 years of integration, despite identical forcing of each ensemble member. Ocean models have previously shown a tipping point where sufficient forcing (e.g. by freshwater inputs) can shut down the AMOC. A variety of fully coupled models have demonstrated this, either through hosing experiments or increased precipitation and warming at high latitudes due to increased concentrations of greenhouse gases. In the GISS simulations, there are no external freshwater perturbations. The bifurcation arises freely in the coupled system and is the result of stochastic variability associated with sea-ice transport and melting in the Irminger Sea following a slowing of the GHG concentration increase at the end of the twenty-first century. We believe this is the first time that a CMIP-class model has shown such a bifurcation across an initial condition ensemble.