ALBERT

Description: The Last Glacial Maximum (LGM, 21,000 years ago) is one of the suite of paleoclimate simulations included in the current phase of the Coupled Model Intercomparison Project (CMIP6). It is an interval when insolation was similar to the present, but global ice volume was at a maximum, eustatic sea level was at or close to a minimum, greenhouse gas concentrations were lower, atmospheric aerosol loadings were higher than today, and vegetation and land-surface characteristics were different from today. The LGM has been a focus for the Paleoclimate Modelling Intercomparison Project (PMIP) since its inception, and thus many of the problems that might be associated with simulating such a radically different climate are well documented. The LGM state provides an ideal case study for evaluating climate model performance because the changes in forcing and temperature between the LGM and pre-industrial are of the same order of magnitude as those projected for the end of the 21st century. Thus, the CMIP6 LGM experiment could provide additional information that can be used to constrain estimates of climate sensitivity. The design of the Tier 1 LGM experiment (lgm) includes an assessment of uncertainties in boundary conditions, in particular through the use of different reconstructions of the ice sheets and of the change in dust forcing. Additional (Tier 2) sensitivity experiments have been designed to quantify feedbacks associated with land-surface changes and aerosol loadings, and to isolate the role of individual forcings. Model analysis and evaluation will capitalize on the relative abundance of paleoenvironmental observations and quantitative climate reconstructions already available for the LGM.

Description: This paper is the first of a series of four GMD (Geoscientific Model Development) papers on the PMIP4-CMIP6 (Paleoclimate Modelling Intercomparison Project - Phase 4 -- Coupled Model Intercomparison Project - Phase 6) experiments. Part 2 (Otto-Bliesner et al., 2017) gives details about the two PMIP4-CMIP6 interglacial experiments, Part 3 (Jungclaus et al., 2017) about the last millennium experiment, and Part 4 (Kageyama et al., 2017) about the Last Glacial Maximum experiment. The mid-Pliocene Warm Period experiment is part of the Pliocene Model Intercomparison Project (PlioMIP) - Phase 2, detailed in Haywood et al. (2016). The goal of the Paleoclimate Modelling Intercomparison Project (PMIP) is to understand the response of the climate system to different climate forcings for documented climatic states very different from the present and historical climates. Through comparison with observations of the environmental impact of these climate changes, or with climate reconstructions based on physical, chemical, or biological records, PMIP also addresses the issue of how well state-of-the-art numerical models simulate climate change. Climate models are usually developed using the present and historical climates as references, but climate projections show that future climates will lie well outside these conditions. Palaeoclimates very different from these reference states therefore provide stringent tests for state-of-the-art models and a way to assess whether their sensitivity to forcings is compatible with palaeoclimatic evidence. Simulations of five different periods have been designed to address the objectives of the sixth phase of the Coupled Model Intercomparison Project (CMIP6): the millennium prior to the industrial epoch (CMIP6 name: past1000); the mid-Holocene, 6000 years ago (midHolocene); the Last Glacial Maximum, 21,000 years ago (lgm); the Last Interglacial, 127,000 years ago (lig127k); and the mid-Pliocene Warm Period, 3.2 million years ago (midPliocene-eoi400). These climatic periods are well documented by palaeoclimatic and palaeoenvironmental records, with climate and environmental changes relevant for the study and projection of future climate changes. This paper describes the motivation for the choice of these periods and the design of the numerical experiments and database requests, with a focus on their novel features compared to the experiments performed in previous phases of PMIP and CMIP. It also outlines the analysis plan that takes advantage of the comparisons of the results across periods and across CMIP6 in collaboration with other MIPs.

Description: The pre-industrial millennium is among the periods selected by the Paleoclimate Model Intercomparison Project (PMIP) for experiments contributing to the sixth phase of the Coupled Model Intercomparison Project (CMIP6) and the fourth phase of the PMIP (PMIP4). The past1000 transient simulations serve to investigate the response to (mainly) natural forcing under background conditions not too different from today, and to discriminate between forced and internally generated variability on interannual to centennial timescales. This paper describes the motivation and the experimental set-ups for the PMIP4-CMIP6 past1000 simulations, and discusses the forcing agents orbital, solar, volcanic, and land use/land cover changes, and variations in greenhouse gas concentrations. The past1000 simulations covering the pre-industrial millennium from 850 Common Era (CE) to 1849 CE have to be complemented by historical simulations (1850 to 2014 CE) following the CMIP6 protocol. The external forcings for the past1000 experiments have been adapted to provide a seamless transition across these time periods. Protocols for the past1000 simulations have been divided into three tiers. A default forcing data set has been defined for the Tier 1 (the CMIP6 past1000) experiment. However, the PMIP community has maintained the flexibility to conduct coordinated sensitivity experiments to explore uncertainty in forcing reconstructions as well as parameter uncertainty in dedicated Tier 2 simulations. Additional experiments (Tier 3) are defined to foster collaborative model experiments focusing on the early instrumental period and to extend the temporal range and the scope of the simulations. This paper outlines current and future research foci and common analyses for collaborative work between the PMIP and the observational communities (reconstructions, instrumental data).

Description: The last extended time period when climate may have been warmer than today was during the Last Interglacial (LIG; ca. 129 to 120 thousand years ago). However, a global view of LIG precipitation is lacking. Here, seven new LIG climate models are compared to the first global database of proxies for LIG precipitation. In this way, models are assessed in their ability to capture important hydroclimatic processes during a different climate. The models can reproduce the proxy-based positive precipitation anomalies from the preindustrial period over much of the boreal continents. Over the Southern Hemisphere, proxy-model agreement is partial. In models, LIG boreal monsoons have 42% wider area than in the preindustrial and produce 55% more precipitation and 50% more extreme precipitation. Austral monsoons are weaker. The mechanisms behind these changes are consistent with stronger summer radiative forcing over boreal high latitudes and with the associated higher temperatures during the LIG.

Description: The Last Glacial Maximum (LGM, 21 000 years ago) is one of the suite of paleoclimate simulations included in the current phase of the Coupled Model Intercomparison Project (CMIP6). It is an interval when insolation was similar to the present, but global ice volume was at a maximum, eustatic sea level was at or close to a minimum, greenhouse gas concentrations were lower, atmospheric aerosol loadings were higher than today, and vegetation and land-surface characteristics were different from today. The LGM has been a focus for the Paleoclimate Modelling Intercomparison Project (PMIP) since its inception, and thus many of the problems that might be associated with simulating such a radically different climate are well documented. The LGM state provides an ideal case study for evaluating climate model performance because the changes in forcing and temperature between the LGM and pre-industrial are of the same order of magnitude as those projected for the end of the 21st century. Thus, the CMIP6 LGM experiment could provide additional information that can be used to constrain estimates of climate sensitivity. The design of the Tier 1 LGM experiment (lgm) includes an assessment of uncertainties in boundary conditions, in particular through the use of different reconstructions of the ice sheets and of the change in dust forcing. Additional (Tier 2) sensitivity experiments have been designed to quantify feedbacks associated with land-surface changes and aerosol loadings, and to isolate the role of individual forcings. Model analysis and evaluation will capitalize on the relative abundance of paleoenvironmental observations and quantitative climate reconstructions already available for the LGM.

Description: Two interglacial epochs are included in the suite of Paleoclimate Modeling Intercomparison Project (PMIP4) simulations in the Coupled Model Intercomparison Project (CMIP6). The experimental protocols for simulations of the mid-Holocene (midHolocene, 6000 years before present) and the Last Interglacial (lig127k, 127 000 years before present) are described here. These equilibrium simulations are designed to examine the impact of changes in orbital forcing at times when atmospheric greenhouse gas levels were similar to those of the preindustrial period and the continental configurations were almost identical to modern ones. These simulations test our understanding of the interplay between radiative forcing and atmospheric circulation, and the connections among large-scale and regional climate changes giving rise to phenomena such as land–sea contrast and high-latitude amplification in temperature changes, and responses of the monsoons, as compared to today. They also provide an opportunity, through carefully designed additional sensitivity experiments, to quantify the strength of atmosphere, ocean, cryosphere, and land-surface feedbacks. Sensitivity experiments are proposed to investigate the role of freshwater forcing in triggering abrupt climate changes within interglacial epochs. These feedback experiments naturally lead to a focus on climate evolution during interglacial periods, which will be examined through transient experiments. Analyses of the sensitivity simulations will also focus on interactions between extratropical and tropical circulation, and the relationship between changes in mean climate state and climate variability on annual to multi-decadal timescales. The comparative abundance of paleoenvironmental data and of quantitative climate reconstructions for the Holocene and Last Interglacial make these two epochs ideal candidates for systematic evaluation of model performance, and such comparisons will shed new light on the importance of external feedbacks (e.g., vegetation, dust) and the ability of state-of-the-art models to simulate climate changes realistically.

Description: In this paper we introduce a Bayesian framework, which is explicit about prior assumptions, for using model ensembles and observations together to constrain future climate change. The emergent constraint approach has seen broad application in recent years, including studies constraining the equilibrium climate sensitivity (ECS) using the Last Glacial Maximum (LGM) and the mid-Pliocene Warm Period (mPWP). Most of these studies were based on ordinary least squares (OLS) fits between a variable of the climate state, such as tropical temperature, and climate sensitivity. Using our Bayesian method, and considering the LGM and mPWP separately, we obtain values of ECS of 2.7 K (0.6–5.2, 5th–95th percentiles) using the PMIP2, PMIP3, and PMIP4 datasets for the LGM and 2.3 K (0.5–4.4) with the PlioMIP1 and PlioMIP2 datasets for the mPWP. Restricting the ensembles to include only the most recent version of each model, we obtain 2.7 K (0.7–5.2) using the LGM and 2.3 K (0.4–4.5) using the mPWP. An advantage of the Bayesian framework is that it is possible to combine the two periods assuming they are independent, whereby we obtain a tighter constraint of 2.5 K (0.8–4.0) using the restricted ensemble. We have explored the sensitivity to our assumptions in the method, including considering structural uncertainty, and in the choice of models, and this leads to 95 % probability of climate sensitivity mostly below 5 K and only exceeding 6 K in a single and most uncertain case assuming a large structural uncertainty. The approach is compared with other approaches based on OLS, a Kalman filter method, and an alternative Bayesian method. An interesting implication of this work is that OLS-based emergent constraints on ECS generate tighter uncertainty estimates, in particular at the lower end, an artefact due to a flatter regression line in the case of lack of correlation. Although some fundamental challenges related to the use of emergent constraints remain, this paper provides a step towards a better foundation for their potential use in future probabilistic estimations of climate sensitivity.

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

Description: Surface ocean pH is declining due to anthropogenic atmospheric CO2 uptake with a global decline of ~0.3 possible by 2100. Extracellular pH influences a range of biological processes, including nutrient uptake, calcification and silicification. However, there are poor constraints on how pH levels in the extracellular microenvironment surrounding phytoplankton cells (the phycosphere) differ from bulk seawater. This adds uncertainty to biological impacts of environmental change. Furthermore, previous modelling work suggests that phycosphere pH of small cells is close to bulk seawater, and this has not been experimentally verified. Here we observe under 140 μmol photons·m−2·s−1 the phycosphere pH of Chlamydomonas concordia (5 µm diameter), Emiliania huxleyi (5 µm), Coscinodiscus radiatus (50 µm) and C. wailesii (100 µm) are 0.11 ± 0.07, 0.20 ± 0.09, 0.41 ± 0.04 and 0.15 ± 0.20 (mean ± SD) higher than bulk seawater (pH 8.00), respectively. Thickness of the pH boundary layer of C. wailesii increases from 18 ± 4 to 122 ± 17 µm when bulk seawater pH decreases from 8.00 to 7.78. Phycosphere pH is regulated by photosynthesis and extracellular enzymatic transformation of bicarbonate, as well as being influenced by light intensity and seawater pH and buffering capacity. The pH change alters Fe speciation in the phycosphere, and hence Fe availability to phytoplankton is likely better predicted by the phycosphere, rather than bulk seawater. Overall, the precise quantification of chemical conditions in the phycosphere is crucial for assessing the sensitivity of marine phytoplankton to ongoing ocean acidification and Fe limitation in surface oceans.

Description: Highlights: • Multi-centennial oscillation with 100–200 years periods is evident in proxy data and model simulations during the Holocene. • Multi-centennial oscillation is a global signal and is more significant in the Northern Hemisphere high latitudes. • None of the external forcings is found to be the sole driver of the multi-centennial variability. • It indicates the multi-centennial oscillation may be due to potential internal drivers and essential feedbacks. Abstract: Variability on centennial to multi-centennial timescales is mentioned as a feature in reconstructions of the Holocene climate. As more long transient model simulations with complex climate models become available and efforts have been made to compile large proxy databases, there is now a unique opportunity to study multi-centennial variability with greater detail and a large amount of data than earlier. This paper presents a spectral analysis of transient Holocene simulations from 9 models and 120 proxy records to find the common signals related to oscillation periods and geographic dependencies and discuss the implications for the potential driving mechanisms. Multi-centennial variability is significant in most proxy records, with the dominant oscillation periods around 120–130 years and an average of 240 years. Spectra of model-based global mean temperature (GMT) agree well with proxy evidence with significant multi-centennial variability in all simulations with the dominant oscillation periods around 120–150 years. It indicates a comparatively good agreement between model and proxy data. A lack of latitudinal dependencies in terms of oscillation period is found in both the model and proxy data. However, all model simulations have the highest spectral density distributed over the Northern hemisphere high latitudes, which could indicate a particular variability sensitivity or potential driving mechanisms in this region. Five models also have differentiated forcings simulations with various combinations of forcing agents. Significant multi-centennial variability with oscillation periods between 100 and 200 years is found in all forcing scenarios, including those with only orbital forcing. The different forcings induce some variability in the system. Yet, none appear to be the predominant driver based on the spectral analysis. Solar irradiance has long been hypothesized to be a primary driver of multi-centennial variability. However, all the simulations without this forcing have shown significant multi-centennial variability. The results then indicate that internal mechanisms operate on multi-centennial timescales, and the North Atlantic-Arctic is a region of interest for this aspect.