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.

The Mid-Piacenzian Warm Period (MPWP; 3.0–3.3 Ma), a warm geological period about three million years ago, has been deemed as a good past analog for understanding the current and future climate change. Based on 12 climate model outputs from Pliocene Model Intercomparison Project Phase 2 (PlioMIP2), we investigate tropical atmospheric circulation (TAC) changes under the warm MPWP and associated underlying mechanisms by diagnosing both atmospheric static stability and diabatic processes. Our findings underscore the advantage of analyzing atmospheric diabatic processes in elucidating seasonal variations of TAC compared to static stability assessments. Specifically, by diagnosing alterations in diabatic processes, we achieve a quantitative understanding and explanation the following TAC changes (incl. Strength and edge) during the MPWP: the weakened (annual, DJF, JJA) Northern Hemisphere and (DJF) Southern Hemisphere Hadley circulation (HC), reduced (annual, DJF) Pacific Walker circulation (PWC) and enhanced (annual, JJA) Southern Hemisphere HC and (JJA) PWC, and westward shifted (annual, DJF, JJA) PWC. We further addressed that the increasing bulk subtropical static stability and/or decreasing vertical shear of subtropical zonal wind - two crucial control factors for changes in subtropical baroclinicity - may promote HC widening, and vice versa. Consequently, our study of spatial diabatic heating and cooling, corresponding to upward and downward motions within the TAC, respectively, provides a new perspective for understanding the processes controlling seasonal TAC changes in response to surface warming.

The mid-Pliocene (∼3 Ma) is one of the most recent warm periods with high CO₂ concentrations in the atmosphere and resulting high temperatures, and it is often cited as an analog for near-term future climate change. Here, we apply a moisture budget analysis to investigate the response of the large-scale hydrological cycle at low latitudes within a 13-model ensemble from the Pliocene Model Intercomparison Project Phase 2 (PlioMIP2). The results show that increased atmospheric moisture content within the mid-Pliocene ensemble (due to the thermodynamic effect) results in wetter conditions over the deep tropics, i.e., the Pacific intertropical convergence zone (ITCZ) and the Maritime Continent, and drier conditions over the subtropics. Note that the dynamic effect plays a more important role than the thermodynamic effect in regional precipitation minus evaporation (PmE) changes (i.e., northward ITCZ shift and wetter northern Indian Ocean). The thermodynamic effect is offset to some extent by a dynamic effect involving a northward shift of the Hadley circulation that dries the deep tropics and moistens the subtropics in the Northern Hemisphere (i.e., the subtropical Pacific). From the perspective of Earth's energy budget, the enhanced southward cross-equatorial atmospheric transport (0.22 PW), induced by the hemispheric asymmetries of the atmospheric energy, favors an approximately 1^∘ northward shift of the ITCZ. The shift of the ITCZ reorganizes atmospheric circulation, favoring a northward shift of the Hadley circulation. In addition, the Walker circulation consistently shifts westward within PlioMIP2 models, leading to wetter conditions over the northern Indian Ocean. The PlioMIP2 ensemble highlights that an imbalance of interhemispheric atmospheric energy during the mid-Pliocene could have led to changes in the dynamic effect, offsetting the thermodynamic effect and, hence, altering mid-Pliocene hydroclimate.

In the Pliocene Model Intercomparison Project Phase 2 (PlioMIP2), coupled climate models have been used to simulate an interglacial climate during the mid-Piacenzian warm period (mPWP; 3.264 to 3.025 Ma). Here, we compare the Atlantic Meridional Overturning Circulation (AMOC), poleward ocean heat transport and sea surface warming in the Atlantic simulated with these models. In PlioMIP2, all models simulate an intensified mid-Pliocene AMOC. How- ever, there is no consistent response in the simulated Atlantic ocean heat transport nor in the depth of the Atlantic overturning cell. The models show a large spread in the simulated AMOC maximum, the Atlantic ocean heat transport and the surface warming in the North Atlantic. Although a few models simulate a surface warming of similar to 8-12 degrees C in the North Atlantic, similar to the reconstruction from Pliocene Research, Interpretation and Synoptic Mapping (PRISM) version 4, most models appear to underestimate this warming. The large model spread and model-data discrepancies in the PlioMIP2 ensemble do not support the hypothesis that an intensification of the AMOC, together with an increase in northward ocean heat transport, is the dominant mechanism for the mid-Pliocene warm climate over the North Atlantic.

The mid-Pliocene warm period (mPWP; ∼3.2 million years ago) is seen as the most recent time period characterized by a warm climate state, with similar to modern geography and ∼400 ppmv atmospheric CO₂ concentration, and is therefore often considered an interesting analogue for near-future climate projections. Paleoenvironmental reconstructions indicate higher surface temperatures, decreasing tropical deserts, and a more humid climate in West Africa characterized by a strengthened West African Monsoon (WAM). Using model results from the second phase of the Pliocene Modelling Intercomparison Project (PlioMIP2) ensemble, we analyse changes of the WAM rainfall during the mPWP by comparing them with the control simulations for the pre-industrial period. The ensemble shows a robust increase in the summer rainfall over West Africa and the Sahara region, with an average increase of 2.5 mm/d, contrasted by a rainfall decrease over the equatorial Atlantic. An anomalous warming of the Sahara and deepening of the Saharan Heat Low, seen in >90 % of the models, leads to a strengthening of the WAM and an increased monsoonal flow into the continent. A similar warming of the Sahara is seen in future projections using both phase 3 and 5 of the Coupled Model Intercomparison Project (CMIP3 and CMIP5). Though previous studies of future projections indicate a west–east drying–wetting contrast over the Sahel, PlioMIP2 simulations indicate a uniform rainfall increase in that region in warm climates characterized by increasing greenhouse gas forcing. We note that this effect will further depend on the long-term response of the vegetation to the CO₂ forcing.

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.

As global warming is proceeding due to rising greenhouse gas concentrations, the Earth system moves towards climate states that challenge adaptation. Past Earth system states are offering possible modelling systems for the global warming of the coming decades. These include the climate of the mid-Pliocene (similar to 3 Ma), the last interglacial (similar to 129-116 ka) and the mid-Holocene (similar to 6 ka). The simulations for these past warm periods are the key experiments in the Paleoclimate Model Intercomparison Project (PMIP) phase 4, contributing to phase 6 of the Coupled Model Intercomparison Project (CMIP6). Paleoclimate modelling has long been regarded as a robust out-of-sample test bed of the climate models used to project future climate changes. Here, we document the model setup for PMIP4 experiments with EC-Earth3-LR and present the large-scale features from the simulations for the mid-Holocene, the last interglacial and the mid-Pliocene. Using the pre-industrial climate as a reference state, we show global temperature changes, large-scale Hadley circulation and Walker circulation, polar warming, global monsoons and the climate variability modes - El Nino-Southern Oscillation (ENSO), the Pacific Decadal Oscillation (PDO) and the Atlantic Multidecadal Oscillation (AMO). EC-Earth3-LR simulates reasonable climate responses during past warm periods, as shown in the other PMIP4-CMIP6 model ensemble. The systematic comparison of these climate changes in past three warm periods in an individual model demonstrates the model's ability to capture the climate response under different climate forcings, providing potential implications for confidence in future projections with the EC-Earth model.

Future climate biogeochemical projections indicate large changes in the ocean with environmental conditions not experienced at present referred to as novel, or may even disappear. These climate-induced changes will most likely affect species distribution via changes in growth, behavior, evolution, dispersal, and species interactions. However, the future risk of novel and disappearing environmental conditions in the ocean is poorly understood, in particular for compound effects of climate and nutrient management changes. We map the compound risk of the occurrence of future novel and disappearing environmental conditions, analyze the outcome of climate and nutrient management scenarios for the world’s largest estuary, the Baltic Sea, and the potential consequences for three charismatic species. Overall, the future projections show, as expected, an increase in environmental novelty over time. The future nutrient reduction management that improves the eutrophication status of the Baltic Sea contributes to large novel and disappearing conditions. We show the consequences of novel and disappearing environmental conditions for fundamental niches of three charismatic species under different scenarios. This first step toward comprehensively analyzing environmental novelty and disappearing conditions for a marine system illustrates the urgent need to include novelty and disappearing projection outputs in Earth System Models. Our results further illustrate that adaptive management is needed to account for the emergence of novelty related to the interplay of multiple drivers. Overall, our analysis provides strong support for the expectation of novel ecological communities in marine systems, which may affect ecosystem services, and needs to be accounted for in sustainable future management plans of our oceans.

There is a well-known mode of rainfall variability associating opposite hydrological conditions over the Sahel region and the Gulf of Guinea, forming a dipole pattern. Previous meteorological observations show that the dipole pattern varies at interannual timescales. Using an EC-Earth climate model simulation for last millennium (850-1850 CE), we investigate the rainfall variability in West Africa over longer timescales. The 1000-year-long simulation data show that this rainfall dipole presents at decadal to multidecadal and centennial variability and long-term trend. Using the singular value decomposition (SVD) analysis, we identified that the rainfall dipole present in the first SVD mode with 60% explained variance and associated with the variabilities in tropical Atlantic sea surface temperature (SST). The second SVD mode shows a monopole rainfall variability pattern centred over the Sahel, associated with the extra-tropical Atlantic SST variability. We conclude that the rainfall dipole-like pattern is a natural variability mode originated from the local ocean-atmosphere-land coupling in the tropical Atlantic basin. The warm SST anomalies in the equatorial Atlantic Ocean favour an anomalous low pressure at the tropics. This low pressure weakens the meridional pressure gradient between the Saharan Heat Low and the tropical Atlantic. It leads to anomalous northeasterly, reduces the southwesterly moisture flux into the Sahel and confines the Gulf of Guinea's moisture convergence. The influence from extra-tropical climate variability, such as Atlantic multidecadal oscillation, tends to modify the rainfall dipole pattern to a monopole pattern from the Gulf of Guinea to Sahara through influencing the Sahara heat low. External forcing-such as orbital forcing, solar radiation, volcanic and land-use-can amplify/dampen the dipole mode through thermal forcing and atmosphere dynamical feedback.

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.7K (0.6-5.2, 5th-95th percentiles) using the PMIP2, PMIP3, and PMIP4 datasets for the LGM and 2.3K (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.7K (0.7-5.2) using the LGM and 2.3K (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.5K (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 5K and only exceeding 6K 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.