Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Plant Sci., 10 February 2026

Sec. Functional Plant Ecology

Volume 17 - 2026 | https://doi.org/10.3389/fpls.2026.1670208

Climate and biocrust types jointly regulate soil multifunctionality and quality in drylands: evidence from the Gurbantunggut Desert

Yonggang Li*Yonggang Li1*Yingjie GaoYingjie Gao1Yunjie HuangYunjie Huang2Yongxing LuYongxing Lu2Benfeng YinBenfeng Yin2Xiaobing ZhouXiaobing Zhou2Hao YuHao Yu1Yuanming Zhang*Yuanming Zhang2*
  • 1School of Plant Protection and Environment, Henan Institute of Science and Technology, Xinxiang, China
  • 2Key Laboratory of Ecological Safety and Sustainable Development in Arid Lands, Xinjiang Institute of Ecology and Geography, Chinese Academy of Sciences, Urumqi, China

Soil multifunctionality (SMF) and the soil quality index (SQI) are essential indicators of soil function, productivity, and health. Additionally, the spatial variability of soil multifunctionality (SVM) signifies soil heterogeneity. Biological soil crusts (Biocrusts) can affect these indicators. However, there is little information about the role of biocrusts in regulating the response of multiple ecosystem functions to climate change. We evaluated the relative importance of climate, soil environment, and biocrusts variables as drivers of SMF, SQI, and SVM at 74 sites in the Gurbantunggut Desert. Soil SMF, and SQI increase with the coverage of lichen and moss crust. Biocrusts index, SMF and SQI increase with an increase in the mean annual temperature. Biocrusts index, SMF and SQI increase first with an increase in mean annual precipitation (MAP)< 163 mm and then decrease. SVM display a significant decreasing trend with the increase of MAP. The structural equation model (SEM) demonstrate that the spatial distribution can significantly influence the biocrusts, soil SQI and SVM. Biocrusts has a significant positive influence on soil SMF (0.47)and SQI (0.31). Soil SMF has a significant negative effect on SVM (-0.50), and SQI (0.59) has a significant positive effect. We provide the first quantitative evidence that biocrust type and a 163 mm precipitation threshold govern SMF through opposing direct vs. indirect temperature pathways, offering a predictive rule-of-thumb for dryland management under climate change. The findings contribute decidedly to our understanding of the patterns and mechanisms driving SMF, SQI, and SVM in drylands, which is important for predicting changes in ecosystem function under climate change.

1 Introduction

When considering terrestrial ecosystems, most ecosystem functions are provided by soils. Soils (abiotic and biotic components) directly contribute to several ecosystem functions such as nutrient and water cycling and biodiversity. Therefore, soil carries immense biodiversity and support key ecosystem processes essential for life (Eldridge et al., 2020; Guerra et al., 2020). These processes encompass nutrient cycling (Hu et al., 2023; Ochoa-Hueso et al., 2018), trophic interactions (Hagh-Doust et al., 2023), individual plant performance, and competitive ability (Day et al., 2003; Huang et al., 2023), as well as community-level productivity (Garcia-Callejas et al., 2023). Soils can be characterized by different physical, chemical and biological properties (Bouma, 2002). These soil properties influence soil multifunctionality (SMF), soil quality index (SQI), and, subsequently, ecosystem processes.

SQI is a comprehensive metric used to evaluate the health status of soil and its capacity to support specific functions (e.g., agricultural production, environmental protection). It quantifies soil quality by integrating multiple physical, chemical, and biological indicators of soil. SQI can assess the capacity of the soil to maintain and support its productivity and health in the relevant ecosystem by combining several physical, chemical, and biological properties (Lal et al., 2021). The SQI assessment would have important implications for resource management and provides a quantitative method for evaluating the soil quality at the soil type-scale (Kafei et al., 2023; Khosravi et al., 2023). Changes in SQI are linked to interactions among ecosystem processes, potentially disrupting the equilibrium between soil physiochemical, microbial, and biochemical processes (Hedo de Santiago et al., 2016). Studies indicate that external materials, like fertilizers, can modify microbial metabolic traits by changing soil nutrient levels, subsequently influencing soil SMF (Huo et al., 2024; Yang et al., 2023). Therefore, the SMF and SQI are calculated from soil properties that respond quickly to grazing, cultivation, or restoration, any change in land use is immediately reflected in the index values (Zornoza et al., 2015). Using multivariate indicators is the best way to estimate SQI because they provide a broader view of the situation than individual property analyses.

In terrestrial ecosystems, there exists a widespread pattern of uneven soil characteristic distribution, which plays a pivotal role in regulating diverse ecosystem processes. Notably, this pattern is particularly pronounced in dryland ecosystems, will directly affect the spatial heterogeneity of desert soil. The spatial variation in soil multifunctionality (SVM) plays a key role in the sustainability and stability of ecosystem functions (Juhos et al., 2016). SVM refers to the non-uniform distribution of SMF (i.e., the ability of soil to simultaneously support multiple ecosystem functions) across space. This variation can occur at different scales. The SVM in drylands is likely to be caused by the significant multifunctional differences between vegetated patches, in which plants largely drive biological processes such as litter decomposition, nutrient cycling. In unvegetated areas physical processes such as wind and water play a larger role in SMF than biotic processes (Durán et al., 2018; Li et al., 2007; Zheng et al., 2019). The SVM is influenced by the complex interplay of biotic and abiotic processes (Allington and Valone, 2014; Ochoa-Hueso et al., 2018). Furthermore, climate changes can also affect the variability in soil functions among different vegetation patch types by altering patch patterns (Durán et al., 2018). Due to the sensitivity of drylands to climate change, increasing aridity and temperature would lead to change in vegetation and soil properties that could adversely affect SVM in these areas worldwide. Therefore, climate change and plant cover are associated with alterations of the SVM in drylands (Ding and Eldridge, 2021; Durán et al., 2018). According to recent studies, increasing aridity is correlated with decreasing plant cover and richness, and increasing woody vegetation encroachment rates across the globe (Delgado-Baquerizo et al., 2013; Vicente-Serrano et al., 2012), which would likely lead to an increase in SVM. Evaluating the effect of environmental factors on desert soil SQI and SVM is crucial not only for understanding the processes of desert ecosystems but also for developing practical and sensitive multivariate indices that can be used as management tools in desert regions.

Drylands are crucial for global sustainability, as they constitute 45% of the earth’s land surface (Reynolds et al., 2007), where many vascular plants are restricted due to the shortage of precipitation; however, biocrusts are widespread. Recent estimates indicate that biocrusts currently cover approximately 30% of dryland soils, constituting around 12% of the earth’s terrestrial surface (Rodriguez-Caballero et al., 2018). Biocrusts are composed of cyanobacteria, green algae, lichens, mosses, and other organisms related to soil particles that play essential fundamental roles in arid and semiarid regions (Xiao et al., 2022), including carbon and nitrogen cycling (Bowker et al., 2013; Li et al., 2019), surface energy balance (Couradeau et al., 2016; Rodriguez-Caballero et al., 2018; Rutherford et al., 2017), erosion (Cantón et al., 2014; Chamizo et al., 2017) and water redistribution (Bowker et al., 2013; Kidron and Büdel, 2014), affecting the colonization and development of vascular plants (Langhans et al., 2009; Li et al., 2005), and supplying habitats for other microorganisms and protozoa (Liu et al., 2017).

The biocrusts types can significantly influence the characteristics of soil nutrient cycling (Gao et al., 2018). There is a distinct difference between soil physicochemical properties, microbial activities, and community compositions of different types of biocrusts (Niu et al., 2017; Xu et al., 2021, Xu et al., 2020; Zhang et al., 2018, Zhang et al., 2016). In different types of biocrusts, soil fungal and bacterial communities have different responses and undertake important roles in maintaining nutrient cycling and SMF (Xu et al., 2021, Xu et al., 2020). There are significant differences in the physical and chemical characteristics of the underlying soil under different types of biocrusts, which will influence soil microbial diversity and activity (Lan et al., 2012; Li et al., 2020). Biocrusts notably enhance the content of C, N, and P in surface soil, and change soil nutrient stoichiometry (Baumann et al., 2021). Cyanobacteria vary significantly in diversity, biomass, and species composition due to resource accumulation among different biocrusts types (Lan et al., 2012; Wu et al., 2009). The moss crusts show higher photosynthesis and N2 fixation activities than the cyanobacterial and lichen crusts (Housman et al., 2006; Lan et al., 2012). Similarly, nitrogenase activity is higher in cyanobacterial and lichen crusts than in moss crusts (Pushkareva et al., 2017; Wu et al., 2009). Biocrust development exerted direct effects on SMF in arid regions, but only indirect effects through changing soil microbial biomass carbon in semi-arid regions (Su et al., 2021). Moss crusts significantly influence soil SMF. For instance, moss crusts alter phosphorus functionality and enhance phosphorus cycling in the ecosystem (Li et al., 2024). Therefore, biocrusts actively participate in soil surface heterogeneity dynamics in terms of biological diversity, soil function and physicochemical properties associated with their spatial structure (Ettema and Wardle, 2002).

Biocrusts are generally located within the uppermost millimeters of the soil surface (Weber et al., 2022) where they are more likely to be affected by environmental factors, such as soil water, nutrients, temperature, and radiation intensity (Grote et al., 2010; Navas Romero et al., 2020). As a result, the soil nutrients in the soil under biocrusts patches are also affected by environmental changes and are considerably higher than those in soil without biocrusts (Li et al., 2019; Tao et al., 2020). Some studies have explored the effects of biocrusts type on soil SMF (Li et al., 2019, Li et al., 2024; Su et al., 2021). However, it is not clear whether the SVM and SQI in desert soil are affected by climate and biocrusts type in the temperate desert. Additionally, does climate affect the SVM, SQI and SMF of desert soils by influencing biocrusts type. To evaluate the role and relative importance of climate and biocrusts type on the soil SMF, SQI, and SVM in desert ecosystems, we collected 74 sites from the Gurbantunggut Desert. Furthermore, previous studies also found that the climate can significantly influence the growth and distribution of different biocrusts types, as well as the physical and chemical properties and functions of soil (Li et al., 2023, Li et al., 2019). The biomass of moss crusts is 3 to 5 times greater than that of algae-lichen crusts (Zhang et al., 2022). Additionally, their thickness and nifH nitrogenase activity are 2 times and 2.3 times higher, respectively. The pseudoroots and polysaccharide carbon input from moss crusts significantly enhance surface SOC and available nitrogen, directly increasing the importance of “nutrient cycling” in SMF. Thus, we hypothesized that 1) Moss crusts increase SMF and SQI more than algal-lichen crusts, with the largest gains in nutrient-cycling functions; 2) Climate can directly influence soil SVM, SMF, and SQI, and can also indirectly affect SMF, SVM, and SQI by altering the type of biocrusts.

2 Materials and methods

2.1 Study area

This study was conducted in the Gurbantunggut Desert (44°11′-46°20′ N, 84°31′-90°00′ E, 300–600 m a.s.l.), which is situated in the center of the Jungger Basin, Central Asia. With a total area of 4.88 × 104 km2, it is the largest fixed and semi-fixed desert in China. Moist air currents from the Indian Ocean are blocked by the Himalayas and fail to reach this area, causing a vast expanse of arid terrain. The mean annual precipitation (MAP) in the desert is around 103–229 mm (Supplementary Figure S2) with gradient distribution (Zhang et al., 2007). The mean annual temperature (MAT) is 6-8°C, while the potential mean annual evaporation is estimated at 2606.6 mm (Zhang et al., 2007). The mean annual wind speed (MAW) of 2–4 m/s. With a mean wind speed of 11.17 m/s, late spring is the windiest time of year. There is prevailing wind from the northwest, northwest and north direction. Soil moisture is about 0.5-2% and the temperature is about 50-60°C. According to nature-reserve management records (1994–2022), zero-grazing policy, and remote-sensing evidence of negligible livestock traces, the study area has experienced no significant anthropogenic disturbance for at least the past three decades. The soil type is an Arenic soil and subclasses of Medium sand (Li et al., 2024).

The distribution of different types of biocrusts was found to be uneven in this desert, with no observed seasonal differences. On the dunes, the species composition varied in different biocrusts types, which could be classified into three groups: algal, lichen, and moss crusts (Supplementary Figure S2). Algal crusts are weakly consolidated, soft, and readily decomposed; they often include fungal or cyanobacterial hyphae and are composed mainly of fungi and/or cyanobacteria. Lichens consist of a symbiotic association between heterotrophic fungi and autotrophic partners such as cyanobacteria or green algae. Moss crusts are characterized by short, hairlike extensions and appear brown when dry and green to brown−green when moist. The type and coverage of biocrusts also vary in different regions with climate change. In this study, algal crust and lichen crust generally coexist in biocrusts, so we consider them as one soil type. Previous studies have shown that algae and lichen crusts have no significant differences in key functional indicators such as soil carbon and nitrogen sequestration, surface aggregate stability, and enzyme activity, and their functional contributions are significantly lower than moss crusts; Recently, large-scale biological crust research in arid regions has generally adopted the “algae lichen crust” merging treatment to avoid classification uncertainty caused by visual discrimination errors. Therefore, we selected two biocrusts types: algal, lichen mixed crust and moss crust in our study (Supplementary Figures S2S4).

2.2 Sample collection and processing

Sites were established at 10 km intervals from the southeast to the northwest of the Gurbantunggut Desert. Each site was divided into five 10 m × 10 m plots (Figure 1). In GIS, divide the east-west and north-south main roads of the Gurbantunggut Desert into grids, randomly select one grid as the starting point, and avoid obvious human interference areas (oil wells, 500 meter buffer zones on highways). Starting from a random starting point, set up one main sampling point every 10 km, for a total of 74 points. Within a range of 100 m × 100 m for each main sampling point, five 10 m × 10 m subplots are set using the system grid method (with 10 m intervals) to ensure spatial representativeness. Each plot was established in a typical interdune area which is flat with no slope. In this study, five 1 m × 10 m transect was established in a subplot. The transect was divided into 10-square grid plots with the aim of biocrusts type, coverage, distribution and patch size more conveniently and precisely (Figure 1). Within each 1 m × 1 m grid, first use a 30 cm × 30 cm sub grid to determine the patches with a biological crust coverage of ≥ 30%. All subsequent soil samples are strictly located directly below the patch and do not deviate laterally. Soil samples were selected from expose area with a minimum distance of 50 cm from surrounding shrubs. Soil samples were randomly collected within subplots. One soil sample in each subplot, which is a mixture of five randomly selected samples within the plot. If there is no biocrusts distribution, bare sand was collected. To precisely compare the properties of desert surface soils, we uniformly collected soil samples from 0-5cm of desert surface soil using a 5 cm height and diameter cutting ring. In total, 370 soil samples were collected from 74 sites, including 300 soils with biocrusts (80 alga-lichen and 220 moss) and 70 bare sand samples (Figure 1). For soil of biocrusts patches, we also collected the top 0–5 cm of soil. However, after collection, we removed the plant components from the biocrust’s soil samples, retaining only the soil part for analysis. The samples were collected in August 2018. Firstly, soil sand content was measured by weighing 50 g of dry soil, crushing, and passing through sieves with different mesh diameters to obtain soil particles of different sizes. The different soil particles were weighted and calculated the weight ratio. Soil samples were stored in a cool, dry place at an average room temperature of 25°C for 2 weeks, until the soil was dry.

Figure 1
Map panel A displays a specific study area in China, marked with red, green, and blue dots indicating moss, lichens and algae, and sand. Panel B details a subplot arrangement with five black squares. Panel C shows transects with alternating black and white vertical stripes. Panel D illustrates a quadrant with a horizontal pattern of alternating black and white squares.

Figure 1. Site distributed of sample collection in Gurbantünggüt Desert, Northwest China. (A) displays a specific study area in China, marked with red, green, and blue dots indicating moss, lichens and algae, and sand. (B) details a subplot arrangement with five black squares. (C) shows transects with alternating black and white vertical stripes. (D) illustrates a quadrant with a horizontal pattern of alternating black and white squares.

2.3 Soil analysis

Use ASTM E11 standard copper sieve, oscillating sieve (Retsch AS200) with an amplitude of 1.5 mm and a screening time of 10 minutes. The sieve apertures (mm) are 2.0, 1.0, 0.5, 0.25, 0.15, and 0.075, respectively. Weigh them step by step to obtain the distribution of gravel and sand particles with an accuracy of 0.01 g, which are larger than 2 mm. A calibrated pH meter (PHS-4, Jiangsu Manufactory of Electrical Analysis Instruments, Jiangyin, China) was used to measure soil pH in a 1:5 soil:water suspension. DDS-11A (LEIC, Shanghai, China) was used to measure electrical conductivity (EC). Soil nutrient levels, i.e., organic carbon (SOC), total nitrogen (TN), total phosphorus (TP), total potassium (TK), NO3 + NO2-N (NO3-N), NH4-N, available phosphorus (AP), and available potassium (AK) were determined (Li et al., 2019). The dichromate oxidation method was performed to determine SOC contents. Soil organic carbon uses Chinese national standard substance GBW-07414 (standard value 1.05 ± 0.04%), with a recovery rate of 97-103%. soil TN was estimated by the Kjeldahl procedure after digestion with concentrated H2SO4 on a distillation unit, and TP was determined by the H2SO4 ammonium molybdate-ascorbic acid method (Barrett et al., 2007). TN and TP use GBW-07415, with a recovery rate of 96-104%. TK was quantified by using inductively coupled plasma spectrometry (Perkin Elmer Optima 3000-DV ICP, Perkin Elmer Inc., Shelton, Connecticut, USA). In order to obtain extractable inorganic nitrogen (NO3 + NO2-N and NH4-N), 10 g soil was extracted in 50 mL 2 M KCl, filtered, and frozen until run on a Lachat auto-analyzer (Barrett et al., 2007). For the analysis of the AP content, we used the molybdenum-antimony colorimetric method (Li et al., 2019). To measure the AK, 5 g soil was extracted in 50 mL 1 mol/L NH4Ac, filtered, and then frozen until run on an atomic absorption spectrophotometer (Perkin Elmer model 2380, Perkin Elmer Inc., USA).

Urease activity (µmol PNP g⁻¹ h⁻¹) was measured by the buffered hydrolysis reaction method modified by Li et al. (2019). In alkaline environments, soil alkaline phosphatase catalyzes the production of yellow p-nitrophenol from p-nitrophenyl phosphate, which has a maximum absorption peak at 405 nm. By detecting the rate of increase of PNP at 405 nm, the activity (µmol PNP g⁻¹ h⁻¹) of the alkaline phosphatase enzyme can be obtained. soil β-glucosidase can catalyze p-nitrobenzene-β-D-glucopyranoside generates a yellow substance, p-nitrophenol (PNP), which has characteristic light absorption at 405 nm, thereby obtaining the activity (µmol PNP g⁻¹ h⁻¹) of soil β-glucosidase (Li et al., 2019). Three technical repetitions were performed for each sample, RSD< 5%. β - glucosidase: LOD = 0.004 µmol PNP g-1 h-1, LOQ = 0.012 µmol PNP g-1 h-1; Alkaline phosphatase: LOD = 0.005 µmol PNP g-1 h-1, LOQ = 0.015 µmol PNP g-1 h-1; All measured values are>LOQ (Supplementary Table S2).

2.4 Soil multifunctionality

We estimated soil SMF using crucial functional variables, including soil SOC, TN, TP, TK, NH4-N, NO3-N, AP, AK, soil carbon, nitrogen, and phosphorus cycling enzymes. The selected functions are closely associated with soil fertility and soil nutrient cycling. They are commonly used as indicators to estimate soil nutrient cycling and soil fertility in drylands. Individual functions were standardized using the min-max transformation with a range from 0 to 1 (Equation 1). The minimum and maximum values of a soil function were estimated for the 74 studied sites.

f(Xi)=(Ximin(X))/(max(X)min(X))(1)

where Xi represents the value observed for the functional variables, min(X) denotes the minimum value, and max(X) denotes the maximum value observed. The index’s range is transformed to 0–1 through data normalization. Then, the soil SMF was calculated using the “averaging approach” (Equation 2) which is widely used in multifunctionality studies (Bowker et al., 2013; Su et al., 2021). The formula for calculating the P-SMF is presented as follows:

SMFi=1Ni= 1Nfi(2)

In this equation, SMFi indicates the multifunctionality for sample i. N represents the overall functional variables count integrated into the computation, while fi denotes the standardized index of the functional variables for sample i. The arithmetic mean SMF was double checked by PCA weighting and threshold method, and the effect direction and significance were both robust.

2.5 Soil quality index

A scoring function analysis framework was used to calculate the SQI (Andrews et al., 2004; Karlen et al., 2008). The SQI was assessed by following a three-step procedure: (1) identification of the minimum dataset of indicators, (2) indicator interpretation, and (3) integration of all indicator scores into one overall SQI value (Andrews et al., 2004; Wang et al., 2024). Principal component analysis was performed on the standardized data matrix of the total dataset to determine potential soil indicators representing the minimum data set. Select principal components with eigenvalue > 1 and principal components with a total variation of at least 5% in the interpretation data set for minimum dataset identification. For each principal components, only the high load index whose load value is within 10% of the maximum weighted load is retained as an important index for indexing principal components.

Kaiser-Meyer-Olkin and Bartlett’s tests were used to determine whether the data were appropriate for principal component analysis. The Kaiser-Meyer-Olkin test was used to assess the correlation between the input variables. If the value of the Bartlett’s tests is significant (p< 0.05), it shows the correlation between the variables and the appropriateness of principal component analysis. In this study, the Kaiser-Meyer-Olkin statistic was > 0.6, and the BTS test was significant. In the minimum dataset, each highly loaded indicator was retained if they were not correlated. Otherwise, only the indicator with the highest weighted loading was chosen for the minimum dataset. After defining the minimum dataset for the SQI, each soil indicator was transformed into unitless combinable scores varying from 0.00 to 1.00 using linear and nonlinear scoring function methods (Andrews et al., 2004, Andrews et al., 2002; Raiesi, 2017). The appropriate scoring algorithms for “scoring indicators” values were selected and interpreted for soil productivity and sustainability. “More is better” and “less is better” scoring curves were used to indicators when a soil indicator was considered good for SQI in increasing order (more is better) or in decreasing order (less is better). For linear scoring, “more is better” (Equation 3) or “less is better” (Equation 4) functions were used as follows:

Si=0.1+((x-b)(a-b))×0.9(3)
Si=1-((x-b)(a-b))×0.9(4)

Where Si was the transformed x variable, b and a are the minimum and maximum threshold values of the x (non-transformed or true variable), respectively.

Nonlinear scoring functions (NLSF), proposed by Masto et al. (2008), were used to normalize SQI (Equation 5)as below:

Si=1[1+eβ(xα)](5)

Where a is the baseline value of the soil variable, where the score is 0.5 or close to the mean value of the upper and lower thresholds and b is the slope. Baselines are generally considered the lowest target values. An acceptable condition can be considered for the system if the soil indicator value is within the threshold (control) limits. If the value is not within the control limits, it is considered a degraded system (Masto et al., 2008).

After scoring the selected indicators, we applied a weighted additive approach to integrate them into the indices (Andrews et al., 2002). In this study, the weights (Wi) of the minimum dataset indicators were given by the ratio of the variance of each variable to the total cumulative variance to get a certain weight value under principal component analysis (Tesfahunegn, 2014). The scoring was calculated using the “more is better” method in this study. After weighing the scored minimum dataset indicators, the integrated SQI was calculated using the below equation (Qi et al., 2009):

SQI=i=1nWiSi(6)

where Wi and Si are the weight and score of indicators i, respectively.

Perform PCA (KMO = 0.74, Bartlett p<0.001) after standardizing the 12 original indicators. Retain the first four principal components based on eigenvalues>1, and explain 83.1% of the variance cumulatively. If Pearson | r |>0.70 for both indicators, only retain the one with the highest load (SOC and TN r=0.71, retain SOC); The final minimum dataset (MDS) contains 7 indicators: SOC, TP, AK, ALP, BG, pH, and EC. The final SQI (0-1) is obtained by weighting and summing the variance contribution rate of each indicator’s PC to the cumulative variance ratio. After changing from non-linear to linear scoring, the SQI (Equation 6) and original results have r=0.93, and the direction of the effect remains unchanged. Delete any indicator, and if the site ranking Kendall τ>0.89, it indicates that the weight and function selection are robust.

2.6 Spatial variation of soil multifunctional

SVM is estimated by first calculating the variability of the soil variables measured at each site. To do so, we estimated the site-level coefficient of variation (CV) using the composite soil samples collected at each site. The CV is a relative measure of heterogeneity that can accommodate variance-mean scaling, preventing variances from increasing with the mean. Therefore, it is more useful for comparing variability within biological properties than absolute measures of variability, such as standard deviation (Fraterrigo and Rusak, 2010; Schlesinger et al., 1990). We estimated SVM (Equation 7) as the arithmetic mean for all individual site-level CVs of all soil variables.

SVMi=1Nj= 1NCVj(7)

where CVj represents the value of CV for the functional variables in a site. Then, the soil SVM was calculated using the “averaging approach” which is widely used. SVMi indicates the SVM for site i. N represents the overall functional variables count integrated into the computation.

The latitude, longitude, and altitude data for the sites were acquired using a GPS locator. The environmental site information was sourced from the National Meteorological Science Data Center (http://data.cma.cn/site/index.html) and the Resource and Environmental Science Data Center of the Chinese Academy of Sciences (http://www.resdc.cn) and included variables such as the MAP, MAW, radiation (Rad) and MAT. The data source has a spatial resolution of 0.25° × 0.25° (approximately 25 km) and a 30-year climate standard value spanning from 1981 to 2010. Using ArcGIS 10.8, perform bilinear interpolation between the GPS coordinates (WGS84) of 74 sampling points and the nearest grid point to obtain the station level annual average; Interpolation error ≤ 0.3°C (MAT) and ≤ 3 mm (MAP), meeting the requirements of regional scale research.

2.7 Statistical analyses

We used R 4.1 software to perform all statistical analyses. We used linear mixed-effects models to account for the hierarchical structure (Plot nested within Site) when testing the effect of biocrust type on SMF, SQI and SVM. Here, we calculated the synthesized biocrusts index following the same process as SMF using biocrusts type, thickness, cover and patch characteristics, which are variables that have been used to indicate biocrust type development and characteristics in previous studies (Belnap et al., 2008; Chamizo et al., 2012; Kidron et al., 2010; Lan et al., 2013). The biocrusts index is closely correlated with visual biocrusts type development (Su et al., 2021). We explored the relationships between MAT, MAP, MAW using regression analysis to evaluate the effects of abiotic factors on biocrusts index. We explored the relationships between biocrusts index, MAT, MAP, MAW, soil pH and EC using regression analysis to evaluate the effects of biotic and abiotic factors on soil SMF, SQI and SVM. Only when non-linear regressions were a better fit to the data, thresholds may be present. Collinearity was absent (VIF< 3.2); MAP × MAT interaction and GAMM non-linearity confirmed the 163 mm threshold without significant interaction terms. Therefore, we explored the presence of thresholds only when non-linear models were a better fit to the data. We did so because segmented regressions threshold models force the existence of at least one threshold. The regression analysis of MAP and SMF, MAP and SQI used segmented linear regression to find the threshold and then conducted linear regression analysis. The breakpoint was estimated by segmented regression with 9999 bootstrap replicates, yielding 163 mm (95% CI = 154 - 172).The regression analysis of soil pH with SMF and SQI used segmented linear regression to find the threshold and then performed linear regression analysis.

In this study, we employed structural equation modeling (SEM) to investigate the direct and indirect effects of biotic and abiotic drivers on the SMF, SQI and SVM. Firstly, we establish an initial model. Subsequently, the model was revised based on the initial model to achieve the optimal model. The following criteria were used to evaluate the model fit: Fisher’s p value (model fits well at 0.05< p< 1.00) and Akaike’s information criterion (AIC; model fits well at low AIC). Spatial Moran’s I and d-separation tests showed no residual autocorrelation (p > 0.35), supporting the independence assumption of SEM. Other unnecessary pathways were eliminated, and only meaningful pathways were retained in the final model. To construct the SEM and gain a comprehensive understanding of the primary drivers of the SVM, we utilized the “piecewiseSEM” package.

3 Results

3.1 SMF, SQI and SVM of biocrusts

Soil under moss crusts exhibited the highest soil SMF and SQI compared with bare sand and algae-lichen crusts (Figure 2). The SQI of soil beneath algae-lichen crusts was significantly higher (p < 0.05) than that of sand. There was no significant difference between sand and algae-lichen crusts. Surface soil TN, available nitrogen (NO3-N and NH4-N), AP, AK, pH, and EC differed significantly among biocrust types (p < 0.05) (Supplementary Figures S5S7). The biocrust index increased with rising MAT (p = 0.003, R2 = 0.12) and MAW (p < 0.001, R2 = 0.29) (Figure 3). The biocrust index rose with MAP up to 163 mm (p = 0.007, R2 = 0.26) and then declined (p = 0.002, R2 = 0.18)(Figure 3).

Figure 2
Violin plot with two panels comparing soil properties across three substrates: Sand, Lichen, and Moss. Panel A shows soil pH, with higher variability in Sand. Panel B displays soil electrical conductivity (EC), also showing highest variability in Moss. Each plot includes a box plot within the violin plot to indicate data distribution.

Figure 2. Changes of soil pH (A) and EC (B) in interspaces of different BSC patches (Sand: bare sand, Lich: lichen crust, Moss: moss crust). Different letters (a, b) above the box plot indicate significant differences at p < 0.05. The sample size of sandy soil is 70, the sample size of algae and lichen crusts is 80, and the sample size of moss crusts is 220.

Figure 3
Violin plots labeled A and B comparing soil indices. A displays soil multifunctional index, B shows soil quality indices for sand, lichens, and moss. Colors differentiate categories with annotations: a, b, c. Both plots highlight distribution and median values.

Figure 3. Changes of SMF (A) and SQI (B) in interspaces of different BSC patches (Sand: bare sand, Lich: lichen crust, Moss: moss crust). Different letters (a–c) above the box plot indicate significant differences at p < 0.05. The sample size of sandy soil is 70, the sample size of algae and lichen crusts is 80, and the sample size of moss crusts is 220.

3.2 The influence of climate and biocrusts

The SMF (p < 0.001, R2 = 0.47) and SQI (p < 0.001, R2 = 0.62)increased with an rising biocrusts index (Figure 4). There was no obvious change in SVM with increasing biocrusts index. SMF and SQI rose with increasing MAP up to a peak and then declined (Figure 5). Both metrics were initially positively correlated with MAP (SMF: p = 0.027, R2 = 0.18; SQI: p = 0.023, R2 = 0.18), and both decreased when MAP > 163 mm (SMF: p = 0.013, R2 = 0.15; SQI: p = 0.0015, R2 = 0.24). SVM and MAP were negatively correlated (p = 0.02, R2 = 0.079). SMF and SQI also increased with MAT (SMF: p < 0.001, R2 = 0.16; SQI: p < 0.001, R2 = 0.22) and with MAW (SMF: p < 0.001, R2 = 0.28; SQI: p < 0.001, R2 = 0.33), whereas SVM showed no significant change with MAT (p > 0.05) or MAW (p > 0.05) (Figure 5).

Figure 4
Scatter plots show relationships between variables. Panels A, B, and G illustrate positive correlations of SMF and SQI with independent variables, with R-squared and p-values indicating statistical significance. Panel D presents dual regression lines for two datasets. Panels C, F, I depict less clear associations in SVM data, noting moderate R-squared in panel F. Gray bands indicate confidence intervals.

Figure 4. The trend of soil SMF, SQI and SVM with MAT, MAP and MAW increased. (A–C) The changes of SMF, SQI and SVM with MAT increased; (D–F) The changes of SMF, SQI, SVM with MAP increased; (G–I) The changes of SMF, SQI, SVM with MAW increased. Equations on each plot show the p and R2 of linear regression. The total sample size is 370.

Figure 5
Scatter plots showing relationships between soil metrics and soil pH or soil EC. Plot A: SMF vs Soil EC, \( R^2 = 0.33 \). Plot B: SQI vs Soil EC, \( R^2 = 0.55 \). Plot C: SVM vs Soil EC, \( R^2 = 0.12 \). Plot D: SMF vs Soil pH, \( R^2 = 0.4 \). Plot E: SQI vs Soil pH, \( R^2 = 0.5 \). Plot F: SVM vs Soil pH. Each plot includes a regression line with a confidence interval.

Figure 5. The trend of soil SMF, SQI and SVM with soil pH and EC increased. (A–C) The changes of SMF, SQI and SVM with soil pH increased: (D–F) The changes of SMF, SQI, SVM with soil EC increased. Equations on each plot show the p and R2 of linear regression. The total sample size is 370.

The SEM results showed that the model accounted for 21% of the variance in SVM (Figure 6). Spatial coordinates strongly influenced biocrusts (Lon: –0.58, Lat: –0.97), SQI (Lon: –0.20, Lat: –0.18) and SVM (Lon: –0.38, Lat: –0.23). MAT exerted a significant negative effect on biocrust development and distribution (-0.33). Biocrusts had significant positive effects on soil SMF (0.47) and SQI (0.31). Biocrusts did not significantly affect SVM. Soil SMF had a significant negative effect on SVM (- 0.50), whereas SQI had a significant positive effect on SVM (0.63).

Figure 6
Path diagram illustrating relationships among Spatial, Biocrusts, Climate, SQI, SMF, and SVM variables. Arrows indicate influence direction and correlation strength, with numbers representing correlation coefficients and R-squared values. Red, blue, and black arrows differentiate relationship types. Additional boxes detail correlations within specific parameter sets: Box1, Box2, and Box3, highlighting LAT, LON, MAT, MAP, MAW, and RAD variables.

Figure 6. Final fitted structural equation models representing relative effects of climate (MAT, MAP, MAW, Rad), geographical position (Lat, Lon), BSC (BSC indicate), soil pH, soil EC, soil sand content (Sand), soil SMF, SQI, and SVM. Boxes signify measured variables. Spatial: latitude (Lat) and longitude (Lon); Climate: include MAT, MAP, MAW, Rad. Standardized path coefficients are displayed, with the width of each arrow equivalent to the strength of the path. Red lines denote the negative paths (p < 0.05). Blue lines indicate positive paths (p < 0.05). The total amount of variance (R2) explained for each endogenous variable (those with arrows pointing to them) is given below the variable. Corresponding probability values are included when p< 0.05 (*p < 0.05, **p < 0.01, ***p < 0.001). The fit of the model was statistically tested (Fisheries’C = 39.833, df = 62, p = 0.987).

4 Discussion

4.1 Effect of biocrust types on SMF and SQI

Our findings confirm that biocrust type significantly modulates both SQI and SMF. Compared with algal–lichen crusts, moss-crusted soils exhibited markedly higher TN, NO3-N, NH4-N, AP, AK, pH, EC, SQI and SMF. Previous studies have demonstrated that as the biological crust layer transitions from algal-lichen crust to moss crust, the soil nutrients beneath the biological crust exhibit an increase (Gao et al., 2018; Fan et al., 2021). Thus, our results align with previous studies showing that biocrusts strongly modulate soil nutrients and physical properties within their colonized zones. One reason for this result is that the distribution of biocrusts in desert ecosystems stabilize the soil (Weber et al., 2022; Yang et al., 2022); The increase in soil stability could also lead to an increase in SMF and SQI. Furthermore, biocrusts play a critical role in regulating soil moisture dynamics by enhancing atmospheric water vapor condensation, decreasing evaporation, boosting water retention, and stabilizing soil surfaces. These effects are particularly pronounced in arid and semi-arid ecosystems, where biocrusts can markedly increase soil water retention and availability for plant growth (Li and Giora, 2021; Li and Xiao, 2022; Sun et al., 2021; Xiao et al., 2022, Xiao et al., 2019). The “fertilizer island” effect exhibited by biocrusts can improve nutrient availability and thus improve soil function (Agnelli et al., 2021; Li et al., 2019; Pushkareva et al., 2017; Young et al., 2022). Furthermore, the “fertilizer island” effect of biocrusts contributes to the accumulation of soil nutrients, which may cause nutrient accumulation of desert surface soil and increase the SMF and SQI.

Biocrusts are well-known for driving the nutrient cycle of desert soil (Wang et al., 2021; Young et al., 2022). Algae, lichens and mosses can increase soil stability by enhancing soil aggregation and reducing erosion (Lange and Belnap, 2016). Therefore, we suggest that the changes in soil SMF and SQI were due to the biocrusts type. Chamizo et al. (2017) validated that biocrusts removal enhanced sediment yield by 20 to 60 times more than that of undisturbed soils colonized by lichen. By reducing biocrust, desert ecosystems may experience a reduction in soil carbon, nitrogen fixation and phosphorus activation and mineralization (Gao et al., 2020). The soil SMF and SQI in moss crust were considerably higher than those in algae-lichen crust and sand. Moss crust, has high soil stability and nutrient accumulation, and can enhance soil nutrients, SMF and SQI in moss crust patches. Hence, we consider that biocrusts can influence SMF and SQI, but this effect is moderated by the biocrusts type and mostly occurs in moss crust.

4.2 Soil SMF, SVM, and SQI are regulated by climate and biocrusts

Desert ecosystems are extremely sensitive to climate change, therefore climate can significantly affect the vegetation distribution characteristics in desert ecosystems. Our results confirm this conclusion that climate change can significantly affect the distribution and developmental characteristics of biocrusts. The biocrusts index shows a significant increasing trend with the increase of MAT and MAW. Previous research has found that temperature have a significant impact on the development of biocrusts type (Ferrenberg et al., 2015). Moreover, the SMF, SQI and SVM increase with the increase of MAT. The apparent contradiction arises because simple regression captures the net temperature benefit, whereas SEM partitions it into opposing direct (negative) and indirect (positive) components. Temperature affects the microbial activity and metabolic processes in biocrusts, and a suitable temperature is beneficial for the transformation from algal-lichen crust to moss crust (Zelikova et al., 2012). Changes in temperature can alter soil moisture content, thereby affecting soil aeration, permeability, and cohesion (Li et al., 2018; Wei et al., 2022). Temperature changes can affect the rate of chemical reactions in soil, such as the release and fixation of mineral nutrients in soil. After accounting for other variables, warming exerts a direct inhibitory effect on microbial activity-primarily by accelerating evaporation and reducing soil water potential-which consequently diminishes SMF. At the same time, warming also indirectly enhances SMF through two mechanisms: prolonging the growing season and improving the utilization of condensed water, thereby promoting moss crust development. These opposing mechanisms operate simultaneously within the same plots and across the same climate gradient. Their net effect remains weakly positive due to partial cancellation of positive and negative contributions, which explains why simple regression and SEM both indicate an overall positive influence of warming on SMF, despite exhibiting opposing signs in their direct pathways. Under different temperature conditions, the types and quantities of organisms in the soil will vary, thereby affecting the SMF. This shift contributes to an increase in soil moisture retention, SQI, and soil stability.

Precipitation is an important environmental factor in desert ecosystems. Precipitation provides the necessary water for microorganisms in biocrusts, and usually the more precipitation, the better the development of different biocrusts types (Wu and Zhang, 2018). Therefore, biocrusts type affects SMF and SQI both directly and indirectly across a climatic gradient (Su et al., 2021). It directly alters SMF and SQI by increasing soil carbon and nitrogen fixation, and phosphorus desorption (Bowker et al., 2014; Xu et al., 2021). Moreover, the aggregation and growth of photoautotrophic organisms in surface soil provide effective coverage and soil carbon storage in areas where plant growth is limited owing to sporadic precipitation (Lefcheck et al., 2015; Maier et al., 2018). Biocrusts type development and the accompanying functional traits also indirectly affect SMF and SQI by altering the community abundance and composition of soil microbial communities (Garibotti et al., 2018; Su et al., 2021; Zhou et al., 2023).

We found that climate directly influenced the soil SMF and SQI. The SMF and SQI show a one-peak model with the increase of MAP. The SMF and SQI are the highest when the MAP is 163 mm (Figures 3, 5). The moss crusts increase with the increase of MAP in the desert. Simultaneously, some studies found that the relative positive effects of biocrust-forming mosses on multifunctionality compared with bare soil increased with increasing aridity, and that biocrusts is crucial to buffer negative effects of climate change on multifunctionality in global drylands (Delgado-Baquerizo et al., 2016). Precipitation affects soil moisture and permeability. Moderate water content helps dissolve and transport soil nutrients, but excessive precipitation may lead to soil erosion (Currier and Sala, 2022). The increase in precipitation may increase the leaching and reduce the nutrient concentration on the soil surface. The present study found that with the increase in precipitation, the microbial community structure of desert soil also increased first and then decreased, and the nutrient content of moss also increased. Furthermore, soil microbial communities and biocrusts species could regulate the effects of global change on SMF in drylands (Liu et al., 2017). All water inputs (rainfall, snow-melt, dew) ultimately derive from precipitation, so we treat them collectively as the precipitation driver. Below the 163 mm threshold every extra 1 mm of MAP increases SMF by 0.8% and SQI by 0.6% (Figure 4) because more frequent wet pulses prolong biocrust photosynthesis. Therefore, the response of the SMF and SQI to MAP is a complex process in drylands. This complexity is further compounded by the interactions between various biotic and abiotic factors. The divergence point (where SQI plateaus but SVM continues to rise) coincides with our 163 mm MAP threshold; future work should test whether this decoupling is universal across temperate deserts.

4.3 SMF, SQI is mainly driven by biocrusts, but is influenced by climate

Based on the results of the SEM, it can be concluded that SQI is mainly affected by the biocrusts and soil SMF. There is a significant negative effect of soil SMF on SVM, but a significant positive effect on SQI. Moreover, it is also found that the Lat and Lon of soil samples and climate can markedly influence biocrusts, SMF and SQI. Due to different soil substrates, SMF and SVM are significantly different between regions. In line with the previous studies, SMF and SQI are affected not only by the biocrusts but also by the climate (Ferrenberg et al., 2022; Huang et al., 2020; Hui et al., 2022; Xiao and Bowker, 2020). Biocrusts have been shown to play an essential role in controlling ecosystem responses to climate change (Maestre et al., 2013; Reed et al., 2012). For example, temperature, soil moisture, and nutrient limitation are thought to affect soil carbon metabolisms, such as microbial growth, respiration, carbon use efficiency, and microbial biomass turnover (Takriti et al., 2018; Hagerty et al., 2014; Dijkstra et al., 2015; Schindlbacher et al., 2015; Zheng et al., 2019). Temperature, soil pH, enzyme activity, and substrate quality can influence soil nutrient cycling (Noll et al., 2019; Wallenstein and Weintraub, 2008; Wanek et al., 2011). Therefore, the results demonstrate that SMF and SQI are substantially influenced by biotic factors (biocrusts type) and abiotic factors (climate). We consider that the spatial and climate in the desert ecosystem are the significant factors that directly affect the SMF, SQI, and SVM.

In this study, we focused on the impact of different types of biocrusts on soil SMF, SQI and SVM in the Gurbantunggut Desert. However, there are certain limitations to this study, as we concentrated solely on two types of biocrusts in the Gurbantunggut Desert. This study has several limitations. First, the findings are based on only two types of biocrusts (algal-lichen and moss crusts) within a single desert region, which may constrain the broader applicability of the results. Second, as biocrusts globally encompass a wider variety of forms, including pure algal, lichen, moss, and mixed assemblages, the functional diversity captured in this study may not be fully representative. Third, the lack of data on associated vegetation and soil microbial communities limits our understanding of their interactive effects with biocrusts on soil multifunctionality and quality.

Biocrusts represent some of the earliest pioneer species in vegetation succession, characterized by minimal resource consumption, and they hold significant potential for applications in desertification control, ecological conservation, and restoration. In the desert regions of northwest China, biocrusts exhibit extensive distribution, high species diversity, and a rich variety of types. However, comprehensive and comparative research on the diversity of biocrusts across different desert areas is lacking, as is a systematic understanding of their relationship with soil nutrient functions. Therefore, in the future, we will primarily investigate the composition and diversity of spore plants and microorganisms in various biocrust types in desert regions, particularly under the influence of global climate change, as well as their relationship with soil functions and quality. Moreover, we aim to enhance our understanding of the relationship between soil nutrient cycling characteristics and the physical and chemical properties of biocrusts, along with the sensitivity of different crust types to climatic factors. We will further explore the application of biocrusts and propose new insights into ecological restoration practices in arid regions, thereby providing valuable decision-making support for the management of desert ecosystems and the restoration of degraded ecosystems.

These limitations highlight the need for future studies to expand the scope of biocrust types examined, including pure algal crusts, cyanobacterial crusts, and mixed crustal assemblages, across multiple desert ecosystems with varying climatic conditions and soil substrates. Establish a synchronous observation network in temperate deserts (e.g., Taklamakan and Badain Jaran), integrating micrometeorological measurements and unmanned aerial vehicle remote sensing to build a dynamic, long-term (≥10 years) dataset on soil multifunctionality in biological crusts. This will help validate the universality of the 163 mm precipitation threshold identified in this study. Using metagenomics and 15N/13C isotope labeling techniques, we will quantify the expression levels of key functional genes (nifH, amoA, phoD) involved in carbon fixation, nitrogen fixation, and phosphorus mineralization across different crust types. This approach aims to clarify the microbial mechanisms through which biological crusts drive soil multifunctionality (SMF). Along the natural climate gradient on the eastern edge of the desert, field manipulation platforms will be deployed to simulate warming (+2°C, +4°C) and altered precipitation regimes (± 30%), combined with crust type replacement treatments. These experiments will allow quantitative distinction between the direct effects of climate and the indirect effects mediated by shifts in crust type on SMF, soil quality index (SQI), and soil vulnerability mitigation (SVM). In subsequent studies, metagenomic sequencing and qPCR will be applied to compare the abundance and expression patterns of key functional genes (e.g., nifH, amoA, phoD) related to carbon fixation, nitrogen fixation, and phosphorus transformation in algal-lichen and moss crusts. A random forest model will be employed to quantify the relative contribution of microbial community structural changes to SMF components, thereby elucidating the microbiological mechanisms by which biocrust type regulates SMF.

5 Conclusions

The present study demonstrates that biocrusts type can significantly and directly affect the soil SMF and SQI, and can also indirectly affect SVM by altering SMF and SQI. Moreover, the climate can also indirectly affect SMF, SQI and SVM by altering biocrusts type. This reinforces the notion that biocrusts are major structural and functional components in drylands. Our results emphasizes that it is crucial to assess them in different biocrusts types and climates if we aim to better understand their role in driving ecosystem function in highly heterogeneous ecosystems such as drylands. We, for the first time, quantified a MAP threshold of 163 mm as the peak inflection point between moss crust cover and SMF. This threshold offers a scientific basis for defining the precipitation boundary and the upper limit of artificial precipitation in ecological restoration efforts across temperate deserts in Central Asia. In these ecosystems, the probability of natural moss crust establishment is highest in areas with MAP between 145–180 mm and a MAT of 6-8°C. It is therefore recommended that artificial moss crust inoculation be prioritized within this climatic range. Due to the lack of information on vegetation and soil microbial communities in this study, it is not possible to calculate the relationships between vegetation, biocrusts and soil microorganisms, as well as their impact on soil SMF, SQI and SVM. We suggest using biocrusts as a medium for soil aboveground and underground, and further exploring the structure and function of desert ecosystems by combining vegetation and soil microbial diversity and community structure. Our research findings contribute to understanding the distribution of different biocrusts types in desert ecosystems and their impact on soil nutrients, functions, and quality, providing a scientific basis for ecological restoration and sustainable management of desert ecosystems.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author contributions

YGL: Conceptualization, Data curation, Formal analysis, Funding acquisition, Methodology, Resources, Writing – original draft, Writing – review & editing. YG: Data curation, Formal analysis, Writing – review & editing. YH: Conceptualization, Formal analysis, Investigation, Methodology, Writing – original draft. YXL: Conceptualization, Methodology, Resources, Writing – original draft. BY: Conceptualization, Investigation, Methodology, Resources, Writing – original draft. XZ: Conceptualization, Data curation, Resources, Writing – review & editing. HY: Methodology, Conceptualization, Resources, Wrting – review & editing. YZ: Conceptualization, Funding acquisition, Methodology, Resources, Writing – review & editing.

Funding

The author(s) declared that financial support was received for this work and/or its publication. This research was supported by the National Natural Science Foundation of China (42007099, U2003214) and the Xinjiang Tianshan Youth Talent Top Project (2022TSYCCX0007).

Acknowledgments

We thank YL and Yi-Xiang Sun, for their assistance with sample collecting in the field. We are grateful to the Associate Editor and reviewers for providing valuable comments.

Conflict of interest

The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declared that generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2026.1670208/full#supplementary-material

Abbreviations

SMF, Soil multifunctionality; SQI, Soil quality index; SVM, Spatial variability of soil multifunctionality; Biocrust, Biological soil crust; EC, Electrical conductivity; MAP, Mean annual precipitation; MAT, Mean annual temperature; SOC, Soil organic carbon; TN, Total nitrogen; TP, Total phosphorus; TK, Total potassium; AP, Available phosphorus; AK, Available potassium.

References

Agnelli, A., Corti, G., Massaccesi, L., Ventura, S., and D’Acqui, L. P. (2021). Impact of biological crusts on soil formation in polar ecosystems. Geoderma 401, 115340. doi: 10.1016/j.geoderma.2021.115340

Crossref Full Text | Google Scholar

Allington, G. R. H. and Valone, T. J. (2014). Islands of fertility: A byproduct of grazing? Ecosystems 17, 127–141. doi: 10.1007/s10021-013-9711-y

Crossref Full Text | Google Scholar

Andrews, S. S., Karlen, D. L., and Cambardella, C. A. (2004). The soil management assessment framework. Soil Sci. Soc. America J. 68, 1945–1962. doi: 10.2136/sssaj2004.1945

Crossref Full Text | Google Scholar

Andrews, S. S., Karlen, D. L., and Mitchell, J. P. (2002). A comparison of soil quality indexing methods for vegetable production systems in Northern California. Agriculture Ecosyst. Environ. 90, 25–45. doi: 10.1016/S0167-8809(01)00174-8

Crossref Full Text | Google Scholar

Barrett, J. E., Virginia, R. A., Lyons, W. B., McKnight, D. M., Priscu, J. C., Doran, P. T., et al. (2007). Biogeochemical stoichiometry of Antarctic Dry Valley ecosystems. J. Geophysical Research-Biogeosciences 112. doi: 10.1029/2005JG000141

Crossref Full Text | Google Scholar

Baumann, K., Eckhardt, K.-U., Acksel, A., Gros, P., Glaser, K., Gillespie, A. W., et al. (2021). Contribution of biological soil crusts to soil organic matter composition and stability in temperate forests. Soil Biol. Biochem. 160, 108315. doi: 10.1016/j.soilbio.2021.108315

Crossref Full Text | Google Scholar

Belnap, J., Phillips, S. L., Witwicki, D. L., and Miller, M. E. (2008). Visually assessing the level of development and soil surface stability of cyanobacterially dominated biological soil crusts. J. Arid Environments 72, 1257–1264. doi: 10.1016/j.jaridenv.2008.02.019

Crossref Full Text | Google Scholar

Bouma, J. (2002). Land quality indicators of sustainable land management across scales. Agric. Ecosyst. Environ. 88, 129–136. doi: 10.1016/S0167-8809(01)00248-1

Crossref Full Text | Google Scholar

Bowker, M. A., Maestre, F. T., Eldridge, D., Belnap, J., Castillo-Monroy, A., Escolar, C., et al. (2014). Biological soil crusts (biocrusts) as a model system in community, landscape and ecosystem ecology. Biodiversity Conserv. 23, 1619–1637. doi: 10.1007/s10531-014-0658-x

Crossref Full Text | Google Scholar

Bowker, M. A., Maestre, F. T., and Mau, R. L. (2013). Diversity and patch-size distributions of biological soil crusts regulate dryland ecosystem multifunctionality. Ecosystems 16, 923–933. doi: 10.1007/s10021-013-9644-5

Crossref Full Text | Google Scholar

Cantón, Y., Román, J. R., Chamizo, S., Rodríguez-Caballero, E., and Moro, M. J. (2014). Dynamics of organic carbon losses by water erosion after biocrust removal. J. Hydrology Hydromechanics 62, 258–268. doi: 10.2478/johh-2014-0033

Crossref Full Text | Google Scholar

Chamizo, S., Cantón, Y., Miralles, I., and Domingo, F. (2012). Biological soil crust development affects physicochemical characteristics of soil surface in semiarid ecosystems. Soil Biol. Biochem. 49, 96–105. doi: 10.1016/j.soilbio.2012.02.017

Crossref Full Text | Google Scholar

Chamizo, S., Rodríguez-Caballero, E., Román, J. R., and Cantón, Y. (2017). Effects of biocrust on soil erosion and organic carbon losses under natural rainfall. Catena 148, 117–125. doi: 10.1016/j.catena.2016.06.017

Crossref Full Text | Google Scholar

Couradeau, E., Karaoz, U., Lim, H. C., Ulisses, N. D. R., Northen, T., Brodie, E., et al. (2016). Bacteria increase arid-land soil surface temperature through the production of sunscreens. Nat. Commun. 7, 10373. doi: 10.1038/ncomms10373

PubMed Abstract | Crossref Full Text | Google Scholar

Currier, C. M. and Sala, O. E. (2022). Precipitation versus temperature as phenology controls in drylands. Ecology 103, e3793. doi: 10.1002/ecy.3793

PubMed Abstract | Crossref Full Text | Google Scholar

Day, K. J., Hutchings, M. J., and John, E. A. (2003). The effects of spatial pattern of nutrient supply on yield, structure and mortality in plant populations. J. Ecol. 91, 541–553. doi: 10.1046/j.1365-2745.2003.00799.x

Crossref Full Text | Google Scholar

Delgado-Baquerizo, M., Maestre, F. T., Eldridge, D. J., Bowker, M. A., Ochoa, V., Gozalo, B., et al. (2016). Biocrust-forming mosses mitigate the negative impacts of increasing aridity on ecosystem multifunctionality in drylands. New Phytol. 209, 1540–1552. doi: 10.1111/nph.13688

PubMed Abstract | Crossref Full Text | Google Scholar

Delgado-Baquerizo, M., Maestre, F. T., Gallardo, A., Bowker, M. A., Wallenstein, M. D., Quero, J. L., et al. (2013). Decoupling of soil nutrient cycles as a function of aridity in global drylands. Nature 502, 672–676. doi: 10.1038/nature12670

PubMed Abstract | Crossref Full Text | Google Scholar

Dijkstra, P., Salpas, E., Fairbanks, D., Miller, E. B., and Hagerty, S. (2015). High carbon use efficiency in soil microbial communities is related to balanced growth, not storage compound synthesis. Soil Biol. Biochem 89, 35–43. doi: 10.1016/j.soilbio.2015.06.021

Crossref Full Text | Google Scholar

Ding, J. and Eldridge, D. J. (2021). Climate and plants regulate the spatial variation in soil multifunctionality across a climatic gradient. Catena 201, 105233. doi: 10.1016/j.catena.2021.105233

Crossref Full Text | Google Scholar

Durán, J., Delgado-Baquerizo, M., Dougill, A. J., Guuroh, R. T., Linstdter, A., Thomas, A. D., et al. (2018). Temperature and aridity regulate spatial variability of soil multifunctionality in drylands across the globe. Ecology, 99, 1184–1193. doi: 10.1002/ecy.2199

PubMed Abstract | Crossref Full Text | Google Scholar

Eldridge, D. J., Delgado-Baquerizo, M., Quero, J. L., Ochoa, V., Gozalo, B., García-Palacios, P., et al. (2020). Surface indicators are correlated with soil multifunctionality in global drylands. J. Appl. Ecol. 57, 424–435. doi: 10.1111/1365-2664.13540

Crossref Full Text | Google Scholar

Ettema, C. H. and Wardle, D. A. (2002). Spatial soil ecology. Trends Ecol. Evol. 17, 0–183. doi: 10.1016/S0169-5347(02)02496-5

Crossref Full Text | Google Scholar

Fan, J., Li, S., Yu, H., and Huang, J. (2021). Soil enzyme activity and carbon, nitrogen and phosphorus stoichiometric characteristics under different types of biocrusts and subsoil in Mu Us Sandland. J. Desert Res. 41, 109–120. doi: 10.7522/j.issn.1000-694X.2021.00041

Crossref Full Text | Google Scholar

Ferrenberg, S., Reed, S. C., and Belnap, J. (2015). Climate change and physical disturbance cause similar community shifts in biological soil crusts. Proc. Natl. Acad. Sci. U.S.A. 112, 12116–12121. doi: 10.1073/pnas.1509150112

PubMed Abstract | Crossref Full Text | Google Scholar

Ferrenberg, S., Tucker, C. L., Reibold, R., Howell, A., and Reed, S. C. (2022). Quantifying the influence of different biocrust community states and their responses to warming temperatures on soil biogeochemistry in field and mesocosm studies. Geoderma 409, 115633. doi: 10.1016/j.geoderma.2021.115633

Crossref Full Text | Google Scholar

Fraterrigo, J. M. and Rusak, J. A. (2010). Disturbance-driven changes in the variability of ecological patterns and processes. Ecol. Lett. 11, 756–770. doi: 10.1111/j.1461-0248.2008.01191.x

PubMed Abstract | Crossref Full Text | Google Scholar

Gao, L., Zhao, Y., Xu, M., Sun, H., and Yang, Q. (2018). The effects of biological soil crust succession on soil ecological stoichiometry characteristics. Acta Ecologica Sinica 38, 678–688 doi: 10.5846/stxb201612132559.

Crossref Full Text | Google Scholar

Gao, L., Bowker, M. A., Sun, H., Zhao, J., and Zhao, Y. (2020). Linkages between biocrust development and water erosion and implications for erosion model implementation. Geoderma 357, 113973. doi: 10.1016/j.geoderma.2019.113973

Crossref Full Text | Google Scholar

Garcia-Callejas, D., Godoy, O., Buche, L., Hurtado, M., Lanuza, J. B., Allen-Perkins, A., et al. (2023). Non-random interactions within and across guilds shape the potential to coexist in multi-trophic ecological communities. Ecol. Lett. 26, 831–842. doi: 10.1111/ele.14206

PubMed Abstract | Crossref Full Text | Google Scholar

Garibotti, I. A., Gonzalez Polo, M., Tabeni, S., and Rasmann, S. (2018). Linking biological soil crust attributes to the multifunctionality of vegetated patches and interspaces in a semiarid shrubland. Funct. Ecol. 32, 1065–1078. doi: 10.1111/1365-2435.13044

Crossref Full Text | Google Scholar

Grote, E. E., Belnap, J., Housman, D. C., and Sparks, J. P. (2010). Carbon exchange in biological soil crust communities under differential temperatures and soil water contents: implications for global change. Global Change Biol. 16, 2763–2774. doi: 10.1111/j.1365-2486.2010.02201.x

Crossref Full Text | Google Scholar

Guerra, C. A., Heintz-Buschart, A., Sikorski, J., Chatzinotas, A., Guerrero-Ramírez, N., Cesarz, S., et al. (2020). Blind spots in global soil biodiversity and ecosystem function research. Nat. Commun. 11, 3870. doi: 10.1038/s41467-020-17688-2

PubMed Abstract | Crossref Full Text | Google Scholar

Hagerty, S. B., Van Groenigen, K. J., Allison, S. D., Hungate, B. A., Schwartz, E., Koch, G. W., et al. (2014). Accelerated microbial turnover but constant growth efficiency with warming in soil. Nat. Climate Change 4, 903–906. doi: 10.1038/nclimate2361

Crossref Full Text | Google Scholar

Hagh-Doust, N., Mikryukov, V., Anslan, S., Bahram, M., Puusepp, R., Dulya, O., et al. (2023). Effects of nitrogen deposition on carbon and nutrient cycling along a natural soil acidity gradient as revealed by metagenomics. New Phytol. 238, 2607–2620. doi: 10.1111/nph.18897

PubMed Abstract | Crossref Full Text | Google Scholar

Hedo de Santiago, J., Lucas-Borja, M. E., Wic-Baena, C., Andrés-Abellán, M., and de las Heras, J. (2016). Effects of thinning and induced drought on microbiological soil properties and plant species diversity at dry and semiarid locations. Land Degradation Dev. 27, 1151–1162. doi: 10.1002/ldr.2361

Crossref Full Text | Google Scholar

Housman, D. C., Powers, H. H., Collins, A. D., and Belnap, J. (2006). Carbon and nitrogen fixation differ between successional stages of biological soil crusts in the Colorado Plateau and Chihuahuan Desert. J. Arid Environments 66, 620–634. doi: 10.1016/j.jaridenv.2005.11.014

Crossref Full Text | Google Scholar

Hu, M., Yan, R., Wu, H., Ni, R., Zhang, D., and Zou, S. (2023). Linking soil phosphorus availability and phosphatase functional genes to coastal marsh erosion: Implications for nutrient cycling and wetland restoration. Sci. Total Environ. 898, 165559. doi: 10.1016/j.scitotenv.2023.165559

PubMed Abstract | Crossref Full Text | Google Scholar

Huang, J., Zhang, G., Zhang, Y., Guan, X., and Guo, R. (2020). Global desertification vulnerability to climate change and human activities. Land Degradation Dev. 31, 1380–1391. doi: 10.1002/ldr.3556

Crossref Full Text | Google Scholar

Huang, Y.-J., Li, Y.-G., Zhou, X.-B., Yin, B.-F., Tao, Y., and Zhang, Y.-M. (2023). Moss patch size as a factor profoundly influencing soil nutrient characteristics and multifunctionality of temperate desert in Central Asia. Ecol. Indic. 155. doi: 10.1016/j.ecolind.2023.110975

Crossref Full Text | Google Scholar

Hui, R., Zhao, R., Liu, L., and Li, X. (2022). Effect of snow cover on water content, carbon and nutrient availability, and microbial biomass in complexes of biological soil crusts and subcrust soil in the desert. Geoderma 406, 115505. doi: 10.1016/j.geoderma.2021.115505

Crossref Full Text | Google Scholar

Huo, R., Wang, J., Wang, K., Zhang, Y., Ren, T., Li, X., et al. (2024). Long-term straw return enhanced crop yield by improving ecosystem multifunctionality and soil quality under triple rotation system: An evidence from a 15 years study. Field Crops Res. 312, 109395. doi: 10.1016/j.fcr.2024.109395

Crossref Full Text | Google Scholar

Juhos, K., Szabó, S., and Ladányi, M. (2016). Explore the influence of soil quality on crop yield using statistically-derived pedological indicators. Ecol. Indic. 63, 366–373. doi: 10.1016/j.ecolind.2015.12.029

Crossref Full Text | Google Scholar

Kafei, F., Rezapour, S., Dalalian, M. R., Sabbaghtazeh, E., and Rafieyan, O. (2023). Soil quality index as affected by long-time continuous cultivation in a Mediterranean sub-humid region. Rendiconti Lincei. Sci. Fisiche e Naturali 34, 563–575.

Google Scholar

Karlen, D. L., Tomer, M. D., Neppel, J., and Cambardella, C. A. (2008). A preliminary watershed scale soil quality assessment in north central Iowa, USA. Soil Till Res. 99, 291–299. doi: 10.1016/j.still.2008.03.002

Crossref Full Text | Google Scholar

Khosravi, A. K., Rezapour, S., Asadzadeh, F., and Nouri, A. (2023). An integrated approach for estimating soil health: Incorporating digital elevation models and remote sensing of vegetation. Comput. Electron. Agric. 210, 107922. doi: 10.1016/j.compag.2023.107922

Crossref Full Text | Google Scholar

Kidron, G. J. and Büdel, B. (2014). Contrasting hydrological response of coastal and desert biocrusts. Hydrological Processes 28, 361–371. doi: 10.1002/hyp.9587

Crossref Full Text | Google Scholar

Kidron, G. J., Vonshak, A., and Abeliovich, A. (2010). Microbiotic crusts as biomarkers for surface stability and wetness duration in the Negev Desert. Earth Surface Processes Landforms 34, 1594–1604. doi: 10.1002/esp.1843

Crossref Full Text | Google Scholar

Lal, R., Monger, C., Nave, L., and Smith, P. (2021). The role of soil in regulation of climate. Philos. Trans. R. Soc. B: Biol. Sci. 376, 20210084.

Google Scholar

Lan, S., Wu, L., Zhang, D., and Hu, C. (2012). Composition of photosynthetic organisms and diurnal changes of photosynthetic efficiency in algae and moss crusts. Plant Soil 351, 325–336. doi: 10.1007/s11104-011-0966-9

Crossref Full Text | Google Scholar

Lan, S., Wu, L., Zhang, D., Hu, C., Lan, S., Wu, L., et al. (2013). Assessing level of development and successional stages in biological soil crusts with biological indicators. Microbial Ecol. 66, 394–403. doi: 10.1007/s00248-013-0191-6

PubMed Abstract | Crossref Full Text | Google Scholar

Lange, O. L. and Belnap, J. (2016). How biological soil crusts became recognized as a functional unit: a selective history, Biological Soil Crusts: An Organizing Principle in Drylands (Cham: Springer), 15–33.

Google Scholar

Langhans, T. M., Storm, C., and Schwabe, A. (2009). Biological soil crusts and their microenvironment: Impact on emergence, survival and establishment of seedlings. Flora 204, 157–168. doi: 10.1016/j.flora.2008.01.001

Crossref Full Text | Google Scholar

Lefcheck, J. S., Byrnes, J. E. K., Isbell, F., Gamfeldt, L., Griffin, J. N., Eisenhauer, N., et al. (2015). Biodiversity enhances ecosystem multifunctionality across trophic levels and habitats. Nat. Commun. 6, 6936. doi: 10.1038/ncomms7936

PubMed Abstract | Crossref Full Text | Google Scholar

Li, F. K. and Giora, J. (2021). Moss-dominated biocrusts enhance water vapor sorption capacity of surface soil and increase non-rainfall water deposition in drylands. Geoderma: Int. J. Soil Sci. 388, 114930. doi: 114930.10.1016/j.geoderma.2021.

Google Scholar

Li, Y. G., Huang, Y. J., Yin, B. F., Zhou, X. B., and Zhang, Y. M. (2024). Shrub and patch size of moss crusts regulate soil multifunctionality in a temperate desert of Central Asia. Soil Sci. Soc. America J. 88, 890–904. doi: 10.1002/saj2.20655

Crossref Full Text | Google Scholar

Li, H., Huo, D., Wang, W., Chen, Y., Cheng, X., Yu, G., et al. (2020). Multifunctionality of biocrusts is positively predicted by network topologies consistent with interspecies facilitation. Mol. Ecol. 29, 1560–1573. doi: 10.1111/mec.15424

PubMed Abstract | Crossref Full Text | Google Scholar

Li, X. R., Jia, X. H., Long, L. Q., and Zerbe, S. (2005). Effects of biological soil crusts on seed bank, germination and establishment of two annual plant species in the tengger desert (N China). Plant Soil 277, 375–385. doi: 10.1007/s11104-005-8162-4

Crossref Full Text | Google Scholar

Li, X. R., Jia, R. L., Zhang, Z. S., Zhang, P., and Hui, R. (2018). Hydrological response of biological soil crusts to global warming: A ten-year simulative study. Glob Chang Biol. 24, 4960–4971. doi: 10.1111/gcb.14378

PubMed Abstract | Crossref Full Text | Google Scholar

Li, J., Okin, G. S., Alvarez, L., and Epstein, H. (2007). Quantitative effects of vegetation cover on wind erosion and soil nutrient loss in a desert grassland of southern New Mexico, USA. Biogeochemistry 85, 317–332. doi: 10.1007/s10533-007-9142-y

Crossref Full Text | Google Scholar

Li, S. L. and Xiao, B. (2022). Cyanobacteria and moss biocrusts increase evaporation by regulating surface soil moisture and temperature on the northern Loess Plateau, China. CATENA, 212, 106068. doi: 10.1016/j.catena.2022.106068

Crossref Full Text | Google Scholar

Li, Y.-G., Zhou, X.-B., Lu, Y., Zhang, Y.-M., and Cheng, X. (2023). Moss C, N, P and K stoichiometry and their relationships are related to soil nutrients and environment in a temperate desert of central Asia. J. Plant Ecol. 16, 0–rtac070. doi: 10.1093/jpe/rtac070

Crossref Full Text | Google Scholar

Li, Y. G., Zhou, X. B., and Zhang, Y. M. (2019). Moss patch size and microhabitats influence stoichiometry of moss crusts in a temperate desert, Central Asia. Plant Soil 443, 55–72. doi: 10.1007/s11104-019-04191-x

Crossref Full Text | Google Scholar

Liu, Y. R., Delgado-Baquerizo, M., Trivedi, P., He, J.-Z., Wang, J.-T., and Singh, B. K. (2017). Identity of biocrust species and microbial communities drive the response of soil multifunctionality to simulated global change. Soil Biol. Biochem. 107, 208–217. doi: 10.1016/j.soilbio.2016.12.003

Crossref Full Text | Google Scholar

Maestre, F. T., Escolar, C., de Guevara, M. L., Quero, J. L., Lazaro, R., Delgado-Baquerizo, M., et al. (2013). Changes in biocrust cover drive carbon cycle responses to climate change in drylands. Glob Chang Biol. 19, 3835–3847. doi: 10.1111/gcb.12306

PubMed Abstract | Crossref Full Text | Google Scholar

Maier, S., Tamm, A., Wu, D., Caesar, J., Grube, M., and Weber, B. (2018). Photoautotrophic organisms control microbial abundance, diversity, and physiology in different types of biological soil crusts. Isme J. 12, 1032–1046. doi: 10.1038/s41396-018-0062-8

PubMed Abstract | Crossref Full Text | Google Scholar

Masto, R. E., Chhonkar, P. K., Singh, D., and Patra, A. K. (2008). Alternative soil quality indices for evaluating the effect of intensive cropping, fertilisation and manuring for 31 years in the semi-arid soils of India. Environ. Monit Assess. 136, 419–435. doi: 10.1007/s10661-007-9697-z

PubMed Abstract | Crossref Full Text | Google Scholar

Navas Romero, A. L., Herrera Moratta, M. A., Martinez Carretero, E., Rodriguez, R. A., and Vento, B. (2020). Spatial distribution of biological soil crusts along an aridity gradient in the central-west of Argentina. J. Arid Environments 176, 104099. doi: 10.1016/j.jaridenv.2020.104099

Crossref Full Text | Google Scholar

Niu, J., Yang, K., Tang, Z., and Wang, Y. (2017). Relationships between soil crust development and soil properties in the desert region of north China. Sustainability 9, 725. doi: 10.3390/su9050725

Crossref Full Text | Google Scholar

Noll, L., Zhang, S., Zheng, Q., Hu, Y., and Wanek, W. (2019). Wide-spread limitation of soil organic nitrogen transformations by substrate availability and not by extracellular enzyme content. Soil Biol. Biochem. 133, 37–49. doi: 10.1016/j.soilbio.2019.02.016

PubMed Abstract | Crossref Full Text | Google Scholar

Ochoa-Hueso, R., Eldridge, D. J., Delgado-Baquerizo, M., Soliveres, S., Bowker, M. A., Gross, N., et al. (2018). Soil fungal abundance and plant functional traits drive fertile island formation in global drylands. J. Ecol. 106, 242–253. doi: 10.1111/1365-2745.12871

Crossref Full Text | Google Scholar

Pushkareva, E., Kvíderová, J., Ŝimek, M., and Elster, J. (2017). Nitrogen fixation and diurnal changes of photosynthetic activity in Arctic soil crusts at different development stage. Eur. J. Soil Biol. 79, 21–30. doi: 10.1016/j.ejsobi.2017.02.002

Crossref Full Text | Google Scholar

Qi, Y., Darilek, J. L., Huang, B., Zhao, Y., Sun, W., and Gu, Z. (2009). Evaluating soil quality indices in an agricultural region of Jiangsu Province, China. Geoderma 149, 325–334. doi: 10.1016/j.geoderma.2008.12.015

Crossref Full Text | Google Scholar

Raiesi, F. (2017). A minimum data set and soil quality index to quantify the effect of land use conversion on soil quality and degradation in native rangelands of upland arid and semiarid regions. Ecol. Indic. 75, 307–320. doi: 10.1016/j.ecolind.2016.12.049

Crossref Full Text | Google Scholar

Reed, S. C., Coe, K. K., Sparks, J. P., Housman, D. C., Zelikova, T. J., and Belnap, J. (2012). Changes to dryland rainfall result in rapid moss mortality and altered soil fertility. Nat. Climate Change 2, 752–755. doi: 10.1038/nclimate1596

Crossref Full Text | Google Scholar

Reynolds, J. F., Smith, D. M. S., Lambin, E. F., Ii, B. L. T., Mortimore, M., Batterbury, S. P. J., et al. (2007). Global desertification: building a science for dryland development. Science 316, 847–851. doi: 10.1126/science.1131634

PubMed Abstract | Crossref Full Text | Google Scholar

Rodriguez-Caballero, E., Belnap, J., Büdel, B., Crutzen, P. J., Andreae, M. O., Pschl, U., et al. (2018). Dryland photoautotrophic soil surface communities endangered by global change. Nat. Geosci. 11, 185–189. doi: 10.1038/s41561-018-0072-1

Crossref Full Text | Google Scholar

Rutherford, W. A., Painter, T. H., Ferrenberg, S., Okin, G. S., and Flagg, C. (2017). Albedo feedbacks to future climate via climate change impacts on dryland biocrusts. Sci. Rep. 7, 44188. doi: 10.1038/srep44188

PubMed Abstract | Crossref Full Text | Google Scholar

Schindlbacher, A., Schnecker, J., Takriti, M., Borken, W., and Wanek, W. (2015). Microbial physiology and soil CO2 efflux after 9 years of soil warming in a temperate forest - no indications for thermal adaptations. Global Change Biol. 21, 4265–4277. doi: 10.1111/gcb.12996

PubMed Abstract | Crossref Full Text | Google Scholar

Schlesinger, W. H., Reynolds, J. F., Cunningham, G. L., Huenneke, L. F., Jarrell, W. M., Virginia, R. A., et al. (1990). Biological feedbacks in global desertification. Science 247, 1043–1048. doi: 10.1126/science.247.4946.1043

PubMed Abstract | Crossref Full Text | Google Scholar

Su, Y. G., Liu, J., Zhang, Y. M., and Huang, G. (2021). More drought leads to a greater significance of biocrusts to soil multifunctionality. Funct. Ecol. 35, 989–1000. doi: 10.1111/1365-2435.13761

Crossref Full Text | Google Scholar

Sun, F., Xiao, B., Li, S., and Kidron, G. J. (2021). Towards moss biocrust effects on surface soil water holding capacity: Soil water retention curve analysis and modeling. Geoderma 399, 115120. doi: 10.1016/j.geoderma.2021.115120

Crossref Full Text | Google Scholar

Takriti, M., Wild, B., Schnecker, J., Mooshammer, M., Knoltsch, A., Lashchinskiy, N., et al. (2018). Soil organic matter quality exerts a stronger control than stoichiometry on microbial substrate use efficiency along a latitudinal transect. Soil Biol. Biochem. 121, 212–220. doi: 10.1016/j.soilbio.2018.02.022

Crossref Full Text | Google Scholar

Tao, Y., Zhou, X. B., Zhang, S. H., Lu, H. Y., and Shao, H. (2020). Soil nutrient stoichiometry on linear sand dunes from a temperate desert in Central Asia. Catena 195, 104847. doi: 10.1016/j.catena.2020.104847

Crossref Full Text | Google Scholar

Tesfahunegn, G. B. (2014). Soil quality assessment strategies for evaluating soil degradation in northern Ethiopia. Appl. Environ. Soil Sci. 2014, 646502. doi: 10.1155/2014/646502

Crossref Full Text | Google Scholar

Vicente-Serrano, S. M., Zouber, A., Lasanta, T., and Pueyo, Y. (2012). Dryness is accelerating degradation of vulnerable shrublands in semiarid Mediterranean environments. Ecol. Monogr. 82, 407–428. doi: 10.1890/11-2164.1

Crossref Full Text | Google Scholar

Wallenstein, M. and Weintraub, M. (2008). Emerging tools for measuring and modeling the in situ activity of soil extracellular enzymes. Soil Biol. Biochem. 40, 2098–2106. doi: 10.1016/j.soilbio.2008.01.024

Crossref Full Text | Google Scholar

Wanek, W., Mooshammer, M., Blöchl, A., Hanreich, A., Keiblinger, K., Zechmeister-Boltenstern, S., et al. (2011). Determination of gross rates of amino acid production and immobilization in decomposing leaf litter by a novel 15N isotope pool dilution technique. Soil Biol. Biochem. 42, 1293–1302. doi: 10.1016/j.soilbio.2010.04.001

Crossref Full Text | Google Scholar

Wang, B., Liu, J., Zhang, X., and Wang, C. (2021). Changes in soil carbon sequestration and emission in different succession stages of biological soil crusts in a sand-binding area. Carbon Balance Manag 16, 27. doi: 10.1186/s13021-021-00190-7

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, P., Wang, S., Chen, F., Zhang, T., and Kong, W. (2024). Preparation of two types plant biochars and application in soil quality improvement. Sci. Total Environ. 906, 167334. doi: 10.1016/j.scitotenv.2023.167334

PubMed Abstract | Crossref Full Text | Google Scholar

Weber, B., Belnap, J., Budel, B., Antoninka, A. J., Barger, N. N., Chaudhary, V. B., et al. (2022). What is a biocrust? A refined, contemporary definition for a broadening research community. . Biol. Rev. Camb Philos. Soc 97, 1768–1785. doi: 10.1111/brv.12862

PubMed Abstract | Crossref Full Text | Google Scholar

Wei, X., Qin, F., Han, B., Zhou, H. K., Liu, M., and Shao, X. (2022). Spatial variations of bacterial communities associated with biological soil crusts along a climatic gradient in alpine grassland ecosystems. Plant Soil 480, 493–506. doi: 10.1007/s11104-022-05595-y

Crossref Full Text | Google Scholar

Wu, L. and Zhang, Y. (2018). Precipitation and soil particle size co-determine spatial distribution of biological soil crusts in the Gurbantunggut Desert, China. J. Arid Land 10, 701–711. doi: 10.1007/s40333-018-0065-3

Crossref Full Text | Google Scholar

Wu, N., Zhang, Y. M., and Downing, A. (2009). Comparative study of nitrogenase activity in different types of biological soil crusts in the Gurbantunggut Desert, Northwestern China. J. Arid Environments 73, 828–833. doi: 10.1016/j.jaridenv.2009.04.002

Crossref Full Text | Google Scholar

Xiao, B. and Bowker, M. A. (2020). Moss-biocrusts strongly decrease soil surface albedo, altering land-surface energy balance in a dryland ecosystem. Sci. Total Environ. 741, 140425. doi: 10.1016/j.scitotenv.2020.140425

PubMed Abstract | Crossref Full Text | Google Scholar

Xiao, B., Bowker, M. A., Zhao, Y. E., Chamizo, S., and Issa, O. M. (2022). Biocrusts: Engineers and architects of surface soil properties, functions, and processes in dryland ecosystems. Geoderma 424, 116015. doi: 10.1016/j.geoderma.2022.116015

Crossref Full Text | Google Scholar

Xiao, B., Sun, F., Hu, K., and Kidron, G. J. (2019). Biocrusts reduce surface soil infiltrability and impede soil water infiltration under tension and ponding conditions in dryland ecosystem. J. Hydrology 568, 792–802. doi: 10.1016/j.jhydrol.2018.11.051

Crossref Full Text | Google Scholar

Xu, L., Zhang, B., Wang, E., Zhu, B., Yao, M., Li, C., et al. (2021). Soil total organic carbon/total nitrogen ratio as a key driver deterministically shapes diazotrophic community assemblages during the succession of biological soil crusts. Soil Ecol. Lett. 3, 328–341. doi: 10.1007/s42832-020-0075-x

Crossref Full Text | Google Scholar

Xu, L., Zhu, B., Li, C., Yao, M., Zhang, B., and Li, X. (2020). Development of biological soil crust prompts convergent succession of prokaryotic communities. Catena 187, 104360. doi: 10.1016/j.catena.2019.104360

Crossref Full Text | Google Scholar

Yang, R. P., Yang, Z.-P., Yang, S., Chen, L.-L., Xin, J., Xu, L., et al. (2023). Nitrogen inhibitors improve soil ecosystem multifunctionality by enhancing soil quality and alleviating microbial nitrogen limitation. Sci. total Environ., 880, 163238. doi: 10.1016/j.scitotenv.2023.163238

PubMed Abstract | Crossref Full Text | Google Scholar

Yang, K., Zhao, Y., Gao, L., Sun, H., and Gu, K. (2022). Nonlinear response of hydrodynamic and soil erosive behaviors to biocrust coverage in drylands. Geoderma 405, 115457. doi: 10.1016/j.geoderma.2021.115457

Crossref Full Text | Google Scholar

Young, K. E., Ferrenberg, S., Reibold, R., Reed, S. C., Swenson, T., Northen, T., et al. (2022). Vertical movement of soluble carbon and nutrients from biocrusts to subsurface mineral soils. Geoderma 405, 115495. doi: 10.1016/j.geoderma.2021.115495

Crossref Full Text | Google Scholar

Zelikova, T. J., Housman, D. C., Grote, E. E., Neher, D. A., and Belnap, J. (2012). Warming and increased precipitation frequency on the Colorado Plateau: implications for biological soil crusts and soil processes. Plant Soil 355, 265–282. doi: 10.1007/s11104-011-1097-z

Crossref Full Text | Google Scholar

Zhang, Y. M., Chen, J., Wang, L., Wang, X. Q., and Gu, Z. H. (2007). The spatial distribution patterns of biological soil crusts in the Gurbantunggut Desert, Northern Xinjiang, China. J. Arid Environments 68, 599–610. doi: 10.1016/j.jaridenv.2006.06.012

Crossref Full Text | Google Scholar

Zhang, B., Kong, W., Wu, N., and Zhang, Y. (2016). Bacterial diversity and community along the succession of biological soil crusts in the Gurbantunggut Desert, Northern China. J. Basic Microbiol. 56, 670–679. doi: 10.1002/jobm.201500751

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, B., Li, Y., and Yuanming, X. (2018). Successional changes of fungal communities along the biocrust development stages. Biol. Fertility Soils Cooperating J. Int. Soc. Soil Science. 54, 285–294. doi: 10.1007/s00374-017-1259-0

Crossref Full Text | Google Scholar

Zhang, Y., Gao, M., Yu, C., Zhang, H., Yan, N., Wu, Q., et al. (2022). Soil nutrients, enzyme activities, and microbial communities differ among biocrust types and soil layers in a degraded karst ecosystem. Catena. 212, 106057. doi: 10.1016/j.catena.2022.106057

Crossref Full Text | Google Scholar

Zheng, Q., Hu, Y., Zhang, S., Noll, L., Böckle, T., Richter, A., et al. (2019). Growth explains microbial carbon use efficiency across soils differing in land use and geology. Soil Biol. Biochem. 128, 45–55. doi: 10.1016/j.soilbio.2018.10.006

PubMed Abstract | Crossref Full Text | Google Scholar

Zhou, H., Li, L., and Liu, Y. (2023). Biological soil crust development affects bacterial communities in the Caragana microphylla community in alpine sandy areas. Front. Microbiol. 14, 1106739. doi: 10.3389/fmicb.2023.1106739

PubMed Abstract | Crossref Full Text | Google Scholar

Zornoza, R., Acosta, J. A., Bastida, F., Domínguez, S. G., Toledo, D. M., and Faz, A. (2015). Identification of sensitive indicators to assess the interrelationship between soil quality, management practices and human health. SOIL 1, 173–185. doi: 10.5194/soil-1-173-2015

Crossref Full Text | Google Scholar

Keywords: biocrusts, desert, soil multifunctionality, soil quality index, spatial variability of soil multifunctionality

Citation: Li Y, Gao Y, Huang Y, Lu Y, Yin B, Zhou X, Yu H and Zhang Y (2026) Climate and biocrust types jointly regulate soil multifunctionality and quality in drylands: evidence from the Gurbantunggut Desert. Front. Plant Sci. 17:1670208. doi: 10.3389/fpls.2026.1670208

Received: 21 July 2025; Accepted: 20 January 2026; Revised: 25 December 2025;
Published: 10 February 2026.

Edited by:

David R Elliott, University of Derby, United Kingdom

Reviewed by:

Pengshuai Shao, Shandong University of Aeronautics, China
Guilin Wu, Chinese Academy of Sciences (CAS), China

Copyright © 2026 Li, Gao, Huang, Lu, Yin, Zhou, Yu and Zhang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Yonggang Li, bGl5b25nZ2FuZ0BoaXN0LmVkdS5jbg==; Yuanming Zhang, eW16aGFuZ0Btcy54amIuYWMuY24=

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.