Genetic variation in the small bivalve Nuculana inaequisculpta along a retreating glacier fjord, King George Island, Antarctica

Climate change is strongly influencing regions of Antarctica but the consequences on microevolutionary processes have been little studied. Patterns of population genetic diversity were analysed in the Antarctic bivalve Nuculana inaequisculpta (Protobranchia: Nuculanidae) from a fjord with 70 years of documented climate-forced glacier retreat. Thirty-nine individuals from five sites at different distances from the glacier terminus were collected, and the COI gene was sequenced from each individual. No statistically significant genetic differentiation was found between sites nor a significant correlation between the proximity of glaciers and genetic diversity, suggesting a high dispersal capability and therefore, a planktonic larval stage for this species. Nevertheless, we encourage increasing the sample size and number of loci in future studies to confirm our findings.

Mid-Term Monitoring of Glacier’s Variations with UAVs: The Example of the Belvedere Glacier

Recently, Unmanned Aerial Vehicles (UAV) have opened up unparalleled opportunities for alpine glacier monitoring, as they allow for reconstructing extensive and high-resolution 3D models. In order to evaluate annual ice flow velocities and volume variations, six yearly measurements were carried out between 2015 and 2020 on the debris-covered Belvedere Glacier (Anzasca Valley, Italian Alps) with low-cost fixed-wing UAVs and quadcopters. Every year, ground control points and check points were measured with GNSS. Images acquired from UAV were processed with Structure-from-Motion and Multi-View Stereo algorithms to build photogrammetric models, orthophotos and digital surface models, with decimetric accuracy. Annual glacier velocities were derived by combining manually-tracked features on orthophotos with GNSS measurements. Velocities ranging between 17 m y−1 and 22 my−1 were found in the central part of the glacier, whereas values between 2 m y−1 and 7 my−1 were found in the accumulation area and at the glacier terminus. Between 2 × 106 m3 and 3.5 × 106m3 of ice volume were lost every year. A pair of intra-year measurements (October 2017–July 2018) highlighted that winter and spring volume reduction was ∼1/4 of the average annual ice loss. The Belvedere monitoring activity proved that decimetric-accurate glacier models can be derived with low-cost UAVs and photogrammetry, limiting in-situ operations. Moreover, UAVs require minimal data acquisition costs and allow for great surveying flexibility, compared to traditional techniques. Information about annual flow velocities and ice volume variations of the Belvedere Glacier may have great value for further understanding glacier dynamics, compute mass balances, or it might be used as input for glacier flow modelling.

Automated detection and analysis of surface calving waves with a terrestrial radar interferometer at the front of Eqip Sermia, Greenland

Glacier calving is a key dynamical process of the Greenland Ice Sheet and a major driver of its increasing mass loss. Calving waves, generated by the sudden detachment of ice from the glacier terminus, can reach tens of meters in height and provide very valuable insights into quantifying calving activity. In this study, we present a new method for the detection of source location, timing, and magnitude of calving waves using a terrestrial radar interferometer. This method was applied to 11 500 1 min interval acquisitions from Eqip Sermia, West Greenland, in July 2018. Over 7 d, more than 2000 calving waves were detected, including waves generated by submarine calving, which are difficult to observe with other methods. Quantitative assessment with a wave power index (WPI) yields a higher wave activity (+49 %) and higher temporally cumulated WPI (+34 %) in deep water than under shallow conditions. Subglacial meltwater plumes, occurring 2.3 times more often in the deep sector, increase WPI and the number of waves by a factor of 1.8 and 1.3, respectively, in the deep and shallow sector. We therefore explain the higher calving activity in the deep sector by a combination of more frequent meltwater plumes and more efficient calving enhancement linked with better connections to warm deep ocean water.

Ice dynamics of the Haupapa/Tasman Glacier measured at high spatial and temporal resolution, Aoraki/Mount Cook, New Zealand

<p>Glaciers are among the clearest of signals for anthropogenic climate change and their retreat is considered symptomatic of the observed warming since the start of the 20th century from anthropogenic sources (Mann et al., 2004). New Zealand has 3,100 mountain glaciers, with those in the Southern Alps experiencing losses of 34% since 1977 and a decline in volume of 51 km3 in 1994 to 41 km3 in 2010 (NIWA, 2011). The direct impact of increasing atmospheric temperatures on glaciers is well understood (Chinn, 2012) through its effects on the melt and accumulation rates (Kirkbride, 2010; Purdie, 2011; Chinn, 1997; Oerlemans, 2001). However lake calving glaciers such as the Tasman Glacier exhibit different behaviour and are suggested to be at least partially decoupled from climate forcing (Benn et al., 2007). Here, I present a temporally and spatially complete study of Haupapa/Tasman Glacier, Aoraki/Mt. Cook over three years to investigate the ice dynamics at the terminus. I used oblique photogrammetry at high resolution for data acquisition and adapted computer vision algorithms for correcting this oblique view to a real-world geometry. This technique has been rarely used (Murray et al., 2015; Messerli and Grinsted, 2015; Ahn and Box, 2010; Harrison et al., 1986 and Flotron, 1973) but owing to its cost-effectiveness and high data yields, it is becoming an increasingly powerful methodology favoured by glaciologists. During the 3 year study period, Tasman Glacier terminus retreat rate Ur was 116 ± 19 m a⁻¹ (2013-2014), 83 ± 18 m a⁻¹ (2014-2015) and 204 ± 20 (2015-2016). A strong seasonal pattern was evident in the calving events. Three major calving events occurred over the study, one occurring in the summer of 2013 and two in the summer of 2016. The latter two events are responsible for the elevated Ur in 2015-2016. These events were characterised as distinct large-magnitude calving (usually as a large tabular iceberg) which continued to drift and break up in the lake for weeks to months. Three large calving events accounted for 47% of the total surface area loss for the 38 month study period with the remaining surface area loss from 2nd order calving including notching at the waterline and the spalling of lamallae of ice from surface fractures, and ice-cliff melt. During the spring/summer months of 2014 and 2015 there was no large buoyancy driven calving event such as those seen in 2013 and 2016, but there were many smaller-magnitude calving events. Smaller-magnitude events were less frequent in winter months as compared to summer months. Ice flow in winter has been shown to be less than in summer (Horgan et al, 2015). While seasonal temperatures and changes to the basal water pressure are linked to these observations, it is also likely that the relatively faster ice flow in summer/autumn could be influencing the rate of 1st and 2nd order calving mechanisms. Overall, the calving rates were calculated as 171 ± 18 m a⁻¹ (2013-2014), 136 ± 17 m a⁻¹ (2014-2015) and accelerated to 256 ± 20 m a⁻¹ in the last year (2015-2016). My results show that almost half of the ice loss at the terminus comes from large, infrequent calving events and that retreat rates for 2015-2016 were high compared to the historic record but the area loss is lower than it has been because of the relatively narrow terminus.</p>

Ice dynamics of the Haupapa/Tasman Glacier measured at high spatial and temporal resolution, Aoraki/Mount Cook, New Zealand

<p>Glaciers are among the clearest of signals for anthropogenic climate change and their retreat is considered symptomatic of the observed warming since the start of the 20th century from anthropogenic sources (Mann et al., 2004). New Zealand has 3,100 mountain glaciers, with those in the Southern Alps experiencing losses of 34% since 1977 and a decline in volume of 51 km3 in 1994 to 41 km3 in 2010 (NIWA, 2011). The direct impact of increasing atmospheric temperatures on glaciers is well understood (Chinn, 2012) through its effects on the melt and accumulation rates (Kirkbride, 2010; Purdie, 2011; Chinn, 1997; Oerlemans, 2001). However lake calving glaciers such as the Tasman Glacier exhibit different behaviour and are suggested to be at least partially decoupled from climate forcing (Benn et al., 2007). Here, I present a temporally and spatially complete study of Haupapa/Tasman Glacier, Aoraki/Mt. Cook over three years to investigate the ice dynamics at the terminus. I used oblique photogrammetry at high resolution for data acquisition and adapted computer vision algorithms for correcting this oblique view to a real-world geometry. This technique has been rarely used (Murray et al., 2015; Messerli and Grinsted, 2015; Ahn and Box, 2010; Harrison et al., 1986 and Flotron, 1973) but owing to its cost-effectiveness and high data yields, it is becoming an increasingly powerful methodology favoured by glaciologists. During the 3 year study period, Tasman Glacier terminus retreat rate Ur was 116 ± 19 m a⁻¹ (2013-2014), 83 ± 18 m a⁻¹ (2014-2015) and 204 ± 20 (2015-2016). A strong seasonal pattern was evident in the calving events. Three major calving events occurred over the study, one occurring in the summer of 2013 and two in the summer of 2016. The latter two events are responsible for the elevated Ur in 2015-2016. These events were characterised as distinct large-magnitude calving (usually as a large tabular iceberg) which continued to drift and break up in the lake for weeks to months. Three large calving events accounted for 47% of the total surface area loss for the 38 month study period with the remaining surface area loss from 2nd order calving including notching at the waterline and the spalling of lamallae of ice from surface fractures, and ice-cliff melt. During the spring/summer months of 2014 and 2015 there was no large buoyancy driven calving event such as those seen in 2013 and 2016, but there were many smaller-magnitude calving events. Smaller-magnitude events were less frequent in winter months as compared to summer months. Ice flow in winter has been shown to be less than in summer (Horgan et al, 2015). While seasonal temperatures and changes to the basal water pressure are linked to these observations, it is also likely that the relatively faster ice flow in summer/autumn could be influencing the rate of 1st and 2nd order calving mechanisms. Overall, the calving rates were calculated as 171 ± 18 m a⁻¹ (2013-2014), 136 ± 17 m a⁻¹ (2014-2015) and accelerated to 256 ± 20 m a⁻¹ in the last year (2015-2016). My results show that almost half of the ice loss at the terminus comes from large, infrequent calving events and that retreat rates for 2015-2016 were high compared to the historic record but the area loss is lower than it has been because of the relatively narrow terminus.</p>

Inferring Past Climate from Moraine Evidence Using Glacier Modelling

<p>Glacier length fluctuations reflect changes in climate, most notably temperature and precipitation. By this reasoning, moraines, which represent former glacier extent, can be used to estimate past climate. However, estimating palaeoclimate from moraines is not a straight-forward process and involves several assumptions. For example, recent studies have suggested that interannual stochastic variability in temperature in a steady-state climate can cause a glacier to experience kilometre-scale fluctuations. Such studies cast doubt on the usefulness of moraines as climate proxy indicators. Detailed glacial geomorphological maps and moraine chronologies have improved our understanding of the spatial and temporal extent of past glacial events in New Zealand. Palaeoclimate estimates associated with these moraines have thus-far come from simple methods, such as the accumulation area ratio, with unquantifiable uncertainties. I used a numerical modelling approach to approximate the present-day glacier mass balance pattern, which includes the effects of snow avalanching on glacier mass balance. I then used the models to reconstruct palaeoclimate for Lateglacial and Holocene glacial events in New Zealand, and to better understand moraine-glacier-climate relationships. The climate reconstructions come from simulating past glacier expansions to specific terminal moraines, but I also simulated glacier fluctuations in response to a previously derived temperature reconstruction, and to interannual stochastic variability in temperature. The purpose behind each simulation was to identify the drivers of significant glacier fluctuations. The modelling results support the hypothesis that New Zealand moraine records reflect past climate, especially changes in temperature. Lateglacial climate was reconstructed to be 2-3 C lower than the present day. This temperature range agrees well with previous estimates from moraines and other climate proxy records in New Zealand. Modelled temperature estimates for the Holocene moraines are slightly colder than those derived from simpler methods, due to a non-linear relationship found between snowline lowering and glacier length. This relationship results from the specific valley shape and glacier geometry, and is likely to occur in other, similarly-shaped glacier valleys. The simulations forced by interannual stochastic variability in temperature do not show significant (>300 m) fluctuations in the glacier terminus. Such fluctuations can not explain the Holocene moraine sequence that I examined, which extends >2 km beyond the present-day glacier terminus. Stochastic temperature change could, however, in part, cause fluctuations in glacier extent during an overall glacier recession. Modelling shows that it is also unlikely that glaciers advanced to Holocene and Lateglacial moraine positions as a result of precipitation changes alone. For these reasons, temperature changes are a necessary part of explaining past glacier extents, especially during the Lateglacial, and the moraines examined here likely reflect changes in mean climate in New Zealand. The glacier modelling studies indicate that simpler methods, such as the accumulation area ratio, can be used to appropriately reconstruct past climate from glacial evidence, as long as the glacier catchment has a straight forward geometry, shallow bed slope and no tributary glaciers. Non-linear relationships between climate change and glacier length develop when valley shape is more complex, and glaciers within these systems are probably better simulated using a modelling approach. Using a numerical modelling approach, it is also possible to gain a greater understanding of glacier response time, length sensitivities, and estimates of ice extent in valleys within the model domain where geomorphic evidence is not available. In this manner, numerical models can be used as a tool for understanding past climate and glacier sensitivity, thus improving the confidence in the palaeoclimate interpretations.</p>

Inferring Past Climate from Moraine Evidence Using Glacier Modelling

<p>Glacier length fluctuations reflect changes in climate, most notably temperature and precipitation. By this reasoning, moraines, which represent former glacier extent, can be used to estimate past climate. However, estimating palaeoclimate from moraines is not a straight-forward process and involves several assumptions. For example, recent studies have suggested that interannual stochastic variability in temperature in a steady-state climate can cause a glacier to experience kilometre-scale fluctuations. Such studies cast doubt on the usefulness of moraines as climate proxy indicators. Detailed glacial geomorphological maps and moraine chronologies have improved our understanding of the spatial and temporal extent of past glacial events in New Zealand. Palaeoclimate estimates associated with these moraines have thus-far come from simple methods, such as the accumulation area ratio, with unquantifiable uncertainties. I used a numerical modelling approach to approximate the present-day glacier mass balance pattern, which includes the effects of snow avalanching on glacier mass balance. I then used the models to reconstruct palaeoclimate for Lateglacial and Holocene glacial events in New Zealand, and to better understand moraine-glacier-climate relationships. The climate reconstructions come from simulating past glacier expansions to specific terminal moraines, but I also simulated glacier fluctuations in response to a previously derived temperature reconstruction, and to interannual stochastic variability in temperature. The purpose behind each simulation was to identify the drivers of significant glacier fluctuations. The modelling results support the hypothesis that New Zealand moraine records reflect past climate, especially changes in temperature. Lateglacial climate was reconstructed to be 2-3 C lower than the present day. This temperature range agrees well with previous estimates from moraines and other climate proxy records in New Zealand. Modelled temperature estimates for the Holocene moraines are slightly colder than those derived from simpler methods, due to a non-linear relationship found between snowline lowering and glacier length. This relationship results from the specific valley shape and glacier geometry, and is likely to occur in other, similarly-shaped glacier valleys. The simulations forced by interannual stochastic variability in temperature do not show significant (>300 m) fluctuations in the glacier terminus. Such fluctuations can not explain the Holocene moraine sequence that I examined, which extends >2 km beyond the present-day glacier terminus. Stochastic temperature change could, however, in part, cause fluctuations in glacier extent during an overall glacier recession. Modelling shows that it is also unlikely that glaciers advanced to Holocene and Lateglacial moraine positions as a result of precipitation changes alone. For these reasons, temperature changes are a necessary part of explaining past glacier extents, especially during the Lateglacial, and the moraines examined here likely reflect changes in mean climate in New Zealand. The glacier modelling studies indicate that simpler methods, such as the accumulation area ratio, can be used to appropriately reconstruct past climate from glacial evidence, as long as the glacier catchment has a straight forward geometry, shallow bed slope and no tributary glaciers. Non-linear relationships between climate change and glacier length develop when valley shape is more complex, and glaciers within these systems are probably better simulated using a modelling approach. Using a numerical modelling approach, it is also possible to gain a greater understanding of glacier response time, length sensitivities, and estimates of ice extent in valleys within the model domain where geomorphic evidence is not available. In this manner, numerical models can be used as a tool for understanding past climate and glacier sensitivity, thus improving the confidence in the palaeoclimate interpretations.</p>

TermPicks: A century of Greenland glacier terminus data for use in machine learning applications

Marine-terminating outlet glacier terminus traces, mapped from satellite and aerial imagery, have been used extensively in understanding how outlet glaciers adjust to climate change variability over a range of time scales. Numerous studies have digitized termini manually, but this process is labor-intensive, and no consistent approach exists. A lack of coordination leads to duplication of efforts, particularly for Greenland, which is a major scientific research focus. At the same time, machine learning techniques are rapidly making progress in their ability to automate accurate extraction of glacier termini, with promising developments across a number of optical and SAR satellite sensors. These techniques rely on high quality, manually digitized terminus traces to be used as training data for robust automatic traces. Here we present a database of manually digitized terminus traces for machine learning and scientific applications. These data have been collected, cleaned, assigned with appropriate metadata including image scenes, and compiled so they can be easily accessed by scientists. The TermPicks data set includes 39,060 individual terminus traces for 278 glaciers with a mean and median number of traces per glacier of 136 ± 190 and 93, respectively. Across all glaciers, 32,567 dates have been picked, of which 4,467 have traces from more than one author (duplication of 14 %). We find a median error of ∼100 m among manually-traced termini. Most traces are obtained after 1999, when Landsat 7 was launched. We also provide an overview of an updated version of The Google Earth Engine Digitization Tool (GEEDiT), which has been developed specifically for future manual picking of the Greenland Ice Sheet.

Prokaryotic Community Succession in Bulk and Rhizosphere Soils Along a High-Elevation Glacier Retreat Chronosequence on the Tibetan Plateau

Early colonization and succession of soil microbial communities are essential for soil development and nutrient accumulation. Herein we focused on the changes in pioneer prokaryotic communities in rhizosphere and bulk soils along the high-elevation glacier retreat chronosequence, the northern Himalayas, Tibetan Plateau. Rhizosphere soils showed substantially higher levels of total organic carbon, total nitrogen, ammonium, and nitrate than bulk soils. The dominant prokaryotes were Proteobacteria, Actinobacteria, Acidobacteria, Chloroflexi, Crenarchaeota, Bacteroidetes, and Planctomycetes, which totally accounted for more than 75% in relative abundance. The dominant genus Candidatus Nitrososphaera occurred at each stage of the microbial succession. The richness and evenness of soil prokaryotes displayed mild succession along chronosequene. Linear discriminant analysis effect size (LEfSe) analysis demonstrated that Proteobacteria (especially Alphaproteobacteria) and Actinobacteria were significantly enriched in rhizosphere soils compared with bulk soils. Actinobacteria, SHA_109, and Thermoleophilia; Betaproteobacteria and OP1.MSBL6; and Planctomycetia and Verrucomicrobia were separately enriched at each of the three sample sites. The compositions of prokaryotic communities were substantially changed with bulk and rhizosphere soils and sampling sites, indicating that the communities were dominantly driven by plants and habitat-specific effects in the deglaciated soils. Additionally, the distance to the glacier terminus also played a significant role in driving the change of prokaryotic communities in both bulk and rhizosphere soils. Soil C/N ratio exhibited a greater effect on prokaryotic communities in bulk soils than rhizosphere soils. These results indicate that plants, habitat, and glacier retreat chronosequence collectively control prokaryotic community composition and succession.

Layered seawater intrusion and melt under grounded ice

Increasing melt of ice sheets at their floating or vertical interface with the ocean is a major driver of marine ice sheet retreat and sea level rise. However, the extent to which warm, salty seawater may drive melting under the grounded portions of ice sheets is still not well understood. Previous work has explored the possibility that dense seawater intrudes beneath relatively light subglacial freshwater discharge, similar to the salt wedge observed in many estuarine systems. In this study, we develop a generalized theory of layered seawater intrusion under grounded ice, including where subglacial hydrology occurs as a macroporous water sheet over impermeable beds or as microporous Darcy flow through permeable till. Using predictions from this theory, we show that seawater intrusion over hard beds may feasibly occur up to tens of kilometers upstream of a glacier terminus or grounding line. On the other hand, seawater is unlikely to intrude more than tens of meters through subglacial till. High-resolution simulations using the Ice-Sheet and Sea-Level System Model (ISSM) show that even just a few hundred meters of basal melt caused by seawater intrusion upstream of marine ice sheet grounding lines can cause projections of marine ice sheet volume loss to be 10–50 % higher or 100 % higher for kilometers of intrusion-induced basal melt. These results suggest that further observational, experimental and numerical investigations are needed to determine whether the conditions under which extensive seawater intrusion occurs and whether it will indeed drive rapid marine ice sheet retreat and sea level rise in the future.