The chemical composition of meltwater-draining Himalayan glacierized basins reflects the dominance of carbonic acid in weathering of silicate and carbonate minerals, yet the role of sulfuric acid-mediated reactions in the mineral weathering and ionic release is still unclear. Here, we present a long-term study (1992–2018) of chemical weathering characteristics of a precipitation-dominated glacierized basin (Dokriani glacier) of central Himalaya. By using new and reprocessed datasets of major ions from the glacial/subglacial zones of the glacier, we suggest that two-thirds of the dissolved load of the meltwater derives from sulfuric acid-mediated weathering of minerals and rocks. We observed a clear control of carbonic acid-mediated reactions in the early ablation periods, while sulfuric acid-mediated reactions dominate in peak and late ablation periods. The slopes and intercepts in best-fit regressions of [*Ca2+ + *Mg2+ vs *SO42− and HCO3−] and [HCO3− vs *SO42−] in meltwater were following the stoichiometric parameters of sulfide oxidation coupled to carbonate dissolution reactions. The glaciers of the central and western Himalaya are in good agreement with the present estimates. We contend that the bedrock lithology has limited or second-order effects over the ionic release from Himalayan glaciers and surmise that these patterns are broadly applicable to the other orogenic systems of the world.

Although protected areas (PAs) play an important role in ecosystem conservation and climate change adaptation, no systematic information is available on PA protection of high-elevation freshwater ecosystems (e.g., lakes and watersheds with glaciers), their biodiversity and their ecosystem services in the tropical Andes. We therefore combined a literature review and map analysis of PAs of International Union for Conservation of Nature (IUCN) and national systems of PAs and freshwater ecosystems. We found that seven national parks were created for water resources protection but were not designed for freshwater conservation (i.e., larger watersheds). High-value biodiversity sites have not been protected, and new local PAs were created due to water resource needs. We quantified 31 Ramsar sites and observed that PAs cover 12% of lakes, 31% of glacial lakes and 12% of the total stream length in the tropical Andes. Additionally, 120 watersheds (average area 631 km2) with glaciers and 40% of the total glacier surface area were covered by PAs. Future research into the role of PAs in ecosystem services provision and more detailed freshwater inventories within and around PAs, especially for those dependent on glacier runoff, will fill key knowledge gaps for freshwater conservation and climate change adaptation in the tropical Andes.

The response of glaciers to climate change has major implications for sea-level change and water resources around the globe. Large-scale glacier evolution models are used to project glacier runoff and mass loss, but are constrained by limited observations, which result in models being over-parameterized. Recent systematic geodetic mass-balance observations provide an opportunity to improve the calibration of glacier evolution models. In this study, we develop a calibration scheme for a glacier evolution model using a Bayesian inverse model and geodetic mass-balance observations, which enable us to quantify model parameter uncertainty. The Bayesian model is applied to each glacier in High Mountain Asia using Markov chain Monte Carlo methods. After 10,000 steps, the chains generate a sufficient number of independent samples to estimate the properties of the model parameters from the joint posterior distribution. Their spatial distribution shows a clear orographic effect indicating the resolution of climate data is too coarse to resolve temperature and precipitation at high altitudes. Given the glacier evolution model is over-parameterized, particular attention is given to identifiability and the need for future work to integrate additional observations in order to better constrain the plausible sets of model parameters.

Information about the ice volume stored in glaciers is of high importance for sustainable water management in many arid regions of Central Asia. Several methods to estimate the ice volume exist. However, none of them take the specific characteristics of flat terminus debris-covered glaciers into account. We present a method for deriving spatially-distributed ice thickness for debris-covered dendritic glaciers, which are common not only in Central Tien Shan but also in several other mountain ranges in High Asia. The method relies on automatically generated branch lines, observed surface velocities and surface topographic parameters as basic input. Branch lines were generated using Thiessen polygons and Dijkstra's path algorithm. Ice thicknesses for four debris-covered glaciers – South Inylchek, Kaindy, Tomur and Koxkar glaciers – have been estimated along the branch line network solving the equation of laminar flow. For Koxkar and South Inylchek glaciers, respectively, maximum thicknesses of ~250 and 380 m were estimated. These results differ by ~50 m compared with GPR measurements with an uncertainty for the debris-covered parts of ~40%. Based on geodetic mass balances, we estimate that the investigated glaciers lost between 6 and 28% of their volume from 1975 to the early 2000s.

Ice volume during the last ten 100 ka glacial cycles was driven by solar radiation flux in the Northern Hemisphere. Early minima in solar radiation combined with critical levels of atmospheric CO2 drove initial glacier expansion. Glacial cycles between Marine Isotope Stage (MIS) 24 and MIS 13, whilst at 100 ka periodicity, were irregular in amplitude, and the shift to the largest amplitude 100 ka glacial cycles occurred after MIS 16. Mountain glaciers in the mid-latitudes and Asia reached their maximum extents early in glacial cycles, then retreated as global climate became increasingly arid. In contrast, larger ice masses close to maritime moisture sources continued to build up and dominated global glacial maxima reflected in marine isotope and sea-level records. The effect of this pattern of glaciation on the state of the global atmosphere is evident in dust records from Antarctic ice cores, where pronounced double peaks in dust flux occur in all of the last eight glacial cycles. Glacier growth is strongly modulated by variations in solar radiation, especially in glacial inceptions. This external control accounts for ~50–60% of ice volume change through glacial cycles. Internal global glacier–climate dynamics account for the rest of the change, which is controlled by the geographic distributions of glaciers.

Icebergs calved from tidewater glaciers represent about one third to one half of the freshwater flux from the Greenland ice sheet to the surrounding ocean. Using multiple satellite datasets, we quantify the first fjord-wide distributions of iceberg sizes and characteristics for three fjords with distinct hydrography and geometry: Sermilik Fjord, Rink Isbræ Fjord and Kangerlussuup Sermia Fjord. We estimate average total iceberg volumes in summer in the three fjords to be 6.4 ± 1.5, 1.7 ± 0.40 and 0.16 ± 0.09 km3, respectively. Iceberg properties are influenced by glacier calving style and grounding line depth, with variations in size distribution represented by exponents of power law distributions that are −1.95 ± 0.06, −1.87 ± 0.05 and −1.62 ± 0.04, respectively. The underwater surface area of icebergs exceeds the subsurface area of glacial termini by at least one order of magnitude in all three fjords, underscoring the need to include iceberg melt in fjord freshwater budgets. Indeed, in Sermilik Fjord, we calculate summertime freshwater flux from iceberg melt of 620 m3 s−1 (±140 m3 s−1), similar in magnitude to subglacial discharge. The method developed here can be extended across Greenland to assess relationships between glacier calving, iceberg discharge and freshwater production.

We explore the possibility of building a continuous glacier reconstruction by analyzing the integrated sedimentary response of a large (440 km2) glacierized catchment in western Norway, as recorded in the downstream lake Nerfloen (N61°56’, E6°52’). A multi-proxy numerical analysis demonstrates that it is possible to distinguish a glacier component in the ~ 8000-yr-long record, based on distinct changes in grain size, geochemistry, and magnetic composition. Principal Component Analysis (PCA) reveals a strong common signal in the 15 investigated sedimentary parameters, with the first principal component explaining 77% of the total variability. This signal is interpreted to reflect glacier activity in the upstream catchment, an interpretation that is independently tested through a mineral magnetic provenance analysis of catchment samples. Minimum glacier input is indicated between 6700–5700 cal yr BP, probably reflecting a situation when most glaciers in the catchment had melted away, whereas the highest glacier activity is observed around 600 and 200 cal yr BP. During the local Neoglacial interval (~ 4200 cal yr BP until present), five individual periods of significantly reduced glacier extent are identified at ~ 3400, 3000–2700, 2100–2000, 1700–1500, and ~ 900 cal yr BP.

The sediments present in some areas of the Orco Valley provide indications on climatic variations that occurred during the last 6000 years on the southern slopes of the Alps. In particular, distribution and ages of peat layers help define periods and extent of glacial fluctuation in the last 2200 years. Sampling of soils involved in periglacial processes provided a basis for development of a chronological framework of late Holocene environmental change. The data indicate a trend toward cooler climate in the second half of the Holocene. A strong relationship exists between phases of River Po flooding and expansion/retreat phases of the Swiss glaciers: major glacial advances were coeval with periods of intense flooding of the River Po, whereas the phases of glacial retreat coincided with periods of little flooding of the Po. Only in three cases do relationships between glacier activity and floods show weak correlations; two of the cases relate to the warmest periods in approximately the last 2200 years, while the third is the present period. Paleoclimatic evidence from the study region indicates the relatively warm Roman Period between about 2200 and 1900 cal yr BP appears to better represent modern conditions than does the Medieval Warm Period.

Stagnation of the terminal region of a glacier occurs in response to sufficiently large and rapid climatic warming, so the presence of stagnation deposits provides quantitative information on the climate change that forced retreat. Here we use a simple flow-line glacier model to investigate the relationship between stagnation, climate forcing and aspects of the glacier bed. For climatic warming greater than the threshold to cause stagnation, larger or faster warming events cause longer regions of a glacier to stagnate. Smaller or slower warming episodes, below the threshold for stagnation, cause retreat while active flow persists along the entire glacier length. The threshold for stagnation depends not only on the climatic forcing but also on many other aspects of the glacier, with stagnation favored by factors including a lower mean bed slope with greater roughness. Quantitative determination of the climatic forcing consistent with the occurrence or absence of stagnation deposits requires that these site-specific characteristics be incorporated in modeling.

The buzzsaw hypothesis refers to the potential for glacial and periglacial processes to rapidly denude mountains at and above glacier Equilibrium Line Altitudes (ELAs), irrespective of uplift rates, rock type or pre-existing topography. Here the appropriateness of the buzzsaw metaphor is examined alongside questions of the links between glacial erosion and ELAs, and whether the glacial system can produce low-relief surfaces or limit summit heights. Plateau fragments in mountains on both active orogens and passive margins that have been cited as products of glacial and periglacial buzzsaw erosion instead generally represent dissected remnants of largely inherited, pre-glacial relief. Summit heights may correlate with ELAs but no causal link need be implied as summit erosion rates are low, cirque headwalls may not directly abut summits and on passive margins, cirques are cut into pre-existing mountain topography. Any simple links between ELAs and glacial erosion break down on passive margins due to topographic forcing of ice-sheet growth, and to the km-scale vertical swaths through which ELAs have shifted through the Quaternary. Glaciers destroy rather than create low-relief rock surfaces through the innate tendency for ice flow to be faster, thicker and warmer along valleys. The glacial buzzsaw cuts down.

Large glaciers descended western valleys of the Olympic Mountains six times during the last (Wisconsin) glaciation, terminating in the Pacific coastal lowlands. The glaciers constructed extensive landforms and thick stratigraphic sequences, which commonly contain wood and other organic detritus. The organic material, coupled with stratigraphic data, provides a detailed radiocarbon chronology of late Pleistocene ice-margin fluctuations. The early Wisconsin Lyman Rapids advance, which terminated prior to ca. 54,000 14C yr B.P., represented the most extensive ice cover. Subsequent glacier expansions included the Hoh Oxbow 1 advance, which commenced between ca. 42,000 and 35,000 14C yr B.P.; the Hoh Oxbow 2 advance, ca. 30,800 to 26,300 14C yr B.P.; the Hoh Oxbow 3 advance, ca. 22,000–19,300 14C yr B.P.; the Twin Creeks 1 advance, 19,100–18,300 14C yr B.P.; and the subsequent, undated Twin Creeks 2 advance. The Hoh Oxbow 2 advance represents the greatest ice extent of the last 50,000 yr, with the glacier extending 22 km further downvalley than during the Twin Creeks 1 advance, which is correlative with the global last glacial maximum. Local pollen data indicate intensified summer cooling during successive stadial events. Because ice extent was diminished during colder stadial events, precipitation—not summer temperature—influenced the magnitude of glaciation most strongly. Regional aridity, independently documented by extensive pollen evidence, limited ice extent during the last glacial maximum. The timing of glacier advances suggests causal links with North Atlantic Bond cycles and Heinrich events.

Glaciers in the southern Himalayas advanced in the early Holocene despite an increase in incoming summer solar insolation at the top of the atmosphere. These glacier advances are in contrast to the smaller alpine glaciers in the western and northern regions of Central Asia. Two different glacier mass-balance models are used to reconcile this Holocene glacier history with climate by quantifying the change in equilibrium-line altitudes (ELA) for simulated changes in Holocene climate. Both ELA models clearly show that the lowering of ELAs in the southern Himalayas is largely due to a decrease in summer temperatures, and that an increase in monsoonal precipitation accounts for less than 30% of the total ELA changes. The decrease in summer temperatures is a dynamic response to the changes in solar insolation, resulting in both a decrease in incoming shortwave radiation at the surface due to an increase in cloudiness and an increase in evaporative cooling. In the western and northern zones of Central Asia, both ELA models show a rise in ELAs in response to a general increase in summer temperatures. This increase in temperatures in the more northern regions is a direct radiative response to the increase in summer solar insolation.

The McMurdo Dry Valleys are a cold, dry polar desert and the alpine glaciers therein exhibit small annual and seasonal mass balances, often <±0.06 m w.e. Typically, winter is the accumulation season, but significant snow storms can occur any time of year occasionally making summer the accumulation season. The yearly equilibrium line altitude is poorly correlated with mass balance because the elevation gradient of mass balance on each glacier can change dramatically from year to year. Most likely, winds redistribute the light snowfall disrupting the normal gradient of increasing mass balance with elevation. Reconstructed cumulative mass balance shows that the glaciers have lost <2 m w.e. over the past half century and area changes show minimal retreat. In most cases these changes are less than the uncertainty and the glaciers are considered in equilibrium. Since 2000, however, the glaciers have lost mass despite relatively stable summer air temperatures suggesting a different mechanism in play. Whether this trend is a harbinger of future changes or a temporary excursion is unclear.

A numerical model to compute the dynamics of glaciers is presented. Ice damage due tocracks or crevasses can be taken into account whenever needed. This model allowssimulations of the past and future retreat of glaciers, the calving process or thebreak-off of hanging glaciers. All these phenomena are strongly affected by climatechange.

Ice is assumed to behave as an incompressible fluid with nonlinear viscosity, so that thevelocity and pressure in the ice domain satisfy a nonlinear Stokes problem. The shape ofthe ice domain is defined using the volume fraction of ice, that is one in the ice regionand zero elsewhere. The volume fraction of ice satisfies a transport equation with asource term on the upper ice-air free surface accounting for ice accumulation or melting.If local effects due to ice damage must be taken into account, the damage functionD is introduced, ranging between zero if no damage occurs and one.Then, the ice viscosity μ in the momentum equation must be replaced by(1 − D)μ. The damage function Dsatisfies a transport equation with nonlinear source terms to model cracks formation orhealing.

A splitting scheme allows transport and diffusion phenomena to be decoupled. Two fixedgrids are used. The transport equations are solved on an unstructured grid of small cubiccells, thus allowing numerical diffusion of the volume fraction of ice to be reduced asmuch as possible. The nonlinear Stokes problem is solved on an unstructured mesh oftetrahedrons, larger than the cells, using stabilized finite elements.

Two computations are presented at different time scales. First, the dynamics ofRhonegletscher, Swiss Alps, are investigated in 3D from 2007 to 2100 using severalclimatic scenarios and without considering ice damage. Second, ice damage is taken intoaccount in order to reproduce the calving process of a 2D glacier tongue submerged bywater.

The meteorological characteristics and hydrological response of an extreme warm, and cold summer in the McMurdo Dry Valleys are compared. The driver behind the warmer summer conditions was the occurrence of down-valley winds, which were not present during the colder summer. Occurrence of the summer down-valley winds coincided with lower than typical mean sea level pressure in the Ross Sea region. There was no significant difference in the amount of solar radiation received during the two summers. Compared to the cold summer, glaciological and hydrological response to the warm summer in Taylor Valley included significant glacier mass loss, and 3- to nearly 6000-fold increase in annual streamflow. Lake levels decreased slightly during the cold summer, and increased between 0.54 and 1.01 m during the warm summer, effectively erasing the prior 14 years of lake level lowering in a period of three months. Lake level rise during the warm summer was shown to be strongly associated with and increase in degree days above freezing at higher elevations. We suggest that strong summer down-valley winds may have been responsible for the generation of large glacial lakes during the Last Glacial Maximum when ice core records recorded annual temperatures significantly colder than present.