Interpretation of short-term ice-sheet elevation changes inferred from satellite altimetry

Satellite altimetry offers means of directly measuring changes in surface elevation over the polar ice sheets of Greenland and Antarctica. By relating these changes to variations in ice mass, it becomes possible to detect short-term changes in the Earth's ice sheets. However, it is not immediately obvious that short-term changes in surface elevation are indicative of any (long-term) trend in ice mass. An increase in ice thickness may very well reflect the response of the glacier to random fluctuations in precipitation. The spectrum of this response is dominated by low frequencies, with the majority of the variance contained in the longer time scales. As a result, the ice-thickness record may exhibit trends that have no climatic significance, but are due to a low-frequency response to random forcing. A simple model for the interpretation of observed elevation changes is developed and applied to measurements made over the Greenland Ice Sheet. It appears to be unlikely that the difference between the rate of thickening derived by Zwally and others (1989) using repeat satellite altimetry, and significantly smaller previous estimates, can be explained as being the response of the ice sheet to random climatic forcing or that this difference can be attributed to a recent increase in accumulation rate.     

A continuous simulation of global ice volume over the past 1 million years with 3-D ice-sheet models

Sea-level records show large glacial-interglacial changes over the past million years, which on these time scales are related to changes of ice volume on land. During the Pleistocene, sea-level changes induced by ice volume are largely caused by the waxing and waning of the large ice sheets in the Northern Hemisphere. However, the individual contributions of ice in the Northern and Southern Hemisphere are poorly constrained. In this study, for the first time a fully coupled system of four 3-D ice-sheet models is used, simulating glaciations on Eurasia, North America, Greenland and Antarctica. The ice-sheet models use a combination of the shallow ice and shelf approximations to determine sheet, shelf and sliding velocities. The framework consists of an inverse forward modelling approach to derive a self-consistent record of temperature and ice volume from deep-sea benthic δ¹⁸O data over the past 1 million years, a proxy for ice volume and temperature. It is shown that for both eustatic sea level and sea water δ¹⁸O changes, the Eurasian and North American ice sheets are responsible for the largest part of the variability. The combined contribution of the Antarctic and Greenland ice sheets is about 10 % for sea level and about 20 % for sea water δ¹⁸O during glacial maxima. However, changes in interglacials are mainly caused by melt of the Greenland and Antarctic ice sheets, with an average time lag of 4 kyr between melt and temperature. Furthermore, we have tested the separate response to changes in temperature and sea level for each ice sheet, indicating that ice volume can be significantly influenced by changes in eustatic sea level alone. Hence, showing the importance of a simultaneous simulation of all four ice sheets. This paper describes the first complete simulation of global ice-volume variations over the late Pleistocene with the possibility to model changes above and below present-day ice volume, constrained by observations of benthic δ¹⁸O proxy data.     

GCM simulations of the Last Glacial Maximum surface climate of Greenland and Antarctica

 The LMDz variable grid GCM was used to simulate the Last Glacial Maximum (LGM, 21 ky Bp.) climate of Greenland and Antarctica at a spatial resolution of about 100 km.The high spatial resolution allows to investigate the spatial variability of surface climate change signals, and thus to address the question whether the sparse ice core data can be viewed as representative for the regional scale climate change. This study addresses primarily surface climate parameters because these can be checked against the, limited, ice core record. The changes are generally stronger for Greenland than for Antarctica, as the imposed changes of the forcing boundary conditions (e.g., sea surface temperatures) are more important in the vicinity of Greenland. Over Greenland, and to a limited extent also in Antarctica, the climate shows stronger changes in winter than in summer. The model suggests that the linear relationship between the surface temperature and inversion strength is modified during the LGM. The temperature dependency of the moisture holding capacity of the atmosphere alone cannot explain the strong reduction in snowfall over central Greenland; atmospheric circulation changes also play a crucial role. Changes in the high frequency variability of snowfall, atmospheric pressure and temperature are investigated and possible consequences for the interpretation of ice core records are discussed. Using an objective cyclone tracking scheme, the importance of changes of the atmospheric dynamics off the coasts of the ice sheets, especially for the high frequency variability of surface climate parameters, is illustrated. The importance of the choice of the LGM ice sheet topography is illustrated for Greenland, where two different topographies have been used, yielding results that differ quite strongly in certain nontrivial respects. This means that the paleo-topography is a significant source of uncertainty for the modelled paleoclimate. The sensitivity of the Greenland LGM climate to the prescribed sea surface conditions is examined by using two different LGM North Atlantic data sets. Received: 23 October 1997 / Accepted: 17 March 1998     

The CO2-induced thickening/thinning of the Greenland and Antarctic ice sheets as simulated by a GCM (CCMI) and an ice-sheet model

Scaling analysis shows that the mean thickness of an ice sheet depends on the product of two poorly known quantities, the ice viscosity and the net snow accumulation rate. We adjust the viscosity of an ice sheet in order to get a consistent value of this product for the present-day ice sheet volume and area given the net snow accumulation rate calculated by an atmospheric general circulation model (GCM). We then hold this artificial rheology constant in further numerical experiments. We hope that in doing so we can partially compensate for systematic GCM errors in simulating the snow accumulation rate, and, therefore, thickening/thinning of ice sheets will depend mostly on the tendency in the net accumulation change rather than on its absolute value. Using this approach, the response of the Greenland and Antarctic ice sheets to doubling CO₂ concentration is simulated and the horizontal distribution of possible thickening/thinning of polar ice obtained. We find that, initially, the region of thickening ice is close to the area of increased snowfall rate, but later it significantly changes under the influence of internal ice flow dynamics. The sea-level changes predicted by our experiments agree with some empirical estimates. The sensitivity experiment with assigned basal sliding does not show significant changes in the large-scale ice topography, meaning, for example, that there is no indication of a possible disintegration of the West Antarctic ice sheet. At the same time, the regional thickening/thinning of ice (and consequently the sea-level change) depends strongly on processes at the ice sheet bottom.     

Climate-model induced differences in the 21st century global and regional glacier contributions to sea-level rise

The large uncertainty in future global glacier volume projections partly results from a substantial range in future climate conditions projected by global climate models. This study addresses the effect of global and regional differences in climate input data on the projected twenty-first century glacier contribution to sea-level rise. Glacier volume changes are calculated with a surface mass balance model combined with volume-area scaling, applied to 89 glaciers in different climatic regions. The mass balance model is based on a simplified energy balance approach, with separated contributions by net solar radiation and the combined other fluxes. Future mass balance is calculated from anomalies in air temperature, precipitation and atmospheric transmissivity, taken from eight global climate models forced with the A1B emission scenario. Regional and global sea-level contributions are obtained by scaling the volume changes at the modelled glaciers to all glaciers larger than 0.1 km² outside the Greenland and Antarctic ice sheets. This results in a global value of 0.102 ± 0.028 m (multi-model mean and standard deviation) relative sea-level equivalent for the period 2012–2099, corresponding to 18 ± 5 % of the estimated total volume of glaciers. Glaciers in the Antarctic, Alaska, Central Asia and Greenland together account for 65 ± 4 % of the total multi-model mean projected sea-level rise. The projected sea-level contribution is 35 ± 17 % larger when only anomalies in air temperature are taken into account, demonstrating an important compensating effect by increased precipitation and possibly reduced atmospheric transmissivity. The variability in projected precipitation and atmospheric transmissivity changes is especially large in the Arctic regions, making the sea-level contribution for these regions particularly sensitive to the climate model used. Including additional uncertainties in the modelling procedure and the input data, the total uncertainty estimate for the future projections becomes ±0.063 m.     

A new method to estimate ice age temperatures

On glacial time scales, the waxing and waning of the Eurasian and North American ice sheets depend largely on variations in atmospheric temperature. As global sea level is primarily determined by the volume of these ice sheets, there is a direct (yet complex) relation between global sea level and the northern hemispheric (NH) temperature. This relation is essentially represented by a model of the NH ice sheets. We use a thermomechanical ice-sheet–ice shelf–bedrock model in conjunction with an inverse method to deduce a time series of NH temperature (from 120 kyr BP until present) that is consistent with the observed global sea level record. The advantage of this method is that it provides the annual mean surface air temperature averaged over the NH continents north of 40°N. The results reveal that ice age temperatures were 4–10°C lower than today, which agrees with other temperature reconstructions. However, reconstructed temperatures are comparitively low during the early stages of the glacial, a feature that is consistent with the rapid growth of the ice sheets. The sensitivity of the results for uncertainties in precipitation rate, in observed sea level and in some other model parameters is examined to quantify the error in reconstructed temperature. During the glacial period (120–15 kyr BP), surface air temperatures in the NH (north of 40°N) were 7.2±1.5°C lower than todays (interglacial) temperatures.     

A projection of future sea level

Evidence is reviewed that suggests faster sea-level rise when climate gets warmer. Four processes appear as dominating on a time scale of decades to centuries: melting of mountain glaciers and small ice caps, changes in the mass balance of the large polar ice sheets (Greenland, Antarctica), possible ice-flow instabilities (in particular on the West Antarctic Ice Sheet), and thermal expansion of ocean water.For a given temperature scenario, an attempt is made to estimate the different contributions. The calculation yields a figure of 9.5 cm of sea-level rise since 1850 AD, which is within the uncertainty range of estimates of the observed rise.A further 33 cm rise is found as most likely for the year 2050, but the uncertainty is very large ( = 32 cm). The contribution from melting of land ice is of the same order of magnitude as thermal expansion. The mass-balance effects of the major ice sheets tend to cancel to some extent (increasing accumulation on Antarctica, increasing ablation on Greenland). For the year 2100 a value of 66 cm above the present-day stand is found ( = 57 cm). The estimates of the standard deviation include uncertainty in the temperature scenario, as presented elsewhere in this volume.     

Linear response functions to project contributions to future sea level

We propose linear response functions to separately estimate the sea-level contributions of thermal expansion and solid ice discharge from Greenland and Antarctica. The response function formalism introduces a time-dependence which allows for future rates of sea-level rise to be influenced by past climate variations. We find that this time-dependence is of the same functional type, R(t) ～ t ^α, for each of the three subsystems considered here. The validity of the approach is assessed by comparing the sea-level estimates obtained via the response functions to projections from comprehensive models. The pure vertical diffusion case in one dimension, corresponding to α =  ?0.5, is a valid approximation for thermal expansion within the ocean up to the middle of the twenty first century for all Representative Concentration Pathways. The approximation is significantly improved for α =  ? 0.7. For the solid ice discharge from Greenland we find an optimal value of α =  ?0.7. Different from earlier studies we conclude that solid ice discharge from Greenland due to dynamic thinning is bounded by 0.42 m sea-level equivalent. Ice discharge induced by surface warming on Antarctica is best captured by a positive value of α = 0.1 which reflects the fact that ice loss increases with the cumulative amount of heat available for softening the ice in our model.     

Response in atmospheric circulation and sources of Greenland precipitation to glacial boundary conditions

The response in northern hemisphere atmospheric circulation and the resulting changes in moisture sources for Greenland precipitation to glacial boundary conditions are studied in NCAR’s CCM3 atmospheric general circulation model fitted with a moisture tracking functionality. We employ both the CLIMAP SST reconstruction and a modification thereto with reconstructions of glacial ice sheets and land masks. The individual components of the boundary conditions are added first one at a time and, finally, together. These steps show the atmospheric circulation to respond approximately linearly to the boundary condition changes, and the full glacial change may thus be decomposed into contributions from SST and topography changes, respectively. We find that using the CLIMAP SST reconstruction leads to a shift from Atlantic toward Pacific source regions not found with the modified reconstruction having cooler tropics and less sea ice. The occurrence of such a shift depends chiefly on the SST reconstruction and not on the existence of the large northern hemisphere glacial ice sheets. The influence of these circulation changes on important factors for ice core interpretation such as precipitation seasonality, condensation temperatures and source temperatures are assessed.     

Mid latitude winter climate variability in the South Indian and southwest Pacific regions since 1300 AD

Mid-latitude winter atmospheric variability in the South Indian Ocean and southwest Pacific Ocean regions of the circum-Antarctic are reconstructed using sea-salt aerosol concentrations measured in the high resolution Law Dome (DSS) ice core from East Antarctica. The sea-salt aerosol concentration data, as sodium (Na), were measured at approximately monthly resolution spanning the past 700 years. Analyses of covariations between Na concentrations in Law Dome ice, and mean sea-level pressure (MSLP) and wind field data were conducted to define the mid-latitude and sub-Antarctic atmospheric circulation patterns associated with variations in Na delivery. High Na concentrations in Law Dome snow are associated with increased meridional aerosol transport from mid-latitude sources. The seasonal average Na concentration for early winter (May, June, July (MJJ)) is strongly correlated to the mid-latitude MSLP field in the South Indian and southwest Pacific Oceans, and southern Australian regions. In addition, the average MJJ Na concentrations display a strong association with the stationary Rossby wave number 3 circulation, and are anti-correlated to the Southern Annular Mode (SAM) index of climate variability: high (low) Na concentrations occurring during negative (positive) SAM phases. This observed relationship is used to derive a proxy record for early-winter MSLP anomalies and the SAM in the South Indian and southwest Pacific Ocean regions over the period 1300–1995 AD. The proxy SAM index from 1300 to 1995 AD shows pronounced decadal-scale variability throughout. The period after 1500 AD is marked by a tendency toward slower variations and a weakly-positive mean SAM (enhanced westerlies in the 50° to 65°S zone) compared to the early part of the record.     

Evaluation of a high-resolution regional climate simulation over Greenland

A simulation of the 1991 summer has been performed over south Greenland with a coupled atmosphere–snow regional climate model (RCM) forced by the ECMWF re-analysis. The simulation is evaluated with in-situ coastal and ice-sheet atmospheric and glaciological observations. Modelled air temperature, specific humidity, wind speed and radiative fluxes are in good agreement with the available observations, although uncertainties in the radiative transfer scheme need further investigation to improve the model’s performance. In the sub-surface snow-ice model, surface albedo is calculated from the simulated snow grain shape and size, snow depth, meltwater accumulation, cloudiness and ice albedo. The use of snow metamorphism processes allows a realistic modelling of the temporal variations in the surface albedo during both melting periods and accumulation events. Concerning the surface albedo, the main finding is that an accurate albedo simulation during the melting season strongly depends on a proper initialization of the surface conditions which mainly result from winter accumulation processes. Furthermore, in a sensitivity experiment with a constant 0.8 albedo over the whole ice sheet, the average amount of melt decreased by more than 60%, which highlights the importance of a correctly simulated surface albedo. The use of this coupled atmosphere–snow RCM offers new perspectives in the study of the Greenland surface mass balance due to the represented feedback between the surface climate and the surface albedo, which is the most sensitive parameter in energy-balance-based ablation calculations.     

Present and future climates of the Greenland ice sheet according to the IPCC AR4 models

The atmosphere?Cocean general circulation models (AOGCMs) used for the IPCC 4th Assessment Report (IPCC AR4) are evaluated for the Greenland ice sheet (GrIS) current climate modelling. The most suited AOGCMs for Greenland climate simulation are then selected on the basis of comparison between the 1970?C1999 outputs of the Climate of the twentieth Century experiment (20C3M) and reanalyses (ECMWF, NCEP/NCAR). This comparison indicates that the representation quality of surface parameters such as temperature and precipitation are highly correlated to the atmospheric circulation (500?hPa geopotential height) and its interannual variability (North Atlantic oscillation). The outputs of the three most suitable AOGCMs for present-day climate simulation are then used to assess the changes estimated by three IPCC greenhouse gas emissions scenarios (SRES) over the GrIS for the 2070?C2099 period. Future atmospheric circulation changes are projected to dampen the zonal flow, enhance the meridional fluxes and therefore provide additional heat and moisture to the GrIS, increasing temperature over the whole ice sheet and precipitation over its northeastern area. We also show that the GrIS surface mass balance anomalies from the SRES A1B scenario amount to ?300?km³/year with respect to the 1970?C1999 period, leading to a global sea-level rise of 5?cm by the end of the 21st century. This work can help to select the boundaries conditions for AOGCMs-based downscaled future projections.     

Simulations of the Last Glacial Maximum climates using a general circulation model: prescribed versus computed sea surface temperatures

 The climate during the Last Glacial Maximum (LGM) has been simulated using the UK Universities Global Atmospheric Modelling Programme (UGAMP) general circulation model (GCM) with both prescribed sea surface temperatures (SSTs) based on the CLIMAP reconstruction and computed SSTs with a simple thermodynamic slab ocean. Consistent with the Paleoclimate Modelling Intercomparison Project (PMIP), the other boundary conditions include the large changes in ice-sheet topography and geography, a lower sea level, a lower concentration of CO₂ in the atmosphere, and a slightly different insolation pattern at the top of the atmosphere. The results are analysed in terms of changes in atmospheric circulation. Emphasis is given to the changes in surface temperatures, planetary waves, storm tracks and the associated changes in distribution of precipitation. The model responds in a similar manner to the changes in boundary conditions to previous studies in global mean statistics, but differs in its treatment of regional climates. Results also suggest that both the land ice sheets and sea ice introduce significant changes in planetary waves and transient eddy activity, which in turn affect regional climates. The computed SST simulations predict less sea ice and cooler tropical temperatures than those based on CLIMAP SSTs. It is unclear as to whether this is a model and/or a data problem, but the resulting changes in land temperatures and precipitation can be large. Snow mass budget analysis suggests that there is net ice loss along the southern edges of the Laurentide and Fennoscandian ice sheets and net ice gain over other parts of the two ice sheets. The net accumulation is mainly due to the decrease in ablation in the cold climate rather than to the changes in snowfall. The characteristics of the Greenland ice-sheet mass balance in the LGM simulations is also quite different from those in the present-day (PD) simulations. The ablation in the LGM simulations is negligible while it is a very important process in the ice mass budget in the PD simulations. Received: 10 January 1997 / Accepted: 11 December 1997     

Snowpack model estimates of the mass balance of the Greenland ice sheet and its changes over the twentyfirst century

The performance of a snow cover model in capturing the ablation on the Greenland ice sheet is evaluated. This model allows an explicit calculation of the formation of melt water, of the fraction of melt water which re-freezes, and of runoff in the ablation region. The input climate variables to the snowpack model come from two climate models. While the higher resolution general circulation model (ECHAM 4), is closest to observations in its estimate of accumulation, it fails to give accurate results in its predictions of runoff, primarily in the southern half of the ice sheet. The two-dimensional low-resolution climate model (MIT 2D LO) produces estimates of runoff from the Greenland ice sheet within the range of uncertainty of the Inter governmental Panel on Climate Change (IPCC1) 1995 estimates. Both models reproduce some of the characteristics of the extent of the wet snow zone observed with satellite remote sensing; the MIT model is closer to observations in terms of areal extent and intensity of the melting in the southern half of the ice-sheet in July and August while the ECHAM model reproduces melting in the northern half of the ice sheet well. Changes in runoff from Greenland and Antarctica are often cited as one of the major concerns linked to anthropogenic changes in climate. Because it is based on physical principles and relies on the surface energy balance as input, the snow cover model can respond to the current climatic forcing as well as to future changes in climate on the century time scale without the limitations inherent in empirical parametrizations. For a reference climate scenario similar to the IPCC's IS92a, the model projects that the Greenland ice sheet does not contribute significantly to changes in the level of the ocean over the twenty-first century. Increases in accumulation over the central portion of the ice sheet offset most of the increase in melting and runoff, which takes place along the margins of the ice sheet. The range of uncertainty in the predictions of sea-level rise is estimated by repeating the calculation with the MIT model for seven climate change scenarios. The range is –0.5 to 1.7 cm.     

Present And Past Ice-Sheet Mass Balance Simulations For Greenland And The Tibetan Plateau

Net annual mass balance was evaluated for Greenland and the Tibetan Plateau using the meteorological forcings from the NCEP reanalysis and two GCMs (FOAM1.0 and CSM1.4) for modern climate and for different time periods extending back to the beginning of the Holocene (11,000 years ago) for the climate models. The ice-sheet budget calculations, using the degree day methodology, were performed on a finer grid than the model output by interpolating monthly precipitation and surface temperature and correcting the latter to account for the GCMs smoothed topography. The computed net mass balance for Greenland in the present day is positive and it ranges between 290–300 mm water equivalent (w.e.)/year for the two models, values close to the NCEP estimate of 250 mm/year. The past climate simulations show that the Greenland mass balance has become slightly more positive since the beginning of the Holocene. The Tibetan Plateaus present-day area average net mass balance is negative and ranges between –1200 and –2000 mm w.e. /year for the two models, values bracketing the NCEP estimate of 1700 mm/year, although the balance is positive over small regions of the plateau consistent with the existence of small ice caps and glaciers. The calculated past mass balance shows an increasingly less negative value for FOAM from 11,000 years ago towards the present and expansion of the positive mass balance areas, mainly due to decreased snow ablation as the summertime insolation decreases with the changes in orbital forcing; in CSM the opposite trend occurs but changes are smaller and less systematic. The result from FOAM shows that the likelihood of ice sheets developing on the Tibetan Plateau may have increased since 11000 years ago, which is consistent with some glacial records.     

Climate sensitivity studies of the Greenland ice sheet using satellite AVHRR,SMMR, SSM/I and in situ data

Summary The feasibility of using satellite data for climate research over the Greenland ice sheet is discussed. In particular, we demonstrate the usefulness of Advanced Very High Resolution Radiometer (AVHRR) Local Area Coverage (LAC) and Global Area Coverage (GAC) data for narrow-band albedo retrieval. Our study supports the use of lower resolution AVHRR (GAC) data for process studies over most of the Greenland ice sheet. Based on LAC data time series analysis, we can resolve relative albedo changes on the order of 2–5%. In addition, we examine Scanning Multichannel Microwave Radiometer (SMMR) and Special Sensor Microwave Imager (SSM/I) passive microwave data for snow typing and other signals of climatological significance. Based on relationships between in situ measurements and horizontally polarized 19 and 37 GHz observations, wet snow regions are identified. The wet snow regions increase in aerial percentage from 9% of the total ice surface in June to a maximum of 26% in August 1990. Furthermore, the relationship between brightness temperatures and accumulation rates in the northeastern part of Greenland is described. We found a consistent increase in accumulation rate for the northeastern part of the ice sheet from 1981 to 1986.With 16 Figures     

Sea-level pressure variability around Antarctica since A.D. 1750 inferred from subantarctic tree-ring records

 A tree-ring chronology network recently developed from the subantarctic forests provides an opportunity to study long-term climatic variability at higher latitudes in the Southern Hemisphere. Fifty long (1911–1985), homogeneous records of monthly mean sea-level pressure (MSLP) from the southern latitudes (15–65^°S) were intercorrelated on a seasonal basis to establish the most consistent, long-term Trans-Polar teleconnections during this century. Variations in summer MSLP between the South America-Antarctic Peninsula and the New Zealand sectors of the Southern Ocean are significantly correlated in a negative sense (r=−0.53, P<0.001). Climatically sensitive chronologies from Tierra del Fuego (54–55^°) and New Zealand (39–47^°) were used to develop verifiable reconstructions of summer (November to February) MSLP for both sectors of the Southern Ocean. These reconstructions, which explain between 37 and 43% of the instrumentally recorded pressure variance, indicate that inverse trends in MSLP from diametrically opposite sides of Antarctica have prevailed during the past two centuries. However, the strength of this relationship varies over time. Differences in normalized MSLP between the New Zealand and the South America-Antarctic Peninsula sectors were used to develop a Summer Trans-Polar Index (STPI), which represents an index of sea-level pressure wavenumber one in the Southern Hemisphere higher latitudes. Tree-ring based reconstructions of STPI show significant differences in large-scale atmospheric circulation between the nineteenth and the twentieth centuries. Predominantly-negative STPI values during the nineteenth century are consistent with more cyclonic activity and lower summer temperatures in the New Zealand sector during the 1800s. In contrast, cyclonic activity appears to have been stronger in the mid-twentieth than previously for the South American sector of the Southern Ocean. Recent variations in MSLP in both regions are seen as part of the long-term dynamics of the atmosphere connecting opposite sides of Antarctica. A detailed analysis of the MSLP and STPI reconstructions in the time and frequency domains indicates that much of the interannual variability is principally confined to frequency bands with a period around 3.3–3.6 y. Cross spectral analysis between the STPI reconstruction and the Southern Oscillation Index suggests that teleconnections between the tropical ocean and extra-tropical MSLP variations may be influencing climate fluctuations at southern latitudes. Received: 18 December 1996/Accepted: 10 January 1997     

Effects of a melted greenland ice sheet on climate, vegetation, and the cryosphere

This paper investigates the possible implications for the earth-system of a melting of the Greenland ice-sheet. Such a melting is a possible result of increased high latitude temperatures due to increasing anthropogenic greenhouse gas emissions. Using an atmosphere-ocean general circulation model (AOGCM), we investigate the effects of the removal of the ice sheet on atmospheric temperatures, circulation, and precipitation. We find that locally over Greenland, there is a warming associated directly with the altitude change in winter, and the altitude and albedo change in summer. Outside of Greenland, the largest signal is a cooling over the Barents sea in winter. We attribute this cooling to a decrease in poleward heat transport in the region due to changes to the time mean circulation and eddies, and interaction with sea-ice. The simulated climate is used to force a vegetation model and an ice-sheet model. We find that the Greenland climate in the absence of an ice sheet supports the growth of trees in southern Greenland, and grass in central Greenland. We find that the ice sheet is likely to regrow following a melting of the Greenland ice sheet, the subsequent rebound of its bedrock, and a return to present day atmospheric CO₂ concentrations. This regrowth is due to the high altitude bedrock in eastern Greenland which allows the growth of glaciers which develop into an ice sheet.     

Regional and global projections of twenty-first century glacier mass changes in response to climate scenarios from global climate models

A large component of present-day sea-level rise is due to the melt of glaciers other than the ice sheets. Recent projections of their contribution to global sea-level rise for the twenty-first century range between 70 and 180 mm, but bear significant uncertainty due to poor glacier inventory and lack of hypsometric data. Here, we aim to update the projections and improve quantification of their uncertainties by using a recently released global inventory containing outlines of almost every glacier in the world. We model volume change for each glacier in response to transient spatially-differentiated temperature and precipitation projections from 14 global climate models with two emission scenarios (RCP4.5 and RCP8.5) prepared for the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. The multi-model mean suggests sea-level rise of 155 ± 41 mm (RCP4.5) and 216 ± 44 mm (RCP8.5) over the period 2006–2100, reducing the current global glacier volume by 29 or 41 %. The largest contributors to projected global volume loss are the glaciers in the Canadian and Russian Arctic, Alaska, and glaciers peripheral to the Antarctic and Greenland ice sheets. Although small contributors to global volume loss, glaciers in Central Europe, low-latitude South America, Caucasus, North Asia, and Western Canada and US are projected to lose more than 80 % of their volume by 2100. However, large uncertainties in the projections remain due to the choice of global climate model and emission scenario. With a series of sensitivity tests we quantify additional uncertainties due to the calibration of our model with sparsely observed glacier mass changes. This gives an upper bound for the uncertainty range of ±84 mm sea-level rise by 2100 for each projection.     

Equilibrium ice sheet scaling in climate modeling

A set of simple scaling formulas related to ice sheet evolution is derived from the dynamic and thermodynamic equations for ice and is used to consider two common situations: (a) when we wish to estimate potential ice sheet characteristics given the prescribed net snow accumulation over an area; and (b) when we wish to reconstruct net snow accumulation and vertical temperature difference within the ice sheet given empirical data only concerning ice sheet area and volume. The scaling formulas are applied to the present day Antarctic and Greenland ice sheets, as well as to some ancient ice sheets, and are used to estimate the potential global sea level change due to greenhouse warming.