Rain-on-snow events can cause wet snow avalanches. Laboratory experiments were carried out to investigate the change in snow strength with increasing water content through rainwater percolation. Snowpack was artificially prepared consisting of a thin ice layer and fine compacted snow, and rainfall (2mmh–1) was artificially applied 22–25.5 and 49–52 hours after the snowpack was formed. Snow hardness was measured with a push–pull force gauge to indicate the snow strength before and after each rain-on-snow event. After the first rainfall, the upper half of the snowpack became wet and a rapid decrease in snow hardness was observed. After the second rainfall the rainwater penetrated the ice layer, high water content was observed above the ice layer but the hardness exceeded that estimated from an empirical relationship between hardness and water content. Micrographs of the snow particles suggest that the delay in grain coarsening observed near the wetting front induces the harder than estimated snow condition.

This study investigates the impact of sea-ice conditions on the production of iceberg plumes from two tidewater glaciers on Devon Ice Cap, Nunavut, Canada. These effects are quantified using a 12 year RADARSAT-1 satellite record from 1997–2008 that contains imagery from approximately every 1–2 weeks in the winter and every 1–4 days in the summer. Iceberg plumes identified in this record are verified against terrestrial time-lapse photography of Belcher Glacier from 2007–08. Results suggest a strong relationship between iceberg plumes and the retreat of sea ice from the glacier termini, with the plumes caused by both the release of previously calved icebergs (ice melange) and new glacier calving. Iceberg plumes are also sometimes observed at other times in the summer and in midwinter (occasionally on both glaciers simultaneously), with these events likely due to new glacier calving alone. Analysis of tides and air temperatures suggests that they provide a minor influence on the timing of iceberg plumes. Instead, it appears that changes in the presence of sea ice are dominant on seasonal timescales, although internal glacier dynamics likely play a significant role for winter plume events that occur when substantial thicknesses of landfast sea ice are present.

The water retention curve (WRC), which shows the relationship between the volumetric liquid water content, θv, and suction, h, is a fundamental part of the characterization of hydraulic properties. Therefore, the formulation of the WRC as a function of snow characteristics is essential for establishing a model of water movement through the snow cover. In this study, we measured the WRC of several snow samples, which had different characteristics (grain size, bulk dry density and grain type), using a gravity drainage column experiment and then analysed these data using the Van Genuchten soil physics model (VG model). The shape of the WRC depended strongly on both the sample grain size, d, and bulk dry density, ρ. Therefore, we introduced the parameter ρ/d to model the WRC of snow. The relationships between the parameters α and n of the VG model and ρ/d change with grain type. For melt forms, α, which is related to the inverse value of the air-entry suction, increases quickly as ρ/d decreases, whereas n, which is related to the gradient of θv vs h, increases with ρ/d. Conversely, neither of these parameters of the VG model for rounded grains showed obvious dependence on ρ/d. These results suggest that water movement through snow cover can be modelled using grain size, bulk dry density and grain type based on the soil physics model.

We identify a series of basal crevasses along a 31 km transect across the northern sector of the Larsen C ice shelf, Antarctica, using in situ ground-penetrating radar. The basal crevasses propagate from a region of multiple, shallow basal fractures to form widely spaced (0.5–2.0 km) but deeply incised (70–134 m) features. Surface troughs, observed in visible imagery, exist above the basal crevasses as the ice vertically shears to reach hydrostatic equilibrium, while widespread surface crevassing occurs along the crests and on the flanks of the undulations, primarily aligned with the topography. We suggest, based on the location of the surface crevasses and the along-flow evolution of the basal crevasses, that the former are induced by a bending stress created by gradients in hydrostatic forces. Using a linear elastic fracture mechanics model, we investigate the sensitivity of basal crevasse propagation to observed trends of ice-shelf thinning and acceleration. Basal crevasses are large-scale structural weaknesses that can both control meltwater ponding and induce surface crevassing. Together, these features may represent an important mechanism in both past and future ice-shelf disintegration events on the Antarctic Peninsula.

Outlet glaciers draining the Antarctic ice sheet into the Amundsen Sea Embayment (ASE) have accelerated in recent decades, most likely as a result of increased melting of their ice-shelf termini by warm Circumpolar Deep Water (CDW). An ocean model forced with climate reanalysis data shows that, beginning in the early 1990s, an increase in westerly wind stress near the continental shelf edge drove an increase in CDW inflow onto the shelf. The change in local wind stress occurred predominantly in fall and early winter, associated with anomalous high sea-level pressure (SLP) to the north of the ASE and an increase in sea surface temperature (SST) in the central tropical Pacific. The SLP change is associated with geopotential height anomalies in the middle and upper troposphere, characteristic of a stationary Rossby wave response to tropical SST forcing, rather than with changes in the zonally symmetric circulation. Tropical Pacific warming similar to that of the 1990s occurred in the 1940s, and thus is a candidate for initiating the current period of ASE glacier retreat.

This paper considers the possible commercial viability of applying the moving particle semi-implicit (MPS) method to avalanches. The MPS method is a powerful tool for reproducing the flow phenomenon with large-scale surface deformation. In order to apply this method to snow avalanches, we modified the original model to introduce constitutive equations of Bingham fluid, dilatant fluid and the erosion–deposition process. The modified model was applied to some cases and evaluated through comparison with experimental results and observed data.

Radar power returned from the basal interface along a 42 km long profile over an ice-rise promontory and the adjacent Roi Baudouin ice shelf, Dronning Maud Land, East Antarctica, is analyzed to infer spatial variations in basal reflectivity and hence the basal environment. Extracting basal reflectivity from basal returned power requires an englacial attenuation model. We estimate attenuation in two ways: (1) using a temperature-dependent model with input from thermomechanical ice-flow models; and (2) using a radar method that linearly approximates the geometrically corrected returned power with ice thickness. The two methods give different results. We argue that attenuation calculated using a modeled temperature profile is more robust than the widely used radar method, especially in locations where depth-averaged attenuation varies spatially or where the patterns of basal reflectivity correlate with the patterns of the ice thickness.

The results of snow-cover chemical composition research in Lefortovo District, Moscow, Russia, are considered. The technogenic impact on snow-cover contamination was evaluated using snow sampling, chemical analysis of liquid and solid phases of snow and data processing. Snow meltwaters were slightly acidic, low-mineralized and contaminated with oil products. Contrasting anomalies of chloride, sodium, nitrites, nitrates, oil products and 3,4-benzopyrene are revealed and the main causes of their formation analyzed. The technogenic load on the air and snow cover was mainly due to an extensive road network and the impact of industrial enterprise. Ecogeochemical associations of elements, polluting snow cover by means of the dust component of air, formed the specific types of geochemical zones due to industrial and road-transport emissions. The common type of man-made contamination in snow cover is characterized by the following general association of contaminants: Hg W Mo Cr Ni Co Pb Cu Zn Cd Ag Bi Sn Ba Nb.

Summer 2000/01 ASTER and Landsat 7 scenes and semi-automated digitization were used to compile a glacier inventory for local glaciers of the Geikie Plateau region, central East Greenland. Of the 332 glaciers (41 591 km2), 120 are tidewater-terminating and drain 90% of the glacierized area. Differencing the 2000/01 tidewater margins from the 1980s GEUS map database ice polygons (113 glaciers) and from 2004/05 ASTER tidewater margins (78 glaciers) shows a cumulative tidewater terminus width decrease from 196 km to 183 km between the 1980s and 2000s, with a corresponding areal loss of ~31 km2 and an effective length change of –14.3 km. Between 2000/01 and 2004/05, areal loss was 26 km2. Average margin retreat rate increased two- to threefold, from 1.7–2.1km2 a–1 (1980s– 2000) to 3.9–5.7km2 a–1 (2000–05). Advances are negligible, apart from two surges, of which one was previously undetected. Coastal, ‘outer’ fjord-terminating, glaciers have two to three times larger areal and effective length retreat rates than ‘inner’ fjord-terminating glaciers. We investigate possible causes and hypothesize that, in addition to ocean temperature and sea ice, changes in sea fog may affect coastal-terminating more than inner fjord-terminating glaciers.

The results of monitoring the rate of recrystallization of snow cover at sites with different degrees of humidity are presented. Observations were carried out at two sites within the same stratigraphic complex of snow cover in the south of Sakhalin island, Sea of Okhotsk. The basic structural changes in snowpack were traced. We determined the rate of growth of the diameter of the crystals in the snow layer. The structure of the snowpack was studied using E.G. Kolomyts’s technique. On the basis of stratigraphic observations, diagrams of the stratification of the snow cover were created, to determine the rate of evolution of ice crystals in a separate layer from its inception.

An experiment using artificial snow was conducted to clarify the influences of temperature and normal load on temporal variations in the shear strength of snow. Artificial snow consisting of dendritic crystals was allowed to accumulate to ∼60 cm depth for the experiment. The shear strength, temperature and weight of the overlying snow were measured at three different depths in the accumulated snow. For the measurement of shear strength, the shear frame index (SFI) was found using a shear frame by placing weights with different masses on the snow contained within the frame. The measured SFI values were treated with the Mohr–Coulomb failure criterion to find the snow cohesion factor C and the internal friction factor tanϕ. The results highlighted similar trends for SFI and C values by which their rate of increase over time was greater with higher snow temperatures and normal load caused by overlying snow. This indicates that C contributes significantly to increased SFI values. tanϕ decreased over time with higher snow temperatures and increased with lower snow temperatures. In low-temperature conditions, in particular, it is likely that snow crystals are compacted but maintain their dendritic morphology.

We compute the eigenmodes (seiches) of the barotropic and baroclinic hydrodynamic equations for an idealized fjord having length and depth scales similar to those of Ilulissat Icefjord, Greenland, into which Jakobshavn Isbræ (also known as Sermeq Kujalleq) discharges. The purpose of the computation is to determine the fjord’s seiche behavior when forced by iceberg calving, capsize and melange movement. Poorly constrained bathymetry and stratification details are an acknowledged obstacle. We are, nevertheless, able to make general statements about the spectra of external and internal seiches using numerical simulations of ideal one-dimensional channel geometry. Of particular significance in our computation is the role of weakly coupled ice melange, which we idealize as a simple array of 20 icebergs of uniform dimensions equally spaced within the fjord. We find that the presence of these icebergs acts to (1) slow down the propagation of both external and internal seiches and (2) introduce band gaps where energy propagation (group velocity) vanishes. If energy is introduced into the fjord within the period range covered by a band gap, it will remain trapped as an evanescent oscillatory mode near its source, thus contributing to localized energy dissipation and ice/melange fragmentation.

Submarine melting at the ice–ocean interface is a significant term in the mass balance of marine-terminating outlet glaciers. However, obtaining direct measurements of the submarine melt rate, or the ocean heat transport towards the glacier that drives this melting, has been difficult due to the scarcity of observations, as well as the complexity of oceanic flows. Here we present a method that uses synoptic velocity and temperature profiles, but accounts for the dominant mode of velocity variability, to obtain representative heat transport estimates. We apply this method to the Sermilik Fjord–Helheim Glacier system in southeastern Greenland. Using lowered acoustic Doppler current profiler (LADCP) and hydrographic data collected in summer 2009, we find a mean heat transport towards the glacier of 29 × 109W, implying a submarine melt rate at the glacier face of 650 ma–1. The resulting adjusted velocity profile is indicative of a multilayer residual circulation, where the meltwater mixture flows out of the fjord at the surface and at the stratification maximum.

Snowpacks can be described as monogenetic rock, where ice is the basic mineral. They can thus be described in terms of a lithologic complex, that is a cluster of different rocks related by formation process or geological era. The stratigraphic complex of snow cover can be cast as a lithologic complex formed in a certain landscape, passing through stages of sedimentation and diagenesis in similar conditions and forming in the same winter a similar spectrum of stratigraphic columns. The stratigraphic complex possesses similar stratification, structure and texture alongside physical, chemical and mechanical characteristics. In landscapes of the same type located in different regions, the same type of stratigraphic complex is formed. The structural transformation of the snowpack is quantitatively described by coefficients of snowpack recrystallization, secondary stratification and texture. Landscape-indicative properties of snowpacks allow the construction of a unified taxonomic scale of stratigraphic complex with a hierarchy of levels: class, subclass, type, subtype and kind.

Snow cover contains various impurities from natural and anthropogenic sources. These impurities percolate into the soil, groundwater and surface water during snowmelt, significantly influencing their composition. We consider the metals content in snow cover in south Sakhalin, Sea of Okhotsk. Samples were taken in February–April 2011 by the Laboratory of Avalanche and Debris-flow Processes of the Sakhalin Department of the FEGI FEB RAS. Sample analysis was performed by the Laboratory of Analytical Chemistry, FEGI FEB RAS. Overall, the metals present in the highest concentrations were Ca, K, Mg and Na. Other metals have low concentrations at the level of total mineralization of surface water and groundwater. Thus, the main sources of metals in south Sakhalin snow cover are the water surface of the sea and anthropogenic sources.

We investigate the sensitivity of sub-ice-shelf melt rates under Pine Island Ice Shelf, West Antarctica, to changes in the oceanic state using an adjoint ocean model that is capable of representing the flow in sub-ice-shelf cavities. The adjoint code is based on algorithmic differentiation (AD) of the Massachusetts Institute of Technology’s ocean general circulation model (MITgcm). The adjoint model was extended by adding into the AD process the corresponding sub-ice-shelf cavity code, which implements a three-equation thermodynamic melt-rate parameterization to infer heat and freshwater fluxes at the ice-shelf/ocean boundary. The inferred sensitivities reveal dominant timescales of 30–60 days over which the shelf exit is connected to the deep interior via advective processes. They exhibit rich three-dimensional time-evolving patterns that can be understood in terms of a combination of the buoyancy forcing by inflowing water masses, the cavity geometry and the effect of rotation and topography in steering the flow in the presence of prominent features in the bedrock bathymetry. Dominant sensitivity pathways are found over a sill, as well as ‘shadow regions’ of very low sensitivities. To the extent that these transient patterns are robust they carry important information for decision-making in observation deployment and monitoring.

Dry-snow metamorphism was the focus of extensive investigations in Russia over several decades. Large amounts of data were collected and used to construct a conceptual model based on crystallographic laws. The difference in terminology compared with the conventional model made these results less accessible to the wider international snow science community. This paper attempts to explain this terminology and its background so that the reported findings are not overlooked in ongoing investigations of the driving forces and consequences of snow metamorphism.

The mass balance of Antarctica is one of the crucial factors for determining sea-level change in a warming climate. The marginal zones of the continent, namely the ice shelves, are most sensitive to climate change. During the 2009/10 austral summer an extensive glaciological field campaign was carried out on Fimbulisen, an ice shelf in East Antarctica, to investigate its recent surface mass balance. Shallow (10–18 m) firn cores were drilled and accumulation rates and stable-isotope ratios determined. For firn-core dating, two different methods were compared: (1) seasonal variations of stable oxygen isotope ratios (δ18O), and (2) dielectric profiling, including using the volcanic eruptions of Pinatubo, Philippines (1991), and El Chichόn, Mexico (1982), as time markers. The mean annual accumulation for the period 1992–2009 ranges from 298 to 349 mmw.e. a–1. The interannual variability at the drilling sites is >30%. Accumulation rates show a weak decreasing trend during the past 20–30 years, which is statistically significant only for one of the cores. Stable-isotope ratios were compared to the snowfall temperature of Neumayer station. Neither the temperatures nor the δ18O values show any trend for the investigated time period.