On the Kamchatka Peninsula in the Russian Far East, 405 glaciers with an estimated total mass of 49 Gt were reported in the 1970s. These have been retreating at an accelerated rate since the start of the 21st century. Because glacier studies in this region are scarce, ice loss and its influence on sea level rise and regional environments is poorly understood. In this study, we analyzed satellite data to quantify glacier mass change from 2000 to 2016 in six major glacier-covered regions on the peninsula. The mean rate of the glacier mass change over the study period was −0.46 ± 0.01 m w.e. a−1 (total mass change was −4.9 ± 0.1 Gt, −304.2 ± 9.1 Mt a−1), which is slightly lower than other regions in mid-latitude and subarctic zones. The mass loss accelerated from >−0.33 ± 0.02 m w.e. a−1 in the period 2000–2006/2010 to <−1.65 ± 0.12 m w.e. a−1 in 2006/2010–2015/16. The increase in mass loss is attributed to a rise in average decadal summer temperatures observed in the region (+0.68°C from 1987–99 to 2000–13). Moreover, a recent trend in Pacific decadal oscillation suggests future acceleration of mass loss due to a decline in winter precipitation.

A 45.97 m long ice core was recovered in the accumulation area of Glaciar Tyndall (50˚59’05’’ S, 73˚31’12’’W; 1756ma.s.l.), Campo de Hielo Patagόnico Sur (southern Patagonia icefield), during December 1999. the firn core was subjected to visual stratigraphic observation and bulk density measurements in the field, and later to analyses of water isotopes (δ18O, δD), major dissolved ions and snow algal biomass. the drillhole remained dry down to about 43 m depth, where a water-soaked layer appeared. Seasonal cycles were found for δ18O, δD and the D-excess, although the amplitudes of the cycles decreased with depth. Major dissolved ions (Na+, K+, Mg2+, Ca2+, Cl–, SO42–) and algal biomass exhibit rapid decreases in the upper 3 m, probably due to meltwater elution. Annual increments defined by the δ18O and D-excess peaks suggest that the minimum net accumulation rates at this location were 17.8ma–1 in 1997/98–1998/99 and 411.0 ma–1 in 1998/99–1999/2000. These are much higher values than those previously obtained from past ice-core studies in Patagonia, but are of the same order of magnitude as those predicted from various observations in ablation areas of Patagonian glaciers.

The Sea of Okhotsk is the southernmost area in the Northern Hemisphere where seasonal sea ice is produced every year. The formation of sea ice drives thermohaline circulation in the Sea of Okhotsk, and this circulation supports the high productivity in the region. However, recent reports have indicated that sea-ice production in the Sea of Okhotsk is decreasing, raising concern that the decreased sea ice will affect not only circulation but also biological productivity in the sea. To reconstruct climatic changes in the Sea of Okhotsk region, we analyzed an ice core obtained from Ichinskaya Sopka (Mount Ichinsky), Kamchatka. We assumed that the remarkable negative peaks of δD in the ice core were caused by expansion of sea ice in the Sea of Okhotsk. Melt feature percentage (MFP), which indicates summer snowmelt, showed high values in the 1950–60s and the mid-1990s–2000s. The high MFP in the 1950–60s was assumed to be caused by an increase in cyclone activity reaching Kamchatka during a negative period of the Pacific Decadal Oscillation index, and that in the 1990–2000s may reflect the increase in solar irradiation during a positive period of the summer Arctic Oscillation index.

The Gorshkov crater glacier at Ushkovsky volcano, Kamchatka, is characterized by a large aspect ratio and special thermodynamic conditions at the bedrock caused by a locally enhanced and spatially varying geothermal heat flux. Furthermore, large parts of this glacier consist of firn rather than pure ice, which alters the rheological properties (such as viscosity and compressibility) of the glacier. We present a newly developed, thermo-mechanically coupled, three-dimensional flow model based on the finite-element (FE) modeling software Elmer, and apply it to the Gorshkov crater glacier. By assuming steady-state conditions, the present-day velocity field, temperature field, basal melting rate and age distribution are simulated. We find that flow velocities are generally small (tens of centimeters per year). Horizontal and vertical velocities are of comparable magnitude, which shows that the shallow-ice approximation is not applicable. Owing to the spatially variable volcanic heat flux, the thermal regime at the ice base is cold in the deeper parts of the glacier and temperate in the shallower parts. The measured temperature profile and age horizons at the K2 borehole are reproduced quite well, and remaining discrepancies may be attributed to transient (non-steady-state) conditions. Firn compressibility is identified as a crucial element for the modeling approach.

As part of a deep ice-core drilling project, GPS and ice-radar measurements were made in the vicinity of Aurora Peak, Alaska, USA. Surface and bed elevation maps were obtained over an area of 400 m×700 m, which includes the drilling site. Surface flow velocity was measured at 18 locations in the same area. The results showed that the drilling site was situated on a gentle north–south ridge with saddle-shaped surface geometry. The horizontal ice-flow velocity at the drilling site was <0.5ma–1 and the ice generally flowed eastward and westward along the surface slope. The glacier bed slopes to the east and has more complex topography than the ice surface. The drilling site was situated above a bedrock dip, with an ice thickness of 252±10 m. These data constitute important information for estimating the effect of ice flow to the ice core retrieved from the depth of up to 180.17 m.

We recorded the burial times of temperature sensors mounted on a specially constructed tower to determine snow accumulation during individual storms in the summit caldera of Mount Wrangell, Alaska, USA, (62° N, 144° W; 4100 m a.s.l.) during the accumulation year June 2005 to June 2006. The experiment showed most of the accumulation occurred in episodic large storms, and half of the total accumulation was delivered in late summer. The timing of individual events correlated well with storms recorded upwind, at Cordova, the closest Pacific coastal weather station (200 km south-southeast), although the magnitude of events showed only poor correlation. Hence, snow accumulation at Mount Wrangell appears to be a reflection of synoptic-scale regional weather systems. The accumulation at Mount Wrangell’s summit (>2.5 m w.e.) exceeded the precipitation at Cordova. Although the direct relationship between accumulation of individual storms at the summit of Mount Wrangell and precipitation events at Cordova may be unique in the region, it is useful for interpreting ice cores obtained on Mount Wrangell. This is especially the case here because the high rate of accumulation allows high time resolution within the core.

General equations of ice dynamics are re-examined, using scale analysis, in order to derive a simplified thermomechanically coupled model for ice flow and heat transfer in ice caps filling volcanic craters. Relatively large aspect ratios between crater depths and diameters, low surface temperatures and intense volcanic heating are the principal characteristics of such craters. The conventional boundary-layer (shallow-ice) approximation is revised to account for these conditions and, in addition, the variable density of the snow, firn and bubbly ice. Large crater depths and intense bottom melting result in low longitudinal balance velocities, controlled by both shear and longitudinal stresses, and hence small surface slopes. In such situations ice can be assumed to be linearly viscous. A flowline model of the glacier dynamics is developed using this assumption. Explicit predictive formulas for ice-particle trajectories and age–depth relations, thus obtained, suggest that the age of ice at the bottom of glaciers in volcanic craters on Kamchatka Peninsula, Russia, may reach hundreds or thousands of years. Ice cores from these glaciers represent unique climatic and volcanic archives.

Debris-laden basal ice is exposed along an ice cliff near Hamna Glacier, Sôya Coast, East Antarctica. The basal ice is about 6.8 m thick and shows conspicuous stratigraphic features. The upper 5.5 m consists of alternating layers of bubble-free and bubbly ice. δ values of the bubble-free ice layers are enriched by 2.4 ±1.0‰ (standard deviation) for δ18O compared to values of neighboring bubbly ice layers above, and slopes of δ18O vs δD are close to 8. Such layers are suggested to have been formed by refreezing of meltwater in an open system. In contrast, part of the bubbly ice layers shows neutral profiles for stable isotopes, suggesting that these ice masses are undisturbed ice-sheet ice which was not affected by melting and freezing. The massive alternating layers are thus considered to have been formed by folding of refrozen and non-melted layers. The lower 1.3 m consists predominantly of bubble-free massive ice. The profile of co-isotopic values shows a change of about 3.0‰ for δ18O at the interface between bubble-free and bubbly ice. Since the isotopic change occurred over a wider thickness than the upper 5.5 m, the basal ice is suggested to have been formed by refreezing of meltwater on a larger scale than the upper 5.5 m.

A glacier at the summit of Ushkovsky volcano, Kamchatka peninsula, Russia, was studied in order to obtain information about the physical characteristics of a glacier that fills a volcanic crater. The glacier has a gentle surface and a concave basal profile with a maximum measured depth of 240 m at site K2. The annual accumulation rate was 0.54 m a−1 w.e., and the 10 m depth temperature was −15.8°C. A 211.70 m long ice core drilled at K2 indicates that (1) the site is categorized as a percolation zone, (2) the stress field in the glacier changes at 180 m depth from vertical and longitudinal compression with transversal extension, which is divergent flow, to a shear-dominated stress field, and (3) the frequent occurrence of ash layers can be a good tool for dating the ice core. The borehole temperature profiles were considered to be non-stationary, but the linear profile made it possible to estimate the basal temperature and the geothermal heat flux at K2. Assuming constant surface and the basal boundary conditions, we constructed two depth–age relationships at K2. These predicted that the bottom ages of the ice core were about 511 or 603 years.

Based on the field data at Koryto glacier, Kamchatka Peninsula, Russia, we constructed a one-dimensional numerical glacier model which fits the behaviour of the glacier. The analysis of meteorological data from the nearby station suggests that the recent rapid retreat of the glacier since the mid-20th century is likely to be due to a decrease in winter precipitation. Using the geographical data of the glacier terminus variations from 1711 to 1930, we reconstructed the fluctuation in the equilibrium-line altitude by means of the glacier model. With summer temperatures inferred from tree-ring data, the model suggests that the winter precipitation from the mid-19th to the early 20th century was about 10% less than that at present. This trend is close to consistent with ice-core results from the nearby ice cap in the central Kamchatka Peninsula.

The Institute of Low Temperature Science at Hokkaido University conducted ice-core drilling in alpine glaciers in the northern North Pacific region to reconstruct climate change in this region for the past few hundred years. We drilled two ice cores in the Kamchatka Peninsula, Russia. An ice core with a length of 211 m was drilled on a glacier at the summit caldera of Ushkovsky mountain in 1998. A second core, with a length of 115m (until bedrock was reached), was drilled on a glacier at the summit caldera of Ichinsky mountain in 2006. We drilled three further ice cores in Alaska, USA. Two ice cores with lengths of 50 and 212m were drilled on a glacier at the summit caldera of Mount Wrangell in 2003 and 2004. The third ice core was drilled on the ice divide among three glaciers, Black Rapids, Trident and Susitna glaciers, which represent a flat saddle north of Aurora Peak in the Alaska Range. This paper details the field operations and characteristics of the different ice-drilling systems used and the problems encountered.

Lake-ice properties at Lake Abashiri, Hokkaido, Japan, were examined using field observations and a 1-D thermodynamic model to clarify formation processes at mid-latitudes subject to significant snowfall as well as moderate air temperature. At all lake sites examined, the ice comprised two distinct layers: a snow ice (SI) layer on top and a congelation ice (CI) layer below. The SI layer occupied as much as 29–73% of the total ice thickness, a much greater fraction than that reported for lakes at Arctic high latitudes. In the model, the CI growth rate was estimated using the traditional heat budget method, while the SI growth rate was calculated assuming the excessive snowfall from the isostatic balance is converted to SI by a snow compression rate (β) with the surface melting rate added when the surface heat budget becomes positive. By tuning the value of β to the observational results of SI thickness, the model outcome successfully reproduced the observational thicknesses of CI and SI, and the break-up date of the lake. Essentially, the model findings show how snow and its formation into SI reduce, by about half, the seasonal variability of total ice thickness.