Modelling the mid-Pliocene Warm Period climate with the IPSL coupled model and its atmospheric component LMDZ5A

This paper describes the experimental design and model results of the climate simulations of the mid-Pliocene Warm Period (mPWP, ca. 3.3–3 Ma) using the Institut Pierre Simon Laplace model (IPSLCM5A), in the framework of the Pliocene Model Intercomparison Project (PlioMIP). We use the IPSL atmosphere ocean general circulation model (AOGCM), and its atmospheric component alone (AGCM), to simulate the climate of the mPWP. Boundary conditions such as sea surface temperatures (SSTs), topography, ice-sheet extent and vegetation are derived from the ones imposed by the Pliocene Model Intercomparison Project (PlioMIP), described in Haywood et al. (2010, 2011). We first describe the IPSL model main features, and then give a full description of the boundary conditions used for atmospheric model and coupled model experiments. The climatic outputs of the mPWP simulations are detailed and compared to the corresponding control simulations. The simulated warming relative to the control simulation is 1.94 °C in the atmospheric and 2.07 °C in the coupled model experiments. In both experiments, warming is larger at high latitudes. Mechanisms governing the simulated precipitation patterns are different in the coupled model than in the atmospheric model alone, because of the reduced gradients in imposed SSTs, which impacts the Hadley and Walker circulations. In addition, a sensitivity test to the change of land-sea mask in the atmospheric model, representing a sea-level change from present-day to 25 m higher during the mid-Pliocene, is described. We find that surface temperature differences can be large (several degrees Celsius) but are restricted to the areas that were changed from ocean to land or vice versa. In terms of precipitation, impact on polar regions is minor although the change in land-sea mask is significant in these areas

Sea Surface Temperature of the mid-Piacenzian Ocean:A Data-Model Comparison

The mid-Piacenzian climate represents the most geologically recent interval of long-term average warmth relative to the last million years, and shares similarities with the climate projected for the end of the 21st century. As such, it represents a natural experiment from which we can gain insight into potential climate change impacts, enabling more informed policy decisions for mitigation and adaptation. Here, we present the first systematic comparison of Pliocene sea surface temperature (SST) between an ensemble of eight climate model simulations produced as part of PlioMIP (Pliocene Model Intercomparison Project) with the PRISM (Pliocene Research, Interpretation and Synoptic Mapping) Project mean annual SST field. Our results highlight key regional and dynamic situations where there is discord between the palaeoenvironmental reconstruction and the climate model simulations. These differences have led to improved strategies for both experimental design and temporal refinement of the palaeoenvironmental reconstruction

Arctic sea ice simulation in the PlioMIP ensemble

Eight general circulation models have simulated the mid-Pliocene warm period (mid-Pliocene, 3.264 to 3.025 Ma) as part of the Pliocene Modelling Intercomparison Project (PlioMIP). Here, we analyse and compare their simulation of Arctic sea ice for both the pre-industrial period and the mid-Pliocene. Mid-Pliocene sea ice thickness and extent is reduced, and the model spread of extent is more than twice the pre-industrial spread in some summer months. Half of the PlioMIP models simulate ice-free conditions in the mid-Pliocene. This spread amongst the ensemble is in line with the uncertainties amongst proxy reconstructions for mid-Pliocene sea ice extent. Correlations between mid-Pliocene Arctic temperatures and sea ice extents are almost twice as strong as the equivalent correlations for the pre-industrial simulations. The need for more comprehensive sea ice proxy data is highlighted, in order to better compare model performances

Evaluation of Arctic warming in mid-Pliocene climate simulations

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×10^{6} km^{2}, with a MMM anomaly of −5.6×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 sea-ice-free conditions (≤1×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

A return to large-scale features of Pliocene climate: the Pliocene Model Intercomparison Project Phase 2

The Pliocene epoch has great potential to improve our understanding of the long-term climatic and environmental consequences of an atmospheric CO2 concentration near ~ 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 and 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.4 and 4.7 °C relative to pre-industrial with a multi-model mean value of 2.8 °C. Annual mean total precipitation rates increase by 6 % (range: 2 %–13 %). On average, surface air temperature (SAT) increases are 1.3 °C greater over the land than over the oceans, and there is a clear pattern of polar amplification with warming polewards of 60° N and 60° S exceeding the global mean warming by a factor of 2.4. In the Atlantic and Pacific Oceans, meridional temperature gradients are reduced, while tropical zonal gradients remain largely unchanged. Although there are some modelling constraints, 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 doubling of CO2 (including ice sheet feedbacks) is approximately 50 % greater than ECS, consistent with results from the PlioMIP1 ensemble. Proxy-derived estimates of Pliocene sea-surface temperatures are used to assess model estimates of ECS and indicate a range in ECS from 2.5 to 4.3 °C. This result is in general accord with the range in ECS presented by previous IPCC Assessment Reports

The Pliocene Model Intercomparison Project Phase 2: large-scale climate features and climate sensitivity

The Pliocene epoch has great potential to improve our understanding of the long-term climatic and environmental consequences of an atmospheric CO2 concentration near ∼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 ∘C relative to the pre-industrial era with a multi-model mean value of 3.2 ∘C. Annual mean total precipitation rates increase by 7 % (range: 2 %–13 %). On average, surface air temperature (SAT) increases by 4.3 ∘C over land and 2.8 ∘C over the oceans. There is a clear pattern of polar amplification with warming polewards of 60∘ N and 60∘ 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 ∘C. This result is in general accord with the ECS range presented by previous Intergovernmental Panel on Climate Change (IPCC) Assessment Reports

Unraveling the mechanisms and implications of a stronger mid-Pliocene Atlantic Meridional Overturning Circulation (AMOC) in PlioMIP2

The mid-Pliocene warm period (3.264–3.025 Ma) is the most recent geological period in which the atmospheric CO2 concentration was approximately equal to the concentration we measure today (ca. 400 ppm). Sea surface temperature (SST) proxies indicate above-average warming over the North Atlantic in the mid-Pliocene with respect to the pre-industrial period, which may be linked to an intensified Atlantic Meridional Overturning Circulation (AMOC). Earlier results from the Pliocene Model Intercomparison Project Phase 2 (PlioMIP2) show that the ensemble simulates a stronger AMOC in the mid-Pliocene than in the pre-industrial. However, no consistent relationship between the stronger mid-Pliocene AMOC and either the Atlantic northward ocean heat transport (OHT) or average North Atlantic SSTs has been found. In this study, we look further into the drivers and consequences of a stronger AMOC in mid-Pliocene compared to pre-industrial simulations in PlioMIP2. We find that all model simulations with a closed Bering Strait and Canadian Archipelago show reduced freshwater transport from the Arctic Ocean into the North Atlantic. This contributes to an increase in salinity in the subpolar North Atlantic and Labrador Sea that can be linked to the stronger AMOC in the mid-Pliocene. To investigate the dynamics behind the ensemble's variable response of the total Atlantic OHT to the stronger AMOC, we separate the Atlantic OHT into two components associated with either the overturning circulation or the wind-driven gyre circulation. While the ensemble mean of the overturning component is increased significantly in magnitude in the mid-Pliocene, it is partly compensated by a reduction in the gyre component in the northern subtropical gyre region. This indicates that the lack of relationship between the total OHT and AMOC is due to changes in OHT by the subtropical gyre. The overturning and gyre components should therefore be considered separately to gain a more complete understanding of the OHT response to a stronger mid-Pliocene AMOC. In addition, we show that the AMOC exerts a stronger influence on North Atlantic SSTs in the mid-Pliocene than in the pre-industrial, providing a possible explanation for the improved agreement of the PlioMIP2 ensemble mean SSTs with reconstructions in the North Atlantic

Reduced El Niño variability in the mid-Pliocene according to the PlioMIP2 ensemble

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

Decomposition of physical processes controlling EASM precipitation changes during the mid-Piacenzian: new insights into data–model integration

The mid-Piacenzian warm period (MPWP, ~3.264–3.025 Ma) has gained widespread interest due to its partial analogy with future climate. However, quantitative data–model comparison of East Asian Summer Monsoon (EASM) precipitation during the MPWP is relatively rare, especially due to problems in decoding the imprint of physical processes to climate signals in the records. In this study, pollen-based precipitation records are reconstructed and compared to the multi-model ensemble mean of the Pliocene Model Intercomparison Project Phase 2 (PlioMIP2). We find spatially consistent precipitation increase in most simulations but a spatially divergent change in MPWP records. We reconcile proxy data and simulation by decomposing physical processes that control precipitation. Our results 1) reveal thermodynamic control of an overall enhancement of EASM precipitation and 2) highlight a distinct control of thermodynamic and dynamical processes on increases of tropical and subtropical EASM precipitation, reflecting the two pathways of water vapor supply that enhance EASM precipitation, respectively