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 degrees N) annual mean surface air temperature (SAT) increases of 3.7 to 11.6 degrees C compared to the pre-industrial period, with a multimodel mean (MMM) increase of 7.2 degrees 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 x 10(6) km(2), with a MMM anomaly of -5.6 x 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 seaice-free conditions (<= 1 x 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.

Marine outlet glaciers on Greenland are retreating, yet it is unclear if the recent fast retreat will persist, and how atmosphere and ocean warming will impact future retreat. We show how a marine outlet glacier in Hardangerfjorden retreated rapidly in response to the abrupt warming following the Younger Dryas cold period (approximately 11,600 years before present). This almost 1000 m deep fjord, with several sills at 300-500 m depth, hosted a 175 km long outlet glacier at the western rim of the Scandinavian Ice Sheet. We use a dynamic ice-flow model constrained by well-dated terminal and lateral moraines to simulate the reconstructed 500-year retreat of Hardangerfjorden glacier. The model includes an idealized oceanic and atmospheric forcing based on reconstructions, but excludes the surface mass balance-elevation feedback. Our simulations show a highly episodic retreat driven by surface melt and warming fjord waters, paced by the fjord bathymetry. Warming air and ocean temperatures by 4-5 degrees C during the period of retreat result in a 125-km retreat of Hardangerfjorden glacier in 500 years. Retreat rates throughout the deglaciation vary by an order of magnitude from 50 to 2500 m a(-1), generally close to 200 m a(-1), punctuated by brief events of swift retreat exceeding 500 m a(-1), each event lasting a few decades. We show that the fastest retreat rates occur in regions of the bed with the largest retrograde slopes; ice shelf length and fjord water depth is less important. Our results have implications for modern glacial fjord settings similar to Hardangerfjorden, where high retreat rates have been observed. Our findings imply that increasing air temperatures and warming subsurface waters in Greenland fjords will continue to drive extensive retreat of marine outlet glaciers. However, the recent high retreat rates are not expected to be sustained for longer than a few decades due to constraints by the fjord bathymetry.