Large-Scale Surface Mass Balance of Ice Sheets from a Comprehensive Atmospheric Model

The surface mass balance for Greenland and Antarctica has been calculated using model data from an AMIP-type experiment for the period 1979?C2001 using the ECHAM5 spectral transform model at different triangular truncations. There is a significant reduction in the calculated ablation for the highest model resolution, T319 with an equivalent grid distance of ca 40?km. As a consequence the T319 model has a positive surface mass balance for both ice sheets during the period. For Greenland, the models at lower resolution, T106 and T63, on the other hand, have a much stronger ablation leading to a negative surface mass balance. Calculations have also been undertaken for a climate change experiment using the IPCC scenario A1B, with a T213 resolution (corresponding to a grid distance of some 60?km) and comparing two 30-year periods from the end of the twentieth century and the end of the twenty-first century, respectively. For Greenland there is change of 495?km³/year, going from a positive to a negative surface mass balance corresponding to a sea level rise of 1.4?mm/year. For Antarctica there is an increase in the positive surface mass balance of 285?km³/year corresponding to a sea level fall by 0.8?mm/year. The surface mass balance changes of the two ice sheets lead to a sea level rise of 7?cm at the end of this century compared to end of the twentieth century. Other possible mass losses such as due to changes in the calving of icebergs are not considered. It appears that such changes must increase significantly, and several times more than the surface mass balance changes, if the ice sheets are to make a major contribution to sea level rise this century. The model calculations indicate large inter-annual variations in all relevant parameters making it impossible to identify robust trends from the examined periods at the end of the twentieth century. The calculated inter-annual variations are similar in magnitude to observations. The 30-year trend in SMB at the end of the twenty-first century is significant. The increase in precipitation on the ice sheets follows closely the Clausius-Clapeyron relation and is the main reason for the increase in the surface mass balance of Antarctica. On Greenland precipitation in the form of snow is gradually starting to decrease and cannot compensate for the increase in ablation. Another factor is the proportionally higher temperature increase on Greenland leading to a larger ablation. It follows that a modest increase in temperature will not be sufficient to compensate for the increase in accumulation, but this will change when temperature increases go beyond any critical limit. Calculations show that such a limit for Greenland might well be passed during this century. For Antarctica this will take much longer and probably well into following centuries.     

冰盖消融的海平面指纹变化及其对GRACE监测结果的影响

全球变暖背景下的冰盖消融以及由此带来海平面上升日益明显,直接影响地球表面的陆地水质量平衡,以及固体地球瞬间弹性响应,研究冰盖质量变化的海平面指纹能够帮助深入了解未来海平面区域变化的驱动因素.本文基于海平面变化方程并考虑负荷自吸效应（SAL）与地球极移反馈的影响,借助美国德克萨斯大学空间研究中心（Center for Space Research,CSR）发布的2003年到2012年十年期间的GRACE重力场月模型数据（RL05）,结合加权高斯平滑的区域核函数,反演得到格陵兰与南极地区冰盖质量变化的时空分布,并利用海平面变化方程计算得到了相对海平面的空间变化,结果表明：格陵兰与南极冰盖质量整体呈明显的消融趋势,变化速率分别为-273.31 Gt/a及-155.56 Gt/a,由此导致整个北极圈相对海平面降低,最高可达约-0.6 cm·a^-1;而南极地区冰盖质量变化趋势分布不一,导致西南极近海相对海平面下降,而东南极地区近海相对海平面上升,最高可达约0.2 cm·a^-1.远离质量负荷区域的全球海平面以上升趋势为主,平均全球相对海平面上升0.71 mm·a^-1,部分远海地区相对海平面上升更加突出（例如北美与澳大利亚）,高出全球平均海平面上升速率将近30%.此外,本文也重点探讨了GRACE监测冰盖消融结果中由于极地近海海平面变化导致的泄漏影响,经此项影响校正后的结果表明：海平面指纹效应对GRACE监测格陵兰与南极地区2003-2012期间整体冰盖消融速率的贡献分别为约3%与9%,建议在后期利用GRACE更精确地估算研究区冰盖质量变化时,应考虑海平面指纹效应的渗透影响.     

Greenland and Antarctica Ice Sheet Mass Changes and Effects on Global Sea Level

Thirteen years of GRACE data provide an excellent picture of the current mass changes of Greenland and Antarctica, with mass loss in the GRACE period 2002–2015 amounting to 265 ± 25 GT/year for Greenland (including peripheral ice caps), and 95 ± 50 GT/year for Antarctica, corresponding to 0.72 and 0.26 mm/year average global sea level change. A significant acceleration in mass loss rate is found, especially for Antarctica, while Greenland mass loss, after a corresponding acceleration period, and a record mass loss in the summer of 2012, has seen a slight decrease in short-term mass loss trend. The yearly mass balance estimates, based on point mass inversion methods, have relatively large errors, both due to uncertainties in the glacial isostatic adjustment processes, especially for Antarctica, leakage from unmodelled ocean mass changes, and (for Greenland) difficulties in separating mass signals from the Greenland ice sheet and the adjacent Canadian ice caps. The limited resolution of GRACE affects the uncertainty of total mass loss to a smaller degree; we illustrate the “real” sources of mass changes by including satellite altimetry elevation change results in a joint inversion with GRACE, showing that mass change occurs primarily associated with major outlet glaciers, as well as a narrow coastal band. For Antarctica, the primary changes are associated with the major outlet glaciers in West Antarctica (Pine Island and Thwaites Glacier systems), as well as on the Antarctic Peninsula, where major glacier accelerations have been observed after the 2002 collapse of the Larsen B Ice Shelf.     

MASS BALANCE OF POLAR ICE FROM LONG WAVELENGTH FEATURES OF THE EARTH'S GRAVITATIONAL FIELD

We have used satellite solutions to the low degree zonal harmonics of the Earth's gravitational potential, and rates of surface accumulation to partially constrain, by means of repeated forward solution, the time rates of thickness change over the Antarctic and Greenland Ice Sheets (dTA and dTG respectively). In addition to the observed zonal coefficients j2 through j5 we impose only one other constraint: That dTA and dTG are proportional to surface accumulation. The lagged response of the Earth to secular changes in ice thickness spanning recent time periods (up to 2000 years before present) and the late Pleistocene is accounted for by means of two viscoelastic rebound models. The sea level contributions from the ice sheets, calculated from dTA and dTG, lower mantle viscosity, and the start time of present-day thickness change are all variables subject to the constraints. For a given set of post glacial rebound inputs, a family of solutions that have similar characteristics and that agree well with observation are obtained from the large number of forward solutions. The off axis position of the Greenland ice sheet makes its contribution to the low degree zonal coefficients less sensitive to the spatial details of the mass balance than to the overall sea level contribution. dTG is therefore modeled as surface mass balance offset by a uniform and constant mass loss. Though dTA varies widely with choices of input parameters, the combined sea level contribution from both ice sheets is reasonably well constrained by the gravity coefficients, and is predicted to range from -0.9 to +1.6 mm yr-1. The sign of the slope of the low degree zonal coefficients versus sea level contribution for Greenland is positive, but for Antarctica, the sign of the slope is positive for even degree and negative for odd degree harmonics. By using this property of the zonal coefficients, it is possible to determine the individual sea level contributions for Greenland and Antarctica. They vary from -0.6 to +0.3 mm yr-1 for the Greenland Ice Sheet, and from -0.3 to +1.3 mm yr-1 for the Antarctic Ice Sheet.     

基于GRACE/GRACE-FO和Swarm卫星研究2002-2020年南极和格陵兰岛冰盖质量时空变化

两极冰盖消融及其质量变化作为全球气候变化的重要指标之一,一直是联合国政府间专门气候委员会IPCC(Intergovernmental Panel on Climate Change)报告的重点关注内容.GRACE(Gravity Recovery and Climate Experiment,2002年4月-2017年...     

Impact of precipitation seasonality changes on isotopic signals in polar ice cores: a multi-model analysis

For Central Greenland, water isotope analysis indicates a temperature difference of about 10°C since the Last Glacial Maximum (LGM). However, borehole thermometry and gas diffusion thermometry indicate that LGM surface temperatures were about 20°C colder than today. Two general circulation model studies have shown that changes in the seasonal precipitation timing in Central Greenland might have caused a warm bias in the LGM water isotope proxy temperatures, and that this bias could explain the difference in the estimated paleotemperatures. Here we present an analysis of a number of atmospheric general circulation model simulations mostly done within the framework of the Paleoclimate Modeling Intercomparison Project. The models suggest that the seasonal cycle of precipitation and surface mass balance over Central Greenland at the LGM might have been very different from today. This supports the idea that the accuracy of the water isotope thermometry at the LGM in Greenland might be compromised as a result of a modified surface mass balance seasonality. However, the models disagree on the amplitude and sign of the bias. For Central East Antarctica, a strong seasonality effect on the LGM isotopic signal is not simulated by any of the analyzed models. For the mid-Holocene (6 kyr BP) the models suggest relatively weak isotope paleothermometry biases linked to changes in the surface mass balance seasonality over both ice sheets.     

ANNUAL SEA LEVEL VARIABILITY INDUCED BY CHANGES IN SEA ICE EXTENT AND ACCUMULATION ON ICE SHEETS: AN ASSESSMENT BASED ON REMOTELY SENSED DATA

Changes of mean annual net accumulation at the surface on the grounded ice sheets of East Antarctica, West Antarctica and Greenland in response to variations in sea ice extent are estimated using grid-point values 100 km apart. The data bases are assembled principally by bilinear interpolation of remotely sensed brightness temperature (Nimbus-5 ESMR, Nimbus-7 SMMR), surface temperature (Nimbus-7 THIR), and surface elevation (ERS-1 radar altimeter). These data, complemented by field data where remotely sensed data are not available, are used in multivariate analyses in which mean annual accumulation (derived from firn emissivity) is the dependent variable; the independent variables are latitude, surface elevation, mean annual surface temperature, and mean annual distance to open ocean (as a source of energy and moisture). The last is the shortest distance measured between a grid point and the mean annual position of the 10% sea ice concentration boundary, and is used as an index of changes in sea ice extent as well as of mean concentration. Stepwise correlation analyses indicate that variations in sea ice extent of ± 50 km would lead to changes in accumulation inversely of ± 4% on East Antarctica, ± 10% on West Antarctica, and ±4% on Greenland. These results are compared with those obtained in a previous study using visually interpolated values from contoured compilations of field data; they substantiate the findings for the Antarctic ice sheets (±4% on East Antarctica, ±9% in West Antarctica), and suggest a reduction by one half of the probable change of accumulation on Greenland (from ±8%). The results also suggest a reduction of the combined contribution to sea level variability to ±0.19 mm a-1 (from ±0.22 mm a-1).     

Effect of ice sheet growth and melting on the slip evolution of thrust faults

Field investigations suggest that postglacial unloading and rebound led to the formation or re-activation of reverse faults even in continental shields like Scandinavia. Here we use finite-element models including a thrust fault embedded in a rheologically layered lithosphere to investigate its slip evolution during glacial loading and subsequent postglacial unloading. The model results show that the rate of thrusting decreases during the presence of an ice sheet and strongly increases during deglaciation. The magnitude of the slip acceleration is primarily controlled by the thickness of the ice sheet, the viscosity of the lithospheric layers and the long-term shortening rate. In contrast, the width of the ice sheet, the rate of deglaciation or the fault dip have an only minor influence on the slip evolution. In all experiments, the slip rate variations are caused by changes in the differential stress. The modelled deglaciation-induced slip acceleration agrees well with the occurrence of large earthquakes soon after the melting of the Fennoscandian ice sheet, which led to the formation of spectacular fault scarps in particular in the Lapland Fault Province. Furthermore, our model results support the idea that the low level of seismicity in currently glaciated regions like Greenland and Antarctica is caused by the presence of the ice sheets. Based on our models we expect that the decay of the Greenland and Antarctica ice sheets in the course of global warming will ultimately lead to an increase in earthquake frequency in these regions.     

Improved Constraints on Models of Glacial Isostatic Adjustment: A Review of the Contribution of Ground-Based Geodetic Observations

The provision of accurate models of Glacial Isostatic Adjustment (GIA) is presently a priority need in climate studies, largely due to the potential of the Gravity Recovery and Climate Experiment (GRACE) data to be used to determine accurate and continent-wide assessments of ice mass change and hydrology. However, modelled GIA is uncertain due to insufficient constraints on our knowledge of past glacial changes and to large simplifications in the underlying Earth models. Consequently, we show differences between models that exceed several mm/year in terms of surface displacement for the two major ice sheets: Greenland and Antarctica. Geodetic measurements of surface displacement offer the potential for new constraints to be made on GIA models, especially when they are used to improve structural features of the Earth’s interior as to allow for a more realistic reconstruction of the glaciation history. We present the distribution of presently available campaign and continuous geodetic measurements in Greenland and Antarctica and summarise surface velocities published to date, showing substantial disagreement between techniques and GIA models alike. We review the current state-of-the-art in ground-based geodesy (GPS, VLBI, DORIS, SLR) in determining accurate and precise surface velocities. In particular, we focus on known areas of need in GPS observation level models and the terrestrial reference frame in order to advance geodetic observation precision/accuracy toward 0.1 mm/year and therefore further constrain models of GIA and subsequent present-day ice mass change estimates.     

Measurement of ice flow velocities from GPS positions logged by short-period seismographs in East Antarctica

The ice flow velocity is a basic feature of glaciers and ice sheets. Measuring ice flow velocities is very important for estimating the mass balance of ice sheets in the Arctic and Antarctic. Traditional methods for measuring ice flow velocity include the use of stakes, snow pits and on-site geodetic GPS and remote sensing measurement methods. Geodetic GPS measurements have high accuracy, but geodetic GPS monitoring points only sparsely cover the Antarctic ice sheets. Moreover, the resolution and accuracy of ice flow velocities based on remote sensing measurements are low. Although the accuracy of the location data recorded by the navigation-grade GPS receivers embedded in short-period seismographs is not as good as that of geodetic GPS,the ice flow velocity can be accurately measured by these navigation-grade GPS data collected over a sufficiently long period. In this paper, navigation-grade GPS location data obtained by passive seismic observations during the 36 th Chinese National Antarctic Research Expedition were used to accurately track the movement characteristics of the ice sheet in the Larsemann Hills of East Antarctica and the Taishan Station area. The results showed that the ice sheet in the two study areas is basically moving northwestward with an average ice flow velocity of approximately 1 m mon-1. The results in the Taishan Station area are basically consistent with the geodetic GPS results, indicating that it is feasible to use the embedded GPS location data from shortperiod seismographs to track the movement characteristics of ice sheets. The ice flow characteristics in the Larsemann Hills are more complex. The measured ice flow velocities in the Larsemann Hills with a resolution of 200 m help to understand its characteristics. In summary, the ice flow velocities derived from GPS location data are of great significance for studying ice sheet dynamics and glacier mass balance and for evaluating the systematic errors caused by ice sheet movements in seismic imaging.     

Measured and modelled absolute gravity changes in Greenland

In glaciated areas, the Earth is responding to the ongoing changes of the ice sheets, a response known as glacial isostatic adjustment (GIA). GIA can be investigated through observations of gravity change. For the ongoing assessment of the ice sheets mass balance, where satellite data are used, the study of GIA is important since it acts as an error source. GIA consists of three signals as seen by a gravimeter on the surface of the Earth. These signals are investigated in this study. The ICE-5G ice history and recently developed ice models of present day changes are used to model the gravity change in Greenland. The result is compared with the initial measurements of absolute gravity (AG) change at selected Greenland Network (GNET) sites.We find that observations are highly influenced by the direct attraction from the ice and ocean. This is especially evident in the measurements conducted at the GNET station near the Helheim Glacier. The effect of the direct attraction diminishes at sites that are more than one degree from the source. Here, the dominant signal is the effect of the elastic signal from present day ice mass changes. We find agreement between the measured and modelled gravity changes at all but one site. This agreement only holds when the direct attraction is considered. For one site, there is no agreement, indicating that some improvements to the modelling results or the processing of the gravity data are needed. In addition, more AG measurements are needed to strengthen the time series of gravity change.     

Ice sheets and the CO2 problem

In this review, the carbon dioxide problem is discussed, with special reference to the possible effects of a global warming on the ice sheets of Greenland and Antarctica. Instead of detailed projections of future climate and the consequences, the basic mechanisms are explained and illustrated with results described in the literature.It is concluded that a doubling of the atmospheric CO₂ content (most likely to occur somewhere in the second half of the next century) will result in a globally-averaged warming of 2–4°C, and an intensification of the hydrological cycle. In the polar regions, this warming will be a few degrees larger and as a consequence the Greenland Ice Sheet will decrease in size. Antarctica, on the other hand, is expected to grow because of the increased snowfall. The instability of the West Antarctic Ice Sheet is also discussed and, although no conclusive prediction to its long-term response can be made, it is argued that on a short time scale (less than about 100 y) nothing dramatically wil happen to this part of Antarctica.     

冰川均衡调整对南极冰质量平衡监测的影响及其不确定性

本文研究了新的全球冰川均衡调整(GIA)模型对南极冰盖质量平衡监测的影响,考虑现有冰川负荷模型和地幔黏滞度模型的差异,完整评估了结果的不确定性,最后结合GRACE和卫星测高的结果进行了对比分析.结果表明,GIA对GRACE监测的等效水柱变化有重大影响,较大的GIA影响出现在西南极,沿罗斯冰架-卡姆布冰流-罗尼冰架-南极...     

Present Day Regional Mass Loss of Greenland Observed with Satellite Gravimetry

This paper summarizes results obtained for Greenland??s mass balance observed with NASA??s GRACE mission. We estimate a Greenland ice sheet mass loss at ?201 ± 19 Gt/year including a discernible acceleration of ?8 ± 7 Gt/year² between March 2003 and February 2010. The mass loss of glacier systems on the South East of Greenland has slowed down while the mass loss increases toward the North along the West side of Greenland. The mass balance can be compared with results obtained by a regional climate model of the Greenland system and ice sheet altimeter data obtained from NASA??s ICEsat mission. Our GRACE-only results differ to within 15% from these independently calculated values; we will comment on the possible causes and the quality of the glacial isostatic adjustment model which is used to correct geodetic datasets.     

基于GRACE RL06数据监测和分析南极冰盖27个流域质量变化

本文基于CSR最新公布的GRACE RL06版本数据,采用Slepian空域反演法估算了南极冰盖27个流域的质量变化.Slepian空域反演法结合了Slepian空间谱集中法和空域反演法的技术优势,能够有效降低GRACE在小区域反演时信号出现的严重泄漏和衰减,进而精确获得南极冰盖在每个流域的质量变化.相对于GRACE RL05版本数据,RL06在条带误差的控制上要更加优化,获得的南极冰盖质量变化时间序列也更加平滑,但在趋势估算上差别并不明显(小于10Gt/a).本文的估算结果显示:在2002年4月至2016年8月期间,整个南极冰盖质量变化速率为-118.6±16.3Gt/a,其中西南极为-142.4±10.5Gt/a,南极半岛为-29.2±2.1Gt/a,东南极则为52.9±8.6Gt/a.南极冰盖损失最大的区域集中在西南极Amundsen Sea Embayment(流域20-23),该地区质量变化速率为-203.5±4.1Gt/a,其次为南极半岛(流域24-27)以及东南极Victoria-Wilkes Land (流域13-15),质量变化速率分别为-29.2±2.1Gt/a和-19.0±4.7Gt/a,其中Amundsen Sea Embayment和南极半岛南部两个地区的冰排放呈现加速状态.南极冰盖质量显著增加的区域主要有西南极的Ellsworth Land(流域1)和Siple Coast(流域18和19)以及东南极的Coats-Queen Maud-Enderby Land (流域3-8),三个地区质量变化速率分别为17.2±2.4Gt/a、43.9±1.9Gt/a和62.7±3.8Gt/a,质量增加大多来自降雪累积,比如:Coats-Queen Maud-Enderby Land在2009年和2011年发生的大规模降雪事件,但也有来自冰川的增厚,如:Siple Coast地区Kamb冰流的持续加厚.此外,对GRACE估算的南极冰盖质量变化年际信号进行初步分析发现,GRACE年际信号与气候模型估算的冰盖表面质量平衡年际信号存在显著的线性相关关系,但与主要影响南极气候年际变化的气候事件之间却不存在线性相关关系,这说明南极冰盖质量变化的年际信号主要受冰盖表面质量平衡的支配,而气候事件对冰盖表面质量平衡的影响可能是复杂的非线性耦合过程.     

利用ICESat数据解算南极冰盖冰雪质量变化

南极冰盖冰雪质量变化反映了全球气候变化,并且直接影响着全球海平面变化.ICESat测高卫星的主要任务之一就是要确定南北两极冰盖的质量变化情况并评估其对全球海平面变化的影响.本文利用2003年10月至2008年12月的ICESat测高数据,针对南极DEM分辨率有限的特殊性,通过求解坡度改正值,解决重复轨道地面脚点不重合的问题,计算了南极大陆(86°S以北区域,后文所述南极冰盖均不包括86°S以南区域)在这5年里的冰雪质量变化情况,得到东南极冰盖的质量变化为-18±20Gt/a,西南极-26±6Gt/a,南极冰盖的冰雪质量变化为-44±21Gt/a,对全球海平面上升的影响约为0.12mm·a~(-1).解算结果表明,南极冰盖质量亏损主要集中在西南极阿蒙森海岸附近冰川以及东南极波因塞特角区域.     

Climates of the Earth and Cryosphere Evolution

The interrelationship between the cryosphere and the climate is not always operating on Earth over a scale of billions or millions of years. Indeed, most of the time, the Earth is regulated at temperatures such that no ice sheet exists. Nevertheless, it is very fruitful to understand the conditions where and when ice sheets were triggered during the Earth??s history. This paper deals with the paleoclimate and the cryosphere in the last 4.6 Ga and explains the different processes that make the climate of the first 4 billion years warm despite the weaker solar luminosity. We also describe the more recent evolution in the last 65 million years when a global decrease in atmospheric CO₂ from around 4 PAL to 1 PAL was associated with a global cooling (1 PAL present atmospheric level = 280 ppm). It is in this context that the Quaternary occurred characterized by low atmospheric CO₂ and the presence of two perennial ice sheets in Greenland and Antarctica. The last million years are certainly the most documented since direct and reliable CO₂ measurements are available. They are characterized by a complex climate/cryosphere dynamics leading to oscillations between long glacial periods with four ice sheets and shorter ones with only two ice sheets (interglacial). We are currently living in one of those interglacials, generally associated with a CO₂ level of 280 ppm. Presently, anthropogenic activities are seriously perturbing the carbon cycle and the atmospheric CO₂ content and therefore the climate. The last but not least question raised in this paper is to investigate whether the anthropogenic perturbation may lead to a melting of the ice sheets.     

Mass balance of the Lambert Glacier basin, East Antarctica

Since it is the largest glacier system in Antarctica, the Lambert Glacier basin plays an important role in the mass balance of the overall Antarctic ice sheet. The observed data and shallow core studies from the inland traverse investigations in recent years show that there are noticeable differences in the distribution and variability of the snow accumulation rate between east and west sides. On the east side, the accumulation is higher on the average and has increased in the past decades, while on the west side it is contrary. The ice movement measurement and the ice flux calculation indicate that the ice velocity and the flux are larger in east than in west, meaning that the major part of mass supply for the glacier is from the east side. The mass budget estimate with the latest data gives that the integrated accumulation over the upstream area of the investigation traverse route is larger than the outflow ice flux by 13%, suggesting that the glacier basin is in a positive mass balance state and the ice thickness will increase if the present climate is keeping.     

Dating a 109.9 m ice core from Dome A(East Antarctica) with volcanic records and a firn densification model

A 109.9 m ice core was extracted at a location about 300 m away from the Dome A summit(80°00′S,77°21″E) by the Chinese team of the International Trans-Antarctic Science Expedition(ITASE) during the 21st Chinese National Antarctica Research Expedition(CHINARE) in January 2005.Two independent methods were used for dating the ice core,volcanic event markers shown by prominent non-sea-salt sulfate(nss-SO42) and the Herron and Langway(H-L) firn densification model.Six promi-nent volcanic events(Agung 1963 AD,Tambora 1815 AD,Kuwae 1453 AD,Unknown 1259 AD,Taupo 186 AD and Pinatubo 1050 BC) were identified by comparison with other Antarctic ice cores.Based on the mean accumulation rates be-tween adjacent events,we estimate the age at the firn pore close-off depth(102 m) was 3516±100 a BP.This is the oldest close-off age ever reported from the Antarctic and the Greenland ice sheets.Calculations using the H-L model show that the age at the same depth is 3581±100 a BP.The two dating techniques differ by 65 years,or～1.8% of the record.We calculated the bottom age of the ice core as 4009±150 a BP using the volcanic dating method and 4115±150 a BP using the H-L model method.     

Observed Mass Balance of Mountain Glaciers and Greenland Ice Sheet in the 20th Century and the Present Trends

Glacier mass balance and secular changes in mountain glaciers and ice caps are evaluated from the annual net balance of 137 glaciers from 17 glacierized regions of the world. Further, the winter and summer balances for 35 glaciers in 11 glacierized regions are analyzed. The global means are calculated by weighting glacier and regional surface areas. The area-weighted global mean net balance for the period 1960?C2000 is ?270 ± 34 mm a^?1 w.e. (water equivalent, in mm per year) or (?149 ± 19 km³ a^?1 w.e.), with a winter balance of 890 ± 24 mm a^?1 w.e. (490 ± 13 km³ a^?1 w.e.) and a summer balance of ?1,175 ± 24 mm a^?1 w.e. (?647 ± 13 km³ a^?1 w.e.). The linear-fitted global net balance is accelerating at a rate of ?9 ± 2.1 mm a^?2. The main driving force behind this change is the summer balance with an acceleration of ?10 ± 2.0 mm a^?2. The decadal balance, however, shows significant fluctuations: summer melt reached its peak around 1945, followed by a decrease. The negative trend in the annual net balance is interrupted by a period of stagnation from 1960s to 1980s. Some regions experienced a period of positive net balance during this time, for example, Europe. The balance has become strongly negative since the early 1990s. These decadal fluctuations correspond to periods of global dimming (for smaller melt) and global brightening (for larger melt). The total radiation at the surface changed as a result of an imbalance between steadily increasing greenhouse gases and fluctuating aerosol emissions. The mass balance of the Greenland ice sheet and the surrounding small glaciers, averaged for the period of 1950?C2000, is negative at ?74 ± 10 mm a^?1 w.e. (?128 ± 18 km³ a^?1 w.e.) with an accumulation of 297 ± 33 mm a^?1 w.e. (519 ± 58 km³ a^?1 w.e.), melt ablation ?169 ± 18 mm a^?1 w.e. (?296 ± 31 km³ a^?1 w.e.), calving ablation ?181 ± 19 mm a^?1 w.e. (?316 ± 33 km³ a^?1 w.e.) and the bottom melt-21 ± 2 mm a^?1 w.e. (?35 ± 4 km³ a^?1 w.e.). Almost half (?60 ± 3 km³ a^?1) of the net mass loss comes from mountain glaciers and ice caps around the ice sheet. At present, it is difficult to detect any statistically significant trends for these components. The total mass balance of the Antarctic ice sheet is considered to be too premature to evaluate. The estimated sea-level contributions in the twentieth Century are 5.7 ± 0.5 cm by mountain glaciers and ice caps outside Antarctica, 1.9 ± 0.5 cm by the Greenland ice sheet, and 2 cm by ocean thermal expansion. The difference of 7 cm between these components and the estimated value with tide-gage networks (17 cm) must result from other sources such as the mass balance of glaciers of Antarctica, especially small glaciers separated from the ice sheet.