Subglacial valleys preserved in the highlands of south and east Greenland record restricted ice extent during past warmer climates

The Greenland Ice Sheet is a key contributor to contemporary global sea level rise, but its long-term history and response to episodes of warming in Earth's geological past remain uncertain. The terrain covered by the ice sheet comprises ∼ 79 % of Greenland and ∼ 1.1 % of the Earth's land surface and contains geomorphological records that may provide valuable insights into past ice-sheet behaviour. Here we use ice surface morphology and radio-echo sounding data to identify ice-covered valleys within the highlands of southern and eastern Greenland and use numerical ice-sheet modelling to constrain the climatological and glaciological conditions responsible for valley incision. Our mapping reveals intricate subglacial valley networks with morphologies that are indicative of substantial glacial modification of an inherited fluvial landscape, yet many of these valleys are presently situated beneath cold-based, slow-moving (i.e. non-erosive) ice. We use the morphology of the valleys and our simple ice-sheet model experiments to infer that incision likely occurred beneath erosive mountain valley glaciers during one or more phases of Greenland's glacial history when ice was restricted to the southern and eastern highlands and when Greenland's contribution to barystatic sea level was up to +7 m relative to today. We infer that this valley incision primarily occurred prior to the growth of a continental-scale ice sheet, most likely during the late Miocene (ca. 7–5 Ma) and/or late Pliocene (ca. 3.6–2.6 Ma). Our findings therefore provide new data-based constraints on early Greenland Ice Sheet extent and dynamics that can serve as valuable boundary conditions in models of regional and global palaeoclimate during past warm periods that are important analogues for climate change in the 21st century and beyond

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

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

Orbital, tectonic and oceanographic controls on Pliocene climate and atmospheric circulation in Arctic Norway

During the Pliocene Epoch, a stronger-than-present overturning circulation has been invoked to explain the enhanced warming in the Nordic Seas region in comparison to low to mid-latitude regions. While marine records are indicative of changes in the northward heat transport via the North Atlantic Current (NAC) during the Pliocene, the long-term terrestrial climate evolution and its driving mechanisms are poorly understood. We present the first two-million-year-long Pliocene pollen record for the Nordic Seas region from Ocean Drilling Program (ODP) Hole 642B, reflecting vegetation and climate in Arctic Norway, to assess the influence of oceanographic and atmospheric controls on Pliocene climate evolution. The vegetation record reveals a long-term cooling trend in northern Norway, which might be linked to a general decline in atmospheric CO2 concentrations over the studied interval, and climate oscillations primarily controlled by precession (23 kyr), obliquity (54 kyr) and eccentricity (100 kyr) forcing. In addition, the record identifies four major shifts in Pliocene vegetation and climate mainly controlled by changes in northward heat transport via the NAC. Cool temperate (warmer than present) conditions prevailed between 5.03–4.30 Ma, 3.90–3.47 Ma and 3.29–3.16 Ma and boreal (similar to present) conditions predominated between 4.30–3.90 Ma, 3.47–3.29 and after 3.16 Ma. A distinct decline in sediment and pollen accumulation rates at c. 4.65 Ma is probably linked to changes in ocean currents, marine productivity and atmospheric circulation. Climate model simulations suggest that changes in the strength of the Atlantic Meridional Overturning Circulation during the Early Pliocene could have affected atmospheric circulation in the Nordic Seas region, which would have affected the direction of pollen transport from Scandinavia to ODP Hole 642B

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

PMIP4-CMIP6: the contribution of the Paleoclimate Modelling Intercomparison Project to CMIP6

The goal of the Palaeoclimate Modelling Intercomparison Project (PMIP) is to understand the response of the climate system to changes in different climate forcings and to feedbacks. Through comparison with observations of the environmental impacts of these climate changes, or with climate reconstructions based on physical, chemical or biological records, PMIP also addresses the issue of how well state-of-the-art models simulate climate changes. Palaeoclimate states are radically different from those of the recent past documented by the instrumental record and thus provide an out-of-sample test of the models used for future climate projections and a way to assess whether they have the correct sensitivity to forcings and feedbacks. Five distinctly different periods have been selected as focus for the core palaeoclimate experiments that are designed to contribute to the objectives of the sixth phase of the Coupled Model Intercomparison Project (CMIP6). This manuscript describes the motivation for the choice of these periods and the design of the numerical experiments, with a focus upon their novel features compared to the experiments performed in previous phases of PMIP and CMIP as well as the benefits of common analyses of the models across multiple climate states. It also describes the information needed to document each experiment and the model outputs required for analysis and benchmarking

Sensitivity of the Greenland Ice Sheet to Pliocene sea surface temperatures

The history of theGrIS (Greenland Ice Sheet), particularly in warm climates of the pre-Quaternary, is poorly known. IRD (ice-rafted debris) records suggest that the ice sheet has existed, at least transiently, since theMiocene and potentially since as long ago as the Eocene. As melting of the GrIS is a key uncertainty in future predictions of climate and sea-level, understanding its behaviour and role within the climate system during pastwarm periods could provide important constraints. The Pliocene has been identified as a key period for understanding warmer than modern climates. Detailed micropalaeontological analyses of the mid-Piacenzian Warm Period (3.264 - 3.025 Ma) have produced a series of SST (sea-surface temperature) reconstructions (PRISM2-AVE, PRISM2-MAX, PRISM2-MIN and PRISM3).Use of these different SSTswithin theHadley CentreGCM(GeneralCirculationModel) and BASISM (BritishAntarctic Survey Ice Sheet Model), consistently show large reductions of Pliocene Greenland ice volumes compared to modern. The changes in climate introduced by the use of different SST reconstructions do change the predicted ice volumes, mainly through precipitation feedbacks. However, the models show a relatively low sensitivity of modelled Greenland ice volumes to different mid-Piacenzian SST reconstructions, with the largest SST induced changes being 20% of Pliocene ice volume or less than a metre of sea-level rise

Late Pliocene Presence / Absence of IRD in the North Atlantic

The late Pliocene was the most recent time in Earth's history to have greater than Pre-Industrial atmospheric concentrations of CO2 and global mean annual temperatures, similar to climate model predictions for the end of this century, and smaller than modern ice sheets. The nature of Northern Hemisphere ice sheets during this period remains poorly constrained. For the first time, we combine outputs derived from a climate model with a thermodynamic iceberg model in order to decipher source regions of ice rafted debris (IRD) found in sediment cores from the North Atlantic Ocean in order to elucidate the nature of ice sheet variability during the late Pliocene. Specifically, we compare the geographical extent of iceberg melt predicted using four scenarios of late Pliocene Northern Hemisphere ice volume (three representing different glacial scenarios approximating Marine Isotope Stage (MIS) M2, and one representing a typical interglacial condition within the mid-Piacenzian Warm Period (mPWP)) with IRD data from core sites in the North Atlantic. Our results suggest ice volume during MIS M2 was greater than during the Pre-Industrial era but less than has been previously suggested. In addition, marine-terminating glaciers were still present during the mPWP. Iceberg model trajectories suggest East Greenland as a potential source of icebergs in accordance with the presence/absence of IRD at sediment core sites within the Labrador Sea. Other potential source regions are also suggested by the model, such as West Greenland and East Canada, and could be differentiated by geochemical provenance testing of late Pliocene North Atlantic IRD