ALBERT

Description: The modeling of paleoclimate, using physically based tools, is increasingly seen as a strong out-of-sample test of the models that are used for the projection of future climate changes. New to the Coupled Model Intercomparison Project (CMIP6) is the Tier 1 Last Interglacial experiment for 127 000 years ago (lig127k), designed to address the climate responses to stronger orbital forcing than the midHolocene experiment, using the same state-of-the-art models as for the future and following a common experimental protocol. Here we present a first analysis of a multi-model ensemble of 17 climate models, all of which have completed the CMIP6 DECK (Diagnostic, Evaluation and Characterization of Klima) experiments. The equilibrium climate sensitivity (ECS) of these models varies from 1.8 to 5.6 ∘C. The seasonal character of the insolation anomalies results in strong summer warming over the Northern Hemisphere continents in the lig127k ensemble as compared to the CMIP6 piControl and much-reduced minimum sea ice in the Arctic. The multi-model results indicate enhanced summer monsoonal precipitation in the Northern Hemisphere and reductions in the Southern Hemisphere. These responses are greater in the lig127k than the CMIP6 midHolocene simulations as expected from the larger insolation anomalies at 127 than 6 ka. New synthesis for surface temperature and precipitation, targeted for 127 ka, have been developed for comparison to the multi-model ensemble. The lig127k model ensemble and data reconstructions are in good agreement for summer temperature anomalies over Canada, Scandinavia, and the North Atlantic and for precipitation over the Northern Hemisphere continents. The model–data comparisons and mismatches point to further study of the sensitivity of the simulations to uncertainties in the boundary conditions and of the uncertainties and sparse coverage in current proxy reconstructions. The CMIP6–Paleoclimate Modeling Intercomparison Project (PMIP4) lig127k simulations, in combination with the proxy record, improve our confidence in future projections of monsoons, surface temperature, and Arctic sea ice, thus providing a key target for model evaluation and optimization.

Description: Based on the fifth phase of the Coupled Model Intercomparison Project (CMIP5)-generation previous Institut Pierre Simon Laplace (IPSL) Earth system model, we designed a new version, IPSL-CM5A2, aiming at running multi-millennial simulations typical of deep-time paleoclimate studies. Three priorities were followed during the setup of the model: (1) improving the overall model computing performance, (2) overcoming a persistent cold bias depicted in the previous model generation and (3) making the model able to handle the specific continental configurations of the geological past. These developments include the integration of hybrid parallelization Message Passing Interface – Open Multi-Processing (MPI-OpenMP) in the atmospheric model of the Laboratoire de Météorologie Dynamique (LMDZ), the use of a new library to perform parallel asynchronous input/output by using computing cores as “I/O servers” and the use of a parallel coupling library between the ocean and the atmospheric components. The model, which runs with an atmospheric resolution of 3.75∘×1.875∘ and 2 to 0.5∘ in the ocean, can now simulate ∼100 years per day, opening new possibilities towards the production of multi-millennial simulations with a full Earth system model. The tuning strategy employed to overcome a persistent cold bias is detailed. The confrontation of a historical simulation to climatological observations shows overall improved ocean meridional overturning circulation, marine productivity and latitudinal position of zonal wind patterns. We also present the numerous steps required to run IPSL-CM5A2 for deep-time paleoclimates through a preliminary case study for the Cretaceous. Namely, specific work on the ocean model grid was required to run the model for specific continental configurations in which continents are relocated according to past paleogeographic reconstructions. By briefly discussing the spin-up of such a simulation, we elaborate on the requirements and challenges awaiting paleoclimate modeling in the next years, namely finding the best trade-off between the level of description of the processes and the computing cost on supercomputers.

Description: Climate change has far-reaching implications in permafrost-underlain landscapes with respect to hydrology, ecosystems, and the population's traditional livelihoods. In the Lena River catchment, eastern Siberia, changing climatic conditions and the associated impacts are already observed or expected. However, as climate change progresses the question remains as to how far we are along this track and when these changes will constitute a significant emergence from natural variability. Here we present an approach to investigate temperature and precipitation time series from observational records, reanalysis, and an ensemble of 65 climate model simulations forced by the RCP8.5 emission scenario. We developed a novel non-parametric statistical method to identify the time of emergence (ToE) of climate change signals, i.e. the time when a climate signal permanently exceeds its natural variability. The method is based on the Hellinger distance metric that measures the similarity of probability density functions (PDFs) roughly corresponding to their geometrical overlap. Natural variability is estimated as a PDF for the earliest period common to all datasets used in the study (1901–1921) and is then compared to PDFs of target periods with moving windows of 21 years at annual and seasonal scales. The method yields dissimilarities or emergence levels ranging from 0 % to 100 % and the direction of change as a continuous time series itself. First, we showcase the method's advantage over the Kolmogorov–Smirnov metric using a synthetic dataset that resembles signals observed in the utilized climate models. Then, we focus on the Lena River catchment, where significant environmental changes are already apparent. On average, the emergence of temperature has a strong onset in the 1970s with a monotonic increase thereafter for validated reanalysis data. At the end of the reanalysis dataset (2004), temperature distributions have emerged by 50 %–60 %. Climate model projections suggest the same evolution on average and 90 % emergence by 2040. For precipitation the analysis is less conclusive because of high uncertainties in existing reanalysis datasets that also impede an evaluation of the climate models. Model projections suggest hardly any emergence by 2000 but a strong emergence thereafter, reaching 60 % by the end of the investigated period (2089). The presented ToE method provides more versatility than traditional parametric approaches and allows for a detailed temporal analysis of climate signal evolutions. An original strategy to select the most realistic model simulations based on the available observational data significantly reduces the uncertainties resulting from the spread in the 65 climate models used. The method comes as a toolbox available at https://github.com/pohleric/toe_tools (last access: 19 May 2020).

Description: The mid-Holocene (6000 years ago) is a standard time period for the evaluation of the simulated response of global climate models using palaeoclimate reconstructions. The latest mid-Holocene simulations are a palaeoclimate entry card for the Palaeoclimate Model Intercomparison Project (PMIP4) component of the current phase of the Coupled Model Intercomparison Project (CMIP6) – hereafter referred to as PMIP4-CMIP6. Here we provide an initial analysis and evaluation of the results of the experiment for the mid-Holocene. We show that state-of-the-art models produce climate changes that are broadly consistent with theory and observations, including increased summer warming of the Northern Hemisphere and associated shifts in tropical rainfall. Many features of the PMIP4-CMIP6 simulations were present in the previous generation (PMIP3-CMIP5) of simulations. The PMIP4-CMIP6 ensemble for the mid-Holocene has a global mean temperature change of −0.3 K, which is −0.2 K cooler than the PMIP3-CMIP5 simulations predominantly as a result of the prescription of realistic greenhouse gas concentrations in PMIP4-CMIP6. Biases in the magnitude and the sign of regional responses identified in PMIP3-CMIP5, such as the amplification of the northern African monsoon, precipitation changes over Europe, and simulated aridity in mid-Eurasia, are still present in the PMIP4-CMIP6 simulations. Despite these issues, PMIP4-CMIP6 and the mid-Holocene provide an opportunity both for quantitative evaluation and derivation of emergent constraints on the hydrological cycle, feedback strength, and potentially climate sensitivity.

Description: At present, variations of primary productivity (PP) in the Bay of Bengal (BoB) are driven by salinity-related stratification, which is controlled by the Indian summer monsoon (ISM). The relationships between PP, precipitation, and more generally climate in the past are not clearly understood. Here, we present a new record of PP based on the examination of coccolithophore assemblages in a 26 000-year sedimentary series retrieved in the northeastern BoB (core MD77-176). We compare our PP records to published climate and monsoon records, as well as outputs from numerical experiments obtained with the Earth system model IPSL-CM5A-LR, including the marine biogeochemical component PISCES, and with the transient climate simulation TraCE-21. Our results show that PP was most probably controlled by nutrient contents and distribution within the upper water column, which were predominantly influenced by (i) regional river systems between 26 and 19 ka, i.e. when sea level was relatively low and climate was relatively dry, and (ii) salinity-related stratification over the last 19 kyr, i.e. when sea level rose and more humid conditions prevailed. During that period, salinity and stratification were directly related to monsoon precipitation dynamics, which were chiefly forced by both insolation and Atlantic meridional overturning circulation (AMOC) strength. During Heinrich Stadial 1 and the Younger Dryas, i.e. when the AMOC collapsed, weaker South Asian precipitation diminished stratification and enhanced PP. During Bølling–Allerød, i.e. when the AMOC recovered, stronger South Asian precipitation increased stratification and subdued PP. Similarly, the precipitation peak recorded around the middle–early Holocene is consistent with a stronger stratification that drives PP minima.

Description: The current-climate Indian monsoon is known to boost biological productivity in the Arabian Sea. This paradigm has been extensively used to reconstruct past monsoon variability from palaeo-proxies indicative of changes in surface productivity. Here, we test this paradigm by simulating changes in marine primary productivity for eight contrasted climates from the last glacial–interglacial cycle. We show that there is no straightforward correlation between boreal summer productivity of the Arabian Sea and summer monsoon strength across the different simulated climates. Locally, productivity is fuelled by nutrient supply driven by Ekman dynamics. Upward transport of nutrients is modulated by a combination of alongshore wind stress intensity, which drives coastal upwelling, and by a positive wind stress curl to the west of the jet axis resulting in upward Ekman pumping. To the east of the jet axis there is however a strong downward Ekman pumping due to a negative wind stress curl. Consequently, changes in coastal alongshore stress and/or curl depend on both the jet intensity and position. The jet position is constrained by the Indian summer monsoon pattern, which in turn is influenced by the astronomical parameters and the ice sheet cover. The astronomical parameters are indeed shown to impact wind stress intensity in the Arabian Sea through large-scale changes in the meridional gradient of upper-tropospheric temperature. However, both the astronomical parameters and the ice sheets affect the pattern of wind stress curl through the position of the sea level depression barycentre over the monsoon region (20–150° W, 30° S–60° N). The combined changes in monsoon intensity and pattern lead to some higher glacial productivity during the summer season, in agreement with some palaeo-productivity reconstructions.

Description: The last deglaciation, which marked the transition between the last glacial and present interglacial periods, was punctuated by a series of rapid (centennial and decadal) climate changes. Numerical climate models are useful for investigating mechanisms that underpin the climate change events, especially now that some of the complex models can be run for multiple millennia. We have set up a Paleoclimate Modelling Intercomparison Project (PMIP) working group to coordinate efforts to run transient simulations of the last deglaciation, and to facilitate the dissemination of expertise between modellers and those engaged with reconstructing the climate of the last 21 000 years. Here, we present the design of a coordinated Core experiment over the period 21–9 thousand years before present (ka) with time-varying orbital forcing, greenhouse gases, ice sheets and other geographical changes. A choice of two ice sheet reconstructions is given, and we make recommendations for prescribing ice meltwater (or not) in the Core experiment. Additional focussed simulations will also be coordinated on an ad hoc basis by the working group, for example to investigate more thoroughly the effect of ice meltwater on climate system evolution, and to examine the uncertainty in other forcings. Some of these focussed simulations will target shorter durations around specific events in order to understand them in more detail and allow for the more computationally expensive models to take part.

Description: Climate change has far-reaching implications in permafrost-underlain landscapes with respect to hydrology, ecosystems and the population’s traditional livelihoods. In the Lena River catchment, Eastern Siberia, changing climatic conditions and the associated impacts are already observed or expected. However, as climate change progresses the question remains as to how far we are along this track and when these changes will constitute a significant emergence from natural variability. Here we present an approach to investigate temperature and precipitation time series from observational records, reanalysis, and an ensemble of 65 climate model simulations forced by the RCP8.5 emission scenario. We focus on the Lena River catchment, where significant environmental changes are already apparent. We developed a novel non-parametric statistical method to identify the time of emergence (ToE) of climate change signals, i.e. the time when a climate signal permanently exceeds its natural variability. The method is based on the Hellinger distance metric that measures the similarity of probability density functions (PDFs) roughly corresponding to their geometrical overlap. Natural variability is estimated as PDF for the earliest period common to all datasets used in the study (1901–1921) and is then compared to PDFs of target periods with moving windows of 21 years at annual and seasonal scale. The method yields dissimilarities or emergence levels ranging from 0 to 100 % and the direction of change as continuous time series itself. For the Lena River catchment, on average, emergence of temperature has a strong onset in the 1970s with a monotonic increase thereafter for validated reanalysis data. At the end of the reanalysis dataset (2004), temperature distributions have emerged by 50–60 %. Climate model projections suggest the same evolution on average and 90 % emergence by 2040. For precipitation the analysis is less conclusive because of high uncertainties in existing reanalysis datasets that also impede an evaluation of the climate models. Model projections suggest hardly any emergence by 2000 but a strong emergence thereafter, reaching 60 % by the end of the investigated period (2089). The presented ToE method provides more versatility than traditional parametric approaches and allows for a detailed temporal analysis of climate signal evolutions. An original strategy to select the most realistic model simulations based on the available observational data significantly reduces the uncertainties resulting from the spread in the 65 climate models used. The method comes as a toolbox available at https://github.com/pohleric/toe_tools.

Description: The penultimate deglaciation (PDG, ∼138–128 thousand years before present, hereafter ka) is the transition from the penultimate glacial maximum (PGM) to the Last Interglacial (LIG, ∼129–116 ka). The LIG stands out as one of the warmest interglacials of the last 800 000 years (hereafter kyr), with high-latitude temperature warmer than today and global sea level likely higher by at least 6 m. Considering the transient nature of the Earth system, the LIG climate and ice-sheet evolution were certainly influenced by the changes occurring during the penultimate deglaciation. It is thus important to investigate, with coupled atmosphere–ocean general circulation models (AOGCMs), the climate and environmental response to the large changes in boundary conditions (i.e. orbital configuration, atmospheric greenhouse gas concentrations, ice-sheet geometry and associated meltwater fluxes) occurring during the penultimate deglaciation. A deglaciation working group has recently been set up as part of the Paleoclimate Modelling Intercomparison Project (PMIP) phase 4, with a protocol to perform transient simulations of the last deglaciation (19–11 ka; although the protocol covers 26–0 ka). Similar to the last deglaciation, the disintegration of continental ice sheets during the penultimate deglaciation led to significant changes in the oceanic circulation during Heinrich Stadial 11 (∼136–129 ka). However, the two deglaciations bear significant differences in magnitude and temporal evolution of climate and environmental changes. Here, as part of the Past Global Changes (PAGES)-PMIP working group on Quaternary interglacials (QUIGS), we propose a protocol to perform transient simulations of the penultimate deglaciation under the auspices of PMIP4. This design includes time-varying changes in orbital forcing, greenhouse gas concentrations, continental ice sheets as well as freshwater input from the disintegration of continental ice sheets. This experiment is designed for AOGCMs to assess the coupled response of the climate system to all forcings. Additional sensitivity experiments are proposed to evaluate the response to each forcing. Finally, a selection of paleo-records representing different parts of the climate system is presented, providing an appropriate benchmark for upcoming model–data comparisons across the penultimate deglaciation.

Description: The Indian monsoon is known to boost biological productivity in the Arabian Sea. This paradigm has been extensively used to reconstruct past monsoon variability from paleo-proxies indicative of changes in surface productivity. Here, we test this paradigm by simulating changes in marine primary productivity for 8 contrasted climates from the last glacial-interglacial cycle. We show that there is no straightforward correlation between boreal summer productivity of the western and central Arabian Sea and summer monsoon strength across the different simulated climates. Locally, productivity is fueled by nutrient supply driven by Ekman dynamics depending on both wind stress and wind stress curl. In our simulations, a stronger monsoon with intensified wind stress on the Arabian Sea can lead to either increased or reduced productivity depending on the exact ocean surface on which increased wind stress and a positive wind stress curl are acting. The effects of winds onto the ocean are modulated by the Indian summer monsoon intensity and pattern (e.g. position of the low-level jet over the Arabian Sea), which in turn are driven by the orbital parameters and the ice sheet cover. The orbital parameters are indeed shown to impact wind stress intensity in the Arabian Sea through large scale changes in the meridional gradient of upper tropospheric temperature. But both the orbital parameters and the ice sheets affects the pattern of wind stress curl through the position of the sea level depression barycentre over the monsoon region (20° W–150° W, 30° S–60° N). The combined changes in monsoon intensity and pattern lead to higher glacial productivity during the summer season, in agreement with some paleo-productivity reconstructions.