In the current study, cryoconite samples were collected from six glaciers on the Tibetan Plateau to analyze n-alkanes and n-alkenes. The findings revealed that the concentrations of n-alkanes and n-alkenes varied from 40.1 to 496.1 μg g−1 and from 4.6 to 13.8 μg g−1, respectively. The carbon preference index of the long-chain n-alkanes ranged from 3.3 to 8.4, and the average chain length ranged from 28.7 to 29.3. Moreover, the δ13C of the n-alkanes in cryoconites were within the range of C3 plants, demonstrating that the n-alkanes in cryoconites were only derived from vascular plants. However, the δDmean were more negative than that of C3 plants, which could be caused by dry and humid conditions of glaciers. Unlike n-alkanes, n-alkenes ranged from C17:1 to C30:1 and showed a weak even-over-odd carbon number preference in the Dongkemadi, Yuzhufeng, Laohugou and Tianshan glacier, but a weak odd carbon preference in the Qiyi glacier. The n-alkenes in the YL Snow Mountains showed an obvious odd-over-even carbon number predominance from C17:1 to C22:1 with Cmax at C19:1, and the even-over-odd carbon number preference from C23:1 to C30:1 with Cmax at C28:1. This demonstrated that the n-alkenes of cryoconites may be mainly derived from in situ production in glaciers.

To further elucidate the Late Pleistocene glacial history of mid-elevation mountainous regions in Central Europe, 10Be cosmic-ray exposure (CRE) dating was applied to moraines in the Zastler Tal Valley in the Southern Black Forest. Periods of glacier recession from moraines in this valley began no later than 16 ka, 15 ka, and 13 ka. CRE ages of moraines in this and other parts of the Southern Black Forest cluster around 17–16 ka and 15–14 ka, thus suggesting a common forcing of glacier recession. Equilibrium-line altitudes (ELAs) during moraine formation were calculated for precipitation reconstruction. Observed spatial discrepancies in ELAs at ca. 15–14 ka are explained best by the size of snow-contributing areas. The reconstructed annual precipitation at the ELA for ca. 16 ka and ca. 15 ka is affected by large uncertainties, representing a wide range from ~50% to ~150% of present-day values. Due to various factors, such as drifting snow, the lower bounds of the estimates appear most realistic, thus concurring with the common hypothesis of less precipitation during the last glacial termination than today in Central Europe. Further research is needed before ELAs of small ice masses can be employed for precise precipitation estimates.

Glacier phototroph blooms on the surfaces of ice sheets and glaciers cause albedo reduction, leading to increased melting rates. We observed seasonal changes in the abundance of phototrophs on the Qaanaaq Ice Cap in northwestern Greenland from June to August 2014, and reproduced these changes using numerical and empirical models. The phototroph community on the ice surface mainly consisted of the glacier alga Ancylonema nordenskioldii and the cyanobacterium Phormidesmis priestleyi. The glacier alga appeared on the ice surface in late June, after which its abundance increased exponentially throughout the melting period. A logistic growth model designed for snow algal growth reproduced the measured exponential increases, suggesting that growth could be explained using the model as a function of the ice melting duration. Cyanobacteria appeared and their abundance increased in late July but did not change exponentially thereafter. The abundance of cyanobacteria was explained with an empirical model expressed as a function of the amount of mineral dust on the bare ice surface. Our numerical and empirical models for reproducing glacier algae and cyanobacteria could be useful for quantifying the albedo reduction caused by their growth and the melt rates of the Greenland ice sheet and glaciers in the future.

The extent of glaciers in the western French Alps foreland during the last glacial maximum (LGM, 26.5–19 ka) has not yet been determined, so understanding glacial paleogeography during the LGM remains an open question. This study focuses on the glacial chronology in the western French Alps piedmont using 10Be surface exposure ages on nine glacial boulders and 12 erratic boulders. Results indicate an LGM glacier advance between ca. 24 and 21 ka. During the late LGM, a smaller glacier readvance or stabilization phase occurred at ca. 19 ka, which was followed by a withdrawal phase between ca. 19 and 16.5 ka. Our outcomes show that the LGM extent in the western French Alps was similar or slightly less extensive than the pre-LGM ice extents during the last glacial. Such paleogeography has also been suggested in the western Italian Alps, which share the same accumulation zone with the western French Alps glaciers. The LGM dynamic of the western French Alps foreland glaciers highlighted by our exposure ages is consistent with the timing of the LGM glacier advances and deglaciation with the western Italian ice lobes.

We use PRISM climatic data (1981–2010) and Landsat images (2012–2013) to establish an empirical relationship linking annual temperature and precipitation to the equilibrium line altitude (ELA) of glaciers in the Sierra Nevada (36–41°N, California, USA). For this, we determined the present-day ELAs of 57 glaciers and the local 0°C isotherms elevation Iso0, averaged over the 1981–2010 period. The difference, for each glacier, is Y, the normalized snowline altitude (Y = ELA – Iso0). We then empirically calibrated a logarithmic relationship between this normalized snowline altitude and mean annual precipitation using data from partially covered glaciers. Our calibration is statistically distinct from that previously established for the tropical and midlatitude Andes (Fox and Bloom [1994], Journal of Geography (Chigaku Zasshi), 103, 867–885; Condom et al. [2007], Global and Planetary Change, 59, 189–202). This new relationship for North America is an easy-to-use tool to permit paleoclimatic reconstructions from paleo-ELAs. For a specific paleoglacial site, paleotemperature can be computed knowing the paleoprecipitation range, and vice versa. We also performed a test showing that, if precipitation is well known, the uncertainty associated with paleotemperature is about 1°C (1σ).

We used Sentinel-2 satellite imagery at 10 m resolution to map the extent of Norway's glaciers and ice-marginal lakes over 2018–19. We applied a standardized semi-automated band ratio method to derive glacier outlines and ice-marginal lakes. To optimise the results, we manually edited the ice-lake interfaces, debris, snow and parts of the glaciers situated under shadow. We compared our Sentinel-2 derived outlines with very high-resolution aerial orthophotos and Pléiades satellite orthoimages. Glaciers larger than 0.3 km2 have area differences within 7%, whereas values are larger for smaller glaciers. The orthophotos and orthoimages provide more details and a higher mapping accuracy for individual glaciers, but require manual digitisation, have smaller spatial and temporal coverage and can have adverse snow conditions. We found a total glacier area of 2328 ± 70 km2 of which the ten largest glaciers accounted for 52%. The glacier area decreased 15% since the previous inventory (Landsat data from 1999 to 2006), the reduction being largest in northern Norway (22%) compared to southern Norway (10%). We detected more than 2000 previously undetected smaller glaciers and ice patches (covering 37 km2) and 360 new ice-marginal lakes.

In the Cordillera Darwin, southernmost South America, we used 10Be and 14C dating, dendrochronology, and historical observations to reconstruct the glacial history of the Dalla Vedova valley from deglacial time to the present. After deglacial recession into northeastern Darwin and Dalla Vedova, by ~16 ka, evidence indicates a glacial advance at ~13 ka coeval with the Antarctic Cold Reversal. The next robustly dated glacial expansion occurred at 870 ± 60 calendar yr ago (approximately AD 1150), followed by less-extensive dendrochronologically constrained advances from shortly before AD 1836 to the mid-twentieth century. Our record is consistent with most studies within the Cordillera Darwin that show that the Holocene glacial maximum occurred during the last millennium. This pattern contrasts with the extensive early- and mid-Holocene glacier expansions farther north in Patagonia; furthermore, an advance at 870 ± 60 yr ago may suggest out-of-phase glacial advances occurred within the Cordillera Darwin relative to Patagonia. We speculate that a southward shift of westerlies and associated climate regimes toward the southernmost tip of the continent, about 900–800 yr ago, provides a mechanism by which some glaciers advanced in the Cordillera Darwin during what is generally considered a warm and dry period to the north in Patagonia.

The long-term refrigerated storage of melted snow and/or ice samples for analyses of insoluble microparticles (hereafter, microparticles) may be limited by increases in the biological particle concentration caused by microbial growth after ~1–2 weeks. In this study, we examined an ultraviolet (UV) disinfection method for the storage of melted snow and/or ice samples and determined the effects of this method on microparticles. Surface snow obtained from Glacier No. 31 in the Suntar-Khayata Range, eastern Siberia, Russia was divided into two portions for UV treatment and untreated controls. Microparticle concentrations and size distributions (in the range of 0.52–12.0 μm) in the samples were measured using a Coulter counter. Whereas the microparticle concentration in untreated samples increased, no obvious increase was observed over 53 d in the samples subjected to UV treatment. Microbial growth was detected in only untreated samples using a viable particle counter. In addition, the original microparticle concentrations and size distributions were unaffected by UV treatment. Our results demonstrated that the microparticle size distribution in untreated melted water samples reflects the growth, decomposition and succession of microorganisms over time and further indicate that UV irradiation is effective for long-term storage for microparticle analysis.

Surface albedo typically dominates the mass balance of mountain glaciers, though long-term trends and patterns of glacier albedo are seldom explored. We calculated broadband shortwave albedo for glaciers in the central Chilean Andes (33–34°S) using end-of-summer Landsat scenes between 1986 and 2020. We found a high inter-annual variability of glacier-wide albedo that is largely a function of the glacier fractional snow-covered area and the total precipitation of the preceding hydrological year (up to 69% of the inter-annual variance explained). Under the 2010–2020 ‘Mega Drought’ period, the mean albedo, regionally averaged ranging from ~0.25–0.5, decreased by −0.05 on average relative to 1986–2009, with the greatest reduction occurring 3500–5000 m a.s.l. In 2020, differences relative to 1986–2009 were −0.14 on average as a result of near-complete absence of late summer snow cover and the driest hydrological year since the Landsat observation period began (~90% reduction of annual precipitation relative to the 1986–2009 period). We found statistically significant, negative trends in glacier ice albedo of up to −0.03 per decade, a trend that would have serious implications for the future water security of the region, because glacier ice melt acts to buffer streamflow shortages under severe drought conditions.

Climate change affects various components of the polar environment, including lacustrine systems in many regions of the Antarctic Peninsula. However, the recent status of these life-preserving habitats remains poorly known. We performed field geomorphological mapping and limnological characterization of lakes distributed across Devil's Bay in the northern area of Vega Island, James Ross Archipelago, and described some significant changes in their physical and chemical properties during the 2013 summer season. Fifty lakes were described in the area, even though the area is relatively small (~12 km2). Six major lake types were distinguished based on their geomorphological settings and their origins, and their hydrochemical properties were compared with the neighbouring lake districts. The geomorphic characteristics of the lakes range from stable lakes on till plains (21 out of 50 lakes), to kettle lakes on fresh moraines with an unstable shoreline, to very dynamic ice-contact lakes with fluctuating water levels and inundated areas. The ice-contact lake dynamics were observed using ‘Lake Payer’ as an example. At the interface between hydrology and geosciences, this study provides new and original data from a yet largely unexplored region and emphasizes the importance of cross-disciplinary research in polar sciences.

We present the first-ever mass-balance (MB) observation (2014–19), reconstruction (between 1978 and 2019) and sensitivity of debris-free Stok glacier (33.98°N, 77.45°E), Ladakh Region, India. In-situ MB was negative throughout the study period except in 2018/19 when the glacier witnessed a balanced condition. For MB modelling, three periods were considered based on the available data. Period I (1978–87, 1988/89) witnessed a near balance condition (−0.03 ± 0.35 m w.e. a−1) with five positive MB years. Whereas Period II (1998–2002, 2003–09) and III (2011–19) experienced high (−0.9 ± 0.35 m w.e. a−1) and moderate (−0.46 ± 0.35 m w.e. a−1) negative MBs, respectively. Glacier area for these periods was derived from the Corona, Landsat and PlanetScope imageries using a semi-automatic approach. The in-situ and modelled MBs were in good agreement with RMSE of 0.23 m w.e. a−1, R2 = 0.92, P < 0.05. The average mass loss was moderate (−0.47 ± 0.35 m w.e. a−1) over 28 hydrological years between 1978 and 2019. Sensitivity analysis showed that the glacier was more sensitive to summer temperature (−0.32 m w.e. a−1 °C−1) and winter precipitation (0.12 m w.e. a−1 for ± 10%). It was estimated that ~27% increase in precipitation is required on Stok glacier to compensate for the mass loss due to 1°C rise in temperature.

Glacier mass loss in Alaska has implications for global sea level rise, fresh water input into the Gulf of Alaska and terrestrial fresh water resources. We map all glaciers (>4000 km2) on the Kenai Peninsula, south central Alaska, for the years 1986, 1995, 2005 and 2016, using satellite images. Changes in surface elevation and volume are determined by differencing a digital elevation model (DEM) derived from Advanced Spaceborne Thermal Emission and Reflection Radiometer stereo images in 2005 from the Interferometric Synthetic Aperture Radar DEM of 2014. The glacier area shrunk by 543 ± 123 km2 (12 ± 3%) between 1986 and 2016. The region-wide mass-balance rate between 2005 and 2014 was −0.94 ± 0.12 m w.e. a−1 (−3.84 ± 0.50 Gt a−1), which is almost twice as negative than found for earlier periods in previous studies indicating an acceleration in glacier mass loss in this region. Area-averaged mass changes were most negative for lake-terminating glaciers (−1.37 ± 0.13 m w.e. a−1), followed by land-terminating glaciers (−1.02 ± 0.13 m w.e. a−1) and tidewater glaciers (−0.45 ± 0.14 m w.e. a−1). Unambiguous attribution of the observed acceleration in mass loss over the last decades is hampered by the scarcity of observational data, especially at high elevation, and by large interannual variability.

Under the influence of a rapidly warming climate, abrupt changes have been observed along the coast of Greenland. This commentary is based on a Japanese research project initiated in 2012, in which we examined the recent changes in the coastal environment and their impacts on human society in Qaanaaq, a village in northwestern Greenland. Initially, our research sought to quantify the mass loss of glaciers and its interaction with the ocean in the Qaanaaq region. Over the course of the project in collaboration with local communities, we soon realised that the changes in glaciers and the ocean directly impacted the ~600 residents of Qaanaaq. We observed natural disasters triggered by climate change. Environmental changes are also important for local economy and industry because loss of sea ice may lead to growth in transportation, tourism and mineral resource exploration. In order to share the results of our study with the Qaanaaq community, and to gain understanding of local and traditional knowledge, we organised an annual meeting in the village every summer since 2016. Our experience demonstrates the critical importance of performing a long-term multidisciplinary study, including participation of the local communities to understand the changing environment, and to contribute to a sustainable future in Qaanaaq.

Microfossil records from ice archives allow vegetation, fire and land-use activity reconstructions on broad spatial scales. Samples typically contain low microfossil concentrations. Therefore, large ice volumes are often needed for palynology. Hence, it is crucial to extract maximum microfossil numbers through appropriate physical-chemical treatments. We compare six methods covering the main water reduction procedures: evaporation, filtration and centrifugation with snow samples. Adding a known number of Lycopodium marker spores prior to sample treatment and a second marker (Eucalyptus) after laboratory processing allows a quantitative microfossil loss assessment during pollen extraction. We applied the best-performing method (average loss of 22%) to high-alpine firn cores from Colle Gnifetti glacier for validation with a natural archive containing extremely low microfossil concentrations. We conclude that samples processed with different microfossil extraction protocols may give different results for pollen concentrations, percentages and ratios between different pollen types, especially if vesiculate conifer pollen is an important pollen assemblage component. We recommend a new evaporation-based method which delivers the smallest and least variable losses among the tested approaches. Since microfossil losses are inevitable during laboratory procedure, adding markers prior to sample processing is mandatory to achieve reliable microfossil concentration and influx estimates.

Postglacial emergence curves are used to infer mantle rheology, delimit ice extent, and test models of the solid Earth response to changing ice and water loads. Such curves are rarely produced by direct dating of land emergence; rather, most rely on the presence of radiocarbon-datable organic material and inferences made between the age of sedimentary deposits and landforms indicative of former sea level. Here, we demonstrate a new approach, 10Be dating, to determine rates of postglacial land emergence in two different settings. In southern Greenland (Narsarsuaq/Igaliku), we date directly the exposure, as relative sea level fell, of gravel beaches and rocky outcrops allowing determination of rapid, post–Younger Dryas emergence. In western Greenland (Kangerlussuaq), we constrain Holocene isostatic response by dating the sequential stripping of terrace sediment driven by land-surface uplift, relative sea-level fall, and resulting fluvial incision. The technique we employ provides high temporal and elevation resolution important for quantifying rapid emergence immediately after deglaciation and less rapid uplift during the middle Holocene. 10Be-constrained emergence curves can improve knowledge of relative sea-level change by dating land emergence along rocky coasts, at elevations and locations where radiocarbon-datable sediments are not present, and without the lag time needed for organic material to accumulate.

We analyzed a 2-year time series of meteorological data (January 2011–December 2012) from three automatic weather stations on Laohugou glacier No. 12, western Qilian Mountains, China. Air temperature, humidity and incoming radiation were significantly correlated between the three sites, while wind speed and direction were not. In this work, we focus on the effects of clouds on other meteorological parameters and on glacier melt. On an average, ~18% of top-of-atmosphere shortwave radiation was attenuated by the clear-sky atmosphere, and clouds attenuated a further 12%. Most of the time the monthly average increases in net longwave radiation caused by clouds were larger than decreases in net shortwave radiation but there was a tendency to lose energy during the daytime when melting was most intense. Air temperature and wind speed related to turbulent heat flux were found to suppress glacier melt during cloudy periods, while increased water vapor pressure during cloudy days could enhance glacier melt by reducing energy loss by latent heat. From these results, we have increased the physical understanding of the significance of cloud effects on continental glaciers.

We present an exhaustive spatial analysis using the geographic, geometric, and hypsometric characteristics of 742 North Cascades glaciers to evaluate changes in their areal extents over a half-century period. Our results indicate that, contrary to our initial expectations, glacier change throughout the study region cannot be explained readily by correlations in glacier location, size, or shape. Because of the large error attributable to annual variations in glacier area due to snowpack, no statistically reliable change could be detected for 444 glaciers in our study (a slight majority). Of the North Cascades glaciers that do exhibit detectable change, a majority decreased in area, but nevertheless, some were detectably growing. These findings suggest that the integration of weather patterns over time does not neatly translate into correlations with natural variations in the geometry of glaciers. Our statistical analyses of the changes observed indicate that geometric data from a large number of glaciers, as well as a surprisingly large amount of spatial change, are required for a credible statistical detection of glacier-length and area changes over a short (multidecadal) period of time.

Changes in the altitudes of glacier snowlines (ELAs) and the altitudes of montane vegetation belts (VBAs) measure Quaternary climatic change. An accepted ‘correction’ to such changes by deducting the amount of contemporary sea level fall is wrong, since the air displaced by the ice sheets approximately fills the space left by the falling sea level and so there is no overall downward movement of the troposphere. This also causes a reduced cooling at the lowered sea level relative to that at the former inter-glacial sea level, about 1°C at the Las Glacial Maximum, which reduces the discrepancies previously noted by others between terrestrial and marine estimates of sea-level cooling. The change in temperature is indicated by the product of the ELA or VBA lowering and the environmental lapse rate (ELR). Prior estimates of ΔELA (−900 ± 135 m) and ELR (−6° ± 0.1°C km−1) would indicate a cooling of −5.4°C at interglacial sea level and −4.4°C at glacial sea level, although glacial-period ELRs are not known reliably. Established ELA corrections for local epeirogenic uplift or subsidence are appropriate.