Burial diagenesis of shales of the Frio Formation resulted in an increase in the abundance of mixed-layer illite-smectite (I-S), albite and chlorite, and a decrease in the abundance of K-feldspar, illite and kaolinite. Some of the mineralogic trends determined in this study contrast with the results of Hower et al. (1976) and other studies of Frio shales. The differences are due to improvements in laboratory and clay quantification techniques since the time of the earlier research. I-S composition changed from ~20% to ≥80% illite, and mineralogic and chemical reaction of I-S continued throughout burial. Shale diagenesis was an open-system process that required addition of K2O and A12O3, and resulted in loss of SiO2. The amount of SiO2 made available by shale diagenesis is sufficient to be the source of the quartz-overgrowth cements in the associated Frio sandstones. The relationships between I-S diagenesis and fluid flow from shales into sandstones, generation of abnormal formation-water fluid pressure, onset of sandstone diagenesis and distribution of authigenic phases in sandstones indicate that reaction of the I-S in shales is one of the most important components of the sandstone/shale/formation water diagenetic system.

Interferometric synthetic aperture radar (InSAR) data suffer from an elevation bias due to signal penetration into the firn and ice surface, rendering the height information unusable for elevation and mass-change detection. This study estimates the penetration bias in X-band InSAR data to quantify its impact on elevation and mass-change detection and to demonstrate the applicability of TanDEM-X digital elevation models (DEMs) for cryosphere research. To achieve this, a multiple linear regression model is applied to a time series of four TanDEM-X DEMs acquired between 2010 and 2018 over the Sverdrup Glacier basin (SGB), Devon Ice Cap, Canada. The resulting penetration corrected TanDEM-X DEMs agreed to within ±14 cm of spatially and temporally coincident precise in situ kinematic dGPS data (±10 cm RMSE). Additionally, multi-year estimations of mass change for the SGB derived from differencing TanDEM-X DEMs over multi-year periods between 2010 and 2018, showed good agreement with mean deviation of 338 ± 166 mm w.e. with independent measurements of mass change derived from annual in situ surface mass balance over the same time periods. The results show that the penetration bias can vary significantly, leading to random under- and overestimations in the detection of elevation and mass changes.

A number of different types of bentonite deposits formed by hydrothermal alteration and diagenetic processes are to be found in the Ordu area of the Eastern Black Sea region. The Ca- and Na-bentonite deposits are related to Upper Cretaceous tholeitic to calc-alkaline volcanites, predominantly dacite and andesite, and also include rhyodacite with lesser basalt and their pyroclastic equivalents. In the present study, dacite (PR1), perlite (PR2), moderately altered rocks (MPR), and Na- and Ca-bentonites were studied to describe and compare their mineralogical and geochemical properties and their conditions of formation by means of X-ray diffraction, optical microscopy, scanning electron microscopy, and chemical analytical techniques.

Ca-bentonites, except for smectite, contain opal-CT, feldspar, biotite, and rarely pyrite, while Na-bentonites contain smectite and less feldspar, opal-CT, kaolinite, and illite.

Progressive alteration of the PR2 caused depletion in K2O and Na2O and enrichment in MgO and CaO in all of the Ca-bentonite samples. Na2O was depleted in all of the Na-bentonites and in most of the MPR. The medium and heavy rare earth elements (MREE and HREE) show mass gain or mass loss in the Na-bentonites. The HREE show nearly immobile behavior in the Ca-bentonites. The rare earth elements (REE) and transition elements (TRE) mostly gained mass in the Ca-bentonites in contrast to Na-bentonites. Large-ion lithophile elements (LILE) are strongly depleted in all of the bentonites. The LREE, MREE, and HREE were strongly depleted in most of the MPR samples. TiO2, Lu, Tm, and Tb show immobile behavior in all samples.

PR1 exhibits a slightly positive Eu anomaly. Two MPR samples show slightly positive Eu anomalies (1.03, 1.13), and one Na-bentonite sample displays a slightly positive Eu anomaly (1.04). Most of the Nabentonites have weakly negative Eu anomalies, whereas perlite and the Ca-bentonite have a strongly negative Eu anomaly. The PR1, PR2, MPR, and Na-bentonite present a positive Ce anomaly, and the Ca-bentonite shows a moderately negative Ce anomaly. The Ca-montmorillonites are mainly hydrothermal in origin and derived from alteration of volcanoclastic material in situ and/or in the subaerial environment. The Na-montmorillonites formed by alteration and diagenesis of volcanoclastic material in the sedimentary basin.

Ice patches have implications for landscape and ecosystem dynamics in polar deserts, however, the understanding of the driving factors that control their spatio-temporal variability is limited. This study aims to assess the seasonal and long-term evolution of ice patches on Ward Hunt Island (WHI; 83°N, Canadian High Arctic) based on field measurements of surface mass and energy balance. Results show that mass gains of the ice patch systems occur mostly through drifting snow, making them highly linked to the topography as well as the frequency and magnitude of wind events. Summer ablation is primarily driven by net radiation, but the short-term variability in melt rate is driven by sensible heat fluxes. The highest ablation rates occur during the passage of warm fronts that combine strong winds and mild temperatures. Conversely, foggy days reduce fluxes of solar radiation and sensible heat to the snow/ice surface, thereby suppressing ablation. Ice patches are less climate-sensitive than other cryospheric elements due to a feedback between snow accumulation and topography, however, summer ablation is strongly influenced by micrometeorology. Model projections of these factors suggest that conditions will become critical for preserving ice patches at WHI and along the northern coast of Ellesmere Island as early as in the next decades.

A distributed mass-balance model is used over a 10-year period for the re-analysis of a glaciological mass-balance time series obtained from Brewster Glacier, New Zealand. Mass-balance modelling reveals glaciological mass balance has been overestimated, with an average mass loss of −516 mm w.e. a−1 not captured by observations at the end of the ablation season, which represents 35% of the annual mass balance. While the average length of the accumulation season (199 days) remains longer than the ablation season (166 days), melting is not uncommon in the core part of the accumulation season, with 2–32% of total snowfall being melted. Refreezing of meltwater is also important, with 10% of surface and subsurface melt being refrozen in the present climate. Net radiation, driven primarily by net shortwave radiation, is the main contributor to melt energy, with melt variability mainly influenced by the turbulent heat fluxes, net longwave radiation and the heat flux from precipitation in the ablation season. Snowfalls in summer are an important moderator of melt, highlighting the critical role of the ice-albedo feedback and phase of precipitation on seasonal mass balance. A complete homogenisation of the long-term glaciological mass balance for Brewster Glacier is still required.

The increased availability of high-resolution and high-quality digital elevation models (DEMs) allows for the investigation of small-scale glaciological changes and improved precision in geodetic mass-balance estimates. However, high precision and careful methodological choices are required to retrieve glacier-wide mass changes at annual to seasonal time scales. Here, we used a 7-year time series of 12 DEMs of the Glacier d'Argentière, in the French Alps, derived from the Pléiades optical satellites to assess the ability of sub-metre stereoscopic satellite images to retrieve annual-to-seasonal mass balances. We found good agreement between the five annual and the four winter mass-balance values estimated using a geodetic method and those of in situ glaciological measurements: mean values via the geodetic method are −0.66 m w.e. and 1.47 m w.e. for annual and winter balances, respectively; mean absolute discrepancies are 0.25 m w.e. (annual) and 0.36 m w.e. (winter). Our study identified three main limitations of this methodology: (i) the intrinsic DEM precision; (ii) the lack of control over the satellite acquisition dates; and (iii) the density assumption. The consistency between the methods demonstrates the potential of short time-scale glacier mass-balance monitoring using very high-resolution satellite images.

Under global warming, many glaciers worldwide are receding. However, recent studies have suggested the extension of the Karakoram Anomaly, a region of anomalous glacier mass gain, into the western Kunlun and eastern Pamir mountains. However, the eastern limit of this anomaly in the Kunlun Mountains is unclear. This study, using changes in glacier area and surface elevation, estimates the eastern limit of the Kunlun-Pamir-Karakoram anomaly at ~85°E. Over the past 50 years, glaciers west of 85°E in the Kunlun Mountains decreased in area from 8401 to 7945 km2 at a rate of −0.12 ± 0.07% a−1, showed a reduction in the rate of retreat through time and have recently gained mass, with surface elevation changes of 0.15 ± 0.35 m a−1 over the period of 2000–2013. Glaciers east of 85°E have experienced greater rates of area change (−61 ± 12 km2 and −0.43 ± 0.13% a−1) over the past 50 years, accelerated area loss in recent years and elevation change rate of −0.51 ± 0.18 m a−1 between 2000 and 2013. These patterns of elevation and area change are consistent with regional increases in summer temperature in the eastern Kunlun Mountains and slight cooling in the western Kunlun Mountains.

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.

The Neogene volcanic region of Cabo de Gata, Almería, SE Spain, is dotted with many outcrops of bentonite, some of them of significant economic interest. The bentonites have their origin in the hydrothermal alteration of pyroclastic rocks (15-7 Ma). The deposits are usually associated with fractures. The major mineral is a dioctahedral Fe- and Mg-smectite (89-75%) and this is accompanied by minor amounts of feldspars, quartz, amphiboles, pyroxenes, biotite, zeolites, disordered tridymite, calcite, etc. This paper describes the geological background, the general characteristics of the bentonites and major aspects of their formation, e.g. type of low-temperature hydrothermal solutions, mass balance, chemical evolutions of the smectites and geochemistry of trace elements. Finally, the characteristics of three of the most important deposits are described.

Rock weathering has been investigated from atomic to global scales through the different but complementary approaches of mineralogy, petrography, geomorphology and geochemistry. The sequences of mineral reactions involved in the alteration process are now well known. They explain the global trend of weathering phenomena but do not account for the actual rock transformation dynamics. In particular, they ignore the intimate relation of the mineral reaction progress with the increase in connected porosity. At the hand specimen scale, heterogeneity is the rule: mineral reactions are controlled by local physicochemical conditions. Alteration processes depend largely on the rock microstructure properties. They proceed through nearly-closed, semi- and completely open microsystems which are interconnected by fractures or pores. Before being leached out by the solutions which flow in the large fractures (flux), the soluble elements migrate inside the connected porosity through chemical diffusion. The dissolution of the primary minerals is mediated through local gradients of chemical potential. With increasing alteration, the rock porosity increases, as does the length of the fluid passageways and their constrictivity and tortuosity. Consequently, the apparent diffusion coefficient for the most soluble elements decreases. The amplitude of the chemical potential gradients for the most soluble elements is reduced by the progressive coating of the reactive surfaces by clays and Fe oxyhydroxides. The residence time of these elements inside the weathered rock increases as alteration progresses; an effect enhanced by their temporary adsorption on the exchangeable sites of clays and Fe oxyhydroxides. Consequently, the weathering rate decreases with time. A possible new way to calculate weathering rates could be to measure the residence time of soluble elements inside the different microsystems during their migration towards the diluted solution which occurs in the large fractures.

Evolution of cold dry snow and firn plays important roles in glaciology; however, the physical formulation of a densification law is still an active research topic. We forced eight firn-densification models and one seasonal-snow model in six different experiments by imposing step changes in temperature and accumulation-rate boundary conditions; all of the boundary conditions were chosen to simulate firn densification in cold, dry environments. While the intended application of the participating models varies, they are describing the same physical system and should in principle yield the same solutions. The firn models all produce plausible depth-density profiles, but the model outputs in both steady state and transient modes differ for quantities that are of interest in ice core and altimetry research. These differences demonstrate that firn-densification models are incorrectly or incompletely representing physical processes. We quantitatively characterize the differences among the results from the various models. For example, we find depth-integrated porosity is unlikely to be inferred with confidence from a firn model to better than 2 m in steady state at a specific site with known accumulation rate and temperature. Firn Model Intercomparison Experiment can provide a benchmark of results for future models, provide a basis to quantify model uncertainties and guide future directions of firn-densification modeling.

Tree rings were used to reconstruct mass balance for Peyto Glacier in the Canadian Rocky Mountains from A.D. 1673 to 1994. Summer balance was reconstructed from tree-ring estimates of summer temperature and precipitation in the Canadian Rockies. Winter balance was derived from tree-ring data from sites bordering the Gulf of Alaska and in western British Columbia. The models for winter and summer balance each explain over 40% of the variance in the appropriate mass-balance series. Over the period 1966–1994 the correlation between the reconstructed and measured net balances is 0.71. Strong positive mass balances are reconstructed for 1695–1720 and 1810–1825, when higher winter precipitation coincided with reduced ablation. Periods of reconstructed positive mass balance precede construction of terminal moraines throughout the Canadian Rockies ca. 1700–1725 and 1825–1850. Positive mass balances in the period 1845–1880 also correspond to intervals of glacier readvance. Mass balances were generally negative between 1760 and 1805. From 1673 to 1883 the mean annual net balance was +70 mm water equivalent per year (w.e./yr.), but it averaged −317 mm w.e./yr from 1884 to 1994. This reconstructed mass balance history provides a continuous record of glacier change that appears regionally representative and consistent with moraine and other proxy climate records.

Most Icelandic glaciers show high-accumulation rates during winter and strong surface melting during summer. Although it is difficult to establish and maintain mass-balance programs on these glaciers, mass-balance series do exist for several of the ice caps (Björnsson and others, 2013). We make use of the frequent volcanic eruptions in Iceland, which cause widespread internal tephra layers in the ice caps, to reconstruct the surface mass balance (SMB) in the ablation zone. This method requires information about surface geometry and ice velocity, derived from remote-sensing information. In addition, the emergence angle of the tephra layer needs to be known. As a proof-of-concept, we utilize a prominent tephra layer of the Mýrdalsjökull Ice Cap to infer local SMB estimates in the ablation area back to 1988. Using tephra-layer outcrop locations across the glacier at different points in time it is possible to determine local mass changes (loss and redistribution) for a large part of the ablation zone, without the use of historic elevation models, which often are not available.

This study presents the first reanalysis of a long-term glacier mass-balance record in the Canadian Arctic. The reanalysis is accomplished through comparison of the 1960–2014 glaciological mass-balance record of White Glacier, Axel Heiberg Island, Nunavut, with a geodetically derived mass change over the same period. The corrections applied to homogenize the two datasets, including adjusting for changes in hypsometry over the period of record and the generic differences between methods, are discussed along with the associated systematic and random errors of the two forms of mass-balance measurement. Statistical comparison of the two datasets reveals that within the error margin there is no significant difference between the average annual glaciological balance (–213 ± 28 mm w.e. a−1) and geodetic balance (–178 ± 16 mm w.e. a−1) at White Glacier over the 54 year record. The validity of this result, and the assumptions made in implementing the glaciological method, are critically assessed.

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.

Thermodynamics of a seasonal supraglacial lake were examined based on field data from three summers. At maximum, the lake body consisted of an upper layer with thin ice on top, and a lower layer with slush, hard ice and sediment at the bottom. Sublimation from the upper ice surface averaged to 0.7 mm d−1, and melting in the ice interior averaged to 9.1 mm d−1 during summer. Albedo was on average 0.6 and light attenuation coefficient was ~1 m−1. Averaged over December and January, and over 3 different years, we found that the net solar heating was 137 W m−2, while the losses averaged to 62 W m−2 for the longwave radiation, 16 W m−2 for the sensible heat flux, 24 W m−2 for the latent heat flux and 3 W m−2 for the bottom flux. The depth scale is determined by the light attenuation distance and thermal diffusion coefficient, and the net liquid water volume ranged from 0.5 to 1.0 m in different years. The potential winter growth is more than summer melting, and thus the lake freezes up completely in winter in the present climate.

The study of glaciers and ice caps in remote and cloudy regions remains difficult using current remote sensing tools. Here the potential of stereo radargrammetry (SRG) with RADARSAT-2 Wide Ultra-Fine images is explored for DEM extraction, elevation changes and mass-balance calculations on Barnes Ice Cap (Nunavut, Canada). Over low-relief terrain surrounding Barnes, a vertical precision of ~7 m (1σ confidence level) is measured, as well as an average vertical bias of ~4 m. Moreover, we show that the C-band penetration depth over the ice cap is insignificant at this time of the year (i.e. late ablation season). This is likely due to a wet surface and the presence of superimposed ice that leads to a surface radar response. Comparing the SRG DEMs with other datasets, an historical glacier-wide mass balance of −0.52 ± 0.19 m w.e. a−1 is estimated for 1960–2013, whereas it decreases to −1.06 ± 0.84 m w.e. a−1 between 2005 and 2013. This clear acceleration of mass loss is in agreement with other recent studies. Given its all-weather functionality and its possible use without ground control points, the RADARSAT-2 SRG technology represents an appropriate alternative for glacier monitoring in cloudy and remote regions.