The mid-Pliocene warm period (3.264–3.025 Ma) is the most recent geological period during which atmospheric CO2 levels were similar to recent historical values (∼400 ppm). Several proxy reconstructions for the mid-Pliocene show highly reduced zonal sea surface temperature (SST) gradients in the tropical Pacific Ocean, indicating an El Niño-like mean state. However, past modelling studies do not show these highly reduced gradients. Efforts to understand mid-Pliocene climate dynamics have led to the Pliocene Model Intercomparison Project (PlioMIP). Results from the first phase (PlioMIP1) showed clear El Niño variability (albeit significantly reduced) and did not show the greatly reduced time-mean zonal SST gradient suggested by some of the proxies.

In this work, we study El Niño–Southern Oscillation (ENSO) variability in the PlioMIP2 ensemble, which consists of additional global coupled climate models and updated boundary conditions compared to PlioMIP1. We quantify ENSO amplitude, period, spatial structure and “flavour”, as well as the tropical Pacific annual mean state in mid-Pliocene and pre-industrial simulations. Results show a reduced ENSO amplitude in the model-ensemble mean (−24 %) with respect to the pre-industrial, with 15 out of 17 individual models showing such a reduction. Furthermore, the spectral power of this variability considerably decreases in the 3–4-year band. The spatial structure of the dominant empirical orthogonal function shows no particular change in the patterns of tropical Pacific variability in the model-ensemble mean, compared to the pre-industrial. Although the time-mean zonal SST gradient in the equatorial Pacific decreases for 14 out of 17 models (0.2 ∘C reduction in the ensemble mean), there does not seem to be a correlation with the decrease in ENSO amplitude. The models showing the most “El Niño-like” mean state changes show a similar ENSO amplitude to that in the pre-industrial reference, while models showing more “La Niña-like” mean state changes generally show a large reduction in ENSO variability. The PlioMIP2 results show a reasonable agreement with both time-mean proxies indicating a reduced zonal SST gradient and reconstructions indicating a reduced, or similar, ENSO variability.

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 (similar to 3 Ma), a time when temperatures were 2-3 degrees 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.

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 degrees N) annual mean surface air temperature (SAT) increases of 3.7 to 11.6 degrees C compared to the pre-industrial period, with a multimodel mean (MMM) increase of 7.2 degrees 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 x 10(6) km(2), with a MMM anomaly of -5.6 x 10 6 km(2), which constitutes a decrease of 53 % compared to the pre-industrial period. The majority (11 out of 16) of models simulate summer seaice-free conditions (<= 1 x 10(6) km(2)) 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.

We investigate hydrology during a past climate slightly warmer than the present: the last interglacial (LIG). With daily output of preindustrial and LIG simulations from eight new climate models we force hydrological model PCR-GLOBWB and in turn hydrodynamic model CaMa-Flood. Compared to preindustrial, annual mean LIG runoff, discharge, and 100-yr 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 to the higher temperatures and evaporation. LIG discharge is much higher for the Niger, Congo, Nile, Ganges, Irrawaddy, and Pearl and lower for the Mississippi, Saint Lawrence, Amazon, Parana, Orange, Zambesi, Danube, and 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.

The Pliocene epoch has great potential to improve our understanding of the long-term climatic and environmental consequences of an atmospheric CO2 concentration near similar to 400 parts per million by volume. Here we present the large-scale features of Pliocene climate as simulated by a new ensemble of climate models of varying complexity and spatial resolution based on new reconstructions of boundary conditions (the Pliocene Model Intercomparison Project Phase 2; PlioMIP2). As a global annual average, modelled surface air temperatures increase by between 1.7 and 5.2 degrees C relative to the pre-industrial era with a multi-model mean value of 3.2 degrees C. Annual mean total precipitation rates increase by 7 % (range: 2 %-13 %). On average, surface air temperature (SAT) increases by 4.3 degrees C over land and 2.8 degrees C over the oceans. There is a clear pattern of polar amplification with warming polewards of 60 degrees N and 60 degrees S exceeding the global mean warming by a factor of 2.3. In the Atlantic and Pacific oceans, meridional temperature gradients are reduced, while tropical zonal gradients remain largely unchanged. There is a statistically significant relationship between a model's climate response associated with a doubling in CO2 (equilibrium climate sensitivity; ECS) and its simulated Pliocene surface temperature response. The mean ensemble Earth system response to a doubling of CO2 (including ice sheet feedbacks) is 67 % greater than ECS; this is larger than the increase of 47 % obtained from the PlioMIP1 ensemble. Proxy-derived estimates of Pliocene sea surface temperatures are used to assess model estimates of ECS and give an ECS range of 2.6-4.8 degrees C. This result is in general accord with the ECS range presented by previous Intergovernmental Panel on Climate Change (IPCC) Assessment Reports.