Genesis and soil environmental implications of intact in-situ rhizoliths in dunes of the Badain Jaran Desert, northwestern China

Qingfeng  Sun · Kazem  Zamanian · Arnaud  Huguet · Omid  Bayat · Hong  Wang · Hanan S.  Badawy

Abstract Desert rhizoliths are generally found as weathered, broken and scattered samples on dune f ield surface, but rarely insitu in their initial states buried under the soil of desert in the Badain Jaran Desert, northwest China. This study off ers an assessment of the morphological, mineralogical, and chemical properties of intact and in-situ rhizoliths found in soils of swales and depressions among dune chains.The characteristics of these rare and precious objects were assessed using optical polarizing microscopy, cathodoluminescence, scanning electronic microscopy, radiocarbon dating, and stable isotopic analyses, providing the opportunity for discussion of the rhizolith formation mechanisms and associated environmental conditions. Field and laboratory investigations showed that the in-situ intact rhizoliths were formed only in the places where Artemisia shrubs are living, and the remaining root relicts within rhizoliths belong to this species. The spatial distribution of rhizoliths also suggested that low topographic positions on a landscape provided soil moisture, and redox environments favored rhizolith formation. A semi-closed redox environment in the subsoil at swales and depressions, where water is always present, along with the sandy soil texture, facilitated fast water percolation to deeper depths and condensation. Such a soil environment not only provides water for Artemisia growth, but also for the weathering of minerals such as felspars and calcite from primary carbonates, and for the decomposition of root relicts. Furthermore, harsh climatic conditions, such as strong winds and solar radiation, led to water evaporation through dead root channels and triggered the calcif ication along the root relicts. The entrapped lithogenic carbonates and to a lesser extent the decomposition of Artemisia roots provided the carbon sources for the rhizoliths formation, while the weathering of soil minerals,particularly feldspars and carbonates, was the main source of Ca. Rhizoliths in the Badain Jaran desert formed relatively quickly, probably over a few soil drying episodes.This led to the entrapment of a large quantity of lithogenic carbonates (more than 90% of carbon) within rhizolith cement. The re-dissolution of the entrapped lithogenic carbonates in rhizolith tubes should be taken into account in the paleoenvironmental interpretation of 14 C ages, the latter suggesting that rhizoliths formed during the Holocene(～ 2053 years cal BP, based on root organic relicts).

Keywords Rhizoliths · Calcif ication · Leptic regosols ·Artemisia roots · Decomposition · Soil moisture

1 Introduction

Soil is at the heart of the Earth’s critical zone (Chorover et al. 2007). Biomineralization, including pedogenic carbonate formation, is a common feature of soil evolution. It occurs not only in soil but also around or within roots, separately or concomitantly. Carbonate rhizoliths are one type of pedogenic carbonates (Lambers et al. 2009) where mineralization takes place around roots in modern (Holocene) soils and paleosols (Sun et al. 2019a, b, 2020; Huguet et al. 2020;Brazier et al. 2020).

Rhizoliths are generally def ined as ichnofossils (root traces) or fossils (petrif ied roots) in geology, taphonomy,and paleontology (Kraus and Hasiotis 2006; Barta, 2011).Phytocretions (minerals precipitate form an external mold of a plant, other than the root. Liutkus 2009) are similar to rhizoliths and formed around plant stem in standing water above sediments. In arid environments, the f luctuations in water table levels, pedogenesis, and erosive dynamics expose the plant remains to subaerial weathering, resulting in low rates of preservation of plant macrofossils (Gastaldo and Demko 2011) and explaining the rare presence of intact rhizoliths in such settings. However, rhizoliths are important indicators of paleo- and modern plants (Nascimento et al. 2019), especially in desert soils. Rhizoliths can be used as environmental markers of humidity, drainage, and local topography (Kraus and Hasiotis 2006; Li et al. 2015).

The formation mechanisms of rhizoliths have been widely discussed (e.g., Liutkus et al. 2005; Kraus and Hasiotis 2006; Owen et al. 2008; Gocke et al. 2010, 2011,2014; Bojanowski et al., 2016 ; Li et al. 2015a, b; Nascimento et al. 2019; Golubtsov et al., 2019) and remain open to question. This is partly due to the diffi culties involved in probing the rhizosphere – the microscale environment immediately surrounding root tissue, and also due to the numerous factors which can aff ect rhizoliths formation.Such factors include, for example, the nature of the cementing mineral, soil humidity and microbial activity (McLaren 1995), root exudation of organic acids (Albalasmeh and Ghezzehei 2014; Marschner 1995), mass-f low (Cramer and Hawkins 2009) and evaporative concentration and evapotranspiration (Owen et al. 2008). Moreover, other factors like the interplay of organic matter, bio-chemicals of production and/or decomposition of living and/or deceased organic tissues, mineral weathering, mineral-water-root interfaces, redox conditions, microorganism, gas exchange,pH, as well as soil physical properties can also aff ect rhizolith formation (Zhao et al. 2020; Dontsova et al. 2020; Liang et al. 2018; Dwivedi et al., 2017; Spohn et al. 2013; Sanaullah et al. 2011; Jones et al. 2009; Smits et al. 2005; Rasse et al. 2005; Hinsinger et al. 2003; Der Hoven and Quade 2002; Van Breemen et al. 1983). It is therefore challenging to determine the formation mechanisms of rhizoliths, which form in diff erent sedimentary and diagenetic environments(Sun et al. 2019b).

Rhizoliths in the deserts of northwest China have been widely studied (Li et al. 2015a, b, 2017; Sun et al. 2019a, b)and have been observed either at the soil surface or sub-surface. The surf icial rhizoliths are most commonly eroded out of dune soil and weathered at the soil surface, where they are subjected to late or epi-diagenesis like wind erosion,dissolution, radiation, and physical thermal expansion and cold contraction. This results in the fragmentation of rhizoliths at the soil surface (Li et al. 2015a, b, 2017; Yang 2000;Chen et al. 2004; Gao et al. 1993), modifying the original characteristics of the rhizoliths associated with their formation underground. Here, we term “in-siturhizoliths (IR)”the rhizoliths preserved vertically and deeply underground and which were not subject to any weathering and epidiagenetic process, in contrast with the reworked and weathered rhizoliths (EWR). To date, pristine rhizoliths from deserts have only been investigated at one site (Tengeri Desert,NW China) (Sun et al. 2020), where they were observed to be preserved horizontally within shallow subsurface dune soils. The pristine rhizoliths from the Badain Jaran deserts diff er from those of the Tengeri Desert, as they are positioned vertically within the deep Holocene soils stratigraphy(Fig. 1 a, b). These two types of pristine rhizoliths – horizontal- or vertical- – might have diff erent mechanisms of formation, potentially related to contrasting environmental settings. In the present study, the pristine rhizoliths from the Badain Jaran Desert were characterized using f ield observations as well as microscopic and isotopic techniques to elucidate their formation mechanisms.

2 Geographical setting

The Badain Jaran desert, at an altitude of 1,200–1,700 m a.s.l., is characterized by the co-existence of mega-dunes and lakes (Ma and Edmunds 2006; Yang et al. 2011; Jiao et al. 2015; Shao et al. 2015; Zhang et al. 2017). The heights of the mega-dunes are usually 200–300 m and can reach over 500 m in its southeastern part (Dong et al., 2013). The annual precipitation is 40–120 mm and mainly occurs in summer. The annual potential evaporation is over 2500 mm with mean annual air temperatures from 9.8 to 10.2 °C. The mean annual surface wind speeds range from 2.8 to 4.6 m·s-1, and the dominant wind directions are northwest and west (Dong et al. 2004; Zhang et al. 2017). The precipitations which are ＜ 5 mm account for approximately 90% of all rain/snow events in the desert and evaporate from the mega-dune surface in 1–3 days (data from the recent f ifty years; Ma et al. 2014 and 2017). Our investigations over the last years showed that the vegetation in the southeastern margin of the desert is mainly constituted ofArtemisia,Psammochloa, andPhragmitesspecies. They are mainly present in the lowermost portions of the landscapes, i.e.swales and depressions between the dunes in the deserts of NW China.

Fig. 1 Locations of the study area and the sampling site

The dune soils contain primary and secondary rock and mineral debris from outside of the deserts and the local weathered bedrocks by wind transportation and mixture,and the average carbonate content of the soils is 6.47 %(Sun et al. 2019a). Previous works also conf irmed that carbonate minerals are present in low amounts in the silicate desert sands based on XRD analysis (Gates et al. 2008; Bai 2011). The calcium carbonate content is between 0.5 % and 2.5 % by weight in the dune sands of the Badain Jaran Desert (Li and Yang 2004). The average carbonate content of the Badain Jaran desert is 4.64% (Wang et al. 2004), while it can be as high as 23% in the calcareous layers in these deserts (Yang et al. 2003). It is suffi ciently large to impact the radiocarbon ages of dune sediment (Yang et al. 2010) and rhizolith carbonate cement. Soil carbonates formed through cycling generations with a mixture of lithogenic and pedogenic pre-carbonate, as well as new-formed carbonate from weathered silicates (Monger 2014; Monger et al. 2015).

3 Materials and methods

In the vast area of the Badain Jaran desert, EWR are commonly scattered due to erosion and weathering on dune soil surfaces. In contrast, IRs were only found at one place in the southeastern peripheral margin of the desert. Therefore,rhizoliths were investigated and collected from f ive sites in an area of the margin (Fig. 2 a, ～39°28′N, ～102°26′E). EWR occurred in the swales between small dunes (Fig. 2 b). A pit of about 50 cm depth was carefully excavated around IR(Fig. 3 a) until it was entirely exposed (Fig. 3 b, c, d). The petrographic and mineralogic analyses of one IR specimen were carried out in the laboratory of the Northwest Petroleum Institute, Lanzhou, China. The specimens were f irst impregnated with resin, cut into transverse sections, and polished to 4.8 cm × 2.8 cm slices. A Zeiss Scope A1 microscope was used to examine the mineral crystal morphology under transmitted plain light. The technique of Dickson( 1966) was applied to stain the calcareous cement of the rhizoliths, with alizarin red-S allowing to discriminate between calcite and dolomite, and then a mixture of potassium ferricyanide and alizarin red-S indicating the presence of ferroan or non-ferroan calcite and/or ferroan dolomite. The luminescence pattern of the grains and cement forming the studied rhizoliths was investigated by cathodoluminescence(CL) with a CITL CL8200 MK5 instrument. A fragment of IR rhizolith was gold coated, and then analyzed using an FEI Quanta 450 scanning electron microscope (SEM) coupled with an energy dispersive X-ray spectrometer (EDX)to examine the ultra-microscale features and the chemical composition of the cement.

The radiocarbon dating and isotope analyses were performed at the Guangzhou Institute of Geochemistry, Chinese Academy of Sciences. The EWR tubes with root hairs (Fig. 3 e), the IRs with root relicts, and the root hairs(Fig. 3 e) were radiocarbon dated. They were treated using the standard acid-alkali-acid method (Shen et al. 2010).The samples were f irstly smashed into small pieces, and then cleaned and washed with ultrapure water (MilliQ,Millipore) in an ultrasonic bath before being treated with phosphoric acid (～ 100%). The CO2gas resulting from this treatment was purif ied in a vacuum line and f inally sealed into a reaction tube for graphitization (Xu et al. 2007). The graphite was measured with a 0.5 MeV compact accelerator mass spectrometry (National Electrostatics Corporation,0.5MV 1.5SDH-1 AMS). The precision of the measurement is close to 2‰. The data were corrected for isotopic fractionation using online AMS δ13C values. All the radiocarbon ages were calibrated into cal BP using Intcal20.14C(Reimer et al. 2020) (Table 1). Root samples and other parts of fresh plant tissue were treated in the traditional way (acidalkali-acid treatment). Samples were dried in a freeze dryer.Subsamples were split up into quartz tubes in the presence of copper oxide and silver wire. The quartz tubes were pumped to 1.0 × 10-3torr in a vacuum line and then sealed using a torch. The sealed tubes were placed into an oven at 900 °C for organic carbon oxidation. The δ13C and δ18O values of the carbonate cement of f ive rhizoliths and threeArtemisiaplant specimens were also determined separately using an isotope ratio mass spectrometer (Thermo Finigan Delta Plus XL; Table 2).

Fig. 2 Geomorphology of sampling sites of the studied rhizoliths

Soil paleotemperatures derived from rhizolith carbonates were calculated using the following equation developed by Dworkin et al. ( 2005) based on the database of the Holocene carbonates (Cerling and Quade 1993):

where δ18O and MAAT represent the oxygen isotopic composition of the carbonate (‰), Pee Dee Belemnite(PDB) and the mean annual air temperature (°C), respectively. This equation is suitable for carbonates with δ18O values comprised between -13‰ and -1‰ (PDB). The values outside this range ref lect post-formation processes (e.g.diagenesis, recrystallization) or evaporative enrichment(Dworkin et al. 2005).

4 Results

4.1 Field observations on rhizoliths and host soils

The eroded and weathered rhizoliths (EWRs) were found in the swales and depressions between dune chains (Fig. 2 a).The whitish broken rhizotubes (Fig. 3 a) were exhumed by wind erosion and looked like frost or snow from distance(Fig. 2 b).Psammochloaherbs (Fig. 2 a) andArtemisia spp.shrubs (Fig. 3 a) occurred sparsely in the swales and depressions (Fig. 2 a).

Only the vertical rhizoliths preserved in sand soil were expected to be pristine,in-siturhizoliths (IR). Only f ive IRs were found in an area of about 20 × 20 m2in a swale where EWRs were also present (Fig. 3 a). One vertical rhizolith was extracted for detailed examination in the present study.Wind striations were visible across the pit wall (Fig. 3 b) and on the rhizoliths (Fig. 3 c). The soil was dry at the surface(Fig. 3 a, b, c) and moist below 10–15 cm depth (Fig. 3 b) and can be classif ied as Leptic Regosols (FAO, 2014).

The underground part of the IR was much thicker and more fragile than the aboveground part (Fig. 3 c). The central part of the IR corresponds to a hollow channel f illed with root residuals and sand (Fig. 3 d). The halo around this central part is whitish to yellow-brown and becomes greybeige with distance from the center (Fig. 3 d).The EWRs appeared as small whitish broken tubes densely and randomly scattered at the soil surface (Fig. 3 a).Hair-like root relicts were found within one EWR (Fig. 3 e).The hollow channel of the EWR was also made of brown root remains mixed with sand (Fig. 3 f).

4.2 Petrography

Microscopically, the IR consists of poorly cemented grains,mainly of quartz, but also feldspars, and carbonates that are stained brownish-red with Alizarin red-S (Fig. 4 a).The lithogenic carbonates or primary carbonates particles originate from windblown materials from around the desert bedrocks (Sun et al. 2019a) and appear as patches or lumps(Fig. 4 a). The grains are entrapped by micritic calcareous cement with no sign of post-diagenetic alteration (Fig. 4 a).

Cathodoluminescence analyses (Fig. 4 b) conf irmed the nature of the clastic grains of the rhizoliths, with quartz grains in dark brown-blue (Mavris et al. 2012; Omer 2015),feldspars in blue (alkali feldspar) or green (plagioclase)(Scholonek and Augustsson 2016), and lithogenic carbonates in orange ( Omer et al. 2014; Kolchugin et al. 2016).The carbonate cement is recognizable as bright orange dots scattered among clastic grains (Fig. 4 b). The cement consists of low magnesium calcite (Fig. 4 c, d).

Fig. 3 Field characteristics of the rhizoliths

4.3 Radiocarbon ages

The results of radiocarbon dating of the IR (Fig. 3 d) and EWR (Fig. 3 e) specimens are shown in Table 1. The age of the organic relicts in the central hollow channel of the IR is2053 ± 16 years cal BP. The root hair relicts within the EWR are modern in age (2009 C.E.). The age of the EWR cement is 5232 ± 19 years cal BP and the IR cement is 8628 ± 32 years cal BP.

Table 1 AMS radiocarbon ages of rhizoliths in the Badain Jaran Desert

4.4 Stable carbon and oxygen isotopes

The δ13C and δ18O values of the carbonate cement and the carbonate particles of rhizoliths andArtemisia spp. tissues were shown in Table 2. The mean δ13C (VPDB) and δ18O(VPDB) values of the EWR cement are -2.6‰ and -6.8‰,respectively. Similar values were observed for the IR cement, with δ13C and δ18O values of -2.5‰ and -6.2‰.Artemisiatissues are characterized by δ13C and δ18O values of -26.7‰ to -28.5‰ and - 27.4‰ to -32.1‰, respectively.

5 Discussion

5.1 Ages of rhizoliths, stable isotopic composition and environmental implications

Modern analogs of initial biomineralization around plant roots are rare (Alonso-Zarza 2018). However, some root relicts were radiocarbon dated to Holocene in modern loess (Gocke et al. 2011), aeolian sand dunes (Cramer and Hawkins 2009), coastal/inland sands (Joseph & Thrivikramaji, 2005), and bay dunes (Rao and Thamban 1997 ).Apparent carbonate radiocarbon ages (Table 1) indicate that the rhizolith formation in the Badain Jaran Desert occurred during the Holocene (8–2 ka ago).

The δ13C values of theArtemisiatissues (-26.74‰ to-28‰, Table 2) conf irm the C3nature ofArtemisia(Liu et al. 2017) because the C3photosynthesis pathway occurs in most shrubs, herbs, and (cool-season) grasses with δ13C values ranging between ～ -25‰ and ～ -32‰ (Vogel 1993;Basum et al. 2015), with a mean value of -27‰ (Smith and White 2004; Basum et al. 2015). Generally, C3plants in arid environments present slightly higher δ13C values than those from temperate regions (Cerling and Quade 1993;Zamanian et al., 2021). When the f luxes of CO2derived from organic matter decomposition in soils are higher than the rates of pedogenic carbonate formation, the δ13C values of secondary carbonates are mainly controlled by those of soil CO2and they can be used as an indicator of the local vegetation cover (Cerling 1984; Cerling and Quade 1993;Quade 2014), after adjustment for isotopic fractionation and diff usion. Rhizoliths form in isotopic equilibrium with root-derived C (Gocke et al. 2011) and their δ13C values are generally enriched by about 14‰ to 16‰ related to the root-derived C (Cerling 1984; Cerling and Quade 1993; Zamanian et al. 2016 a). In the present study, the difference between the δ13C values ofArtemisiatissues (ca.-27‰) and rhizoliths (about -2‰ to -3‰) (Table 2) is much larger than 14-16 ‰. Logically and theoretically,there are three possible explanations for this diff erence: (1)rhizoliths originate from other plant species thanArtemisia,with a C4metabolic pathway; (2) rhizoliths were subjectedto diagenetic processes, such as wind erosion, dissolution,radiation, physical thermal expansion, and cold contraction on the soil surface; (3) there was a large contribution of lithogenic carbonates due to relatively fast formation of rhizoliths in the Badain Jaran desert.

Table 2 δ 13 C and δ 18 O stable isotope values of rhizoliths cement and Artemisia tissues in the Badain Jaran Desert

Fig. 4 Microscopic features of the studied rhizoliths

The f irst hypothesis can be rejected, as f ield observations clearly showed the close association between rhizoliths andArtemisiaplants (Figs. 2 a and 3 a). Similarly, the second one can be excluded, as microscopic analyses revealed no diagenetic changes in the rhizoliths carbonates (Fig. 4 a, b, c).Therefore, the third hypothesis is the most likely and is also supported by a large age discrepancy between root remains and the rhizolith cement (Table 1). A mixture of secondary/pedogenic and lithogenic carbonates in rhizolith cement and also diff usion of atmospheric CO2in soils with low rates of respiration can lead to enrichment in δ13C values. Amundson et al. ( 1989) found that a low density of vegetation cover was the main reason for the entrance of atmospheric CO2(with δ13C values ～ -6.5‰ for pre-industrial CO2) into soils and enrichment of δ13C values in the pedogenic carbonates of the Mojave Desert, California. The enrichment of δ13C values of secondary carbonates was also observed in the deserts of central Iran, which was related to a decline in the vegetation cover and inf lux of atmospheric CO2over time (Bayat et al. 2018). Therefore, the lithogenic carbonate residues inside the rhizoliths of the desert can explain the enrichment of δ13C values. Furthermore, sparse vegetation density in the Badain Jaran desert (Dong et al. 2004; Zhang et al. 2017), low respiration rates, and a sandy texture of the soil (Figs. 2 b and 3 b) lead to inf iltration of atmospheric CO2and can further enhance this enrichment.

The δ18O value of secondary carbonates was shown to be related to the δ18O value of local meteoric water which is controlled by air temperature and evaporative conditions (Cerling 1984; Cerling and Quade 1993; Quade 2014; Zamanian et al., 2021). As14C dating of carbonates shows that rhizoliths were formed during the Holocene(Table 1), we applied an empirical model for paleotemperature reconstruction based on Holocene carbonates (Dworkin et al. 2005). The reconstructed temperatures from rhizoliths of the Badain Jaran desert (Table 2) are ca. 2 to 3 °C higher than the mean annual temperature recorded in the region(Dong et al. 2004; Zhang et al. 2017). However, the analyzed rhizoliths consist of a large contribution (more than 90%) of lithogenic carbonates from our pre-studies (Sun et al. 2019a) that challenges the reliability of any paleo-environmental reconstructions based on rhizoliths for Badain Jaran desert, even though the determined paleo-environmental information looks reasonable.

5.2 Mechanism of rhizolith formation in dune soil of Badain Jaran desert

There are two potential sources of Ca2+for calcif ication of the initial intact rhizoliths in the Badain Jaran desert: (i) dissolution of the lithogenic carbonates of dune clastic particles transported by the wind; (ii) weathering of other minerals like feldspars (e.g. plagioclases) (Fig. 4 b; Gudbrandsson et al. 2014). The release of organic acids via organisms during root decomposition increases the dissolution rate and solubility of silicate minerals like feldspars (Welch and Ullman 1996), which releases calcium cations. The mechanisms of feldspar dissolution and weathering are subjects of discussion (Yuan et al. 2019). Nevertheless, the weathering rates of feldspars were shown to be proportional to their exposed surface areas (Holdren and Speyer 1985) and the population of surface complexes with H+, OH-or ligands (Bloom and Nater 1991).

When the root dies, it can still provide a conduit for downward and upward percolation of solutions and air.This leads to the dissolution of carbonates and minerals weathering as well as carbonates re-precipitation within or around the decaying roots, which forms a petrif ied root tubule (Sarjeant 1975; Nascimento et al. 2019). Calcite can precipitate along the root hollow channel due to the steep chemical gradient occurring around the dead root. The formation of calcite in and around roots is favored by the episodic drying phases of the soil, which often occurs in deserts due to strong solar radiation and wind (Cohen 1982;Kraus and Hasiotis 2006). These episodic drying events prevent the complete dissolution of lithogenic carbonates, but the entrapment in the rhizolith cement (i.e. the pedogenic carbonate). This formation procedure points to the fact that rhizoliths form relatively fast, over a few drying events,and even after root death. Furthermore, the sandy structure of the soil facilitated cementation with a relatively small amount of precipitated pedogenic carbonates. The specif ic fast formation and incomplete recrystallization of carbonates (Zamanian et al. 2016b and c) in equilibrium with respired CO2make the studied rhizoliths in Badain Jaran Desert inappropriate for paleo-environmental reconstructions and radiocarbon dating. This is in contrast with many recently published works that suggest the suitability of rhizoliths for paleo-environmental reconstruction (Gocke et al.2010, 2011, 2014; Bojanowski et al. 2015; Li et al. 2015a,b; Nascimento et al. 2019).

Soil moisture is a key factor in the rhizolith formation process. Rhizolith formation around decaying roots needs a semi-closed, weak redox environment (Fig. 5). Indeed,both shallow dry, and deeper waterlogged soils are not favorable environments for root decomposition and associated CO2production, as an open and dry environment prevents root decay, while closed reducing environments lead to organics fermentation and CH4production rather than HCO3-(Sun et al. 2019a). The homogenous porous sandy soil of the Badain Jaran Desert (Fig. 3 a, b) and the harsh climatic condition of the region provide a favorable condition in terms of soil moisture, i.e., episodic wetting and drying cycles. Soil moisture is further controlled by landscape geomorphology. The desert vegetation is solely distributed at the foot of dunes (Fig. 2 a, b) on the beds of swales and depressions, where rainwater accumulates and underground pore water migrates from higher landscape positions. In the depressions, the Ca2+concentration also increases (Bai 2011), which further favors rhizolith formation.

In the Badain Jaran Desert, IRs are progressively eroded out of dune soils and become thinner due to wind and weathering, forming the EWRs. Some features of the former roots, such as root hairs within the hollow channels of the EWR (Fig. 3 d, e), can still be preserved because of extreme aridity at the surface and also protection off ered by the rhizolith carbonate cement. Nevertheless, f ine roots are rarely found within the EWRs, as they are, mostly, lost by wind erosion (Fig. 5).

Fig. 5 Conceptual diagram showing formation mechanism of the IR in dune sandy soils. After root death, rhizolith formation starts in semiclosed redox conditions at about 30–50 cm depth below the surface. In these conditions, water is available for chemical reactions (i.e., mineral dissolution and carbonate crystallization) leading to root decomposition, CO 2 release, and mineral dissolution. At the same time, diff erent types of organic acids are produced, inducing a decrease in pH around the root. A steep chemical gradient is generated between the root and surrounding soil. Ca 2+ ions are released via the dissolution of minerals such as feldspar and lithogenic carbonates. They may later crystallize as calcite around the decomposed roots. At the same time, the channels left after root decomposition provide a conduit for water transportation to the surface, leading then to evaporation and further carbonate precipitation in the channel and/or around the decomposed root relicts.The water transportation to the surface (thick yellow arrows) is due to extreme weather conditions in deserts such as strong winds and solar irradiation. The rhizoliths can be preserved in underground soil. Then,wind may remove the surface soil layer and expose the rhizoliths on the surface. The latter can become thinner and thinner by chemical and mechanical erosion and weathering on the soil surface, and at last broken, lying on the surface

6 Conclusions

The f ield observations and laboratory analyses of rhizoliths from dune soils of Badain Jaran Desert constrain the main factors inf luencing their formation:

a. Vegetation type. This is the f irst key factor for rhizolith formation in this area as the rhizoliths are only formed around the modern dead roots ofArtemisiarather than roots of any other species such asPsammochloa villosa(Trin.) Bor, andPhragmites Australis.

b. Soil moisture. Rhizoliths are only present in moist,deeper soil horizons of the dune swales. These sandy soils with suitable water content provide weak oxidation-reduction conditions for gas exchange as well as water migration and evaporation leading to continuous root decomposition, solutes transport, and calcif ication.

c. Topographic position. Soil pore-water migrates from upper parts of dunes to swales and depressions,where rainfall and dune condensation water is collected. Water is permanently available forArtemisiagrowth in soils of swales and depressions despite the extreme evaporative condition of the region.

d. High aridity and evaporation potential. The high aridity and evaporation due to strong solar radiation and strong winds favor the condensation and crystallization of dissolved carbonates around roots and the formation of rhizoliths.Thein-siturhizoliths are shown to be formed over a very short time after root death, leading to the incorporation of“old” carbon from lithogenic carbonates (as much as 90% of rhizolith carbonates) into the rhizolith structure (Sun et al.2019a), as evidenced by the14C date discrepancies between rhizolith carbonates and organic relicts of the roots. Isotopic analyses also highlighted the eff ects of the low-density vegetation cover and high rates of evaporation in deserts on13C and18O isotopic enrichment of rhizoliths. Altogether,these results show that rhizoliths are the ref lection of the soil complex and specif ic environmental and climatic conditions in which they are formed.

Acknowledgements We highly appreciate the constructive comments from Professor Zhuolun Li. We give our thanks to Professor Ping Ding for radiocarbon dating and isotopic work, Pu Wang and Shiming Zhang for their petrographic and mineralogical work; Wenhui Xue and Wentao Pei for their assistance during f ield exploration. The authors also acknowledge Professor Eric Verrecchia, Professor Xin Wang, Professor Yinglong Chen, Professor Zhiqiang Lei, and Dr. Keyu Fa for their comments. The research was supported by the Natural Science Foundation of China (41561046), Chinese-German Centre (Sino-German Mobility M-0069) and the German Research Foundation (DFG)(ZA 1068/4-1).

Funding Open Access funding enabled and organized by Projekt DEAL.

Availability of data and material All data generated and analyzed during this study are included in this published article.

Declarations

Conf lict of interest On behalf of all authors, the corresponding author states that there is no conf lict of interest.

Open AccessThis article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing,adaptation, distribution and reproduction in any medium or format,as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.