Applying the energy- and water balance model for incorporation of the cryospheric component into a climate model. Part II. Modeled mass balance on the green land ice sheet surface

The Greenland ice sheet is a very important potential source of fresh water inflow to the World Ocean under warming climate conditions. Apparently, it was the same during the Last Interglacial 130-115 thousand years ago. In order to quantify input of the Greenland ice sheet to the rise of the global mean sea level in the past or in the future, we include a surface mass balance model block into the Earth System Model. The computational algorithm is based on the calculation of energy balance on the ice sheet surface. The key tuning parameter of the model is the daily amplitude of air surface temperature. It defines the area and the rate of snow or ice melting. The range of possible values of this parameter is determined during a series of numerical experiments. High sensitivity of meltwater runoff volume to surface air temperature amplitude is revealed.     

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.     

Greenland surface mass balance simulated by a regional climate model and comparison with satellite-derived data in 1990–1991

The 1990 and 1991 ablation seasons over Greenland are simulated with a coupled atmosphere-snow regional climate model with a 25-km horizontal resolution. The simulated snow water content allows a direct comparison with the satellite-derived melt signal. The model is forced with 6-hourly ERA-40 reanalysis at its boundaries. An evaluation of the simulated precipitation and a comparison of the modelled melt zone and the surface albedo with remote sensing observations are presented. Both the distribution and quantity of the simulated precipitation agree with observations from coastal weather stations, estimates from other models and the ERA-40 reanalysis. There are overestimations along the steep eastern coast, which are most likely due to the “topographic barrier effect”. The simulated extent and time evolution of the wet snow zone compare generally well with satellite-derived data, except during rainfall events on the ice sheet and because of a bias in the passive microwave retrieved melt signal. Although satellite-based surface albedo retrieval is only valid in the case of clear sky, the interpolation and the correction of these data enable us to validate the simulated albedo on the scale of the whole Greenland. These two comparisons highlight a large sensitivity of the remote sensing observations to weather conditions. Our high-resolution climate model was used to improve the retrieval algorithms by taking more fully into account the atmosphere variability. Finally, the good agreement of the simulated melting surface with the improved satellite signal allows a detailed estimation of the melting volume from the simulation.     

The 1988–2003 Greenland ice sheet melt extent using passive microwave satellite data and a regional climate model

Measurements from ETH-Camp and JAR1 AWS (West Greenland) as well as coupled atmosphere-snow regional climate simulations have highlighted flaws in the cross-polarized gradient ratio (XPGR) technique used to identify melt from passive microwave satellite data. It was found that dense clouds (causing notably rainfall) on the ice sheet severely perturb the XPGR melt signal. Therefore, the original XPGR melt detection algorithm has been adapted to better incorporate atmospheric variability over the ice sheet and an updated melt trend for the 1988–2003 period has been calculated. Compared to the original algorithm, the melt zone area increase is eight times higher (from 0.2 to 1.7% year⁻¹). The increase is higher with the improved XPGR technique because rainfall also increased during this period. It is correlated to higher atmospheric temperatures. Finally, the model shows that the total ice sheet runoff is directly proportional to the melt extent surface detected by satellites. These results are important for the understanding of the effect of Greenland melting on the stability of the thermohaline circulation.     

Short‐term volume changes of the Greenland ice sheet in response to doubled CO2 conditions

This paper focuses on the rôle of accumulation and cloudiness changes in the response of the Greenland ice sheet to global warming. Changes in accumulation or cloudiness were often neglected, or coupled to temperature changes. We used model output on temperature, precipitation and cloudiness from a GCM (ECHAM4 T106). The GCM output was used to drive the Greenland model that exists of a vertically averaged ice flow model, coupled to a 1D surface energy balance model that calculates the ablation. Variables are temperature, accumulation and cloudiness. Sensitivity experiments with this model show that changes in accumulation are very important for the ice sheet mass balance, whereas cloudiness is of secondary importance. If the Greenland model is forced by the GCM output, the Greenland model is found to contribute 70% less to sea level rise after 70 years than is indicated by the results presented in the IPCC report. This large discrepancy is mainly due to the fact that the enhanced ablation is strongly compensated by increased accumulation. Comparing the result obtained here with changes in mass balance derived directly from the same general circulation model, indicates a 20% larger contribution to sea level. This increase is due to changes in ice flow, and a different method for the ablation calculation.     

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.     

Temporal and Spatial Variations of the Aerodynamic Roughness Length in the Ablation Zone of the Greenland Ice Sheet

To understand the response of the Greenland ice sheet to climate change the so-called ablation zone is of particular importance, since it accommodates the yearly net surface ice loss. In numerical models and for data analysis, the bulk aerodynamic method is often used to calculate the turbulent surface fluxes, for which the aerodynamic roughness length (z ₀) is a key parameter. We present, for the first time, spatial and temporal variations of z ₀ in the ablation area of the Greenland ice sheet using year-round data from three automatic weather stations and one eddy-correlation mast. The temporal variation of z ₀ is found to be very high in the lower ablation area (factor 500) with, at the end of the summer melt, a maximum in spatial variation for the whole ablation area of a factor 1000. The variation in time matches the onset of the accumulation and ablation season as recovered by sonic height rangers. During winter, snow accumulation and redistribution by snow drift lead to a uniform value of z ₀≈ 10⁻⁴ m throughout the ablation area. At the beginning of summer, snow melt uncovers ice hummocks and z ₀ quickly increases well above 10⁻² m in the lower ablation area. At the end of summer melt, hummocky ice dominates the surface with z ₀ > 5  ×  10⁻³ m up to 60 km from the ice edge. At the same time, the area close to the equilibrium line (about 90 km from the ice edge) remains very smooth with z ₀ = 10⁻⁵ m. At the beginning of winter, we observed that single snow events have the potential to lower z ₀ for a very rough ice surface by a factor of 20 to 50. The total surface drag of the abundant small-scale ice hummocks apparently dominates over the less frequent large domes and deep gullies. The latter results are verified by studying the individual drag contributions of hummocks and domes with a drag partition model.     

Greenland��s contribution to global sea-level rise by the end of the 21st century

The Greenland ice sheet holds enough water to raise the global sea level with ??7 m. Over the last few decades, observations manifest a substantial increase of the mass loss of this ice sheet. Both enhanced melting and increase of the dynamical discharge, associated with calving at the outlet-glacier fronts, are contributing to the mass imbalance. Using a dynamical and thermodynamical ice-sheet model, and taking into account speed up of outlet glaciers, we estimate Greenland??s contribution to the 21st-century global sea-level rise and the uncertainty of this estimate. Boundary fields of temperature and precipitation extracted from coupled climate-model projections used for the IPCC Fourth Assessment Report, are applied to the ice-sheet model. We implement a simple parameterization for increased flow of outlet glaciers, which decreases the bias of the modeled present-day surface height. It also allows for taking into account the observed recent increase in dynamical discharge, and it can be used for future projections associated with outlet-glacier speed up. Greenland contributes 0?C17?cm to global sea-level rise by the end of the 21st century. This range includes the uncertainties in climate-model projections, the uncertainty associated with scenarios of greenhouse-gas emissions, as well as the uncertainties in future outlet-glacier discharge. In addition, the range takes into account the uncertainty of the ice-sheet model and its boundary fields.     

Late Pliocene to Pleistocene sensitivity of the Greenland Ice Sheet in response to external forcing and internal feedbacks

The timing and nature of ice sheet variations on Greenland over the last ～5 million years remain largely uncertain. Here, we use a coupled climate-vegetation-ice sheet model to determine the climatic sensitivity of Greenland to combined sets of external forcings and internal feedbacks operating on glacial-interglacial timescales. In particular, we assess the role of atmospheric pCO₂, orbital forcing, and vegetation dynamics in modifying thresholds for the onset of glaciation in late Pliocene and Pleistocene. The response of circum-Arctic vegetation to declining levels of pCO₂ (from 400 to 200 ppmv) and decreasing summer insolation includes a shift from boreal forest to tundra biomes, with implications for the surface energy balance. The expansion of tundra amplifies summer surface cooling and heat loss from the ground, leading to an expanded summer snow cover over Greenland. Atmospheric and land surface fields respond to forcing most prominently in late spring-summer and are more sensitive at lower Pleistocene-like levels of pCO₂. We find cold boreal summer orbits produce favorable conditions for ice sheet growth, however simulated ice sheet extents are highly dependent on both background pCO₂ levels and land-surface characteristics. As a result, late Pliocene ice sheet configurations on Greenland differ considerably from late Pleistocene, with smaller ice caps on high elevations of southern and eastern Greenland, even when orbital forcing is favorable for ice sheet growth.     

Climate modification by future ice sheet changes and consequences for ice sheet mass balance

The future evolution of global ice sheets under anthropogenic greenhouse forcing and its impact on the climate system, including the regional climate of the ice sheets, are investigated with a comprehensive earth system model consisting of a coupled Atmosphere–Ocean General Circulation Model, a dynamic vegetation model and an ice sheet model. The simulated control climate is realistic enough to permit a direct coupling of the atmosphere and ice sheet components, avoiding the use of anomaly coupling, which represents a strong improvement with respect to previous modelling studies. Glacier ablation is calculated with an energy-balance scheme, a more physical approach than the commonly used degree-day method. Modifications of glacier mask, topographic height and freshwater fluxes by the ice sheets influence the atmosphere and ocean via dynamical and thermodynamical processes. Several simulations under idealized scenarios of greenhouse forcing have been performed, where the atmospheric carbon dioxide stabilizes at two and four times pre-industrial levels. The evolution of the climate system and the ice sheets in the simulations with interactive ice sheets is compared with the simulations with passively coupled ice sheets. For a four-times CO₂ scenario forcing, a faster decay rate of the Greenland ice sheet is found in the non-interactive case, where melting rates are higher. This is caused by overestimation of the increase in near-surface temperature that follows the reduction in topographic height. In areas close to retreating margins, melting rates are stronger in the interactive case, due to changes in local albedo. Our results call for careful consideration of the feedbacks operating between ice sheets and climate after substantial decay of the ice sheets.     

Present and future near-surface wind climate of Greenland from high resolution regional climate modelling

The present and twenty-first century near-surface wind climate of Greenland is presented using output from the regional atmospheric climate model RACMO2. The modelled wind variability and wind distribution compare favourably to observations from three automatic weather stations in the ablation zone of southwest Greenland. The Weibull shape parameter is used to classify the wind climate. High values (κ > 4) are found in northern Greenland, indicative of uniform winds and a dominant katabatic forcing, while lower values (κ < 3) are found over the ocean and southern Greenland, where the synoptic forcing dominates. Very high values of the shape parameter are found over concave topography where confluence strengthens the katabatic circulation, while very low values are found in a narrow band along the coast due to barrier winds. To simulate the future (2081–2098) wind climate RACMO2 was forced with the HadGEM2-ES general circulation model using a scenario of mid-range radiative forcing of +4.5 W m^?2 by 2100. For the future simulated climate, the near-surface potential temperature deficit reduces in all seasons in regions where the surface temperature is below the freezing point, indicating a reduction in strength of the near-surface temperature inversion layer. This leads to a wind speed reduction over the central ice sheet where katabatic forcing dominates, and a wind speed increase over steep coastal topography due to counteracting effects of thermal and katabatic forcing. Thermally forced winds over the seasonally sea ice covered region of the Greenland Sea are reduced by up to 2.5 m s^?1.     

Impact of Greenland and Antarctic ice sheet interactions on climate sensitivity

We use the Earth system model of intermediate complexity LOVECLIM to show the effect of coupling interactive ice sheets on the climate sensitivity of the model on a millennial time scale. We compare the response to a 2×CO₂ warming scenario between fully coupled model versions including interactive Greenland and Antarctic ice sheet models and model versions with fixed ice sheets. For this purpose an ensemble of different parameter sets have been defined for LOVECLIM, covering a wide range of the model??s sensitivity to greenhouse warming, while still simulating the present-day climate and the climate evolution over the last millennium within observational uncertainties. Additional freshwater fluxes from the melting ice sheets have a mitigating effect on the model??s temperature response, leading to generally lower climate sensitivities of the fully coupled model versions. The mitigation is effectuated by changes in heat exchange within the ocean and at the sea?Cair interface, driven by freshening of the surface ocean and amplified by sea?Cice-related feedbacks. The strength of the effect depends on the response of the ice sheets to the warming and on the model??s climate sensitivity itself. The effect is relatively strong in model versions with higher climate sensitivity due to the relatively large polar amplification of LOVECLIM. With the ensemble approach in this study we cover a wide range of possible model responses.     

The effect of snow on Antarctic sea ice simulations in a coupled atmosphere-sea ice model

 The effect of a snow cover on sea ice accretion and ablation is estimated based on the ‘zero-layer’ version sea ice model of Semtner, and is examined using a coupled atmosphere-sea ice model including feedbacks and ice dynamics effects. When snow is disregarded in the coupled model the averaged Antarctic sea ice becomes thicker. When only half of the snowfall predicted by the atmospheric model is allowed to land on the ice surface sea ice gets thicker in most of the Weddell and Ross Seas but thinner in East Antarctic in winter, with the average slightly thicker. When twice as much snowfall as predicted by the atmospheric model is assumed to land on the ice surface sea ice also gets much thicker due to the large increase of snow-ice formation. These results indicate the importance of the correct simulation of the snow cover over sea ice and snow-ice formation in the Antarctic. Our results also illustrate the complex feedback effects of the snow cover in global climate models. In this study we have also tested the use of a mean value of 0.16 Wm^-1 K^-1 instead of 0.31 for the thermal conductivity of snow in the coupled model, based on the most recent observations in the eastern Antarctic and Bellingshausen and Amundsen Seas, and have found that the sea ice distribution changes greatly, with the ice becoming much thinner by about 0.2 m in the Antarctic and about 0.4 m in the Arctic on average. This implies that the magnitude of the thermal conductivity of snow is of considerable importance for the simulation of the sea ice distribution. An appropriate value of the thermal conductivity of snow is as crucial as the depth of the snow layer and the snowfall rate in a sea ice model. The coupled climate models require accurate values of the effective thermal conductivity of snow from observations for validating the simulated sea ice distribution under the present climate conditions. Received: 20 November 1997/Accepted: 27 July 1998     

Global Warming and the Greenland Ice Sheet

The Greenland coastal temperatures have followed the early 20th century global warming trend. Since 1940, however, the Greenland coastal stations data have undergone predominantly a cooling trend. At the summit of the Greenland ice sheet the summer average temperature has decreased at the rate of 2.2 °C per decade since the beginning of the measurements in 1987. This suggests that the Greenland ice sheet and coastal regions are not following the current global warming trend. A considerable and rapid warming over all of coastal Greenland occurred in the 1920s when the average annual surface air temperature rose between 2 and 4 °C in less than ten years (at some stations the increase in winter temperature was as high as 6 °C). This rapid warming, at a time when the change in anthropogenic production of greenhouse gases was well below the current level, suggests a high natural variability in the regional climate. High anticorrelations (r = ?0.84 to?0.93) between the NAO (North Atlantic Oscillation) index and Greenland temperature time series suggest a physical connection between these processes. Therefore, the future changes in the NAO and Northern Annular Mode may be of critical consequence to the future temperature forcing of the Greenland ice sheet melt rates.     

Modelling the response of glaciers to a doubling in atmospheric CO2: a case study of Storglaciären, northern Sweden

To predict the evolution of glaciers in an enhanced greenhouse climate, results from a global climate model, a glacier melt/accumulation model, and a glacier flow model were combined. The method was applied to Storglaciären, a small well-studied glacier in northern Sweden. The difference between the present climate and a 2 × CO₂ climate around the year 2050 was extracted from a model experiment with the ECHAM4-T106 high resolution climate model for time slices at present and in 2050, using prescribed boundary conditions of sea surface temperature and sea-ice distribution, which are derived from a lower resolution transient run of the ECHAM4-T42/OPIC-coupled atmosphere ocean model between present and 2050. The local climatic conditions on the glacier for 2050 were obtained by adding the modelled local climate changes to the observed local present-day climate. The combination of the comprehensive models presented offers a tool to test and calibrate simplified models which are applicable to a much larger sample of glaciers. For the region of Storglaciären, the GCM projected temperature is found to increase most strongly during the winter months, but also shows a warming during the transition from spring to summer, and again between summer and fall, thus extending the melt season by three to four weeks. Precipitation, on the other hand, decreases by approximately 5% during May to September while there is a stronger increase of approximately 14% for the rest of the year. The consequent increase in winter accumulation on Storglaciären is more than compensated by the increase in ablation during the melt season. The glacier flow model predicts a 300 m retreat of the glacier terminus by the middle of the next century, and a loss of 30% of the present ice mass.     

On the Response of the Greenland Ice Sheet to Greenhouse Climate Change

Numerical computations are performed with the three-dimensional polythermal ice-sheet model SICOPOLIS in order to investigate the possible impact of a greenhouse-gas-induced climate change on the Greenland ice sheet. The assumed increase of the mean annual air temperature above the ice covers a range from T = 1°C to 12°C, and several parameterizations for the snowfall and the surface melting are considered. The simulated shrinking of the ice sheet is a smooth function of the temperature rise, indications for the existence of critical thresholds of the climate input are not found. Within 1000 model years, the ice-volume decrease is limited to 10% of the present volume for T 3°C, whereas the most extreme scenario, T = 12°C, leads to an almost entire disintegration, which corresponds to a sea-level equivalent of 7 m. The different snowfall and melting parameterizations yield an uncertainty range of up to 20% of the present ice volume after 1000 model years.     

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     

Future surface mass balance of the Antarctic ice sheet and its influence on sea level change, simulated by a regional atmospheric climate model

A regional atmospheric climate model with multi-layer snow module (RACMO2) is forced at the lateral boundaries by global climate model (GCM) data to assess the future climate and surface mass balance (SMB) of the Antarctic ice sheet (AIS). Two different GCMs (ECHAM5 until 2100 and HadCM3 until 2200) and two different emission scenarios (A1B and E1) are used as forcing to capture a realistic range in future climate states. Simulated ice sheet averaged 2 m air temperature (T_2m) increases (1.8–3.0 K in 2100 and 2.4–5.3 K in 2200), simultaneously and with the same magnitude as GCM simulated T_2m. The SMB and its components increase in magnitude, as they are directly influenced by the temperature increase. Changes in atmospheric circulation around Antarctica play a minor role in future SMB changes. During the next two centuries, the projected increase in liquid water flux from rainfall and snowmelt, together 60–200 Gt year^?1, will mostly refreeze in the snow pack, so runoff remains small (10–40 Gt year^?1). Sublimation increases by 25–50 %, but remains an order of magnitude smaller than snowfall. The increase in snowfall mainly determines future changes in SMB on the AIS: 6–16 % in 2100 and 8–25 % in 2200. Without any ice dynamical response, this would result in an eustatic sea level drop of 20–43 mm in 2100 and 73–163 mm in 2200, compared to the twentieth century. Averaged over the AIS, a strong relation between $\Updelta$ SMB and $\Updelta\hbox{T}_{2{\rm m}}$ of 98 ± 5 Gt w.e. year^?1 K^?1 is found.