Retreat of Mendenhall Glacier nearJuneau, Alaska, U.S.A., has exposed a bedrock ridge spotted with “siltskins”, patchy coatings of calcite-cemented clay-to sandsized lithic grains. Coatings 0.5–20 mm thick occur in two distinct morphologies. Thin, striated siltskins coat mainly stoss faces. Thicker, corrugated siltskins on lee faces consist of parallel micro-ridges elongated downslope. Thin-section analysis shows that siltskins consist of a basal, calcite-rich layer overlain by microlaminated layers of calcite-cemented lithic grains. Scanning electron microscope (SEM) energy-dispersive spectrometer (EDS) analysis of laminae and surfaces shows laterally persistent Ca/Si differences. Isotopic values for δO18 and δO13 ranged from −19.52‰ to −12.74‰ and −6.18‰ to −3.44‰, respectively, consistent with deposition from subglacial waters of varying isotopic composition and with derivation of carbon from inorganic sources. Corrugated siltskins are complex depositional features modified by erosional processes. Parallel micro-ridges spaced 1–10 mm apparently formed as sediment-rich water dripped down lee-slope rock faces. Ice–rock separation, flow energy and the transported sediment controlled the layering and depositional forms. Siltskins probably formed when a subglacial cavity system was active on the rock ridge and provide clues about how microscale hydrologic processes interact with larger-scale subglacial systems.

Twice-daily satellite observations from the Special Sensor Microwave Imager (SSM/I) indicate melt onset and refreeze on southeast-Alaskan icefields. Melt and refreeze are based on 37 GHz vertically polarized brightness temperatures (Tb) and diurnal-amplitude variations (DAV). Two types of melt regime have different summer characteristics. Onset is characterized by increasing average daily Tb and a switch from low- to high-amplitude DAV. Melt timing, calibrated using Juneau Icefield temperatures, correlates well with nearby stream hydrographs. Some pixels maintain high Tb throughout the melt season and return to low-amplitude DAV after melt onset. Refreeze on these pixels is identified by decrease in Tb and accompanying high-amplitude DAV. Other pixels maintain high DAV throughout the summer, indicating nocturnal refreeze. Fall refreeze is determined by the end of high-amplitude DAV. Interannual variability in melt timing and ablation-season length is high. Melt onset and refreeze timing show a regional tendency toward earlier glacier-melt onset and longer ablation seasons from 1988–98.

We apply a time-dependent distributed glaciohydraulic model to Vatnajökull ice cap, Iceland, aiming to determine the large-scale subglacial drainage structure, the importance of basally derived meltwater, the influence of a permeable glacier bed and Vatnajökull’s discharge contribution to major rivers in Iceland. The model comprises two coupled layers that represent the subglacial horizon perched on a subsurface aquifer in the western sector and bedrock in the eastern sector. To initialize and drive the simulations, we use digital elevation models of the ice surface and bed, the 1999/2000 measured mass balance and an estimate of subglacial geothermal heat fluxes. The modelled subglacial flow field differs substantially from that derived by hydraulic-potential calculations, and the corresponding distribution of basal effective pressure shows a strong correlation between low effective pressure and surge-prone areas in northeastern and southern sectors of Vatnajökull. Simulations suggest that geothermally derived basal melt may account for up to ∼5% of the annual glacial discharge, and buried aquifers may evacuate up to ∼30% of subglacialwater.Time-dependent tests yield estimates of the glacial discharge component in various outlet rivers and suggest a possible seasonal migration of subglacial hydraulic divides. This study of present-day Vatnajökull hydrology forms the starting point for investigations of its future evolution.

Glacier responses to tephra deposition are shown to be highly variable where wind-transported eruption plumes produce narrow distal fallout zones with steep lateral thickness gradients. Significant but short-lived advances of faster-flowing glaciers can be triggered by deposition from modest eruptions. The 1947 eruption of Hekla, south Iceland, covered nearby glaciers with variable thicknesses of tephra, causing dramatic spatial variations in ablation rate. Relative snow ablation increased by an estimated >80% at the eastern margin of Eyjafjallajökull, but decreased by >54% at the western margin. Relative ice ablation increased by only 4% in the east but decreased by >75% at the western margin, only 15 km distant. The effect on mass balance therefore depends on tephra distribution as well as on the nature of the glacier surface. On Gígjökull, retardation of ablationwas greatest at the terminus, associated with an anomalous advance of ≥328 m between 1947 and 1954. Other glaciers in the fallout zone show either no recognizable response, or accelerated retreat. Advances will be short-lived due to the rapid redistribution of supraglacial tephra in a maritime climate, and the potential for discriminating between volcanically and climatically forced advances in the glacio-geomorphological record is low. The long-term glaciological effect of volcanism is to create thin, low-albedo covers after reworking, which make conditions less favourable for glaciation.

Using the finite-element we have modeled the stress field near the calving face of an idealized tidewater glacier under a variety of assumptions about submarine calving-face height, subaerial calving-face height, and ice rheology. These simulations all suggest that a speed maximum should be present at the calving face near the waterline. In experiments without crevassing, the decrease in horizontal velocity above this maximum culminates in a zone of longitudinal compression at the surface somewhat up-glacier from the face. This zone of compression appears to be a consequence of the non-linear rheology of ice. It disappears when a linear rheology is assumed. Explorations of the near-surface stress field indicate that when pervasive crevassing of the surface ice is accounted for in the simulations (by rheological softening), the zone of compressive strain rates does not develop. Variations in the pattern of horizontal velocity with glacier thickness support the contention that calving rates should increase with water depth at the calving face. In addition, the height of the subaerial calving face may have an importance that is not visible in current field data owing to the lack of variation in height of such faces in nature. Glaciers with lower calving faces may not have sufficient tensile stress to calve actively, while tensile stresses in simulated higher faces are sufficiently high that such faces will be unlikely to build in nature.

Based on aerial photographs, topographical maps and the Landsat-5 image data, we have analyzed fluctuations of glaciers in the western Qilian Shan, north-west China, from the Little Ice Age (LIA) to 1990. The areas and volumes of glaciers in the whole considered region decreased 15% and 18%, respectively, from the LIA maximum to 1956. This trend of glacier shrinkage continued and accelerated between 1956 and 1990. These latest decreases in area and volume were about 10% in 34 years. The recent shrinkage may be due either to a combination of higher temperatures and lower precipitation during the period 1956–66, or to continuous warming in the high glacierized mountains from 1956 to 1990. As a consequence, glacier runoff from ice wastage between 1956 and 1990 has increased river runoff by 6.2 km3 in the four river basins under consideration. Besides, the equilibrium-line altitude (ELA) rise estimated from the mean terminus retreat of small glaciers <1 km long is 46 m, which corresponds to a 0.3°C increase of mean temperatures in warm seasons from the LIA to the 1950s.

The transition of inland to complex ice flow is investigated in the downstream area of Shirase Glacier catchment, East Antarctica. The 250 km long Mizuho Triangulation Chain (MTC; Naruse, 1978), 230 km upstream of the grounding line of Shirase Glacier, is re-analyzed by a force-balance analysis, in combination with balance-velocity estimates. Numerical model experiments were carried out with a two-dimensional thermomechanical higher-order ice-sheet model along a central flowline extending from the ice divide (Dome Fuji) to the grounding line of Shirase Glacier, passing through the MTC. Results show that complex flow originates near the MTC, where longitudinal normal stresses and normal drag in particular play a decisive role in the force balance. The fast-flow area, however, is geologically controlled and confined to the “bottleneck” or ice-flow convergence area, a few tens of kilometres upstream of the grounding line. In between the MTC and the grounding zone, longitudinal normal stress and basal sliding are dominant features of the ice-flow regime.

Discharge measurements, climate observations and hydrochemical samples gathered monthly (1998/99) in the Yanamarey and Uruashraju glacier-fed catchments of the Cordillera Blanca, Peru, permit an analysis of the glacier meltwater contribution to stream-flow. These glacier catchments feed the Río Santa, which discharges into the Pacific Ocean. Based on a water-balance computation, glacier melt contributes an estimated 35% of the average discharge from the catchments. For comparison, a volumetric end-member mixing model of oxygen isotopes shows glacier melt contributes 30–45% to the total annual discharge. Based on stream geochemistry, discharge from the Yanamarey glacier catchment provides 30% of the annual volume discharged from the Querococha watershed, which is <10% glacierized. By analogy, the larger Río Santa watershed, also <10% glacierized, receives at least 12% of its annual discharge from melting glacier ice. Tributary watersheds to the Río Santa with larger fractions of glacier cover have less variable runoff and enhanced discharge, demonstrating that the glaciers effectively buffer stream discharge seasonally. With continued glacier melting, stream-flow will likely become more variable, and there will be less dry-season runoff.

There is no discernible river discharge from Austre Okstindbreen, Norway, in winter; any water formed by basal melting is likely to be stored in isolated cavities or sediments at the bed. In summer, a baseflow component of discharge, relatively depleted of 18O, is diluted by water more enriched in 18O. Glacier ice meltwater with a high δ18O value passes rapidly through channelized systems, and emerges in the glacier river within a few hours. Headward extension of the systems in late summer may tap water stored within the glacier since the previous summer’s close-down. Stored water also is released during early melt-season events, when the low-capacity drainage systems cannot accommodate inputs. The high δ18O value indicates that the probable source of the released water is glacier ice or low-altitude snowfall. Stratigraphic variations of δ18O values, resulting from changing weather conditions, characterize the accumulating winter snowpack. During the melt season, the mean δ18O value of the residual pack rises. The water leaving it, depleted of 18O, passes slowly through the glacier, smoothing out variations of supply. δ18O values of rain vary between and within events, and their effect on glacier-river δ18O values is unpredictable.

We present an analysis of many European Remote-sensing Satellite (ERS) synthetic aperture radar (SAR) images of Vatnajökull, Iceland, by comparing them with Advanced Very High Resolution Radiometer (AVHRR) images, mass-balance observations and modelled firn stratigraphy. Summer SAR and AVHRR images both detect the surface firn line (or the transient snow-line when it lies below the firn line) as a distinct boundary.Winter SAR images of Vatnajökull generally display the late-summer surface firn line, and often also a transition at a higher elevation that possibly reflects a subsurface continuation of the firn line. On some winter images the firn line is not adequately detected, possibly due to a shallow firn pack. We found no differences in back-scatter between melting snow and firn, and for most years no differences in reflectance either. For Vatnajökull, equilibrium-line altitude is therefore not a particularly useful estimator of the mean specific mass balance (Bm). Unlike albedo images, SAR images do not display interannual variations of the signal within the accumulation area that are clearly related to Bm. Hence, SAR images of Vatnajökull contain less information about Bm than albedo images. However, for some drainage basins, mean transient snow-line/firn-line altitude during the melting season can be used to infer Bm.

Marine-geological and -geophysical data collected from the continental shelf in Pine Island Bay, Antarctica, reveal a complex paleo-subglacial drainage system controlled by bedrock topography and subglacial meltwater discharge. Significant amounts of freely flowing meltwater existed beneath former ice sheets in Pine Island Bay. Subglacial drainage is characterized by descriptions of glacial landforms imaged on the sea floor and sedimentary deposits collected in piston cores. Bedrock geology is characterized using seismic data. Large-scale landforms on the shelf include channels and cavities incised into impermeable crystalline bedrock. There is a transition from randomly oriented channels on the inner shelf to a dendritic pattern of elongate channels on the middle shelf. On the outer shelf, a change in basal conditions occurs where sedimentary deposits bury crystalline bedrock. No evidence for flowing meltwater exists on sedimentary substrates. Instead, meltwater formed at the ice–sediment contact was incorporated into the sediments, contributing to development of a deforming bed, which was sampled in piston cores. Characterization of subglacial meltwater processes that occurred in the past may aid in understanding the role meltwater plays in stability of the West Antarctic ice sheet today.

Knowledge of iceberg calving is important for understanding instabilities of tidewater glaciers and ice sheets. Since 1995 the terminus of LeConte Glacier, Alaska, U.S.A., has retreated about 2 km and the glacier has thinned approximately 120 m at its 1999 terminus position. Our focus is short-term (hours to weeks) variability of the frequency and magnitude of calving events and calving flux. Both photogrammetric and visual observations are employed in a temporal analysis over a several-week period. We combined these data with measurements of ice speed, tide level, surface water input and water-storage estimates in an attempt to better understand the calving process. Contrary to results obtained over longer time-scales on other glaciers, our results show no correlation between ice speed and the frequency of calving. However, calving events do not appear to occur randomly; often they are a response to measurable changes in other parameters within the terminus region. Caclving can often be attributed to buoyancy perturbations and possibly flexure of the nearly floating terminus. Given the multiple possibilities for buoyancy perturbations, we have found no simple relationship between any specific forcing parameter and calving at short time-scales.

The structural and stable-oxygen-isotopic properties of sea-ice samples collected on land-fast sea ice in the northern Baltic Sea, at 60.2–65.7° N, 21.3–26.9° E, were used to estimate the contribution of snow ice to the total sea-ice thickness, and the fraction of snow in the snow-ice layers and in the total sea-ice thickness. The textural analysis of the sea-ice cores revealed a general structure of a granular surface layer and a columnar basal layer. Using different isotopic criteria to distinguish snow-ice layers, we estimate that snow-ice layers contributed on average 24–32% of the total ice thickness in the study region. Transition ice (intermediate granular–columnar ice) also contributed significantly to the total ice thickness. The snow in the snow-ice layers contributed on average 18–21% of the total sea-ice thickness (by mass), a significant contribution when compared to other regions with sea-ice cover. The methods used could be improved to give more accurate estimates, especially by improving knowledge of the temporal variability of the isotopic properties of the snow incorporated into the sea-ice cover. However, our results can be regarded as plausible estimates for the contribution of snow ice and snow to sea-ice growth in the study area.

Glaciar Coronas is one of the small cirque glaciers located in the Maladeta Range, central Spanish Pyrenees. The ice extends from 3220 to 3065 m a.s.l. Using early and modern photographs and aerial photographs, we study the evolution of the glacier during the 20th century, with particular focus on evaluating the variations in area, ice thickness and the withdrawal of the ice front. This information is compared with some basic climatic data (quarterly and annual evolution of temperatures and precipitation) recorded in the nearby Seira observatory. The results show that ice wastage, clearly associated with the main regional climatic evolution, has been dominant during the 20th century at Glaciar Coronas, but the rate has varied and has not been continuous. Three periods have been detected: (1) 1920s–early 1960s: ice depletion period; (2) early 1960s–mid-1980s: stabilization of the deglaciation process; (3) mid-1980s–1998: marked ice depletion period. The study shows a decrease in area of Glaciar Coronas from 19 ha in 1928 to 2.4 ha in 1998. The decrease in length from 1928 to 1998 has been 325 m.

Using a slightly modified form of the Spring–Hutter equations, glacial outburst floods are simulated from three classic sites, “Hazard Lake”, Yukon, Canada, Summit Lake, British Columbia, Canada, and Grímsvötn, Iceland, in order to calibrate the hydraulic roughness associated with subglacial conduits. Previous work has suggested that the Manning roughness of the conduits is remarkably high, but the new calibration yields substantially lower values that are representative of those for natural streams and rivers. The discrepancy can be traced to a poor assumption about the effectiveness of heat transfer at the conduit walls. The simulations reveal behaviour that cannot be inferred from simplified theories: (1) During flood onset, water pressure over much of the conduit can exceed the confining pressure of surrounding ice. (2) Local values of fluid potential gradient can differ substantially from the value averaged over the length of the conduit, contradicting an assumption of simple theories. (3) As the flood progresses, the location of flow constrictions that effectively control the flood magnitude can jump rapidly over large distances. (4) Predicted water temperature at the conduit outlet exceeds that suggested by measurements of exit water temperature.

This paper describes the use of a neutron probe to measure detailed stratigraphy in ice and snow. The Wallingford neutron probe, developed for measurement of soil moisture, consists of an annular radioactive source of fast neutrons around the centre of a cylindrical detector for slow (thermal) neutrons. In snow and ice, the fast neutrons lose energy by scattering from hydrogen atoms, and the number of slow neutrons arriving at the detector (the count rate) is related to the density of the medium. Calibration equations for count rate as a function of snow density and borehole diameter have been derived. Snow-density profiles from boreholes obtained using the probe show that, despite the smoothing produced by the neutron-scattering process, annual variations in density can be resolved. The potential contribution of the neutron probe to improvements in mass-balance monitoring is discussed.

The response of Arctic ice masses to climate change is studied using ice cores containing information on past climatic and environmental features. Interpretation of this information requires accurate chronological data. Absolute dating of ice cores from sub-polar Arctic glaciers is possible using well-known radioactive layers deposited by atmospheric nuclear tests (maximum fallout in 1963) and the Chernobyl accident (1986). Analysis of several isotopes (3H, 137Cs) shows that 3H provides the most accurate dating of the 1963 maximum, as indicated also in comparison with results from total-beta measurements (90Sr and 137Cs). Mean annual net mass balances are derived from the dated ice cores from 1963 up to the date of the drillings. The 137Cs and 3H deposited by nuclear tests, after decay correction, are used to define a melt index for all 13 ice cores studied. The relative strength of melting and percolation post-depositional processes is studied on the basis of these 137Cs and 3H deposits.

A simple model uses once-daily meteorological values in the U.S. National Centers for Environmental Prediction and U.S. National Center for Atmospheric Research (NCEP–NCAR) re-analysis database to estimate summer balance of South Cascade Glacier, Washington, U.S.A., each year over 1959–99. The rms error, 0.30 m w.e. (r2 = 0.71), is comparable to measurement error. The model relates summer balance linearly to temperature T > 0°C at 2000 m andto snow flux at 1650 m, the altitudes in recent years of the equilibrium line and terminus. The snow flux is the product of the humidity and westerly wind component at 850 hPa when temperature T <+2°C at 1650 m. Temperatures are interpolated linearly between the 850 and 700 hPa levels. Both the positive 2000 m temperature and the snow flux are summed from 26 April to 4 October. When the summer estimates are combined with those from a winter balance model using the same database, the rms error in estimating net balance is 0.40 m (r2 = 0.81). The indicated sensitivities of balance to warming of 1°C are −0.51 m for summer and −0.24 m for winter. On the assumption that the total −0.75 m °C−1 sensitivity exists at all altitudes, a warming of only 0.7°C would be sufficient to overcome the 1986–98 average net balance +0.5 m at the top of the glacier.

We assess the accuracy and precision of the U.S. National Science Foundation’s Support Office for Aerogeophysical Research (SOAR) laser profiling system for mapping topography and detecting surface elevation changes of West Antarctic ice streams. The procedures used to process, calibrate and validate the laser, navigation and global positioning system (GPS) data are presented. The primary objective is to produce surface elevations with the best possible resolution. Repeat surveys of a grid of lines over Whillans Ice Stream and Ice Streams C and E were conducted in the 1997/98 and 1999/2000 seasons. The procedure has been calibrated using special test flights conducted over areas that have been surveyed with precise geodetic GPS equipment mounted on snow-mobiles. After calibration, agreement between the two surfaces is ±10 cm rms. The accuracy and precision of the procedure have been evaluated at points where laser flight-lines cross over one another. The accuracy of the system is found to range from 0.09 to 0.22 m.