Transient Climate Simulations Research Articles

Contrasting latitudinal evolution of East Asian monsoonal precipitation during the Last Interglacial (130–120 ka)

Evolution of East Asian monsoonal precipitation across the Last Interglacial (LIG) remains controversial, owing to the discrepancies between various proxies and their low temporal resolution. Through a transient high-resolution global climate simulation covering the interval of 130–120 ka, we illustrate a long-term increasing (decreasing) trend in summer precipitation over south China (northeast Asia) during the LIG (i.e. 130–120 ka). The out-of-phase precipitation evolution across latitudes were coherently regulated by the weakened monsoonal circulation, southward moved western North Pacific high, and southward displaced East Asian westerly jet from the early to late LIG. These atmospheric circulation variations were in turn determined by sea surface temperature anomalies over the Pacific and the propagation of extratropical Rossby waves originating from North Africa. Our results may provide important insights for reconciling discrepancies between precipitation proxies during the LIG and for precipitation behavior in a warmer-than-present world.

Fundamental Shift From Summer to Winter of Holocene Rainfall Regime in the Tropics

AbstractThe seasonal rainfall regime is a key factor control on local ecological and social processes and is commonly thought to be stable under long‐term climate changes. Here we present a unique high‐resolution rainfall record from the Thai‐Malay Peninsula, combined with a state‐of‐the‐art transient climate simulation, demonstrating a fundamental rainfall regime shift from summer to winter during the Holocene. Transient model simulation and new sensitivity experiments further reveal that westward migration of the boundary between summer and winter rainfall regimes results in a summer to winter rainfall regime shift forced by distinct changes in summer and winter monsoons. Our findings suggest that the seasonal rainfall regime could be unstable under climate change around the boundaries of rainfall regimes in the tropics and possibly worldwide, which might be more critical for shaping both past and future ecological environments.

Revisiting Western United States Hydroclimate During the Last Deglaciation

AbstractDuring the last ice age, the western United States was covered by large lakes, sustained partly by higher levels of precipitation. Increased rainfall was driven by the atmospheric circulation associated with the presence of large North American ice sheets, yet Pleistocene lakes generally reached their highstands not at glacial maximum but during deglaciation. Prior modeling studies, however, showed nearly monotonic drying since the last glacial maximum. Here I show that iTraCE, a new transient climate simulation of the last deglaciation, reproduces a robust peak in winter rainfall over the Great Basin near 16 ka. The simulated peak is driven by a transient strengthening and southward shift of the midlatitude jet. While meltwater forcing is an important driver of changes to the North Pacific Jet, changing orbital conditions and rising atmospheric CO2 also shift the jet south and contribute to wetter conditions over the western US during deglaciation.

Reversals in Temperature‐Precipitation Correlations in the Northern Hemisphere Extratropics During the Holocene

AbstractFuture precipitation levels remain uncertain because climate models have struggled to reproduce observed variations in temperature‐precipitation correlations. Our analyses of Holocene proxy‐based temperature‐precipitation correlations and hydrological sensitivities from 2,237 Northern Hemisphere extratropical pollen records reveal a significant latitudinal dependence and temporal variations among the early, middle, and late Holocene. These proxy‐based variations are largely consistent with patterns obtained from transient climate simulations (TraCE21k). While high latitudes and subtropical monsoon areas show mainly stable positive correlations throughout the Holocene, the mid‐latitude pattern is temporally and spatially more variable. In particular, we identified a reversal from positive to negative temperature‐precipitation correlations in the eastern North American and European mid‐latitudes from the early to mid‐Holocene that mainly related to slowed down westerlies and a switch to moisture‐limited convection under a warm climate. Our palaeoevidence of past temperature‐precipitation correlation shifts identifies those regions where simulating past and future precipitation levels might be particularly challenging.

Warm pool ocean heat content regulates ocean-continent moisture transport.

The Indo-Pacific Warm Pool (IPWP) exerts a dominant role in global climate by releasing huge amounts of water vapour and latent heat to the atmosphere and modulating upper ocean heat content (OHC), which has been implicated in modern climate change1. The long-term variations of IPWP OHC and their effect on monsoonal hydroclimate are, however, not fully explored. Here, by combining geochemical proxies and transient climate simulations, we show that changes of IPWP upper (0-200 m) OHC over the past 360,000 years exhibit dominant precession and weaker obliquity cycles and follow changes in meridional insolation gradients, and that only 30%-40% of the deglacial increases are related to changes in ice volume. On the precessional band, higher upper OHC correlates with oxygen isotope enrichments in IPWP surface water and concomitant depletion in East Asian precipitation as recorded in Chinese speleothems. Using an isotope-enabled air-sea coupled model, we suggest that on precessional timescales, variations in IPWP upper OHC, more than surface temperature, act to amplify the ocean-continent hydrological cycle via the convergence of moisture and latent heat. From an energetic viewpoint, the coupling of upper OHC and monsoon variations, both coordinated by insolation changes on orbital timescales, is critical for regulating the global hydroclimate.

Influence of Hillslope Flow on Hydroclimatic Evolution Under Climate Change

AbstractWe analyzed the influence of hillslope flow on projections of climate change by comparing two transient climate simulations with the IPSL climate model between 1980 and 2100. Hillslope flow induces a reorganization and increment of soil moisture (+10%), which increases evapotranspiration (+4%) and precipitation (+1%) and decreases total runoff (−3%) and air temperature (−0.1 °C) on an annual average over land for 1980–2010 when compared to simulation not representing hillslope flow. These changes in land/atmosphere fluxes are not homogenous and depend on regional climate and surface conditions. Hillslope flow also influences climate change projections. On average over land, it amplifies the positive trend of soil moisture (+23%), evapotranspiration (+50%), and precipitation (+7%) and slightly attenuates global warming (−1%), especially for daily maximum air temperature. The role of hillslope flow in supporting surface/atmosphere fluxes is more evident at a regional scale. Where precipitation is projected to decrease, hillslope flow is shown to attenuate the related declines in evapotranspiration, precipitation, and total runoff, regardless of aridity conditions and mean air temperature. Where precipitation is projected to increase, hillslope flow amplifies evapotranspiration enhancement but attenuates the increase in precipitation and total runoff. Warming is generally attenuated, especially in semiarid and cold areas, and humid and warm/temperate regions, but the signal is weak. These results demonstrate the role of hillslope flow in enhancing water and energy fluxes between the surface and the atmosphere. They also suggest that including hillslope flow in climate models would weaken the projected intensification of hydrological extreme events.

Transient climate simulations of the Holocene (version 1) – experimental design and boundary conditions

Abstract. The Holocene, which started approximately 11.5 ka, is the latest interglacial period with several rapid climate changes with timescales, from decades to centuries, superimposed on the millennium-scale mean climate trend. Climate models provide useful tools to investigate the underlying dynamic mechanisms for the climate change during this well-studied time period. Thanks to the improvements in the climate model and computational power, transient simulation of the Holocene offers an opportunity to investigate the climate evolution in response to time-varying external forcings and feedbacks. Here, we present the design of a new set of transient experiments for the whole Holocene from 11.5 ka to the preindustrial period (1850; HT-11.5 ka) to investigate both the combined and separated effects of the main external forcing of orbital insolation, atmospheric greenhouse gas (GHG) concentrations, and ice sheets on the climate evolution over the Holocene. The HT-11.5 ka simulations are performed with a relatively high-resolution version of the comprehensive Earth system model CESM1.2.1 without acceleration, both fully and singly forced by time-varying boundary conditions of orbital configurations, atmospheric GHGs, and ice sheets. Preliminary simulation results show a slight decrease in the global annual mean surface air temperature from 11.5 to 7.5 ka due to both changes in orbital insolation and GHG concentrations, with an abrupt cooling at approximately 7.5 ka, which is followed by a continuous warming until the preindustrial period, mainly due to increased GHG concentrations. Both at global and zonal mean scales, the simulated annual and seasonal temperature changes at 6 ka lie within the range of the 14 PMIP4 model results and are overall stronger than their arithmetic mean results for the Middle Holocene simulations. Further analyses on the HT-11.5 ka transient simulation results will be covered by follow-up studies.

In-phase millennial-scale glacier changes in the tropics and North Atlantic regions during the Holocene

Based on new and published cosmic-ray exposure chronologies, we show that glacier extent in the tropical Andes and the north Atlantic regions (TANAR) varied in-phase on millennial timescales during the Holocene, distinct from other regions. Glaciers experienced an early Holocene maximum extent, followed by a strong mid-Holocene retreat and a re-advance in the late Holocene. We further explore the potential forcing of TANAR glacier variations using transient climate simulations. Since the Atlantic Meridional Overturning Circulation (AMOC) evolution is poorly represented in these transient simulations, we develop a semi-empirical model to estimate the “AMOC-corrected” temperature and precipitation footprint at regional scales. We show that variations in the AMOC strength during the Holocene are consistent with the observed glacier changes. Our findings highlight the need to better constrain past AMOC behavior, as it may be an important driver of TANAR glacier variations during the Holocene, superimposed on other forcing mechanisms.

Identifying Global‐Scale Patterns of Vegetation Change During the Last Deglaciation From Paleoclimate Networks

AbstractDuring the last deglaciation (∼19–11 ka before present), the global mean temperature increased by 3–8 K. The concurrent hydroclimate and land cover changes are not well constrained. Here, we use a pollen database to quantify global‐scale vegetation changes during this transitional period at orbital (∼104 years) and millennial timescales (∼103 years). We focus on the proportion of tree and shrub pollen, the arboreal pollen (AP) fraction. Temporal similarities over long distances are identified by a paleoclimate network approach. At the orbital scale, we find coherent AP variations in the low and mid‐latitudes which we attribute to the global climate forcing. While AP fractions predominantly increased through the deglaciation, we identify regions where AP fractions decreased. For millennial timescales, we do not observe spatially coherent similarity structures. We compare our results with networks computed from three deglacial climate simulations with the CCSM3, HadCM3, and LOVECLIM models. Networks based on simulated precipitation patterns reproduce the characteristics of the AP network. Sensitivity experiments with statistical emulators indicate that, indeed, precipitation variations explain the diagnosed patterns of vegetation change better than temperature and CO2 variations. Our findings support previous interpretations of deglacial forest evolution in the mid‐latitudes being the result of atmospheric circulation changes. The network analysis identifies differences in the vegetation‐climate‐CO2 relationship simulated by CCSM3 and HadCM3. We conclude that network analyses are a promising tool to benchmark transient climate simulations with dynamical vegetation changes. This may result in stronger constraints of future hydroclimate and land cover changes.

Retreat and Regrowth of the Greenland Ice Sheet During the Last Interglacial as Simulated by the CESM2‐CISM2 Coupled Climate–Ice Sheet Model

AbstractDuring the Last Interglacial, approximately 129 to 116 ka (thousand years ago), the Arctic summer climate was warmer than the present, and the Greenland Ice Sheet retreated to a smaller extent than its current state. Previous model‐derived and geological reconstruction estimates of the sea‐level contribution of the Greenland Ice Sheet during the Last Interglacial vary widely. Here, we conduct a transient climate simulation from 127 to 119 ka using the Community Earth System Model (CESM2), which includes a dynamic ice sheet component (the Community Ice Sheet Model, CISM2) that is interactively coupled to the atmosphere, land, ocean, and sea ice components. Vegetation distribution is updated every 500 years based on biomes simulated using a monthly climatology to force the BIOME4 equilibrium vegetation model. Results show a substantial retreat of the Greenland Ice Sheet, reaching a minimum extent at 121.9 ka, equivalent to a 3.0 m rise in sea level relative to the present day, followed by gradual regrowth. In contrast, a companion simulation employing static vegetation based on pre‐industrial conditions shows a much smaller ice‐sheet retreat, highlighting the importance of the changes in high‐latitude vegetation distribution for amplifying the ice‐sheet response.

Strong Southern African Monsoon and weak Mozambique Channel throughflow during Heinrich events: Implication for Agulhas leakage

The Mozambique Channel is a conduit of trade wind-driven throughflow that is a key component of the Agulhas Current and Agulhas leakage, a flux of warm and salty water from the tropical Indo-Pacific to the Atlantic Ocean. Agulhas leakage is thought to modulate Atlantic meridional overturning circulation variability. Previous studies from the Cape Basin suggest that enhanced Agulhas leakage played an important role in accelerating glacial terminations. The southern African monsoon response to abrupt climate changes associated with meltwater-induced reorganizations of the North Atlantic meridional overturning circulation, and its impact on the Mozambique Channel throughflow and, by extension, on the Agulhas leakage is not well understood. Here we present a high-resolution 26,000 year-long hydroclimate record of northern Madagascar, a core region of the southern hemisphere monsoon domain, and a mixed layer temperature reconstruction using sediment cores collected from the runoff-influenced eastern Mozambique Channel. The record indicates precipitation increases centered at 11.7-12.5 thousand years before present (kyr BP), 14.5-19 kyr BP, 23-24 kyr BP, 25-26 kyr BP. Considering age model uncertainties, this is the first strong evidence for southern African monsoon strengthening in response to meltwater-induced northern high latitude climate instabilities during the Younger Dryas (YD), Heinrich Stadial 1 (HS1), HS2 and the HS-like event prior to HS2, in agreement with the results of transient climate simulations. Furthermore, our study shows a reversal of the mixed layer temperature gradient between the western and eastern Mozambique Channel during Heinrich event 1 (HE1). We posit that the gradient reversal indicates a weakening of the trade wind-driven South Equatorial Current and Mozambique Channel throughflow that likely weakened the Agulhas leakage, potentially creating a feedback that may have contributed to the sustained weakening of the AMOC during HE1 by reducing the amount of heat and salt leakage into the Atlantic.

Influence of Arctic sea-ice loss on the Greenland ice sheet climate

The Arctic is the region on Earth that is warming the fastest. At the same time, Arctic sea ice is reducing while the Greenland ice sheet (GrIS) is losing mass at an accelerated pace. Here, we study the seasonal impact of reduced Arctic sea ice on GrIS surface mass balance (SMB), using the Community Earth System Model version 2.1 (CESM2), which features an advanced, interactive calculation of SMB. Addressing the impact of sea-ice reductions on the GrIS SMB from observations is difficult due to the short observational records. Also, signals detected using transient climate simulations may be aliases of other forcings. Here, we analyze dedicated simulations from the Polar Amplification Model Intercomparison Project with reduced Arctic sea ice and compare them with preindustrial sea ice simulations while keeping all other forcings constant. In response to reduced sea ice, the GrIS SMB increases in winter due to increased precipitation, driven by the more humid atmosphere and increasing cyclones. In summer, surface melt increases due to a warmer, more humid atmosphere providing increased energy transfer to the surface through the sensible and latent heat fluxes, which triggers the melt-albedo feedback. Further, warming occurs throughout the entire troposphere over Baffin Bay. This deep warming results in regional enhancement of the 500 hPa geopotential heights over the Baffin Bay and Greenland, increasing blocking and heat advection over the GrIS’ surface. This anomalous circulation pattern has been linked to recent increases in the surface melt of the GrIS.

First assessment of the earth heat inventory within CMIP5 historical simulations

Abstract. The energy imbalance at the top of the atmosphere over the last century has caused an accumulation of heat within the ocean, the continental subsurface, the atmosphere and the cryosphere. Although ∼90 % of the energy gained by the climate system has been stored in the ocean, the other components of the Earth heat inventory cannot be neglected due to their influence on associated climate processes dependent on heat storage, such as sea level rise and permafrost stability. However, there has not been a comprehensive assessment of the heat inventory within global climate simulations yet. Here, we explore the ability of 30 advanced general circulation models (GCMs) from the fifth phase of the Coupled Model Intercomparison Project (CMIP5) to simulate the distribution of heat within the Earth's energy reservoirs for the period 1972–2005 of the Common Era. CMIP5 GCMs simulate an average heat storage of 247±172 ZJ (96±4 % of total heat content) in the ocean, 5±9 ZJ (2±3 %) in the continental subsurface, 2±3 ZJ (1±1 %) in the cryosphere and 2±2 ZJ (1±1 %) in the atmosphere. However, the CMIP5 ensemble overestimates the ocean heat content by 83 ZJ and underestimates the continental heat storage by 9 ZJ and the cryosphere heat content by 5 ZJ, in comparison with recent observations. The representation of terrestrial ice masses and the continental subsurface, as well as the response of each model to the external forcing, should be improved in order to obtain better representations of the Earth heat inventory and the partition of heat among climate subsystems in global transient climate simulations.

Possible obliquity-forced warmth in southern Asia during the last glacial stage

Orbital-scale global climatic changes during the late Quaternary are dominated by high-latitude influenced ~100,000-year global ice-age cycles and monsoon influenced ~23,000-year low-latitude hydroclimate variations. However, the shortage of highly-resolved land temperature records remains a limiting factor for achieving a comprehensive understanding of long-term low-latitude terrestrial climatic changes. Here, we report paired mean annual air temperature (MAAT) and monsoon intensity proxy records over the past 88,000 years from Lake Tengchongqinghai in southwestern China. While summer monsoon intensity follows the ~23,000-year precession beat found also in previous studies, we identify previously unrecognized warm periods at 88,000–71,000 and 45,000–22,000 years ago, with 2–3 °C amplitudes that are close to our recorded full glacial-interglacial range. Using advanced transient climate simulations and comparing with forcing factors, we find that these warm periods in our MAAT record probably depends on local annual mean insolation, which is controlled by Earth’s ~41,000-year obliquity cycles and is anti-phased to annual mean insolation at high latitudes. The coincidence of our identified warm periods and intervals of high-frequent dated archaeological evidence highlights the importance of temperature on anatomically modern humans in Asia during the last glacial stage.

Holocene precipitation changes in northeastern China from CCSM3 transient climate simulations

To date, climate records have mainly shown three different trends of Holocene precipitation evolution in northeastern (NE) China, and the underlying mechanisms remain unclear. Here, we use model results from Holocene transient climate simulations conducted by the Community Climate System Model version 3 to investigate the evolution of precipitation in NE China and the associated mechanisms. The model results indicate that precipitation changes within NE China show obvious spatial discrepancies. In particular, the annual precipitation maximum occurs in the early Holocene for the western subregion, while it occurs in the mid-late Holocene for the eastern subregion. These results show two different trends of Holocene precipitation within NE China capturing the large-scale precipitation changes appearing in climate records. These spatial features are closely related to the gradual weakening of the East Asian summer monsoon during the Holocene and are mainly attributed to the combined effects of orbital forcing and the ice sheet. Changes in orbital parameters play a major role in the decreased precipitation in the western subregion, while changes in the ice sheet contribute more to the increased precipitation in the eastern subregion. The observed model-data discrepancy partly relates to the low horizontal resolution employed and the physical processes and parameterizations of the model used.

Overestimated acceleration of the advective Brewer–Dobson circulation due to stratospheric cooling

AbstractTropospheric warming and stratospheric cooling influence the vertical structure of the atmosphere. Numerous studies have analysed the thermal expansion of the troposphere, however, stratospheric cooling reverses the sign of this shift in the middle stratosphere, causing a downward shift in the upper stratosphere and mesosphere. This is a robust feature in transient climate model simulations, but its impact is commonly unappreciated. Here, we quantify the trend difference of the residual mean vertical velocity (), a proxy for diagnosing the advective Brewer–Dobson circulation (BDC) strength, which arises from implicit neglect of the shrinking distance between stratospheric pressure levels in the CCMI‐1 (Chemistry‐Climate Model Initiative part 1) data request. There, a log‐pressure formula with constant scale height is recommended to compute . However, stratospheric cooling in transient climate simulations causes a reduction of the geometrical distance between pressure levels and thereby also the scale height significantly decreases over time. Using the general scale height definition for the transformation, the trends are therefore smaller. In both cases, the units are m·s−1, but in the latter case it is the constant measure of length geopotential metres and not log‐pressure metres. We quantify that, due to the temperature dependence of log‐pressure metres, past studies that based trend analyses on log‐pressure overestimated the advective BDC acceleration by ∼ 20%. This result is consistent among the CCMI‐1 projections over the 1960–2100 period. We highlight that other diagnostics can also be affected by the neglect of the declining stratospheric pressure level distance. A detailed description of the diagnostics is necessary for consistent assessments of trends. Data processing tools should generally not include the constant scale height assumption if the data are used for trend analyses.

Modulation of Indian monsoon by water vapor and cloud feedback over the past 22,000 years

To predict how monsoons will evolve in the 21st century, we need to understand how they have changed in the past. In paleoclimate literature, the major focus has been on the role of solar forcing on monsoons but not on the amplification by feedbacks internal to the climate system. Here we have used the results from a transient climate simulation to show that feedbacks amplify the effect of change in insolation on the Indian summer monsoon. We show that during the deglacial (22 ka to 10 ka) monsoons were predominantly influenced by rising water vapor due to increasing sea surface temperature, whereas in the Holocene (10 ka to 0 ka) cloud feedback was more important. These results are consistent with another transient simulation, thus increasing confidence despite potential model biases. We have demonstrated that insolation drives monsoon through different pathways during cold and warm periods, thereby highlighting the changing role of internal factors.

Evolution of the Asian–African Monsoonal Precipitation over the last 21 kyr and the Associated Dynamic Mechanisms

Abstract The Asian–African monsoonal precipitation (AAMP) has a significant impact on the water availability, biodiversity, and livelihoods of billions of people. A comprehensive understanding of the AAMP behavior over Earth’s history will help to make better future projections. Using a set of transient climate simulations over the last 21 000 years (21 ka), the variation of the AAMP and its responses to various external forcings, including orbital insolation, greenhouse gases (GHGs), and ice sheets, are explored. The precipitation evolutions in the individual monsoon domains have the characteristic of hemispheric synchrony over the last 21 ka. Specifically, the AAMP increased from the Last Glacial Maximum to the early Holocene with several abrupt events and then decreased subsequently. The raised orbital insolation and GHGs lead to an overall AAMP increase, but the enhanced insolation tends to induce a systematic northward shift of the Asian–African monsoon domain. Decreased meltwater discharge could promote the African and Indian monsoonal precipitation through strengthening the Atlantic Ocean meridional overturning circulation. However, the lowering of ice sheets (i.e., orographic effect) results in an anomalous dipole precipitation pattern between North China and India. An analysis of the moisture budget suggests that, although different external forcings may lead to the same sign of precipitation change (e.g., both increased insolation and GHGs can cause the enhanced AAMP), the thermodynamic and dynamic contributions to precipitation could vary greatly by region and forcing. This study provides a reference for the long-term behavior of the AAMP with rising GHGs, higher insolation, and potential melting of the Greenland Ice Sheet.

Paleo calendar-effect adjustments in time-slice and transient climate-model simulations (PaleoCalAdjust v1.0): impact and strategies for data analysis

Abstract. The “paleo calendar effect” is a common expression for the impact that changes in the length of months or seasons over time, related to changes in the eccentricity of Earth's orbit and precession, have on the analysis or summarization of climate-model output. This effect can have significant implications for paleoclimate analyses. In particular, using a “fixed-length” definition of months (i.e., defined by a fixed number of days), as opposed to a “fixed-angular” definition (i.e., defined by a fixed number of degrees of the Earth's orbit), leads to comparisons of data from different positions along the Earth's orbit when comparing paleo with modern simulations. This effect can impart characteristic spatial patterns or signals in comparisons of time-slice simulations that otherwise might be interpreted in terms of specific paleoclimatic mechanisms, and we provide examples for 6, 97, 116, and 127 ka. The calendar effect is exacerbated in transient climate simulations in which, in addition to spatial or map-pattern effects, it can influence the apparent timing of extrema in individual time series and the characterization of phase relationships among series. We outline an approach for adjusting paleo simulations that have been summarized using a modern fixed-length definition of months and that can also be used for summarizing and comparing data archived as daily data. We describe the implementation of this approach in a set of Fortran 90 programs and modules (PaleoCalAdjust v1.0).

Varying Sensitivity of East Asia Summer Monsoon Circulation to Temperature Change Since Last Glacial Maximum

AbstractProxy records of the East Asia summer monsoon (EASM) and temperature reveal an in‐phase relationship of the two since the Last Glacial Maximum (LGM), which is consistent with their expected physical relationship; but the response amplitude of EASM to temperature variation seemed larger over the Holocene than over the last deglaciation. Using a state‐of‐the‐art transient climate simulation, we confirm this proxy‐inferred phenomenon, and estimate the “sensitivity of EASM circulation to summer temperature over adjacent region” (defined as SMT here) is about five times over the Holocene of that over the last deglaciation. This varying SMT attributed to the diverse SMT under the solo forcing of orbital insolation and CO2 (eight times of the former vs. the latter), and their varying confounding effect along time. The results presented here should improve our understanding of the past change of EASM, and help project its variation magnitude under ongoing CO2‐forced global warming.