Parallel Ice Sheet Model Research Articles

Geochronology and ice‐flow modelling of the Late Quaternary glaciers on Mt. Soğanlı, Türkiye

ABSTRACTUnderstanding Earth's climate history through the chronology and reconstruction of palaeoglaciers is a central topic of palaeoclimatology. Examining the remnants of past glaciations preserved in isolated alpine environments provides insight into critical palaeoclimatic conditions. Accurate dating techniques, such as terrestrial cosmogenic nuclides, are crucial for understanding the timing of these changes. The Anatolian Peninsula, characterized by diverse topographical and climatic factors, witnessed the presence of glaciers during the Late Pleistocene, shaping its alpine landscape significantly. While mountain ranges such as the Taurus and Pontic favoured glacier development due to their elevation and lower temperatures, isolated mountains such as Mount Soğanlı in the interior hinterland also supported significant palaeoglaciers. Despite previous studies on Anatolia's glaciers, numerical dating of glacial deposits on Mount Soğanlı has remained elusive until now, hindering precise palaeoclimatic interpretations. This study employs cosmogenic 36Cl surface exposure dating and physical‐based ice‐flow modelling using an open‐source framework, Parallel Ice Sheet Model (PISM), to determine the timing of glaciations and to reconstruct past glacier extents on Mount Soğanlı. The results reveal glacier retreat ages of 48.3 ± 13.4 and 18.3 ± 4.4 ka based on two well‐preserved moraines. Additionally, palaeoclimate simulations matched with the field observations indicate considerably colder conditions (6.65–8.15°C colder than today) necessary to sustain glaciers on Mount Soğanlı, when precipitation amounts were the same as today. These findings contribute to a better understanding of Anatolia's palaeoclimatic fluctuations and align with similar studies in neighbouring regions, enhancing our knowledge of past climate dynamics of the Eastern Mediterranean and the Middle East.

Feedback mechanisms controlling Antarctic glacial-cycle dynamics simulated with a coupled ice sheet–solid Earth model

Abstract. The dynamics of the ice sheets on glacial timescales are highly controlled by interactions with the solid Earth, i.e., the glacial isostatic adjustment (GIA). Particularly at marine ice sheets, competing feedback mechanisms govern the migration of the ice sheet's grounding line (GL) and hence the ice sheet stability. For this study, we developed a coupling scheme and performed a suite of coupled ice sheet–solid Earth simulations over the last two glacial cycles. To represent ice sheet dynamics we apply the Parallel Ice Sheet Model (PISM), and to represent the solid Earth response we apply the 3D VIscoelastic Lithosphere and MAntle model (VILMA), which, in addition to load deformation and rotation changes, considers the gravitationally consistent redistribution of water (the sea-level equation). We decided on an offline coupling between the two model components. By convergence of trajectories of the Antarctic Ice Sheet deglaciation we determine optimal coupling time step and spatial resolution of the GIA model and compare patterns of inferred relative sea-level change since the Last Glacial Maximum with the results from previous studies. With our coupling setup we evaluate the relevance of feedback mechanisms for the glaciation and deglaciation phases in Antarctica considering different 3D Earth structures resulting in a range of load-response timescales. For rather long timescales, in a glacial climate associated with the far-field sea-level low stand, we find GL advance up to the edge of the continental shelf mainly in West Antarctica, dominated by a self-amplifying GIA feedback, which we call the “forebulge feedback”. For the much shorter timescale of deglaciation, dominated by the marine ice sheet instability, our simulations suggest that the stabilizing sea-level feedback can significantly slow down GL retreat in the Ross sector, which is dominated by a very weak Earth structure (i.e., low mantle viscosity and thin lithosphere). This delaying effect prevents a Holocene GL retreat beyond its present-day position, which is discussed in the scientific community and supported by observational evidence at the Siple Coast and by previous model simulations. The applied coupled framework, PISM–VILMA, allows for defining restart states to run multiple sensitivity simulations from. It can be easily implemented in Earth system models (ESMs) and provides the tools to gain a better understanding of ice sheet stability on glacial timescales as well as in a warmer future climate.

Simulation of a former ice field with Parallel Ice Sheet Model – Snežnik study case

Abstract. In this paper, we present a reconstruction of climate conditions during the Last Glacial Maximum on a karst plateau Snežnik, which lies in Dinaric Mountains (southern Slovenia) and bears evidence of glaciation. The reconstruction merges geomorphological ice limits, classified as either clear or unclear, and a computer modelling approach based on the Parallel Ice Sheet Model (PISM). Based on extensive numerical experiments where we studied the agreements between simulated and geomorphological ice extent, we propose using a combination of a high-resolution precipitation model that accounts for orographic precipitation combined with a simple elevation-based temperature model. The geomorphological ice extent can be simulated with climate to be around 6 °C colder than the modern day and with a lower-than-modern-day amount of precipitation, which matches other state-of-the art climate reconstructions for the era. The results indicate that an orographic precipitation model is essential for the accurate simulation of the study area, with moist southern winds from the nearby Adriatic Sea having a predominant effect on the precipitation patterns. Finally, this study shows that transforming climate conditions towards wetter and warmer or drier and colder does not significantly change the conditions for glacier formation.

Coupling MAR (Modèle Atmosphérique Régional) with PISM (Parallel Ice Sheet Model) mitigates the positive melt–elevation feedback

Abstract. The Greenland Ice Sheet is a key contributor to sea level rise. By melting, the ice sheet thins, inducing higher surface melt due to lower surface elevations, accelerating the melt coming from global warming. This process is called the melt–elevation feedback and can be considered by using two types of models: either (1) atmospheric models, which can represent the surface mass balance (SMB), or SMB estimates resulting from simpler models such as positive degree day models or (2) ice sheet models representing the surface elevation evolution. The latter ones do not represent the surface mass balance explicitly as well as polar-oriented climate models. A new coupling between the MAR (Modèle Atmosphérique Régional) regional climate model and the PISM (Parallel Ice Sheet Model) ice sheet model is presented here following the CESM2 (Community Earth System Model; SSP5-8.5, Shared Socioeconomic Pathway) scenario until 2100 at the MAR lateral boundaries. The coupling is extended to 2200 with a stabilised climate (+7 ∘C compared to 1961–1990) by randomly sampling the last 10 years of CESM2 to force MAR and reaches a sea level rise contribution of 64 cm. The fully coupled simulation is compared to a one-way experiment where surface topography remains fixed in MAR. However, the surface mass balance is corrected for the melt–elevation feedback when interpolated on the PISM grid by using surface mass balance vertical gradients as a function of local elevation variations (offline correction). This method is often used to represent the melt–elevation feedback and prevents a coupling which is too expensive in computation time. In the fully coupled MAR simulation, the ice sheet morphology evolution (changing slope and reducing the orographic barrier) induces changes in local atmospheric patterns. More specifically, wind regimes are modified, as well as temperature lapse rates, influencing the melt rate through modification of sensible heat fluxes at the ice sheet margins. We highlight mitigation of the melt lapse rate on the margins by modifying the surface morphology. The lapse rates considered by the offline correction are no longer valid at the ice sheet margins. If used (one-way simulation), this correction implies an overestimation of the sea level rise contribution of 2.5 %. The mitigation of the melt lapse rate on the margins can only be corrected by using a full coupling between an ice sheet model and an atmospheric model.

The evolution of future Antarctic surface melt using PISM-dEBM-simple

Abstract. It is virtually certain that Antarctica's contribution to sea-level rise will increase with future warming, although competing mass balance processes hamper accurate quantification of the exact magnitudes. Today, ocean-induced melting underneath the floating ice shelves dominates mass losses, but melting at the surface will gain importance as global warming continues. Meltwater at the ice surface has crucial implications for the ice sheet's stability, as it increases the risk of hydrofracturing and ice-shelf collapse that could cause enhanced glacier outflow into the ocean. Simultaneously, positive feedbacks between ice and atmosphere can accelerate mass losses and increase the ice sheet's sensitivity to warming. However, due to long response times, it may take hundreds to thousands of years until the ice sheet fully adjusts to the environmental changes. Therefore, ice-sheet model simulations must be computationally fast and capture the relevant feedbacks, including the ones at the ice–atmosphere interface. Here we use the novel surface melt module dEBM-simple (a slightly modified version of the “simple” diurnal Energy Balance Model) coupled to the Parallel Ice Sheet Model (PISM, together referred to as PISM-dEBM-simple) to estimate the impact of 21st-century atmospheric warming on Antarctic surface melt and ice dynamics. As an enhancement compared to the widely adopted positive degree-day (PDD) scheme, dEBM-simple includes an implicit diurnal cycle and computes melt not only from the temperature, but also from the influence of solar radiation and changes in ice albedo, thus accounting for the melt–albedo feedback. We calibrate PISM-dEBM-simple to reproduce historical and present-day Antarctic surface melt rates given by the regional atmospheric climate model RACMO2.3p2 and use the calibrated model to assess the range of possible future surface melt trajectories under Shared Socioeconomic Pathway SSP5-8.5 warming projections until the year 2100. To investigate the committed impacts of the enhanced surface melting on the ice-sheet dynamics, we extend the simulations under fixed climatological conditions until the ice sheet has reached a state close to equilibrium with its environment. Our findings reveal a substantial surface-melt-induced speed-up in ice flow associated with large-scale elevation reductions in sensitive ice-sheet regions, underscoring the critical role of self-reinforcing ice-sheet–atmosphere feedbacks in future mass losses and sea-level contribution from the Antarctic Ice Sheet on centennial to millennial timescales.

Modeling sensitivities of thermally and hydraulically driven ice stream surge cycling

Abstract. Modeling ice sheet instabilities is a numerical challenge of potentially high real-world relevance. Yet, differentiating between the impacts of model physics, numerical implementation choices, and numerical errors is not straightforward. Here, we use an idealized North American geometry and climate representation (similarly to the HEINO (Heinrich Event INtercOmparison) experiments – Calov et al., 2010) to examine the process and numerical sensitivity of ice stream surge cycling in ice flow models. Through sensitivity tests, we identify some numerical requirements for a more robust model configuration for such contexts. To partly address model-specific dependencies, we use both the Glacial Systems Model (GSM) and the Parallel Ice Sheet Model (PISM). We show that modeled surge characteristics are resolution dependent, though they converge (decreased differences between resolutions) at finer horizontal grid resolutions. Discrepancies between fine and coarse horizontal grid resolutions can be reduced by incorporating sliding at sub-freezing temperatures. The inclusion of basal hydrology increases the ice volume lost during surges, whereas the dampening of basal-temperature changes due to a bed thermal model leads to a decrease.

Coupled climate-glacier modelling of the last glaciation in the Alps

Abstract Our limited knowledge of the climate prevailing over Europe during former glaciations is the main obstacle to reconstruct the past evolution of the ice coverage over the Alps by numerical modelling. To address this challenge, we perform a two-step modelling approach: First, a regional climate model is used to downscale the time slice simulations of a global earth system model in high resolution, leading to climate snapshots during the Last Glacial Maximum (LGM) and the Marine Isotope Stage 4 (MIS4). Second, we combine these snapshots and a climate signal proxy to build a transient climate over the last glacial period and force the Parallel Ice Sheet Model to simulate the dynamical evolution of glaciers in the Alps. The results show that the extent of modelled glaciers during the LGM agrees with several independent key geological imprints, including moraine-based maximal reconstructed glacial extents, known ice transfluences and trajectories of erratic boulders of known origin and deposition. Our results highlight the benefit of multiphysical coupled climate and glacier transient modelling over simpler approaches to help reconstruct paleo glacier fluctuations in agreement with traces they have left on the landscape.

The stability of present-day Antarctic grounding lines – Part 1: No indication of marine ice sheet instability in the current geometry

Abstract. Theoretical and numerical work has shown that under certain circumstances grounding lines of marine-type ice sheets can enter phases of irreversible advance and retreat driven by the marine ice sheet instability (MISI). Instances of such irreversible retreat have been found in several simulations of the Antarctic Ice Sheet. However, it has not been assessed whether the Antarctic grounding lines are already undergoing MISI in their current position. Here, we conduct a systematic numerical stability analysis using three state-of-the-art ice sheet models: Úa, Elmer/Ice, and the Parallel Ice Sheet Model (PISM). For the first two models, we construct steady-state initial configurations whereby the simulated grounding lines remain at the observed present-day positions through time. The third model, PISM, uses a spin-up procedure and historical forcing such that its transient state is close to the observed one. To assess the stability of these simulated states, we apply short-term perturbations to submarine melting. Our results show that the grounding lines around Antarctica migrate slightly away from their initial position while the perturbation is applied, and they revert once the perturbation is removed. This indicates that present-day retreat of Antarctic grounding lines is not yet irreversible or self-sustained. However, our accompanying paper (Part 2, Reese et al., 2023a) shows that if the grounding lines retreated further inland, under present-day climate forcing, it may lead to the eventual irreversible collapse of some marine regions of West Antarctica.

Effects of extreme melt events on ice flow and sea level rise of the Greenland Ice Sheet

Abstract. Over the past decade, Greenland has experienced several extreme melt events, the most pronounced ones in the years 2010, 2012 and 2019. With progressing climate change, such extreme melt events can be expected to occur more frequently and potentially become more severe and persistent. So far, however, projections of ice loss and sea level change from Greenland typically rely on scenarios which only take gradual changes in the climate into account. Using the Parallel Ice Sheet Model (PISM), we investigate the effect of extreme melt events on the overall mass balance of the Greenland Ice Sheet and the changes in ice flow, invoked by the altered surface topography. As a first constraint, this study estimates the overall effect of extreme melt events on the cumulative mass loss of the Greenland Ice Sheet. We find that the sea level contribution from Greenland might increase by 2 to 45 cm (0.2 % to 14 %) by the year 2300 if extreme events occur more frequently in the future under a Representative Concentration Pathway 8.5 (RCP8.5) scenario, and the ice sheet area might be reduced by an additional 6000 to 26 000 km2 by 2300 in comparison to future warming scenarios without extremes. In conclusion, projecting the future sea level contribution from the Greenland Ice Sheet requires consideration of the changes in both the frequency and intensity of extreme events. It is crucial to individually address these extremes at a monthly resolution as temperature forcing with the same excess temperature but evenly distributed over longer timescales (e.g., seasonal) leads to less sea level rise than for the simulations of the resolved extremes. Extremes lead to additional mass loss and thinning. This, in turn, reduces the driving stress and surface velocities, ultimately dampening the ice loss attributed to ice flow and discharge. Overall, we find that the surface elevation feedback largely amplifies melting for scenarios with and without extremes, with additional mass loss attributed to this feedback having the greatest impact on projected sea level.

The impact of spatially varying ice sheet basal conditions on sliding at glacial time scales

AbstractSpatially variable basal conditions are thought to govern how ice sheets behave at glacial time scales (>1000 years) and responsible for changes in dynamics between the core and peripheral regions of the Laurentide and Fennoscandian ice sheets. Basal motion is accomplished via the deformation of unconsolidated sediments, or via sliding of the ice over an undeformable bed. We present an ice sheet sliding module for the Parallel Ice Sheet Model (PISM) that takes into account changes in sediment cover and incorporates surface meltwater. This model routes meltwater, produced at the surface and base of the ice sheet, toward the margin of the ice sheet. Basal sliding is accomplished through the deformation of water saturated sediments, or sliding at the ice-bed interface. In areas with continuous, water saturated sediments, sliding is almost always accomplished through sediment deformation. In areas with incomplete cover, sliding has a stronger dependence on the supply of water. We find that the addition of surface meltwater to the base is a more important factor for ice sheet evolution than the style of sliding. In a glacial cycle simulation, our model causes a more rapid buildup of the Laurentide Ice Sheet.

An ensemble of Antarctic deglacial simulations constrained by geological observations

The Antarctic ice sheet has the potential to make a significant contribution to future sea-level rise. Understanding this potential and making projections of future ice sheet mass change requires use of numerical models. Confidence in model projections can be improved by hindcasting: testing the model against past ice sheet changes. Robust deglacial model reconstructions are also used to correct satellite gravimetric measurements of current change for glacial isostatic adjustment processes. Here we present a new model ensemble of post-Last Glacial Maximum Antarctic deglaciation reconstructions. The ensemble is generated using the Parallel Ice Sheet Model (PISM) in which we vary a range of parameters including key glaciological controls, basal sediment strength, mantle viscosity, and climate forcing. We test the ensemble results against a database of geological constraints, and develop a new scoring scheme that allows us to determine metrics of model performance against past ice sheet grounding line extent, ice sheet thickness, and thinning rates, as well as present-day ice sheet configuration. We discuss the parameter combinations that lead to the highest-scoring simulations and we also compare our ensemble performance with existing published deglacial models, using the same scoring scheme. Exploring the characteristics of the highest-scoring ensemble members highlights some key features of deglacial behaviour including a relatively narrow range of past excess ice volumes at the LGM, Holocene retreat behind present-day grounding lines with commensurate volume minima, and readvance behaviours. The comparison also allows us to identify areas where more geological data would have high constraining power for ice sheet models.

Modelled sensitivity of Monte San Lorenzo ice cap, Patagonian Andes, to past and present climate

Sparse measurements of glacier mass balance, velocity and ice thickness in Patagonia challenge our ability to understand glacier sensitivity to climate change and relate past glacier fluctuations to palaeoclimate change. Small ice caps, such as Monte San Lorenzo, have short response times and high climate sensitivity, making well-dated moraines in their glacier foregrounds an important tool for exploring glacier response to rapid changes in palaeoclimate. Here, the Parallel Ice Sheet Model (PISM) is used to model ice flow across a domain centred on the Monte San Lorenzo ice cap. Ice-flow parameters are calibrated to match present-day ice extent, velocity and thickness. Our aim is, firstly, to quantify present-day physical glacier properties, and ice cap dynamics and sensitivities, and secondarily, to evaluate the controls on the deglaciation of the ice cap within the context of the Southern Hemisphere palaeoclimate system during the Last Glacial-Interglacial Transition (LGIT). The simulated present-day ice cap shows high surface mass flux, with ablation at outlet glacier tongues up to 18 m w. e. a−1, accumulation at the highest elevations of up to 5.5 m w. e. a−1 and a simulated Equilibrium Line Altitude (ELA) of 1750–2000 m asl. The ice cap is more sensitive to changes in precipitation relative to changes in temperature. We provide envelopes with likely ranges of palaeotemperature and palaeoprecipitation for glacial advances to moraines formed during the Last Glacial-Interglacial Transition and Holocene. Our numerical model predicts that cooling and an increase in precipitation is required to force glacial advance to mapped moraine limits at 12.1 ka (2°C cooler, 50% more precipitation), 5.6 ka (0°C cooler, 50% more precipitation) and 0.2 ka (1°C cooler, 25% more precipitation). Our modelling results thus provide insights into the present-day mass balance, thermal regime and velocity of the ice cap, explores the sensitivities of this ice cap to various model and climatic parameters, and provide palaeoclimatic envelopes for readvances during the LGIT and Holocene in Patagonia.

Reconstruction of Last Glacial Maximum glaciers and palaeoclimate in the central Taurus Range, Mt. Karanfil, of the Eastern Mediterranean

We report the first glacial chronology of Mt. Karanfil (3059 m above sea level; a.s.l.), a small mountain in the south-central Taurus Range in Turkey. Well-preserved terminal and lateral moraines are clustered in four groups; M1, M2, M3, and M4, between ∼2000 and ∼2400 m a.s.l., composed of limestone blocks. Each moraine group originated from one to three cirque areas (C1, C2, C3, and C4 cirques) on the mountain's northern face between 2400 and 2600 m a.s.l. The best estimated mean ages of moraines M1 (17.3 ± 2.4 ka), M2 (18.6 ± 1.6 ka), and M3 (20.0 ± 2.6 ka) indicate that glaciers reached their maximum extents on the global Last Glacial Maximum (LGM), and they were persistent till the Late Glacial. The last moraine group (M4), more oriented towards the north, is dated to 14.6 ± 2.4 ka, implying that the M4 glacier lasted a few thousand years more. Also, a relict rock glacier in the C1 cirque floor was dated to 9.9 ± 0.9 ka and testifies to the development of permafrost conditions at the onset of Holocene. The timing of glaciation on Mt. Karanfil is close to the global LGM, as similar terrestrial cosmogenic nuclide chronologies were obtained from other mountains in the eastern Mediterranean. We also used the Parallel Ice Sheet Model (PISM) to reconstruct the glaciers and climate during the LGM on Mt. Karanfil. Twenty-one palaeoclimate simulations were run as a function of present temperature and precipitation to reach the steady-state glacier extents, and the modelled glacial areas were compared with the field-observed ice extents. The best-fit scenarios indicate that the LGM climate on Mt. Karanfil was 8.3 ± 0.5 °C colder than today if the precipitation levels were the same as today. More humid (20% wetter) or arid (20% drier) conditions than today bring the palaeotemperature estimates to 6.9 ± 0.4 °C or 10.4 ± 0.6 °C lower than the present, respectively.

Stabilizing effect of mélange buttressing on the marine ice-cliff instability of the West Antarctic Ice Sheet

Abstract. Owing to global warming and particularly high regional ocean warming, both Thwaites and Pine Island Glaciers in the Amundsen region of the Antarctic Ice Sheet could lose their buttressing ice shelves over time. We analyse the possible consequences using the parallel ice sheet model (PISM), applying a simple cliff-calving parameterization and an ice mélange-buttressing model. We find that the instantaneous loss of ice-shelf buttressing, due to enforced ice-shelf melting, initiates grounding-line retreat and triggers marine ice sheet instability (MISI). As a consequence, the grounding line progresses into the interior of the West Antarctic Ice Sheet and leads to a sea level contribution of 0.6 m within 100 a. By subjecting the exposed ice cliffs to cliff calving using our simplified parameterization, we also analyse marine ice cliff instability (MICI). In our simulations it can double or even triple the sea level contribution depending on the only loosely constrained parameter that determines the maximum cliff-calving rate. The speed of MICI depends on this upper bound of the calving rate, which is given by the ice mélange buttressing the glacier. However, stabilization of MICI may occur for geometric reasons. Because the embayment geometry changes as MICI advances into the interior of the ice sheet, the upper bound on calving rates is reduced and the progress of MICI is slowed down. Although we cannot claim that our simulations bear relevant quantitative estimates of the effect of ice-mélange buttressing on MICI, the mechanism has the potential to stop the instability. Further research is needed to evaluate its role for the past and future evolution of the Antarctic Ice Sheet.

Brief communication: Estimating the ice thickness of the Müller Ice Cap to support selection of a drill site

Abstract. The Müller Ice Cap will soon set the scene for a new drilling project. Therefore, ice thickness estimates are necessary for planning, since thickness measurements of the ice cap are sparse. Here, three models are presented and compared: (i) a simple Semi-Empirical Ice Thickness Model (SEITMo) based on an inversion of the shallow-ice approximation by the use of a single radar line in combination with the glacier outline, surface slope, and elevation; (ii) an iterative inverse method using the Parallel Ice Sheet Model (PISM), and (iii) a velocity-based inversion of the shallow-ice approximation. The velocity-based inversion underestimates the ice thickness at the ice cap top, making the model less useful to aid in drill site selection, whereas PISM and the SEITMo mostly agree about a good drill site candidate. However, the new SEITMo is insensitive to mass balance, computationally fast, and provides as good fits as PISM.

PISM-LakeCC: Implementing an adaptive proglacial lake boundary in an ice sheet model

Abstract. During the Late Pleistocene and Holocene retreat of paleo-ice sheets in North America and Europe, vast proglacial lakes existed along the land terminating margins. These proglacial lakes impacted ice sheet dynamics by imposing boundary conditions analogous to a marine terminating margin. Such lacustrine boundary conditions cause changes in the ice sheet geometry, stress balance and frontal ablation and therefore affect the mass balance of the entire ice sheet. Despite this, dynamically evolving proglacial lakes have rarely been considered in detail in ice sheet modeling endeavors. In this study, we describe the implementation of an adaptive lake boundary in the Parallel Ice Sheet Model (PISM), which we call PISM-LakeCC. We test our model with a simplified glacial retreat setup of the Laurentide Ice Sheet (LIS). By comparing the experiments with lakes to control runs with no lakes, we show that the presence of proglacial lakes locally enhances the ice flow, which leads to a lowering of the ice sheet surface. In some cases, this also results in an advance of the ice margin and the emergence of ice lobes. In the warming climate, increased melting on the lowered ice surface drives the glacial retreat. For the LIS, the presence of lakes triggers a process similar to marine ice sheet instability, which caused the collapse of the ice saddle over Hudson Bay. In the control experiments without lakes, Hudson Bay is still glaciated when the climate reaches present-day (PD) conditions. The results of our study demonstrate that glacio-lacustrine interactions play a significant role in the retreat of land terminating ice sheet margins.

The role of an interactive Greenland ice sheet in the coupled climate-ice sheet model EC-Earth-PISM

Ice sheet processes are often simplified in global climate models as changes in ice sheets have been assumed to occur over long time scales compared to ocean and atmospheric changes. However, numerous observations show an increasing rate of mass loss from the Greenland Ice Sheet and call for comprehensive process-based models to explore its role in climate change. Here, we present a new model system, EC-Earth-PISM, that includes an interactive Greenland Ice Sheet. The model is based on the EC-Earth v2.3 global climate model in which ice sheet surface processes are introduced. This model interacts with the Parallel Ice Sheet Model (PISM) without anomaly or flux corrections. Under pre-industrial climate conditions, the modeled climate and ice sheet are stable while keeping a realistic interannual variability. In model simulations forced into a warmer climate of four times the pre-industrial CO2 concentration, the total surface mass balance decreases and the ice sheet loses mass at a rate of about 500 Gt/year. In the climate warming experiments, the resulting freshwater flux from the Greenland Ice Sheet increases 55% more in the experiments with the interactive ice sheet and the climate response is significantly different: the Arctic near-surface air temperature is lower, substantially more winter sea ice covers the northern hemisphere, and the ocean circulation is weaker. Our results indicate that the melt-albedo feedback plays a key role for the response of the ice sheet and its influence on the changing climate in the Arctic. This emphasizes the importance of including interactive ice sheets in climate change projections.

Impact of paleoclimate on present and future evolution of the Greenland Ice Sheet.

Using transient climate forcing based on simulations from the Alfred Wegener Institute Earth System Model (AWI-ESM), we simulate the evolution of the Greenland Ice Sheet (GrIS) from the last interglacial (125 ka, kiloyear before present) to 2100 AD with the Parallel Ice Sheet Model (PISM). The impact of paleoclimate, especially Holocene climate, on the present and future evolution of the GrIS is explored. Our simulations of the past show close agreement with reconstructions with respect to the recent timing of the peaks in ice volume and the climate of Greenland. The maximum and minimum ice volume at around 18-17 ka and 6-5 ka lag the respective extremes in climate by several thousand years, implying that the ice volume response of the GrIS strongly lags climatic changes. Given that Greenland's climate was getting colder from the Holocene Thermal Maximum (i.e., 8 ka) to the Pre-Industrial era, our simulation implies that the GrIS experienced growth from the mid-Holocene to the industrial era. Due to this background trend, the GrIS still gains mass until the second half of the 20th century, even though anthropogenic warming begins around 1850 AD. This is also in agreement with observational evidence showing mass loss of the GrIS does not begin earlier than the late 20th century. Our results highlight that the present evolution of the GrIS is not only controlled by the recent climate changes, but is also affected by paleoclimate, especially the relatively warm Holocene climate. We propose that the GrIS was not in equilibrium throughout the entire Holocene and that the slow response to Holocene climate needs to be represented in ice sheet simulations in order to predict ice mass loss, and therefore sea level rise, accurately.

Transient Increase in Arctic Deep-Water Formation and Ocean Circulation under Sea Ice Retreat

AbstractWhile a rapid sea ice retreat in the Arctic has become ubiquitous, the potential weakening of the Atlantic meridional overturning circulation (AMOC) in response to global warming is still under debate. As deep mixing occurs in the open ocean close to the sea ice edge, the strength and vertical extent of the AMOC is likely to respond to ongoing and future sea ice retreat. Here, we investigate the link between changes in Arctic sea ice cover and AMOC strength in a long simulation with the EC-Earth–Parallel Ice Sheet Model (PISM) climate model under the emission scenario RCP8.5. The extended duration of the experiment (years 1850–2300) captures the disappearance of summer sea ice in 2060 and the removal of winter sea ice in 2165. By introducing a new metric, the Arctic meridional overturning circulation (ArMOC), we document changes beyond the Greenland–Scotland ridge and into the central Arctic. We find an ArMOC strengthening as the areas of deep mixing move north, following the retreating winter sea ice edge into the Nansen Basin. At the same time, mixing in the Labrador and Greenland Seas reduces and the AMOC weakens. As the winter sea ice edge retreats farther into the regions with high surface freshwater content in the central Arctic Basin, the mixing becomes shallower and the ArMOC weakens. Our results suggest that the location of deep-water formation plays a decisive role in the structure and strength of the ArMOC; however, the intermittent strengthening of the ArMOC and convection north of the Greenland–Scotland ridge cannot compensate for the progressive weakening of the AMOC.

Impact of the melt–albedo feedback on the future evolution of the Greenland Ice Sheet with PISM-dEBM-simple

Abstract. Surface melting of the Greenland Ice Sheet contributes a large amount to current and future sea level rise. Increased surface melt may lower the reflectivity of the ice sheet surface and thereby increase melt rates: the so-called melt–albedo feedback describes this self-sustaining increase in surface melting. In order to test the effect of the melt–albedo feedback in a prognostic ice sheet model, we implement dEBM-simple, a simplified version of the diurnal Energy Balance Model dEBM, in the Parallel Ice Sheet Model (PISM). The implementation includes a simple representation of the melt–albedo feedback and can thereby replace the positive-degree-day melt scheme. Using PISM-dEBM-simple, we find that this feedback increases ice loss through surface warming by 60 % until 2300 for the high-emission scenario RCP8.5 when compared to a scenario in which the albedo remains constant at its present-day values. With an increase of 90 % compared to a fixed-albedo scenario, the effect is more pronounced for lower surface warming under RCP2.6. Furthermore, assuming an immediate darkening of the ice surface over all summer months, we estimate an upper bound for this effect to be 70 % in the RCP8.5 scenario and a more than 4-fold increase under RCP2.6. With dEBM-simple implemented in PISM, we find that the melt–albedo feedback is an essential contributor to mass loss in dynamic simulations of the Greenland Ice Sheet under future warming.