Proglacial Area Research Articles

Catchment-scale network analysis of functional sediment connectivity during an extreme rainfall event in the Grastal catchment, Austrian Central Alps

Global warming significantly impacts sediment dynamics in glaciated catchments, affecting water resource operations, water quality, recreational activities, and ecological systems. The propagation of climate-change-induced geomorphic changes and the catchment's sediment yield are moderated by sediment connectivity, defined as the degree to which a geomorphic system facilitates sediment transfer. Quantifying functional sediment connectivity at the catchment scale remains a challenge. To address this, we propose a novel approach combining graph theory with the morphological method. This approach is exemplified through a detailed case study of a 2022 thunderstorm event in the Grastal valley, Tyrol, Austria. First, a graph of potential sediment cascades is constructed using a geomorphological map, a digital elevation model and a flow routing algorithm. A short-term Digital Elevation Model of Difference (DoD) from consecutive ALS surveys is then used to infer sediment fluxes and calculate the Sediment Delivery Ratio (SDR) for each landform. The primary sediment mobilising processes were debris flows and fluvial erosion, with a significant proportion of debris flow material being deposited on slopes, not reaching the fluvial corridor. Strong fluvial erosion was observed in the proglacial area, but the propagation of these geomorphic changes is halted by an alluvial fan and a lake. Most landforms can be clearly categorised as connecting or disconnecting features based on their SDR. In total, a maximum of 12 % of mobilised sediments exited the catchment. Our findings demonstrate that (i) short-term, catchment-wide DoDs are valuable for assessing functional connectivity at an event temporal scale, (ii) using landforms as fundamental spatial units allows for the identification and in-depth analysis of critical sediment sinks and sources, and (iii) graph analysis facilitates the catchment-wide calculation of sediment delivery ratios between meaningful fundamental units and the delineation of significant sediment cascades.

Radiological and elemental composition of cryoconite and glacier mice from Vatnajökull, Iceland

Cryoconite has been demonstrated to be an efficient accumulator of some classes of contaminants on glaciers in both mountain and polar environments, however the accumulation of contaminants in cryoconite in Iceland has received very little attention to date. To understand the spatial variability of natural and anthropogenic fallout radionuclides and metals on glaciers in Iceland, we present the first study of this region including both cryoconite from three glaciers: Virkisjökull; Skaftafellsjökull; and Falljökull, together with moss balls (‘glacier mice’) from Falljökull. The cryoconite samples and glacier mice were analysed using XRF spectrometry to assess their elemental composition and gamma spectrometry to identify, and quantify, fallout radionuclides, primarily 7Be, 137Cs, 241Am, excess 210Pb, and 40K. The results revealed that the cryoconite samples had similar compositions, influenced by local geology and natural sources of volcanic ash and dust. Higher concentrations of radionuclides and heavy metals were found in both cryoconite and glacier mice compared to control samples comprising nearby proglacial sediments. In comparison to other glaciers in the Northern Hemisphere, however, cryoconite from Icelandic glaciers contains some of the lowest activity concentrations of key radionuclides. Consequently, cryoconite deposits that are released and diluted during the melt and retreat of Icelandic glaciers are unlikely to be of environmental concern following transport to proglacial areas.

A hydrogeological conceptual model of aquifers in catchments headed by temperate glaciers

Abstract. For reliable forecasting of the evolution of critical water resources, as well as of potential flood and landslide hazards and their response to climate change, it is necessary to improve the understanding and quantification of unknown aquifer systems in glacierized catchments. We focus on four southeastern outlet glaciers of the main Icelandic ice cap, Vatnajökull. A multidisciplinary approach is carried out, including the acquisition of new in situ data to characterize aquifers and their groundwater dynamics. Moreover, the recharge to aquifers from glacial melt and effective rainfall is estimated. From a detailed analysis of all available data and the determination of the dynamic characteristics of the aquifers, a hydrogeological conceptual model of glacierized catchments is constructed: (i) two distinct aquifers, their hydraulic conductivities and their hydrodynamic responses to climate forcing are identified; (ii) a comprehensive water balance for the whole catchment is obtained; (iii) the subglacial recharge to the aquifers is shown to be 4 times higher than in the proglacial area; and (v) the importance of the impact of the glacial melt recharge on the groundwater system is demonstrated. Thus, we highlight the major role that the groundwater component has in the hydrodynamic functioning of glacierized catchments.

Multitemporal characterization of a proglacial system: a multidisciplinary approach

Abstract. The recession of Alpine glaciers causes an increase in the extent of proglacial areas and leads to changes in the water discharge and sediment balance (morphodynamics and sediment transport). Although the processes occurring in proglacial areas are relevant not only from a scientific point of view but also for the purpose of climate change adaptation, there is a lack of work on the continuous monitoring and multitemporal characterization of these areas. This study offers a multidisciplinary approach that merges the contributions of different scientific disciplines, such as hydrology, geophysics, geomatics, and water engineering, to characterize the Rutor Glacier and its proglacial area. Since 2020, we have surveyed the glacier and its proglacial area using both uncrewed and crewed aerial surveys (https://doi.org/10.5281/zenodo.8089499, Corte et al., 2023c; https://doi.org/10.5281/zenodo.10100968, Corte et al., 2023f; https://doi.org/10.5281/zenodo.10074530, Corte et al., 2023g; https://doi.org/10.5281/zenodo.10101236, Corte et al., 2023h; https://doi.org/10.5281/zenodo.7713146, Corte et al., 2023b). We have determined the bathymetry of the most downstream proglacial lake and the thickness of the sediments deposited on its bottom (https://doi.org/10.5281/zenodo.7682072, Corte et al., 2023a). The water depth at four different locations within the hydrographic network of the proglacial area (https://doi.org/10.5281/zenodo.7697100, Corte et al., 2023d) and the bedload at the glacier snout (https://doi.org/10.5281/zenodo.7708800, Corte et al., 2023e) have also been continuously monitored. The synergy of our approach enables the characterization, monitoring, and understanding of a set of complex and interconnected processes occurring in a proglacial area.

Accelerating Glacier Area Loss Across the Andes Since the Little Ice Age

AbstractAndean glaciers are losing mass rapidly but a centennial‐scale context to those rates is lacking. Here we show the extent of >5,500 glaciers during the Little Ice Age chronozone (LIA; c. 1,400 to c. 1,850) and compute an overall area change of −25% from then to year 2000 at an average rate of −36.5 km2 yr−1 or −0.11% yr−1. Glaciers in the Tropical Andes (Peru, Bolivia) have depleted the most; median −56% of LIA area, and the fastest; median −0.16% yr−1. Up to 10 × acceleration in glacier area loss has occurred in Tropical mountain sub‐regions comparing LIA to 2,000 rates to post‐2000 rates. Regional climate controls inter‐regional variability, whereas local factors affect intra‐region glacier response time. Analyzing glacier area change by river basins and by protected areas leads us to suggest that conservation and environmental management strategies should be re‐visited as proglacial areas expand.

The impact of glacier retreat on Andean high wetlands: Assessing the geochemical transfer and sediment provenance in the proglacial area of Huayna-Potosí (Bolivia)

The accelerated decline of tropical Andean glaciers is affecting the water cycle of mountains in the region. In the Cordillera Real (Bolivia) on the Huayna-Potosí peak (6088 m a.s.l.), the rapid retreat of the Oeste Glacier is exposing new rock outcrops and glacial materials. Changes in the hydrological regime by glacier retreat are likely to modify the supply of sediments and subsequently affect the geochemical transfer of particle-bound elements within proglacial and downstream ecosystems. To address this issue, a sampling campaign in the proglacial area of Huayna-Potosí aimed to identify and collect sediments from the primary active landforms. To characterise elemental transfer and identify sediment provenance using fingerprinting techniques, streambed sediments were collected along the main stream from the glacier tongue at 5002 m a.s.l. to a wetland shallow lake at 4700 m a.s.l.. Geomorphological features control the production of fine particles and the distribution of particle-bound geochemical elements. Complex processes of solubilization, mobilisation, and retention govern the elemental transfer. Significant increases in 238U, as well as some major elements, were observed in the surface lake sediments, while Pb, Cr, and Cu decreased, and 137Cs was not detected. The spatial distribution patterns of geochemical elements and the assessment of the sediment provenance indicate that deglaciation fosters geochemical processes and the elemental transfer dynamics in the wetlands. Moraines and colluvium were found to be the primary contributors to sediment yield (93 %) while the rich organic soils in swamps contributed the least (7 %). The enhanced supply of sediments to the shallow lake in parallel to glacier decline may affect water quality and bring about further changes in the wetlands. Identifying primary sediment sources is essential for water management. Our findings enhance understanding of the compositional characteristics of mobilised sediment and their variations, aiding in more effective resource management to preserve these high-altitude wetlands.

Controls on Sediment Transport From a Glacierized Catchment in the Swiss Alps Established Through Inverse Modeling of Geomorphic Processes

AbstractAs climate warms, hydrology and geomorphology in glacierized catchments are evolving, changing sediment export from these catchments, thus impacting downstream ecosystems and communities. Currently, much uncertainty persists regarding interactions among geomorphic processes that evacuate sediment from glacierized catchments. Here, we present a catchment‐scale numerical model of subglacial and proglacial sediment transport with debris meltout processes. We apply the model in a Monte Carlo framework to suspended sediment data collected over 2014–2020 from the Fieschertal catchment in the Swiss Alps, assessing possible combinations of geomorphic processes responsible for suspended sediment discharge. The ensemble of model outputs quantifies the interaction of different geomorphic processes, including the trade‐off between bedrock erosion and sediment previously stored in the catchment. Model runs suggest that, at some periods, up to 20% of the sediment leaving the glacier is deposited in the proglacial area relative to the catchment's total sediment discharge. This shows that the model captures the proglacial area change from sediment source to sediment sink in different hydrological and glaciological conditions. Furthermore, the findings highlight the impact of glacier retreat on the catchment's sediment dynamics, which both reduces the proglacial area's slope and introduces additional sediment to the proglacial area through debris meltout. The model outputs and the parameters quantify the interaction among geomorphic processes, yielding key insights into the drivers of sediment transport in alpine regions. The implications of these interactions are discussed in the context of interpreting processes responsible for controlling erosion rates from glacierized regions and modeling sediment transport in glacierized catchments.

Assessment of plot-scale sediment transport on young moraines in the Swiss Alps using a fluorescent sand tracer

Abstract. Glacial retreat uncovers large bodies of unconsolidated sediment that are prone to erosion. However, our knowledge of overland flow (OF) generation and sediment transport on moraines that have recently become ice-free is still limited. To investigate how the surface characteristics of young moraines affect OF and sediment transport, we installed five bounded runoff plots on two moraines of different ages in a proglacial area of the Swiss Alps. On each plot we conducted three sprinkling experiments to determine OF characteristics (i.e., total OF and peak OF flow rate) and measured sediment transport (turbidity, sediment concentrations, and total sediment yield). To determine and visualize where sediment transport takes place, we used a fluorescent sand tracer with an afterglow as well as ultraviolet (UV) and light-emitting diode (LED) lamps and a high-resolution camera. The results highlight the ability of this field setup to detect sand movement, even for individual fluorescent sand particles (300–500 µm grain size), and to distinguish between the two main mechanisms of sediment transport: OF-driven erosion and splash erosion. The higher rock cover on the younger moraine resulted in longer sediment transport distances and a higher sediment yield. In contrast, the higher vegetation cover on the older moraine promoted infiltration and reduced the length of the sediment transport pathways. Thus, this study demonstrates the potential of the use of fluorescent sand with an afterglow to determine sediment transport pathways as well as the fact that these observations can help to improve our understanding of OF and sediment transport processes on complex natural hillslopes.

Supervised Geomorphic Mapping of Himalayan Rivers Based on Sentinel-2 Data

The Himalayan region is a hotspot in terms of expected future hydrological and geomorphological variations induced by climate change on proglacial areas and the related implications for human societies established along the downstream rivers. Due to the remoteness of the proglacial zones in the Himalayas and the associated logistical problems in carrying out traditional field and UAV-based morphological monitoring activities, remote sensing here plays a crucial role to monitor past and current fluvial dynamics, which could be used to anticipate future changes; however, there has been, so far, limited research on morphological changes in Himalayan proglacial rivers. To address this gap, a morphological classification model was designed to classify recent changes in Himalayan proglacial rivers using the Google Earth Engine platform. The model is the first of its kind developed for the Himalayan region and uses multispectral S-2 satellite data to delineate submerged water channels, vegetated surfaces, and emerged, unvegetated sediment bars, and then to track their variations over time. The study focused on three training sites: Langtang-Khola (Nepal), Saltoro (Pakistan), and Nubra (Jammu and Kashmir) rivers, and one testing site, the Ganga-Bhagirathi River (India). A total of 900 polygons were used as training samples for the random forest classifier, which were further divided into 70% calibration and 30% validation datasets for the training sites, and a separate validation dataset was acquired from the testing site to assess the model performance. The model achieved high accuracy, with an average overall accuracy of 96% and a kappa index of 0.94, indicating the reliability of the S2 data for modeling proglacial geomorphic features in the Himalayan region. Therefore, this study provides a reliable tool to detect past and current morphological changes occurring in the Himalayan proglacial rivers, which will be of great value for both research and river management purposes.

Modeling of surface energy balance for Icelandic glaciers using remote-sensing albedo

Abstract. During the melt season, absorbed solar energy, modulated at the surface by albedo, is one of the main governing factors controlling surface melt variability for glaciers in Iceland. An energy balance model was applied with the possibility of utilizing spatiotemporal Moderate Resolution Imaging Spectroradiometer (MODIS) satellite-derived daily surface albedo driven by high-resolution climate forcing data to reconstruct the surface energy balance (SEB) for all Icelandic glaciers for the period 2000–2021. The SEB was reconstructed from April through September for 2000–2021 at a daily time step with a 500 m spatial resolution. Validation was performed using observations from various glaciers spanning distinct locations and elevations with good visual and statistical agreement. The results show that spatiotemporal patterns for the melt season have high annual and interannual variability for Icelandic glaciers. The variability was influenced by high climate variability, deposition of light-absorbing particles (LAPs) from volcanic eruptions and dust hotspots in pro-glacial areas close to the glaciers. Impacts of LAPs can lead to significant melt enhancement due to lowering of albedo and increased short-wave radiative energy forced at the surface. Large impacts on the SEB were observed for years with high LAP deposits, such as the volcanic eruption years of 2004, 2010 and 2011 and the sand- and dust-rich year of 2019. The impacts of volcanic eruptions and other LAP events were estimated using historical mean albedo under the same climatology forcing to provide estimations of melt energy enhancements. The impact of LAPs was often significant even though the glaciers were far away from the eruption location. On average, the melt enhancements due to LAPs were ∼27 % in 2010, ∼16 % in 2011 and ∼14 % in 2019 for Vatnajökull, Hofsjökull and Langjökull.

Lateral moraine failure in the valley of the Djankuat Catchment (Central Caucasus) and subsequent morphodynamics

Moraine failures are severe events in proglacial areas that lead to major surface changes. On July 1, 2015, a lateral moraine collapsed within the 9.1 km2 glacierized Djankuat Catchment in the Northern Caucasus. This study quantifies the surface changes over a 7-year period after the collapse and proposes the main event triggers using hydrometeorological measurements and remote sensing data. Digital Elevation Models (DEMs), derived from Unmanned Aerial Vehicle (UAV) surveys and satellite images for the 2010–2022 period, were used to estimate erosion/deposition associated with the event and further surface changes in several zones: the breach, the fan, the headward erosion area and gully erosion on lateral moraine slopes. The failure occurred due to an extreme precipitation event (7-day cumulative sum: 227 mm, 24 h prior: 60 mm), simultaneous snowmelt within the whole upstream catchment, and continuous positive temperatures. During the 2015 event, approximately 200,000 m3 of material was instantly eroded, corresponding to a surface lowering of 1.91 m a−1 (over 52,227 m2). Only 25 % of the eroded material redeposited on the fan. Later in 2017–2019, the breach underwent infilling with a 10 m deposition layer in the bottom (net change ca. 14,600 m3), while the stream in the headward erosion area continued to incise (ca. −12,900 m3). From 2019 to 2022, the dynamics decreased, and a positive sediment budget was assessed for all zones except the intensive local erosion on the fan in 2021–2022 (ca. −15,200 m3). Gully changes suggested a constant erosion of 90–100 mm a−1. Based on daily rainfall measurements for the 2008–2022 period, we proposed a 3-day and prior 15-day cumulative precipitation threshold that outlines major erosional events in the Djankuat River valley. A cumulative 3-day precipitation of 100 mm could be a sufficient but not necessary condition since some significant sediment export was recorded below that threshold.

Polycyclic aromatic hydrocarbon dynamics in soils along proglacial chronosequences in the Alps

Polycyclic aromatic hydrocarbons (PAHs) were studied in the soils of three proglacial areas in France (Noir and Chardon Glaciers) and Italy (Miage Glacier). PAH contents, PAH stocks and PAH contents normalized to the total organic carbon contents (PAHs/TOC ratio) were investigated along proglacial soil chronosequences to infer their evolutions with soil age (from 3 to 4200 years), where the PAH contamination was only related to long-range atmospheric transport. Evolutions of PAH and TOC contents, PAHs/TOC ratio and PAH stock were fitted with exponential and logarithmic relations. For the three proglacial areas, PAH contents increased rapidly during the first 150 years of soil development, ranged from 4 to 152 ng·g−1, and showed a strong relationship with total organic carbon (TOC) contents (r = 0.83, p < 0.05). The joint increase of PAH and TOC contents suggested that PAH accumulation in soils were not only driven by PAH inputs but also by the capacity of soils to store these contaminants. PAH contents in the oldest soils (from 1200 BCE and 2200 BCE) were similar than for soils from 1850 CE. The period 1850–2019 CE corresponded to a decrease in the PAHs/TOC ratio suggesting both a faster accumulation of TOC than PAHs and a dilution effect of PAHs already present in soils. For the oldest soils, the PAHs/TOC ratio appeared similar to those for soils from 1850 CE, with values ranging from 0.48 to 2.06 ng·mg−1, suggesting an equilibrium between both parameters for soils older than 170 years. Finally, PAH stocks ranged from 0.41 mg·m−2 to 6.80 mg·m−2 in the youngest and oldest soils, respectively. These results do not allow us to identify the same period of greatest emission as other studies (estimated ~1960), but they revealed changes in the capacity of soils to store these pollutants.

Primary succession and its driving variables – a sphere-spanning approach applied in proglacial areas in the upper Martell Valley (Eastern Italian Alps)

Abstract. Climate change and the associated glacier retreat lead to considerable enlargement and alterations of the proglacial systems. The colonisation of plants in this ecosystem was found to be highly dependent on terrain age, initial site conditions and geomorphic disturbances. Although the explanatory variables are generally well understood, there is little knowledge on their collinearities and resulting influence on proglacial primary succession. To develop a sphere-spanning understanding of vegetation development, a more interdisciplinary approach was adopted. In the proglacial areas of Fürkeleferner, Zufallferner and Langenferner (Martell Valley, Eastern Italian Alps), in total 65 plots of 5×2 m were installed to perform the vegetation analysis on vegetation cover, species number and species composition. For each of those, 39 potential explanatory variables were collected, selected through an extensive literature review. To analyse and further avoid multicollinearities, 33 of the explanatory variables were clustered via principal component analysis (PCA) to five components. Subsequently, generalised additive models (GAMs) were used to analyse the potential explanatory factors of primary succession. The results showed that primary succession patterns were highly related to the first component (elevation and time), the second component (solar radiation), the third component (soil chemistry), the fifth component (soil physics) and landforms. In summary, the analysis of all explanatory variables together provides an overview of the most important influencing variables and their interactions; thus it provides a basis for the debate on future vegetation development in a changing climate.

Estimating the Evolution of a Post-Little Ice Age Deglaciated Alpine Valley through the DEM of Difference (DoD)

Since the end of the Little Ice Age (LIA, ~1830), the accelerated glaciers’ shrinkage along mid-latitude high mountain areas promoted a quick readjustment of geomorphological processes with the onset of the paraglacial dynamic, making proglacial areas among the most sensitive Earth’s landscapes to ongoing climate change. A potentially useful remote-sensing method for investigating such dynamic areas is the DEM (Digital Elevation Model) of Difference (DoD) technique, which quantifies volumetric changes in a territory between successive topographic surveys. After a detailed geomorphological analysis and comparison with historical maps of the Martello Valley (central Italian Alps), we applied the DoD for reconstructing post-LIA deglaciation dynamics and reported on the surface effects of freshly-onset paraglacial processes. The head of the valley is still glacierized, with three main ice bodies resulting from the huge reduction of a single glacier present at the apogee of the LIA. Aftermath: the glaciers lose 60% of their initial surface area, largely modifying local landforms and expanding the surface of the proglacial areas. The DoD analysis of the 2006–2015 timeframe (based on registered DEM derived from LiDAR—Light Detection and Ranging—data) highlights deep surface elevation changes ranging from +38 ± 4.01 m along the foot of rock walls, where gravitative processes increased their intensity, to −47 ± 4.01 m where the melting of buried ice caused collapses of the proglacial surface. This approach permits estimating the volume of sediments mobilized and reworked by paraglacial processes. Here, in less than 10 years, −23,675 ± 1165 m3 of sediment were removed along the proglacial area and transported down valley, highlighting the dynamicity of proglacial areas.

Cooling glaciers in a warming climate since the Little Ice Age at Qaanaaq, northwest Kalaallit Nunaat (Greenland)

AbstractThe centennial response of land‐terminating glaciers in Greenland to climate change is largely unknown. Yet, such information is important to understand ongoing changes and for projecting the future evolution of Arctic subpolar glaciers, meltwater runoff, and sediment fluxes. This paper analyses the topography, geomorphology, and sedimentology of prominent moraine ridges and the proglacial areas of ice cap outlet glaciers on the Qaanaaq peninsula (Piulip Nunaa). We determine geometric changes of glaciers since the neoglacial maximum; the Little Ice Age (LIA), and we compare glacier behaviour during the LIA with that of the present day. There has been very little change in the rate of volume loss of each outlet glacier since the LIA compared with the rate between 2000 and 2019. However, the percentage of each glacier that is likely composed of cold‐based ice has increased since the LIA, typically by 20%. The LIA moraines comprise subrounded, striated, and faceted clasts that evidence subglacial transport, and outwash plains, flutes, kames, and eskers that evidence subglacial motion and meltwater within temperate ice. Contrastingly, contemporary ice margins and their convex ice surfaces comprise pronounced primary foliation, ephemeral supraglacial drainage, sediment drapes from thrust plane fractures, and an absence of open crevasses and moulins. These calculations and observations together lead us to interpret that these outlet glaciers have transitioned towards an increasingly cold‐based thermal regime despite a warming regional climate. Thermal regime transitions control glacier dynamics and therefore should be incorporated into glacier evolution models, especially where polythermal glaciers prevail and where climate is changing rapidly.

Quantitative Long-Term Monitoring (1890–2020) of Morphodynamic and Land-Cover Changes of a LIA Lateral Moraine Section

Aerial photographs of the European Alps usually only reach back to the middle of the 20th century, which limits the time span of corresponding studies that quantitatively analyse long-term surface changes of proglacial areas using georeferenced orthophotos. To the end of the Little Ice Age, this leads to a gap of about 100 years. Using digital monoplotting and several historical terrestrial photographs, we show the quantification of surface changes of a Little Ice Age lateral moraine section until the late second half of the 19th century, reaching a total study period of 130 years (1890–2020). The (initial) gully system expands (almost) continuously over the entire study period from 1890 to 2020. Until 1953, the vegetation-covered areas also expanded (mainly scree communities, alpine grasslands and dwarf shrub communities), before decreasing again, especially between 1990 and 2003, due to large-scale erosion within the gully system. Furthermore, our results show that the land-cover development was impacted by temperature and precipitation changes. With the 130-year study period, we contribute to a substantial improvement in the understanding of the processes in the proglacial by analysing the early phase and thus the immediate response of the lateral moraine to the ice exposure.

The vertical atmospheric structure of the partially glacierised Mittivakkat valley, southeast Greenland

Air temperature inversions, a situation in which atmospheric temperature increases with height, are key components of the Arctic planetary boundary layer. The present study investigates the spatial and temporal variations of temperature inversions over different surface types (rock, gravel, snow, ice) along the Mittivakkat valley (southeast Greenland). For this purpose, 113 vertical profiles with high spatio-temporal resolution of air temperature and relative humidity were collected with unoccupied aerial vehicles (UAVs) during a 13-day field campaign in summer 2019. Air temperature inversions were present in 83% of the profiles, of which 24% were surface-based inversions and 76% were elevated inversions. The proglacial area covered with bare rock and gravel induces surface heating and convection during the day and, through interaction with local circulation patterns, leads to the frequent formation of elevated inversions. In contrast, the glacier surface itself acts as a persistent cooling surface and leads to the formation of surface-based inversions. A low-level fog layer that forms under the inversion layer may be causing non-linear vertical ablation gradients on Mittivakkat Gletsjer. Furthermore, we demonstrate that atmospheric measurements using UAVs can better capture small-scale processes than other products like radiosonde or modeled reanalysis data.

Hydrological connections in a glaciated Andean catchment under permafrost conditions (33°S)

Fresh water supply is critical along the Andes, where drought conditions over the past decade are projected to persist. At high Andean headwater catchments, frozen ground conditions are assumed to modulate groundwater flow paths and their hydrological signals at different timescales. However, knowledge of hydrological connections in subtropical Andean catchments is still very sparse. This study assessed hydrological connections and their impacts on groundwater contribution to baseflow in a headwater proglacial aquifer located in central Chile at 33° S and 3600 m a.s.l. We collected and analyzed snow, glacial stream, and groundwater spring water samples between 2019 and 2021. We combined of water isotope and metagenomic proxies with the hydraulic parameterization of the catchment to deliver mean transit time distributions through the proglacial aquifer. The new hydrological insights for the region include the finding that groundwater spring signals delivered sub-decadal transit times, implying likely origins from glacial or interstitial ice. Additionally, the stable isotope signature showed that groundwater consistently differs from snow and surface runoff. The 16S rRNA metabarcoding analyses demonstrated the presence of psychrophilic microorganisms in groundwater springs, supporting the idea of a late warm-season activation of interstitial ice due to thawing events associated with a differential relative-abundance of specific cryophilic bacteria. Finally, our results suggest hydrological connections and dampening timeframes between glaciers, proglacial areas, and groundwater springs, most likely from thawing sources.

Catchment characteristics and seasonality control the composition of microbial assemblages exported from three outlet glaciers of the Greenland Ice Sheet.

Glacial meltwater drains into proglacial rivers where it interacts with the surrounding landscape, collecting microbial cells as it travels downstream. Characterizing the composition of the resulting microbial assemblages in transport can inform us about intra-annual changes in meltwater flowpaths beneath the glacier as well as hydrological connectivity with proglacial areas. Here, we investigated how the structure of suspended microbial assemblages evolves over the course of a melt season for three proglacial catchments of the Greenland Ice Sheet (GrIS), reasoning that differences in glacier size and the proportion of glacierized versus non-glacierized catchment areas will influence both the identity and relative abundance of microbial taxa in transport. Streamwater samples were taken at the same time each day over a period of 3 weeks (summer 2018) to identify temporal patterns in microbial assemblages for three outlet glaciers of the GrIS, which differed in glacier size (smallest to largest; Russell, Leverett, and Isunnguata Sermia [IS]) and their glacierized: proglacial catchment area ratio (Leverett, 76; Isunnguata Sermia, 25; Russell, 2). DNA was extracted from samples, and 16S rRNA gene amplicons sequenced to characterize the structure of assemblages. We found that microbial diversity was significantly greater in Isunnguata Sermia and Russell Glacier rivers compared to Leverett Glacier, the latter of which having the smallest relative proglacial catchment area. Furthermore, the microbial diversity of the former two catchments continued to increase over monitored period, presumably due to increasing hydrologic connectivity with proglacial habitats. Meanwhile, diversity decreased over the monitored period in Leverett, which may have resulted from the evolution of an efficient subglacial drainage system. Linear discriminant analysis further revealed that bacteria characteristic to soils were disproportionately represented in the Isunnguata Sermia river, while putative methylotrophs were disproportionately abundant in Russell Glacier. Meanwhile, taxa typical for glacierized habitats (i.e., Rhodoferax and Polaromonas) dominated in the Leverett Glacier river. Our findings suggest that the proportion of deglaciated catchment area is more influential to suspended microbial assemblage structure than absolute glacier size, and improve our understanding of hydrological flowpaths, particulate entrainment, and transport.

The relative imprint of forming factors on soil characteristics in a recently deglaciated area: concerns about chronosequences approach

ABSTRACT In high-altitude environments, most of the debris-free glacier forelands are increasingly widening and the bare surfaces left by retreating glaciers offer the opportunity to investigate the evolution of soil through time. The main aim of this study is to discuss the applicability of a chronosequence approach in a deglaciated area, by considering the role of different soil forming factors. A study case has been selected (Alpe Veglia, Lepontine Alps), where field and laboratory characterizations were performed along a transect of 12 soil profiles from the proglacial area of the Aurona Glacier to the Alpe Veglia hollow, crossing different age glacial deposits. The results of soil physical and chemical analyses underline a time-trend of soil properties. On the other hand, even if the soil properties variability along the transect can be mainly explained according to the soil chronosequence approach, the data seem to highlight how the various soil forming factors (e.g. parent material, relief, vegetation) influence soil features, partly masking the effect of the time factor inducing a divergence from a traditional chronosequence. In particular, the morphology of the surrounding reliefs, depending on bedrock lithologies and structures, and the geomorphic dynamics seem to affect soil formation and evolution.