Description: Numerical simulations employing prognostic stable water isotopes can not only facilitate our understanding of hydrological processes and climate change but also allow for a direct comparison between isotope signals obtained from models and various archives. In the current work, we describe the performance and explore the potential of a new version of the Earth system model AWI-ESM (Alfred Wegener Institute Earth System Model), labeled AWI-ESM-2.1-wiso, in which we incorporated three isotope tracers into all relevant components of the water cycle. We present here the results of pre-industrial (PI) and mid-Holocene (MH) simulations. The model reproduces the observed PI isotope compositions in both precipitation and seawater well and captures their major differences from the MH conditions. The simulated relationship between the isotope composition in precipitation (δ18Op) and surface air temperature is very similar between the PI and MH conditions, and it is largely consistent with modern observations despite some regional model biases. The ratio of the MH–PI difference in δ18Op to the MH–PI difference in surface air temperature is comparable to proxy records over Greenland and Antarctica only when summertime air temperature is considered. An amount effect is evident over the North African monsoon domain, where a negative correlation between δ18Op and the amount of precipitation is simulated. As an example of model applications, we studied the onset and withdrawal date of the MH West African summer monsoon (WASM) using daily variables. We find that defining the WASM onset based on precipitation alone may yield erroneous results due to the substantial daily variations in precipitation, which may obscure the distinction between pre-monsoon and monsoon seasons. Combining precipitation and isotope indicators, we suggest in this work a novel method for identifying the commencement of the WASM. Moreover, we do not find an obvious difference between the MH and PI periods in terms of the mean onset of the WASM. However, an advancement in the WASM withdrawal is found in the MH compared to the PI period due to an earlier decline in insolation over the northern location of Intertropical Convergence Zone (ITCZ).

Description: In order to quantify the relative importance of individual boundary conditions and forcings, including greenhouse gases, ice sheets, and Earth's orbital parameters, on determining Last Glacial Maximum (LGM) climate, we have performed a series of LGM experiments using a state-of-the-art climate model AWI-ESM, in which different combinations of boundary conditions and forcings have been applied following the protocol of Paleoclimate Modelling Intercomparison Project phase 4 (PMIP4). In good agreement with observational proxy records, a general colder and drier climate is simulated in our full-forced LGM experiment as compared to the present-day simulation. Our simulated results from non-full-forced sensitivity simulations reveal that both the greenhouse gases and ice sheets play a major role in defining the anomalous LGM surface temperature compared to today. Decreased greenhouse gases in LGM as compared to present day leads to a non-uniform global cooling with polar amplification effect. The presence of LGM ice sheets favors a warming over the Arctic and northern Atlantic oceans in boreal winter, as well as a cooling over regions with the presence of ice sheets. The former is induced by a strengthening in the Atlantic meridional overturning circulation (AMOC) transporting more heat to high latitudes, whilst the latter is due to the increased surface albedo and elevation of ice sheets. We find that the Northern Hemisphere monsoon precipitation is influenced by the opposing effects of LGM greenhouse gases and ice sheets. Specifically, the presence of ice sheets leads to significant drying in the Northern Hemisphere monsoon regions, while a reduction in greenhouse gases results in increased monsoon rainfall. Based on our model results, continental ice sheets exert a major control on atmospheric dynamics and the variability of El Niño–Southern Oscillation (ENSO). Moreover, our analysis also implies a nonlinearity in climate response to LGM boundary conditions and forcings.

Description: Recent global warming is pronounced in high-latitude regions (e.g. northern Asia), and will cause the vegetation to change. Future vegetation trends (e.g. the “arctic greening”) will feed back into atmospheric circulation and the global climate system. Understanding the nature and causes of past vegetation changes is important for predicting the composition and distribution of future vegetation communities. Fossil pollen records from 468 sites in northern and eastern Asia were biomised at selected times between 40 cal ka bp and today. Biomes were also simulated using a climate-driven biome model and results from the two approaches compared in order to help understand the mechanisms behind the observed vegetation changes. The consistent biome results inferred by both approaches reveal that long-term and broad-scale vegetation patterns reflect global to hemispheric-scale climate changes. Forest biomes increase around the beginning of the late deglaciation, become more widespread during the early and middle Holocene, and decrease in the late Holocene in fringe areas of the Asian Summer Monsoon. At the southern and southwestern margins of the taiga, forest increases in the early Holocene and shows notable species succession, which may have been caused by winter warming at ca. 7 cal ka bp. At the northeastern taiga margin (central Yakutia and northeastern Siberia), shrub expansion during the last deglaciation appears to prevent the permafrost from thawing and hinders the northward expansion of evergreen needle-leaved species until ca. 7 cal ka bp. The vegetation climate disequilibrium during the early Holocene in the taiga-tundra transition zone suggests that projected climate warming will not cause a northward expansion of evergreen needle-leaved species.

Description: The equilibrium climate sensitivity (ECS) of climate models is calculated as the equilibrium global mean surface air warming resulting from a simulated doubling of the atmospheric CO2 concentration. In these simulations, long-term processes in the climate system, such as land ice changes, are not incorporated. Hence, climate sensitivity derived from paleodata has to be compensated for these processes, when comparing it to the ECS of climate models. Several recent studies found that the impact these long-term processes have on global temperature cannot be quantified directly through the global radiative forcing they induce. This renders the prevailing approach of deconvoluting paleotemperatures through a partitioning based on radiative forcings inaccurate. Here, we therefore implement an efficacy factor ε[LI] that relates the impact of land ice changes on global temperature to that of CO2 changes in our calculation of climate sensitivity from paleodata. We apply our refined approach to a proxy-inferred paleoclimate dataset, using ε[LI]=0.45+0.34−0.20 based on a multi-model assemblage of simulated relative influences of land ice changes on the Last Glacial Maximum temperature anomaly. The implemented ε[LI] is smaller than unity, meaning that per unit of radiative, forcing the impact on global temperature is less strong for land ice changes than for CO2 changes. Consequently, our obtained ECS estimate of 5.8±1.3 K, where the uncertainty reflects the implemented range in ε[LI], is ∼50 % higher than when differences in efficacy are not considered.

Description: Teleconnections refer to links between regions that are distant to each other, but nevertheless exhibit some relation. The study of such teleconnections is a well-known task in climate research. Climate simulation shall model known teleconnections. Detecting teleconnections in climate simulations is a crucial aspect in judging the quality of the simulation output. It is common practice to run scripts to execute a sequence of analysis steps on the climate simulations to search for teleconnections. Such a scripting approach is not flexible and targeted towards one specific goal. It is desirable to have one tool that allows for a flexible analysis of all teleconnection patterns with a dataset. We present such a tool, where the extracted information is provided in an intuitive visual form to users, who then can interactively explore the data. We developed an analysis workflow that is modeled around four views showing different facets of the data with coordinated interaction. We present a teleconnection study with simulation ensembles and reanalysis data obtained by data assimilation to observe how well the teleconnectivity patterns match and to demonstrate the effectiveness of our tool.

Description: CO2-induced global warming will affect seasonal to decadal temperature patterns. Expected changes will be particularly strong in extratropical regions where temperatures will increase at faster rates than at lower latitudes. Despite that, it is still poorly constrained how precisely short-term climate dynamics will change in a generally warmer world, particularly in nearshore surface waters in the extratropics, i.e., the ecologically most productive regions of the ocean on which many human societies depend. Specifically, a detailed knowledge of the relationship between pCO2 and seasonal SST is crucial to understand interactions between the ocean and the atmosphere. In the present investigation, we have studied for the first time how rising atmospheric pCO2 levels forced surface temperature changes in Central Europe (paleolatitude ~45 °N) during the mid-Oligocene (fromca. 31 to 25Ma), a time interval of Earth history during which global conditions were comparable to those predicted for the next few centuries. For this purpose, we computed numerical climate models for the Oligocene (winter, summer, annual average) assuming an atmospheric carbon dioxide rise from 400 to 560 ppm (current level to two times pre-industrial levels, PAL) and from 400 to 840 ppm (= three times PAL), respectively. These models were compared to seasonally resolved sea surface temperatures (SST) reconstructed from δ18O values of fossil bivalve shells (Glycymeris planicostalis, G. obovata, Palliolum pictum, Arctica islandica and Isognomon maxillata sandbergeri) and shark teeth (Carcharias cuspidata, C. acutissimaand Physogaleus latus) collected fromthe shallow water deposits of the Mainz and Kassel Basins (Germany). Multi-taxon oxygen isotope-based reconstructions suggest a gradual rise of temperatures in surface waters (upper 30 to 40m), on average, by asmuch as 4 °C during the Rupelian stage followed by a 4 °C cooling during the Chattian stage. Seasonal temperature amplitudes increased by ca. 2 °C during the warmest time interval of the Rupelian stage,withwarming beingmore pronounced during summer (5 °C) than during winter (3 °C). According to numerical climate simulations, the warming of surface waters during the early Oligocene required a CO2 increase by at least 160 ppm, i.e., 400 ppm to 560 ppm. Given that atmospheric carbon dioxide levels predicted for the near future will likely exceed this value significantly, the Early Oligocene warming gives a hint of the possible future climate in Central Europe under elevated CO2 levels.

Description: Tree-ring maximum latewood density (MXD) records from Fennoscandia have been widely used to infer regional- and hemispheric-scale mean temperature variability. Here, we explore if MXD records can also be used to infer past variability of summer temperature extremes across Fennoscandia. The first principal component (PC1) based on 34 MXD chronologies in Fennoscandia explains 50% of the total variance in the observed warm-day extremes over the period 1901–1978. Variations in both observed summer warm-day extremes and PC1 are influenced by the frequency of anomalous anticyclonic pattern over the region, summer sea surface temperatures over the Baltic, North and Norwegian Seas, and the strength of the westerly zonal wind at 200 hPa across Fennoscandia. Both time series are associated with nearly identical atmospheric circulation and SST patterns according to composite map analysis. In a longer context, the first PC based on 3 millennium-long MXD chronologies in central and northern Fennoscandia explains 83% of the total variance of PC1 from the 34 MXD chronologies over the period 1901–1978, 48% of the total variance of the summer warm-day extreme variability over the period 1901–2006, and 36% of the total variance in the frequency of a summer anticyclonic pattern centered over eastern-central Fennoscandia in the period 1948–2006. The frequency of summer warm-day extremes in Fennoscandia is likely linked to a meridional shift of the northern mid-latitude jet stream. This study shows that the MXD network can be used to infer the variability of past summer warm-day extremes and the frequency of the associated summer anticyclonic circulation pattern over Fennoscandia.