Topsoil organic matter build-up in glacier forelands around the world

Grassland degraded patchiness reduces microbial necromass content but increases contribution to soil organic carbon accumulation

Plant and microbially derived carbon (C) are the two major sources of soil organic carbon (SOC), and their ratio impacts SOC composition, accumulation, stability, and turnover. The contributions of and the key factors defining the plant and microbial C in SOC with grassland patches are not well known. Here, we aim to address this issue by analyzing lignin phenols, amino sugars, glomalin-related soil proteins (GRSP), enzyme activities, particulate organic carbon (POC), and mineral-associated organic carbon (MAOC). Shrubby patches showed increased SOC and POC due to higher plant inputs, thereby stimulating plant-derived C (e.g., lignin phenol) accumulation. While degraded and exposed patches exhibited higher microbially derived C due to reduced plant input. After grassland degradation, POC content decreased faster than MAOC, and plant biomarkers (lignin phenols) declined faster than microbial biomarkers (amino sugars). As grassland degradation intensified, microbial necromass C and GRSP (gelling agents) increased their contribution to SOC formation. Grassland degradation stimulated the stabilization of microbially derived C in the form of MAOC. Further analyses revealed that microorganisms have a C and P co-limitation, stimulating the recycling of necromass, resulting in the proportion of microbial necromass C in the SOC remaining essentially stable with grassland degradation. Overall, with the grassland degradation, the relative proportion of the plant component decreases while than of the microbial component increases and existed in the form of MAOC. This is attributed to the physical protection of SOC by GRSP cementation. Therefore, different sources of SOC should be considered in evaluating SOC responses to grassland degradation, which has important implications for predicting dynamics in SOC under climate change and anthropogenic factors.

Topsoil dilution by subsoil admixture had less impact on soil organic carbon stock development than fertilizer form and erosion state

Enhancing the agroecosystems carbon (C) sink function for climate mitigation faced challenges, particularly with traditional measures with limited suitability for increasing soil organic carbon (SOC) stocks. Inducing a SOC undersaturation in the topsoil by abrupt subsoil admixture is a way to create an additional C sink. However, the deep tillage traditionally used for this topsoil dilution was not always successful. It was due to a lack of knowledge and suitable approaches to record the effect of all relevant factors in SOC recovery, including soil conditions and fertilizer forms. We addressed these problems by establishing a three-factorial experiment: I) "moderate topsoil dilution," II) "N fertilization form," and III) "soil erosion state," representing three soil types in the hummocky ground moraine landscape of NE Germany. SOC dynamics were determined over a year of winter rye cropping using a novel robotic chamber system capable of measuring CO2 exchange on 36 experimental plots with a reduced methodological bias than previous measuring systems. The averaged net ecosystem carbon balance, a proxy for SOC stock change, indicated that topsoil dilution only reduced further SOC losses. The N fertilizer form had a significantly stronger and more differentiated effect. While the mineral N fertilization consistently produced only C sources, the organic fertilization, in combination with the diluted topsoil, led to a C sink. This C-sink function was, however, more pronounced in the eroded soil than in the non-eroded soil. Overall, the results have made clear that the impact of topsoil dilution on the further development of the SOC stock is only possible if the effect of other relevant factors, such as N fertilizer form and erosion state, are taken into account.

Chronological distribution and potential sources of persistent toxic substances in soils from the glacier foreland of Midtre Lovénbreen, Svalbard

Despite its reputation as one of the cleanest regions globally, recent studies have identified the presence of various persistent toxic substances (PTSs) in the environmental matrices collected from Svalbard. This study investigated the chronological distribution and potential sources of 81 PTSs in soils from the glacier foreland of Midtre Lovénbreen. Soil samples (n = 45) were categorized by age based on exposure to the atmosphere due to glacier retreat in July 2014 into five age groups: 80-100 years (n = 7), 60-80 years (n = 12), 40-60 years (n = 16), 20-40 years (n = 7), and <20 years (n = 3). Concentrations of polychlorinated biphenyls (PCBs, n = 32) in soils varied with age, ranging from 0.29 to 0.74 ng g-1 dw. In addition, the concentrations of polycyclic aromatic hydrocarbons (PAHs, n = 28), perylene, and alkyl-PAHs (n = 20) in soils ranged from 21 to 80 ng g-1 dw, 2.9-62 ng g-1 dw, and 73-420 ng g-1 dw, respectively. The concentrations of PTSs were observed to be greater in older soils. Principal component analysis revealed that PCBs in soils originated from various product sources. Positive matrix factorization modeling estimated the association of PAHs in soils with potential origins, such as diesel emissions, petroleum and coal combustion, and coal. Potential sources of PAHs were mainly coal in younger soils and diesel emissions and petroleum combustion in older soils. Alkyl-PAH compositions in the soil were similar to those of bituminous coal, with a noteworthy degree of weathering observed in older soils. The accumulation rate and flux of PTSs in soils exhibited compound-specific patterns, reflecting factors such as long-range transport, fate, origin, and recent inputs. These findings can serve as baseline data for protecting and preserving polar environments.

The development of terrestrial ecosystems emerging after glacier retreat

The global retreat of glaciers is dramatically altering mountain and high-latitude landscapes, with new ecosystems developing from apparently barren substrates1-4. The study of these emerging ecosystems is critical to understanding how climate change interacts with microhabitat and biotic communities and determines the future of ice-free terrains1,5. Here, using a comprehensive characterization of ecosystems (soil properties, microclimate, productivity and biodiversity by environmental DNA metabarcoding6) across 46 proglacial landscapes worldwide, we found that all the environmental properties change with time since glaciers retreated, and that temperature modulates the accumulation of soil nutrients. The richness of bacteria, fungi, plants and animals increases with time since deglaciation, but their temporal patterns differ. Microorganisms colonized most rapidly in the first decades after glacier retreat, whereas most macroorganisms took longer. Increased habitat suitability, growing complexity of biotic interactions and temporal colonization all contribute to the increase in biodiversity over time. These processes also modify community composition for all the groups of organisms. Plant communities show positive links with all other biodiversity components and have a key role in ecosystem development. These unifying patterns provide new insights into the early dynamics of deglaciated terrains and highlight the need for integrated surveillance of their multiple environmental properties5.

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.

Local climate modulates the development of soil nematode communities after glacier retreat

The worldwide retreat of glaciers is causing a faster than ever increase in ice-free areas that are leading to the emergence of new ecosystems. Understanding the dynamics of these environments is critical to predicting the consequences of climate change on mountains and at high latitudes. Climatic differences between regions of the world could modulate the emergence of biodiversity and functionality after glacier retreat, yet global tests of this hypothesis are lacking. Nematodes are the most abundant soil animals, with keystone roles in ecosystem functioning, but the lack of global-scale studies limits our understanding of how the taxonomic and functional diversity of nematodes changes during the colonization of proglacial landscapes. We used environmental DNA metabarcoding to characterize nematode communities of 48 glacier forelands from five continents. We assessed how different facets of biodiversity change with the age of deglaciated terrains and tested the hypothesis that colonization patterns are different across forelands with different climatic conditions. Nematodes colonized ice-free areas almost immediately. Both taxonomic and functional richness quickly increased over time, but the increase in nematode diversity was modulated by climate, so that colonization started earlier in forelands with mild summer temperatures. Colder forelands initially hosted poor communities, but the colonization rate then accelerated, eventually leveling biodiversity differences between climatic regimes in the long term. Immediately after glacier retreat, communities were dominated by colonizer taxa with short generation time and r-ecological strategy but community composition shifted through time, with increased frequency of more persister taxa with K-ecological strategy. These changes mostly occurred through the addition of new traits instead of their replacement during succession. The effects of local climate on nematode colonization led to heterogeneous but predictable patterns around the world that likely affect soil communities and overall ecosystem development.

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.

Microbial dynamics in soils of the Damma glacier forefield show succession in the functional genetic potential

Glacier retreat is a visible consequence of climate change worldwide. Although taxonomic change of the soil microbiomes in glacier forefields have been widely documented, how microbial genetic potential changes along succession is little known. Here, we used shotgun metagenomics to analyse whether the soil microbial genetic potential differed between four stages of soil development (SSD) sampled along three transects in the Damma glacier forefield (Switzerland). The SSDs were characterized by an increasing vegetation cover, from barren soil, to biological soil crust, to sparsely vegetated soil and finally to vegetated soil. Results suggested that SSD significantly influenced microbial genetic potential, with the lowest functional diversity surprisingly occurring in the vegetated soils. Overall, carbohydrate metabolism and secondary metabolite biosynthesis genes overrepresented in vegetated soils, which could be partly attributed to plant-soil feedbacks. For C degradation, glycoside hydrolase genes enriched in vegetated soils, while auxiliary activity and carbohydrate esterases genes overrepresented in barren soils, suggested high labile C degradation potential in vegetated, and high recalcitrant C degradation potential in barren soils. For N-cycling, organic N degradation and synthesis genes dominated along succession, and gene families involved in nitrification were overrepresented in barren soils. Our study provides new insights into how the microbial genetic potential changes during soil formation along the Damma glacier forefield.

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.

Llamas (Llama glama) enhance proglacial ecosystem development in Cordillera Blanca, Peru

Worldwide, mountain glaciers are shrinking rapidly. Consequently, large areas are becoming available for the development of novel alpine ecosystems. These harsh environments, however, delay primary succession. In this study with a local community, we conducted an inclusion experiment to investigate whether Llama glama influences soils and vegetation primary succession following glacial retreat. At the foot of the Uruashraju glacier in the Cordillera Blanca, Peru (~ 4680 m.a.s.l.), we established four llama inclusion plots and four control plots that we studied from 2019 to 2022, 24-40 years after deglacierization. After three years, the llama plots had significantly increased soil organic carbon and soil nitrogen. In the llama plots, we found a large, significant increase in vascular plant cover (+ 57%) between the second and third years of experimentation, and we identified four new species that were not present in 2019. Our results suggest that Llama glama, through their latrine behavior and role as a seed disperser, enhances the primary succession and novel ecosystem formation in recently deglacierized landscapes. Our study provides scientific support that rewilding of native Andean camelids may favor adaptation to glacier retreat and inform conservation and management strategies in proglacial landscapes.

Heterogeneous changes of soil microclimate in high mountains and glacier forelands

Landscapes nearby glaciers are disproportionally affected by climate change, but we lack detailed information on microclimate variations that can modulate the impacts of global warming on proglacial ecosystems and their biodiversity. Here, we use near-subsurface soil temperatures in 175 stations from polar, equatorial and alpine glacier forelands to generate high-resolution temperature reconstructions, assess spatial variability in microclimate change from 2001 to 2020, and estimate whether microclimate heterogeneity might buffer the severity of warming trends. Temporal changes in microclimate are tightly linked to broad-scale conditions, but the rate of local warming shows great spatial heterogeneity, with faster warming nearby glaciers and during the warm season, and an extension of the snow-free season. Still, most of the fine-scale spatial variability of microclimate is one-to-ten times larger than the temporal change experienced during the past 20 years, indicating the potential for microclimate to buffer climate change, possibly allowing organisms to withstand, at least temporarily, the effects of warming.

Future emergence of new ecosystems caused by glacial retreat

Glacier shrinkage and the development of post-glacial ecosystems related to anthropogenic climate change are some of the fastest ongoing ecosystem shifts, with marked ecological and societal cascading consequences1-6. Yet, no complete spatial analysis exists, to our knowledge, to quantify or anticipate this important changeover7,8. Here we show that by 2100, the decline of all glaciers outside the Antarctic and Greenland ice sheets may produce new terrestrial, marine and freshwater ecosystems over an area ranging from the size of Nepal (149,000 ± 55,000 km2) to that of Finland (339,000  ±  99,000  km2). Our analysis shows that the loss of glacier area will range from 22 ± 8% to 51 ± 15%, depending on the climate scenario. In deglaciated areas, the emerging ecosystems will be characterized by extreme to mild ecological conditions, offering refuge for cold-adapted species or favouring primary productivity and generalist species. Exploring the future of glacierized areas highlights the importance of glaciers and emerging post-glacial ecosystems in the face of climate change, biodiversity loss and freshwater scarcity. We find that less than half of glacial areas are located in protected areas. Echoing the recent United Nations resolution declaring 2025 as the International Year of Glaciers' Preservation9 and the Global Biodiversity Framework10, we emphasize the need to urgently and simultaneously enhance climate-change mitigation and the in situ protection of these ecosystems to secure their existence, functioning and values.

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.

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.

Soil Formation and Mass Redistribution during the Holocene Using Meteoric 10Be, Soil Chemistry and Mineralogy

Soil development and erosion are important and opposing processes in the evolution of high-mountainous landscapes, though their dynamics are not fully understood. We compared soil development between a calcareous and a siliceous chronosequence in the central Swiss Alps at high altitudes, which both cover soil formation over the Holocene. We calculated element mass balances, long-term erosion rates based on meteoric 10Be and we determined the rates of soil formation. We also analyzed the shifts in the mineralogical composition, weathering indices, the particle size distribution, carbon stocks and oxalate extractable Fe, Al, and Mn. The siliceous soils had high chemical weathering rates at the early stage of soil formation that strongly decreased after a few millennia. The development of calcareous soil was characterized by high carbonate losses and a shift to finer soil texture. Soil erosion hampered the upbuilding of soil horizons in the early stages of soil development, which led to a delay in soil and vegetation development. This study shows how soil formation drivers change over time. In the early stages of soil development, the parent material predominantly drives soil formation while at later stages the vegetation becomes more dominant as it influences surface stability, hydrological pathways, and chemical weathering that determine water drainage and retention.

Dynamics of Ecological Communities Following Current Retreat of Glaciers

Glaciers are retreating globally, and the resulting ice-free areas provide an experimental system for understanding species colonization patterns, community formation, and dynamics. The last several years have seen crucial advances in our understanding of biotic colonization after glacier retreats, resulting from the integration of methodological innovations and ecological theories. Recent empirical studies have demonstrated how multiple factors can speed up or slow down the velocity of colonization and have helped scientists develop theoretical models that describe spatiotemporalchanges in community structure. There is a growing awareness of how different processes (e.g., time since glacier retreat, onset or interruption of surface processes, abiotic factors, dispersal, biotic interactions) interact to shape community formation and, ultimately, their functional structure through succession. Here, we examine how these studies address key theoretical questions about community dynamics and show how classical approaches are increasingly being combined with environmental DNA metabarcoding and functional trait analysis to document the formation of multitrophic communities, revolutionizing our understanding of the biotic processes that occur following glacier retreat.

The Retreat of Mountain Glaciers since the Little Ice Age: A Spatially Explicit Database

Most of the world’s mountain glaciers have been retreating for more than a century in response to climate change. Glacier retreat is evident on all continents, and the rate of retreat has accelerated during recent decades. Accurate, spatially explicit information on the position of glacier margins over time is useful for analyzing patterns of glacier retreat and measuring reductions in glacier surface area. This information is also essential for evaluating how mountain ecosystems are evolving due to climate warming and the attendant glacier retreat. Here, we present a non-comprehensive spatially explicit dataset showing multiple positions of glacier fronts since the Little Ice Age (LIA) maxima, including many data from the pre-satellite era. The dataset is based on multiple historical archival records including topographical maps; repeated photographs, paintings, and aerial or satellite images with a supplement of geochronology; and own field data. We provide ESRI shapefiles showing 728 past positions of 94 glacier fronts from all continents, except Antarctica, covering the period between the Little Ice Age maxima and the present. On average, the time series span the past 190 years. From 2 to 46 past positions per glacier are depicted (on average: 7.8).