Glacier Dynamics Research Articles

Understanding drivers of glacier-length variability over the last millennium

Abstract. Changes in glacier length reflect the integrated response to local fluctuations in temperature and precipitation resulting from both external forcing (e.g., volcanic eruptions or anthropogenic CO2) and internal climate variability. In order to interpret the climate history reflected in the glacier moraine record, the influence of both sources of climate variability must therefore be considered. Here we study the last millennium of glacier-length variability across the globe using a simple dynamic glacier model, which we force with temperature and precipitation time series from a 13-member ensemble of simulations from a global climate model. The ensemble allows us to quantify the contributions to glacier-length variability from external forcing (given by the ensemble mean) and internal variability (given by the ensemble spread). Within this framework, we find that internal variability is the predominant source of length fluctuations for glaciers with a shorter response time (less than a few decades). However, for glaciers with longer response timescales (more than a few decades) external forcing has a greater influence than internal variability. We further find that external forcing also dominates when the response of glaciers from widely separated regions is averaged. Single-forcing simulations indicate that, for this climate model, most of the forced response over the last millennium, pre-anthropogenic warming, has been driven by global-scale temperature change associated with volcanic aerosols.

Rock glacier inventory, permafrost probability distribution modeling and associated hazards in the Hunza River Basin, Western Karakoram, Pakistan

The destabilization of rock glaciers and permafrost variations is of great importance to the safety of the population and infrastructure in the Karakoram region because of their effects on land stability and river obstructions. In this study, we compiled the first complete rock glacier inventory for the Hunza Basin, western Karakoram, of 616 rock glaciers with an area of 194 km2 between 2800 and 5700 m a.s.l. We categorized the rock glaciers as intact or relict, and their distributions and destabilization were further analyzed and used along with in situ climate and elevation dataset to model the permafrost probability distribution. The modeled areas where the permafrost zonation index (PZI) is 0.5–1.00 indicate that permafrost occurs over 85% of the catchment area and lies above 3525 m a.s.l., which closely matches the zero-degree isotherm of 3800 m a.s.l. Based on the sensitivity analysis of the independent variables, elevation is the most sensitive variable, followed by net radiation, for predicting the probabilities of the presence and absence of permafrost. The model distributions are quite precise, with median posterior areas under the curve of 0.98 and 0.96 for model training and testing, respectively. We analyzed the rock glacier destabilization for 68 rock glaciers that interacted with river channels, of which 50 blocked or diverted river channels. Destabilized rock glaciers can be closely linked to the 0 °C isotherm between 3400 and 4600 m a.s.l. The significant damage caused by periodic floods from the subsequent blockage of river channels by landslides can be attributed to variations in permafrost. Which demolished infrastructure, including a hydropower plant, suspension bridge and water supply system in Hassan-abad catchment. Quantification of rock glacier dynamics and permafrost in the region can further improve policies related to the reduction in disaster risk and mitigation of associated hazards.

Remote Detection of Surge-Related Glacier Terminus Change across High Mountain Asia

High Mountain Asia (HMA) hosts the largest glacier concentration outside of polar regions. It is also distinct glaciologically as it forms one of two major surge clusters globally, and many glaciers there contradict the globally observed glacier recession trend. Surging glaciers are critical to our understanding of HMA glacier dynamics, threshold behaviour and flow instability, and hence have been the subject of extensive research, yet many dynamical uncertainties remain. Using the cloud-based geospatial data platform, Google Earth Engine (GEE) and GEE-developed tool, GEEDiT, to identify and quantify trends in the distribution and phenomenological characteristics of surging glaciers synoptically across HMA, we identified 137 glaciers as surging between 1987–2019. Of these, 55 were newly identified, 15 glaciers underwent repeat surges, and 18 were identified with enhanced glaciological hazard potential, most notably from Glacier Lake Outburst Floods (GLOFs). Terminus position time series analysis from 1987–2019 facilitated the development of a six-part phenomenological classification of glacier behaviour, as well as quantification of surge variables including active phase duration, terminus advance distance and rate, and surge periodicity. This research demonstrates the application of remote sensing techniques and the GEE platform to develop our understanding of surging glacier distribution and terminus phenomenology across large areas, as well as their ability to highlight potential geohazard locations, which can subsequently be used to focus monitoring efforts.

A theory of glacier dynamics and instabilities Part 1: Topographically confined glaciers

AbstractWe present a theoretical framework that integrates the dynamics of glaciers with and without the topographic confinement. This Part 1 paper concerns the former, which may exhibit surge cycles when subjected to thermal switches associated with the bed condition. With the topographic trough setting the glacier width and curbing the lateral drainage of the meltwater, the problem falls under the purview of the undrained plastic bed (UPB) formalism. Employing the UPB, we shall examine the external controls of the glacial behavior and test them against observations. Through our non-dimensionalization scheme, we construct a 2-D regime diagram, which allows a ready prognosis of the glacial properties over the full range of the external conditions, both climate- and size-related. We first discern the boundaries separating the glacial regimes of steady-creep, cyclic-surging and steady-sliding. We then apply the regime diagram to observed glaciers for quantitative comparisons. These include the Svalbard glaciers of both normal and surge types, Northeast Greenland Ice Stream characterized by steady-sliding, and Hudson Strait Ice Stream exhibiting cyclic surges. The quantitative validation of our model containing no free parameters suggests that the thermal switch may unify the dynamics of these diverse glaciers.

Automatic delineation of glacier grounding lines in differential interferometric synthetic-aperture radar data using deep learning

Delineating the grounding line of marine-terminating glaciers—where ice starts to become afloat in ocean waters—is crucial for measuring and understanding ice sheet mass balance, glacier dynamics, and their contributions to sea level rise. This task has been previously done using time-consuming, mostly-manual digitizations of differential interferometric synthetic-aperture radar interferograms by human experts. This approach is no longer viable with a fast-growing set of satellite observations and the need to establish time series over entire continents with quantified uncertainties. We present a fully-convolutional neural network with parallel atrous convolutional layers and asymmetric encoder/decoder components that automatically delineates grounding lines at a large scale, efficiently, and accompanied by uncertainty estimates. Our procedure detects grounding lines within 232 m in 100-m posting interferograms, which is comparable to the performance achieved by human experts. We also find value in the machine learning approach in situations that even challenge human experts. We use this approach to map the tidal-induced variability in grounding line position around Antarctica in 22,935 interferograms from year 2018. Along the Getz Ice Shelf, in West Antarctica, we demonstrate that grounding zones are one order magnitude (13.3 ± 3.9) wider than expected from hydrostatic equilibrium, which justifies the need to map grounding lines repeatedly and comprehensively to inform numerical models.

Seismic observations of a complex firn structure across the Amery Ice Shelf, East Antarctica

AbstractWe use seismic refraction data to investigate the firn structure across a suture zone on the Amery Ice Shelf, East Antarctica, and the possible role of glacier dynamics in firn evolution. In the downstream direction, the data reveal decreasing compressional-wave velocities and increasing penetration depth of the propagating wave in the firn layer, consistent with$\sim$1 m firn thickening every 6 km. The boundary between the Lambert Glacier unit to the west and a major suture zone and the Mawson Escarpment Ice Stream unit to the east, is marked by differences in firn thicknesses, compressional-wave velocities and seismic anisotropy in the across-flow direction. The latter does not contradict the presence of a single-maximum crystal orientation fabric oriented 45–$90^{\circ }$away from the flow direction. This is consistent with the presence of transverse simple shear governing the region's underlying ice flow regime, in association with elevated strain along the suture zone. The confirmation and quantification of the implied dynamic coupling between firn and the underlying ice requires integration of future seismic refraction, coring and modelling studies. Because firn is estimated to cover$\sim$98% of the Antarctic continent any such coupling may have widespread relevance to ice-sheet evolution and flow.

Seasonal ice flow velocity variations of Polar Record Glacier, East Antarctica during 2016–2019 using Sentinel-1 data

Remote sensing-based investigation of ice flow dynamics of Polar Record Glacier (PRG) during December–April of 2016–2019 has been conducted in this work. Using offset tracking method on the Sentinel-1 Synthetic Aperture Radar (SAR) images, we estimated the glacier ice flow velocity. Ice flow velocity near the glacier terminus indicated higher velocity during January and subsequently showed lower values in December during 2016–2019. The maximum and minimum velocity of the glacier was found to be ∼2.54 and 0.03 m/d for 2016–2017, ∼2.49 and 0.03 m/d for 2017–2018, and ∼2.47 and 0.03 m/d for 2018–2019. Results indicate that the maximum velocity occurred at the terminus of the glacier and minimum flow was detected on the ice sheet portion of the glacier. The average ice front position receded from ∼900 m in 2016–2017 to ∼650 m in 2017–2018 to ∼200 m in 2018–2019. A slowdown near the glacier terminus is observed for the period 2018–2019, while the total glacier advancement was observed to be 28 km2 from 877 km2 to 905 km2 during the period of 2016–2019. The enhanced velocity at the terminus in the early winter (March and April) was attributed to ice-free water surrounding the PRG terminus which eliminated the buttressing effect.

Greenland-wide inventory of ice marginal lakes using a multi-method approach

Ice marginal lakes are a dynamic component of terrestrial meltwater storage at the margin of the Greenland Ice Sheet. Despite their significance to the sea level budget, local flood hazards and bigeochemical fluxes, there is a lack of Greenland-wide research into ice marginal lakes. Here, a detailed multi-sensor inventory of Greenland’s ice marginal lakes is presented based on three well-established detection methods to form a unified remote sensing approach. The inventory consists of 3347 (pm 8%) ice marginal lakes (>0.05,{{text{ km }}^{2}}) detected for the year 2017. The greatest proportion of lakes lie around Greenland’s ice caps and mountain glaciers, and the southwest margin of the ice sheet. Through comparison to previous studies, a sim 75% increase in lake frequency is evident over the west margin of the ice sheet since 1985. This suggests it is becoming increasingly important to include ice marginal lakes in future sea level projections, where these lakes will form a dynamic storage of meltwater that can influence outlet glacier dynamics. Comparison to existing global glacial lake inventories demonstrate that up to 56% of ice marginal lakes could be unaccounted for in global estimates of ice marginal lake change, likely due to the reliance on a single lake detection method.

First ice thickness measurements in Tierra del Fuego at Schiaparelli Glacier, Chile

Abstract. Cordillera Darwin in Tierra del Fuego (Chile) remains one of the least studied glaciated regions in the world. However, this region being one of very few terrestrial sites at this latitude in the Southern Hemisphere has the potential to provide key information on the effect of climate variability and climate change on the cryosphere at sub-polar mid-latitudes of the Southern Hemisphere. Schiaparelli Glacier is located at the northern side of the Cordillera Darwin draining the north side of Monte Sarmiento (2187 m a.s.l.). Despite being one of the largest glaciers in the Cordillera Darwin, no previous in situ observation of its ice thickness had been made either at this glacier or at any other location in the Cordillera Darwin. Ice thickness is one of the fundamental parameters to understand glacier dynamics, constrain ice dynamical modelling, and predict glacier evolution. In April 2016 we performed the first successful ice thickness measurements using terrestrial ground-penetrating radar in the ablation area of Schiaparelli Glacier (Gacitúa et al., 2020, https://doi.org/10.1594/PANGAEA.919331). The measurements were made along a transect line perpendicular to the ice flow. Results show a valley-shaped bedrock with a maximum ice thickness of 324 m within a distinct glacier trough. The bedrock is located below current sea level for 51 % of the transect measurements with a minimum of −158 m, which illustrates that the local topography is subject to considerable glacier-related over-deepening.

Rockwall Slope Erosion in the Northwestern Himalaya

AbstractRockwall slope erosion is an important component of alpine landscape evolution, yet the role of climate and tectonics in driving this erosion remains unclear. We define the distribution and magnitude of periglacial rockwall slope erosion across 12 catchments in Himachal Pradesh and Jammu and Kashmir in the Himalaya of northern India using cosmogenic 10Be concentrations in sediment from medial moraines. Beryllium‐10 concentrations range from 0.5 ± 0.04 × 104 to 260.0 ± 12.5 × 104 at/g SiO2, which yield erosion rates between 7.6 ± 1.0 and 0.02 ± 0.004 mm/a. Between ∼0.02 and ∼8 m of rockwall slope erosion would be possible in this setting across a single millennium, and >2 km when extrapolated for the Quaternary period. This erosion affects catchment sediment flux and glacier dynamics, and helps to establish the pace of topographic change at the headwaters of catchments. We combine rockwall erosion records from the Himalaya of Himachal Pradesh, Jammu and Kashmir, and Uttarakhand in India and Baltistan in Pakistan to create a regional erosion data set. Rockwall slope erosion rates progressively decrease with distance north from the Main Central Thrust and into the interior of the orogen. The distribution and magnitude of this erosion is most closely associated with records of Himalayan denudation and rock uplift, where the highest rates of change are recorded in the Greater Himalaya sequences. This suggests that tectonically driven uplift, rather than climate, is a first order control on patterns of rockwall slope erosion in the northwestern Himalaya. Precipitation and temperature would therefore come as secondary controls.

Quantitative estimates of Holocene glacier meltwater variations on the Western Tibetan Plateau

Knowledge of the alpine glacier meltwater variations is fundamental prerequisite for understanding glacier dynamics and assessing the availability of freshwater resources. Glaciers on the Tibetan Plateau (TP) are sources of water for most major Asian rivers, but their melting history remains unclear, preventing in-depth understanding of their mechanisms. Here, we propose the authigenic carbonate δ18O from glacial lakes as a quantitative proxy to estimate variations of glacier meltwater. In the Western Kunlun Mountain, δ18O record at Guozha Co indicates that maximum glacier meltwater (−28.62±25.76Gt) occurred at 9.5–8.5 ka BP, and minimum glacier meltwater (24.53±25.02Gt) at 1.3–0.5 ka BP. Nearly 20% of regional glaciers melted from the Early to Late Holocene, likely controlled by the summer temperature and accumulation of melting potential estimated by positive degree-day. Based on the projected temperature, this study suggests the TP glaciers likely face severe threats at the current rates of global warming.

Longbasaba Glacier recession and contribution to its proglacial lake volume between 1988 and 2018

AbstractDuring the last few decades, the lake-terminating glaciers in the Himalaya have receded faster than the land-terminating glaciers as proglacial lakes have exacerbated the mass loss of their host glaciers. Monitoring the impacts of glacier recession and dynamics on lake extent and water volume provides an approach to assess the mass interplay between glaciers and proglacial lakes. We describe the recession of Longbasaba Glacier and estimate the mass wastage and its contribution to the water volume of its proglacial lake. The results show that the glacier area has decreased by 3% during 1988–2018, with a more variable recession prior to 2008 than in the last decade. Longbasaba Lake has expanded by 164% in area and 237% in water volume, primarily as a result of meltwater inflow produced from surface lowering of the glacier. Over the periods 1988–2000 and 2000–18, the mass loss contributed by glacier thinning has decreased from 81 to 61% of the total mass loss, accompanied by a nearly doubled contribution from terminus retreat. With the current rate of retreat, Longbasaba glacier is expected to terminate in its proglacial lake for another four decades. The hazard risk of this lake is expected to continue to increase in the near future because of the projected continued glacier mass loss and related lake expansion.

Observing traveling waves in glaciers with remote sensing: new flexible time series methods and application to Sermeq Kujalleq (Jakobshavn Isbræ), Greenland

Abstract. The recent influx of remote sensing data provides new opportunities for quantifying spatiotemporal variations in glacier surface velocity and elevation fields. Here, we introduce a flexible time series reconstruction and decomposition technique for forming continuous, time-dependent surface velocity and elevation fields from discontinuous data and partitioning these time series into short- and long-term variations. The time series reconstruction consists of a sparsity-regularized least-squares regression for modeling time series as a linear combination of generic basis functions of multiple temporal scales, allowing us to capture complex variations in the data using simple functions. We apply this method to the multitemporal evolution of Sermeq Kujalleq (Jakobshavn Isbræ), Greenland. Using 555 ice velocity maps generated by the Greenland Ice Mapping Project and covering the period 2009–2019, we show that the amplification in seasonal velocity variations in 2012–2016 was coincident with a longer-term speedup initiating in 2012. Similarly, the reduction in post-2017 seasonal velocity variations was coincident with a longer-term slowdown initiating around 2017. To understand how these perturbations propagate through the glacier, we introduce an approach for quantifying the spatially varying and frequency-dependent phase velocities and attenuation length scales of the resulting traveling waves. We hypothesize that these traveling waves are predominantly kinematic waves based on their long periods, coincident changes in surface velocity and elevation, and connection with variations in the terminus position. This ability to quantify wave propagation enables an entirely new framework for studying glacier dynamics using remote sensing data.

Evolution of debris cover on glaciers of the Eastern Alps, Austria, between 1996 and 2015

AbstractDebris cover on glaciers is an important component of glacial systems as it influences climate–glacier dynamics and thus the lifespan of glaciers. Increasing air temperatures, permafrost thaw and rock faces freshly exposed by glacier downwasting in accumulation zones result in increased rockfall activity and debris input. In the ablation zone, negative mass balances result in an enhanced melt‐out of englacial debris. Glacier debris cover thus represents a clear signal of climate warming in mountain areas. To assess the temporal development of debris on glaciers of the Eastern Alps, Austria, we mapped debris cover on 255 glaciers using Landsat data at three time steps. We applied a ratio‐based threshold classification technique and analysed glacier catchment characteristics to understand debris sources better. Across the Austrian Alps, debris cover increased by more than 10% between 1996 and 2015 while glaciers retreated in response to climate warming. Debris cover distribution shows significant regional variability, with some mountain ranges being characterised by mean debris cover on glaciers of up to 75%. We also observed a general rise of the mean elevation of debris cover on glaciers in Austria. The debris cover distribution and dynamics are highly variable due to topographic, lithological and structural settings that determine the amount of debris delivered to and stored in the glacier system. Despite strong variation in debris cover, all glaciers investigated melted at increasing rates. We conclude that the retarding effects of debris cover on the mass balance and melt rate of Austrian glaciers is strongly subdued compared with other mountain areas. The study indicates that, if this trend continues, many glaciers in Austria may become fully debris covered. However, since debris cover seems to have little impact on melt rates, this would not lead to prolonged existence of debris‐covered ice compared with clean ice glaciers.

Morphology and evolution of supraglacial hummocks on debris‐covered Himalayan glaciers

AbstractThick supraglacial debris layers often have an undulating, hummocky topography that influences the lateral transport of debris and meltwater and provides basins for supraglacial ponds. The role of ablation and other processes associated with supraglacial debris in giving rise to this hummocky topography is poorly understood. Characterizing hummocky topography is a first step towards understanding the feedbacks driving the evolution of debris‐covered glacier surfaces and their potential impacts on mass balance, hydrology and glacier dynamics. Here we undertake a geomorphological assessment of the hummocky topography on five debris‐covered glaciers in the Everest region of the central Himalaya. We characterize supraglacial hummocks through statistical analyses of their vertical relief and horizontal geometry. Our results establish supraglacial hummocks as a distinct landform. We find that a typical hummock has an elongation ratio of 1.1:1 in the direction of ice flow, length of 214 ± 109 m and width of 192 ± 88 m. Hummocky topography has a greater amplitude across‐glacier (15.4 ± 10.9 m) compared to along the glacier flow line (12.6 ± 8.3 m). Consequently, hummock slopes are steeper in the across‐glacier direction (8.7 ± 4.3°) than in the direction of ice flow (5.6 ± 4.0°). Longer, wider and higher‐amplitude hummocks are found on larger glaciers. We postulate that directional anisotropy in the hummock topography arises because, while the pattern of differential ablation driving topography evolution is moderated by processes including the gravitational redistribution of debris across the glacier surface, it also inherits an orientation preference from the distribution of englacial debris in the underlying ice. Our morphometric data inform future efforts to model these interactions, which should account for additional factors such as the genesis of supraglacial ponds and ice cliffs and their impact on differential ablation.

Drivers of Pine Island Glacier speed-up between 1996 and 2016

Abstract. Pine Island Glacier in West Antarctica is among the fastest changing glaciers worldwide. Over the last 2 decades, the glacier has lost in excess of a trillion tons of ice, or the equivalent of 3 mm of sea level rise. The ongoing changes are thought to have been triggered by ocean-induced thinning of its floating ice shelf, grounding line retreat, and the associated reduction in buttressing forces. However, other drivers of change, such as large-scale calving and changes in ice rheology and basal slipperiness, could play a vital, yet unquantified, role in controlling the ongoing and future evolution of the glacier. In addition, recent studies have shown that mechanical properties of the bed are key to explaining the observed speed-up. Here we used a combination of the latest remote sensing datasets between 1996 and 2016, data assimilation tools, and numerical perturbation experiments to quantify the relative importance of all processes in driving the recent changes in Pine Island Glacier dynamics. We show that (1) calving and ice shelf thinning have caused a comparable reduction in ice shelf buttressing over the past 2 decades; that (2) simulated changes in ice flow over a viscously deforming bed are only compatible with observations if large and widespread changes in ice viscosity and/or basal slipperiness are taken into account; and that (3) a spatially varying, predominantly plastic bed rheology can closely reproduce observed changes in flow without marked variations in ice-internal and basal properties. Our results demonstrate that, in addition to its evolving ice thickness, calving processes and a heterogeneous bed rheology play a key role in the contemporary evolution of Pine Island Glacier.

TanDEM-X:Deriving InSAR Height Changes and Velocity Dynamics of Great Aletsch Glacier

Great Aletsch Glacier (Grosser Aletschgletscher), the largest glacier in the European Alps, contains 20% of the entire Swiss ice mass. Therefore, it has been selected as a super-testsite for the TanDEM-X satellite mission. Dense time series with a repeat interval down to 11 d were acquired between 2011 and 2019 using two polarizations (HH and VV) and across-track baselines of 0-1 km. To evaluate the use of interferometric single-pass synthetic aperture radar mission for glaciological applications, we implemented a processing pipeline in interactive data language (IDL) and computed 130 digital elevation models (DEMs) from bistatic radar interferograms. We present a method to circumvent a common pitfall during orthorectification of radar DEM differences. Regression analysis of DEM time series shows a height loss of up to 8 m a$^{-1}$ on the tongue and 1.5 m a$^{-1}$ when averaged over the whole glacier area. In spring 2013, we observed X-band penetration depths of 4±2 m in the accumulation area. For strongly crevassed areas, the coherence drops already to 0.5 for across-track baselines <b xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">B</b> <sub xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">⊥</sub> > 200 m. With patch-based incoherent offset tracking, we obtained an almost complete 200 m resolution velocity map. Velocities reach up to 0.8 m d$^{-1}$, show a seasonal variability of ±0.05 m d$^{-1}$, and agree within 0.04 m d$^{-1}$(root mean square error) with field measurements. Copolar phase differences suggest an approximation of snow accumulation. We demonstrate orthorectification of the backscatter intensity using simultaneously acquired TanDEM-X interferograms, which allows for the decoupling of horizontal velocity estimates from phase-center height changes due to penetration and ice melt.

Оценка состояния малого оледенения в Хибинах

The presence of very small glaciation in the Khibiny mountains was discovered by V.F. Perov in 1958. Since then and until the beginning of the 21st study of glaciers in this area has not been carried out. In general, warming is recorded on the Kola Peninsula. The average annual increase in temperature is 2.3°C ± 1°C during the last 50 years, which correlates with the data of other researchers. The glaciers in the Khibiny mountains should have significantly reduced their size in this connection. In other Arctic regions degradation of glaciation is noted. Some researchers predicted to disappear of Khibiny glaciation. In order to assess the change in glacier areas, we carried out a comprehensive study of small glaciers in the Khibiny mountains. We studied dynamics of glaciers over the past 15 years using field observations and the use of GIS and remote sensing methods. In our study, we used aerial photographs of 1958, digital models built on the basis of field measurements, and satellite images for years different in meteorological parameters. Our research shows, despite the slight degradation of glaciation, this tendency is modest and glaciers’ area remains relatively stable. Our research shows that glaciers shrank significantly in the early 2000s and showed a tendency towards complete disappearance. Their area has decreased by more than two times compared to the figures obtained in 1958. This fact may be caused by the increase in snow precipitation through in recent years. After 2007, there has been an increase in snow precipitation. In subsequent years, starting in 2008, the glaciers began to regain their size. Over the past 15 years, the average area of two glaciers (No. 3 and No. 4) exceeds the size of 1958. Glaciers No. 1 and No. 2 have shrinkage.

High diversity of planktonic prokaryotes in Arctic Kongsfjorden seawaters in summer 2015

Kongsfjorden is a typical glacier fjord in the European Arctic. In order to know the distribution and diversity of planktonic prokaryotes in the fjord subjected to strong physical gradient related with local glacier dynamics, prokaryote communities along a transect extending from the outer to the inner fjord was investigated in the summer of 2015 using 454 pyrosequencing of archaeal and bacterial 16S rRNA genes. Higher diversity and species richness were detected in Bacteria than in Archaea in each sample. Archaeal 16S rRNA gene sequences were not detected in the outer area of the fjord. Nitrososphaeria within the phylum Crenarchaeota and Poseidoniia within the Thermoplasmatota were two dominant archaeal groups in the fjord. At the genus level, Nitrosopumilus, Nitrosopelagicus, MGIIa-L1, and MGIIb–O2 dominated seawater archaeal communities. Higher abundances of Poseidoniale-related archaea were observed in the deep water than in the surface water. The Proteobacteria, Bacteroidota, and Actinobacteriota dominated the planktonic bacterial communities, with Verrucomicrobiota (represented by Verrucomicrobiae) dominating the middle area of the fjord. At the class level, Bacteroidia, Alphaproteobacteria, Gammaproteobacteria, and Actinobacteria dominated the bacterial communities within Kongsfjorden. Sequences affiliated with genera Polaribacter and Yoonia were frequently detected in the fjord. Results suggest that planktonic archaeal and bacterial community compositions in Kongsfjorden can be driven by different major environmental factors.

Contribution of Crevasse Advection and Mixed Mode Calving to Glacier Dynamics

Contribution of Crevasse Advection and Mixed Mode Calving to Glacier Dynamics