scholarly journals Inferred basal friction and surface mass balance of the Northeast Greenland Ice Stream using data assimilation of ICESat (Ice Cloud and land Elevation Satellite) surface altimetry and ISSM (Ice Sheet System Model)

2014 ◽  
Vol 8 (6) ◽  
pp. 2335-2351 ◽  
Author(s):  
E. Larour ◽  
J. Utke ◽  
B. Csatho ◽  
A. Schenk ◽  
H. Seroussi ◽  
...  

Abstract. We present a new data assimilation method within the Ice Sheet System Model (ISSM) framework that is capable of assimilating surface altimetry data from missions such as ICESat (Ice Cloud and land Elevation Satellite) into reconstructions of transient ice flow. The new method relies on algorithmic differentiation to compute gradients of objective functions with respect to model forcings. It is applied to the Northeast Greenland Ice Stream, where surface mass balance and basal friction forcings are temporally inverted, resulting in adjusted modeled surface heights that best fit existing altimetry. This new approach allows for a better quantification of basal and surface processes and a better understanding of the physical processes currently missing in transient ice-flow models to better capture the important intra- and interannual variability in surface altimetry. It also demonstrates that large spatial and temporal variability is required in model forcings such as surface mass balance and basal friction, variability that can only be explained by including more complex processes such as snowpack compaction at the surface and basal hydrology at the bottom of the ice sheet. This approach is indeed a first step towards assimilating the wealth of high spatial resolution altimetry data available from EnviSat, ICESat, Operation IceBridge and CryoSat-2, and that which will be available in the near future with the launch of ICESat-2.

2014 ◽  
Vol 8 (3) ◽  
pp. 2331-2373 ◽  
Author(s):  
E. Larour ◽  
J. Utke ◽  
B. Csatho ◽  
A. Schenk ◽  
H. Seroussi ◽  
...  

Abstract. We present a new data assimilation method within the ISSM framework that is capable of assimilating surface altimetry data from missions such as ICESat-1 into reconstructions of transient ice flow. The new method relies on algorithmic differentiation to compute gradients of diagnostics with respect to model forcings. It is applied to the North East Greenland Ice Stream where surface mass balance and basal friction forcings are temporally inverted, resulting in significantly improved modeled surface heights that match existing altimetry. This new approach allows for a better quantification of basal and surface processes, and a better understanding of the physical processes currently missing in transient ice flow models to better capture the important intra and inter-annual variability in surface altimetry. It also demonstrates that large spatial and temporal variability is required in model forcings such as surface mass balance and basal friction, variability that can only be explained by including more complex processes such as snowpack compaction at the surface and basal hydrology at the bottom of the ice sheet. This approach is indeed a first step towards assimilating the wealth of high spatial resolution altimetry data available from EnviSat, ICESat-1, Operation IceBridge and CryoSat-2, and that will be available in the near future with the launch of ICESat-2.


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.


2014 ◽  
Vol 27 (1) ◽  
pp. 215-226 ◽  
Author(s):  
Miren Vizcaíno ◽  
William H. Lipscomb ◽  
William J. Sacks ◽  
Michiel van den Broeke

Abstract This study presents the first twenty-first-century projections of surface mass balance (SMB) changes for the Greenland Ice Sheet (GIS) with the Community Earth System Model (CESM), which includes a new ice sheet component. For glaciated surfaces, CESM includes a sophisticated calculation of energy fluxes, surface albedo, and snowpack hydrology (melt, percolation, refreezing, etc.). To efficiently resolve the high SMB gradients at the ice sheet margins and provide surface forcing at the scale needed by ice sheet models, the SMB is calculated at multiple elevations and interpolated to a finer 5-km ice sheet grid. During a twenty-first-century simulation driven by representative concentration pathway 8.5 (RCP8.5) forcing, the SMB decreases from 372 ± 100 Gt yr−1 in 1980–99 to −78 ± 143 Gt yr−1 in 2080–99. The 2080–99 near-surface temperatures over the GIS increase by 4.7 K (annual mean) with respect to 1980–99, only 1.3 times the global increase (+3.7 K). Snowfall increases by 18%, while surface melt doubles. The ablation area increases from 9% of the GIS in 1980–99 to 28% in 2080–99. Over the ablation areas, summer downward longwave radiation and turbulent fluxes increase, while incoming shortwave radiation decreases owing to increased cloud cover. The reduction in GIS-averaged July albedo from 0.78 in 1980–99 to 0.75 in 2080–99 increases the absorbed solar radiation in this month by 12%. Summer warming is strongest in the north and east of Greenland owing to reduced sea ice cover. In the ablation area, summer temperature increases are smaller due to frequent periods of surface melt.


2018 ◽  
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. 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 to global sea level rise between 1.9 and 13.0 cm 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 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.


1977 ◽  
Vol 18 (80) ◽  
pp. 359-371 ◽  
Author(s):  
I. M. Whillans

Abstract The continuity relationship that is often used in the study of ice sheets and ice shelves is developed by integrating the equation of continuity through the ice thickness. This equation is then integrated again with respect to horizontal distance from an ice divide, showing that the difference between the true ice velocity and the balance velocity, which is defined, is a measure of the time chance of the mass of a column through the ice thickness. The relationship is applied using data from along the “Byrd” station strain network, Antarctica. This region is found to be thinning slowly (0.03 m a−1 of ice of mean density) and uniformly, but it is still close to steady-state. The calculations would show a larger thinning rate if bottom sliding contributed more to the ice movement and integral shear contributed less, but the “Byrd” station bore-hole tilting results of Garfield and Ueda (1975, 1976), together with surface velocity measurements at “Byrd” station, indicate that most of the ice flow is by deformation within the ice mass. This large amount of internal deformation is more than that predicted by most “flow laws”, probably because of the strongly oriented ice-crystal fabric in the ice sheet. The cause of ice thinning is probably decreased surface mass balance beginning before A.D. 1550. The consistent relationship between measured velocity and balance velocity indicates that the ice flow is simple and that flow lines are in the same direction at depth as at the surface when considered smoothed over a distance of 10 km. Because the ice sheet is at present thinning, the balance velocity, calculated only from flow line and surface mass-balance data, and the somewhat mistaken assumption of steady-state is 15% less than the true ice velocity. This rather small difference confirms the use of balance-velocity estimates where velocity measurements are not available.


2012 ◽  
Vol 6 (6) ◽  
pp. 1561-1576 ◽  
Author(s):  
F. Gillet-Chaulet ◽  
O. Gagliardini ◽  
H. Seddik ◽  
M. Nodet ◽  
G. Durand ◽  
...  

Abstract. Over the last two decades, the Greenland ice sheet (GrIS) has been losing mass at an increasing rate, enhancing its contribution to sea-level rise (SLR). The recent increases in ice loss appear to be due to changes in both the surface mass balance of the ice sheet and ice discharge (ice flux to the ocean). Rapid ice flow directly affects the discharge, but also alters ice-sheet geometry and so affects climate and surface mass balance. Present-day ice-sheet models only represent rapid ice flow in an approximate fashion and, as a consequence, have never explicitly addressed the role of ice discharge on the total GrIS mass balance, especially at the scale of individual outlet glaciers. Here, we present a new-generation prognostic ice-sheet model which reproduces the current patterns of rapid ice flow. This requires three essential developments: the complete solution of the full system of equations governing ice deformation; a variable resolution unstructured mesh to resolve outlet glaciers and the use of inverse methods to better constrain poorly known parameters using observations. The modelled ice discharge is in good agreement with observations on the continental scale and for individual outlets. From this initial state, we investigate possible bounds for the next century ice-sheet mass loss. We run sensitivity experiments of the GrIS dynamical response to perturbations in climate and basal lubrication, assuming a fixed position of the marine termini. We find that increasing ablation tends to reduce outflow and thus decreases the ice-sheet imbalance. In our experiments, the GrIS initial mass (im)balance is preserved throughout the whole century in the absence of reinforced forcing, allowing us to estimate a lower bound of 75 mm for the GrIS contribution to SLR by 2100. In one experiment, we show that the current increase in the rate of ice loss can be reproduced and maintained throughout the whole century. However, this requires a very unlikely perturbation of basal lubrication. From this result we are able to estimate an upper bound of 140 mm from dynamics only for the GrIS contribution to SLR by 2100.


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.


2020 ◽  
Vol 14 (4) ◽  
pp. 1425-1435 ◽  
Author(s):  
Brice Noël ◽  
Leonardus van Kampenhout ◽  
Willem Jan van de Berg ◽  
Jan T. M. Lenaerts ◽  
Bert Wouters ◽  
...  

Abstract. We present a reconstruction of historical (1950–2014) surface mass balance (SMB) of the Greenland ice sheet (GrIS) using a high-resolution regional climate model (RACMO2; ∼11 km) to dynamically downscale the climate of the Community Earth System Model version 2 (CESM2; ∼111 km). After further statistical downscaling to 1 km spatial resolution, evaluation using in situ SMB measurements and remotely sensed GrIS mass change shows good agreement. Comparison with an ensemble of previously conducted RACMO2 simulations forced by climate reanalysis demonstrates that the current product realistically represents the long-term average and variability of individual SMB components and captures the recent increase in meltwater runoff that accelerated GrIS mass loss. This means that, for the first time, climate forcing from an Earth system model (CESM2), which assimilates no observations, can be used without additional corrections to reconstruct the historical GrIS SMB and its recent decline that initiated mass loss in the 1990s. This paves the way for attribution studies of future GrIS mass loss projections and contribution to sea level rise.


2019 ◽  
Vol 13 (12) ◽  
pp. 3193-3208 ◽  
Author(s):  
Raymond Sellevold ◽  
Leonardus van Kampenhout ◽  
Jan T. M. Lenaerts ◽  
Brice Noël ◽  
William H. Lipscomb ◽  
...  

Abstract. The modeling of ice sheets in Earth system models (ESMs) is an active area of research with applications to future sea level rise projections and paleoclimate studies. A major challenge for surface mass balance (SMB) modeling with ESMs arises from their coarse resolution. This paper evaluates the elevation class (EC) method as an SMB downscaling alternative to the dynamical downscaling of regional climate models. To this end, we compare EC-simulated elevation-dependent surface energy and mass balance gradients from the Community Earth System Model 1.0 (CESM1.0) with those from the regional climate model RACMO2.3. The EC implementation in CESM1.0 combines prognostic snow albedo, a multilayer snow model, and elevation corrections for two atmospheric forcing variables: temperature and humidity. Despite making no corrections for incoming radiation and precipitation, we find that the EC method in CESM1.0 yields similar SMB gradients to RACMO2.3, in part due to compensating biases in snowfall, surface melt, and refreezing gradients. We discuss the sensitivity of the results to the lapse rate used for the temperature correction. We also evaluate the impact of the EC method on the climate simulated by the ESM and find minor cooling over the Greenland ice sheet and Barents and Greenland seas, which compensates for a warm bias in the ESM due to topographic smoothing. Based on our diagnostic procedure to evaluate the EC method, we make several recommendations for future implementations.


2013 ◽  
Vol 26 (20) ◽  
pp. 7793-7812 ◽  
Author(s):  
Miren Vizcaíno ◽  
William H. Lipscomb ◽  
William J. Sacks ◽  
Jan H. van Angelen ◽  
Bert Wouters ◽  
...  

Abstract The modeling of the surface mass balance (SMB) of the Greenland Ice Sheet (GIS) requires high-resolution models in order to capture the observed large gradients in the steep marginal areas. Until now, global climate models have not been considered suitable to model ice sheet SMB owing to model biases and insufficient resolution. This study analyzes the GIS SMB simulated for the period 1850–2005 by the Community Earth System Model (CESM), which includes a new ice sheet component with multiple elevation classes for SMB calculations. The model is evaluated against observational data and output from the regional model Regional Atmospheric Climate Model version 2 (RACMO2). Because of a lack of major climate biases, a sophisticated calculation of snow processes (including surface albedo evolution) and an adequate downscaling technique, CESM is able to realistically simulate GIS surface climate and SMB. CESM SMB agrees reasonably well with in situ data from 475 locations (r = 0.80) and output from RACMO2 (r = 0.79). The simulated mean SMB for 1960–2005 is 359 ± 120 Gt yr−1 in the range of estimates from regional climate models. The simulated seasonal mass variability is comparable with mass observations from the Gravity Recovery and Climate Experiment (GRACE), with synchronous annual maximum (May) and minimum (August–September) and similar amplitudes of the seasonal cycle. CESM is able to simulate the bands of precipitation maxima along the southeast and northwest margins, but absolute precipitation rates are underestimated along the southeastern margin and overestimated in the high interior. The model correctly simulates the major ablation areas. Total refreezing represents 35% of the available liquid water (the sum of rain and melt).


Sign in / Sign up

Export Citation Format

Share Document