scholarly journals Emulation of long-term changes in global climate: Application to the late Pliocene and future

2017 ◽  
Author(s):  
Natalie S. Lord ◽  
Michel Crucifix ◽  
Dan J. Lunt ◽  
Mike C. Thorne ◽  
Nabila Bounceur ◽  
...  

Abstract. Multi-millennial transient simulations of climate changes have a range of important applications, such as for investigating key geologic events and transitions for which high resolution palaeoenvironmental proxy data are available, or for projecting the long-term impacts of future climate evolution on the performance of geological repositories for the disposal of radioactive wastes. However, due to the high computational requirements of current fully coupled General Circulation Models (GCMs), long-term simulations can generally only be performed with less complex models and/or at lower spatial resolution. In this study, we present novel long-term “continuous” projections of climate evolution based on the output from GCMs, via the use of a statistical emulator. The emulator is calibrated using ensembles of GCM simulations which have varying orbital configurations and atmospheric CO2 concentrations and enables a variety of investigations of long-term climate change to be conducted which would not be possible with other modelling techniques at the same temporal and spatial scales. To illustrate the potential applications, we apply the emulator to the late Pliocene (by modelling SAT), comparing its results with palaeo-proxy data for a number of global sites, and to the next 200 thousand years (kyr) (by modelling SAT and precipitation). A range of CO2 scenarios are modelled for each period. During the late Pliocene, we find that emulated SAT varies on an approximately precessional timescale, with evidence of increased obliquity response at times. A comparison of atmospheric CO2 concentration for this period, estimated using the proxy data and emulator results and using proxy CO2 records, finds that relatively similar concentrations are produced at lower latitudes, although higher latitude sites show larger discrepancies. In our second illustrative application, spanning the next 200 kyr into the future, we find that SAT oscillations appear to be primarily influenced by obliquity for the first ~ 120 kyr, whilst eccentricity is relatively low, after which precession plays a more dominant role. Conversely, variations in precipitation over the entire period demonstrate a strong precessional signal. Overall, we find that the emulator provides a useful and powerful tool for rapidly simulating the long-term evolution of climate, both past and future, due to its relatively high spatial resolution and relatively low computational cost.

2017 ◽  
Vol 13 (11) ◽  
pp. 1539-1571 ◽  
Author(s):  
Natalie S. Lord ◽  
Michel Crucifix ◽  
Dan J. Lunt ◽  
Mike C. Thorne ◽  
Nabila Bounceur ◽  
...  

Abstract. Multi-millennial transient simulations of climate changes have a range of important applications, such as for investigating key geologic events and transitions for which high-resolution palaeoenvironmental proxy data are available, or for projecting the long-term impacts of future climate evolution on the performance of geological repositories for the disposal of radioactive wastes. However, due to the high computational requirements of current fully coupled general circulation models (GCMs), long-term simulations can generally only be performed with less complex models and/or at lower spatial resolution. In this study, we present novel long-term continuous projections of climate evolution based on the output from GCMs, via the use of a statistical emulator. The emulator is calibrated using ensembles of GCM simulations, which have varying orbital configurations and atmospheric CO2 concentrations and enables a variety of investigations of long-term climate change to be conducted, which would not be possible with other modelling techniques on the same temporal and spatial scales. To illustrate the potential applications, we apply the emulator to the late Pliocene (by modelling surface air temperature – SAT), comparing its results with palaeo-proxy data for a number of global sites, and to the next 200 kyr (thousand years) (by modelling SAT and precipitation). A range of CO2 scenarios are prescribed for each period. During the late Pliocene, we find that emulated SAT varies on an approximately precessional timescale, with evidence of increased obliquity response at times. A comparison of atmospheric CO2 concentration for this period, estimated using the proxy sea surface temperature (SST) data from different sites and emulator results, finds that relatively similar CO2 concentrations are estimated based on sites at lower latitudes, whereas higher-latitude sites show larger discrepancies. In our second illustrative application, spanning the next 200 kyr into the future, we find that SAT oscillations appear to be primarily influenced by obliquity for the first ∼ 120 kyr, whilst eccentricity is relatively low, after which precession plays a more dominant role. Conversely, variations in precipitation over the entire period demonstrate a strong precessional signal. Overall, we find that the emulator provides a useful and powerful tool for rapidly simulating the long-term evolution of climate, both past and future, due to its relatively high spatial resolution and relatively low computational cost. However, there are uncertainties associated with the approach used, including the inability of the emulator to capture deviations from a quasi-stationary response to the forcing, such as transient adjustments of the deep-ocean temperature and circulation, in addition to its limited range of fixed ice sheet configurations and its requirement for prescribed atmospheric CO2 concentrations.


2014 ◽  
Vol 7 (2) ◽  
pp. 2173-2216
Author(s):  
H. Wan ◽  
P. J. Rasch ◽  
K. Zhang ◽  
Y. Qian ◽  
H. Yan ◽  
...  

Abstract. This paper explores the feasibility of an experimentation strategy for investigating sensitivities in fast components of atmospheric general circulation models. The basic idea is to replace the traditional serial-in-time long-term climate integrations by representative ensembles of shorter simulations. The key advantage of the proposed method lies in its efficiency: since fewer days of simulation are needed, the computational cost is less, and because individual realizations are independent and can be integrated simultaneously, the new dimension of parallelism can dramatically reduce the turnaround time in benchmark tests, sensitivities studies, and model tuning exercises. The strategy is not appropriate for exploring sensitivity of all model features, but it is very effective in many situations. Two examples are presented using the Community Atmosphere Model version 5. The first example demonstrates that the method is capable of characterizing the model cloud and precipitation sensitivity to time step length. A nudging technique is also applied to an additional set of simulations to help understand the contribution of physics-dynamics interaction to the detected time step sensitivity. In the second example, multiple empirical parameters related to cloud microphysics and aerosol lifecycle are perturbed simultaneously in order to explore which parameters have the largest impact on the simulated global mean top-of-atmosphere radiation balance. Results show that in both examples, short ensembles are able to correctly reproduce the main signals of model sensitivities revealed by traditional long-term climate simulations for fast processes in the climate system. The efficiency of the ensemble method makes it particularly useful for the development of high-resolution, costly and complex climate models.


2021 ◽  
Author(s):  
Gunter Stober ◽  
Ales Kuchar ◽  
Dimitry Pokhotelov ◽  
Huixin Liu ◽  
Han-Li Liu ◽  
...  

Abstract. Long-term and continuous observations of mesospheric/lower thermospheric winds are rare, but they are important to investigate climatological changes at these altitudes on time scales of several years, covering a solar cycle and longer. Such long time series are a natural heritage of the mesosphere/lower thermosphere climate, and they are valuable to compare climate models or long term runs of general circulation models (GCMs). Here we present a climatological comparison of wind observations from six meteor radars at two conjugate latitudes to validate the corresponding mean winds and atmospheric diurnal and semidiurnal tides from three GCMs, namely Ground-to-Topside Model of Atmosphere and Ionosphere for Aeronomy (GAIA), Whole Atmosphere Community Climate Model Extension (Specified Dynamics) (WACCM-X(SD)) and Upper Atmosphere ICOsahedral Non-hydrostatic (UA-ICON) model. Our results indicate that there are interhemispheric differences in the seasonal characteristics of the diurnal and semidiurnal tide. There also are some differences in the mean wind climatologies of the models and the observations. Our results indicate that GAIA shows a reasonable agreement with the meteor radar observations during the winter season, whereas WACCM-X(SD) shows a better agreement with the radars for the hemispheric zonal summer wind reversal, which is more consistent with the meteor radar observations. The free running UA-ICON tends to show similar winds and tides compared to WACCM-X(SD).


2021 ◽  
Author(s):  
Peter Nooteboom ◽  
Michiel Baatsen ◽  
Peter Bijl ◽  
Erik van Sebille ◽  
Appy Sluijs ◽  
...  

<p>Simulations of the geological past using General Circulation Models (GCMs) are computationally expensive. Mainly because of the long equilibration time scales, most of these GCMs have ocean components with a horizontal resolution of 1° or coarser. Such models are non-eddying and the effects of mesoscale ocean eddies on the transport of heat and salt are parameterized. However, from present-day ocean modeling studies, it is known that eddying ocean models better represent regional and time-mean ocean flows compared to non-eddying models. At the same time, proxy data from sediment sample sites represent climate at specific locations. Hence, the coarse ocean resolution of typical palaeo-GCMs lead to a challenge for model-data comparison in past climates.</p><p>Here we present the first simulations of a global eddying Eocene ocean with a 0.1° (horizontal) resolution model, which are initialized and forced with data from a coarser resolution (1° horizontally) equilibrated coupled ocean-atmosphere GCM. We investigate the response of the model equilibrium state to the change in ocean resolution and the consequences this has for model-data comparison in the middle-late Eocene (38Ma). We find that, compared to the non-eddying model, the eddying ocean resolution of palaeomodels reduce the biases in both sea surface temperatures and biogeographic patterns which are derived from proxy data.</p>


2019 ◽  
Author(s):  
Allison C. Michaelis ◽  
Gary M. Lackmann ◽  
Walter A. Robinson

Abstract. We present multi-seasonal simulations representative of present-day and future thermodynamic environments using the global Model for Prediction Across Scales-Atmosphere (MPAS) version 5.1 with high resolution (15 km) throughout the Northern Hemisphere. We select ten simulation years with varying phases of El Niño-Southern Oscillation (ENSO) and integrate each for 14.5 months. We use analysed sea surface temperature (SST) patterns for present-day simulations. For the future climate simulations, we alter present-day SSTs by applying monthly-averaged temperature changes derived from a 20-member ensemble of Coupled Model Intercomparison Project Phase 5 (CMIP5) general circulation models (GCMs) following the Representative Concentration Pathway (RCP) 8.5 emissions scenario. Daily sea ice fields, obtained from the monthly-averaged CMIP5 ensemble mean sea ice, are used for present-day and future simulations. The present-day simulations provide a reasonable reproduction of large-scale atmospheric features in the Northern Hemisphere such as the wintertime midlatitude storm tracks, upper-tropospheric jets, and maritime sea-level pressure features as well as annual precipitation patterns across the tropics. The simulations also adequately represent tropical cyclone (TC) characteristics such as strength, spatial distribution, and seasonal cycles for most of Northern Hemispheric basins. These results demonstrate the applicability of these model simulations for future studies examining climate change effects on various Northern Hemispheric phenomena, and, more generally, the utility of MPAS for studying climate change at spatial scales generally unachievable in GCMs.


2019 ◽  
Vol 12 (8) ◽  
pp. 3725-3743 ◽  
Author(s):  
Allison C. Michaelis ◽  
Gary M. Lackmann ◽  
Walter A. Robinson

Abstract. We present multi-seasonal simulations representative of present-day and future environments using the global Model for Prediction Across Scales – Atmosphere (MPAS-A) version 5.1 with high resolution (15 km) throughout the Northern Hemisphere. We select 10 simulation years with varying phases of El Niño–Southern Oscillation (ENSO) and integrate each for 14.5 months. We use analyzed sea surface temperature (SST) patterns for present-day simulations. For the future climate simulations, we alter present-day SSTs by applying monthly-averaged temperature changes derived from a 20-member ensemble of Coupled Model Intercomparison Project phase 5 (CMIP5) general circulation models (GCMs) following the Representative Concentration Pathway (RCP) 8.5 emissions scenario. Daily sea ice fields, obtained from the monthly-averaged CMIP5 ensemble mean sea ice, are used for present-day and future simulations. The present-day simulations provide a reasonable reproduction of large-scale atmospheric features in the Northern Hemisphere such as the wintertime midlatitude storm tracks, upper-tropospheric jets, and maritime sea-level pressure features as well as annual precipitation patterns across the tropics. The simulations also adequately represent tropical cyclone (TC) characteristics such as strength, spatial distribution, and seasonal cycles for most Northern Hemisphere basins. These results demonstrate the applicability of these model simulations for future studies examining climate change effects on various Northern Hemisphere phenomena, and, more generally, the utility of MPAS-A for studying climate change at spatial scales generally unachievable in GCMs.


2016 ◽  
Vol 33 (2) ◽  
pp. 303-311 ◽  
Author(s):  
N. C. Privé ◽  
R. M. Errico

AbstractGeneral circulation models can now be run at very high spatial resolutions to capture finescale features, but saving the full-spatial-resolution output at every model time step is usually not practical because of storage limitations. To reduce storage requirements, the model output may be produced at reduced temporal and/or spatial resolutions. When this reduced-resolution output is then used in situations where spatiotemporal interpolation is required, such as the generation of synthetic observations for observing system simulation experiments, interpolation errors can significantly affect the quality and usefulness of the reduced-resolution model output. Although it is common in practice to record model output at the highest possible spatial resolution with relatively infrequent temporal output, this may not be the best option to minimize interpolation errors. In this study, two examples using a high-resolution global run of the Goddard Earth Observing System Model, version 5 (GEOS-5), are presented to illustrate cases in which the optimal output dataset configurations for interpolation have high temporal frequency but reduced spatial resolutions. Interpolation errors of tropospheric temperature, specific humidity, and wind fields are investigated. The relationship between spatial and temporal output resolutions and interpolation errors is also characterized for the example model.


2002 ◽  
Vol 713 ◽  
Author(s):  
Mostafa Fayek ◽  
Keld A. Jensen ◽  
Rodney C. Ewing ◽  
Lee R. Riciputi

ABSTRACTUranium deposits can provide important information on the long-term performance of radioactive waste forms because uraninite (UO2+X) is similar to the UO2 in spent nuclear fuel. The Oklo-Okélobondo U-deposits, Gabon, serve as natural laboratory where the long-term (hundreds to billions of years) migration of uranium and other radionuclides can be studied over large spatial scales (nm to km). The natural fission reactors associated with the Oklo- Okélobondo U-deposits occur over a range of depths (100 to 400 m) and provide a unique opportunity to study the behavior of uraninite in near surface oxidizing environments versus more reducing conditions at depth. Previously, it has been difficult to constrain the timing of interaction between U-rich minerals and post-depositional fluids. These problems are magnified because uraninite is susceptible to alteration, it continuously self-anneals radiation damage, and because these processes are manifested at the nm to μm scale. Uranium, lead and oxygen isotopes can be used to study fluid-uraninite interaction, provided that the analyses are obtained on the micro-scale. Secondary ionization mass spectrometry (SIMS) permits in situ measurement of isotopic ratios with a spatial resolution on the scale of a few μm. Preliminary U-Pb results show that uraninite from all reactor zones are highly discordant with ages aaproaching the timing of fission chain reactions (1945±50 Ma) and resetting events at 1180±47 Ma and 898±46 Ma. Oxygen isotopic analyses show that uraninite from reactors that occur in near surface environments (δ18O= −14.4‰ to −8.5‰) have reacted more extensively with groundwater of meteoric origin relative to reactors located at greater depths (μ18O= −10.2‰ to −7.3‰). This study emphasizes the importance of using in situ high spatial resolution analysis techniques for natural analogue studies.


2020 ◽  
Author(s):  
Traute Crueger ◽  
Hauke Schmidt ◽  
Bjorn Stevens

<p>Under present day conditions the observations approximately show a hemispheric symmetry of the top of atmosphere (TOA)  short wave (SW) reflection despite the asymmetry of surface SW reflection. This has been confirmed by climate models. With models in an aqua planet setup, Voigt et al. (2014) found that tropical clouds largely compensate surface SW hemispheric asymmetries, however to a different degree in dependence on the convection scheme.</p><p>In this study, we question, whether there is also a hemispheric symmetry of TOA SW radiation under changed atmospheric radiation conditions. For that reason, we analyze experiments performed with a set of fully coupled general circulation models. The experiments were performed with either a) hemispheric asymmetric incoming radiation, b) increased atmospheric CO2 concentrations, c) increased atmospheric CO2 concentrations combined with increased stratospheric aerosol burden, or d) increased atmospheric CO2 concentration in conjunction with increased ocean albedo.</p><p>We show that generally, a hemispheric symmetry of TOA SW radiation does not occur. Overall, among the group of models, the hemispheric TOA SW radiation budgets are roughly similar for the distinct experiments, although the models utilyze different convection schemes.  We discuss the role of surface and atmospheric feedbacks in the different experiments, especially of tropical and extratropical clouds.</p><p>Reference:<br>Voigt, A., B. Stevens, J. Bader, and T. Mauritsen, 2014: Compensation of Hemispheric Albedo Asymmetries by Shifts of the ITCZ and Tropical Clouds. J. Climate, 27, 1029–1045, https://doi.org/10.1175/JCLI-D-13-00205.1.</p>


Author(s):  
R. Philbin ◽  
M. Jun

Abstract. This study validates the near-surface temperature and precipitation output from decadal runs of eight atmospheric ocean general circulation models (AOGCMs) against observational proxy data from the National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) reanalysis temperatures and Global Precipitation Climatology Project (GPCP) precipitation data. We model the joint distribution of these two fields with a parsimonious bivariate Matérn spatial covariance model, accounting for the two fields' spatial cross-correlation as well as their own smoothnesses. We fit output from each AOGCM (30-year seasonal averages from 1981 to 2010) to a statistical model on each of 21 land regions. Both variance and smoothness values agree for both fields over all latitude bands except southern mid-latitudes. Our results imply that temperature fields have smaller smoothness coefficients than precipitation fields, while both have decreasing smoothness coefficients with increasing latitude. Models predict fields with smaller smoothness coefficients than observational proxy data for the tropics. The estimated spatial cross-correlations of these two fields, however, are quite different for most GCMs in mid-latitudes. Model correlation estimates agree well with those for observational proxy data for Australia, at high northern latitudes across North America, Europe and Asia, as well as across the Sahara, India, and Southeast Asia, but elsewhere, little consistent agreement exists.


Sign in / Sign up

Export Citation Format

Share Document