The past two thousand years is a key interval for climate science. This period encompasses both the era of human-induced global warming and a much longer interval when changes in Earth’s climate were governed principally by natural drivers and unforced variability. Since 2009, the Past Global Changes (PAGES) 2k Network has brought together hundreds of scientists from around the world to reconstruct and understand the climate of the Common Era using open and collaborative approaches to palaeoclimate science, including virtual meetings. The third phase of the network will end in December 2021. Here we highlight some key outputs of PAGES 2k and present the major themes and scientific questions emerging from recent surveys of the community. We explore how these might boost a new phase of PAGES 2k or a successor project(s). This year we will further reach out to the community through Town Hall consultations, vEGU and other meetings, and a PAGES 2k global webinar series. The aim of these activities is to foster development of post-2021 community-led PAGES initiatives that connect past and present climate.

How to cite: Eggleston, S. S., Phipps, S., Bothe, O., McGregor, H. V., Martrat, B., Linderholm, H., Konecky, B., Abram, N., and St. George, S.: Advancing community-led research into the climate of the Common Era, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-10979, https://doi.org/10.5194/egusphere-egu21-10979, 2021.

Better understanding of current climate changes needs a full knowledge about regional specific of thermal conditions at the end of Little Ice Age. So, the earliest available meteorological data are important. First regular daily qualitative meteorological observations were taken in Moscow city from 1657 to 1675. Episodic short series of instrumental measurements were made there for the first time in 1731; regular daily measurements started in 1779 when one of Mannheim network stations was founded in Moscow.

         All known old data series of the air temperature T measurements in Moscow since 1779 were collected and analyzed. Mannheim station existed there from 1779 to 1797 but average values of T are available from issues of Ephemerides Societatis Meteorologicae Palatinae only for the period 1779–1792. High accuracy of measurements at Mannheim network is confirmed by high correlation co-efficient between monthly-averaged T values in Moscow and at closest stations (Warsaw and St. Petersburg): up to 0.82-0.84 on separate months.

         Different methodical questions (unknown location of the station, unknown conditions of thermometer installation, its height and shading, an accuracy of its calibration, etc.) were studied. As a result it was found that the most probable error due to thermometer installation close to the northern building wall is ±0.1÷0.2 ºС; the error of daily-averaged T due to unknown height of measurements is ±0.1 ºС; the calibration accuracy in Mannheim was about ±0.1 ºС. Thus, a total error of T on average of a day in the 18^th century was not higher than ±0.3÷0.4 ºС. Probably it was even less because separate components of the error may be multidirectional. For the first time mean-annual T in Moscow was received for 1783, and the most probable values were estimated for 1784 and 1785 using the data of the closest Mannheim station (Saint-Petersburg) for separate months with data gaps. The end of Little Ice Age manifeted at extremely low minimal values of T: up to –31 ˚R (–38.8 ˚С) in December 17^th, 1788. However, thermal conditions from June to September changed only a bit since the 18^th century till nowadays (differences are not statistically significant with the 0.95 confidence probability).  

         Later measurements in Moscow were renewed since 1808 and broken again in August of 1812 due to Napoleon’s invasion and terrible Moscow fire. For the first time unknown data series of everyday measurements which were made by Ivan Lange in 1816–1817 were found and studied. As is known the famous 1816 ‘Year Without a Summer’ was noted almost all over the World by extremely cold summer as probable result of Mount Tambora eruption in 1815. Nevertheless, it was found that summer of 1816 in Moscow was comparatively cool but not extremely cold: monthly-averaged T there was 15.7, 17.3 and 14.5 ˚С in June, July and August, respectively, and 15.8 ˚С on average of the summer. Thus, 1816 occupies only 27^th place among the coldest summers in the city during 216 years.

         Author is thankful to the memory of his late PhD student Ekaterina L. Vasilenko.

How to cite: Lokoshchenko, M. A.: Climate of Moscow at the end of Little Ice Age, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-6663, https://doi.org/10.5194/egusphere-egu21-6663, 2021.

Based on copies of the original data (source: Oeschger Center for Climate Change Research) we perform climate reconstructions for Paris between 1665 - 1709. The focus lies on the following meteorological variables: temperature, cloudiness, direction of movement of the clouds, precipitation and humidity. Apart from humidity, these meteorological variables were measured three times a day over the entire period from Louis Morin. Temperature and humidity were measured with instruments, whereas cloud cover, direction of movement of the clouds and precipitation were measured in a descriptive manner. In addition to the quantitative temperature measurements, conclusions about synoptic air movements over Europe are possible due to the additional meteorological variables. The Late Maunder Minimum is characterised by cold winters and moderate summers. Winter is characterised by a lower frequency of westerly direction of movement of the clouds. This reduction of advection from the ocean leads to cooling in Paris and also to less precipitation in winter. This can be seen very strongly between the last decade of the 17^th century (cold) and the first decade of the 18^th century (warm). A lower frequency of westerly direction of movement of the clouds can also be seen in summer, but the influence is stronger in winter than in summer. However, this reduction leads to moderate/warm temperatures in summer. So unusually cold winters in the Late Maunder Minimum can be attributed to a lower frequency of westerly direction of movement of the clouds.

How to cite: Pliemon, T., Foelsche, U., Rohr, C., and Pfister, C.: Analysis of Subdaily Meteorological Measurements by Louis Morin in the Late Maunder Minimum 1665 – 1709 in Paris, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-11953, https://doi.org/10.5194/egusphere-egu21-11953, 2021.

This work describes the compilation of global instrumental climate data with a focus on the 18th and early 19th centuries. This database provides early instrumental data recovered for thousands of locations around the world. Instrumental meteorological measurements from periods prior to the start of national weather services are designated “early instrumental data”. Much of the data is taken from repositories we know (GHCN, ISTI, CRUTEM, Berkeley Earth, HISTALP). In addition, many of these stations have not been digitized before. Therefore,  we provide a new global collection of monthly averages of multivariable meteorological parameters before 1890 based on land-based meteorological station data. The product will be form as the most comprehensive global monthly climate data set, encompassing temperature, pressure, and precipitation as ever done. These data will be quality controlled and analyzed with respect to climate variability and they be assimilated into global climate model simulations to provide monthly global reconstructions. The collection has resulted in a completely new database that is uniform, where no interpolations are included. Therefore, we are left with climate reconstruction that becomes very authentic. This compilation will describe the procedure and various challenges we have encountered by creating a unified database that can later be used for e.g. models. It will also describe the strategy for quality control that has been adopted is a sequence of tests.

How to cite: Lundstad, E., Brugnera, Y., and Brönnimann, S.: Global Instrumental Meteorological Database Before 1890 – a useful overview, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-15903, https://doi.org/10.5194/egusphere-egu21-15903, 2021.

Energy exchanges among climate subsystems are of critical importance to determine the climate sensitivity of the Earth's system to changes in external forcing, to quantify the magnitude and evolution of the Earth's energy imbalance, and to make projections of future climate. Additionally, climate phenomena sensitive to land heat storage, such as permafrost stability and sea level rise, are important due to their impacts on society and ecosystems. Thus, ascertaining the magnitude and change of the Earth's energy partition within climate subsystems has become urgent in recent years. 

Here, we provide new global estimates of changes in ground surface temperature, ground surface heat flux and continental heat storage derived from geothermal data using an expanded database and new techniques developed in the last two decades. This new dataset contains 253 recent borehole profiles that were not included in previous estimates of global continental heat storage. In addition, our analysis considers additional sources of uncertainty that were not included in previous borehole studies. Results reveal markedly higher changes in ground heat flux and heat storage within the continental subsurface during the second half of the 20th century than previously reported, with a land mean temperature increase of 1 K and continental heat gains of around 12 ZJ relative to preindustrial times. Half of the heat gained by the continental subsurface since 1960 have occurred in the last twenty years. These results may be important for estimates of climate sensitivity based on energy budget constrains, as well as for the evaluation of global transient climate simulations in terms of the Earth’s heat inventory and energy-dependent subsurface processes. Our estimate of land heat storage is included in the new assessment of the components of the Earth’s heat inventory recently released (von Schuckmann et al. 2020), together with the oceans, the atmosphere and the cryosphere.

How to cite: Cuesta-Valero, F. J., García-García, A., Beltrami, H., González-Rouco, J. F., and García-Bustamante, E.: Long-Term Global Ground Heat Flux and Continental Heat Storage from Geothermal Data, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-589, https://doi.org/10.5194/egusphere-egu21-589, 2021.

The surface temperature response to changes in our planet’s external forcing is larger at higher latitudes, a phenomenon known as polar amplification. The Arctic amplification has been particularly intense during the last century, with arctic-wide paleoclimatic reconstructions and state-of-the-art model simulations revealing a twofold arctic warming in comparison with the average global temperature increase. As a consequence, Arctic ground temperatures respond with rapid warming, but this response varies with snow cover and permafrost processes. Thus, changes in arctic ground temperatures are difficult to reconstruct from data, and to simulate in climate models.

Here, we reconstruct the ground surface temperature histories of 120 borehole temperature profiles above 60ºN for the last 400 years. Past surface temperature evolution from each profile was estimated using a Perturbed Parameter Inversion approach based on a singular value decomposition method. Long-term surface temperature climatologies (circa 1300 and 1700 CE) and quasi-steady state heat flow are also estimated from linear regression through the depth range 200 to 300 m of each borehole temperature profile. The retrieved temperatures are assessed against simulated ground surface temperatures from five Past Millennium and five Historical experiments from the Paleoclimate Modelling Intercomparison Project Phase III (PMIP3), and the fifth phase of the Coupled Model Intercomparison Project (CMIP5) archives, respectively.

Preliminary results from borehole estimates and PMIP3/CMIP5 simulations reveal that changes in recent Arctic ground temperatures vary spatially and are related to each site’s earlier thermal state of the surface. The magnitudes of ground warming from data and simulations differ with large discrepancies among models. As a consequence, a better understanding of freezing processes at and below the air-ground interface is necessary to interpret subsurface temperature records and global climate model simulations in the Arctic.

How to cite: Beltrami, H., Cuesta-Valero, F. J., García-García, A., Gruber, S., and Jaume-Santero, F.: Assessing Arctic Ground Surface Temperatures from Borehole Temperatures and Paleoclimatic Model Simulations, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-9222, https://doi.org/10.5194/egusphere-egu21-9222, 2021.

Historical solar irradiance is a critical input to climate models. As no direct measurements are available before 1978, reconstructions of past irradiance changes are employed instead. Such reconstructions are based on the knowledge that solar irradiance on time scales of interest to climate studies is modulated by the evolution of the solar surface magnetic structures, such as sunspots and faculae. This calls for historical records or proxies of such features. The longest direct, and thus mostly used, record is the sunspot number. It allows a reasonable description of the emergence and evolution of active regions, which are larger magnetic regions containing sunspots. At the same time, a significant amount of the magnetic flux on the Sun emerges in the form of the so-called ephemeral magnetic regions, which are weaker short-lived bipolar regions that do not contain sunspots. Due to their high frequency, ephemeral regions are an important source of the irradiance variability, especially on time scales longer than the solar cycle. Difficulties in their proper accounting are a main reason for the high uncertainty in the secular irradiance variability. Existing models either do not account for their evolution at all or link them linearly to active regions. We use a new, more realistic model of the ephemeral region emergence, relying on recent independent solar observations, as input to a surface flux transport model (SFTM) to simulate the evolution of the magnetic field in such regions. The latter can then be used to reconstruct the solar irradiance since the Maunder minimum.

How to cite: Hofer, B., Krivova, N. A., Solanki, S. K., Cameron, R., Wu, C.-J., and Usoskin, I. G.: Accounting for small bipolar magnetic regions in solar irradiance reconstructions, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-8323, https://doi.org/10.5194/egusphere-egu21-8323, 2021.

This study presents a reconstruction of climate change in central Korea during the last 3,000 years, using a core from a montane peatland of Yongneup. Multiple proxies of pollen, macrocharcoal, and geochemistry were analysed to provide three findings as follows: First, abrupt climate events at ca. 2.8 and 2.3 ka BP possibly accompanied dry summer as well as cold and arid winter seasons on the Korean peninsula. The first macrocharcoal analysis on the peninsula indicates increased wildfire activities during these dry periods. Next, a weakening of summer monsoon during El Niño-like phases was clearly found during the late Holocene. This confirms previous findings of a dominant oceanic influence on hydroclimate variability on the Korean peninsula. Finally, changes in temperature were likely synchronous with a global trend, indicated by the total organic content (TOC) and arboreal pollen percentages. Due to its location at a high-altitude, the environment of Yongneup has possibly sensitively responded to fluctuations in temperature. Altogether, these findings suggest that temperature and precipitation changes on the Korean peninsula have been separately influenced by insolation and oceanic circulations, respectively.

How to cite: Park, J., Jin, Q., Choi, J., and Park, J.: Late-Holocene climate variability of coastal East Asia reconstructed from the Yongneup fen in central Korea, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-1893, https://doi.org/10.5194/egusphere-egu21-1893, 2021.

The increase in anthropogenic induced warming over the last two centuries is impacting marine environments. Marine planktic calcifying organisms interact sensitively to changes in sea surface temperatures (SST), and the food web structure. Here, we study two high resolution multicore records from two western Mediterranean Sea regions (Alboran and Balearic basins), areas highly affected by both natural climate change and anthropogenic warming. Cores cover the time interval from the Medieval Climate Anomaly (MCA) to present. Reconstructed SSTs are in good agreement with other results, tracing temperature changes through the Common Era, and show a clear 20^th century warming signal. Both cores show opposite abundance fluctuations of planktic foraminiferal species (Globigerina bulloides, Globorotalia inflata and Globorotalia truncatulinoides) a common group of marine calcifying zooplankton. The abundance ratios between these species show the switch between winter / spring surface productivity and deep winter mixing in the Balearic basin. In the Alboran Sea, Globigerina bulloides and Globorotalia inflata instead respond to local upwelling dynamics. In the pre-industrial era, changes in planktic foraminiferal productivity and species composition can be explained mainly by the natural variability of North Atlantic Oscillation (NAO), and, to lesser extent, by the Atlantic Multidecadal Oscillation (AMO). In the industrial era, starting from about 1800 Common Era (CE), this variability is affected by anthropogenic surface warming, leading to enhanced vertical stratification of the upper water column, and resulting in a decrease of surface productivity at both sites. We found that natural planktic foraminiferal population dynamics in the western Mediterranean is already altered by enhanced anthropogenic impact in the industrial era, suggesting that in this region natural cycles and influences are being overprinted by human influences.

How to cite: Pallacks, S., Ziveri, P., Martrat, B., Mortyn, G. P., Grelaud, M., Schiebel, R., Incarbona, A., Garcia-Orellana, J., and Anglada-Ortiz, G.: Planktic Foraminifera changes in the western Mediterranean Anthropocene, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-13653, https://doi.org/10.5194/egusphere-egu21-13653, 2021.

There is a major knowledge gap in the past climate oscillation of the Arabian desert, especially during the past two millennium. Reliable continuous continental records that archives at high resolution past environmental variability are useful sentinels of paleoclimate changes. Reliable interpretation from climatic proxies retrieved from lake records are crucial for identifying periodicities and the onset of climatic events and evaluating inter-annual and decadal trends driven by shifting of the Intertropical Convergence Zone (ITCZ). A multiproxy approach is presented for a ~3.3 m composite core from a karst lake located in Gayal el Bazal, southern Yemen. Sedimentary proxies, including grain size distribution and magnetic susceptibility (MS) coupled with geochemistry (XRF), provide an initial picture of centennial-scale environmental changes over the southern Arabian desert. The chronology of the core was anchored by five radiocarbon (¹⁴C) dates of terrestrial plants (wood) extracted from sediment samples and indicates the core extends to ~800 AD. Our data provides a snapshot for better understanding the impact of Indian Ocean monsoon variability at an exceptional resolution for a region that lacks sufficient information. Our data indicates that during the ‘Little Ice Age’ (~1500-1800 AD) was arid relative to the warm conditions that prevailed during the Medieval Warming Period (~800 to 1200 AD). The arid phase was marked by high Ca/(Al, Fe, Ti) values, increased inorganic carbon content, decreased MS values, and gypsum precipitation. Furthermore, end-member mixing analyses (EMMA) derived from the grain-size distribution corroborates the production of carbonate sand probably due to an increase in flash floods occurring concurrently with low lake levels under generally dry conditions. Aridity during the Little Ice Age is consistent with evidence and theory for weakened boreal summer monsoons during intervals of northern hemisphere cooling. Overall, this study will provide insight into the monsoon variability and a record for understanding the interactions between northward migrations of the ITCZ and tropical monsoonal dynamics during the late Holocene. In the context of current climate change and increasing population pressure, a deeper understanding of their long-term hydrological variability, this study is highly essential to satisfactorily forecast the sustainability of lakes as a resource in a warming world.

How to cite: Parth, S., Russell, J. M., and Waldmann, N.: Reconstructing 1200 years of hydroclimate variability in the southern margins of the Arabian Desert, inferred from an ancient lake in southern Yemen, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-14756, https://doi.org/10.5194/egusphere-egu21-14756, 2021.

Recent Antarctic surface climate change has been characterized by greater warming trends in West Antarctica than in East Antarctica. Although the changes over recent decades are well studied, the short instrumental record limits our ability to determine if such asymmetric patterns are common for Antarctica and the processes at their origin. Here, we will focus on the years 0-1000 CE as some ice core records display very contrasted trends during this period. Furthermore, the climate models are unable to reproduce the warming displayed in some reconstructions from 1 to 500 CE over East Antarctica. In order to understand the origin of these apparent incompatibilities and investigate the effect of proxy selection on regional reconstructions over 0-1000 CE, we performed several offline data assimilation experiments based on different groups of d¹⁸O records and the isotope-enabled general circulation models (iCESM). When assimilating different d18O data sets, large differences appear in the pattern of temperature trend over 0-500 CE, but the patterns over 500-1000 CE are more consistent among the various experiments. This implies that the spatial pattern of temperature trend over 0-500 CE is still uncertain because of this high sensitivity on the choice of the proxies to constrain the model results, while the pattern over 500-1000 is more robust, with the greater cooling over West Antarctica than East Antarctica. This pattern over 500-1000 CE relates to the intensifying of the low pressure centered in the Amundsen Sea, which induces enhanced southerly flow through most of WAIS.

How to cite: Lyu, Z., Goosse, H., and Dalaiden, Q.: Spatial patterns of multi-centennial temperature trend in Antarctica over 0-1000 CE: insights from ice core records and modeling, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-571, https://doi.org/10.5194/egusphere-egu21-571, 2021.

Weather systems in the southern Indian Ocean drive synoptic-scale precipitation, temperature and wind variability in East Antarctica, sub-Antarctic islands and southern Australia.  Over seasonal to decadal timescales, the mean condition associated with combinations of these synoptic weather patterns (e.g., extratropical cyclones, fronts and regions of high pressure) is often referred to as variability in the westerly wind belt or the Southern Annular Mode (SAM). The westerly wind belt is generally considered to be zonally symmetric around Antarctica however, on a daily timescale this is not the case. To capture the daily variability of regional weather systems, we used synoptic typing (Self-Organising Maps) to group weather patterns based on similar features, which are often lost when using monthly or seasonal mean fields. We identified nine key regional weather types based on anomaly pattern and strength. These include four meridional nodes, three mixed nodes, one zonal node and one transitional node. The meridional nodes are favourable for transporting warm, moist air masses to the subantarctic and Antarctic region, and are associated with increased precipitation and temperature where the systems interact with the Antarctic coastline.  These nodes have limited association with the SAM, especially during austral spring.  In contrast, the zonal and mixed nodes were strongly correlated with the SAM however, the regional synoptic representation of SAM positive conditions is not zonally symmetric and is represented by three separate nodes.  These different types of SAM positive conditions mean that the commonly used hemispheric Marshall index often fails to capture the regional variability in surface weather conditions in the southern Indian Ocean. Our results show the importance of considering different synoptic set ups of SAM conditions, particularly SAM positive, and identify conditions that are potentially missed by SAM variability (e.g., extreme precipitation events). Our results are particularly important to consider when interpreting SAM or westerly wind belt reconstructions in the study region (from ice cores, tree rings, or lake sediments).  Here we present a case study using the synoptic typing results to enhance our understanding of the Law Dome (East Antarctica) ice core record, focussing on links to large scale modes of climate variability and Australian hydroclimate.  These results enhance the usefulness of ice core proxies in coastal East Antarctica and assist with determining where and how it is appropriate to use coastal East Antarctic ice core records for reconstructions of large scale modes of climate variability (e.g. SAM and ENSO) and remote hydroclimate conditions.

How to cite: Udy, D., Vance, T., Kiem, A., Holbrook, N., and Curran, M.: Synoptic climatology of southern Indian Ocean and paleoclimate proxy interpretation, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-3756, https://doi.org/10.5194/egusphere-egu21-3756, 2021.

Great Basin Bristlecone pine (Pinus longaeva) is known for its longevity. The longest continuous tree-ring width chronology covers more than 9000 years. Tree-ring width of upper treeline bristlecone pine trees is influenced by summer temperature variability at decadal to centennial scales, but to infer a temperature signal on interannual scales, Maximum Latewood Density (MXD) is a better proxy. Here, we present a preliminary MXD chronology to investigate the temperature signal in upper treeline and lower elevation bristlecone pines. MXD was measured with an X-ray Computed Tomography toolchain in 24 dated cores,  with the oldest sample dating back to 776 CE. Ring and fibre angles were corrected and two MXD chronologies for different elevations were developed, which will be used to study climate-growth relationships and the effect of elevation on them. Future scanning will allow constructing a 5000+ year-long MXD chronology from upper treeline sites, which will provide an annual-resolution North American temperature record covering the mid-to-late Holocene.

How to cite: De Mil, T., Salzer, M., Pearson, C., Trouet, V., and Van den Bulcke, J.: Maximum latewood density records of Great Basin Bristlecone pine (Pinus Longaeva) from the White Mountains, California, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-7550, https://doi.org/10.5194/egusphere-egu21-7550, 2021.

We use an interdisciplinary approach combining stable isotopes in tree rings, pollen data, ice cores from temperature-limited environment in the Siberian north and developed a comprehensive description of the climatic changes over the past 1500 years. We found that the Climatic Optimum Period was warmer and drier compared to the Medieval one, but rather similar to the recent period. Our results indicate that the Medieval Warm period in the Taimyr Peninsula started earlier and was wetter compared to the northeastern part of Siberia (northeastern Yakutia). Summer precipitation reconstruction obtained from carbon isotopes in tree-ring cellulose from Taimyr Peninsula significantly correlated with the pollen data of the Lama Lake (Andreev et al. 2004) and oxygen isotopes of the ice core from Severnaya Zemlya (Opel et al. 2013) recording wetter climate conditions during the Medieval Warm period compared to the northeastern part of Siberia. Common large-scale climate variability was confirmed by significant relationship between oxygen isotope data in tree-ring cellulose from the Taimyr Peninsula and northeastern Yakutia, and oxygen isotope ice core data from Severnaya Zemlja during the Medieval Warm period and the recent one. Finally, we showed that the recent warming on the Taimyr Peninsula is not unprecedented in the Siberian north. Similar climate conditions were recorded by stable isotopes in tree rings, pollen, and ice core data 6000 years ago. On the northeastern part of Siberia newly developed a 1500-year summer vapor pressure deficit (VPD) reconstruction showed, that VPD increased recently, but does not yet exceed the maximum values reconstructed during the Medieval Warm period. The most humid conditions in the northeastern part of Siberia were recorded in the Early Medieval period and during the Little Ice Age. However, the increasing VPD under elevated air temperature in the last decades affects the hydrological regime of these sensitive ecosystems by greater evapotranspiration rates. Further VPD increase will significantly affect Siberian forests most likely leading to drought even under additional access of thawed permafrost water.

This work was supported by the FP7-PEOPLE-IIF-2008 - Marie Curie Action: "International Incoming Fellowships" 235122 and "Reintegration Fellowships" 909122 “Climatic and environmental changes in the Eurasian Subarctic inferred from tree-ring and stable isotope chronologies for the past and recent periods” and the Government of Krasnoyarsk Kray and Russian Foundation for Basic Research and Krasnoyarsk Foundation 20-44-240001 “Adaptation of conifer forests on the north of the Krasnoyarsk region (Taimyr Peninsula) to climatic changes after extreme events over the past 1500 years“ awarded to Olga V. Churakova (Sidorova).

How to cite: Churakova (Sidorova), O., Fonti, M., Siegwolf, R., Trushkina, T., Vaganov, E., and Saurer, M.: Response of Siberian trees to climatic changes over the past 1500 years, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-8995, https://doi.org/10.5194/egusphere-egu21-8995, 2021.

How to cite: Kim, W. M., Blender, R., Sigl, M., Raible, C., and Messmer, M.: Statistical Characteristics of Global Daily Extreme Precipitation during the last 3350 years (1501BC – 1849 AD) , EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-5405, https://doi.org/10.5194/egusphere-egu21-5405, 2021.

For more than a decade TEX₈₆ and U^K’₃₇, derived from ratios of biomarker lipids have widely been used as organic paleotemperature proxies. Yet, these proxies, especially TEX₈₆, have several uncertainties associated with factors such as depth and seasonal biases which are complicating its application as an annual mean sea-surface temperature (SST) proxy. To constrain this impact, we performed a relatively simple modelling exercise where we use instrumental temperature and nutrient data from 40 locations across the globe to predict theoretical proxy values and compare them with measured core-top proxy values.

The model first uses instrumental nutrient and temperature data, and probability density functions to predict the theoretical depth occurrence of the source organisms of the two proxies. Additionally, seasonal bias was introduced by predicting seasonal occurrences using instrumental nutrient and chlorophyll data. This was used to calculate the depth- and season weighed temperature signal annually deposited in the sediment, which in turn was converted to theoretical proxy values using culture or mesocosm calibrations. This showed, as expected, that depth and seasonal biases introduced scatter in the correlation between theoretical proxy values and annual mean SST but still highly significant for both U^K’₃₇ (r²= 0.96), and TEX₈₆ (r²= 0.77). We find that the theoretical proxy values are much lower than measured proxy value for TEX₈₆, which tentatively suggests that TEX₈₆ might in fact be coming from shallower depths or that the mesocosm calibration is incorrect. Our model for U^K’₃₇ results in theoretical values similar to measured values except for low temperature locations. This might suggest an influence of seasonal bias towards more warmer summer seasons which is more pronounced in high latitudes than in tropics.

How to cite: Varma, D., Reichart, G.-J., and Schouten, S.: Depth and seasonal biases in organic temperature proxies: a modelling study, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-1692, https://doi.org/10.5194/egusphere-egu21-1692, 2021.

How to cite: Cuesta-Valero, F. J., García-García, A., Beltrami, H., Zorita, E., and Jaume Santero, F.: Long-term Surface Temperature (LoST) Database as a Complement for Transient and Control Preindustrial Simulations, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-16519, https://doi.org/10.5194/egusphere-egu21-16519, 2021.

The Last Interglacial (LIG) period (~129,000–116,000 years BP) and the mid-Holocene (MH) (~6,000 years BP) are the two most recent intervals with temperatures comparable to low emissions scenarios for the end of the 21^st century. During the LIG and the MH differences in the seasonal and latitudinal distribution of insolation led to enhanced northern hemisphere high-latitude warmth relative to the pre-industrial, despite similar greenhouse gas concentrations, marking these intervals as potentially useful analogs for future change in regions like North America. Further, the inclusion of both LIG (127 ka) and MH (6 ka) experiments in the CMIP6-PMIP4 effort provides an opportunity to better understand the regional hydroclimate responses to radiative forcing during these two intervals. The dense coverage of paleoclimate proxy records for North America during the MH (N=260 sites) reveals a pattern of relative aridity in the Pacific Northwest and Western Canada and wetness in the southern Great Basin and Mexico. However, the seasonality and driving mechanisms of rainfall patterns across the continent remain poorly understood. Our understanding of terrestrial hydroclimate in North America during the LIG is more limited (N=39 sites), largely because the LIG is beyond the range of radiocarbon dating.

Here we present spatial comparisons between output from 14 PMIP4 global circulation models and LIG and MH networks of moisture-sensitive proxies compiled for the North American continent. We utilize two statistical measures of agreement – weighted Cohen’s Kappa and Gwet’s AC2 – to assess the degree of categorical agreement between moisture patterns produced by the models and the proxy networks for each time-slice. PMIP4 models produce variable precipitation anomalies relative to the pre-industrial for both the LIG and MH experiments, often disagreeing on both the sign and magnitude of precipitation changes across much of North America. The models showing the best agreement with the proxy network are similar but not identical for the two measures, with Gwet’s AC2 values tending to be larger than Cohen’s Kappa values for all models. This pattern is enhanced for the much larger MH proxy network and is likely related to the fact that Gwet’s AC2 is a more predictable statistic in the presence of high agreement. Overall agreement is lower for the mid-Holocene than for the LIG, reflecting smaller MH rainfall anomalies in the models. The models with the highest agreement scores during the LIG produce aridity in the Rocky Mountains and Pacific Northwest and wetness in Alaska, the Yukon, the Great Basin, and parts of the Mid-West and Eastern US, although spatial coverage of the proxies in these latter two regions is poor. The models with the highest agreement score for the mid-Holocene tend to produce aridity across Canada and the northern US with dry conditions extending down the US Pacific coast and increased wetness in the American Southeast and across the North American Monsoon region. Our analyses help elucidate the driving mechanisms of rainfall patterns during past warm states and can inform which models may be the most useful for predictions of near-future hydroclimate change across North America.

How to cite: de Wet, C., Oster, J., Ibarra, D., and Belanger, B.: North American rainfall patterns during past warm states: A proxy network-model comparison for the Last Interglacial and the mid-Holocene, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-6574, https://doi.org/10.5194/egusphere-egu21-6574, 2021.

The Atacama Desert in Northern Chile is considered to be the driest desert on Earth. At present, the annual rainfall amount is less than 1mm for parts of the hyper-arid core of the desert. The processes controlling this hyper-aridity are known, but the mean state and variability of the regional climate on geological time scales is not well understood. In this study, we aim to analyse climate conditions in the Atacama Desert from PMIP4 simulations. Our focus is on the Last Glacial Maximum (LGM), when climate records from the Central Atacama point to a substantially different climate with wetter conditions than at present (Diederich et al., 2020). We statistically analyse and evaluate PMIP4 historical simulations with respect to circulation patterns over the Southeast Pacific and Western South America which are associated with rare rainfall events in the Atacama Desert. For the evaluation, PMIP4 simulations for the historical period are compared to Reanalysis data, and we will focus on troughs and cutoff lows over the subtropical Southeast Pacific, and on the Bolivian High (Reyers et al., 2020). We then assess changes of the characteristics, e.g., the frequency of occurrence, of such circulation patterns for Paleo-climate conditions compared to the present. In the framework of our study, we perform km-scale simulations with the regional climate model WRF, using results from PMIP4 experiments for the historical period and for the LGM as boundary conditions. In the future, these simulations will be used to better understand the meso-scale processes, e.g., involved in local wind systems, that contribute to changes in the hydrological cycle and potentially impact the dust-emission activity of the desert. This study is part of the Collaborative Research Centre 1211 “Earth- Evolution at the dry Limit” (https://sfb1211.uni-koeln.de/).

Diederich, J,  Wennrich, V, Bao, R, and co-authors (2020). A 68 ka precipitation record from the hyperarid core of the Atacama Desert in northern Chile. Global and Planetary Change, 184, 103054. DOI:10.1016/j.gloplacha.2019.103054.

Reyers, M, Boehm, C, Knarr, L, Shao, Y, and Crewell, S (2020). Synoptic-to-Regional-Scale Analysis of Rainfall in the Atacama Desert (18°–26°S) Using a Long-Term Simulation with WRF. Monthly Weather Review, 149, 91-112. DOI:10.1175/MWR-D-20-0038.1.

How to cite: Reyers, M., Fiedler, S., and Shao, Y.: Paleo-climate shifts in the Atacama Desert from PMIP4 simulations, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-3185, https://doi.org/10.5194/egusphere-egu21-3185, 2021.

During the last million years, northern Africa has alternated between arid and humid conditions, as recorded by different kinds of climate archives, including fossil pollen, lake sediments, marine sediments and archaeological remains. Variations occur at millennial scale, with dry phases being similar to the current desert state in the region, and with wet phases, known as African Humid Periods (AHPs), characterised by a strong summer monsoon which can carry enough moisture inland to support rivers, lakes and lush vegetation further north than seen today. Recent sediment records from the Mediterranean Sea revealed that the previous five AHPs had different intensities, in relation to rainfall and vegetation extent. Motivated by these findings, our work focuses on explaining what caused such differences in intensity. To this end, we use the CLIMBER-2 climate model to study the AHP response to changes in three drivers of atmospheric dynamics: Earth's orbit variations, atmospheric concentration of CO₂ and inland ice extent. Global transient simulations of the last 190,000 years are used in new factorisation analyses, which allow us to separate the individual contributions of the forcings to the AHP intensity, as well as those of their synergies. We confirm the predominant role of the orbital forcing in the strength of the last five AHPs, and our simulations agree with previous estimates of a threshold in orbital forcing above which an AHP develops. Moreover, we show that atmospheric CO₂ and the extent of ice sheets can also add up to be as important as the orbital parameters. High values of CO₂, past a 205 ppm threshold, and low values of ice sheets extent, below an 8 % of global land surface threshold, yield the AHPs with the most precipitation and vegetation. Additionally, our results show that AHPs differ not only in amplitude, but also in their speed of change, and we find that the non-linear vegetation response of AHPs does not correlate with a single forcing and that the vegetation growth response is faster than its subsequent decline. In regards to future change, an extension of the simulations until the next 50,000 years, shows CO₂ to be the main driver of AHPs, with orbital forcing only setting the pace and their intensities being scenario-dependent.

How to cite: Duque-Villegas, M., Claussen, M., Brovkin, V., and Kleinen, T.: External and internal forcing of African Humid Periods from MIS 6 to MIS 1, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-12167, https://doi.org/10.5194/egusphere-egu21-12167, 2021.

During the mid-Holocene, an expansion of vegetation, lakes and wetlands over North Africa reinforced the West African monsoon precipitation increase that was initiated by changes in the orbital forcing. Sedimentary records reflect these surface changes, however, they provide only limited spatial and temporal information about the size and distribution of mid-Holocene lakes and wetlands. Previous simulation studies that investigated the influence of mid-Holocene lakes and wetlands on the West African monsoon precipitation, prescribed either a small lake and wetland extent or focusing on mega-lakes only. In contrast to these simulation studies, we investigate the range of simulated West African monsoon precipitation changes caused by a small and a potential maximum lake and wetland extent during the mid-Holocene.

Therefore, four mid-Holocene sensitivity experiments are conducted using the atmosphere model ICON-A and the land model JSBACH4 at 160 km resolution. The simulations have a 30-year evaluation period and only differ in their lake and wetland extent over North Africa: (1) pre-industrial lakes, (2) small lake extent, (3) maximum lake extent and (4) maximum wetland extent. The small lake extent is given by the reconstruction map of Hoelzmann et al. (1998) and the potential maximum lake and wetland extent is given by a model derived map of Tegen et al. (2002).

The simulation results reveal that the maximum lake extent shifts the Sahel precipitation threshold (> 200 mm/year) about 3 ° further northward than the small lake extent. The major precipitation differences between the small and maximum lake extent results from the lakes over the West Sahara. Additionally, the maximum wetland extent causes a stronger West African monsoon precipitation increase than the equally large maximum lake extent, particularly at higher latitudes.

How to cite: Specht, N., Claußen, M., and Kleinen, T.: Influence of a small and maximum lake and wetland extent on the simulated West African monsoon precipitation during the mid-Holocene, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-15050, https://doi.org/10.5194/egusphere-egu21-15050, 2021.

We give an overview on the global change in mid-to late Holocene vegetation pattern derived from a transient MPI-ESM1.2 simulation and discuss the vegetation trend in the context of the simulated Holocene climate change. The model captures the main trends found in reconstructions. Most prominent are the southward retreat of the northern treeline, coinciding the strong reduction of forest cover in the high northern latitudes during the Holocene, and the vast increase of the Sahara desert that is embedded in a general decrease and equator-ward retreat of the vegetation in the northern hemispheric monsoon margin regions. In contrast, large parts of the extratropical North American continent experience a greening during the Holocene, caused by an increase in forest and grass cover.

While the broad forest decline in the high northern latitudes can mainly be explained by the cooling of the warm season climate, precipitation is the driving factor for the tropical and extratropical vegetation trends on the northern hemisphere south of 60°N. The model indicates that most of the changes in rainfall can be related to the weakening of the northern hemispheric monsoon systems and the response of the global atmospheric circulation to this weakening.

The southern hemisphere is less affected by changes in total vegetation cover during the last 8000 years, but the monsoon related increase in precipitation and the insolation-induced cooling of the winter climate lead to shifts in the vegetation composition, mainly in between the woody plant functional types (PFTs).

The simulated large-scale global vegetation pattern almost linearly follow the subtle, approximately linear orbital forcing. Non-linear and more rapid changes in vegetation cover occur only on a regional level. The most striking area is the western Sahel-Sahara domain that experiences a rapid vegetation decline to a rather desertic state, in line with a strong decrease in moisture availability. The model also indicates rapid shifts in the vegetation composition in some regions in the high northern latitudes, in South Asia and in the monsoon margins of the southern hemisphere. These rapid transitions are mainly triggered by changes in the winter temperatures, which go into, or move out of, the bioclimatic tolerance range of the individual PFTs defined in the model and therefore have to be interpreted differently.

In summary, our model results identify the global monsoon system as the key player in Holocene climate and vegetation history and point to a far greater importance of the monsoon systems on the extra-monsoonal regions than previously assumed.

How to cite: Dallmeyer, A., Claussen, M., and Herzschuh, U.: Dominant role of the global monsoon intensity on large-scale Holocene vegetation transitions, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-10640, https://doi.org/10.5194/egusphere-egu21-10640, 2021.

We carry out three sets of last interglacial (LIG) experiments, named lig127k, and of pre-industrial experiments, named piControl, both as part of PMIP4/CMIP6 using three versions of the MIROC model: MIROC4m, MIROC4m-LPJ, and MIROC-ES2L. The results are compared with reconstructions from climate proxy data. All models show summer warming over northern high-latitude land, reflecting the differences between the distributions of the LIG and present-day solar irradiance. Globally averaged temperature changes are −0.94 K (MIROC4m), −0.39 K (MIROC4m-LPJ), and −0.43 K (MIROC-ES2L).
Only MIROC4m-LPJ, which includes dynamical vegetation feedback from the change in vegetation distribution, shows annual mean warming signals at northern high latitudes, as indicated by proxy data. In contrast, the latest Earth system model (ESM) of MIROC, MIROC-ES2L, which considers only a partial vegetation effect through the leaf area index, shows no change or even annual cooling over large parts of the Northern Hemisphere. Results from the series of experiments show that the inclusion of full vegetation feedback is necessary for the reproduction of the strong annual warming over land at northern high latitudes. The LIG experimental results show that the warming predicted by models is still underestimated, even with dynamical vegetation, compared to reconstructions from proxy data, suggesting that further investigation and improvement to the climate feedback mechanism are needed.

How to cite: O'ishi, R., Chan, W.-L., Abe-Ouchi, A., Sherriff-Tadano, S., Ohgaito, R., and Yoshimori, M.: PMIP4/CMIP6 last interglacial simulations using three different versions of MIROC: importance of vegetation, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-3792, https://doi.org/10.5194/egusphere-egu21-3792, 2021.

The climate evolution of the past few thousand years is essential for understanding the context in which civilisation arose and for understanding the natural background of anthropogenic influence. Proxy-inferred records show a complex picture of earlier warming and later cooling during the Holocene depending on region and reconstruction method. In contrast climate model simulations almost uniformly show warming throughout the past 10,000 years and for example also fail to reproduce a major advance of rainbelt over the Sahara.  These discrepancies raise questions about the reliability of climate models on longer-time scales.

We present a suite of four new transient Holocene simulations covering the last 8500 years using the HadCM3B-M21aD coupled general circulation. We use an optimised version of this model which is able to replicate the greening of the Sahara through changes to the atmospheric convection and vegetation schemes. We apply transient changes in Earth’s orbit, ice-sheets and sea-level and greenhouse gases, and optionally solar output, volcanic eruptions and anthropogenic land-use change.  The simulations without land-use show a warming throughout the Holocene, albeit with significantly higher variability once volcanic eruptions are included. With the inclusion of land-use change temperature trends in Northern Hemisphere are reversed from around 4000 years before present.

We explore the contribution of different forcings to the regional trends in the model ensemble and we compare the simulations against the Holocene reconstructions to evaluate the relative importance of each forcing. We also use the model ensemble to quantify the terrestrial coverage of proxy locations that is required to reliably infer global mean temperature variations.

How to cite: Hopcroft, P. and Valdes, P.: Contribution of forcings to Holocene climate evolution, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-4304, https://doi.org/10.5194/egusphere-egu21-4304, 2021.

How to cite: Wirths, C., Ziegler, E., Toohey, M., Schindlbeck-Belo, J. C., Kutterolf, S., Anders, H., and Rehfeld, K.: Comparing temperature trends and variability over the Holocene in climate models of low and high complexity, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-1842, https://doi.org/10.5194/egusphere-egu21-1842, 2021.

Multiple evidences from the analysis of satellite, in-situ and proxy data show that the climate is already changing toward a warmer Earth System due to our emissions of CO2 into the atmosphere. However, the magnitude and the extent of changes remain difficult to predict. A change in the ocean thermohaline circulation and its consequences for climate, such as drought, regional sea-level and ocean carbon uptake remain under debate as this circulation has been long thought to be stable during warm Earth periods – Interglacials. However, recent high-resolution reconstructions of carbon isotopes (δ¹³C) from the deep North Atlantic challenge this idea of stability and point toward abrupt modifications in the ocean interior biogeochemistry and/or ocean thermohaline circulation during the Last Interglacial (LIG, 125ka – 115ka).

Our model simulation of the LIG reproduces the observed magnitude and timescale of the reconstructed variations of δ¹³C, highlighting crucial dynamical changes in two regions of the North Atlantic deep-water formation (south of Greenland and south of Svalbard). These regions are found to drive the variations in the strength of the Atlantic Overturning Circulation (AMOC) when the Arctic sea-ice extent is perturbed.

Our study suggests that the AMOC may have experienced great instability phase during some parts of the LIG. The water mass geometry reorganization from the warm onset at 125ka to the glacial inception at 115ka could also have greatly impacted the distribution of carbon in the interior Ocean. Changes in sea-ice cover either south of Svalbard or in the Southern Ocean seem to play a determining role. However, in our global warming context, our study suggests that the mechanisms responsible for the LIG AMOC instability of the LIG may not occur by the end of the century if the Arctic sea-ice retreats from the high latitudes of the North Atlantic as projected by climate models.

How to cite: Kessler, A., Roche, D., Galaasen, E., Tjiputra, J., Bouttes, N., and Ninnemann, U.: AMOC instability during the Last Inerglacial, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-12364, https://doi.org/10.5194/egusphere-egu21-12364, 2021.

To better understand the processes contributing to future climate change, palaeoclimate model simulations are an important tool because they allow testing of the models’ ability to simulate very different climates than that of today.  As part of CMIP6/PMIP4, the latest version of the UK’s physical climate model, HadGEM3-GC31-LL (hereafter, for brevity, HadGEM3), was recently used to simulate the mid-Holocene (~6 ka) and Last Interglacial (~127 ka) simulations and the results were compared to the preindustrial era, previous versions of the same model and proxy data (see Williams et al. 2020, Climate of the Past).  Here, we use the same model to go further back in time, presenting the results from the mid-Pliocene Warm Period (~3.3 to 3 ma, hereafter the “Pliocene” for brevity).  This period is of particular interest when it comes to projections of future climate change under various scenarios of CO₂ emissions, because it is the most recent time in Earth’s history when CO₂ levels were roughly equivalent to today.  In response, albeit due to slower mechanisms than today’s anthropogenic fossil fuel driven-change, during the Pliocene global mean temperatures were 2-3°C higher than today, more so at the poles.

Here, we present results from the HadGEM3 Pliocene simulation.  The model is responding to the Pliocene boundary conditions in a manner consistent with current understanding and existing literature.  When compared to the preindustrial era, global mean temperatures are currently ~5°C higher, with the majority of warming coming from high latitudes due to polar amplification from a lack of sea ice.  Relative to other models within the Pliocene Modelling Intercomparison Project (PlioMIP), this is the 2^nd warmest model, with the majority of others only showing up to a 4.5°C increase and many a lot less.  This is consistent with the relatively high sensitivity of HadGEM3, relative to other CMIP6-class models.  When compared to a previous generation of the same UK model, HadCM3, similar patterns of both surface temperature and precipitation changes are shown (relative to preindustrial).  Moreover, when the simulations are compared to proxy data, the results suggest that the HadGEM3 Pliocene simulation is closer to the reconstructions than its predecessor.

How to cite: Williams, C., Lunt, D., Sellar, A., Roberts, W., Smith, R., Hopcroft, P., and Stone, E.: Simulation of the mid-Pliocene Warm Period using HadGEM3-GC31-LL: Pliocene climate relative to the pre-industrial era, previous model versions, other climate models and proxy data, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-7707, https://doi.org/10.5194/egusphere-egu21-7707, 2021.

Earth System Models (ESMs) project drying of the northern subtropics by the end of the 21^st century. However, geologic evidence from intervals with elevated concentrations of atmospheric carbon dioxide (pCO₂), like the mid-Pliocene, suggest mesic subtropical conditions. Several hypotheses, including an El Niño-like SST pattern and weaker Hadley circulation, have been proposed to explain this mismatch. Here, we show that PlioMIP2 ensemble broadly capture the pattern of proxy reconstructed Pliocene hydroclimate, notably a wetter Sahel and southeast Asia. Sensitivity simulations reveal that this pattern is driven by summertime rainfall increases as a result of lowered albedo and a distinct surface warming pattern, generated by prescribed vegetation and ice sheet changes. The resultant tropospheric moistening and stationary wave pattern enhance moisture convergence into the northern subtropics. Our results suggest that mid-Pliocene hydroclimate is part of the Earth system feedback to sustained CO₂ concentrations similar to today.

How to cite: Feng, R., Bhattacharya, T., Otto-bliesner, B., and Brady, E. and the PlioMIP2: Mid-Pliocene mesic subtropical hydroclimate over continents driven by land surface changes, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-13907, https://doi.org/10.5194/egusphere-egu21-13907, 2021.

The mid-Pliocene climate is the most recent geological period with a greenhouse gas concentration of approximately 400 ppmv, similar to the present day. Proxy reconstructions indicate enhanced warming in the high North Atlantic in the mid-Pliocene, which has been suggested to be a response to a stronger Atlantic Meridional Overturning Circulation (AMOC). PlioMIP2 ensemble results show a stronger AMOC and simulated North Atlantic sea surface temperatures (SSTs) match reconstructions better than PlioMIP1. A major difference between PlioMIP1 and PlioMIP2 is the closure of the Bering Strait and Canadian Archipelago in the Pliocene. Previous studies have shown that closure of these Arctic gateways leads to an enhanced AMOC due to altered freshwater fluxes in the Arctic.

Analysis of our Community Earth System Model (CESM1) simulations shows that the simulated increase in North Atlantic SSTs and strengthened AMOC in the Pliocene is a result of Pliocene boundary conditions rather than CO₂ concentration increase. Here we compare results from two runs with pre-industrial boundary conditions and 280 and 560 ppmv CO₂ concentrations and three runs with PlioMIP2 boundary conditions and 280, 400 and 560 ppmv CO₂ concentrations. Results show a 10-15% stronger AMOC in the Pliocene simulations as well as enhanced warming and saltening of the North Atlantic sea surface. While there is a stronger AMOC, the Atlantic northward ocean heat transport (OHT) in the Pliocene simulations only increases 0-3% with respect to the pre-industrial. Analysis indicates there is an altered relationship between the AMOC and OHT in the Pliocene, pointing to fundamentally different behavior of the AMOC in the Pliocene simulations. This is supported by a specific spatial pattern of deep water formation (DWF) areas in the Pliocene simulations that is significantly different from that of the pre-industrial. In the Pliocene simulations, DWF areas adjacent to south Greenland disappear and new DWF areas appear further southwards in the Labrador Sea off the coast of Newfounland. These results indicate that insight into the effect of the palaeogeographic boundary conditions is crucial to understanding the Pliocene climate and its potential as a geological equivalent to a future greenhouse climate.

How to cite: Weiffenbach, J., Baatsen, M., and von der Heydt, A.: Impact of Arctic gateways closure on the Atlantic Meridional Overturning Circulation in the Pliocene, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-9760, https://doi.org/10.5194/egusphere-egu21-9760, 2021.

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.

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.

How to cite: Nooteboom, P., Baatsen, M., Bijl, P., van Sebille, E., Sluijs, A., Dijkstra, H., and von der Heydt, A.: Model-data comparison in a strongly eddying Eocene ocean, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-8383, https://doi.org/10.5194/egusphere-egu21-8383, 2021.

Early Eocene Climatic Optimum (EECO, ~53-51 million years) is one of the past warm periods, associated with high CO₂ concentrations (~900-2500 ppmv), which can serve as an analogue for our possible future, high C0₂ climate. One notable feature of this hothouse climate state is the weaker meridional temperature gradient relative to pre-industrial values. This have been confirmed by both proxies and models, but the extent of the temperature gradient still requires more research. Models are challenged to reproduce the stronger than present day polar amplification signal, and it is also shown that high latitude proxy data are often influenced by seasonal bias. Thus, there is an uncertainty regarding both the observed and modelled meridional gradient and the mentioned issues complicate also the comparison between modeled and proxy data.

In our work we aim to investigate the EECO period with a simple energy balance box model and apply the maximum entropy production principle to explore the possible scenarios of meridional temperature gradients. We find that the maximum entropy production principle could be beneficial in the paleoclimate context since it has the utility to give an accurate prediction for non-equilibrium systems with the minimal amount of information. We also assess the heat transport signaled by proxy data and by state-of-the-art model outputs in accordance to our theoretical constrains based on the idealized test case.

How to cite: Kelemen, F. D. and Ahrens, B.: Exploring the possible meridional temperature gradient of Early Eocene Climatic Optimum with an energy balance model, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-7365, https://doi.org/10.5194/egusphere-egu21-7365, 2021.

Past climates contain precious information about the workings of the climate system, and about what can be expected in a changed climate. The Last Interglacial (LIG; ca. 125,000 years ago) is the most recent period of climate warmer than modern, at least in the Northern Hemisphere. Because of this, it has been often proposed that the LIG holds a partial analogy with a future warmer climate forced by enhanced greenhouse effect. Still, such analogy has never been examined in a quantitative manner. Here we address the question: for which scenario, time horizon, regions and season is the climate of the LIG a useful analogue of the future? We use the results of 13 climate models that performed the standard experiments of PMIP4 and CMIP6, and present a comparison of hemispheric temperature and precipitation between the LIG and SSP scenarios of the future. We also two independent assessments of models performance, by comparing their temperature and precipitation to climate reanalysis of the last decades and to proxies of the LIG. Insights gained from this comparison can inform studies in disciplines beyond climate studies, such as hydrology and ecology.

How to cite: Scussolini, P., Bakker, P., De Luca, P., Coumou, D., Bosmans, J., Lohmann, G., Thomas, Z., Turney, C., Menviel, L., Obase, T., Abe-Ouchi, A., Braconnot, P., Otto-Bliesner, B., Yin, Q., Prange, M., Tzedakis, C., Capron, E., Renssen, H., Ward, P., and Aerts, J.: Global temperature and hydroclimate in warmer climates of the past and future: the Last Interglacial versus greenhouse scenarios, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-15417, https://doi.org/10.5194/egusphere-egu21-15417, 2021.

We present new insights from the project PalMod, which started in 2016 and is envisioned to run for a decade. The modelling initiative PalMod aims at filling the long-standing scientific gaps in our understanding of the dynamics and variability of the climate system during the last glacial-interglacial cycle. One of the grand challenges in this context is to quantify the processes that determine the spectrum of climate variability on timescales that range from seasons to millennia. Climatic processes are intimately coupled across these timescales. Understanding variability at any one timescale requires understanding of the whole spectrum. If we could successfully simulate the spectrum of climate variability during the last glacial cycle in Earth system models, would this enable us to more reliably assess the future climate change? Such simulations are necessary to deduce, for example, if a regime shift in climate variability could occur during the next centuries and millennia in response to global warming. PalMod is specifically designed to enhance our understanding of the Earth system dynamics and its variability on timescales up to the multimillennial with complex Earth System Models.

The following major goals were achieved up to now:

• Full coupling of atmosphere, ocean and ice-sheet models, enabling investigation of Heinrich Events and bi-stability of the AMOC, and millennial-scale transient climate-ice sheet simulations.
• Implementation of a coupled ocean and land biogeochemistry enabling simulations with prognostic atmospheric CO₂ concentrations and including improved representation of methane (CH₄) in transient deglaciation runs.
• Systematic comparison of newly compiled proxy data with model simulations.

The major goal for the next two years is to set up the fully coupled physical-biogeochemical model which will be tested for three time periods: deglaciation, glacial inception and Marine Isotope Stage 3 (MIS3). This fully coupled model will be eventually used to simulate the complete glacial cycle and project the climate over the next few millennia.

How to cite: Fieg, K., Latif, M., Schulz, M., and Ilyina, T.: From the last interglacial to the future – new insights from modeling the last glacial-interglacial cycle in PalMod, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-11083, https://doi.org/10.5194/egusphere-egu21-11083, 2021.

We propose a reduced-complexity process-based model for the long-term evolution of the global ice volume, atmospheric CO₂ concentration and global mean temperature. The model only external forcings are the orbital forcing and anthropogenic CO₂ cumulative emissions. The model consists of a system of three coupled non-linear differential equations, representing physical mechanisms relevant for the evolution of the climate – ice sheets – Carbon cycle system in timescales longer than thousands of years. The model is successful in reproducing the glacial-interglacial fluctuations of the last 800 kyr, in good agreement with paleorecords both in terms of timing and amplitude, with a correlation between modelled and paleo global ice volume of up to 0.86.

Using different model realisations, we generate a probabilistic forecast of the evolution of the Earth system over the next 1 million years under natural and several fossil-fuel CO₂ release scenarios. In the natural scenario, the model assigns high probability of occurrence of long interglacials in the periods between present and 50 kyr after present, and between 400 kyr and 500 kyr after present. The next full glacial conditions are most likely to occur 90 kyr after present. The model shows that even already achieved cumulative CO₂ anthropogenic emissions (500 PgC) are capable of affecting the climate evolution for up to half million years, indicating that the beginning of the next glaciation is highly unlikely in the next 150 kyr. If cumulative fossil-fuel CO₂ emissions reach 3000 PgC, or higher, the model predicts with high probability ice-free Northern Hemisphere landmass conditions will prevail in the next half million years, postponing the natural occurrence of the next glacial inception to 600 kyr after present.

How to cite: Talento, S. and Ganopolski, A.: Evolution of the climate in the next million years: A reduced-complexity model for glacial cycles and impact of fossil fuel CO2, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-8512, https://doi.org/10.5194/egusphere-egu21-8512, 2021.

Between the last glacial maximum (LGM) and preindustrial times (PI), the atmospheric concentration of CH₄, as shown by reconstructions from ice cores, roughly doubled. It then doubled again from PI to the present. Ice cores, however, cannot tell us how that development will continue in the future, and ice cores also cannot shed light on the causes of the rise in methane, as well as the rapid fluctuations during periods such as the Bolling-Allerod and Younger Dryas.

We use a methane-enabled version of MPI-ESM, the Max Planck Institute for Meteorology Earth System Model, to investigate changes in methane cycling in a transient ESM experiment from the LGM to the present, continuing onwards into the future for the next millennium. The model is driven by prescribed orbit, greenhouse gases and ice sheets, with all other changes to the climate system determined internally. Methane cycling is modelled by modules representing the atmospheric transport and sink of methane, as well as terrestrial sources and sinks from soils, termites, and fires. Thus, the full natural methane cycle – with the exception of geological and animal emissions – is represented in the model. For historical and future climate, anthropogenic emissions of methane are considered, too.

We show that the methane increase since the LGM is largely driven by source changes, with LGM emissions substantially reduced in comparison to the early Holocene and preindustrial states due to lower temperature, CO₂, and soil carbon. Depending on the future climate scenario, these dependencies then lead to further increases in CH₄, with a further doubling of atmospheric CH₄ easily possible if one of the higher radiative forcing scenarios is followed. Furthermore, the future increases in CH₄ will persist for a long time, as CH₄ only decreases when the climate system cools again.

How to cite: Kleinen, T., Gromov, S., Steil, B., and Brovkin, V.: Methane in the climate system -- from the last glacial to the future, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-12440, https://doi.org/10.5194/egusphere-egu21-12440, 2021.

Much of the inter-model spread in equilibrium climate sensitivity (ECS) estimates is attributed to cloud and convective parameterizations which model cloud and water vapor feedbacks. These parameterizations also directly influence water isotopes, which may be retrieved not only from modern observations, but also a plethora of paleoclimate archives that represent a much broader range of variability than is available in modern measurements. And thus, these water isotope tracers can be used to constrain ECS by flagging unrealistic parts of the parameterization phase space via model biases in a perturbed parameterization ensemble (PPE) of paleoclimate simulations. In this proof-of-concept study, we evaluate a suite of isotope-enabled atmosphere-only GISS-E2.1 simulations, each with varying cloud and convective perturbations, against speleothem and ice core δ¹⁸O for the Last Glacial Maximum (LGM, 21000 years ago), mid-Holocene (MH, 6000 years ago) and pre-Industrial periods. The first-order spatial pattern of δ¹⁸O of precipitation (δ¹⁸O_p) is in excellent agreement between proxy data and all parameterizations across all time periods. While the simulations generally capture large scale δ¹⁸O_p patterns, the magnitude of change is consistently smaller in all simulations than those of the proxies, highlighting uncertainties in both models and proxies. Not a single set of parameterizations worked well in all climate states, indicating that improving future simulations requires determining all plausible parameter combinations critical in refining ECS. Further, it may be that certain parameterization choices represent certain types of variability better than others, and there may be a non-unique solution to ideal clouds/convection parameterization choices that is modulated by the question asked.

How to cite: Ramos, R. D., LeGrande, A. N., Griffiths, M. L., Elsaesser, G. S., Litchmore, D. T., Tierney, J. E., Pausata, F. S. R., and Nusbaumer, J.: Using paleoclimate data to constrain cloud parameterizations in GISS-E2.1, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-14136, https://doi.org/10.5194/egusphere-egu21-14136, 2021.

In simulations of the climate during the Last Glacial Maximum (LGM), we employ two different isotope-enabled atmospheric general circulation models (NCAR iCAM3 and MPI ECHAM6-wiso) and use simulated (by coupled climate models) as well as reconstructed (from a new global climatology of the ocean surface duing the LGM, GLOMAP) surface conditions.

The resulting atmospheric fields reflect the more pronounced structure and gradients in the reconstructions, for example, the precipitation is more depleted in oxygen-18 in the high latitudes and more enriched in low latitudes, especially in the tropical convective regions over the maritime continent in the equatorial Pacific and Indian Oceans and over the equatorial Atlantic Ocean. Furthermore, at the sites of ice cores and speleothems, the model-data fit improves in terms of the coefficients of determination and root-mean square errors.

In additional sensitivity experiments, we also use the climatologies by Annan and Hargreaves (2013) and Tierney et al. (2020) and consider the impact of changes in reconstructed sea-ice extent and the global-mean sea-surface temperature.

Our findings imply that the correct simulation or reconstruction of patterns and gradients in sea-surface conditions are crucial for a successful comparison to oxygen-isotope data from ice cores and speleothems.

How to cite: Paul, A., Cauquoin, A., Mulitza, S., Tharammal, T., and Werner, M.: Sensitivity of simulated oxygen isotopes in ice cores and speleothems to Last Glacial Maximum surface conditions, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-2305, https://doi.org/10.5194/egusphere-egu21-2305, 2021.

We present new reconstructions of global climatological temperature fields for the Last Glacial Maximum and the mid-Pliocene Warm Period.

The method is based on an Ensemble Kalman Smoother which combines globally complete modelled temperature fields, with sparse proxy-based estimates of local temperature anomalies. This ensures spatially coherent fields which respect physical principles and which are also tied closely to observational estimates. 

For the Last Glacial Maximum, we use the full set of PMIP2/3/4 model simulations, and we combine this with a wide range of proxy-based SST and SAT estimates of local temperature to ensure the best possible global coverage. Our reconstruction has a global mean surface air temperature anomaly of -5.3 +- 0.9C relative to the pre-industrial climate, and thus lies roughly half-way between the estimates of Annan and Hargreaves (2013) and Tierney et al (2020). We examine the reasons for these differences and discuss their implications.

For the mid-Pliocene Warm Period, we use the PlioMIP 1 and 2 model simulations and the PRISM proxy estimates for the 3.2 Ma time slice. These data are considerably more sparse and uncertain than for the LGM and our reconstruction is correspondingly more uncertain. We obtain an estimate of 5.6 +- 1.6C which is considerably warmer than most previous estimates, suggesting a significant discrepancy between the models and the data. We investigate the reasons for this and discuss the implications.

How to cite: Annan, J., Hargreaves, J., and Mauritsen, T.: Reconstructing the surface temperature fields of the Last Glacial Maximum and mid-Pliocene Warm Period using climate models and data., EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-11047, https://doi.org/10.5194/egusphere-egu21-11047, 2021.

Comparing simulations and data from paleoclimate archives such as speleothems can test the capability of climate models to capture past climate changes. In past, present, and future, the hydrologic response to radiative forcing changes is far less understood and more uncertain than thermal changes.
 
Speleothems store terrestrial climate information in the form of isotopic oxygen in mineral and are found mostly in the low-to mid-latitudes of the landmasses. Their usually well preserved (semi-)continuous time series of oxygen isotope ratio δ¹⁸O can cover full Glacial-Interglacial cycles and are used for past climate reconstructions. However, the measured δ¹⁸O in the mineral is influenced by multiple climate and cave-related variables and does, therefore, not directly represent past temperature or precipitation. 

We assess the capability of the isotope-enabled models HadCM3 and ECHAM5-MPI/OM to simulate decadal to centennial climate variability beyond the instrumental period. In particular, we investigate the relationship between simulated δ¹⁸O and precipitation variability under different background conditions. By comparing simulated δ¹⁸O values at cave locations to the large global speleothem database SISALv2 (Comas-Bru et al. 2020), we also examine the consistency between modeled and archived temporal changes in δ¹⁸O in the mean state and variability. Our strategy involves forward-modeling of cave processes such as temperature-dependent fractionation and transit times to constrain a simple speleothem proxy model for the simulation output. For the late Holocene, we observe a strongly underestimated simulated isotopic variability on decadal to centennial timescales. We further test how much this underestimation depends on the background radiative forcing conditions by comparing the Last Glacial Maximum, the mid-Holocene, and the late Holocene. This provides deeper insight on low to mid-latitude state-dependent climate variability on decadal to centennial time scales. 

Reference:

Comas-Bru, L., et. al. SISALv2: a comprehensive speleothem isotope database with multiple age-depth models. Earth System Science Data 12, 2579-2606 (2020) https://essd.copernicus.org/articles/12/2579/2020/

How to cite: Buehler, J., Weitzel, N., Baudouin, J.-P., Werner, M., and Rehfeld, K.: Last Glacial to present-day variability of surface climate from oxygen isotope signatures in speleothems and model simulations, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-9847, https://doi.org/10.5194/egusphere-egu21-9847, 2021.

Experiment outputs are now available from the Coupled Model Intercomparison Project’s 6^th phase (CMIP6) and the past climate experiments defined in the Model Intercomparison Project’s 4^th phase (PMIP4). All of this output is freely available from the Earth System Grid Federation (ESGF). Yet there are overheads in analysing this resource that may prove complicated or prohibitive. Here we document the steps taken by ourselves to produce ensemble analyses covering past and future simulations. We outline the strategy used to curate, adjust the monthly calendar aggregation and process the information downloaded from the ESGF. The results of these steps were used to perform analysis for several of the initial publications arising from PMIP4. We provide post-processed fields for each simulation, such as climatologies and common measures of variability. Example scripts used to visualise and analyse these fields is provided for several important case studies.

How to cite: Zhao, A. and Brierley, C.: Workflow and tools to analyse the PMIP4-CMIP6 ensemble, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-2476, https://doi.org/10.5194/egusphere-egu21-2476, 2021.

How to cite: Weitzel, N., Andres, H., Baudouin, J.-P., Bothe, O., Dolman, A., Jonkers, L., Kapsch, M., Kleinen, T., May, M., Mikolajewicz, U., Paul, A., and Rehfeld, K.: Towards model-data comparison of the deglacial temperature evolution in space and time, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-9683, https://doi.org/10.5194/egusphere-egu21-9683, 2021.

Modelling studies are evaluated by comparing the simulation outputs to an observational reference. In climate science, the number and complexity of the models and the mass of data have led the community to develop standardised methods and automated tools, such as the Climate Variability Diagnostic Package or the ESMValTool. However, these tools are mostly designed to evaluate simulations of the instrumental period. Different methods are required to compare paleoclimate simulations to palaeodata. For example, new variables are being modelled, such as vegetation, ice sheet extent, or isotopic ratio, and are used for the evaluation. Changing boundary conditions in transient simulations further complicate the evaluation process: traditional indices that characterise circulation (e.g. monsoon) or modes of variability (e.g. NAO, ENSO) need to be adapted, while new ones are needed to investigate modes of longer timescale and abrupt events. Finally, the palaeodata also present challenges: various type of uncertainty, complex relation to climate variables, and different spatio-temporal representativeness compared to model outputs. Here, we summarise the challenges of model-data comparison in paleo-climate studies. We then review some of the different methods and tools already developed by the community, such as biome comparison and Bayesian approaches to quantify model-data deviation. We finally discuss the implementation of an evaluation framework which aims to provide both adaptable tools to the community and automated standardised analyses.

How to cite: Baudouin, J.-P., Bothe, O., Chevalier, M., Weitzel, N., Dallmeyer, A., Brierley, C., and Rehfeld, K.: Model-data comparison challenges in paleo-climate analyses: Towards an evaluation toolbox for transient climate model simulations , EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-5716, https://doi.org/10.5194/egusphere-egu21-5716, 2021.

How to cite: Ziegler, E., Andres, H., Ellerhoff, B., Kapsch, M.-L., Kutterolf, S., Mikolajewicz, U., Schindlbeck-Belo, J. C., Toohey, M., Wirths, C., Weitzel, N., and Rehfeld, K.: State-dependency of temperature variability in transient simulations of the last Deglaciation from models of varying complexity, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-11006, https://doi.org/10.5194/egusphere-egu21-11006, 2021.

Fossil pollen datasets can help to understand the temporal and spatial distribution patterns and driving forces of the past terrestrial biomes in high northern latitudes. Here we present a global pollen dataset since the Last Glacial Maximum, synthesized from 2821 palynological records from the Neotoma Paleoecology Database and additional literature. All terrestrial pollen taxa were taxonomically harmonized on genus (woody taxa) or family level (herb taxa) and temporally standardized by using a defined parameter setting for Bayesian age-depth modeling based on 14C dating. The age-depth models were statistically compared with existing models for each record. With a biomization approach, we reconstructed biomes for several time-slices throughout the last 22000 years with a temporal resolution of roughly 500 years. The reconstructed biome distributions are compared to simulated biome distributions inferred from a transient simulation for the last 25000 years, performed in the comprehensive Earth System Model of the Max Planck Institute (MPI-ESM). The overall biome trend agrees well, but the simulation shows lower forest cover in the high northern latitudes and reaches the maximum forest cover in the Holocene much earlier than the reconstructions indicate.

How to cite: Li, C., Dallmeyer, A., Böhmer, T., Postl, A., and Herzschuh, U.: Northern hemispheric biome changes synthesized from taxonomically harmonized and temporally standardized fossil pollen record since the Last Glacial Maximum in comparison to MPI-ESM simulations, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-12224, https://doi.org/10.5194/egusphere-egu21-12224, 2021.

Southern Ocean sea ice and oceanic fronts are known to play an important role on the climate system, carbon cycles, bottom ocean circulation, and Antarctic ice sheet. However, many models of the previous Past-climate Model Intercomparison Project (PMIP) underestimated sea-ice extent (SIE) for the Last Glacial Maximum (LGM)(Roche et al., 2012; Marzocchi and Jensen, 2017), mainly because of surface bias (Flato et al., 2013) that may have an impact on mean ocean temperature (MOT). Indeed, recent studies further suggest an important link between Southern Ocean sea ice and mean ocean temperature (Ferrari et al., 2014; Bereiter et al., 2018 among others). Misrepresent the Antarctic sea-ice extent could highly impact deep ocean circulation, the heat transport and thus the MOT. In this study, we will stress the relationship between the distribution of Antarctic sea-ice extent and the MOT through the analysis of the PMIP3 and PMIP4 exercise and by using a set of MIROC models. To date, the latest version of MIROC improve its representation of the LGM Antarctic sea-ice extent, affecting the deep circulation and the MOT distribution (Sherriff-Tadano et al., under review).

Our results show that available PMIP4 models have an overall improvement in term of LGM sea-ice extent compared to PMIP3, associated to colder deep and bottom ocean temperature. Focusing on MIROC (4m) models, we show that models accounting for Southern Ocean sea-surface temperature (SST) bias correction reproduce an Antarctic sea-ice extent, 2D-distribution, and seasonal amplitude in good agreement with proxy-based data. Finally, using PMIP-MIROC analyze, we show that it exists a relationship between the maximum SIE and the MOT, modulated by the Antarctic intermediate and bottom waters.

How to cite: Vadsaria, T., Sherriff-Tadano, S., Abe-Ouchi, A., Obase, T., Chan, W.-L., and Crosta, X.: Last Glacial Maximum Antarctic sea ice linked with global mean ocean temperature: evidence from PMIP3, PMIP4 and MIROC-4m simulations, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-3550, https://doi.org/10.5194/egusphere-egu21-3550, 2021.

The last deglaciation provides as unique a framework to investigate the processes of ice sheet and climate interaction during periods of mass loss as in the current climate. Here we simulate the Last Glacial Maximum (LGM) northern hemisphere ice sheets climate, surface mass balance (SMB), and dynamics with the Community Earth System Model version 2 (CESM2, Danabasoglu et al., 2020)) and the Community Ice Sheet Model version 2 (CISM2, Lipscomb et al., 2019). This LGM simulation will be later used as starting point for coupled CESM2-CISM2 simulations of the last deglaciation.

CESM2 is run at the nominal resolution used for IPCC-type projections (approx. 1 degree for all components). The model includes an advanced snow/firn and SMB calculation (van Kampenhout et al, 2019; Sellevold et al, 2019) the land component (CLM, cite) that has been evaluated and applied to the simulation of the future Greenland melt (van Kampenhout et al, 2020, Muntjewerf et al., 2020a,b, Sellevold & Vizcaino, 2020).

Our analysis examines how the global, Arctic, and North Atlantic climate result in the simulated radiative and turbulent heat fluxes over the ice sheets, and the mass fluxes from precipitation, refreezing, runoff, and sublimation. We also examine the simulated ice streams in CISM2, which is run at 8 km under a higher-order approximation for ice flow.

How to cite: Bradley, S. L., Petrini, M., Sellevold, R., Vizcaino, M., Lipscomb, W. H., and Georgiou, S.: Simulation of LGM ice sheets with the Community Earth System Model version 2 , EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-15463, https://doi.org/10.5194/egusphere-egu21-15463, 2021.

The DO events of the last ice age represent one of the best studied abrupt climate transitions, yet we still lack a comprehensive explanation for them. There is uncertainty whether current IPCC-relevant models can effectively represent the processes that cause DO events. Current Earth system models (ESMs) seem overly stable against external perturbations and incapable of reproducing most abrupt climate changes of the past (Valdes, 2011). If this holds true, this could noticeably influence their capability to predict future abrupt transitions, with significant consequences for the delivery of precise climate change projections.  In this task, the objectives of this study are (1) to cross compare existing simulations that show spontaneous DO-type oscillations using a common set of diagnostics so we can compare the mechanisms and the characteristics of the oscillations, and (2) to formulate possible pathways to a DO PMIP protocol that could help investigate cold-period instabilities through a range of insolation-, freshwater-, GHG-, and NH ice sheet-related forcings, as well as evaluating the possibility of spontaneous internal oscillations.

Although most abrupt DO events happened during MIS3, only few studies investigate DO events in coupled general circulation models under MIS 3 conditions (e.g., Kawamura et al., 2017; Zhang and Prange, 2020). Here, we thus propose that the MIS3 period could be the focus of such a DO-event modelling protocol. More specific sensitivity experiments performed under MIS 3 boundary conditions are needed in order to (1) better understand the mechanisms behind millennial-scale climate variability, (2) explore AMOC variability under intermediate glacial conditions, and (3) help answer the question: “are models too stable?”.

How to cite: Malmierca-Vallet, I., Sime, L. C., and Valdes, P. J.: Possible pathways to a Dansgaard-Oeschger (DO) PMIP protocol, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-8730, https://doi.org/10.5194/egusphere-egu21-8730, 2021.

Understanding the processes causing variations in the carbon cycle is critical to accurately simulate the future carbon cycle and climate. Paleoclimate models can provide insights about these processes since they are used under different conditions than present-day’s and evaluated against paleoproxy data. In particular, the Last Glacial Maximum (LGM) has been a focus of the Paleoclimate Modelling Intercomparison Project (PMIP) as it is well-documented thanks to numerous paleoclimate archives. Around 21,000 years ago, the LGM was a colder period with extensive ice sheets in the Northern Hemisphere and a resulting lower sea-level. Although this period has been studied for years, the causes of the lower atmospheric CO₂ concentration at the time (around 186 ppm, against 280 ppm at the pre-industrial) remain unclear, and models struggle to simulate this low CO₂ value. The ocean is thought to have played a significant role due to different processes (through changes of the biological pump efficiency, ocean circulation, sea-ice, and CO₂ solubility due to colder temperatures), but no consensus has been reached yet as to their contribution (Khatiwala et al. [2019], Yu et al. [2016], Marzocchi and Jansen [2019]).

Despite the carbon cycle being simulated by more and more climate models, it has not been systematically analysed within the framework of PMIP multimodel comparisons. In this context, the ongoing PMIP-carbon project aims at comparing climate-carbon interactions in LGM simulations, and includes results from both intermediate complexity models and general circulation models. The PMIP protocol proposes standardized forcing parameters and boundary conditions (Kageyama et al. [2017]) and specifies a few recommendations for ocean biogeochemistry models (adjustment of salinity, dissolved inorganic carbon, alkalinity, and nutrients to account for the change in ocean volume). Indeed, the bathymetric changes associated with a sea-level drop of 133 m entail a change of the reservoir size and potential technical issues concerning the conservation of carbon.

In this study, we use outputs from PMIP-carbon models and other models available on the ESGF (MIROC4m-COCO, MIROC-ES2L, CESM, IPSL-CM5A2, UVic, LOVECLIM, iLOVECLIM, CLIMBER_2P ; GISS-E2-R, MRI-CGCM3, MPI-ESM-P, CNRM-CM5, MIROC-ESM) to compute total ocean volumes and compare them to high resolution topographic data (etopo1 for the PI, GLAC-1D and ICE-6G-C for the LGM). We show that the deglacial volume change is rarely accurate. We then use the iLOVECLIM model with a new bathymetry implementation method (Lhardy et al. [in review, 2020]) to demonstrate the effect of an improved ocean volume on the simulated oceanic carbon content, resulting in an increase of the already overestimated atmospheric CO₂ concentration. We also quantify the effect of the mentioned adjustments of salinity, alkalinity, and carbon. The results reinforce the idea that a realistic ocean volume is needed, as well as consistency between models in dealing with large changes in bathymetry.

How to cite: Lhardy, F., Bouttes, N., Roche, D., Abe-Ouchi, A., Chase, Z., Ivanovic, R., Jochum, M., Kageyama, M., Kobayashi, H., Menviel, L., Muglia, J., Nuterman, R., Oka, A., Schmittner, A., Vettoretti, G., and Yamamoto, A.: The impact of bathymetry on the simulated carbon at the Last Glacial Maximum, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-7297, https://doi.org/10.5194/egusphere-egu21-7297, 2021.

The last deglaciation is a time of large climate transition from a cold Last Glacial Maximum at 21,000 years BP with extensive ice sheets, to the warmer Holocene 9,000 years BP onwards with reduced ice sheets. Despite more and more proxy data documenting this transition, the evolution of climate is not fully understood and difficult to simulate. The PMIP4 protocol (Ivanovic et al., 2016) has indicated which boundary conditions to use in model simulations during this transition. The common boundary conditions should enable consistent multi model and model-data comparisons. While the greenhouse gas concentration evolution and orbital forcing are well known and easy to prescribe, the evolution of ice sheets is less well constrained and several choices can be made by modelling groups. First, two ice sheet reconstructions are available: ICE-6G (Peltier et al., 2015) and GLAC-1D (Tarasov et al., 2014). On top of topographic changes, it is left to modelling groups to decide whether to account for the associated bathymetry and land-sea mask changes, which is technically more demanding. These choices could potentially lead to differences in the climate evolution, making model comparisons more complicated.

We use the iLOVECLIM model of intermediate complexity (Goosse et al., 2010) to evaluate the impact of different ice sheet reconstructions and the effect of bathymetry changes on the global climate evolution during the Last deglaciation. We test the two ice sheet reconstructions (ICE-6G and GLAC-1D), and have implemented changes of bathymetry and land-sea mask. In addition, we also evaluate the impact of accounting for the Antarctic ice sheet evolution compared to the Northern ice sheets only.

We show that despite showing the same long-term changes, the two reconstructions lead to different evolutions. The bathymetry plays a role, although only few changes take place before ~14ka. Finally, the impact of the Antarctic ice sheet is important during the deglaciation and should not be neglected.

References

Goosse, H., et al., Description of the Earth system model of intermediate complexity LOVECLIM version 1.2, Geosci. Model Dev., 3, 603–633, https://doi.org/10.5194/gmd-3-603-2010, 2010

Ivanovic, R. F., et al., Transient climate simulations of the deglaciation 21–9 thousand years before present (version 1) – PMIP4 Core experiment design and boundary conditions, Geosci. Model Dev., 9, 2563–2587, https://doi.org/10.5194/gmd-9-2563-2016, 2016

Peltier, W. R., Argus, D. F., and Drummond, R., Space geodesy constrains ice age terminal deglaciation: The global ICE-6G_C (VM5a) model, J. Geophys. Res.-Sol. Ea., 120, 450–487, doi:10.1002/2014JB011176, 2015

Tarasov,L.,  et al., The global GLAC-1c deglaciation chronology, melwater pulse 1-a, and a question of missing ice, IGS Symposium on Contribution of Glaciers and Ice Sheets to Sea-Level Change, 2014

How to cite: Bouttes, N., Roche, D., Lhardy, F., Quiquet, A., Paillard, D., and Goosse, H.: Impact of ice sheet reconstructions on the deglacial climate in iLOVECLIM, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-8776, https://doi.org/10.5194/egusphere-egu21-8776, 2021.

How to cite: Romé, Y., Ivanovic, R., and Gregoire, L.: Meltwater driven abrupt climate changes in the North Atlantic simulated in a Heinrich Stadial 1 background, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-8724, https://doi.org/10.5194/egusphere-egu21-8724, 2021.

The global mean sea level in the last interglacial (LIG, about 130,000 to 115,000 years before present) was very likely higher than the present level, driven mainly by mass loss of the Antarctic ice sheet. Some studies have suggested that this mass loss may have been caused by the warmer temperature over the Southern Ocean in the LIG compared with the present interglacial. However, the ultimate cause of the difference in Antarctic warming between the last and current interglacials has not been explained. Here, based on transient simulations of the last deglaciation using a fully coupled ocean–atmosphere model, we show that greater meltwater (by a factor of 1.5 relative to the last deglaciation) during the middle and later stages of the deglaciation could have produced the difference in Antarctic warmth. Northern Hemisphere ice sheet model experiments suggest that the difference in meltwater was caused by slightly smaller orbital eccentricity in our current interglacial than in the LIG, indicating that mass loss of the Antarctic ice sheet is influenced by the preceding northern summer insolation and disintegration of Northern Hemisphere ice sheets.

How to cite: Obase, T., Abe-Ouchi, A., and Saito, F.: Antarctic warmth in the last interglacial driven by Northern insolation and deglaciation, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-8145, https://doi.org/10.5194/egusphere-egu21-8145, 2021.

A popular hypothesis in paleoclimatology posits that the episodic discharges of glacial water from the Laurentide Ice Sheet (LIS) to the North Atlantic caused abrupt changes in ocean circulation and climate during the last (de)glacial periods. Implicit in this hypothesis is that the glacial water spread away from the coast and reached critical sites of deep water formation. Among the processes that could favour the offshore export of glacial water released along the eastern North American coast is the entrainment with the Gulf Stream near Cape Hatteras, where the Stream is observed to detach from the coast in the modern climate, or at other locations between Cape Hatteras and the Grand Banks of Newfoundland.

Here we investigate the fate of glacial water released in the western North Atlantic from the Laurentian Channel, which geologic evidence suggests to have been the main route of ice discharge from the Québec-Labrador Ice Dome of the LIS. To this end, we conduct numerical experiments with an ocean circulation model with eddy-resolving resolution and configured to represent the region north of Bermuda and west of the Grand Banks. Experiments with different regional sea levels are performed which correspond to different estimates of global sea level since the Last Glacial Maximum. In each experiment, glacial water in liquid form is discharged from the Laurentian Channel, providing a paleoceanographic analogue of the dam-break problem. As expected from the action of the Coriolis force and from the properties of the glacial water inflow, the discharged water turns to the right of the Channel and then produces a narrow buoyant current that flows along the coast to the southwest towards Cape Hatteras. Our presentation will focus on the interaction of this current with the Gulf Stream, particularly with its meanders and rings, and on the role of this interaction both in the seaward export of glacial water and in the modification of the Stream itself.

How to cite: Marchal, O. and Condron, A.: Interaction of a Glacial Water Outflow with the Gulf Stream, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-6346, https://doi.org/10.5194/egusphere-egu21-6346, 2021.

Freshwater, in the form of glacial runoff, is hypothesized to play a critical role in centennial to millennial scale climate variability, eg. the Younger Dryas and Dansgaard Oeschger events. Freshwater injection, or hosing, model experiments demonstrate that freshwater has the capability to generate abrupt climate transitions.  However, in an attempt to mitigate the inability of most models to resolve the smaller-scale features relevant to freshwater transport (such as boundary currents and mesoscale eddies), these hosing experiments commonly apply the entirety of the freshwater directly to the regions of deepwater formation (DWF). Our results indicate that this can inflate the freshwater signal in those regions by as much as four times. We propose a novel method of freshwater injection for such low-resolution models that spatially distributes the freshwater in accord with the results of eddy-permitting modelling. Furthermore, this “freshwater fingerprint” method not only impacts the timing of simulated climate transitions but also can allow us to evaluate how much we are overestimating the effects of freshwater when injected directly into sites of DWF.

The freshwater fingerprints we develop are based on a suite of freshwater injection experiments performed using an eddy permitting Younger Dryas configuration of the MITGCM. Freshwater injection locations include the Mackenzie River, Gulf of St. Lawrence, Gulf of Mexico and a location off the coast of Norway, with flux amounts bounded by glacial reconstructions. These simulations indicate that freshwater from the Mackenzie River and Fennoscandia have the largest impact on salinity in most of the conventional sites of DWF (GIN and Labrador Seas, and in these simulations, predominantly the northern North Atlantic due to extensive sea ice), while freshwater from the Gulf of St. Lawrence is effective at freshening only the northern North Atlantic. The Gulf of Mexico has little impact on any DWF region we consider, mostly because the lower but continual flux in our simulations does not allow freshwater to penetrate northward past the Gulf Stream. The dilution of the freshwater signal as it is transported from the site of injection to the DWF zones leads to a reduction in the effective freshwater forcing, making hosing directly over DWF zones even with realistic freshwater amounts unrealistic. Thus, we construct freshwater fingerprints from these simulations by extracting the freshwater anomaly spatial pattern averaged over the last 5 simulation years, vertically integrating the field and normalizing it.

The freshwater fingerprint is then implemented in the COSMOS Earth Systems Model, which is run at resolutions typical for paleoclimate simulations (non-eddy permitting). Initial results show that freshwater from the Mackenzie River using our  fingerprint method leads to a more gradual cooling than if the meltwater is released directly over the hosing region (50-70N). We conclude that hosing over DWF zones, even with realistic freshwater amounts, produces an unrealistically large climate response. Additional results for the remaining injection locations and with the fingerprint implemented in a simpler climate model will be presented.

How to cite: Love, R., Andres, H., Condron, A., Zhang, X., Lohmann, G., and Tarasov, L.: Towards more physically constrained freshwater injection via eddy permitting simulations of the last glacial cycle, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-10357, https://doi.org/10.5194/egusphere-egu21-10357, 2021.

The Dead Sea depression features exceptionally well preserved lacustrine sedimentary sequences and fossil lake-level markers that attest to a much more extensive lake with a maximum highstand water level of more than 200 m above the modern Dead Sea. Lake-level reconstructions based on sedimentary sequences places this highstand phase within the interval of 15-29 ka. Regional paleoclimatic records, however, indicate arid conditions during this time. This apparent contradiction has been explained by spatially heterogeneous moisture delivery resulting from a southward shift of the Westerly wind system and a change in the path and intensity of winter storms. A newly established lake level-chronology based on ¹⁴C- and U/Th-dating of fossil stromatolites has provided contrasting results with respect to previous investigations. Accordingly, the paleolake-highstand was of much shorter duration and occurred at least 10 ka earlier than previously suggested. The new lake-level curve agrees with evidence of arid glacial and humid interglacial periods in the Levant.  In this study we compared these different lake-level reconstructions quantitatively, using a distributed hydrological balance model. This model computes evaporation based on an aerodynamic- /mass-transfer approach. Calibration and validation of this model is achieved by using ~30 years of pre-anthropogenic lake-level observations combined with interpolated climate surfaces based on weather-station records. In the paleo-hydrological reconstruction we account for parameter uncertainties using Monte-Carlo simulations. Our preliminary results show a pronounced sensitivity of the lake-level to precipitation, wind speed, and surface roughness.

How to cite: Mey, J., Jara, J., and Strecker, M. R.: Late Pleistocene paleo-hydrological reconstruction based on a new lake-level curve of the Dead Sea, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-13482, https://doi.org/10.5194/egusphere-egu21-13482, 2021.

The Bond 4 event starting at 7000 yr BP and culminating around 5200 yr BP corresponds to the largest (in magnitude and duration) invasion of drifting ice across the subpolar North Atlantic during the Holocene (Bond et al., 2001). While several studies have focused on other events of the Holocene, such as the 8.2 ka, the 4.2 ka and the Little Ice Age, little is known about the mid-Holocene 5.2 ka event. Here we present a global compilation of carefully selected high-resolution time series of sea surface temperature (SST; N=58) and humidity/precipitation (N=35) to characterize in space and time the 5.2 ka event pattern.

The SST records show the occurrence of cold conditions in the North Atlantic, western Mediterranean as well as in the western Pacific Ocean. However, they indicate warming in the high latitude North Atlantic, the southeastern Atlantic, the eastern Mediterranean and the Arabian and Red seas. Humidity/precipitation data (mainly based on oxygen  isotope records in speleothems) indicate dry conditions in the northern hemisphere subtropical and mid latitude regions of all continents. Based on these data and others from marine and lacustrine records in tropical regions, we suggest a possible weakening of monsoon systems, i.e. in Africa, North America, southwest Asia as well as East Asia. Precipitation reduced as well in most regions of the Mediterranean (i.e. except Iberian Peninsula). All together these data indicate severe climate conditions during the 5.2 ka event.

Based on the recent compilation of sortable silt from the high latitude North Atlantic of McCave and Andrews (2019), the 5.2 ka event coincides with a decrease of the main Shallow and bottom ocean flows (i.e. North Iceland Irminger Current, East Greenland Current, Iceland-Scotland overflow) probably reflecting a weakening of the North Atlantic Deep Water formation. This event also corresponds to the occurrence of several solar minima as well as several tropical volcanic mega-eruptions that could have triggered a global colder and drier climate (Steinhilber et al., 2012; Kobashi et al., 2017). 

How to cite: Jalali, B. and Sicre, M.-A.: The spatiotemporal pattern of the Bond 4 event (5.2 ka) : a global data-based review, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-15576, https://doi.org/10.5194/egusphere-egu21-15576, 2021.

In the winter of 2017 three undisturbed sediment cores were retrieved from the Szemes Basin of Lake Balaton. The sediments were sampled for AMS ¹⁴C dating and we used 8 of the radiocarbon dates for age-depth modelling. Based on this, the investigated sediment sequence covers the entire Holocene and Late Glacial period and the bottom of the sediment is ca. ~16,000 cal yr BP old. X-ray fluorescence spectrometry (XRF) was used to reconstruct rapid changes in the element content of the lake sediment. The evaluation of the measured results makes it possible to reconstruct the changes in the discharge environment and lake water level that can be related to the climate and human impact. Based on the data, two major evaporation events can be observed at 5500 BP and 8100 BP. These results were also verified by oxygen isotope studies. To reconstruct the energy of the deposition environment, particle size analysis was performed. The obtained results confirmed that river sediments are common at the bottom of Lake Balaton sediments, while biogenic carbonate dominates in the upper, Holocene part of the sediment core. To identify each mineral phase in the sediment, X-ray diffraction (XRD) studies were used to determine the ratio of calcite to Mg-calcite. Based on our XRF measurements, focusing primarily on quantitative changes in magnesium and calcium, transmission electron microscopy (TEM) studies were performed, mainly in the Mg enrichment layer around 8100 BP. The precipitation of biogenic carbonate in Lake Balaton is still taking place, mainly in the form of calcite and Mg-calcite. Their relative proportions strongly depend on the Mg saturation of the water and the substrates on which they are separated. From our results we can draw conclusions about the possible previous deeper phases of the lake and the evaporation conditions of the water. The data obtained from transmission electron microscopy shows a good agreement with the results of the XRF measurement, the proportion of Mg-calcite increases around 8100 BP that likely indicate drier climatic conditions connectable to the well know 8.2 ka cal BP climatic reversal.

How to cite: Pálfi, I., Pósfai, M., Kristály, F., Veres, D., Arnaud, F., Szalai, Z., Gresina, F., Szabó, Z., Darabos, G., Korponai, J., Czuppon, G., Demény, A., and Magyari, E.: Geochemical survey of Lake Balaton sediments: holocene paleoenvironment and paleoclimate, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-15654, https://doi.org/10.5194/egusphere-egu21-15654, 2021.

In Austral summer 2017/18 daily surface snow samples were taken (weather allowing) at two depths, 0-1cm and 6-7cm, at Neumayer III Station, Dronning Maud Land DML, Antarctica. Stable isotope ratios (18O, D, d-excess) of the snow samples were analysed in the AWI isotope lab. In parallel, water vapor stable isotopes were measured continuously on a routine base with a Picarro cavity ring-down spectroscope analyser (CRDS). Neumayer III is also a full meteorological observatory measuring all important meteorological variables including upper-air data. Meteorological data were directly compared to both snow and vapor isotope data. The corresponding synoptic situations were analysed using data from AMPS (Antarctic Mesoscale Prediction System), which employs WRF (Weather Research and Forecasting Model), a mesoscale atmospheric model that has been successfully used in earlier studies in DML. AMPS is run operationally at NCAR for Antarctic weather forecasting, particularly for flight operations of the US Antarctic Program (USAP). Additionally, back-trajectory calculations to investigate moisture sources and transport were carried out using FLEXPART, an open-source Lagrangian particle dispersion model. Due to logistical problems, the measuring period during the expedition was too short for statistical analysis, thus we focus on case studies here. In particular, periods with no precipitation were investigated, since earlier studies in Greenland have shown that the interaction of snow surface and atmosphere is important for the stable isotope ratio in the snow, thus in later ice cores that are used to derive paleo temperatures. A better understanding of the highly complex relationship between water vapor stable isotopes and meteorological conditions (including moisture source and transport) as well as the interaction between surface snow and water vapor is necessary for a correct paleoclimatic interpretation of ice cores.

How to cite: Schlosser, E., Bagheri, S., Powers, J. G., Manning, K. W., Hoerhold, M., Behrens, M., and Werner, M.: Atmospheric influences on water stable isotopes in Antarctic water vapor and surface snow – implications for ice core interpretation, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-2225, https://doi.org/10.5194/egusphere-egu21-2225, 2021.

Paleoclimate reconstructions from ice core records can be hampered due to the lack of a reliable chronology, especially when the stratigraphy is disturbed and conventional dating methods cannot be readily applied. The noble-gas radioisotopes ⁸¹Kr and ³⁹Ar can in these cases provide robust constraints as they yield absolute, radiometric ages. ⁸¹Kr (half-life 229 ka) covers the time span of 50-1300 ka, which is particularly relevant for polar ice cores, whereas ³⁹Ar (half-life 269 a) with a dating range of 50-1800 a is suitable for high mountain glaciers. For a long time the use of ⁸¹Kr and ³⁹Ar for dating of ice samples was hampered by the lack of a detection technique that can meet its extremely small abundance at a reasonable sample size.

Here, we present ⁸¹Kr and ³⁹Ar dating of Antarctic and Tibetan ice cores with the detection method Atom Trap Trace Analysis (ATTA), using 5-10 kg of ice for ⁸¹Kr and 2-5 kg for ³⁹Ar. Recent advances in further decreasing the sample size and increasing the dating precision will be discussed. Current studies include ⁸¹Kr dating in shallow ice cores from the Larsen Blue ice area, East Antarctica, in order to retrieve climate signals from the last glacial termination. Moreover, an ³⁹Ar profile from a central Tibetan ice core has been obtained in combination with layer counting based on isotopic and visual stratigraphic signals. The presented studies demonstrate how ⁸¹Kr and ³⁹Ar can constrain the age range of ice cores and complement other methods in developing an ice core chronology.

[1] Z.-T. Lu, Tracer applications of noble gas radionuclides in the geosciences, Earth-Science Reviews 138, 196-214, (2014)
[2] C. Buizert, Radiometric ⁸¹Kr dating identifies 120,000-year-old ice at Taylor Glacier, Antarctica, Proceedings of the National Academy of Sciences, 111, 6876, (2014)

[3] L. Tian, ⁸¹Kr Dating at the Guliya Ice Cap, Tibetan Plateau, Geophysical Research Letters, (2019)

http://atta.ustc.edu.cn

How to cite: Ritterbusch, F., Ahn, J., Gu, J.-Q., Jiang, W., Lee, G., Lu, Z.-T., Shao, L., Tian, L., Tong, A.-L., and Yang, G.-M.: Constraining ice core chronologies with 39Ar and 81Kr, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-3624, https://doi.org/10.5194/egusphere-egu21-3624, 2021.

Seasonal changes of meteorological and ground surface conditions can cause seasonal changes of source area of mineral particles supplied to the Greenland ice sheet. Difference in metal composition and size distribution of them reflects their source area difference. To clarify seasonal changes in the source area of mineral particles, we examined the metal composition and size distribution of them in snow pit samples obtained from EGRIP (East Greenland Ice Core Project).

In summer 2017, we dug a 2.01 m snow pit at EGRIP camp. Snow samples were collected at 0.03 m intervals. The snow samples were melted, then fractions of them were analyzed for the particle size distribution (0.52-12μm) with a Coulter counter (Beckman Coulter: Multisizer 4). Other fractions of the samples were treated with a microwave acid digestion method to decompose particulates. Total concentrations of Al, Ca and Na in samples were measured by inductively coupled plasma mass spectrometry.

Total concentrations of Al (t-Al) and concentrations of non-sea-salt Ca (nssCa) showed peak values in winter to spring layers. In those layers, nssCa/t-Al ratios and volume fractions of fine particles (<4 μm) tended to be relatively high. In some of those layers, the nssCa/t-Al ratios were 2-3 times higher than the crustal average of Ca/Al ratio. This result suggests that fine Ca-rich mineral particles originated from remote desert areas were supplied to EGRIP in those seasons. In contrast, in summer to autumn layers, the nssCa/t-Al ratios and volume fractions tended to be relatively low. This result can be explained by supplies of coarse (≥ 4 μm) Ca-poor mineral particles originated from soil areas near Greenland in those seasons.

How to cite: Komuro, Y., Nakazawa, F., Goto-Azuma, K., Hirabayashi, M., Shigeyama, W., and Nagatsuka, N.: Recent seasonal changes of metal composition and size distribution of mineral particles in snow at EGRIP, Greenland, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-4944, https://doi.org/10.5194/egusphere-egu21-4944, 2021.