Description: A new global climate model setup using FESOM2.0 for the sea ice‐ocean component and ECHAM6.3 for the atmosphere and land surface has been developed. Replacing FESOM1.4 by FESOM2.0 promises a higher efficiency of the new climate setup compared to its predecessor. The new setup allows for long‐term climate integrations using a locally eddy‐resolving ocean. Here it is evaluated in terms of (1) the mean state and long‐term drift under preindustrial climate conditions, (2) the fidelity in simulating the historical warming, and (3) differences between coarse and eddy‐resolving ocean configurations. The results show that the realism of the new climate setup is overall within the range of existing models. In terms of oceanic temperatures, the historical warming signal is of smaller amplitude than the model drift in case of a relatively short spin‐up. However, it is argued that the strategy of “de‐drifting” climate runs after the short spin‐up, proposed by the HighResMIP protocol, allows one to isolate the warming signal. Moreover, the eddy‐permitting/resolving ocean setup shows notable improvements regarding the simulation of oceanic surface temperatures, in particular in the Southern Ocean.

Description: Recent evidence shows that wind‐driven ocean currents, like the western boundary currents, are strongly affected by global warming. However, due to insufficient observations both on temporal and spatial scales, the impact of climate change on large‐scale ocean gyres is still not clear. Here, based on satellite observations of sea surface height and sea surface temperature, we find a consistent poleward shift of the major ocean gyres. Due to strong natural variability, most of the observed ocean gyre shifts are not statistically significant, implying that natural variations may contribute to the observed trends. However, climate model simulations forced with increasing greenhouse gases suggest that the observed shift is most likely to be a response of global warming. The displacement of ocean gyres, which is coupled with the poleward shift of extratropical atmospheric circulation, has broad impacts on ocean heat transport, regional sea level rise, and coastal ocean circulation.

Description: Lack of constraint on spatial and long-term temporal variability of the El Niño southern Oscillation (ENSO) and its sensitivity to external forcing limit our ability to evaluate climate models and ENSO future projections. Current knowledge of Holocene ENSO variability derived from paleoclimate reconstructions does not separate the role of insolation forcing from internal climate variability. Using an updated synthesis of coral and bivalve monthly resolved records, we build composite records of seasonality and interannual variability in four regions of the tropical Pacific: Eastern Pacific (EP), Central Pacific (CP), Western Pacific (WP) and South West Pacific (SWP). An analysis of the uncertainties due to the sampling of chaotic multidecadal to centennial variability by short records allows for an objective comparison with transient simulations (mid-Holocene to present) performed using four different Earth System models. Sea surface temperature and pseudo-δ18O are used in model-data comparisons to assess the potential influence of hydroclimate change on records. We confirm the significance of the Holocene ENSO minimum (HEM) 3-6ka compared to low frequency unforced modulation of ENSO, with a reduction of ENSO variance of ∼50 % in EP and ∼80 % in CP. The approach suggests that the increasing trend of ENSO since 6ka can be attributed to insolation, while models underestimate ENSO sensitivity to orbital forcing by a factor of 4.7 compared to data, even when accounting for the large multidecadal variability. Precession-induced change in seasonal temperature range is positively linked to ENSO variance in EP and to a lesser extent in other regions, in both models and observations. Our regional approach yields insights into the past spatial expression of ENSO across the tropical Pacific. In the SWP, today under the influence of the South Pacific Convergence Zone (SPCZ), interannual variability was increased by ∼200 % during the HEM, indicating that SPCZ variability is independent from ENSO on millennial time scales.

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: Numerical simulations employing prognostic sta- ble water isotopes can not only facilitate our understanding of hydrological processes and climate change but also al- low 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 We- gener 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) simula- tions. The model reproduces the observed PI isotope compo- sitions in both precipitation and seawater well and captures their major differences from the MH conditions. The sim- ulated relationship between the isotope composition in precipitation (d18Op) and surface air temperature is very similar between the PI and MH conditions, and it is largely consis- tent 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 pre- cipitation 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 dis- tinction between pre-monsoon and monsoon seasons. Com- bining 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 with- drawal 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: Improving our understanding of the controls on Antarctic precipitation is critical for gaining insights into past and future polar and global environmental changes. Here we develop innovative water tracing diagnostics in the atmospheric general circulation model ECHAM6. These tracers provide new detailed information on moisture source locations and properties of Antarctic precipitation. In the preindustrial simulation, annual mean Antarctic precipitation originating from the open ocean has a source latitude range of 49–35◦ S, a source sea surface temperature range of 9.8– 16.3 ◦ C, a source 2 m relative humidity range of 75.6 %– 83.3%, and a source 10m wind velocity (vel10) range of 10.1 to 11.3ms−1. These results are consistent with estimates from existing literature. Central Antarctic precipitation is sourced from more equator-ward (distant) sources via elevated transport pathways compared to coastal Antarctic precipitation. This has been attributed to a moist isentropic framework; i.e. poleward vapour transport tends to follow constant equivalent potential temperature. However, we find notable deviations from this tendency especially in the lower troposphere, likely due to radiative cooling. Heavy precipitation is sourced by longer-range moisture transport: it comes from 2.9◦ (300 km, averaged over Antarctica) more equator-ward (distant) sources compared to the rest of precipitation. Precipitation during negative phases of the Southern Annular Mode (SAM) also comes from more equator-ward moisture sources (by 2.4◦, averaged over Antarctica) compared to precipitation during positive SAM phases, likely due to amplified planetary waves during negative SAM phases. Moreover, source vel10 of annual mean precipitation is on average 2.1 m s−1 higher than annual mean vel10 at moisture source locations from which the precipitation originates. This shows that the evaporation of moisture driving Antarctic precipitation occurs under windier conditions than average. We quantified this dynamic control of Southern Ocean surface wind on moisture availability for Antarctic precipitation. Overall, the innovative water tracing diagnostics enhance our under- standing of the controlling factors of Antarctic precipitation.

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: The Atlantic Meridional Overturning Circulation (AMOC) is a key feature of the North Atlantic with global ocean impacts. The AMOC's response to past changes in forcings during the Holocene provides important context for the coming centuries. Here, we investigate AMOC trends using an emerging set of transient simulations using multiple global climate models for the past 6,000 years. Although some models show changes, no consistent trend in overall AMOC strength during the mid-to-late Holocene emerges from the ensemble. We interpret this result to suggest no overall change in AMOC, which fits with our assessment of available proxy reconstructions. The decadal variability of the AMOC does not change in ensemble during the mid- and late-Holocene. There are interesting AMOC changes seen in the early Holocene, but their nature depends a lot on which inputs are used to drive the experiment.

Description: Changes of the Atlantic meridional overturning circulation (AMOC) in the mid‐Holocene compared to the preindustrial state are explored in different coupled climate models. Using time‐slice integrations by a newly developed global finite‐element model ECHAM6‐FESOM with unstructured mesh and high resolution, our simulations show an enhanced mid‐Holocene AMOC, accompanied by an increase in the ocean salinity over regions of deep water formation. We identify two different processes affecting the AMOC: (1) a more positive phase of North Atlantic Oscillation (NAO) increased water density over the Labrador Sea through anomalous net evaporation and surface heat loss; (2) a decreased import of sea ice from the Arctic causes a freshwater reduction in the northern North Atlantic Ocean. Using the coupled model ECHAM6‐MPIOM in T63GR15 and T31GR30 grids, we find that the simulated AMOC has significant discrepancy with different model resolutions. In detail, stronger‐than‐present mid‐Holocene AMOC is revealed by simulations with the T63GR15 grid, which resembles the result of ECHAM6‐FESOM, while a decline of the mid‐Holocene AMOC is simulated by the low resolution model with the T31GR30 grid. Such discrepancy can be attributed to different changes in Labrador Sea density which is mainly affected by (1) NAO‐induced net precipitation and deep water convection, (2) freshwater transport from the Arctic Ocean, and (3) the strength of AMOC itself. Finally, we analyzed available coupled climate models showing a diversity of responses of AMOC to mid‐Holocene forcings, most of which reveal positive AMOC changes related to northern high latitudes salinification.