Description: The Pliocene Model Intercomparison Project is the first coordinated climate model comparison for a warmer palaeoclimate with atmospheric CO2 significantly higher than pre-industrial concentrations. The simulations of the mid-Pliocene warm period show global warming of between 1.8 and 3.6 °C above pre-industrial surface air temperatures, with significant polar amplification. Here we perform energy balance calculations on all eight of the coupled ocean–atmosphere simulations within PlioMIP Experiment 2 to evaluate the causes of the increased temperatures and differences between the models. In the tropics simulated warming is dominated by greenhouse gas increases, with cloud albedo feedbacks enhancing the warming in most of the models, but by widely varying amounts. The responses to mid-Pliocene climate forcing in the Northern Hemisphere mid-latitudes are substantially different between the climate models, with the only consistent response being a warming due to increased greenhouse gases. In the high latitudes all the energy balance components become important, but the dominant warming influence comes from the clear sky albedo. This demonstrates the importance of specified ice sheet and high latitude vegetation boundary conditions and simulated sea ice and snow albedo feedbacks. The largest components in the overall uncertainty are associated with cloud albedo feedbacks in the tropics and polar clear sky albedo, particularly in sea ice regions. These simulations show that high latitude albedo feedbacks provide the most significant enhancements to Pliocene greenhouse warming.

Description: Based on simulations with 15 climate models in the Pliocene Model Intercomparison Project (PlioMIP), the regional climate of East Asia (focusing on China) during the mid-Pliocene is investigated in this study. Compared to the pre-industrial, the multi-model ensemble mean (MMM) of all models shows the East Asian summer winds (EASWs) largely strengthen in monsoon China, and the East Asian winter winds (EAWWs) strengthen in south monsoon China but slightly weaken in north monsoon China in the mid-Pliocene. The MMM of all models also illustrates a warmer and wetter mid-Pliocene climate in China. The simulated weakened mid-Pliocene EAWWs in north monsoon China and intensified EASWs in monsoon China agree well with geological reconstructions. However, there is a large model–model discrepancy in simulating mid-Pliocene EAWW, which should be further addressed in the future work of PlioMIP.

Description: Thermodynamic arguments imply that global mean rainfall increases in a warmer atmosphere; however, dynamical effects may result in more significant diversity of regional precipitation change. Here we investigate rainfall changes in the mid-Pliocene Warm Period (~ 3 Ma), a time when temperatures were 2–3ºC warmer than the pre-industrial era, using output from the Pliocene Model Intercomparison Projects phases 1 and 2 and sensitivity climate model experiments. In the Mid-Pliocene simulations, the higher rates of warming in the northern hemisphere create an interhemispheric temperature gradient that enhances the southward cross-equatorial energy flux by up to 48%. This intensified energy flux reorganizes the atmospheric circulation leading to a northward shift of the Inter-Tropical Convergence Zone and a weakened and poleward displaced Southern Hemisphere Subtropical Convergences Zones. These changes result in drier-than-normal Southern Hemisphere tropics and subtropics. The evaluation of the mid-Pliocene adds a constraint to possible future warmer scenarios associated with differing rates of warming between hemispheres.

Description: The past provides evidence of abrupt climate shifts and changes in the frequency of climate and weather extremes. We explore the non‐linear response to orbital forcing and then consider climate millennial variability down to daily weather events. Orbital changes are translated into regional responses in temperature, where the precessional response is related to nonlinearities and seasonal biases in the system. We question regularities found in climate events by analyzing the distribution of inter‐event waiting times. Periodicities of about 900 and 1150 years are found in ice cores besides the prominent 1500‐years cycle. However, the variability remains indistinguishable from a random process, suggesting that centennial‐to‐millennial variability is stochastic in nature. New numerical techniques are developed allowing for a high resolution in the dynamically relevant regions like coasts, major upwelling regions, and high latitudes. Using this model, we find a strong sensitivity of the Atlantic meridional overturning circulation depending on where the deglacial meltwater is injected into. Meltwater into the Mississippi and near Labrador hardly affect the large‐scale ocean circulation, whereas subpolar hosing mimicking icebergs yields a quasi shutdown. The same multi‐scale approach is applied to radiocarbon simulations enabling a dynamical interpretation of marine sediment cores. Finally, abrupt climate events also have counterparts in the recent climate records, revealing a close link between climate variability, the statistics of North Atlantic weather patterns, and extreme events.

Description: A range of future climate scenarios are projected for high atmospheric CO2 concentrations, given uncertainties over future human actions as well as potential environmental and climatic feedbacks. The geological record offers an opportunity to understand climate system response to a range of forcings and feedbacks which operate over multiple temporal and spatial scales. Here, we examine a single interglacial during the late Pliocene (KM5c, ca. 3.205±0.01 Ma) when atmospheric CO2 exceeded pre-industrial concentrations, but were similar to today and to the lowest emission scenarios for this century. As orbital forcing and continental configurations were almost identical to today, we are able to focus on equilibrium climate system response to modern and near-future CO2. Using proxy data from 32 sites, we demonstrate that global mean sea-surface temperatures were warmer than pre-industrial values, by ∼2.3°C for the combined proxy data (foraminifera Mg∕Ca and alkenones), or by ∼3.2–3.4°C (alkenones only). Compared to the pre-industrial period, reduced meridional gradients and enhanced warming in the North Atlantic are consistently reconstructed. There is broad agreement between data and models at the global scale, with regional differences reflecting ocean circulation and/or proxy signals. An uneven distribution of proxy data in time and space does, however, add uncertainty to our anomaly calculations. The reconstructed global mean sea-surface temperature anomaly for KM5c is warmer than all but three of the PlioMIP2 model outputs, and the reconstructed North Atlantic data tend to align with the warmest KM5c model values. Our results demonstrate that even under low-CO2 emission scenarios, surface ocean warming may be expected to exceed model projections and will be accentuated in the higher latitudes.

Description: We present the Alfred Wegener Institute's contribution to the Pliocene Model Intercomparison Project Phase 2 (PlioMIP2) wherein we employ the Community Earth System Models (COSMOS) that include a dynamic vegetation scheme. This work builds on our contribution to Phase 1 of the Pliocene Model Intercomparison Project (PlioMIP1) wherein we employed the same model without dynamic vegetation. Our input to the PlioMIP2 special issue of Climate of the Past is twofold. In an accompanying paper we compare results derived with COSMOS in the framework of PlioMIP2 and PlioMIP1. With this paper we present details of our contribution with COSMOS to PlioMIP2. We provide a description of the model and of methods employed to transfer reconstructed mid-Pliocene geography, as provided by the Pliocene Reconstruction and Synoptic Mapping Initiative Phase 4 (PRISM4), to model boundary conditions. We describe the spin-up procedure for creating the COSMOS PlioMIP2 simulation ensemble and present large-scale climate patterns of the COSMOS PlioMIP2 mid-Pliocene core simulation. Furthermore, we quantify the contribution of individual components of PRISM4 boundary conditions to characteristics of simulated mid-Pliocene climate and discuss implications for anthropogenic warming. When exposed to PRISM4 boundary conditions, COSMOS provides insight into a mid-Pliocene climate that is characterised by increased rainfall (+0.17 mm d−1) and elevated surface temperature (+3.37 ∘C) in comparison to the pre-industrial (PI). About two-thirds of the mid-Pliocene core temperature anomaly can be directly attributed to carbon dioxide that is elevated with respect to PI. The contribution of topography and ice sheets to mid-Pliocene warmth is much smaller in contrast – about one-quarter and one-eighth, respectively, and nonlinearities are negligible. The simulated mid-Pliocene climate comprises pronounced polar amplification, a reduced meridional temperature gradient, a northwards-shifted tropical rain belt, an Arctic Ocean that is nearly free of sea ice during boreal summer, and muted seasonality at Northern Hemisphere high latitudes. Simulated mid-Pliocene precipitation patterns are defined by both carbon dioxide and PRISM4 paleogeography. Our COSMOS simulations confirm long-standing characteristics of the mid-Pliocene Earth system, among these increased meridional volume transport in the Atlantic Ocean, an extended and intensified equatorial warm pool, and pronounced poleward expansion of vegetation cover. By means of a comparison of our results to a reconstruction of the sea surface temperature (SST) of the mid-Pliocene we find that COSMOS reproduces reconstructed SST best if exposed to a carbon dioxide concentration of 400 ppmv. In the Atlantic to Arctic Ocean the simulated mid-Pliocene core climate state is too cold in comparison to the SST reconstruction. The discord can be mitigated to some extent by increasing carbon dioxide that causes increased mismatch between the model and reconstruction in other regions.

Description: We investigate hydrology during a past climate slightly warmer than the present: the Last Interglacial (LIG). With daily output of pre‐industrial and LIG simulations from eight new climate models we force hydrological model PCR‐GLOBWB, and in turn hydrodynamic model CaMa‐Flood. Compared to pre‐industrial, annual mean LIG runoff, discharge, and 100‐year flood volume are considerably larger in the Northern Hemisphere, by 14%, 25% and 82%, respectively. Anomalies are negative in the Southern Hemisphere. In some boreal regions, LIG runoff and discharge are lower despite higher precipitation, due the higher temperatures and evaporation. LIG discharge is much higher for the Niger, Congo, Nile, Ganges, Irrawaddy, Pearl, and lower for the Mississippi, Saint Lawrence, Amazon, Paraná, Orange, Zambesi, Danube, Ob. Discharge is seasonally postponed in tropical rivers affected by monsoon changes. Results agree with published proxies on the sign of discharge anomaly in 15 of 23 sites where comparison is possible.

Description: We compare results obtained from modeling the mid-Pliocene warm period using the Community Earth System Models (COSMOS, version: COSMOS-landveg r2413, 2009) with the two different modeling methodologies and sets of boundary conditions prescribed for the two phases of the Pliocene Model Intercomparison Project (PlioMIP), tagged PlioMIP1 and PlioMIP2. Here, we bridge the gap between our contributions to PlioMIP1 (Stepanek and Lohmann, 2012) and PlioMIP2 (Stepanek et al., 2020). We highlight some of the effects that differences in the chosen mid-Pliocene model setup (PlioMIP2 vs. PlioMIP1) have on the climate state as derived with COSMOS, as this information will be valuable in the framework of the model–model and model–data comparison within PlioMIP2. We evaluate the model sensitivity to improved mid-Pliocene boundary conditions using PlioMIP's core mid-Pliocene experiments for PlioMIP1 and PlioMIP2 and present further simulations in which we test model sensitivity to variations in paleogeography, orbit, and the concentration of CO2. Firstly, we highlight major changes in boundary conditions from PlioMIP1 to PlioMIP2 and also the challenges recorded from the initial effort. The results derived from our simulations show that COSMOS simulates a mid-Pliocene climate state that is 0.29°C colder in PlioMIP2 if compared to PlioMIP1 (17.82°C in PlioMIP1, 17.53°C in PlioMIP2; values based on simulated surface skin temperature). On the one hand, high-latitude warming, which is supported by proxy evidence of the mid-Pliocene, is underestimated in simulations of both PlioMIP1 and PlioMIP2. On the other hand, spatial variations in surface air temperature (SAT), sea surface temperature (SST), and the distribution of sea ice suggest improvement of simulated SAT and SST in PlioMIP2 if employing the updated paleogeography. Our PlioMIP2 mid-Pliocene simulation produces warmer SSTs in the Arctic and North Atlantic Ocean than those derived from the respective PlioMIP1 climate state. The difference in prescribed CO2 accounts for 0.5°C of temperature difference in the Arctic, leading to an ice-free summer in the PlioMIP1 simulation, and a quasi ice-free summer in PlioMIP2. Beyond the official set of PlioMIP2 simulations, we present further simulations and analyses that sample the phase space of potential alternative orbital forcings that have acted during the Pliocene and may have impacted geological records. Employing orbital forcing, which differs from that proposed for PlioMIP2 (i.e., corresponding to pre-industrial conditions) but falls into the mid-Pliocene time period targeted in PlioMIP, leads to pronounced annual and seasonal temperature variations. Our result identifies the changes in mid-Pliocene paleogeography from PRISM3 to PRISM4 as the major driver of the mid-Pliocene warmth within PlioMIP and not the minor differences in forcings.

Description: Palaeoclimate simulations improve our understanding of the climate, inform us about the performance of climate models in a different climate scenario, and help to identify robust features of the climate system. Here, we analyse Arctic warming in an ensemble of 16 simulations of the mid-Pliocene Warm Period (mPWP), derived from the Pliocene Model Intercomparison Project Phase 2 (PlioMIP2). The PlioMIP2 ensemble simulates Arctic (60–90 °N) annual mean surface air temperature (SAT) increases of 3.7 to 11.6 °C compared to the pre-industrial period, with a multi-model mean (MMM) increase of 7.2 °C. The Arctic warming amplification ratio relative to global SAT anomalies in the ensemble ranges from 1.8 to 3.1 (MMM is 2.3). Sea ice extent anomalies range from −3.0 to −10.4×106 km2, with a MMM anomaly of −5.6×106 km2, which constitutes a decrease of 53 % compared to the pre-industrial period. The majority (11 out of 16) of models simulate summer sea-ice-free conditions (≤1×106 km2) in their mPWP simulation. The ensemble tends to underestimate SAT in the Arctic when compared to available reconstructions, although the degree of underestimation varies strongly between the simulations. The simulations with the highest Arctic SAT anomalies tend to match the proxy dataset in its current form better. The ensemble shows some agreement with reconstructions of sea ice, particularly with regard to seasonal sea ice. Large uncertainties limit the confidence that can be placed in the findings and the compatibility of the different proxy datasets. We show that while reducing uncertainties in the reconstructions could decrease the SAT data–model discord substantially, further improvements are likely to be found in enhanced boundary conditions or model physics. Lastly, we compare the Arctic warming in the mPWP to projections of future Arctic warming and find that the PlioMIP2 ensemble simulates greater Arctic amplification than CMIP5 future climate simulations and an increase instead of a decrease in Atlantic Meridional Overturning Circulation (AMOC) strength compared to pre-industrial period. The results highlight the importance of slow feedbacks in equilibrium climate simulations, and that caution must be taken when using simulations of the mPWP as an analogue for future climate change.