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.     

Multistability of the Greenland ice sheet and the effects of an adaptive mass balance formulation

The effect of a warmer climate on the Greenland ice sheet as well as its ability to regrow from a reduced geometry is important knowledge when studying future climate. Here we use output from a general circulation model to construct adaptive temperature and precipitation patterns to force an ice flow model off-line taking into consideration that the patterns change in a non-uniform way (both spatially and temporally) as the geometry of the ice sheet evolves and as climate changes. In a series of experiments we investigate the retreat from the present day configuration, build-up from ice free conditions of the ice sheet during a warmer-than-present climate and how the ice sheet moves between states. The adaptive temperature and accumulation patterns as well as two different constant-pattern formulations are applied and all experiments are run to steady state. All results fall into four different groups of geometry regardless of the applied accumulation pattern and initial state. We find that the ice sheet is able to survive and build up at higher temperatures using the more realistic adaptive patterns compared to the classic constant patterns. In contrast, decay occurs at considerably higher temperatures than build-up when the other formulations are used. When studying the motion between states it is clear that the initial state is crucial for the result. The ice sheet is thus multistable at least for certain temperature forcings, and this implies that the ice sheet not does not necessarily return to its initial configuration after a temperature excursion.     

Response of the Antarctic ice sheet to future greenhouse warming

Possible future changes in land ice volume are mentioned frequently as an important aspect of the greenhouse problem. This paper deals with the response of the Antarctic ice sheet and presents a tentative projection of changes in global sea level for the next few hundred years, due to changes in its surface mass balance. We imposed a temperature scenario, in which surface air temperature rises to 4.2° C in the year 2100 AD and is kept constant afterwards. As GCM studies seem to indicate a higher temperature increase in polar latitudes, the response to a more extreme scenario (warming doubled) has also been investigated. The mass balance model, driven by these temperature perturbations, consists of two parts: the accumulation rate is derived from present observed values and is consequently perturbed in proportion to the saturated vapour pressure at the temperature above the inversion layer. The ablation model is based on the degree-day method. It accounts for the daily temperature cycle, uses a different degree-day factor for snow and ice melting and treats refreezing of melt water in a simple way. According to this mass balance model, the amount of accumulation over the entire ice sheet is presently 24.06 × 10¹¹ m³ of ice, and no runoff takes place. A 1°C uniform warming is then calculated to increase the overall mass balance by an amount of 1.43 × 10¹¹ m³ of ice, corresponding to a lowering of global sea level with 0.36 mm/yr. A temperature increase of 5.3°C is needed for the increase in ablation to become more important than the increase in accumulation and the temperature would have to rise by as much as 11.4°C to produce a zero surface mass balance. Imposing the Bellagio-scenario and accumulating changes in mass balance forward in time (static response) would then lower global sea level by 9 cm by 2100 AD. In a subsequent run with a high-resolution 3-D thermomechanic model of the ice sheet, it turns out that the dynamic response of the ice sheet (as compared to the direct effect of the changes in surface mass balance) becomes significant after 100 years or so. Ice-discharge across the grounding-line increases, and eventually leads to grounding-line retreat. This is particularly evident in the extreme case scenario and is important along the Antarctic Peninsula and the overdeepened outlet glaciers along the East Antarctic coast. Grounding-line retreat in the Ross and Ronne-Filchner ice shelves, on the other hand, is small or absent.     

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.     

Climate effect of ozone changes caused by present and future air traffic

 The potential of aircraft-induced ozone changes to force a substantial climate impact is investigated by means of simulations with an atmospheric general circulation model, coupled to a mixed layer ocean model. We present results from several numerical experiments that are based on ozone change patterns for 1992 aviation and on a future scenario for the year 2015. In both cases, the climate signal is statistically significant. The strength of the ozone impact is of comparable magnitude to that arising from aircraft CO₂ emissions, thus meaning a non-negligible contribution to the total climate effect of aviation emissions. There are indications of a characteristic signature of the aircraft ozone related temperature response pattern, distinctly different from that associated with the increase of well-mixed greenhouse gases. Likewise, the climate sensitivity to non-uniform ozone changes including a strong concentration perturbation at the tropopause may be higher than the climate sensitivity to uniform changes of a greenhouse gas. In a hierarchy of experiments, for which the spatial structure of an aircraft-related ozone perturbation was left fixed, while the amplitude of the perturbation was artificially increased, the climate signal depends in a non-linear way on the radiative forcing. Received: 10 September 1998 / Accepted: 4 May 1999     

Verification of model simulated mass balance, flow fields and tabular calving events of the Antarctic ice sheet against remotely sensed observations

The Antarctic ice sheet (AIS) has the greatest potential for global sea level rise. This study simulates AIS ice creeping, sliding, tabular calving, and estimates the total mass balances, using a recently developed, advanced ice dynamics model, known as SEGMENT-Ice. SEGMENT-Ice is written in a spherical Earth coordinate system. Because the AIS contains the South Pole, a projection transfer is performed to displace the pole outside of the simulation domain. The AIS also has complex ice-water-granular material-bedrock configurations, requiring sophisticated lateral and basal boundary conditions. Because of the prevalence of ice shelves, a ‘girder yield’ type calving scheme is activated. The simulations of present surface ice flow velocities compare favorably with InSAR measurements, for various ice-water-bedrock configurations. The estimated ice mass loss rate during 2003–2009 agrees with GRACE measurements and provides more spatial details not represented by the latter. The model estimated calving frequencies of the peripheral ice shelves from 1996 (roughly when the 5-km digital elevation and thickness data for the shelves were collected) to 2009 compare well with archived scatterometer images. SEGMENT-Ice’s unique, non-local systematic calving scheme is found to be relevant for tabular calving. However, the exact timing of calving and of iceberg sizes cannot be simulated accurately at present. A projection of the future mass change of the AIS is made, with SEGMENT-Ice forced by atmospheric conditions from three different coupled general circulation models. The entire AIS is estimated to be losing mass steadily at a rate of ~120 km³/a at present and this rate possibly may double by year 2100.     

Sensitivity of a Greenland ice sheet model to ice flow and ablation parameters: consequences for the evolution through the last climatic cycle

Sensitivity experiments are conducted to test the influence of poorly known model parameters on the simulation of the Greenland ice sheet by means of a three dimensional numerical model including the mechanical and thermal processes within the ice. Two types of experiments are performed: steady-state climatic conditions and simulations over the last climatic cycle with a climatic forcing derived from the GRIP record. The experiments show that the maximum altitude of the ice sheet depends on the ice flow parameters (deformation and sliding law coefficients, geothermal flux) and that it is low when the ice flow is fast. On the other hand, the maximum altitude is not sensitive to the ablation strength and consequently during the climatic cycle it is driven by changes in accumulation rate. The ice sheet extension shows the opposite sensitivity: it is barely affected by ice flow velocity and the ice covered area is smaller for large ablation coefficients. For colder climates, when there is no ablation, the ice sheet extension depends on the sea level. An interesting result is that the variations with time of the altitude at the ice divide (Summit) do not depend on the parameters we tested. The present modelled ice sheets resulting from the climatic cycle experiments are compared with the present measured ice sheet in order to find the set of parameters that gives the best fit between modelled and measured geometry. It seems that, compared to the parameter set most commonly used, higher ablation rate coefficents must be used. Received: 19 September 1995 / Accepted: 30 May 1996     

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.     

Lake ice records used to detect historical and future climatic changes

Historical ice records, such as freeze and breakup dates and the total duration of ice cover, can be used as a quantitative indicator of climatic change if long homogeneous records exist and if the records can be calibrated in terms of climatic changes. Lake Mendota, Wisconsin, has the longest uninterrupted ice records available for any lake in North America dating back to 1855. These records extend back prior to any reliable air temperature data in the midwestern region of the U.S. and demonstrate significant warming of approximately 1.5 °C in fall and early winter temperatures and 2.5 °C in winter and spring temperatures during the past 135 years. These changes are not completely monotonie, but rather appear as two shorter periods of climatic change in the longer record. The first change was between 1875 and 1890, when fall, winter, and spring air temperatures increased by approximately 1.5 °C. The second change, earlier ice breakup dates since 1979, was caused by a significant increase in winter and early spring air temperatures of approximately 1.3 °C. This change may be indicative of shifts in regional climatic patterns associated with global warming, possibly associated with the Greenhouse Effect.With the relationships between air temperature and freeze and break up dates, we can project how the ice cover of Lake Mendota should respond to future climatic changes. If warming occurs, the ice cover for Lake Mendota should decrease approximately 11 days per 1 °C increase. With a warming of 4 to 5 °C, years with no ice cover should occur in approximately 1 out of 15 to 30 years.     

Long-term changes of water balance components in the basins of rivers fed by snow and ice

The variability of the main components of the annual water balance (precipitation, evaporation, glacial alimentation, and dynamic water reserves in the basin) for 1935–1990 is, for the first time, determined for the area where the Zeravshan runoff is formed, higher than hydrological post Dupuli is located. Long-term data on the annual Zeravshan River runoff from an area of 10 200 km² were derived from the measurements at Dupuli hydrological post. The other water balance components were determined with the help of computation methods. Comparison of the measured and calculated volumes of the annual runoff demonstrated that a relative difference between them is systematic, and as a whole for a computation period it is in the interval from ?0.31 to ?4.78%. The annual balance of accumulation and thawing of solid precipitation on glaciers and in the extraglacial area is also determined in the Zeravshan River basin. A new method for computing and mapping spatial variability of the maximum snowline altitude is developed.     

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.     

Greenland Ice Sheet Contribution to Future Global Sea Level Rise based on CMIP5 Models简

Sea level rise （SLR） is one of the major socioeconomic risks associated with global warming. Mass losses from the Greenland ice sheet （GrIS） will be partially responsible for future SLR, although there are large uncertainties in modeled climate and ice sheet behavior. We used the ice sheet model SICOPOLIS （Simulation COde for POLythermal Ice Sheets） driven by climate projections from 20 models in the fifth phase of the Coupled Model Intercomparison Project （CMIP5） to estimate the GrlS contribution to global SLR. Based on the outputs of the 20 models, it is estimated that the GrIS will contribute 0-16 （0-27） cm to global SLR by 2100 under the Representative Concentration Pathways （RCP） 4.5 （RCP 8.5） scenarios. The projected SLR increases further to 7-22 （7-33） cm with 2~basal sliding included. In response to the results of the multimodel ensemble mean, the ice sheet model projects a global SLR of 3 cm and 7 cm （10 cm and 13 cm with 2~basal sliding） under the RCP 4.5 and RCP 8.5 scenarios, respectively. In addition, our results suggest that the uncertainty in future sea level projection caused by the large spread in climate projections could be reduced with model-evaluation and the selective use of model outputs.     

Thresholds for irreversible decline of the Greenland ice sheet

The Greenland ice sheet will decline in volume in a warmer climate. If a sufficiently warm climate is maintained for a few thousand years, the ice sheet will be completely melted. This raises the question of whether the decline would be reversible: would the ice sheet regrow if the climate cooled down? To address this question, we conduct a number of experiments using a climate model and a high-resolution ice-sheet model. The experiments are initialised with ice sheet states obtained from various points during its decline as simulated in a high-CO₂ scenario, and they are then forced with a climate simulated for pre-industrial greenhouse gas concentrations, to determine the possible trajectories of subsequent ice sheet evolution. These trajectories are not the reverse of the trajectory during decline. They converge on three different steady states. The original ice-sheet volume can be regained only if the volume has not fallen below a threshold of irreversibility, which lies between 80 and 90% of the original value. Depending on the degree of warming and the sensitivity of the climate and the ice-sheet, this point of no return could be reached within a few hundred years, sooner than CO₂ and global climate could revert to a pre-industrial state, and in that case global sea level rise of at least 1.3 m would be irreversible. An even larger irreversible change to sea level rise of 5 m may occur if ice sheet volume drops below half of its current size. The set of steady states depends on the CO₂ concentration. Since we expect the results to be quantitatively affected by resolution and other aspects of model formulation, we would encourage similar investigations with other models.     

大气河对南极冰盖及海冰的影响

南极冰盖与海冰对全球气候具有重要影响。大气河作为高纬度地区经向水汽输送的重要途径,其对南极冰盖与海冰的影响在近年来愈发受到重视。南极大气河通常形成于高压脊(阻塞高压)与温带气旋之间的强向极经向输送带内。低频的大气河活动为南极带来强降雪,有利于冰盖质量增加。然而,强暖湿水汽侵入同时会导致表面融化、冰架崩解和极端高温,对冰盖质量存在潜在负贡献。大气河携带极端暖湿水汽与强风通过热力与动力过程导致海冰密集度下降。目前,大气河的识别算法仍不完善,其对液态降水的直接影响、与南大洋的相互作用等仍不清楚,需要进一步明晰大气河对南极冰盖与海冰的影响机制,以准确预估未来大气河对南极冰盖物质平衡与海冰变化的作用。     

The sea ice mass budget of the Arctic and its future change as simulated by coupled climate models

Arctic sea ice mass budgets for the twentieth century and projected changes through the twenty-first century are assessed from 14 coupled global climate models. Large inter-model scatter in contemporary mass budgets is strongly related to variations in absorbed solar radiation, due in large part to differences in the surface albedo simulation. Over the twenty-first century, all models simulate a decrease in ice volume resulting from increased annual net melt (melt minus growth), partially compensated by reduced transport to lower latitudes. Despite this general agreement, the models vary considerably regarding the magnitude of ice volume loss and the relative roles of changing melt and growth in driving it. Projected changes in sea ice mass budgets depend in part on the initial (mid twentieth century) ice conditions; models with thicker initial ice generally exhibit larger volume losses. Pointing to the importance of evolving surface albedo and cloud properties, inter-model scatter in changing net ice melt is significantly related to changes in downwelling longwave and absorbed shortwave radiation. These factors, along with the simulated mean and spatial distribution of ice thickness, contribute to a large inter-model scatter in the projected onset of seasonally ice-free conditions.     

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     

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.