A latest Cretaceous to earliest Paleogene dinoflagellate cyst zonation of Antarctica, and implications for phytoprovincialism in the high southern latitudes

The thickest uppermost Cretaceous to lowermost Paleogene (Maastrichtian to Danian) sedimentary succession in the world is exposed on southern Seymour Island (65° South) in the James Ross Basin, Antarctic Peninsula. This fossiliferous shallow marine sequence, which spans the Cretaceous–Paleogene boundary, has allowed a high-resolution analysis of well-preserved marine palynomorphs. Previous correlation of Cretaceous–Paleogene marine palynomorph assemblages in the south polar region relied on dinoflagellate cyst biozonations from New Zealand and southern Australia. The age model of the southern Seymour Island succession is refined and placed within the stratigraphical context of the mid to high southern palaeolatitudes. Quantitative palynological analysis of a new 1102 m continuous stratigraphical section comprising the uppermost Snow Hill Island Formation and the López de Bertodano Formation (Marambio Group) across southern Seymour Island was undertaken. We propose the first formal late Maastrichtian to early Danian dinoflagellate cyst zonation scheme for the Antarctic based on this exceptional succession. Two new late Maastrichtian zones, including three subzones, and one new early Danian zone are defined. The oldest beds correlate well with the late Maastrichtian of New Zealand. In a wider context, a new South Polar Province based on Maastrichtian to Danian dinoflagellate cysts is proposed, which excludes most southern South American marine palynofloras. This interpretation is supported by models of ocean currents around Antarctica and implies an unrestricted oceanic connection across Antarctica between southern South America and the Tasman Sea

The transient response of ice volume to orbital forcing during the warm Late Pliocene

Examining the nature of ice sheet and sea level response to past episodes of enhanced greenhouse gas forcing may help constrain future sea level change. Here, for the first time, we present the transient nature of ice sheets and sea level during the late Pliocene. The transient ice sheet predictions are forced by multiple climate snapshots derived from a climate model set up with late Pliocene boundary conditions, forced with different orbital forcing scenarios appropriate to two Marine Isotope Stages (MISs), MIS KM5c, and K1. Our results indicate that during MIS KM5c both the Antarctic and Greenland ice sheets contributed to sea level rise relative to present and were relatively stable. Insolation forcing between the hemispheres was out of phase during MIS K1 and led to an asynchronous response of ice volume globally. Therefore, when variations of precession were high, inferring the behavior of ice sheets from benthic isotope or sea level records is complex

Modelling the enigmatic Late Pliocene Glacial Event: Marine Isotope Stage M2

The Pliocene Epoch (5.2 to 2.58Ma) has often been targeted to investigate the nature ofwarmclimates. However, climate records for the Pliocene exhibit significant variability and show intervals that apparently experienced a cooler than modern climate. Marine Isotope Stage (MIS) M2 (~3.3 Ma) is a globally recognisable cooling event that disturbs an otherwise relatively (compared to present-day) warm background climate state. It remains unclear whether this event corresponds to significant ice sheet build-up in the Northern and Southern Hemisphere. Estimates of sea level for this interval vary, and range from modern values to estimates of 65 m sea level fall with respect to present day. Here we implement plausibleM2 ice sheet configurations into a coupled atmosphere–ocean climate model to test the hypothesis that larger-than-modern ice sheet configurations may have existed at M2. Climate model results are compared with proxy climate data available for M2 to assess the plausibility of each ice sheet configuration. Whilst the outcomes of our data/model comparisons are not in all cases straight forward to interpret, there is little indication that results from model simulations in which significant ice masses have been prescribed in the Northern Hemisphere are incompatible with proxy data from the North Atlantic, Northeast Arctic Russia, North Africa and the Southern Ocean. Therefore, our model results do not preclude thepossibilityof the existenceof larger icemasses duringM2 in the Northern or SouthernHemisphere. Specifically they are not able to discount the possibility of significant icemasses in the Northern Hemisphere during the M2 event, consistent with a global sea-level fall of between 40 m and 60 m. This study highlights the general need for more focused and coordinated data generation in the future to improve the coverage and consistency in proxy records for M2, which will allow these and future M2 sensitivity tests to be interrogated further

The warm winter paradox in the Pliocene northern high latitudes

Reconciling palaeodata with model simulations of the Pliocene climate is essential for understanding a world with atmospheric CO2 concentration near 400 ppmv (parts per million by volume). Both models and data indicate an amplified warming of the high latitudes during the Pliocene; however, terrestrial data suggest that Pliocene northern high-latitude temperatures were much higher than can be simulated by models. We focus on the mid-Pliocene warm period (mPWP) and show that understanding the northern high-latitude terrestrial temperatures is particularly difficult for the coldest months. Here the temperatures obtained from models and different proxies can vary by more than 20 ∘C. We refer to this mismatch as the “warm winter paradox”. Analysis suggests the warm winter paradox could be due to a number of factors including model structural uncertainty, proxy data not being strongly constrained by winter temperatures, uncertainties in data reconstruction methods, and the fact that the Pliocene northern high-latitude climate does not have a modern analogue. Refinements to model boundary conditions or proxy dating are unlikely to contribute significantly to the resolution of the warm winter paradox. For the Pliocene high-latitude terrestrial summer temperatures, models and different proxies are in good agreement. Those factors which cause uncertainty in winter temperatures are shown to be much less important for the summer. Until some of the uncertainties in winter high-latitude Pliocene temperatures can be reduced, we suggest a data–model comparison should focus on the summer. This is expected to give more meaningful and accurate results than a data–model comparison which focuses on the annual mean

The role of atmospheric CO2 in controlling sea surface temperature change during the Pliocene,

We present the role of CO2 forcing in controlling Late Pliocene sea surface temperature (SST) change using six models from Phase 2 of the Pliocene Model Intercomparison Project (PlioMIP2) and palaeoclimate proxy data from the PlioVAR working group. At a global scale, SST change in the Late Pliocene relative to the pre-industrial is predominantly driven by CO2 forcing in the low and mid-latitudes and non-CO2 forcing in the high latitudes. We find that CO2 is the dominant driver of SST change at the vast majority of proxy data sites assessed (17 out of 19), but the relative dominance of this forcing varies between all proxy sites, with CO2 forcing accounting for between 27 % and 82 % of the total change seen. The dearth of proxy data sites in the high latitudes means that only two sites assessed here are predominantly forced by non-CO2 forcing (such as changes to ice sheets and orography), both of which are in the North Atlantic Ocean.We extend the analysis to show the seasonal patterns of SST change and its drivers at a global scale and at a site-specific level for three chosen proxy data sites. We also present a new estimate of Late Pliocene climate sensitivity using site-specific proxy data values. This is the first assessment of site-specific drivers of SST change in the Late Pliocene and highlights the strengths of using palaeoclimate proxy data alongside model outputs to further develop our understanding of the Late Pliocene. We use the best available proxy and model data, but the sample sizes remain limited, and the confidence in our results would be improved with greater data availability

On the causes of mid-Pliocene warmth and polar amplification

The mid-Pliocene (~ 3 to 3.3 Ma ago), is a period of sustained global warmth in comparison to the late Quaternary (0 to ~ 1 Ma ago), and has potential to inform predictions of long-term future climate change. However, given that several processes potentially contributed, relatively little is understood about the reasons for the observed warmth, or the associated polar amplification. Here, using a modelling approach and a novel factorisation method, we assess the relative contributions to mid-Pliocene warmth from: elevated CO2, lowered orography, and vegetation and ice sheet changes. The results show that on a global scale, the largest contributor to mid-Pliocene warmth is elevated CO2. However, in terms of polar amplification, changes to ice sheets contribute significantly in the Southern Hemisphere, and orographic changes contribute significantly in the Northern Hemisphere. We also carry out an energy balance analysis which indicates that that on a global scale, surface albedo and atmospheric emmissivity changes dominate over cloud changes. We investigate the sensitivity of our results to uncertainties in the prescribed CO2 and orographic changes, to derive uncertainty ranges for the various contributing processes

Accounting for Centennial Scale Variability when Detecting Changes in ENSO: a study of the Pliocene

The El Niño Southern Oscillation (ENSO) is the dominant mode of interannual climate variability. However, climate models are inconsistent in future predictions of ENSO, and long term variations in ENSO cannot be quantified from the short instrumental records available. Here we analyse ENSO behaviour in millennial-scale climate simulations of a warm climate of the past, the mid-Pliocene Warm Period (mPWP; ∼3.3 − 3.0Ma). We consider centennial-scale variability in ENSO for both the mPWP and the preindustrial, and consider which changes between the two climates are detectable above this variability. We find that El Niño typically occurred 12% less frequently in the mPWP but with a 20% longer duration, and with stronger amplitude in precipitation and temperature. However low frequency variability in ENSO meant that Pliocene-preindustrial changes in El Niño temperature amplitude in the NINO3.4 region (5° N-5° S, 170° W-120° W) were not always detectable. The Pliocene-preindustrial El Niño temperature signal in the NINO4 region (5° N-5° S, 160° E-150° W) and the El Niño precipitation signal are usually larger than centennial scale variations of El Niño amplitude, and provide consistent indications of ENSO amplitude change. The enhanced mPWP temperature signal in the NINO4 region is associated with an increase in Central Pacific El Niño events similar to those observed in recent decades and predicted for the future. This study highlights the importance of considering centennial scale variability when comparing ENSO changes between two climate states. If centennial scale variability in ENSO has not first been established, results suggesting changes in ENSO behaviour may not be robust

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

Amplified surface warming in the south-west Pacific during the mid-Pliocene (3.3–3.0 Ma) and future implications

Based on Nationally Determined Contributions concurrent with Shared Socioeconomic Pathways (SSPs) 2-4.5, the IPCC predicts global warming of 2.1–3.5 ∘C (very likely range 10–90th percentile) by 2100 CE. However, global average temperature is a poor indicator of regional warming and global climate models (GCMs) require validation with instrumental or proxy data from geological archives to assess their ability to simulate regional ocean and atmospheric circulation, and thus, to evaluate their performance for regional climate projections. The south-west Pacific is a region that performs poorly when GCMs are evaluated against instrumental observations. The New Zealand Earth System Model (NZESM) was developed from the United Kingdom Earth System Model (UKESM) to better understand south-west Pacific response to global change, by including a nested ocean grid in the south-west Pacific with 80 % greater horizontal resolution than the global-scale host.Here, we reconstruct regional south-west Pacific sea-surface temperatures (SSTs) for the mid-Pliocene warm period (mPWP; 3.3–3.0 Ma), which has been widely considered a past analogue with an equilibrium surface temperature response of +3 ∘C to an atmospheric CO2 concentration of ∼350–400 ppm, in order to assess the warming distribution in the south-west Pacific. This study presents proxy SSTs from seven deep sea sediment cores distributed across the south-west Pacific. Our reconstructed SSTs are derived from molecular biomarkers preserved in the sediment – alkenones (i.e. U index) and isoprenoid glycerol dialkyl glycerol tetraethers (i.e. TEX86 index) – and are compared with SSTs reconstructed from the Last Interglacial (125 ka), Pliocene Model Intercomparison Project (PlioMIP) outputs and transient climate model projections (NZESM and UKESM) of low- to high-range SSPs for 2090–2099 CE.Mean interglacial equilibrium SSTs during the mPWP for the south-west Pacific sites were on average 4.2 ∘C (1.8–6.1 ∘C likely range) above pre-industrial temperatures and show good agreement with model outputs from NZESM and UKESM under mid-range SSP 2–4.6 conditions. These results highlight that not only is the mPWP an appropriate analogue when considering future temperature change in the centuries to come, but they also demonstrate that the south-west Pacific region will experience warming that exceeds that of the global mean if atmospheric CO2 remains above 350 ppm