Optimising a regional Antarctic Ice Sheet model to investigate basal conditions and initialise transient experiments

Author(s):  
Rupert Gladstone ◽  
Yufang Zhang ◽  
Thomas Zwinger ◽  
Fabien Gillet-Chaulet ◽  
Michael Wolovick ◽  
...  

<p>Computer models for ice sheet dynamics are the primary tools for making future predictions of ice sheet behaviour, marine ice sheet instability, and ice sheet contributions to sea level change.  Such modelling studies face a number of challenges, and we consider here two examples.  The dominant mode of flow for ice streams is sliding at the bed, and the physical processes that control sliding are hard to observe. Ice sheet models often prescribe basal resistance as a function of sliding velocity.  But laboratory experiments and real-world observations indicate that basal resistance is also dependent on the water pressure in the sub-glacial hydrologic system, a property that is hard to constrain.  Initialising an ice sheet model for future projections is usually implemented either by a multi-millennial spin up or else by optimisation simulations, both of which have significant drawbacks.  In particular, long spin-up simulations cannot easily ensure a close match to present-day ice geometry, and optimisations cannot easily ensure an overall ice sheet mass balance that matches the present-day mass balance.</p><p>Using a 3D Stokes-flow ice dynamic model, we carry out optimisations for two Antarctic catchments: The Pine Island Glacier (PIG) in West Antarctica and the Lambert-Amery Glacier System (LAGS) in East Antarctica.  We optimise both the basal resistance and flow enhancement in order to minimise discrepancy between modelled and observed (from satellite) horizontal velocities at the ice upper surface.  We use these optimised model configurations to estimate the transient mass trend and also look at the 3D velocity field, its sensitivity to choice of boundary conditions in the normal direction at upper and lower surfaces, and its implications for the 3D temperature structure.  These simulations provide an estimate of the present-day thermo-mechanical state of the PIG and LAGS.</p><p>We demonstrate that constraining only horizontal velocity in the optimisations can lead to unrealistic normal velocities at the upper surface.  We show that this can, in turn, strongly impact on the catchment’s total mass budget (through locally unconstrained thinning/thickening rates) and lead to a large-scale bias in temperatures simulated using the optimised model with the steady state assumption, due to unphysical advection of heat through the ice upper surface.</p><p>We employ the optimised model to estimate basal melt, due mainly to friction heat, and drive a subglacial hydrology model beneath the PIG, providing a model-based estimate of the distribution of basal water pressure.  We use this, along with simulated sliding velocity and basal resistance, to evaluate some commonly used sliding relations.</p>

1984 ◽  
Vol 30 (105) ◽  
pp. 131-139 ◽  
Author(s):  
L. W. Morland ◽  
G. D. Smith ◽  
G. S. Boulton

AbstractThe sliding law is defined as a basal boundary condition for the large-scale bulk ice flow, relating the tangential tractionτb, overburden pressurepb, and tangential velocityubon a smoothed-out mean bed contour. This effective bed is a lower boundary viewed on the scale of the bulk ice flow and is not the physical ice/rock or sediment interface. The sliding relation reflects on the same scale the complex motion taking place in the neighbourhood of the physical interface. The isothermal steady-state ice-sheet analysis of Morland and Johnson (1980, 1982) is applied to known surface profiles from the Greenland ice sheet and Devon Island ice cap, with their corresponding mass-balance distributions, to determineτb,pb, andubfor each case. These basal estimates are used in turn to construct, using least-squares correlation, polynomial representations for an overburden dependenceλ(pb) in the adopted form of sliding lawτb═λ(pb)ub1/mwithm ≥1.The two different data sets determine functionsλ(pb) of very different magnitudes, reflecting very different basal conditions. A universal sliding law must therefore contain more general dependence on basal conditions, but the two relations determined appear to describe the two extremes. Hence use of both relations in turn to determine profiles compatible with given mass-balance distributions can be expected to yield extremes of the possible profiles, and further to show the sensitivity of profile form to variation of the sliding relation. The theory is designed as a basis for reconstruction of former ice sheets and their dynamics which are related to the two fundamental determinants of surface mass balance and basal boundary condition.


2016 ◽  
Vol 57 (73) ◽  
pp. 69-78 ◽  
Author(s):  
Christopher M. Little ◽  
Nathan M. Urban

ABSTRACTProjections of ice-sheet mass balance require regional ocean warming projections derived from atmosphere-ocean general circulation models (AOGCMs). However, the coarse resolution of AOGCMs: (1) may lead to systematic or AOGCM-specific biases and (2) makes it difficult to identify relevant water masses. Here, we employ a large-scale metric of Antarctic Shelf Bottom Water (ASBW) to investigate circum-Antarctic temperature biases and warming projections in 19 different Coupled Model Intercomparison Project Phase 5 (CMIP5) AOGCMs forced with two different ‘representative concentration pathways’ (RCPs). For high-emissions RCP 8.5, the ensemble mean 21st century ASBW warming is 0.66, 0.74 and 0.58°C for the Amundsen, Ross and Weddell Seas (AS, RS and WS), respectively. RCP 2.6 ensemble mean projections are substantially lower: 0.21, 0.26, and 0.19°C. All distributions of regional ASBW warming are positively skewed; for RCP 8.5, four AOGCMs project warming of greater than 1.8°C in the RS. Across the ensemble, there is a strong, RCP-independent, correlation between WS and RS warming. AS warming is more closely linked to warming in the Southern Ocean. We discuss possible physical mechanisms underlying the spatial patterns of warming and highlight implications of these results on strategies for forcing ice-sheet mass balance projections.


2020 ◽  
Author(s):  
Emma Lewington ◽  
Stephen Livingstone ◽  
Chris Clark ◽  
Andrew Sole ◽  
Robert Storrar

<p>Despite being widely studied, subglacial meltwater landforms are typically mapped and investigated individually, thus the drainage system as a whole remains poorly understood. Here, we identify and map all visible traces of subglacial meltwater flow across the Keewatin sector of the former Laurentide Ice Sheet from the ArcticDEM, generating significant new insights into the connectedness of the drainage system.</p><p>Due to similarities in spacing, morphometry and spatial location, we suggest that the 100s-1000s m wide features often flanking and connecting sections of eskers (i.e. tunnel valleys, meltwater tracks and esker splays) are varying expressions of the same phenomena and collectively term these features ‘meltwater corridors’. Based on observations from contemporary ice masses, we propose a new formation model based on the pressure fluctuations surrounding a central conduit, in which the esker records the imprint of the central conduit and the wider meltwater corridors the interactions with the surrounding distributed drainage system, or variable pressure axis (VPA).</p><p>We suggest that the widespread aerial coverage of meltwater corridors across the Keewatin sector provides constraints on the extent of basal uncoupling induced by basal water pressure fluctuation and variations in spatial distribution and evolution of the subglacial drainage system, which have important implications for ice sheet dynamics. </p>


2018 ◽  
Vol 31 (21) ◽  
pp. 8895-8915 ◽  
Author(s):  
Michael R. Gallagher ◽  
Matthew D. Shupe ◽  
Nathaniel B. Miller

The Greenland Ice Sheet (GrIS) plays a crucial role in the Arctic climate, and atmospheric conditions are the primary modifier of mass balance. This analysis establishes the relationship between large-scale atmospheric circulation and principal determinants of GrIS mass balance: moisture, cloud properties, radiative forcing, and temperature. Using self-organizing maps (SOMs), observations from the Integrated Characterization of Energy, Clouds, Atmospheric State, and Precipitation at Summit (ICECAPS) project are categorized by daily sea level pressure (SLP) gradient. The results describe in detail how southerly, northerly, and zonal circulation regimes impact observations at Summit Station, Greenland. This southerly regime is linked to large anomalous increases in low-level liquid cloud formation, cloud radiative forcing (CRF), and surface warming at Summit Station. An individual southerly pattern relates to the largest positive anomalies, with the most extreme 25% of cases leading to CRF anomalies above 21 W m−2 and temperature anomalies beyond 8.5°C. Finally, the July 2012 extreme melt event is analyzed, showing that the prolonged ice sheet warming was related to persistence of these southerly circulation patterns, causing an unusually extended period of anomalous CRF and temperature. These results demonstrate a novel methodology, connecting daily atmospheric circulation to a relatively brief record of observations.


1990 ◽  
Vol 81 (4) ◽  
pp. 371-384 ◽  
Author(s):  
Richard C. A. Hindmarsh

ABSTRACTComprehensible explanations of the operation of earth climate systems should consist of descriptions of the operation of a few degrees of freedoms. Qualitative interpretations of results from large-scale numerical models generally follow this principle, but do not render formal definitions of the precise nature of such degrees of freedom.At its simplest, ice-sheet kinematics requires knowledge of the evolving height and span. Rheology and surface mass-balance impose different requirements upon the co-evolution of these variables, meaning a two-degree of freedom model is over-prescribed. By means of a perturbation expansion about the analytic similarity solution for viscous spreading, eigenfunctions corresponding to degrees of freedom in the ice-sheet profile are obtained, and are used to decompose mass-balance distributions. Only a few eigenfunctions are needed to replicate numerical models, implying that ice-sheets in plane flow may operate with fewer than ten degrees of freedom.Unstable evolution of ice-sheets can occur, when the operation of a very large number of degrees of freedom can be manifested. Previous work is reviewed and new results for the unstable transformation of valley glaciers into ice-sheets are presented. Phasing of initiation may be an unpredictable phenomenon.


2017 ◽  
Author(s):  
Sebastian H. Mernild ◽  
Glen E. Liston ◽  
Andrew P. Beckerman ◽  
Jacob C. Yde

Abstract. Knowledge about variations in runoff from Greenland to adjacent fjords and seas is important for the hydrochemistry and ocean research communities to understand the link between terrestrial and marine Arctic environments. Here, we simulate the Greenland Ice Sheet (GrIS) surface mass balance (SMB), including refreezing and retention, and runoff together with catchment-scale runoff from the entire Greenland landmass (n = 3,272 simulated catchments) throughout the 35-year period 1979–2014. SnowModel/HydroFlow was applied at 3-h intervals to resolve the diurnal cycle and at 5-km horizontal grid increments using ERA-Interim (ERA-I) reanalysis atmospheric forcing. Simulated SMB was low compared to earlier studies, whereas the GrIS surface conditions and precipitation were similar. Variations in meteorological and surface ice and snow cover conditions influenced the seasonal variability in simulated catchment runoff; variations in the GrIS internal drainage system were assumed negligible and a time-invariant digital elevation model was applied. Approximately 80 % of all catchments showed increasing runoff trends over the 35 years, with on average relatively high and low catchment-scale runoff from the SW and N parts of Greenland, respectively. Outputs from an Empirical Orthogonal Function (EOF) analysis were combined with cross-correlations indicating a direct link (zero lag time) between modeled catchment-scale runoff and variations in the large-scale atmospheric circulation indices North Atlantic Oscillation (NAO) and Atlantic Multidecadal Oscillation (AMO). This suggests that natural variabilities in AMO and NAO constitute major controls on catchment-scale runoff variations in Greenland.


2018 ◽  
Author(s):  
Leonardus van Kampenhout ◽  
Alan M. Rhoades ◽  
Adam R. Herrington ◽  
Colin M. Zarzycki ◽  
Jan T. M. Lenaerts ◽  
...  

Abstract. In this study, the resolution dependence of the simulated Greenland Ice Sheet surface mass balance in the variable-resolution Community Earth System Model (VR-CESM) is investigated. Coupled atmosphere-land simulations are performed on three regionally refined grids over Greenland at 1° (~111 km), 0.5°(~55 km), and 0.25° (~28 km), maintaining a quasi-uniform resolution of 1° (~111 km) over the rest of the globe. The SMB in the accumulation zone is significantly improved compared to airborne radar and in-situ observations, with a general wetting at the margins and a drying in the interior GrIS. Total precipitation decreases with resolution, which is in line with best-available regional climate model results. In the ablation zone, VR-CESM starts developing a positive SMB bias in some locations. Potential driving mechanisms are proposed, amongst which are diversions in large scale circulation, changes in cloud cover, and changes in summer snowfall. Overall, our results demonstrate that VR-CESM is a viable new tool in the cryospheric sciences and can be used to dynamically downscale future scenarios and/or be interactively coupled to dynamical ice sheet models.


2021 ◽  
Author(s):  
Michael McPhail ◽  
Ian Hewitt

<p>The presence of subglacial water can have a significant effect on the motion of an ice sheet. The rate at which the ice slides over the bedrock is mediated by subglacial water pressure. Meltwater on the surface of the sheet can drain through cracks and moulins; drastically increasing the amount of water under the sheet. This source of water fluctuates seasonally and diurnally, much faster than the timescale associated with large-scale glacier evolution. We are interested in the effect that this short-term variation in the subglacial hydrology, and therefore water pressure, has on the long-term behaviour of the ice sheet.<span>  </span>In particular, we are interested in how important it is to resolve the short-timescale variations in ice sliding speed.</p><p> </p><p>We use a mathematical model to study the response of the subglacial drainage system to time-varying surface melt input. By coupling this subglacial hydrology through an effective-pressure-dependent sliding law to the momentum equation for the overlying ice sheet, we study the impact of short-term meltwater fluctuations on the ice velocity.<span>  </span>We study these interactions using a one-dimensional (1D) flowline model representing a confined glacier, allowing us to explore a range of couplings between the ice flow and hydrology.<span>  </span>This enables us to assess the influence of the fluctuating meltwater input on the long-term behaviour of the ice sheet. We find that using a time-averaged effective pressure with an asynchronous coupling to the momentum equation gives a reasonable estimate for the time-averaged ice-sheet velocity, despite the nonlinearity of the governing equations. We use the results to suggest how hydrological coupling might be achieved in larger-scale models where resolving the short-term fluctuations is likely to be infeasible. <span> </span></p>


2000 ◽  
Vol 31 ◽  
pp. 133-140 ◽  
Author(s):  
Ives Janssens ◽  
Philippe Huybrechts

AbstractRetention of meltwater runoff by percolation and/or refreezing in the snowpack cannot be neglected when studying the surface mass balance of the Greenland ice sheet. In this paper, we make a detailed comparison of several treatments proposed in the literature to account for this process in large-scale mass-balance parameterizations. The melt on the Greenland ice sheet is calculated with a revised degree-day model using updated datasets of surface elevation and precipitation rate on a 5 km grid. Crucial model parameters are recalibrated by comparing mass-balance characteristics with available observations on a regional basis. We discuss the role of meltwater retention in the light of the overall mass balance of the Greenland ice sheet and its sensitivity to climatic change, and display patterns of effective-retention fractions for the various methods. As a main conclusion it appears that overall results are quite similar for the various models, but that meltwater retention has a large spatial variation not described by the simple treatments. Using the most comprehensive retention model, the sensitivity of the runoff is found to be +0.35 mm ˚C–1 of sea-level change per year. We also present a new map of the different zones (facies) that characterize the accumulation area of the Greenland ice sheet, which is useful for interpreting field data and calibrating satellite observations.


Sign in / Sign up

Export Citation Format

Share Document