Beschreibung: During the Late Pleistocene and Holocene retreat of paleo-ice sheets in North America and Europe, vast proglacial lakes existed along the land terminating margins. These proglacial lakes impacted ice sheet dynamics by imposing boundary conditions analogous to a marine terminating margin. Such lacustrine boundary conditions cause changes in the ice sheet geometry, stress balance and frontal ablation and therefore affect the mass balance of the entire ice sheet. Despite this, dynamically evolving proglacial lakes have rarely been considered in detail in ice sheet modeling endeavors. In this study, we describe the implementation of an adaptive lake boundary in the Parallel Ice Sheet Model (PISM), which we call PISM-LakeCC. We test our model with a simplified glacial retreat setup of the Laurentide Ice Sheet (LIS). By comparing the experiments with lakes to control runs with no lakes, we show that the presence of proglacial lakes locally enhances the ice flow, which leads to a lowering of the ice sheet surface. In some cases, this also results in an advance of the ice margin and the emergence of ice lobes. In the warming climate, increased melting on the lowered ice surface drives the glacial retreat. For the LIS, the presence of lakes triggers a process similar to marine ice sheet instability, which caused the collapse of the ice saddle over Hudson Bay. In the control experiments without lakes, Hudson Bay is still glaciated when the climate reaches present-day (PD) conditions. The results of our study demonstrate that glacio-lacustrine interactions play a significant role in the retreat of land terminating ice sheet margins.

Beschreibung: Climate and environments of the mid-Pliocene warm period (3.264 to 3.025 Ma) have been extensively studied. Whilst numerical models have shed light on the nature of climate at the time, uncertainties in their predictions have not been systematically examined. The Pliocene Model Intercomparison Project quantifies uncertainties in model outputs through a coordinated multi-model and multi-model/data intercomparison. Whilst commonalities in model outputs for the Pliocene are clearly evident, we show substantial variation in the sensitivity of models to the implementation of Pliocene boundary conditions. Models appear able to reproduce many regional changes in temperature reconstructed from geological proxies. However, data/model comparison highlights that models potentially underestimate polar amplification. To assert this conclusion with greater confidence, limitations in the time-averaged proxy data currently available must be addressed. Furthermore, sensitivity tests exploring the known unknowns in modelling Pliocene climate specifically relevant to the high latitudes are essential (e.g. palaeogeography, gateways, orbital forcing and trace gasses). Estimates of longer-term sensitivity to CO2 (also known as Earth System Sensitivity; ESS), support previous work suggesting that ESS is greater than Climate Sensitivity (CS), and suggest that the ratio of ESS to CS is between 1 and 2, with a "best" estimate of 1.5.

Beschreibung: 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.

Beschreibung: 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).

Beschreibung: The response of evapotranspiration to anthropogenic warming is of critical importance for the water and carbon cycle. Contradictory conclusions about evapotranspiration changes are caused primarily by their brevity in time and sparsity in space, as well as the strong influence of internal variability. Here, we present the first gridded reconstruction of the summer (June, July, and August) vapor pressure deficit (VPD) for the past 4 centuries at the European level. This gridded reconstruction is based on 26 European tree ring oxygen isotope records and is obtained using a random forest approach. According to validation scores obtained with the Nash-Sutcliffe model efficiency, our reconstruction is robust over large parts of Europe since 1600, in particular for the westernmost and northernmost regions, where most tree ring records are located. Based on our reconstruction, we show that from the mid-1700s a trend towards higher summer VPD occurred in central Europe and the Mediterranean region that is related to a simultaneous increase in temperature and decrease in precipitation. This increasing summer VPD trend continues throughout the observational period and in recent times. Moreover, our summer VPD reconstruction helps to visualize the local and regional impacts of the current climate change, as well as to minimize statistical uncertainties of historical VPD variability. This paper provides also new insights into the relationship between summer VPD and large-scale atmospheric circulation, and we show that summer VPD has two preferred modes of variability, namely a NW-SE dipole-like mode and a N-S dipole-like mode. Furthermore, the interdisciplinary use of the data should be emphasized, as summer VPD is a crucial parameter for many climatological feedback processes in the Earth's surface system. The reconstructed summer VPD gridded data over the last 400 years are available at the following link: 10.5281/zenodo.5958836 (Balting et al., 2022).

Beschreibung: We compare the main atmospheric drivers of the melt season over the Greenland Ice Sheet (GrIS) in ERA5 and ERA-Interim (ERAI) in their overlapping period 1979–2018. In summer, ERA5 differs significantly from ERAI, especially in the melt regions. Small-scale ERA5 − ERAI differences near the ice sheet’s margins and over steep slopes can be explained by the different resolution, while the large-scale differences indicate a different representation of physical processes in the two reanalyses: averaged over the lower parts of the GrIS, the mean near-surface air temperature is 1 K lower, while the mean downward shortwave radiation at the surface is on average 15 W m−2 higher than in ERAI. Comparison with observational weather station data shows a significant warm bias in ERAI and, for ERA5, a significant positive bias in downward shortwave radiation. Consequently, methods that previously estimated the GrIS surface mass balance from the ERAI surface energy balance need to be carefully recalibrated before converting to ERA5 forcing.