bluuebunny/crossref_metadata_2025 huggingface dataset:
Development and validation of a snow albedo algorithm for the MODIS instrument
Authors: Klein Andrew G. , Stroeve Julienne
Published: November 2007

A prototype snow albedo algorithm has been developed for the Moderate Resolution Imaging Spectroradiometer (MODIS). It complements existing MODIS products by providing albedo measurements for areas mapped as snow on a global daily basis by MODIS. Cloud detection and atmospheric correction are accomplished using existing MODIS products. Models of the bidirectional reflectance of snow created using a discrete-ordinate radiative transfer (DISORT) model are used to correct for anisotropic scattering effects over non-forested surfaces. Initial algorithm validation is undertaken through comparisons with broadband albedo measurements made at the U.S. National Oceanic and Atmospheric Administration (NOAA) Surface Radiation Budget Network (SURFRAD) site in Fort Peck, MT. In situ SURFRAD albedo measurements are compared to daily MODIS snow albedo retrievals for the period 21–26 November 2000 created from five narrow-to-broadband albedo conversion schemes. The prototype MODIS algorithm produces reasonable broadband albedo estimates. Maximum daily differences between the five MODIS broadband albedo retrievals and in situ albedo are 15%. Daily differences between the best MODIS broadband estimate and the measured SURFRAD albedo are 1–8%. However, no single conversion scheme consistently provides the closest albedo estimate. Further validation and algorithm development using data from North America and Greenland is ongoing.
A 21 year record of Arctic sea-ice extents and their regional, seasonal and monthly variability and trends
Authors: Parkinson Claire L. , Cavalieri Donald J.
Published: November 2007

Satellite passive-microwave data have been used to calculate sea-ice extents over the period 1979–99 for the north polar sea-ice cover as a whole and for each of nine regions. Over this 21 year time period, the trend in yearly-average ice extents for the ice cover as a whole is –32 900±6100 km2 a–1 (–2.7 ±0.5% per decade), indicating a statistically significant reduction in sea-ice coverage. Regionally, the reductions are greatest in the Arctic Ocean, the Kara and Barents Seas and the Seas of Okhotsk and Japan; and seasonally, the reductions are greatest in summer, for which season the 1979–99 trend in ice extents is –41600±12 900 km2 a–1 (–4.9±1.5% per decade). On a monthly basis, the reductions are greatest in July and September for the north polar ice cover as a whole, in September for the Arctic Ocean, in June and July for the Kara and Barents Seas, and in April for the Seas of Okhotsk and Japan. Of the nine regions, only the Bering Sea and the Gulf of St Lawrence show positive ice-extent trends on a yearly-average basis. However, the increases in these two regions are not statistically significant. For the north polar region as a whole, and for the Arctic Ocean, the Seas of Okhotsk and Japan, and Hudson Bay, the negative trends in the yearly averages are statistically significant at a 99% confidence level.
Assessment of the relative accuracy of hemispheric-scale snow-cover maps
Authors: Hall Dorothy K. , Kelly Richard E. J. , Riggs George A. , Chang Alfred T. C. , Foster James L.
Published: November 2007

There are several hemispheric-scale satellite-derived snow-cover maps available, but none has been fully validated. For the period 23 October–25 December 2000, we compare snow maps of North America derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) and operational snow maps from the U.S. National Oceanic and Atmospheric Administration (NOAA) National Operational Hydrologic Remote Sensing Center (NOHRSC), both of which rely on satellite data from the visible and near-infrared parts of the spectrum; we also compare MODIS maps with Defense Meteorological Satellite Program (DMSP) Special Sensor Microwave/Imager (SSM/I) passive-microwave snow maps for the same period. The maps derived from visible and near-infrared data are more accurate for mapping snow cover than are the passive-microwave-derived maps, but discrepancies exist as to the location and extent of the snow cover even between operational snow maps. The MODIS snow-cover maps show more snow in each of the 8 day periods than do the NOHRSC maps, in part because MODIS maps the effects of fleeting snowstorms due to its frequent coverage. The large (~30 km) footprint of the SSM/I pixel, and the difficulty in distinguishing wet and shallow snow from wet or snow-free ground, reveal differences up to 5.33 x 106 km2 in the amount of snow mapped using MODIS vs SSM/I data. Algorithms that utilize both visible and passive-microwave data, which would take advantage of the all-weather mapping capability of the passive-microwave data, will be refined following the launch of the Advanced Microwave Scanning Radiometer (AMSR) in the fall of 2001.
Time-series analysis of passive-microwave-derived central North American snow water equivalent imagery
Authors: Derksen C. , Walker A. , LeDrew E. , Goodison B.
Published: November 2007

The Meteorological Service of Canada has developed a series of operational snow water equivalent (SWE) retrieval algorithms for central Canada, based on the vertically polarized difference index for the 19 and 37 GHz channels of the Special Sensor Microwave/Imager (SSM/I). Separate algorithms derive SWE for open environments, deciduous, coniferous and sparse forest cover. A final SWE value represents the area-weighted average based on the proportional land cover within each pixel. In this study, 5 day averaged (pentad) passive-microwave-derived SWE imagery for the winter season (December–February) of 1994/95 is compared to in situ data from central Canada in order to assess algorithm performance. Investigation of regions with varying proportional land cover within the four algorithm classes shows that retrieved SWE remains within ±10–20mm of surface observations, independent of fractional within-pixel land cover. Following algorithm evaluation, ten winter seasons (1988/89 through 1997/98) of pentad central North American SWE imagery are subjected to a rotated principal-component analysis (PCA). Although there are no trends in total study-area SWE, the PCA results identify the interseasonal variability in the SWE accumulation and ablation centers of action through the SSM/I time series.
Mapping of permafrost on Vega Island, Antarctic Peninsula, using satellite images and aerial photography
Authors: Ermolin Evgeniy , Angelis Hernán De , Skvarca Pedro
Published: November 2007

The work presented deals with detailed mapping of permafrost in Devil Bay, Vega Island, northeastern Antarctic Peninsula. Mapping of landforms and permafrost features within a periglacial plain was performed using high resolution visible satellite images and aerial photography. Two maps of permafrost were produced: one based on analysis and interpretation of visible satellite imagery and another on low-altitude aerial photography, both yielding similar results. The principles of morphogenesis were applied to map production, distinguishing both the syncryogenic and epicryogenic formations and each constitutive cryofacies. The interpretation of remote-sensing data allowed areas of occurrence of particular cryogenic processes to be defined. Remote sensing was found to be useful for permafrost mapping at both medium and large scales, and applicable for future extension to other regions in Antarctica.
Snow dunes and glazed surfaces in Antarctica: new field and remote-sensing data
Authors: Frezzotti Massimo , Gandolfi Stefano , Marca Floriana La , Urbini Stefano
Published: November 2007

As part of the International Trans-Antarctic Scientific Expedition project, the Italian Antarctic Programme undertook two traverses from the Terra Nova station to Talos Dome and to Dome C. Along the traverses, the party carried out several tasks (drilling, glaciological and geophysical exploration). The difference in spectral response between glazed surfaces and snow makes it simple to identify these areas on visible/near-infrared satellite images. Integration of field observation and remotely sensed data allows the description of different mega-morphologic features: wide glazed surfaces, sastrugi glazed surface fields, transverse dunes and megadunes. Topography global positioning system, ground penetrating radar and detailed snow-surface surveys have been carried out, providing new information about the formation and evolution of mega-morphologic features. The extensive presence, (up to 30%) of glazed surface caused by a long hiatus in accumulation, with an accumulation rate of nil or slightly negative, has a significant impact on the surface mass balance of a wide area of the interior part of East Antarctica. The aeolian processes creating these features have important implications for the selection of optimum sites for ice coring, because orographic variations of even a few metres per kilometre have a significant impact on the snow-accumulation process. Remote-sensing surveys of aeolian macro-morphology provide a proven, high-quality method for detailed mapping of the interior of the ice sheet’s prevalent wind direction and could provide a relative indication of wind intensity.
The regional distribution of the dry-snow zone on the Antarctic Peninsula north of 70˚ S
Authors: Rau Frank , Braun Matthias
Published: November 2007

In order to estimate the contribution of the Antarctic Peninsula to global sea-level rise as a result of the observed warming in this region, the spatial extent of snow-melt-producing areas needs to be quantified. By using the dry-snow line derived from synthetic aperture radar (SAR) imagery as the uppermost limit of frequent or occasional surface melt, an estimation of the spatial extent of areas with non-zero ablation rates is facilitated. Three calibrated RADARSAT ScanSAR mosaics covering the northern Antarctic Peninsula were analyzed applying a threshold of –14 dB to identify the dry-snow line. The area of the dry-snow radar zone was determined to be 23300±2000 km2.Areas affected frequently or occasionally by snowmelt cover 85 000±9000 km2. In addition, the dry-snow line as derived from multi-temporal ERS-1/-2 imagery serves as an indicator of climate variability in the uppermost areas of polar glaciers and ice sheets. The upward shift of the dry-snow line between 1992 and 1998 on the eastern side of the Antarctic Peninsula at 68˚ S is interpreted as a direct response to the increasing number of high-temperature events during the 1991–2000 decade.
In situ and satellite surface temperature records in Antarctica
Authors: Shuman Christopher A. , Comiso Josefino C.
Published: November 2007

Air-temperature records (TA) during 1992 from five inland Antarctic automatic weather station (AWS) sites were compared with the best available infrared temperatures (TIR) from the Advanced Very High Resolution Radiometer (AVHRR) as well as calibrated passive-microwave temperatures (TC) from the Special Sensor Microwave/Imager (SSM/I). Daily and monthly average TA, TIR, and TC data indicate that each approach captures generally similar trends at each site but each approach also has limitations. AWS TA data are considered the most accurate but represent spatially restricted areas and may have long gaps due to sensor or transmission problems. AVHRR TIR data have daily variability similar to the TA record but have numerous small gaps due to cloud cover or observation interruptions. An offset between TA and TIR (>4 K) at the South Pole site was identified that may be due to the inclusion of data with large satellite scan angles necessary to cover this area. SSM/I TC data have the most continuity but exhibit calibration problems, a significantly damped daily response and do not cover all of Antarctica. Individual daily differences between TA and TIR as well as TA and TC can exceed 17 K, but all sites have mean daily differences of about 1 Kor better, after compensating for the offset at South Pole, and standard deviations of <6K. Monthly temperature differences are typically 5 K or better, with standard deviations generally <3K. And finally, using the available data, the 1992 average temperature differences are <1 K.
RADARSAT-1 Antarctic Mapping Project: change-detection and surface velocity campaign
Authors: Jezek Kenneth C.
Published: November 2007

The RADARSAT-1 Antarctic Mapping Project (RAMP) is a collaboration between NASA and the Canadian Space Agency to map Antarctica using synthetic aperture radar (SAR). The first Antarctic Mapping Mission (AMM-1) was successfully completed in October 1997. Data from the acquisition phase of the 1997 campaign have been used to achieve the primary goal of producing the first high-resolution SAR image map of Antarctica. The Modified Antarctic Mapping Mission (MAMM) occurred during the fall of 2000. The acquisition strategy concentrated on collecting highest- resolution RADARSAT-1 data of Antarctica’s fast glaciers for change detection and feature-retracking estimates of surface velocity. Additionally, extensive data were acquired for interferometric analysis over the entire viewable region, which extends north of 80.1˚ S latitude. This paper summarizes the goals and strategy behind MAMM. It goes on to discuss ice-sheet margin changes observed on several ice shelves around the continent. Margin changes are documented by comparing AMM-1 and MAMM data with earlier datasets including European Remote-sensing Satellite-1 SAR imagery, Landsat imagery, the Antarctic Digital Database (version 1) coastline and Declassified Intelligence Satellite Photography. Analysis reveals a complex pattern of ice-margin advance and retreat without indicating any systematic behavior in ice-sheet extent about the ice-sheet perimeter.
Correlation and trend studies of the sea-ice cover and surface temperatures in the Arctic
Authors: Comiso Josefino C.
Published: November 2007

Co-registered and continuous satellite data of sea-ice concentrations and surface ice temperatures from 1981 to 2000 are analyzed to evaluate relationships between these two critical climate parameters and what they reveal in tandem about the changing Arctic environment. During the 19 year period, the Arctic ice extent and actual ice area are shown to be declining at a rate of –2.0±0.3% dec –1 and 3.1 ±0.4% dec–1, respectively, while the surface ice temperature has been increasing at 0.4 ±0.2 K dec–1, where dec is decade. The extent and area of the perennial ice cover, estimated from summer minimum values, have been declining at a much faster rate of –6.7±2.4% dec–1 and –8.3±2.4% dec–1, respectively, while the surface ice temperature has been increasing at 0.9 ±0.6K dec–1. This unusual rate of decline is accompanied by a very variable summer ice cover in the 1990s compared to the 1980s, suggesting increases in the fraction of the relatively thin second-year, and hence a thinning in the perennial, ice cover during the last two decades. Yearly anomaly maps show that the ice-concentration anomalies are predominantly positive in the 1980s and negative in the 1990s, while surface temperature anomalies were mainly negative in the 1980s and positive in the 1990s. The yearly ice-concentration and surface temperature anomalies are highly correlated, indicating a strong link especially in the seasonal region and around the periphery of the perennial ice cover. The surface temperature anomalies also reveal the spatial scope of each warming (or cooling) phenomenon that usually extends beyond the boundaries of the sea-ice cover.
The influence of temperature on the small-strain viscous deformation mechanics of snow: a comparison with polycrystalline ice
Authors: Scapozza Carlo , Bartelt Perry A.
Published: November 2007

Glen’s law is commonly used to model the viscous deformation of polycrystalline ice. It is a power law that relates stress to viscous strain rate and contains three material parameters: n, a power-law exponent, Q, an activation energy, and A0, a material constant. Because polycrystalline ice is the constituent material of snow, it is to be expected that the viscous deformation mechanics of snow are related to the viscous behaviour of polycrystalline ice, especially under small strains and low strain rates when kinematic effects in the ice matrix like bond breakage, bond formation and grain sliding are of secondary importance. Based on 64 deformation-controlled compression tests on fine-grained snow in the density range 200–430kg m–3 and temperature range T = –20 to –2°C, we show that Glen’s law—with material parameters similar to those for polycrystalline ice—can be applied to model the viscous deformation of high-density snow However, the values of the ice material parameters are valid for densities above a relatively low density of 400 kg m–3; they are not valid for snow with densities below 360 kg m–3. We present the variation of n, Q and A for snow as a function of density and temperature. A possible explanation for this behaviour is that the ice grains in low-density snow are less constrained. Therefore, deformation mechanisms, such as grain-boundary sliding, increase in overall importance, leading to smaller n values and higher activation energies, Q. Although the material behaviour of low-density snow can be accurately modelled using a power law, the power-law parameters depart substantially from those of polycrystalline ice. The large variation of n and Q with temperature and density underscores the difficulty of predicting snow avalanches.
Simulation of fast-flow features of the Fennoscandian ice sheet during the Last Glacial Maximum
Authors: Forsström Pirjo-Leena , Sallasmaa Olli , Greve Ralf , Zwinger Thomas
Published: November 2007

In order to reconstruct the palaeoglaciation in Fennoscandia and northern Asia during the late-Weichselian ice-age phase, simulations with the dynamic and thermodynamic ice-sheet model SICOPOLIS are carried out. Our focus is on the Last Glacial Maximum (LGM) around 20 kyr BP. Climate forcing is based on mean annualsur-face temperature and precipitation derived from present data and Palaeoclimatic Modelling Intercomparison Project (PMIP) UKMO21 results for the LGM. These distributions are interpolated via a glacial index defined by the Greenland Icecore Project (GRIP) δ18O record. The extent of the Scandinavian and the Barents ice sheets is reproduced in good agreement with the Quaternary Environments of the Eurasian North (QUEEN) reconstruction, but the Kara Sea and Taymyr Peninsula areas are excessively glaciated. The fast-flow regions derived from the simulations, which are generally connected to regions with a temperate base and temperate ice above, are compared to hypothesized palaeo-ice-stream locations, especially in the Norwegian Channel and the Baltic area. In the Norwegian Channel, temperate basal conditions with temperate ice above prevail and favour fast flow. In the Baltic area, ice-sheet advance is generally accompanied by slow ice velocities (<200 m a–1). Some temporary fast-flow features occur due to transitional temperate-base conditions, and higher velocities arise in retreat phases.
Laboratory studies of the flow rates of debris-laden ice
Authors: Jacka Tim H. , Shavawn Donoghue  , Li Jun , William F. Budd  , Anderson Ross M.
Published: November 2007

Ice-sheet basal ice is warmer than that above because of the heat from the Earth’s interior. The stresses acting on the basal ice are greatest. In addition, the basal ice often contains debris consisting of silt and small stones picked up from the rock over which the ice flows. Because the base is the warmest part of an ice sheet and the stress there is greatest, flow rates in the basal ice are large and often contribute most of the ice movement. It is therefore important, for accurate modelling of the ice sheets, to know whether the debris within the basal ice enhances or retards the flow of the ice. In this paper, we describe laboratory deformation tests in uniaxial compression and in simple shear, on sand-laden ice. We find no significant dependence of flow rate on sand content (up to 15% volume) in the stress range 0.13–0.5 MPa and temperature range –0.02 to –18.0°C. Further work needs to include laboratory tests on debris-laden ice extracted from the polar ice sheets. This work is underway.
Glacier changes in southern Spitsbergen, Svalbard,1901–2000
Authors: Pälli Anja , Moore John C. , Jania Jacek , Glowacki Piotr
Published: November 2007

High-resolution ground-penetrating radar surveys at 50 MHz on the polythermal glaciers Hornbreen, Hambergbreen and several surrounding glaciers in southern Spitsbergen, Svalbard, are presented and interpreted. Accurate positioning was obtained using differential global positioning system (DGPS). Digital elevation models (DEMs) of the bedrock and surface were constructed. Comparison of DGPS data and surface DEMs with data from the topographic mappings from 1936 oblique stereoscopic aerial photographs and from Mission Russe in 1899–1901 shows that the Hornbreen and Hambergbreen surfaces are about 60–100 m thinner today in the upper part than at the beginning of the 20th century. Hornbreen has retreated by 13.5 km from the central part of the front, and Hambergbreen by 16 km. All the fronts of the nearby east-coast glaciers in this area have retreated. The bedrock DEM shows that the Hornbreen and Hambergbreen beds lie at –25 to 25 m a.s.l. The combination of sub-sea-level fronts and increasing steepness of the glaciers suggests that the low-lying glaciated valley filled by Hornbreen and Hambergbreen may become a partially inundated ice-free isthmus within perhaps 100 years.
Seasonal variability in ice crystal properties at NorthGRIP: a case study around 301 m depth
Authors: Svensson Anders , Baadsager Pauli , Persson Asbjørn , Schøtt Hvidberg Christine , Siggaard-Andersen Marie-Louise
Published: November 2007

The aim of this case study is to quantify the seasonal variability in crystal properties and to discuss the reason for the variability. A continuous 1.10 m long vertical thin-section profile covering approximately five annual cycles has been obtained from the North Greenland Icecore Project (NorthGRIP) ice core at around 301 m depth. The crystal outline and the c-axis orientation of more than 13000 crystals in the profile have been measured on a new Australian automated ice-crystal analyzer. In 2.5 cm resolution we observe a strong seasonal variability in crystal areas of >30%deviation from the average value of 6.7 mm2. Each year, a band of smaller crystals is observed in ice deposited during spring. The area distribution function is found to be close to a lognormal distribution. The crystal areas are compared to the concentration of chemical impurities in the ice; at a 5 cm resolution, the best correlation is found with the concentration of Ca2+. Our results show no seasonal variability of the average c-axis orientation of ice crystals.
Anisotropic ice flow leading to the onset of Ice Stream D, West Antarctica: numerical modelling based on the observations from Byrd Station borehole
Authors: Wang Weili , Jay Zwally H. , Hulbe Christina L. , Siegert Martin J. , Joughin Ias
Published: November 2007

An ice-sheet flowline model is used to simulate the flow of ice along two particle paths toward the onset to Ice Stream D, West Antarctica. One path is near the centre line of the main tributary to the ice stream, while the second passes by the Byrd Station borehole site. In this paper, we analyze the flow of the moderately fast-flowing tributaries in terms of ice-fabric anisotropy and estimate the steady-state ice-flow regions with the compatible developed crystal orientation fabrics along two particle paths. Comparison between modelled isochrones and internal layers detected from radio-echo sounding surveys in the area is used to suggest that flow upstream of the onset to Ice Stream D appears to have been stable since at least the Last Glacial Maximum.
Properties of ice crystals in NorthGRIP late- to middle-Holocene ice
Authors: Svensson Anders , Schmidt Karen G. , Dahl-Jensen Dorthe , Johnsen Sigfús J. , Wang Yun , Kipfstuhl Sepp , Thorsteinn Thorsteinsson
Published: November 2007

Detailed measurements of crystal outlines and fabrics have been performed on 35 000 crystals in fifteen 10 × 20 cm2 vertical thin sections from the North Greenland Icecore Project (NorthGRIP) ice core, evenly distributed in the depth interval 115–880m. The crystals exhibit important changes over this period. As the ice gets older the mean crystal area increases towards a constant value, the shape of the crystals becomes increasingly irregular, and the area distribution of crystals develops from a single log-normal distribution into a bimodal lognormal distribution. The c-axis fabric of the ice shows a smooth development of an increasingly stronger vertical fabric with depth, and the formation of a weak vertical girdle. Already in the younger samples the fabric is rather strongly oriented towards vertical. The fabric and the area of individual crystals are found not to correlate. A simple model, which takes into account the vertical strain of the ice, is applied in an attempt to determine the crystal growth rate at NorthGRIP.
Effects of first-order stress gradients in an ice sheet evaluated by a three-dimensional thermomechanical coupled model
Authors: Saito Fuyuki , Abe-Ouchi Ayako , Blatter Heinz
Published: November 2007

A three-dimensional thermomechanically coupled ice-sheet model including calculation of the first-order stress gradients is developed to perform numerical studies on the effects of stress components neglected in the zeroth-order shallow-ice approximation. Steady-state solutions are obtained with both a shallow-ice and a first-order ice-sheet model for several idealistic, radially symmetric ice-sheet configurations. The results show that the effects of the normal deviatoric-stress gradients on ice thickness are generally small, but the influence on basal temperatures and age profiles at an ice divide is significant.
An ice crystal model for Jupiter’s moon Europa
Authors: K Schmidt Karen Guldbæ , Dahl-Jensen Dorthe
Published: November 2007

A simple model for crystal growth in the ice shell of Europa has been made in order to estimate the size of ice crystals at Europa’s surface. If mass is lost from the surface of Europa due to sputtering processes, and the ice thickness is constant in time, ice crystals will be transported upwards in the ice shell. The crystals will therefore grow under varying conditions through the shell. The model predicts that ice crystals are 4 cm– 80 m across at the surface. For the preferred parameter values, a crystal size of the order of 7 m is calculated.
Lakes and subglacial hydrological networks around Dome C, East Antarctica
Authors: Rémy Frédérique , Testut Laurent , Legrésy Benoît , Forieri Alessandro , Bianchi Cesido , Tabacco Ignazio E.
Published: November 2007

Precise topography from European Remote-sensing Satellite radar altimetry and high density of airborne radio-echo sounding in the area surrounding Dome C, Antarctica, show a link between surface features and subglacial lakes. In this paper, we extend the study to fine structures by computing a curvature-based coefficient (cy) related to surface undulations. These coefficient variations reveal many surface undulations, and some elongated features of this parameter seem to link known subglacial lakes. A population of high values of this coefficient, assumed to correspond to transitions between sliding and non-sliding flow regime, strengthen the appearance of a network which would link most of the lakes in the area. The existence of such a network impacts on ice-flow dynamics and on subglacial-lake studies.
The ice-thickness distribution of Unteraargletscher, Switzerland
Authors: Bauder Andreas , Funk Martin , Gudmundsson G. Hilmar
Published: November 2007

Results from 15 years of work on glacier-bed mapping by radio-echo soundings on Finsteraar-, Lauteraar- and Unteraargletscher, Switzerland, are summarized, and a new and greatly improved map of the ice-thickness distribution presented. In contrast to the tongue of Unteraargletscher, its two main tributaries, Lauteraar and Finsteraar, are both deep and narrow, and ice-thickness determination depends on the detection of more than just the primary reflection. Migrating the data led to considerably improved bed determination. Wherever possible, additional information on ice thicknesses gained from numerous hot-water drillings to the glacier bed is used as an independent verification of the results of the radar measurements, and a fair agreement is found.
Modeling the marine extent of Northern Hemisphere ice sheets during the last glacial cycle
Authors: Zweck Chris , Huybrechts Philippe
Published: November 2007

Mechanisms that determine time-dependent changes of the marine ice margin in dynamic ice-sheet models are important but poorly understood. Here we derive an empirical formulation for changes in the marine extent when modelling the Northern Hemisphere ice sheets over the last glacial cycle in a three-dimensional thermomechanically coupled ice-sheet model. We assume that the strongest control on changes in marine extent is ice calving, and that the variable most crucial to calving is water depth. The empirical marine-extent relationship is tuned so that the major marine-retreat history of the Laurentide and Eurasian ice sheets is modelled accurately in time and space. We find that this empirical treatment relating marine extent to water depth is sufficient to reproduce the observations, and discuss the implications for the physics of marine margin changes and the dynamics of the Northern Hemisphere ice sheets since the Last Glacial Maximum.
On the rheology of till
Authors: Fowler Andrew C.
Published: November 2007

The deformation of subglacial till is instrumental in causing certain glacier surges, the motion of rapid ice streams in ice sheets, and ice-sheet surges which are associated with Heinrich events, and consequent rapid climatic shifts, during the last ice age. It may also be the means whereby drumlins are formed, and these in turn may act as a brake on large-scale ice flow. It is therefore important in building models to understand the rate at which till deforms, and how this controls the basal ice velocity. In recent years, two paradigms have emerged. On the one hand, theoreticians have tended to use a viscous rheology, though this lacks quantitative support. On the other hand, field and laboratory studies suggest that till behaves plastically. In this paper I will examine some of the dynamic consequences of this latter assumption, and show how the dichotomy between viscous and plastic may be less clear-cut than previously thought.
Effects of basal sliding on isochrones and flow near an ice divide
Authors: Pettit Erin C. , Jacobson H. Paul , Waddington Edwin D.
Published: November 2007

If an ice sheet is frozen to its bed, deep ice directly under a divide experiences low deviatoric stress and is relatively hard, because the rheology of polar ice is described by a power-law constitutive relation. In steady state, stratigraphic layers tend to form an arch (“Raymond bump”) in this region. However, when the basal ice can slide, the stresses are redistributed, and longitudinal extension due to sliding is associated with increased deviatoric stress in the deep ice under the divide. This increased deviatoric stress weakens the tendency to form a Raymond bump. To find a realistic spatial distribution of sliding under an ice divide, we incorporate a thin layer of viscous till in a finite-element plane-strain flow model. The resulting basal “sliding” velocity varies approximately linearly with distance from the ice divide. By varying the till viscosity, we can adjust the amount of basal motion. We find that the Raymond bump decays exponentially with the fraction of total ice flux carried by sliding: the arch is 50% smaller when the sliding flux is only 7% of the total ice flux. This implies that the possibility of a wet bed must be considered when inferring past ice-divide locations from radar internal layering.
Spatial stability of Ice Stream D and its tributaries, West Antarctica, revealed by radio-echo sounding and interferometry
Authors: Siegert Martin J. , Payne Antony J. , Joughin Ian
Published: November 2007

It has been shown recently that ice streams are fed by fast-flowing tributaries occupying well-defined subglacial troughs and with shared source areas. Here, ice-penetrating radio-echo sounding (RES) data are analyzed in conjunction with ice surface velocities derived from interferometric synthetic aperture radar (InSAR), to determine the englacial properties of tributaries feeding Ice Stream D, West Antarctica. All of Ice Stream D’s tributaries are coincident with “buckled” internal ice-sheet layers, most probably deformed by the processes responsible for enhanced ice flow. Between the tributaries well-preserved internal layers occur. The data reveal that no lateral migration of the ice-stream tributaries has occurred recently. This is consistent with thermomechanical ice-flow modelling, which indicates that the flow of Ice Stream D is controlled by a subglacial trough and is unaffected by changes to the flow of neighbouring Ice Stream C.
Ice-fabrics study in the upper 1500 m of the Dome C (East Antarctica) deep ice core
Authors: Wang Yun , Kipfstuhl Sepp , Azuma Nobuhiko , Thorsteinsson Thorsteinn , Miller Heinz
Published: November 2007

A study of c-axis orientations in the upper 1500m of the Dome C (East Antarctica) deep ice core has been carried out using an automatic ice-fabric analyzer (AIFA). Twenty-nine vertical and a few horizontal thin sections from different depths in the core have been analyzed. Several statistical parameters describing fabric strength and fabric shapes have been calculated from the c-axis orientation data. The fabric diagrams display a near-random c-axis orientation distribution in the uppermost parts of the ice sheet. A tendency of c-axis rotation towards a broad single-maximum fabric is observed in the lowest part of the studied interval. The fabric development at Dome C thus appears typical for ice-sheet summit and dome sites. The fabric development at Dome C is compared with the fabric evolution in the Dome F and GRIP ice cores, and data on crystal size obtained with image-analysis techniques are presented. Studies of misorientation angles between the c axes of neighbouring crystals reveal little evidence for polygonization, but microscopic observations show that sub-grain boundaries are present in half of the grains at any depth.
Identifying isochrones in GPR profiles from DEP-based forward modeling
Authors: Eisen Olaf , Wilhelms Frank , Nixdorf Uwe , Miller Heinrich
Published: November 2007

Isochronic continuous horizons between 20 and 90 m depth in a ground-penetrating radar (GPR) profile, recorded in Dronning Maud Land, Antarctica, are identified by comparison of synthetic and measured single radar traces. The measured radar-gram is derived from a stacked GPR profile; the synthetic radargram is computed by convolution of the complex reflection coefficient profile, based on dielectric profiling (DEP) data of a 150 m ice core, with a depth-invariant wavelet. It reproduces prominent reflections of the measured radargram to a considerable degree. Analyzing matching peaks in both radargrams enables us to identify isochronic reflections and transfer individual volcanic-event datings to the GPR profile. Reflections are primarily caused by changes in permittivity; changes in conductivity are of minor importance. However, several peaks in permittivity andconductivity show a good correlation and indicate that some reflections are related to acidic layers. The results demonstrate the possibility of reproducing radargrams from ice-core property profiles, a necessary step for the interpretation of remotely sensed radar data and the general significance of connecting ice-core and radar data for correct interpretations. Problems related to forward modeling, data gaps, origin of permittivity peaks, and GPR profiles used for comparison, are discussed.
Distribution of basal melting and freezing beneath tributaries of Ice Stream C: implication for the Holocene decay of the West Antarctic ice sheet
Authors: Vogel Stefan W. , Tulaczyk Slawek , Joughin Ian R.
Published: November 2007

Ice-stream tributaries connect the relatively slow-moving interior of the West Antarctic ice sheet (WAIS) with the fast-flowing Siple Coast ice streams. Basal water underneath these ice streams reduces basal resistance and enables the fast motion of the ice. Basal melting being the only source for this water, it is important to include the distribution of basal melting and freezing into numerical models assessing the stability of the WAIS. However, it is very difficult to constrain its distribution from existing field observations. Past borehole observations confirmed the presence of a wet bed at Byrd Station in the WAIS interior and at different locations within Siple Coast ice streams. However, the recent discovery of a 12–25m thick sediment-laden bubble-free basal ice layer at the UpC boreholes indicates that basal freezing is either currently occurring or had occurred upstream during the last glacial–interglacialcycle.We use a flowline model of ice thermodynamics to assess and quantify the spatial and temporal distribution of basal melting and freezing beneath Ice Stream C tributaries, taking into account the geothermal flux, shear heating and heat conduction away from the bed. Under the assumption that the ice was moving over a weak bed (τb =1–10 kPa) our model is able to reproduce a layer of frozen-on ice similar in thickness to the UpC “sticky spot” basal ice layer. Increased basal melting in the early Holocene possibly could have initiated the Holocene decay of the WAIS, whereas increased freezing rates over the past few thousand years could have decreased the amount of basal water in the system, resulting in a strengthening of the bed. This is consistent with current force-budget calculations for ice-stream tributaries and with observed stoppages and slow-downs of ice streams.
Sediment, glaciohydraulic supercooling, and fast glacier flow
Authors: Alley Richard B. , Lawson Daniel E. , Evenson Edward B. , Larson Grahame J.
Published: November 2007

Glaciers often advance over proglacial sediments, which then may enhance basal motion. For glaciers with abundant meltwater, thermodynamic considerations indicate that the sediment–ice contact in the direction of ice flow tends toward an angle opposed to and somewhat steeper than the surface slope (by slightly more than 50%). A simple model based on this hypothesis yields the extent of over-ridden sediments as a function of sediment thickness and strength, a result that may be useful in guiding additional fieldwork for hypothesis testing. Sediment-floored as well as rock-floored overdeepenings are common features along glacier flow paths and are expected based on theories of glacier erosion, entrainment, transport and deposition.
Thermodynamics of basal freeze-on: predicting basal and subglacial signatures of stopped ice streams and interstream ridges
Authors: Christoffersen Poul , Tulaczyk Slawek
Published: November 2007

We have constructed a numerical model that simulates the response of subglacial sediments to basal freeze-on. The model is set up to emulate the basal zone of drilling sites in the Ross Sea sector of the West Antarctic ice sheet. We treat basal freeze-on at an ice–sediment interface as a thermodynamic process that couples the flow of water, heat and solutes in unfrozen subglacial sediments underlying a freezing ice base. The coupling of these flows occurs through the Clapeyron equation, which specifies the dependence of the basal freezing/melting temperature on ice pressure, water pressure, solute concentration and surface tension effects. Thermally driven water flow is induced when an ice base becomes supercooled below the pressure-melting point because ice–water surface tension inhibits ice growth in small pore spaces of fine-grained subglacial sediments. Our model results show that basal freeze-on is capable of inducing considerable changes in the basal zone of both ice streams and interstream ridges. These changes are associated with specific signatures that compare with borehole observations and geophysical surveys. Water-pressure levels are reduced, and thick layers of debris-laden basal ice develop. These basal ice layers and underlying sediments contain a distinct isotopic signal. The predicted stable-isotope ratios reflect Rayleigh-type isotopic fractionation whose significance increases with increasing freezing rates. Supercooling of the ice base induces also measurable changes in the ice-temperature profile of the glacier. Till porosity represents another quantity whose evolution is influenced strongly by basal freeze-on. In particular, measurements of vertical porosity distribution beneath stopped ice streams could be used to back-calculatethe timing of the onset of basal freezing. Our model results show that the basal zone of ice streams and interstream ridges responds sensitively to changes in basal melting/freezing rates. This sensitivity may allow reconstruction of past conditions beneath ice streams and interstream ridges from measurements made on basal ice samples and subglacial sediment samples. Our model results also indicate that meltwater from fast-flowing ice streams may be driven towards the freezing ice base of interstream ridges.
Surge potential and drainage-basin characteristics in East Greenland
Authors: Jiskoot Hester , Murray Tavi , Luckman Adrian
Published: November 2007

We introduce a new glacier inventory of central East Greenland and use the collected data to test proposed theories on surging. The glacier inventory contains 259 glaciers, of which 10 have observed surges and a further 61 are inferred surge-type. The total glaciated area is 5.5×103km2. The inventory was created from a combination of remote-sensing data and maps, and some 24 glacial and geological inventory parameters were collected for each glacier. A multivariate logistic analysis is used to test which combination of glacial and environmental data is conducive to surging behaviour in East Greenland. Three different models suggest that glaciers with a large complexity, low slope and oriented in a broad arc from northeast to south are most likely to be of surge type. Geological conditions, and hence substrate character, appeared not to be related to surge potential. On the basis of these results and the surge dynamics in this region, we suggest a hydrologically controlled surge mechanism operates in central East Greenland.
Interannual surface velocity variations of Pine Island Glacier, West Antarctica
Authors: Rabus Bernhard T. , Lang Oliver
Published: November 2007

The surface velocity of Pine Island Glacier, West Antarctica, during the period 1992–2000 is measured with synthetic aperture radar feature-tracking techniques. Over the observation period, we find a monotonic acceleration with a spatially uniform amplitude of about 12% of the surface velocity. The acceleration extends > 80 km inland of the grounding line into a zone of prominent arcuate crevasses.The upper limit of these crevasses has migrated up-glacier by 0.2 km a−1 correlated with a velocity increase of similar size in the crevassed zone. On the other hand, there is no clear correlation between the velocity variations and observations of grounding-line migration. These findings suggest ongoing dynamic thinning of Pine Island Glacier, providing independent confirmation of recent interferometric results obtained by Rignot and others (2002).
Tidally driven stick–slip motion in the mouth of Whillans Ice Stream, Antarctica
Authors: Bindschadler Robert A. , Vornberger Patricia L. , King Matt A. , Padman Laurie
Published: November 2007

We show that the ice plain in the mouth of Whillans Ice Stream (formerly Ice Stream B), Antarctica, moves by stick–slip motion. During a spring-tide period, rapid motions regularly occur near high tide and during falling tide. This correlation is weaker during a neap-tide period when the tidal magnitudes are less. Precise timing of these motion events suggests that they propagate through the region with a mean velocity of 88 m s−1.We hypothesize that this speed is associated with the propagation of shear waves through a wet subglacial till. Motion events are also seen on more smoothly flowing floating ice. Event delays are very short between grounded and floating stations, suggesting the events propagate through the ice shelf as an elastic wave. We further hypothesize the events are caused by the interaction of a sticky bed, the accumulation of stored elastic strain through the compression of ice by upstream inflow, and tidal forcing. Motion events seem to be triggered either by reduction of vertical normal stresses at high tide or by the increase of shear stresses from sub-shelf ocean currents during falling tide. Event magnitudes are not related to the length of the preceding quiescent period, suggesting significant viscous dissipation within the till.
The distribution and flow characteristics of surge-type glaciers in the Canadian High Arctic
Authors: Copland Luke , Sharp Martin J. , Dowdeswell Julian A.
Published: November 2007

A systematic review of 1959/60 aerial photography, and 1999/2000 Landsat 7 imagery, has identified 51 surge-type polythermal glaciers in the Canadian High Arctic. These were identified from the presence of features such as looped medial moraines, intense folding visible at the surface, rapid terminus advance, heavy surface crevassing, and high surface velocities. These observations suggest that surging glaciers are much more common than previously believed in the Canadian High Arctic, where only six surge-type glaciers have previously been described. Of the 51 surge-type glaciers identified in this study, 15 were observed in the active phase in the 1959/60 and/or 1999/2000 imagery. The most dramatic advances have occurred on western Axel Heiberg Island, where Iceberg,“Good Friday Bay” and Airdrop Glaciers have all advanced by 4–7 km between1959 and 1999. For glaciers with repeat Landsat 7 coverage from 1999 and 2000, image correlation software was used to determine the magnitude and spatial distribution of surge velocities. For example,“Mittie” Glacier on Manson Icefield was moving at a rate of up to 1 kma–over a distance of at least 25 km back from its terminus. The terminus of this glacier has advanced by at least 4 km since 1959, and the glacier was observed to be heavily crevassed during overflights in April 2000, with clear signs of surface lowering of 10–25 m indicated by a strandline.
Dilatant till layer near the onset of streaming flow of Ice Stream C, West Antarctica, determined by AVO (amplitude vs offset) analysis
Authors: Anandakrishnan Sridhar
Published: November 2007

A powerful seismic technique that exploits the phase of the ice-bottom reflections shows that soft till is widespread beneath a West Antarctic ice stream very close to the onset of streaming flow. The amplitude vs offset (AVO) method measures the change in amplitude of the reflection as a function of increasing angle of incidence. For a decrease in acoustic impedance with depth, the reflection phase is negative at low angles of impedance but positive at intermediate angles. The change in phase by 180° is an obvious and robust measure of the relative acoustic impedance contrasts. This technique is only usable when there is a change in phase vs offset, conditions which obtain for “UpB-type” tills (high water pressures and porosity, low compressional- and shear-wave velocities, similar to those observed at Upstream B camp). I have applied this technique to the far upstream regions of Ice Stream C and find that a dilatant ( and presumably deforming), relatively thick (meters) till layer has formed beneath the ice stream within tens of km of the region identified as the transition from inland flow to ice-stream flow. These results suggest that the onset of rapid basal motion is linked to the formation of this deforming subglacial layer.