scholarly journals Supplementary material to "The impact of model resolution on the simulated Holocene retreat of the Southwestern Greenland Ice Sheet using the Ice Sheet System Model (ISSM)"

Author(s):  
Joshua K. Cuzzone ◽  
Nicole-Jeanne Schlegel ◽  
Mathieu Morlighem ◽  
Eric Larour ◽  
Jason P. Briner ◽  
...  
2018 ◽  
Author(s):  
Joshua K. Cuzzone ◽  
Nicole-Jeanne Schlegel ◽  
Mathieu Morlighem ◽  
Eric Larour ◽  
Jason P. Briner ◽  
...  

Abstract. Geologic archives constraining the variability of the Greenland Ice Sheet (GrIS) during the Holocene provide targets for ice sheet models to test sensitivities to variations in past climate and model formulation. Even as data-model comparisons are becoming more common, many models simulating the behaviour of the GrIS during the past rely on meshes with coarse horizontal resolution (≥ 10 km). In this study, we explore the impact of model resolution on the simulated nature of retreat across Southwestern Greenland during the Holocene. Four simulations are performed using the Ice Sheet System Model (ISSM); three which use a uniform mesh and horizontal mesh resolutions of 20 km, 10 km, and 5 km and one non-uniform mesh with resolution ranging from 2 km to 15 km. We find that the simulated retreat can vary significantly between models with different horizontal resolutions based on how well the bed topography is resolved. In areas of low topographic relief, model resolution plays a negligible role in simulated differences in retreat, with models instead responding similarly to surface mass balance (SMB) driven retreat. Conversely, in areas where the bed topography is complex and high in relief, such as fjords, lower resolution models (10 km and 20 km) simulate unrealistic retreat driven as ice-surface lowering intersects bumps in the bed topography which would otherwise be resolved as troughs in the higher resolution models. Our results highlight the important role that high resolution grids play in simulating retreat in areas of complex bed topography, but also suggest that models using non-uniform grids can save computational resources through coarsening the mesh in areas of non-complex bed topography where the SMB drives retreat. Additionally, these results emphasize that care must be taken with ice sheet models when tuning model parameters to match reconstructed margins, particularly for lower resolution models in regions where complex bed topography is poorly resolved.


2019 ◽  
Vol 13 (3) ◽  
pp. 879-893 ◽  
Author(s):  
Joshua K. Cuzzone ◽  
Nicole-Jeanne Schlegel ◽  
Mathieu Morlighem ◽  
Eric Larour ◽  
Jason P. Briner ◽  
...  

Abstract. Geologic archives constraining the variability of the Greenland ice sheet (GrIS) during the Holocene provide targets for ice sheet models to test sensitivities to variations in past climate and model formulation. Even as data–model comparisons are becoming more common, many models simulating the behavior of the GrIS during the past rely on meshes with coarse horizontal resolutions (≥10 km). In this study, we explore the impact of model resolution on the simulated nature of retreat across southwestern Greenland during the Holocene. Four simulations are performed using the Ice Sheet System Model (ISSM): three that use a uniform mesh and horizontal mesh resolutions of 20, 10, and 5 km, and one that uses a nonuniform mesh with a resolution ranging from 2 to 15 km. We find that the simulated retreat can vary significantly between models with different horizontal resolutions based on how well the bed topography is resolved. In areas of low topographic relief, the horizontal resolution plays a negligible role in simulated differences in retreat, with each model instead responding similarly to retreat driven by surface mass balance (SMB). Conversely, in areas where the bed topography is complex and high in relief, such as fjords, the lower-resolution models (10 and 20 km) simulate unrealistic retreat that occurs as ice surface lowering intersects bumps in the bed topography that would otherwise be resolved as troughs using the higher-resolution grids. Our results highlight the important role that high-resolution grids play in simulating retreat in areas of complex bed topography, but also suggest that models using nonuniform grids can save computational resources through coarsening the mesh in areas of noncomplex bed topography where the SMB predominantly drives retreat. Additionally, these results emphasize that care must be taken with ice sheet models when tuning model parameters to match reconstructed margins, particularly for lower-resolution models in regions where complex bed topography is poorly resolved.


2021 ◽  
Author(s):  
Paul Halas ◽  
Jeremie Mouginot ◽  
Basile de Fleurian ◽  
Petra Langebroek

<div> <p>Ice losses from the Greenland Ice Sheet have been increasing in the last two decades, leading to a larger contribution to the global sea level rise. Roughly 40% of the contribution comes from ice-sheet dynamics, mainly regulated by basal sliding. The sliding component of glaciers has been observed to be strongly related to surface melting, as water can eventually reach the bed and impact the subglacial water pressure, affecting the basal sliding.  </p> </div><div> <p>The link between ice velocities and surface melt on multi-annual time scale is still not totally understood even though it is of major importance with expected increasing surface melting. Several studies showed some correlation between an increase in surface melt and a slowdown in velocities, but there is no consensus on those trends. Moreover those investigations only presented results in a limited area over Southwest Greenland.  </p> </div><div> <p>Here we present the ice motion over many land-terminating glaciers on the Greenland Ice Sheet for the period 2000 - 2020. This type of glacier is ideal for studying processes at the interface between the bed and the ice since they are exempted from interactions with the sea while still being relevant for all glaciers since they share the same basal friction laws. The velocity data was obtained using optical Landsat 7 & 8 imagery and feature-tracking algorithm. We attached importance keeping the starting date of our image pairs similar, and avoided stacking pairs starting before and after melt seasons, resulting in multiple velocity products for each year.  </p> </div><div> <p>Our results show similar velocity trends for previously studied areas with a slowdown until 2012 followed by an acceleration. This trend however does not seem to be observed on the whole ice sheet and is probably specific to this region’s climate forcing. </p> </div><div> <p>Moreover comparison between ice velocities from different parts of Greenland allows us to observe the impact of different climatic trends on ice dynamics.</p> </div>


2021 ◽  
Author(s):  
Joanna Davies ◽  
Anders Møller Mathiasen ◽  
Kristiane Kristensen ◽  
Christof Pearce ◽  
Marit-Solveig Seidenkrantz

<p>The polar regions exhibit some of the most visible signs of climate change globally; annual mass loss from the Greenland Ice Sheet (GrIS) has quadrupled in recent decades, from 51 ± 65 Gt yr<sup>−1</sup> (1992-2001) to 211 ± 37 Gt yr<sup>−1</sup> (2002-2011). This can partly be attributed to the widespread retreat and speed-up of marine-terminating glaciers. The Zachariae Isstrøm (ZI) is an outlet glacier of the Northeast Greenland Ice Steam (NEGIS), one of the largest ice streams of the GrIS (700km), draining approximately 12% of the ice sheet interior. Observations show that the ZI began accelerating in 2000, resulting in the collapse of the floating ice shelf between 2002 and 2003. By 2014, the ice shelf extended over an area of 52km<sup>2</sup>, a 95% decrease in area since 2002, where it extended over 1040km<sup>2</sup>. Paleo-reconstructions provide an opportunity to extend observational records in order to understand the oceanic and climatic processes governing the position of the grounding zone of marine terminating glaciers and the extent of floating ice shelves. Such datasets are thus necessary if we are to constrain the impact of future climate change projections on the Arctic cryosphere.</p><p>A multi-proxy approach, involving grain size, geochemical, foraminiferal and sedimentary analysis was applied to marine sediment core DA17-NG-ST8-92G, collected offshore of the ZI, on  the Northeast Greenland Shelf. The aim was to reconstruct changes in the extent of the ZI and the palaeoceanographic conditions throughout the Early to Mid Holocene (c.a. 12,500-5,000 cal. yrs. BP). Evidence from the analysis of these datasets indicates that whilst there has been no grounded ice at the site over the last 12,500 years, the ice shelf of the ZI extended as a floating ice shelf over the site between 12,500 and 9,200 cal. yrs. BP, with the grounding line further inland from our study site. This was followed by a retreat in the ice shelf extent during the Holocene Thermal Maximum; this was likely to have been governed, in part, by basal melting driven by Atlantic Water (AW) recirculated from Svalbard or from the Arctic Ocean. Evidence from benthic foraminifera suggest that there was a shift from the dominance of AW to Polar Water at around 7,500 cal. yrs. BP, although the ice shelf did not expand again despite of this cooling of subsurface waters.</p>


2018 ◽  
Vol 12 (10) ◽  
pp. 3097-3121 ◽  
Author(s):  
Reinhard Calov ◽  
Sebastian Beyer ◽  
Ralf Greve ◽  
Johanna Beckmann ◽  
Matteo Willeit ◽  
...  

Abstract. We introduce the coupled model of the Greenland glacial system IGLOO 1.0, including the polythermal ice sheet model SICOPOLIS (version 3.3) with hybrid dynamics, the model of basal hydrology HYDRO and a parameterization of submarine melt for marine-terminated outlet glaciers. The aim of this glacial system model is to gain a better understanding of the processes important for the future contribution of the Greenland ice sheet to sea level rise under future climate change scenarios. The ice sheet is initialized via a relaxation towards observed surface elevation, imposing the palaeo-surface temperature over the last glacial cycle. As a present-day reference, we use the 1961–1990 standard climatology derived from simulations of the regional atmosphere model MAR with ERA reanalysis boundary conditions. For the palaeo-part of the spin-up, we add the temperature anomaly derived from the GRIP ice core to the years 1961–1990 average surface temperature field. For our projections, we apply surface temperature and surface mass balance anomalies derived from RCP 4.5 and RCP 8.5 scenarios created by MAR with boundary conditions from simulations with three CMIP5 models. The hybrid ice sheet model is fully coupled with the model of basal hydrology. With this model and the MAR scenarios, we perform simulations to estimate the contribution of the Greenland ice sheet to future sea level rise until the end of the 21st and 23rd centuries. Further on, the impact of elevation–surface mass balance feedback, introduced via the MAR data, on future sea level rise is inspected. In our projections, we found the Greenland ice sheet to contribute between 1.9 and 13.0 cm to global sea level rise until the year 2100 and between 3.5 and 76.4 cm until the year 2300, including our simulated additional sea level rise due to elevation–surface mass balance feedback. Translated into additional sea level rise, the strength of this feedback in the year 2100 varies from 0.4 to 1.7 cm, and in the year 2300 it ranges from 1.7 to 21.8 cm. Additionally, taking the Helheim and Store glaciers as examples, we investigate the role of ocean warming and surface runoff change for the melting of outlet glaciers. It shows that ocean temperature and subglacial discharge are about equally important for the melting of the examined outlet glaciers.


2016 ◽  
Vol 12 (12) ◽  
pp. 2195-2213 ◽  
Author(s):  
Heiko Goelzer ◽  
Philippe Huybrechts ◽  
Marie-France Loutre ◽  
Thierry Fichefet

Abstract. As the most recent warm period in Earth's history with a sea-level stand higher than present, the Last Interglacial (LIG,  ∼  130 to 115 kyr BP) is often considered a prime example to study the impact of a warmer climate on the two polar ice sheets remaining today. Here we simulate the Last Interglacial climate, ice sheet, and sea-level evolution with the Earth system model of intermediate complexity LOVECLIM v.1.3, which includes dynamic and fully coupled components representing the atmosphere, the ocean and sea ice, the terrestrial biosphere, and the Greenland and Antarctic ice sheets. In this setup, sea-level evolution and climate–ice sheet interactions are modelled in a consistent framework.Surface mass balance change governed by changes in surface meltwater runoff is the dominant forcing for the Greenland ice sheet, which shows a peak sea-level contribution of 1.4 m at 123 kyr BP in the reference experiment. Our results indicate that ice sheet–climate feedbacks play an important role to amplify climate and sea-level changes in the Northern Hemisphere. The sensitivity of the Greenland ice sheet to surface temperature changes considerably increases when interactive albedo changes are considered. Southern Hemisphere polar and sub-polar ocean warming is limited throughout the Last Interglacial, and surface and sub-shelf melting exerts only a minor control on the Antarctic sea-level contribution with a peak of 4.4 m at 125 kyr BP. Retreat of the Antarctic ice sheet at the onset of the LIG is mainly forced by rising sea level and to a lesser extent by reduced ice shelf viscosity as the surface temperature increases. Global sea level shows a peak of 5.3 m at 124.5 kyr BP, which includes a minor contribution of 0.35 m from oceanic thermal expansion. Neither the individual contributions nor the total modelled sea-level stand show fast multi-millennial timescale variations as indicated by some reconstructions.


2013 ◽  
Vol 9 (4) ◽  
pp. 1629-1643 ◽  
Author(s):  
M. Blaschek ◽  
H. Renssen

Abstract. The relatively warm early Holocene climate in the Nordic Seas, known as the Holocene thermal maximum (HTM), is often associated with an orbitally forced summer insolation maximum at 10 ka BP. The spatial and temporal response recorded in proxy data in the North Atlantic and the Nordic Seas reveals a complex interaction of mechanisms active in the HTM. Previous studies have investigated the impact of the Laurentide Ice Sheet (LIS), as a remnant from the previous glacial period, altering climate conditions with a continuous supply of melt water to the Labrador Sea and adjacent seas and with a downwind cooling effect from the remnant LIS. In our present work we extend this approach by investigating the impact of the Greenland Ice Sheet (GIS) on the early Holocene climate and the HTM. Reconstructions suggest melt rates of 13 mSv for 9 ka BP, which result in our model in an ocean surface cooling of up to 2 K near Greenland. Reconstructed summer SST gradients agree best with our simulation including GIS melt, confirming that the impact of the early Holocene GIS is crucial for understanding the HTM characteristics in the Nordic Seas area. This implies that modern and near-future GIS melt can be expected to play an active role in the climate system in the centuries to come.


2021 ◽  
Vol 15 (3) ◽  
pp. 1627-1644
Author(s):  
Andrea J. Pain ◽  
Jonathan B. Martin ◽  
Ellen E. Martin ◽  
Åsa K. Rennermalm ◽  
Shaily Rahman

Abstract. Accelerated melting of the Greenland Ice Sheet has increased freshwater delivery to the Arctic Ocean and amplified the need to understand the impact of Greenland Ice Sheet meltwater on Arctic greenhouse gas budgets. We evaluate subglacial discharge from the Greenland Ice Sheet for carbon dioxide (CO2) and methane (CH4) concentrations and δ13C values and use geochemical models to evaluate subglacial CH4 and CO2 sources and sinks. We compare discharge from southwest (a sub-catchment of the Isunnguata Glacier, sub-Isunnguata, and the Russell Glacier) and southern Greenland (Kiattut Sermiat). Meltwater CH4 concentrations vary by orders of magnitude between sites and are saturated with respect to atmospheric concentrations at Kiattut Sermiat. In contrast, meltwaters from southwest sites are supersaturated, even though oxidation reduces CH4 concentrations by up to 50 % during periods of low discharge. CO2 concentrations range from supersaturated at sub-Isunnguata to undersaturated at Kiattut Sermiat. CO2 is consumed by mineral weathering throughout the melt season at all sites; however, differences in the magnitude of subglacial CO2 sources result in meltwaters that are either sources or sinks of atmospheric CO2. At the sub-Isunnguata site, the predominant source of CO2 is organic matter (OM) remineralization. However, multiple or heterogeneous subglacial CO2 sources maintain atmospheric CO2 concentrations at Russell but not at Kiattut Sermiat, where CO2 is undersaturated. These results highlight a previously unrecognized degree of heterogeneity in greenhouse gas dynamics under the Greenland Ice Sheet. Future work should constrain the extent and controls of heterogeneity to improve our understanding of the impact of Greenland Ice Sheet melt on Arctic greenhouse gas budgets, as well as the role of continental ice sheets in greenhouse gas variations over glacial–interglacial timescales.


Sign in / Sign up

Export Citation Format

Share Document