A quasi-three-dimensional (3-D) climate model (Sellers, 1983) was used to simulate the climate of theLast Glacial Maximum (LGM) in order to provide climatic input for themodelling of the Northern Hemisphere ice sheets. The climate model isbasically a coarse-gridded general circulation (GCM) with simplifieddynamics, and was subject to appropriate boundary conditions for ice-sheetelevation, atmospheric CO2 concentration and orbital parameters.When compared with the present-daysimulation, the simulated climate at theLast Glacial Maximum is characterized by a global annual cooling of 3.5°Cand a reduction in global annualprecipitation of 7.5%, which agrees wellwith results from other, more complex GCMs. Also the patterns of temperaturechange compare fairly with mostother GCM results, except for a smallercooling over the North Atlantic and the larger cooling predicted for thesummer rather than for the winter over Eurasia.

The climate model is able to simulate changes in Northern Hemispheretropospheric circulation, yielding enhanced westerlies in the vicinity ofthe Laurentide and Eurasian ice sheets. However, the simulated precipitationpatterns are less convincing, and show a distinct mean precipitationincrease over the Laurentide ice sheet. Nevertheless, when using themean-monthly fields of LGM minus present-day anomalies of temperature andprecipitation rate to drive a three-dimensional thermomechanical ice-sheetmodel, it was demonstrated that within realistic bounds of the ice-flow andmass-balance parameters, veryreasonable reconstructions of the Last GlacialMaximum ice sheets could be obtained.

The frequencies of both linear and angular oscillations in a vertical plane of a floating iceberg are shown to converge as the horizontal dimensions become relatively larger. For icebergs of thickness of 250 m, the calculated period of about 30 s is supported by actual observations which have been reported. Because swell in the Southern Ocean may extend to such long periods, a critical relationship between iceberg dimensions and ocean wavelengths has been formulated. This relates the conditions under which the length and natural period of an iceberg may, respectively, coincide with the wavelength and wave period of the ocean.

Aerial photography and satellite imagery of the Petersen ice shelf, Nunavut, Canada, from 1959 to 2012 show that it was stable until June 2005, after which a series of major calving events in the summers of 2005, 2008, 2011 and 2012 resulted in the loss of ∼61% of the June 2005 ice-shelf area. This recent series of calving events was initiated by the loss of extensive regions of ˃50-year-old multi-year landfast sea ice from the front of the ice shelf in summer 2005. Each subsequent calving event has been preceded by open-water conditions and resulting loss of pack-ice pressure across the front of the ice shelf, and most occurred during record warm summers. Ground-penetrating radar (GPR) ice thickness measurements and RADARSAT-2 derived observations of surface motion indicate that tributary glaciers provided total ice input of 1.19-5.65 Mta–1 to the ice shelf from 2011 to 2012, far below the mean surface loss rate of 28.45 Mta–1. With recent losses due to calving and little evidence for current basal freeze-on, this suggests that the Petersen ice shelf will no longer exist by the 2040s, or sooner if further major calving events occur.

This study describes the morphological and dynamic changes of Parkachik Glacier, Suru River valley, Ladakh Himalaya, India. We used medium-resolution satellite images; CORONA KH-4, Landsat and Sentinel-2A from 1971–2021, and field surveys between 2015 and 2021. In addition, we used the laminar flow-based Himalayan Glacier Thickness Mapper and provide results for recent margin fluctuations, surface ice velocity, ice thickness, and identified glacier-bed overdeepenings. The results revealed that overall the glacier retreated by −210.5 ± 80 m with an average rate of 4 ± 1 m a−1 between 1971 and 2021. Whereas a field study suggested that the glacier retreat increased to −123 ± 72 m at an average rate of −20 ± 12 m a−1 between 2015 and 2021. Surface ice velocity was estimated using COSI-Corr on the Landsat data. Surface ice velocity in the lower ablation zone was 45 ± 2 m a−1 in 1999–2000 and 32 ± 1 m a−1 in 2020–2021, thus reduced by 28%. Further, the maximum thickness of the glacier is estimated to be ~441 m in the accumulation zone, while for glacier tongue it is ~44 m. The simulation results suggest that if the glacier continues to retreat at a similar rate, three lakes of different dimensions may form in subglacial overdeepenings.

Computations over 50 000 years into steady state with Greve’s polythermal ice-sheet model and its numerical code are performed for the Greenland ice sheet with today’s climatological input (surface temperature and accumulation function) and three values of the geothermal heat flux: (42, 54.6, 29.4) mW m−2. It is shown that through the thermomechanical coupling the geometry as well as the thermal regime, in particular that close to the bed, respond surprisingly strongly to the basal thermal heat input. The most sensitive variable is the basal temperature field, but the maximum height of the summit also varies by more than ±100m. Furthermore, some intercomparison of the model outputs with the real ice sheet is carried out, showing that the model provides reasonable results for the ice-sheet geometry as well as for the englacial temperatures.

The water balance of the mountainous Durance river catchment, French Alps, is simulated from 1981 to 1994 with a soil-vegetation-atmosphere transfer (SVAT) model. Particular attention is paid to the snow-cover evolution using a detailed model of the snowpack evolution. The results are validated by comparison of the simulated discharges calculated by the SVAT with daily observations at three gauging stations located in the watershed. Three different spatial resolutions are used (1, 8 and 46 km) in order to evaluate the impact on the surface-water-budget results. Comparison with the finest resolution indicates the need for sub-grid-scale parameterization for the model with larger resolution.

Supraglacial streams form annually during the melt season, transporting dissolved solutes from the melting ice and snowpack to subglacial flow paths and the glacier terminus. Although nutrient and carbon processing has been documented in other supraglacial environments (cryoconite holes, snowpack), little work has examined the potential for in-stream nutrient retention in supraglacial streams. Here we carried out a solute nutrient injection experiment to quantify NH3 +, PO4 3− and labile dissolved organic carbon (DOC) retention in a supraglacial stream. The experiment was performed on a 100 m stream reach on Mendenhall Glacier, an outlet glacier on the Juneau Icefield, southeastern Alaska, USA. The study stream contained two distinct reaches of equal length. The first reach had a lower velocity (0.04 ms−1) and contained abundant gravel sediment lining the ice–water interface, while the second reach was devoid of bedload sediment and had an order-of-magnitude higher velocity. At the end of the second reach, the stream emptied into a moulin, which is typical of supraglacial streams on this and other temperate glaciers. We found that N and P were transported largely conservatively, although NO3 − increased along the reach, suggestive of nitrification. Labile DOC was retained slightly within the stream, although rates were low relative to the travel times observed within the supraglacial stream. Although our findings show that these streams have low processing rates, measurable in-stream nitrification and dissolved organic matter uptake within this biologically unfavorable environment suggests that supraglacial streams with longer residence times and abundant fine substrate have the potential to modify and retain nutrients during transport to the glacier terminus.

A degree-day model extended for surface mass-balance calculations has been applied to derive the sensitivity of Gran Campo Nevado ice cap (GCN), southwest Patagonia, to climate change. Seasonal sensitivity characteristics were computed using automatic weather station data gathered in the period 2000–05. Results indicate pronounced mass-balance sensitivity to temperature during the summer, with monthly values of –0.27±0.01mw.e. K–1. Monthly sensitivity to a 10% precipitation perturbation fluctuates around +0.03mw.e The sensitivity characteristics obtained were used to model the surface mass-balance evolution of GCN during the 20th and 21 st centuries based on monthly means of air temperature and precipitation derived from bias-corrected weather station data and statistically downscaled re-analysis and general climate model data. Surface mass balance shows a persistently negative trend ranging from around +1mw.e. a–1 at the beginning of the 20th century down to almost –1.5mw.e. a–1 during the first years of the 21st century, with only a few positive years occurring occasionally during the second half of the 20th century. The scenario for the end of the 21 st century totals approximately –4.5mw.e. a–1, i.e. an estimated ice volume loss for GCN of 59 km3 during 1900–2099.

Though solar radiation is important for glacier mass-balance simulation, solar radiation data are not always available. As a result of analyzing meteorological data measured in the Himalaya and the Tibetan Plateau, a favorable correlation between precipitation and atmospheric transmissivity of solar radiation is found in terms of monthly values. Monthly mean solar radiation is derived from the relationship between atmospheric transmissivity of solar radiation and precipitation with input of monthly precipitation, latitude, skyline and time. The differences between estimated and observed monthly mean solar radiation are <40Wm−2 in most cases. However, the differences at some sites are significantly large. The error in the estimated solar radiation during the monsoon season can be large when the monthly mean precipitation rate is about 5 mm d−1. Though the error in the estimated solar radiation during the non-monsoon season is generally small due to low precipitation in the Himalaya and the Tibetan Plateau during this season, it can exceed 100 W m−2.

The major results from a comprehensive study of the Amery Ice Shelf are presented, following the work of a wintering expedition in 1968 and supplemented by further measurements during the summer seasons of 1969 to 1971. The Programme included ice-core drilling, oversnow surveys for ice movement and optical levelling, ice-thickness sounding, and measurements of snow accumulation. The new data obtained provide the basis for a more accurate assessment of the mass balance and dynamics of the ice shelf than was possible from the earlier surveys.

The results indicate a substantial growth of basal ice under the ice shelf inland where the ice thickness is greater than 450 m. Further towards the ice front the high strain thinning is approximately balanced by the horizontal ice advection.

The velocity distribution over the ice shelf is primarily governed by a substantial surface slope towards the ice front and high restraining shear stress along the sides.

Basal channels, which form where buoyant plumes of ocean water and meltwater carve troughs upwards into ice-shelf bases, are widespread on Antarctic ice shelves. The formation of these features modulates ice-shelf basal melt by influencing the flow of buoyant plumes, and influences structural stability through concentration of strain and interactions with fractures. Because of these effects, and because basal channels can change rapidly, on timescales similar to those of ice-shelf evolution, constraining the impacts of basal channels on ice shelves is necessary for predicting future ice-shelf destabilization and retreat. We suggest that future research priorities should include constraining patterns and rates of basal channel change, determining mechanisms and detailed patterns of basal melt, and quantifying the influence that channel-related fractures have on ice-shelf stability.

We describe a derivation of surface velocities and associated errors for Black Rapids Glacier, Alaska, U.S.A., using single-orbital-path synthetic aperture radar interferometry (InSAR). The technique described is adapted to small temperate glaciers with complex flow patterns. We also describe a motion anomaly, apparent in the InSAR phase signal, that persisted on Black Rapids Glacier for at least 78 days during winter 1991/92 and recurred in 1996. This anomaly is interpreted using a basal hydrology hypothesis in which a hydraulic head is maintained at the glacier bed at close to the overburden pressure. This permits a cumulative influx of 1.6 × 106 m3 of water under the glacier, a sort of shallow subglacial lake, that migrates downstream at an average rate of 30 m d− 1 over 78 days. The motion anomaly is speculated to be an unsuccessful bid for surge initiation.

Several years of daily microwave satellite ice drift are combined with moored upward-looking sonar (ULS) ice drafts into an ice-volume flux record at points along a flux gate across the Weddell Sea, Antarctica. Monthly ice transport varies at the mooring locations from a maximum export of 0.4 m3 s−1 near joinville Island to −0.4 m3 s−1 imported along the Fimbul and Riiser-Larsen ice-shelf margins. Winter peaks are observed at each end of the flux gate, where high concentrations of deformed ice are advected in and out of the basin along the coastline. The central gyre, in contrast, exhibits negligible seasonality and much smaller volume transports. During the period of overlapping ULS operation, the mean monthly integrated ice export west of the gyre center is 59 × 103 m3 s−1, and the import in the East Wind Drift is −17 x103 m3s -1. ULS data are compared with ERS satellite observation of radar backscatter to obtain an empirical relationship between ice thickness and the rate of change of backscatter with incidence angle. Resulting proxy ice-thickness data are combined with Special Sensor Microwave/Imager-derived ice velocities to obtain seasonally varying estimates of net ice-volume flux for the period 1992−98. Significant interannual variability is observed in ice-volume flux expressed as fresh-water transport. A maximum annual mean of 0.054 Sv is observed in 1992; with a minimum of 0.015 Sv in 1996. A 6 year mean transport of 0.032 Sv is observed. Maximum seasonal ice export occurs in July 1992, with a minimum in November 1996. The 10 year mean area flux is 30 × 103 m2s –1 Interannual variations in net volume flux closely follow variations in area flux, with summer minima in 1990/91 and 1996/97. Maximum area transport occurs in 1991, and although this predates the ERS-1 scatterometer data, ice-thickness estimates by Harms and others confirm 1991 as a decadal peak in net integrated fresh-water transport.

Moraines in five valleys of the Turkestan range of the Pamir-Alay, Kyrgyz Republic, were analyzed based on their geographical position and elevation, morphology and ecological characteristics. The moraines represent four glacial advances: Turkestan Stage I, during the last glacial; Stage II, during the late Holocene; Stage III, during the Little Ice Age (AD 1400−1900); and Stage IV, during the 20th century. In Turkestan Stage I, the glaciers expanded to about 5−10 km from their present termini during marine isotope stage (MIS) 2. Radiocarbon ages for several soil layers buried in a lateral moraine of the Asan-Usin glacier frontal area indicate that glaciers shrank or stagnated between at least ∼48 and 42 kyr BP, during the warm climate of mid-MIS 3. These findings suggest that the glacier expansion in mid-MIS 3 was relatively small compared to the glacier expansion of MIS 2 during the last glacial.

The first topographic and ice-motion maps of the northwestern flank of Hielo Patagónico Norte (HPN, northern Patagonia Icefield), in Chile, were produced using satellite synthetic-aperture interferometric radar data acquired by NASA’s Spaceborne Imaging Radar C instrument in October 1994. The topographic map has a 10 m vertical precision with a 30 m horizontal spacing, which should be sufficient to serve as a reference for monitoring future mass changes of the icefield. The ice-motion map is accurate to within 4 mm d−1 (or 1 m a−1). The radar-derived surface topography and ice velocity are used to estimate the ice discharge from the accumulation area of four outlet glaciers, and the calving flux and mass balance of Glaciar San Rafael. The results demonstrate the use of SAR interferometry for monitoring glaciological parameters on a spatial and temporal scale unattainable by any other means.

A network of automated time-lapse cameras was deployed on Variegated Glacier, Alaska, to establish the temporal and spatial patterns of velocity change at a one-day time resolution. Results from the summers of 1979, 1980, and 1981 are presented; a surge occurred in 1982 and 1983. The principal velocity variations were pulses of increased speed, lasting about one day and referred to as “early”- or “late”-season motion events. The former recurred quasi-periodically on the upper part of the glacier in the early part of the melt season; the latter occurred later in the summer and were correlated with major storms or melting. Supplemental information about the occurrence of motion events was obtained from monitoring of seismic activity. Evidence for several other types of velocity changes was found.

Sca-icc drift and deformation were measured with an array of drifting buoys during a 1995 winter experiment off the East Antarctic continental shelf south of the Antarctic Divergence. The buoys were configured so that deformation of the icefield could be monitored on a range of spatial scales from 2 to 130 km. The mean hourly drift rate during the 3 week-long experiment was 0.21 m s−1, and the mean daily translation of the field was 17.3 km. Differential kinematic parameters calculated from the data show a very high short-term variance, indicating that high-frequency processes are dominant. Spectral analysis of the velocity data shows a major peak of the energy spectrum at the frequency of passage of synoptic weather systems, and a second peak at the inertial frequency. A major storm event occurred during the experiment. Net divergence over this phase of the experiment, as measured by a five-buoy array, is small compared to the short-period variance. This alternating divergence and convergence has a marked effect on the net ice growth. Intense freezing and rapid new ice formation occurs in the open water areas formed during divergence, and this is thickened by rafting and ridge-building during the subsequent convergence. New open water areas equivalent to 10% of the total area formed during the first phase of the experiment. A one-dimensional multilayer thermodynamic model of ice growth shows that this led to an increase of 2.8 cm in the area-averaged ice growth over a 7 day interval, which is equivalent to 40-50% of the total estimated ice growth over the region.

The lowest glacial cirques in Wales are in the South Wales Coalfield and in western mid-Wales: the highest are in the northeast and on the highest mountains. Floor altitudes show great local variability, but in general rise to the northeast in most of Wales and northward in southern Wales, as does the former glaciation level (from 470–710 m a.s.l.). The pattern is similar for reconstructed Younger Dryas glaciers, which occupied the higher and even some of the lower cirques. If cirque development had spread to lower areas as the ice sheet built up, cirques would be expected in peripheral areas such as the northeast and southwest, away from the main ice-sheet sources. This is not seen and there is no clear evidence of time-transgressive cirque glaciation. Cirques relate to phases of glaciation with an ELA a little below that in the Younger Dryas.

Analyses of air trapped in an ice core from the South Pole indicate that the CO2 concentration may have increased by about 10 ppm and that the 13C/12C ratio decreased slightly in the thirteenth century. These changes, if really of atmospheric origin, must be due to a significant input into the atmosphere of CO2, either of biogenic or of oceanic origin.

18O/16O ratios in CO2 from different ice cores are much lower than those which have been observed in atmospheric carbon dioxide. A possible explanation is that the CO2 has equilibrated isotopically with the ice. We have calculated equilibrium isotope-fractionation factors between ice and carbon dioxide and found that the observed 18O/16O ratios of CO2 are indeed near isotopic equilibrium with the ice. This indicates that an exchange of oxygen atoms probably occurs between ice and included CO2.