Organic carbon accumulation in the glacier forelands with regard to variability of environmental conditions in different ecogenesis stages of High Arctic ecosystems

From glacier forelands to human settlements: Patterns, environmental drivers, and risks of antibiotic resistance genes

Antibiotic resistance genes (ARGs) are biological pollutants widely present in glaciers, such as ice, snow, and melt water. However, it remains unclear whether ARGs in glaciers influence their distribution in human settlements within the glacier basins. Therefore, we investigated the distribution pattern and driving factors of ARGs in the Laohugou glacier basins on the Tibetan Plateau. Using high-throughput quantitative PCR, the total abundance of ARGs in the Laohugou glacier basins ranged from 7.53 × 10⁶ to 1.83 × 10⁹ copies/g, including 128 detected ARGs across 11 classes, with aminoglycoside resistance genes being the dominant group. The abundance of ARGs exhibited a U-shaped pattern along the elevational gradient, with higher levels in glacier regions and human settlements, and the lowest abundance at mid-elevations. While glacier melting and anthropogenic disturbance are likely major contributors to this pattern, other potential mechanisms may also be involved, such as elevation-dependent microbial community composition, atmospheric deposition and release of legacy ARGs from melting permafrost and glacial ice. Together, these processes likely interact to shape the observed ARG pattern in this alpine watershed. We further verified that the distribution of ARGs was strongly correlated with microbial community structure, especially bacterial communities (r > 0.50; p < 0.05). Network analysis showed that Nitrolancea negatively correlated with several core ARGs, suggesting its potential role in regulating the spread of ARGs. Random forest analysis and structural equation modeling (SEM) indicated that, after accounting for various driving factors, organic matter and bacterial biomass were the primary drivers of increased ARG abundance. This study provides a foundation for assessing the risks of ARGs in glacier basins under global climate change, offering insights into risk mitigation strategies and guiding future ecological and public health research.

Dynamics and drivers of mycorrhizal fungi after glacier retreat

The development of terrestrial ecosystems depends greatly on plant mutualists such as mycorrhizal fungi. The global retreat of glaciers exposes nutrient-poor substrates in extreme environments and provides a unique opportunity to study early successions of mycorrhizal fungi by assessing their dynamics and drivers. We combined environmental DNA metabarcoding and measurements of local conditions to assess the succession of mycorrhizal communities during soil development in 46 glacier forelands around the globe, testing whether dynamics and drivers differ between mycorrhizal types. Mycorrhizal fungi colonized deglaciated areas very quickly (< 10 yr), with arbuscular mycorrhizal fungi tending to become more diverse through time compared to ectomycorrhizal fungi. Both alpha- and beta-diversity of arbuscular mycorrhizal fungi were significantly related to time since glacier retreat and plant communities, while microclimate and primary productivity were more important for ectomycorrhizal fungi. The richness and composition of mycorrhizal communities were also significantly explained by soil chemistry, highlighting the importance of microhabitat for community dynamics. The acceleration of ice melt and the modifications of microclimate forecasted by climate change scenarios are expected to impact the diversity of mycorrhizal partners. These changes could alter the interactions underlying biotic colonization and belowground-aboveground linkages, with multifaceted impacts on soil development and associated ecological processes.

Metagenomic insights into novel microbial lineages with distinct ecological functions in the Arctic glacier foreland ecosystems

Land-terminating glaciers are retreating globally, resulting in the expansion of the ice-free glacier forelands (GFs). These GFs act as a natural laboratory to study microbial community succession, soil formation, and ecosystem development. Here, we have employed gene-centric and genome-resolved metagenomic approaches to disseminate microbial diversity, community structure, and their associated biogeochemical processes involved in the carbon, nitrogen, and sulfur cycling across three GF ecosystems. Here, we present a compendium of draft Metagenome Assembled Genomes (MAGs) belonging to bacterial (n = 899) and archaeal (n = 4) domains. These MAGs were reconstructed using a total of 27 shotgun metagenomic datasets obtained from three different GFs, including Midtre Lovénbreen glacier (Svalbard), Russell glacier (Greenland), and Storglaciaren (Sweden). The taxonomic classification revealed that 98% of MAGs remained unclassified at species levels, suggesting the presence of novel microbial lineages. The abundance of metabolic genes associated with carbon, nitrogen, and sulfur cycling pathways varied between and within the samples collected across the three GF ecosystems. Our findings indicate that MAGs from different GFs share close phylogenetic relationships but exhibit significant differences in abundance, distribution patterns, and metabolic functions. This compendium of novel MAGs, encompassing autotrophic, phototrophic, and chemolithoautotrophic microbial groups reconstructed from GF ecosystems, represents a valuable resource for further studies.

Metabarcoding data reveal vertical multitaxa variation in topsoil communities during the colonization of deglaciated forelands

Ice-free areas are expanding worldwide due to dramatic glacier shrinkage and are undergoing rapid colonization by multiple lifeforms, thus representing key environments to study ecosystem development. It has been proposed that the colonization dynamics of deglaciated terrains is different between surface and deep soils but that the heterogeneity between communities inhabiting surface and deep soils decreases through time. Nevertheless, tests of this hypothesis remain scarce, and it is unclear whether patterns are consistent among different taxonomic groups. Here, we used environmental DNA metabarcoding to test whether community diversity and composition of six groups (Eukaryota, Bacteria, Mycota, Collembola, Insecta, and Oligochaeta) differ between the surface (0-5 cm) and deeper (7.5-20 cm) soil at different stages of development and across five Alpine glaciers. Taxonomic diversity increased with time since glacier retreat and with soil evolution. The pattern was consistent across groups and soil depths. For Eukaryota and Mycota, alpha-diversity was highest at the surface. Time since glacier retreat explained more variation of community composition than depth. Beta-diversity between surface and deep layers decreased with time since glacier retreat, supporting the hypothesis that the first 20 cm of soil tends to homogenize through time. Several molecular operational taxonomic units of bacteria and fungi were significant indicators of specific depths and/or soil development stages, confirming the strong functional variation of microbial communities through time and depth. The complexity of community patterns highlights the importance of integrating information from multiple taxonomic groups to unravel community variation in response to ongoing global changes.

Changes of Soil Dissolved Organic Matter and Its Relationship with Microbial Community along the Hailuogou Glacier Forefield Chronosequence

Glacier-retreated areas are ideal areas to study soil biogeochemical processes during vegetation succession, because of the limited effect of other environmental and climatic factors. In this study, the changes of soil dissolved organic matter (DOM) and its relationship with microbial communities along the Hailuogou Glacier forefield chronosequence were investigated. Both microbial diversity and DOM molecular chemodiversity recovered rapidly at the initial stage, indicating the pioneering role of microorganisms in soil formation and development. The chemical stability of soil organic matter enhanced with vegetation succession due to the retaining of compounds with high oxidation state and aromaticity. The molecular composition of DOM affected microbial communities, while microorganisms tended to utilize labile components to form refractory components. This complex relationship network between microorganisms and DOM components played an important role in the development of soil organic matter as well as the formation of stable soil carbon pool in glacier-retreated areas.

Vascular plant and cryptogam abundance as well as soil chemical properties shape microbial communities in the successional gradient of glacier foreland soils

In the glacier forelands, microbes play a fundamental role in soil development and shaping the vegetation structure. Such ecosystems represent various stages of soil development and are, therefore, an excellent place to study the interrelationship between soil, plants, and microorganisms. The aim of the study was to assess the effects of vegetation and soil physicochemical properties developing after glacier retreat on soil microbial communities. Specifically, abundance, species richness and the composition of arbuscular mycorrhizal fungi (AMF), as well as microbial biomass and community structure in soils were compared between plots established in 800-meter transects of three glacier forelands in northern Sweden. The cover of vascular plants and cryptogams, soil C content, AMF spore density and species richness, AMF biomass indicators, total microbial biomass, and bacterial phospholipid fatty acids (PLFA) were significantly and positively related to the distance from the glacier terminus. On the other hand, macronutrient concentrations and pH decreased along with increasing distance. No significant impact of the distance from the glacier terminus on the ratio fungal/bacterial PLFA was observed. Moreover, we found a significant effect of both glacier and the distance from the glacier terminus on the microbial community structure. AMF species richness and spore density in the glacier forelands were generally low, which is probably due to a limited supply of inoculum in primary successional ecosystems. Most microbial biochemical markers and AMF parameters were positively associated with the number of arbuscular mycorrhizal plant species and vascular plant and lichen cover as well as C content in soil, whereas negatively with soil macronutrients and pH. This could be related to an increase in plant cover and a decrease in soil nutrient levels as plant succession progresses. Our results showed that vegetation, soil C content, and microbial communities are interlinked and exhibit concordant patterns along successional gradients.

The crucial interactions between climate and soil

Since the birth of soil science, climate has been recognized as a soil-forming factor, along with parent rock, time, topography, and organisms (from which humans were later kept distinct), often prevalent on the other factors on the very long term. But the climate is in turns affected by soils and their management. This paper describes the interrelationships between climate - and its current change - and soil, focusing on each single factor of its formation. Parent material governs, primarily through the particle size distribution, the capacity of soil to retain water and organic matter, which are two main soil-related drivers of the climate. Time is the only unmanageable soil-forming factor; however, extreme climatic phenomena can upset the soil or even dismantle it, so as to slow down the pathway of pedogenesis or even make it start from scratch. Topography, which drives the pedogenesis mostly controlling rainfall distribution - with repercussions also on the climate - is not anymore a given factor because humans have often become a shaper of it. Indeed humans now play a key role in affecting in a plethora of ways those soil properties that most deal with climate. The abundance and diversity of the other organisms are generally positive to soil quality and as a buffer for climate, but there are troubling evidences that climate change is decreasing soil biodiversity. The corpus of researches on mutual feedback between climate and soil has essentially demonstrated that the best soil management in terms of climate change mitigation must aim at promoting vegetation growth and maximizing soil organic matter content and water retention. Some ongoing virtuous initiatives (e.g., the Great Green Wall of Africa) and farming systems (e.g., the conservation agriculture) should be extended as much as possible worldwide to enable the soil to make the greatest contribution to climate change mitigation.

Arctic methylmercury cycling

Anthropogenic mercury (Hg) undergoes long-range transport to the Arctic where some of it is transformed into methylmercury (MeHg), potentially leading to high exposure in some Arctic inhabitants and wildlife. The environmental exposure of Hg is determined not just by the amount of Hg entering the Arctic, but also by biogeochemical and ecological processes occurring in the Arctic. These processes affect MeHg uptake in biota by regulating the bioavailability, methylation and demethylation, bioaccumulation and biomagnification of MeHg in Arctic ecosystems. Here, we present a new budget for pools and fluxes of MeHg in the Arctic and review the scientific advances made in the last decade on processes leading to environmental exposure to Hg. Methylation and demethylation are key processes controlling the pool of MeHg available for bioaccumulation. Methylation of Hg occurs in diverse Arctic environments including permafrost, sediments and the ocean water column, and is primarily a process carried out by microorganisms. While microorganisms carrying the hgcAB gene pair (responsible for Hg methylation) have been identified in Arctic soils and thawing permafrost, the formation pathway of MeHg in oxic marine waters remains less clear. Hotspots for methylation of Hg in terrestrial environments include thermokarst wetlands, ponds and lakes. The shallow sub-surface enrichment of MeHg in the Arctic Ocean, in comparison to other marine systems, is a possible explanation for high MeHg concentrations in some Arctic biota. Bioconcentration of aqueous MeHg in bacteria and algae is a critical step in the transfer of Hg to top predators, which may be dampened or enhanced by the presence of organic matter. Variable trophic position has an important influence on MeHg concentrations among populations of top predator species such as ringed seal and polar bears distributed across the circumpolar Arctic. These scientific advances highlight key processes that affect the fate of anthropogenic Hg deposited to Arctic environments.

Environmental micro-niche filtering shapes bacterial pioneer communities during primary colonization of a Himalayas' glacier forefield

The pedogenesis from the mineral substrate released upon glacier melting has been explained with the succession of consortia of pioneer microorganisms, whose structure and functionality are determined by the environmental conditions developing in the moraine. However, the microbiome variability that can be expected in the environmentally heterogeneous niches occurring in a moraine at a given successional stage is poorly investigated. In a 50 m² area in the forefield of the Lobuche glacier (Himalayas, 5050 m above sea level), we studied six sites of primary colonization presenting different topographical features (orientation, elevation and slope) and harbouring greyish/dark biological soil crusts (BSCs). The spatial vicinity of the sites opposed to their topographical differences, allowed us to examine the effect of environmental conditions independently from the time of deglaciation. The bacterial microbiome diversity and their co-occurrence network, the bacterial metabolisms predicted from 16S rRNA gene high-throughput sequencing, and the microbiome intact polar lipids were investigated in the BSCs and the underlying sediment deep layers (DLs). Different bacterial microbiomes inhabited the BSCs and the DLs, and their composition varied among sites, indicating a niche-specific role of the micro-environmental conditions in the bacterial communities' assembly. In the heterogeneous sediments of glacier moraines, physico-chemical and micro-climatic variations at the site-spatial scale are crucial in shaping the microbiome microvariability and structuring the pioneer bacterial communities during pedogenesis.

SEM-EDS and water chemistry characteristics at the early stages of glacier recession reveal biogeochemical coupling between proglacial sediments and meltwater

Most glaciers worldwide are undergoing climate-forced recession, but the impact of glacier changes on biogeochemical cycles is unclear. This study examines the influence of proglacial sediment weathering on meltwater chemistry at the early stages of glacier recession in the High Arctic of Svalbard. Scanning electron microscopy-energy dispersive X-ray spectroscopy (SEM-EDS) in combination with a wide range of geochemical analyses were used in this study. The SEM-EDS analyses of sediments collected in front of Werenskioldbreen show general degradation of pyrite and carbonate grains with age. The outer parts of pyrite grains have a gradual decrease in sulphur and gradual increase in iron oxides due to pyrite oxidation. This process was less advanced in the proglacial zone younger than 100 years compared to older sites such as the terminal moraine from the Little Ice Age. In both the proglacial zone and the terminal moraine, physical weathering of mineral grains, including formation of microcracks and microfractures, clearly enhanced pyrite oxidation. A consequence of proglacial sediment weathering is that the river chemistry is strongly affected by carbonate dissolution driven by sulphuric acid from sulphide oxidation. Also, reactive iron oxides, a product of sulphide oxidation, are mobilized in the proglacial zone. The results of this study show that proglacial weathering in the High Arctic of Svalbard is strongly coupled to river geochemistry, especially during the early stages of proglacial exposure after glacier recession.

Cryosphere Microbiome Biobanks for Mountain Glaciers in China

The glaciers in China have an important role as one of the most climate-sensitive constituents of the Tibetan Plateau which is known as the Asian Water Tower. Although the cryosphere is one of the most extreme environments for organisms, the soils of the glacier foreland harbor surprisingly rich microbiomes. A large amount of accelerated glacier retreat accompanied by global warming will not only raise the sea level, but it will also lead to the massive release of a considerable amount of carbon stored in these glaciers. The responses of glacier microbiomes could alter the biogeochemical cycle of carbon and have a complex impact on climate change. Thus, understanding present-day and future glacier microbiome changes is crucial to assess the feedback on climate change and the impacts on ecosystems. To this end, we discuss here the diversity and biogeochemical functions of the microbiomes in Chinese mountain glacier ecosystems.

Seasonal Changes in the Photosynthetic Activity of Terrestrial Lichens and Mosses in the Lichen Scots Pine Forest Habitat

Photosynthetic activity is one of the most important metabolic processes that can be quickly and easily studied in the field. It can be used for identifying the environmental factors affecting ecosystem balance, as any stressor influencing metabolic and physiological processes will have a measurable effect on photosynthesis. The aim of this study was to measure the photosynthetic activity of selected lichens and mosses and investigate its changes resulted from diurnal and seasonal variability. We studied two lichens (Cladonia mitis Sandst and Cladonia uncialis (L.) Weber ex F.H. Wigg.) and two mosses (Pleurozium schreberi (Willd. ex Brid.) Mitt. and Dicranum scoparium (L.) Hedw.). Samples were collected in the area of lichen Scots pine forest of the “Bory Tucholskie” National Park. Our study revealed that the photosynthetic activity of cryptogams depended on species, season, time of the day, and water availability. Cladonia species, which are the main component of lichen Scots pine forests, have higher photosynthetic activity than Pleurozium schreberi, which represents species of fresh coniferous forests. Photosynthetic activity increased from spring through summer and reached the highest values in autumn. It was also higher in soaked samples collected in the morning and afternoon compared to noon. Despite the water access, noon samples still showed the lowest activity. This can result from natural changes in humidity during the day to which cryptogams are well-adapted.

Flux of Polynuclear Aromatic Compounds (PAHs) from the Atmosphere and from Reindeer/Bird Feces to Arctic Soils in Ny-Ålesund (Svalbard)

Atmospheric, soil, and feces samples were collected in Ny-Ålesund during July 2015. The concentrations, distributions, congener profiles, and contaminant migration levels were analyzed for 16 polycyclic aromatic hydrocarbons (PAHs) identified by the United States Environmental Protection Agency (US EPA) as priority contaminants (16 PAHs). Mean concentrations in the gas and particle phases were 37.8 ng m^-3 and 2.9 ng m^-3, respectively, and mean concentrations in soil and reindeer/bird feces were 329.1 ng g^-1 and 720.7 ng g^-1, respectively, on a dry weight (d.w.) basis. In more than three phases, naphthalene and phenanthrene dominated the concentrations of the 16 total PAHs (Σ₁₆PAH) and the concentration of PAHs in the gas phase was much higher than in the particle phase. The main sources of local PAHs may be coal combustion and air-surface exchange. There was a volatilization tendency from soil to air for 2-4 ring PAHs, and exchange fluxes were ~ 10⁵ times greater than the deposition fluxes of 5-6 ring PAHs. The underground migration of PAHs was investigated in Ny-Ålesund; the results showed flux values of ~ 0.07% from the initial PAH concentrations.

Patterns and drivers of cryptogam and vascular plant diversity in glacier forelands

Vascular and nonvascular plants are affected by environmental factors determining their distribution and shaping their diversity and cover. Despite the cryptogam commonness in Arctic communities, previous studies have often focused on limited number of factors and their impact on only selected species of vascular plants or cryptogams. Our study aimed to investigate in detail the differences in species diversity and cover of cryptogams and vascular plants in the glacier forelands and mature tundra on Svalbard. Furthermore, we determined the biotic and abiotic factors that affected diversity, cover and distribution of cryptogam and vascular plant species. In 2017, we established 201 plots in eight locations (each including habitat type of foreland and mature tundra) and surveyed species abundance, sampled soils and environmental data. Results revealed that diversity and cover of analysed groups differed significantly between locations and habitat types, except for cryptogam cover in mature tundra in terms of location. Distance to the glacier terminus, slope, soil conductivity, nutrient content, and clay content impacted both plant groups' diversity. In contrast, distance to the glacier terminus, nutrient content and soil pH affected their cover. In addition, for cryptogam diversity and cover, foreland location and vascular plant cover were also important, while for vascular plant cover time elapsed after glacier retreat was significant. Distribution of both groups' species in forelands was associated with time elapsed after glacier retreat, soil pH, and nutrient contents. Soil texture and distance to the glacier terminus additionally influenced cryptogam distribution. The positive impact of vascular plants on cryptogam diversity and cover indicates complex relationships between these groups, even in forelands' relatively simple communities. As the cryptogam diversity in the polar areas is high but still largely unknown, future studies on species ecology and climate change impact on vegetation should consider both vascular plants and cryptogams and interactions between these groups.

Soil development and spatial differentiation in a glacial river valley under cold and extremely arid climate of East Pamir Mountains

Melting glaciers release new ground surfaces, which may be either a source of greenhouse gas emissions or a sink for carbon dioxide. Studies carried out in subpolar and alpine ecosystems confirm the relatively rapid soil development and increase of carbon and nitrogen pools. However, observations from high-mountain glacier forelands in cold and dry climate are very scarce. This study analyses the impact of major environmental factors related to climate, topography, and vegetation, over a time-scale, on soil development and spatial soil differentiation in the foreland of Uisu Glacier, East Pamir Mountains. Moreover, the usefulness of the World Reference Base (WRB) and Soil Taxonomy in the classification of poorly developed soils in the ultracontinental climate was assessed. Geomorphological, pedological, and botanical surveys covered a sequence of terraces, alluvial fans, and end-moraines. Typical characteristics of the soils in the glacier foreland were: very high stoniness, coarse texture, high content of calcium carbonate, alkaline reaction, and low salinity. Soil development has extremely low intensity and was manifested in (a) soil organic carbon pools being among the lowest reported in the world (up to 1.4 kg m^-2 in the layer 0-50 cm), and (b) the presence of cambic/calcic horizons only on landforms older than of Mid-Holocene age (estimated). It was concluded that both the extremely cold and extremely dry climate conditions in the Uisu Glacier foreland limit the water flux and availability, suppress vegetation density and variability, and slow down the rate of soil development. Both WRB and Soil Taxonomy were able to reflect the advances in soil development and spatial soil differentiation (Calcaric Hyperskeletic Leptosols - Calcaric Cambisols - Cambic Calcisols, and Gelifluvents - Haplocambids - Haplocalcids, respectively); however, highlighting different features developed under an extremely cold and dry climate conditions of the East Pamir Mountains.

Topographic Wetness Index calculation guidelines based on measured soil moisture and plant species composition

Soil moisture controls environmental processes and species distributions, but it is difficult to measure and interpolate across space. Topographic Wetness Index (TWI) derived from digital elevation model is therefore often used as a proxy for soil moisture. However, different algorithms can be used to calculate TWI and this potentially affects TWI relationship with soil moisture and species assemblages. To disentangle insufficiently-known effects of different algorithms on TWI relation with soil moisture and plant species composition, we measured the root-zone soil moisture throughout a growing season and recorded vascular plants and bryophytes in 45 temperate forest plots. For each plot, we calculated 26 TWI variants from a LiDAR-based digital terrain model and related these TWI variants to the measured soil moisture and moisture-controlled species assemblages of vascular plants and bryophytes. A flow accumulation algorithm determined the ability of the TWI to predict soil moisture, while the flow width and slope algorithms had only a small effects. The TWI calculated with the most often used single-flow D8 algorithm explained less than half of the variation in soil moisture and species composition explained by the TWI calculated with the multiple-flow FD8 algorithm. Flow dispersion used in the FD8 algorithm strongly affected the TWI performance, and a flow dispersion close to 1.0 resulted in the TWI best related to the soil moisture and species assemblages. Using downslope gradient instead of the local slope gradient can strongly decrease TWI performance. Our results clearly showed that the method used to calculate TWI affects study conclusion. However, TWI calculation is often not specified and thus impossible to reproduce and compare among studies. We therefore provide guidelines for TWI calculation and recommend the FD8 flow algorithm with a flow dispersion close to 1.0, flow width equal to the raster cell size and local slope gradient for TWI calculation.