High resolution climate and vegetation simulations of the Mid-Pliocene, a model-data comparison over western Europe and the Mediterranean region

International audienceHere we perform a detailed comparison between climate model results and climate reconstructions in western Europe and the Mediterranean area for the mid-Piacenzian warm interval (ca 3 Myr ago) of the Late Pliocene epoch. This region is particularly well suited for such a comparison as several quantitative climate estimates from local pollen records are available. They show evidence for temperatures significantly warmer than today over the whole area, mean annual precipitation higher in northwestern Europe and equivalent to modern values in its southwestern part. To improve our comparison, we have performed high resolution simulations of the mid-Piacenzian climate using the LMDz atmospheric general circulation model (AGCM) with a stretched grid which allows a finer resolution over Europe. In a first step, we applied the PRISM2 (Pliocene Research, Interpretation, and Synoptic Mapping) boundary conditions except that we used modern terrestrial vegetation. Second, we simulated the vegetation for this period by forcing the ORCHIDEE (Organizing Carbon and Hydrology in Dynamic Ecosystems) dynamic global vegetation model (DGVM) with the climatic outputs from the AGCM. We then supplied this simulated terrestrial vegetation cover as an additional boundary condition in a second AGCM run. This gives us the opportunity to investigate the model's sensitivity to the simulated vegetation changes in a global warming context

Ability of a dynamical climate sensitive disease model to reproduce historical Rift Valley Fever outbreaks over Africa

AbstractRift Valley Fever (RVF) is a zoonosis transmitted by Aedes and Culex mosquitoes, and is considered a priority pathogen by the WHO. RVF epidemics mostly occur in Africa and can decimate livestock herds, causing significant economic losses and posing health risks for humans. RVF transmission is associated with the occurrence of El Niño events that cause floods in eastern Africa and favour the emergence of mosquitoes in wetlands. Different risk models have been developed to forecast RVF transmission risk but very few studies have validated models at pan-African scale. This study aims to validate the skill of the Liverpool Rift Valley Fever model (LRVF) in reproducing RVF epidemics over Africa and to explore the relationship between simulated climatic suitability for RVF transmission and large-scale climate modes of variability such as the El Niño Southern Oscillation (ENSO) and the Dipole Mode Index (DMI). Our results show that the LRVF model correctly simulates RVF transmission hotspots and reproduces large epidemics that affected African countries. LRVF was able to correctly reproduce major RVF epidemics in Somalia, Kenya, Zambia and to a lesser extent for Mauritania and Senegal. The positive phases of ENSO and DMI are associated with an increased risk of RVF over the Horn of Africa, with important time lags. Following research activities should focus on the development of predictive modelling systems at different time scales.</jats:p

Heinrich event 1: an example of dynamical ice-sheet reaction to oceanic changes

Heinrich events, identified as enhanced ice-rafted detritus (IRD) in North Atlantic deep sea sediments (Heinrich, 1988; Hemming, 2004) have classically been attributed to Laurentide ice-sheet (LIS) instabilities (MacAyeal, 1993; Calov et al., 2002; Hulbe et al., 2004) and assumed to lead to important disruptions of the Atlantic meridional overturning circulation (AMOC) and North Atlantic deep water (NADW) formation. However, recent paleoclimate data have revealed that most of these events probably occurred after the AMOC had already slowed down or/and NADW largely collapsed, within about a thousand years (Hall et al., 2006; Hemming, 2004; Jonkers et al., 2010; Roche et al., 2004), implying that the initial AMOC reduction could not have been caused by the Heinrich events themselves. Here we propose an alternative driving mechanism, specifically for Heinrich event 1 (H1; 18 to 15 ka BP), by which North Atlantic ocean circulation changes are found to have strong impacts on LIS dynamics. By combining simulations with a coupled climate model and a three-dimensional ice sheet model, our study illustrates how reduced NADW and AMOC weakening lead to a subsurface warming in the Nordic and Labrador Seas resulting in rapid melting of the Hudson Strait and Labrador ice shelves. Lack of buttressing by the ice shelves implies a substantial ice-stream acceleration, enhanced ice-discharge and sea level rise, with peak values 500–1500 yr after the initial AMOC reduction. Our scenario modifies the previous paradigm of H1 by solving the paradox of its occurrence during a cold surface period, and highlights the importance of taking into account the effects of oceanic circulation on ice-sheets dynamics in order to elucidate the triggering mechanism of Heinrich events.Peer reviewe

Modeling the mid-piacenzian warm climate using the water isotope-enabled Community Earth System Model (iCESM1.2-ITPCAS)

The mid-Piacenzian Warm Period (MPWP, ~ 3.264–3.025 Ma) is the most recent example of a persistently warmer climate in equilibrium with atmospheric CO2 concentrations similar to today. Towards studying patterns and dynamics of a warming climate the MPWP is often compared to today. Following the Pliocene Model Intercomparison Project, Phase 2 (PlioMIP2) protocol we prepare a water isotope-enabled Community Earth System Model (iCESM1.2) simulation that is warmer and wetter than the PlioMIP2 multi-model ensemble (MME). While our simulation resembles PlioMIP2 MME in many aspects we find added insights. (1) Considerable warmth at high latitudes exceeds previous simulations. Polar amplification (PA) is comparable to proxies, enabled by iCESM1.2’s high climate sensitivity and a distinct method of ocean initialization. (2) Major driver of warmth is the downward component of clear-sky surface long-wave radiation. (3) In iCESM1.2 modulated dominance of dynamic (δDY) processes causes different low-latitude (~ 30 S°–10°N) precipitation response than the PlioMIP2 MME, where thermodynamic processes (δTH) dominate. (4) Modulated local condensation leads to lower δ18O across tropical Indian Ocean and surrounding Asian-African-Australian monsoon regions. (5) We find contrasting changes in tropical atmospheric circulations (Hadley and Walker cells). Anomalous regional meridional (zonal) circulation, forced by changes in tropical-subtropical (tropical) diabatic processes, presents a more comprehensive perspective than explaining weakened and expanded Hadley circulation (strengthened and westward-shifted Walker circulation) via static stability. (6) Enhanced Atlantic meridional overturning circulation owes to a closed Bering Strait

Large-scale features of Pliocene climate: results from the Pliocene Model Intercomparison Project

Climate and environments of the mid-Pliocene warm period (3.264 to 3.025 Ma) have been extensively studied. Whilst numerical models have shed light on the nature of climate at the time, uncertainties in their predictions have not been systematically examined. The Pliocene Model Intercomparison Project quantifies uncertainties in model outputs through a coordinated multi-model and multi-model/data intercomparison. Whilst commonalities in model outputs for the Pliocene are clearly evident, we show substantial variation in the sensitivity of models to the implementation of Pliocene boundary conditions. Models appear able to reproduce many regional changes in temperature reconstructed from geological proxies. However, data/model comparison highlights that models potentially underestimate polar amplification. To assert this conclusion with greater confidence, limitations in the time-averaged proxy data currently available must be addressed. Furthermore, sensitivity tests exploring the known unknowns in modelling Pliocene climate specifically relevant to the high latitudes are essential (e.g. palaeogeography, gateways, orbital forcing and trace gasses). Estimates of longer-term sensitivity to CO2 (also known as Earth System Sensitivity; ESS), support previous work suggesting that ESS is greater than Climate Sensitivity (CS), and suggest that the ratio of ESS to CS is between 1 and 2, with a "best" estimate of 1.5

Snowball Earth climate dynamics and Cryogenian geology-geobiology

Geological evidence indicates that grounded ice sheets reached sea level at all latitudes during two long-lived Cryogenian (58 and ≥5 My) glaciations. Combined uranium-lead and rhenium-osmium dating suggests that the older (Sturtian) glacial onset and both terminations were globally synchronous. Geochemical data imply that CO2 was 102 PAL (present atmospheric level) at the younger termination, consistent with a global ice cover. Sturtian glaciation followed breakup of a tropical supercontinent, and its onset coincided with the equatorial emplacement of a large igneous province. Modeling shows that the small thermal inertia of a globally frozen surface reverses the annual mean tropical atmospheric circulation, producing an equatorial desert and net snow and frost accumulation elsewhere. Oceanic ice thickens, forming a sea glacier that flows gravitationally toward the equator, sustained by the hydrologic cycle and by basal freezing and melting. Tropical ice sheets flow faster as CO2 rises but lose mass and become sensitive to orbital changes. Equatorial dust accumulation engenders supraglacial oligotrophic meltwater ecosystems, favorable for cyanobacteria and certain eukaryotes. Meltwater flushing through cracks enables organic burial and submarine deposition of airborne volcanic ash. The subglacial ocean is turbulent and well mixed, in response to geothermal heating and heat loss through the ice cover, increasing with latitude. Terminal carbonate deposits, unique to Cryogenian glaciations, are products of intense weathering and ocean stratification. Whole-ocean warming and collapsing peripheral bulges allow marine coastal flooding to continue long after ice-sheet disappearance. The evolutionary legacy of Snowball Earth is perceptible in fossils and living organisms

Orbitally forced ice sheet fluctuations during the Marinoan Snowball Earth glaciation

Two global glaciations occurred during the Neoproterozoic. Snowball Earth theory posits that these were terminated after millions of years of frigidity when initial warming from rising atmospheric CO2 concentrations was amplified by the reduction of ice cover and hence a reduction in planetary albedo. This scenario implies that most of the geological record of ice cover was deposited in a brief period of melt-back. However, deposits in low palaeo-latitudes show evidence of glacial–interglacial cycles. Here we analyse the sedimentology and oxygen and sulphur isotopic signatures of Marinoan Snowball glaciation deposits from Svalbard, in the Norwegian High Arctic. The deposits preserve a record of oscillations in glacier extent and hydrologic conditions under uniformly high atmospheric CO2 concentrations. We use simulations from a coupled three-dimensional ice sheet and atmospheric general circulation model to show that such oscillations can be explained by orbital forcing in the late stages of a Snowball glaciation. The simulations suggest that while atmospheric CO2 concentrations were rising, but not yet at the threshold required for complete melt-back, the ice sheets would have been sensitive to orbital forcing. We conclude that a similar dynamic can potentially explain the complex successions observed at other localities

Mid-Pliocene Atlantic Meridional Overturning Circulation simulated in PlioMIP2

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. However, 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 ∼ 8–12 ∘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

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×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