Projecting future sea level

Through the use of fossil fuels as an energy source, mankind is slowly changing the constitution of the atmosphere. The emission of CO₂ and other greenhouse gases changes the radiative properties of the earth/atmosphere system, and as a result climate is expected to become warmer. As a starting point for the sea-level rise scenario discussed here it is assumed that the globally-averaged increase of surface air temperatures will amount to 2 to 4°C in the second half of the next century (i.e. around 2085 AD). One of the consequences of this warming is an accelerated rise in sea level, caused by thermal expansion of ocean water and further retreat of mountain glaciers. The Greenland Ice Sheet will also decrease in size, but on the other hand, Antarctica is expected to grow slightly due to increased snowfall. Taken together, the projection for future sea level presented here suggest that by 2085 AD, global sea-level stand will be 28–66 cm higher than the present level, which implies a rate of sea-level rise of about 2 to 4 times that observed during the last 100 yr. Our scenario does not include a contribution resulting from the possible collapse of the West Antarctic Ice Sheet. If this collapse is indeed likely to occur after the major peripheral ice shelves have thinned considerably, the effects on sea level will be small in the coming 100 yr. First, the oceans surrounding Antarctica must have warmed sufficiently to reduce the winter sea-ice extent to allow circumpolar deep water to penetrate into the sub-shelf cavities, thus increasing basal melt rates on the ice shelves. Of course, on longer time scales, West Antarctica could become the major contributor to rising sea level.     

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.     

冰盖消融的海平面指纹变化及其对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更精确地估算研究区冰盖质量变化时,应考虑海平面指纹效应的渗透影响.

     

利用ICESat数据确定格陵兰冰盖高程和体积变化

两极冰盖消融是造成海平面上升的重要原因,作为世界第二大冰盖,格陵兰冰盖消融速度在进入21世纪以后明显加快,引起了广泛关注.本文利用ICESat卫星激光测高数据,探讨了坡度改正的方法,通过改进平差模型解决了病态问题,并采用重复轨道方法计算了2003年9月至2009年10月间格陵兰冰盖的体积和高程变化趋势,对格陵兰冰盖各冰川流域系统的变化情况进行了详细分析.结果表明,格陵兰冰盖在这6年间平均高程变化趋势为-16.79±0.84 cm·a^-1,体积变化速率为-301.37±15.16 km³·a^-1,体积流失主要发生在冰盖边缘,其中DS1、DS8等流域的体积损失正在加剧,而高程在2000 m以上的冰盖内陆地区表现出高程积聚的状态,但增长速度明显减缓.与现有研究成果的对比表明,算法优化后的本文结果更具可靠性.

     

A Review of Recent Updates of Sea-Level Projections at Global and Regional Scales

Sea-level change (SLC) is a much-studied topic in the area of climate research, integrating a range of climate science disciplines, and is expected to impact coastal communities around the world. As a result, this field is rapidly moving, and the knowledge and understanding of processes contributing to SLC is increasing. Here, we discuss noteworthy recent developments in the projection of SLC contributions and in the global mean and regional sea-level projections. For the Greenland Ice Sheet contribution to SLC, earlier estimates have been confirmed in recent research, but part of the source of this contribution has shifted from dynamics to surface melting. New insights into dynamic discharge processes and the onset of marine ice sheet instability increase the projected range for the Antarctic contribution by the end of the century. The contribution from both ice sheets is projected to increase further in the coming centuries to millennia. Recent updates of the global glacier outline database and new global glacier models have led to slightly lower projections for the glacier contribution to SLC (7–17 cm by 2100), but still project the glaciers to be an important contribution. For global mean sea-level projections, the focus has shifted to better estimating the uncertainty distributions of the projection time series, which may not necessarily follow a normal distribution. Instead, recent studies use skewed distributions with longer tails to higher uncertainties. Regional projections have been used to study regional uncertainty distributions, and regional projections are increasingly being applied to specific regions, countries, and coastal areas.     

Contemporary (1960–2012) Evolution of the Climate and Surface Mass Balance of the Greenland Ice Sheet

We assess the contemporary (1960–2012) surface mass balance (SMB) of the Greenland ice sheet (GrIS), its individual components and trends. We use output of the high-resolution (11 km) regional atmospheric climate model (RACMO2), evaluated with automatic weather stations and GRACE data. A persistent negative North Atlantic oscillation index over the last 6 years resulted in the summertime advection of relatively warm continental air toward the GrIS. Added to the enhanced radiative forcing by increased CO₂ levels, this has resulted in an increase in near-surface temperature of more than 2 K during 2007–2012 compared to 1960–1990. The associated decrease in albedo led to an extra absorption of shortwave radiation of ～6 Wm^?2 (11 %) in the summer months, which is the main driver of enhanced surface melting and runoff in recent years. From 1990 onward, we see a steady increase in meltwater runoff and an associated decrease in the SMB, accelerating after 2005, with the record low SMB year in 2010. Despite the fact that the GrIS was subject to the highest surface melt rates in 2012, relatively high accumulation rates prevented 2012 to set a record low SMB. In 2012, melt occurred relatively high on the ice sheet where melt water refreezes in the porous firn layer. Up to 2005, increased runoff was partly offset by increased accumulation rates. Since then, accumulation rates have decreased, resulting in low SMB values. Other causes of decreased SMB are the loss of firn pore space and decreasing refreezing rates in the higher ablation area. The GrIS has lost in total 1,800 ± 300 Gt of mass from surface processes alone since 1990 and about half of that in the last 6 years.     

Understanding and Modelling Rapid Dynamic Changes of Tidewater Outlet Glaciers: Issues and Implications

Recent dramatic acceleration, thinning and retreat of tidewater outlet glaciers in Greenland raises concern regarding their contribution to future sea-level rise. These dynamic changes seem to be parallel to oceanic and climatic warming but the linking mechanisms and forcings are poorly understood and, furthermore, large-scale ice sheet models are currently unable to realistically simulate such changes which provides a major limitation in our ability to predict dynamic mass losses. In this paper we apply a specifically designed numerical flowband model to Jakobshavn Isbrae (JIB), a major marine outlet glacier of the Greenland ice sheet, and we explore and discuss the basic concepts and emerging issues in our understanding and modelling ability of the dynamics of tidewater outlet glaciers. The modelling demonstrates that enhanced ocean melt is able to trigger the observed dynamic changes of JIB but it heavily relies on the feedback between calving and terminus retreat and therefore the loss of buttressing. Through the same feedback, other forcings such as reduced winter sea-ice duration can produce similar rapid retreat. This highlights the need for a robust representation of the calving process and for improvements in the understanding and implementation of forcings at the marine boundary in predictive ice sheet models. Furthermore, the modelling uncovers high sensitivity and rapid adjustment of marine outlet glaciers to perturbations at their marine boundary implying that care should be taken in interpreting or extrapolating such rapid dynamic changes as recently observed in Greenland.     

Shortening the recurrence periods of extreme water levels under future sea-level rise

Sea-level rise, as a result of global warming, may lead to more natural disasters in coastal regions where there are substantial aggregations of population and property. Thus, this paper focuses on the impact of sea-level rise on the recurrence periods of extreme water levels fitted using the Pearson type III (P-III) model. Current extreme water levels are calculated using observational data, including astronomical high tides and storm surges, while future extreme water levels are determined by superposing scenario data of sea-level rise onto current extreme water levels. On the basis of a case study using data from Shandong Province, China, results indicated that sea-level rise would significantly shorten the recurrence periods of extreme water levels, especially under higher representative concentration pathway (RCP) scenarios. Results showed that by the middle of the century, 100-year current extreme water levels for all stations would translate into once in 15–30 years under RCP 2.6, and once in ten to 25 years under RCP 8.5. Most seriously, the currently low probability event of a 1000-year recurrence would become common, occurring nearly every 10 years by 2100, based on projections under RCP 8.5. Therefore, according to this study, corresponding risk to coastlines could well be increase in future, as the recurrence periods of extreme water levels would be shortened with climate change.     

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.     

Response of a Coupled Ocean�CAtmosphere Model to Greenland Ice Melting

We investigate the transient response of the global coupled ocean?Catmosphere system to enhanced freshwater forcing representative of melting of the Greenland ice sheets. A 50-year long simulation by a coupled atmosphere?Cocean general circulation model (CGCM) is compared with another of the same length in which Greenland melting is prescribed. To highlight the importance of coupled atmosphere?Cocean processes, the CGCM results are compared with those of two other experiments carried out with the oceanic general circulation model (OGCM). In one of these OGCM experiments, the prescribed surface fluxes of heat, momentum and freshwater correspond to the unperturbed simulation by the CGCM; in the other experiment, Greenland melting is added to the freshwater flux. The responses by the CGCM and OGCM to the Greenland melting have similar patterns in the Atlantic, albeit the former having five times larger amplitudes in sea surface height anomalies. The CGCM shows likewise stronger variability in all state variables in all ocean basins because the impact of Greenland melting is quickly communicated to all ocean basins via atmospheric bridges. We conclude that the response of the global climate to Greenland ice melting is highly dependent on coupled atmosphere?Cocean processes. These lead to reduced latent heat flux into the atmosphere and an associated increase in net freshwater flux into the ocean, especially in the subpolar North Atlantic. The combined result is a stronger response of the coupled system to Greenland ice sheet melting.     

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.     

Estimating the Glacier Contribution to Sea-Level Rise for the Period 1800�C2005

In this study, a new estimate of the contribution of glaciers and ice caps to the sea-level rise over the period 1800?C2005 is presented. We exploit the available information on changes in glacier length. Length records form the only direct evidence of glacier change that has potential global coverage before 1950. We calculate a globally representative signal from 349 glacier length records. By means of scaling, we deduce a global glacier volume signal, that is calibrated on the mass-balance and geodetic observations of the period 1950?C2005. We find that the glacier contribution to sea-level rise was 8.4 ± 2.1 cm for the period 1800?C2005 and 9.1 ± 2.3 cm for the period 1850?C2005.     

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.     

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.     

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

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

Representing Grounding Line Dynamics in Numerical Ice Sheet Models: Recent Advances and Outlook

Recent satellite observations of the Antarctic and Greenland ice sheets show accelerated ice flow and associated ice sheet thinning along coastal outlet glaciers in contact with the ocean. Both processes are the result of grounding line retreat due to melting at the grounding line (the grounding line is the contact of the ice sheet with the ocean, where it starts to float and forms an ice shelf or ice tongue). Such rapid ice loss is not yet included in large-scale ice sheet models used for IPCC projections, as most of the complex processes are poorly understood. Here we report on the state-of-the art of grounding line migration in marine ice sheets and address different ways in which grounding line migration can be attributed and represented in ice sheet models. Using one-dimensional ice flow models of the ice sheet/ice shelf system we carried out a number of sensitivity experiments with different spatial resolutions and stress approximations. These are verified with semi-analytical steady state solutions. Results show that, in large-scale finite-difference models, grounding line migration is dependent on the numerical treatment (e.g. staggered/non-staggered grid) and the level of physics involved (e.g. shallow-ice/shallow-shelf approximation).     

Sea-level rise: A review of recent past and near-future trends

Global mean sea level is a potentially sensitive indicator of climate change. Global warming will contribute to worldwide sea-level rise (SLR) from thermal expansion of ocean water, melting of mountain glaciers and polar ice sheets. A number of studies, mostly using tide-gauge data from the Permanent Service for Mean Sea Level, Bidston Observatory, England, have obtained rates of global SLR within the last 100 years that range between 0·3 and 3 mm yr^?1, with most values concentrated between 1 and 2 mm yr^?1. However, the reliability of these results has been questioned because of problems with data quality and physical processes that introduce a high level of spatial and temporal variability. Sources of uncertainty in the sea-level data include variations in winds, ocean currents, river runoff, vertical earth movements, and geographically uneven distribution of long-term records. Crustal motions introduce a major source of error. To a large extent, these can be filtered by employing palaeo-sea-level proxies, and geophysical modelling to remove glacio-isostatic changes. Ultimately, satellite geodesy will help resolve the inherent ambiguity between the land and ocean level changes recorded by tide gauges. Future sea level is expected to rise by ～ 1 m, with a ‘best-guess’ value of 48 cm by the year 2100. Such rates represent an acceleration of four to seven times over present rates. Local land subsidence could substantially increase the apparent SLR. For example, Louisiana is currently experiencing SLR trends nearly 10 times the global mean rate. These recently reduced SLR estimates are based on climate models that predict a zero to negative contribution to SLR from Antarctica. Most global climate models (GCMs) indicate an ice accumulation over Antarctica, because in a warmer world, precipitation will exceed ablation/snow-melt. However, the impacts of attritional processes, such as thinning of the ice shelves, have been downplayed according to some experts. Furthermore, not all climate models are in agreement. Opposite conclusions may be drawn from the results of other GCMs. In addition, the West Antarctic Ice Sheet is potentially subject to dynamic and volcanic instabilities that are difficult to predict. Because of the great uncertainty in SLR projections, careful monitoring of future sea-level trends by upgraded tide-gauge networks and satellite geodesy will become essential. Finally, because of the high spatial variability in crustal subsidence rates, wave climates and tidal regimes, it will be the set of local conditions (especially the relative sea-level rise), rather than a single global mean sea-level trend, that will determine each locality's vulnerability to future SLR.     

The Deformational Response of a Viscoelastic Solid Earth Model Coupled to a Thermomechanical Ice Sheet Model

We apply a coupled thermomechanical ice sheet—self-gravitating viscoelastic solid Earth model (SGVEM), allowing for the dynamic exchange of ice thickness and bedrock deformation, in order to investigate the effect of viscoelastic deformation on ice dynamics and vice versa. In a synthetic glaciation scenario, we investigate the interaction between the ice sheet and the solid Earth deformation, the glacial-isostatic adjustment (GIA), accounting for an atmospheric forcing depending on the ice sheet surface altitude. We compare the results from the coupled model to runs with the common elastic lithosphere/relaxing asthenosphere (ELRA) model, where the lithosphere is represented by a thin plate and the mantle relaxes with one characteristic relaxation time, as well as to a rigid Earth without any deformation. We find that the deformational behaviour of the SGVEM on ice dynamics (i.e. stored ice volume, ice thickness and velocity field) is comparable to the ELRA for an optimal choice of the parameters in steady state, but exhibits differences in the transient behaviour. Beyond the ice sheet, in the region of peripheral forebulge, the differences in the transient surface deformation between ELRA and SGVEM are substantial, demonstrating the inadequacy of the ELRA model for interpreting constraints on GIA in the periphery of the ice sheet, such as sea-level indicators and GPS uplift rates.     

E) Organisations Nationales

Abstract

Sudden rapid advances or surges of glaciers and sections of smail ice caps are well known. After remaining dormant or in retreat over long periods of time these ice masses suddenly move forward rapidly with speeds about 2 orders of magnitude greater than usual. If such surges were to occur in large sections of the Antarctic ice sheet serious consequences could result. These include a significant rise of sea level, a substantial increase in the high-albedo ice cover around the continent especially in summer, and a cooling of the Antarctic ocean by the additional ice melting.

A numerical model has now been developed which simulates surging in certain glaciers and ice sheets in an apparently realistic manner. This model has been found to give close representations to a number of existing real surging temperate ice masses from small mountain glaciers to large sectors of ice caps. The model reproduces realistically many features of these ice masses such as the period of the surge, the duration, the velocity of advance, the magnitude of the advance, and the changes in ice thickness.

The application of the model to the Antarctic ice sheet is made more difficult by the problem caused by the temperature dependence of the flow properties of ice. This means that for a complete study the interaction with the environment needs to be considered. However, at this stage preliminary calculations indicate a number of features that are relevant to the effect of Antarctic ice surges on the global climate. These include the period between surges, the duration of the surge, the amount of ice advanced and the changes in thickness of the ice sheet.     

Evaluating sun–climate relationships since the Little Ice Age

From the coldest period of the Little Ice Age to the present time, the surface temperature of the Earth increased by perhaps 0.8°C. Solar variability may account for part of this warming which, during the past 350 years, generally tracks fluctuating solar activity levels. While increases in greenhouse gas concentrations are widely assumed to be the primary cause of recent climate change, surface temperatures nevertheless varied significantly during pre-industrial periods, under minimal levels of greenhouse gas variations. A climate forcing of 0.3 W m⁻² arising from a speculated 0.13% solar irradiance increase can account for the 0.3°C surface warming evident in the paleoclimate record from 1650 to 1790, assuming that climate sensitivity is 1°C W⁻¹ m⁻² (which is within the IPCC range). The empirical Sun–climate relationship defined by these pre-industrial data suggests that solar variability may have contributed 0.25°C of the 0.6°C subsequent warming from 1900 to 1990, a scenario which time dependent GCM simulations replicate when forced with reconstructed solar irradiance. Thus, while solar variability likely played a dominant role in modulating climate during the Little Ice Age prior to 1850, its influence since 1900 has become an increasingly less significant component of climate change in the industrial epoch. It is unlikely that Sun–climate relationships can account for much of the warming since 1970, not withstanding recent attempts to deduce long term solar irradiance fluctuations from the observational data base, which has notable occurrences of instrumental drifts. Empirical evidence suggests that Sun–climate relationships exist on decadal as well as centennial time scales, but present sensitivities of the climate system are insufficient to explain these short-term relationships. Still to be simulated over the time scale of the Little Ice Age to the present is the combined effect of direct radiative forcing, indirect forcing via solar-induced ozone changes in the atmosphere, and speculated charged particle mechanisms whose pathways and sensitivities are not yet specified.