Introduction

The CO2 concentration in soils is usually much higher than in the atmosphere (Amundson and Davidson 1990) and can be very dynamic (Riveros‐Iregui et al. 2007). However, the diel dynamic of the CO2 concentration in soils in relation to temperature is not yet fully understood (Philipps et al. 2011). The soil CO2 concentration results from the production and transport of CO2. CO2 is produced in soils by microbial and root respiration (Schimel et al. 2001) as well as by the dissolution of carbonates in calcareous soils (Ramnarine et al. 2012). Only at a few, exceptional sites is CO2 derived from geothermal sources (Kämpf et al. 2013). The production of CO2 by soil microorganisms depends on soil temperature and moisture (Wood et al. 1993; Moyano et al. 2013; Schimel 2018). The production of CO2 in soil by plant roots depends additionally on solar radiation (Carbone and Trumbore 2007). Traditionally, it was believed that the transport of CO2 from the soil to the atmosphere is driven only by diffusion (Amundson and Davidson 1990). However, recent studies demonstrated that also advection, i.e., movement of CO2 through movement of air, contributes to the transport of CO2 from the soil to the atmosphere (Maier et al. 2012; Rey 2015; Moya et al. 2019). Advective movement of air in soil pores and from the soil to the atmosphere can be caused by small turbulence-induced pressure fluctuations which can be induced by wind (Takles et al. 2004; Flechard et al. 2007; Maier et al. 2012; Sánchez-Cañete et al. 2013; Moya et al. 2019; Laemmel et al. 2019). In addition, thermal differences between soil and air can cause movement of air from soil pores to the atmosphere, called thermal convection, free convection, or thermal convective venting (TCV) (Rose and Guo 1995; Ganot et al. 2014; Roland et al. 2015; Rey et al. 2015; Levintal et al. 2017; Laemmel et al. 2019). TCV likely occurs mostly in (semi-)desert ecosystems (Weisbrod et al. 2009; Ganot et al. 2014; Levintal et al. 2017) because of the large differences between soil and air temperature at night (Noy-Meir 1973; Nachshon et al. 2011). However, there are so far only very few studies that explored TCV under field conditions.

TCV is caused by a higher temperature of the soil than of the air in combination with high soil permeability. The temperature difference between soil and air creates unstable conditions in the soil profile as warm, less dense gas underlies cold, denser gas (Nachshon et al. 2011; Ganot et al. 2014). If the buoyancy forces overcome the impeding viscous forces, convective gas motion within the soil profile takes place. Whether thermal convective venting (TCV) occurs in homogeneous porous media can be assessed based on the Rayleigh-Darcy number (Ra). Ra is a dimensionless number that shows the ratio of buoyant and viscous forces.

$$Ra=\frac{g\beta \Delta TkH}{\nu {\alpha }_{s}}$$
(1)

where g is the gravitational acceleration, \(\beta\) is the volumetric thermal expansion coefficient, \(\Delta T\) is the temperature difference between the bottom and the top of the soil, k is the permeability, H is the soil depth, \(\nu\) is the kinematic viscosity of the fluid, and \({\alpha }_{s}\) is the thermal diffusivity of the soil (Nield and Bejan 2013). Equation 1 shows that TCV is likely to occur in deep soils and in highly permeable soils, when the soil is warmer than the air. An increase in both \(\Delta T\) and H will increase the buoyancy forces that promote heat convection, and hence Ra. Increasing H will decrease the thermal gradient, and thus the proportion of heat transferred by conduction, which increases the buoyancy forces. Under the assumption that (i) soil is homogeneous and isotropic, (ii) Darcy’s law and the Boussinesq approximation are valid (Nield and Bejan 2013), and (iii) air is the only mobile fluid in the soil, the minimal Ra value that is needed for TCV to occur is approximately 27 (Ribando and Torrance 1976; Tan et al. 2003).

During recent years, new CO2 sensors have made it possible to measure CO2 concentrations at high frequency, which gives rise to the description of hourly and sub-hourly changes in the soil CO2 concentration (Jassal et al. 2005; Riveros‐Iregui et al. 2007; Sánchez-Cañete et al. 2013; Zhang et al. 2015; Moya et al. 2019). Based on these high frequency measurements of the soil CO2 concentration, a hysteresis loop was found when plotting the CO2 concentration as a function of the soil temperature in a high-altitudinal, semi-arid forest soil (Riveros‐Iregui et al. 2007). A very similar relationship between soil CO2 concentration and soil temperature was also observed in a temperate pine plantation (Zhang et al. 2015) and Mediterranean cropland soils (Min et al. 2020). In addition, several authors reported a hysteresis loop when plotting the CO2 efflux as a function of soil temperature (Parkin and Kaspar 2003; Gaumont-Guay et al. 2006; Ruehr et al. 2010; Jia et al. 2013; Wang et al. 2014; Song et al. 2015; Zhang et al. 2015).

The reasons for the hysteresis observed between soil temperature and soil CO2 concentration or CO2 efflux have been debated intensively. Initially, the phenomenon has been explained by changes in soil CO2 production by biota that depend on soil temperature (Riveros‐Iregui et al. 2007). Later, Ruehr et al. (2010) pointed out that the hysteretic relationship between soil temperature and CO2 observed in a forest soil indicated an unusually high temperature sensitivity of soil respiration (Q10 > 150), making it rather unlikely that changes in CO2 are caused by temperature-dependent biotic activity. Philipps et al. (2011) showed based on a model that the hysteresis can result from heat and CO2 transport processes in soil alone without changes in biological CO2 production. In the following years, several authors discussed possible biotic and the abiotic components of the hysteresis (Wang et al. 2014; Song et al. 2015; Zhang et al. 2015). However, despite the fact that temporal dynamics of CO2 are more closely related to the air temperature than to the soil temperature in some ecosystems (Parking and Kasper 2003), there is no study that analyzed and discussed changes in air temperature as the reason for the hysteresis, to our knowledge.

The objective of this study was to determine the soil CO2 concentration and explore its temporal dynamic in relation to soil and air temperature in a deeply weathered desert soil. We hypothesized that the CO2 concentration undergoes diel changes due to the strong changes in air temperature and shows a hysteretic relationship both with soil and air temperature.

Material and methods

Study site

The study site (for a photo see Supplementary Information) is located in the Coastal Cordillera of Chile in the reserve Santa Gracia close to the city La Serena, at 615 m a.s.l. (− 29.759769°, − 70.457333°). The mean annual temperature is 16.1 °C and mean annual precipitation is 87 mm per year (Oeser et al. 2018). The soil is locates on a hillslope (15° slope angle), and many coarse stones are visible on the surface (5% cover). The vegetation is sparse (∼10% cover) and dominated by Proustia cuneifolia, Cordia decandra and Adesmia spp. The parent material of the soil is granodiorite (Oeser et al. 2018). The soil is a Cambisol (according to the world reference base for soil resources (WRB)) and consists of a 30 cm deep A horizon overlying a B horizon, ranging from 30 to 80 cm depth, a B-Cw horizon ranging from 80 to 120 cm, and below 120 cm until > 200 cm a Cw horizon. The bulk density varies between 1.42 and 1.66 g cm−3 (Bernhard et al. 2018). The root density is low; 3–5 fine roots per dm2 in the A horizon, and 1–2 fine roots per dm2 in the B horizon. In the B-Cw and Cw horizon, roots are only present in fractures. The bedrock is deeply weathered, which is typical for the Coastal Cordillera of Chile (Vázquez et al. 2016), and particularly the upper meters of the soil and saprolite contain many fractures (Krone et al. 2021). More information about the soil at the study site can be found in Bernhard et al. (2018), Oeser et al. (2018), Brucker and Spohn (2019), and Krone et al. (2021).

Sampling and study design

In March 2019, we excavated a 2 m deep soil profile at the site, and the excavated material was stored in different heaps according to the depth it was taken from. Soil samples were collected from the following depths: 0–5, 5–10, 10–20 cm, and further in 20 cm intervals until 200 cm soil depth (i.e. 20–40, 40–60 cm, etc.). Subsequently, we dug 25 cm long, flat, horizontal tunnels into the profile wall using a chisel at the following depths; 15, 30, 50, 90, 135, and 185 cm. In these tunnels, the temperature and the CO2 sensors were placed. The tunnels were closed with the material that was removed when building the tunnels and, subsequently, the whole profile was closed by filling the material of the different soil depths increments into the soil pit in the right order. The temperature sensor (DS18B20) and the CO2 sensor (R30, Sense Air) of one soil depth were controlled together by a self-built Arduino clone board with flash storage and real time clock crystal (https://www.bayceer.uni-bayreuth.de/bayeos/en/sensor_net/gru/html.php?id_obj=145275). The two sensors of each soil depth and the corresponding Arduino board were powered by four lithium-ion batteries. The CO2 sensors had been calibrated before deployment in the field using gases with known CO2 concentration in closed glass jars. The temperature and the CO2 concentration were measured and logged every 60 min. Most sensors stopped measuring at the end of January 2020 when the battery charge got low. The sensors in soil depth 130 cm measured continuously over the whole year. Batteries of all sensors were changed on March 8, 2020, and subsequently, temperature as well as CO2 concentration were measured and logged until May 5, 2020.

We also used data from the weather station (Campbell Scientific) installed at the study site (for details see Übernickel et al. (2020)). The weather station measures the air temperature 200 cm above the soil using a temperature sensor (HC2S3, Rotronic) and the precipitation using a tipping bucket rain gauge (ARG 100, Environmental Measurements Limited), at an interval of 60 min.

Soil chemical and physical analyses

One part of each soil sample was dried at 15 °C. Subsequently, the gravel content was determined by sieving (< 60 mm and > 2 mm). The texture of the < 2 mm fraction was determined using the standard soil pipette texture method according to Burt (2004) after destruction of organic matter using H2O2. The pH was measured in a soil:water ratio of 1:2.5 using a pH meter (pH 340, WTW GmbH, Weilheim, Germany) with a Sentix 81 electrode. A subsample of the fine earth fraction was dried at 105 °C and ground in a ball mill for the determination of the carbonate content and the total organic C (TOC) and total N (TN) content. TOC and TN were measured using an element analyzer (Vario Max Elementar, Hanau). The absence of carbonates was confirmed using 30% HCl.

In order to assess the amount of microbial biomass, we determined the DNA content in field moist samples. The samples were passed through a sieve (2 mm). DNA was extracted from 400 mg of field moist soil using a DNA extraction kit (FastDNA™ SPIN Kit for Soil, MP Biomedicals) with small modifications as in Spohn et al. (2016). The DNA concentration was measured with the picogreen assay (Sandaa et al. 1998) using a kit (Quant-iT™ PicoGreen® dsDNA Reagent, Life Technologies).

Data analysis and calculations

For the data analysis, the first six weeks after installation of the sensors were not considered because there might have been artifacts in the CO2 concentration due to the settling of the soil in the soil pit. For the following period from May 1, 2019 until May 5, 2020, we calculated monthly means of temperature and CO2 concentration across all measurements of one month (n = 720 or 744, and 696 in February), daily means across all 24 h of the day (n = 24), and hourly means for specific months across all days of the month (n = 30 or 31 and 29 in February). All data analyses were conducted in R (R Core Team 2018).

We calculated Ra according to Eq. 1, using the following values for the constants; g = 9.806 m s−2, \(\beta\) = 3.41 × 10–3 °C−1, \(\nu\) = 1.53 × 10–5 m2 s−1. We computed \({\alpha }_{s}\) following Ganot et al. (2014) as

$${\alpha }_{s}={\kappa }_{s}{\kappa }_{f}^{-1}{\alpha }_{f}$$
(2)

where \({\kappa }_{s}\) and \({\kappa }_{f}\) are the thermal conductivities of the soil and the fluid (i.e., the air), and \({\alpha }_{f}\) is the thermal diffusivity of the fluid, using the following values: \({\kappa }_{f}\) = 2.58 × 10–2 W m−1 °C−1, \({\alpha }_{f}\) = 2.16 × 10–5 m2 s−1, \({\kappa }_{s}\) = 0.27 W m−1 °C−1 (Ganot et al. 2014). We considered ΔT as the difference between the soil temperature in a given soil depth (specified in the Results section) and the air temperature above the soil, following Levintal et al. (2017). We used different estimates of the permeability (k) based on data about the permeability for soil, sand as well as gravel and sand mixtures (Ganot et al. 2014; Côté et al. 2011; Bear, 2013), as specified in the Results section.

Results

Soil physical and chemical properties

The soil had a high content of fragments of weathered saprolite > 2 mm (called gravel in the following) and a sandy texture (Table 1). In the 100–200 cm depth, the gravel content ranged around 50% of the soil mass. The gravel content and the sand content increased with increasing soil depth, while the silt and clay content decreased with increasing depth (Table 1). The pH decreased from pH 8.3 in the upper 5 cm to pH 6.1 in 180–200 cm depth (Table 1). The soil was completely carbonate free. Both TOC and TN concentrations were highest in the upper 5 cm and decreased strongly with increasing soil depth. The TOC concentration decreased by 80% from the upper 5 cm to 40–60 cm depth. The DNA concentration was much higher in the upper 60 cm than in the lower part of the soil profile. The molar TOC:TN ratio ranged between 9.7 and 11.1 and showed no trend with soil depth. In contrast to the TOC concentration, the DNA concentration was highest from 10 to 20 cm depth (Table 1).

Table 1 Soil properties, including gravel content (mineral particles > 2 mm), soil texture (of the fraction < 2 mm), pH determined in water, total organic carbon (TOC), total nitrogen (TN), and soil DNA content, b.d.l. stands for below detection limit

Annual variation of temperature and CO2 concentrations

The mean air temperature ranged between a monthly mean of 9.2 °C in July and 20.2 °C in December. There were a few small precipitation events between May 30 and mid-July (Fig. 1). During these events, the precipitation rate did not exceed 3–4 mm per hour. The first of these events occurred on May 29.

Fig. 1
figure 1

Air temperature and precipitation at the study site from May 2019 until May 2020

The mean monthly soil temperature in 15 cm depth was 13.5 °C in July and 27.1 °C in December. The variability in the daily means of soil temperature decreased with increasing soil depth (Fig. 2a). Between May and August the daily mean soil temperature was higher in the subsoil than in the topsoil and increased with increasing depth. In contrast, from late October to April the temperature declined with increasing depth (Fig. 2a).

Fig. 2
figure 2

Soil temperature (a) and soil CO2 concentration (b) in six soil depths. Black lines show daily means and colored shades show daily standard deviations (n = 24)

The CO2 concentration increased with depth from 15 to 185 cm soil depth throughout the year (Fig. 2b). The mean of the CO2 concentration for the period from May 1, 2019 until May 5, 2020 was 463 ppm in 15 cm depth, 533 ppm in 30 cm depth, 712 ppm in 50 cm depth, and 696 ppm in 90 cm depth, 1032 ppm in 130 cm depth, and 1578 ppm in 185 cm depth. The mean CO2 concentrations of the six soil depths were positively correlated with soil depth (R2 = 0.93, p < 0.05). The variability in the daily means of the CO2 concentration decreased with increasing soil depth (Fig. 2b).

Diel variation of soil CO2

The soil CO2 concentration in the topsoil changed considerably over the day (Fig. 3). In July, the mean CO2 concentration in 15 cm depth, calculated across all days of the month, increased by 95 ppm between 10:00 h and 18:00 h. In December, the CO2 concentration in 15 cm depth changed more strongly than in July over the course of the day, namely by 211 ppm, from 418 ppm at 9:00 h to 629 ppm at 17:00 h (Fig. 3a). In December, in 30 cm depth, the maximum CO2 concentration was shifted by some hours in comparison to 15 cm soil depth, and was observed at 23:00 h (Fig. 3b). In 50 cm soil depth, the maximum CO2 concentration was time-shifted in comparison to 30 cm soil depth, and was observed at 2:00 h (Fig. 3c).

Fig. 3
figure 3

CO2 concentrations in 15 cm (a), 30 cm (b), and 50 cm (c) depth in July and December as a function of time of the day. Symbols indicate means, and arrows depict standard deviations calculated across all days of the two months in 2019 (n = 24)

Diel variation of temperature and Ra

Similar to the CO2 concentration in the topsoil, the air temperature in the upper 15 cm varied strongly over the course of the day during all seasons. In July, the mean air temperature, calculated across all days of the month, ranged from 3.6 °C at 9:00 h to 17.8 °C at 17:00 h. In December, the mean air temperature ranged from 11.0 °C at 7:00 h to 30.4 °C at 16:00 h (Fig. 3a). In contrast, the soil temperature changed little during the day, even in the uppermost 15 cm (Fig. 4b). In July, the mean soil temperature at 15 cm depth calculated across all days of the month, ranged from 12.5 °C at 13:00 h to 14.4 °C at 23:00 h. In December, the mean temperature at 15 cm soil depth ranged from 25.5 °C at 12:00 h to 28.7 °C at 23:00 h (Fig. 4b). As already shown in Fig. 2b for the daily means, in winter (July), the soil temperature was larger in 185 cm than in 15 cm depth, whereas in summer (December) the soil temperature was inverted compared to winter during all hours of the day (Figs. 4b, c).

Fig. 4
figure 4

Air temperature (a), soil temperature in 15 cm (b) and 185 cm (c) soil depth as well as the difference between soil temperature and air temperature for the uppermost 15 cm (d), 50 cm (e) and 185 cm (f) in July and December as a function of time of the day. Symbols indicate means, and arrows depict standard deviations calculated across all days of the two months in 2019 (n = 31). Please notice that the temperature differences are calculated as the soil temperature minus the air temperature

In summer (December), the temperature difference across the 185 cm ranged between 11 and − 8 °C (Fig. 4f). The temperature difference between the soil in 15 cm depth and the air was 10–19 °C between 21:00 h and 9:00 h in December (Fig. 4d). In contrast, in July, the temperature difference between the soil in 15 cm depth and the air was only 4–6 °C between 23:00 h and 9:00 h, and the soil was colder than the air during the day (Fig. 4d). In Winter (July), the temperature difference between the soil in 185 cm depth and the air above the soil was positive throughout the day and ranged between 2 and 16 °C (Fig. 4f)

Considering the whole soil profile down to 185 cm (H = 1.85 m) in July, Ra increased during the night until 10:00 h in the morning (Fig. 5a). Assuming a permeability of 1 × 10−7 m2, Ra reached Rac in the period between 4:00 and 10:00 h. At a permeability lower than 1 × 10−7 m2, Ra remained below Rac throughout the day. In summer (December), the soil temperature decreased with increasing soil depth (Fig. 2b). Therefore, for December, we only calculated Ra for the upper 15 cm (H = 0.15 m), and found that a permeability of 1 × 10−6 m2 or higher is required to reach a Ra larger than Rac during nighttime (Fig. 5b).

Fig. 5
figure 5

Rayleigh-Darcy number as a function of time of the day for the upper 185 cm of the soil profile in July (a) and the upper 15 cm of the soil in December (b). The black horizontal line indicates the critical Rayleigh-Darcy number (Rac = 27)

Hysteresis

When plotting the CO2 concentration in 15 cm soil depth as a function of the soil temperature at 15 cm depth, a hysteresis loop was observed, both for July and December (Fig. 6). The hysteresis resulted from the fact that the CO2 concentration was at its minimum in the morning, and reached its maximum in the late afternoon, while the soil temperature was at its minimum around noon and reached its maximum at 23:00 h (Fig. 6a). Another hysteresis loop was observed when plotting the CO2 concentration in 15 cm soil depth as a function of the air temperature, for both July and December (Fig. 6b). In contrast to the soil temperature that changed little during the day, air temperature changed to a similar extent as the soil CO2 concentration (in 15 cm depth) over the course of the day (Fig. 6b).

Fig. 6
figure 6

CO2 concentration in 15 cm soil depth as a function of soil temperate (a) and air temperature (b) in December and July 2019. Shown are hourly measurements calculated across all days of the month (n = 31) for December and July. Colors indicate the hour of the day

The hysteresis loops of the soil temperature and the CO2 concentration turn counterclockwise (Fig. 6a) and the hysteresis loops of the air temperature and the CO2 concentration turn clockwise (Fig. 6b). The reason for this is that the soil CO2 concentration in the topsoil had a similar dynamic as the air temperature (Fig. 6a), whereas it had time-shifted minima and maxima compared to the soil temperature in the topsoil over the course of the day (see also Figs. 3a, 4a, b).

CO2 after precipitation event

On May 29, 4 mm precipitation occurred at the study site after there had been no precipitation event with more than 0.25 mm of precipitation for several months (Fig. 1 and Übernickel et al. (2020)). Following this precipitation event in May, increased CO2 concentrations were observed at 15, 30, and 50 cm soil depths (Fig. 7a). The CO2 concentration increased faster following the rain event in 15 cm depth than in 50 cm depth. The increase in the CO2 concentration after the rain event was largest in 15 cm depth (115 ppm). The number of days over which the mean CO2 concentration remained significantly elevated was 8 days in 15 cm depth and 9 days both in 30 and 50 cm depth (Fig. 7b). The daily oscillation of the CO2 concentration described before (see Fig. 3) was still observed with the same amplitude after the precipitation event (Fig. 7a).

Fig. 7
figure 7

Soil CO2 concentrations in the time from May 21 to June 11, covering the rain event which took place on May 29. Hourly measured CO2 concentrations in six different soil depths (a) and daily mean CO2 concentrations and daily standard deviations (n = 24) in 15, 30 and 50 cm depth (b) are shown

Discussion

The unique dataset of high frequency measurements of CO2 and temperature at different depths in a coarse-textured desert soil presented here reveals strong diel oscillation of the CO2 concentration that decreased with soil depth. The oscillation of the CO2 concentration coincided with a strong diel variation of the air temperature, whereas the soil temperature, even in the topsoil, changed very little over the course of one day.

Diel changes in CO2 concentrations

The strong diel oscillation of the soil CO2 concentration in the topsoil could be caused either by variation in the production of CO2 in the soil or variation in the efflux of CO2 from the soil. The soil CO2 concentration changed by 95 to 211 ppm over the course of the day, while the soil temperature changed only by 2 °C during the same time (Fig. 6a). Given the small change in temperature, it is unlikely that such a strong change in the CO2 concentration is caused by the temperature-dependence of the production of CO2 by microorganisms or plant roots, as pointed out earlier by Ruehr et al. (2010) with respect to similar observations. We propose that the diel variation in the CO2 concentration is caused by TCV leading to transport of CO2-rich air from the soil to the atmosphere at night when the soil is substantially warmer than the atmosphere (Fig. 4). Thermal convection has been shown to cause transport of CO2 from soil when the soil has a high porosity and is warmer than the air by Ganot et al. (2014) based on column experiments.

In order to evaluate under which conditions TCV can occur in the soil under study, we calculated the Rayleigh-Darcy number (Fig. 5). In July, the soil temperature increased with increasing soil depths, and thus the soil depth over which the temperature gradient spans (H, see Eq. 1) is large. Given the large soil depth and the relatively large temperature differences between the subsoil and the air at night (Fig. 4), Ra values above Rac were reached at night and in the morning (Fig. 5a), indicating that TCV took place. The soil has high gravel and sand contents that increase with depth (Table 1), and thus a relatively high permeability. However, estimating the mean permeability over the whole soil profile is associated with a large uncertainty, which causes a large uncertainty in the Ra estimate (Fig. 5). It also needs to be considered that the bedrock at the study soil is deeply weathered, and the upper meters of the saprolite have many fractures (Krone et al. 2021). Thus, the depth (H) over which conductive transport can take place is very likely several m, and not only 1.85 m as assumed here. If a larger H is considered, Ra > Rac is reached also with smaller permeability (k). Finally, it needs to be taken into account that the soil profile is located on a hillslope (Fig. S1, Supplementary Information). For bodies with sloping boundaries, the conventional Rayleigh criterion is not directly applicable, and at slopes, TCV is much more likely to occur than on a plane, even in moderately permeable materials (Rose and Guo 1995; Nield and Bejan 2013).

At the second and third soil depth (30 and 50 cm), the CO2 concentrations also underwent diel oscillation that was time-shifted compared to the topsoil (Fig. 3). The dynamic of the CO2 concentration in the different depths of the soil can be explained by a convection cell model, where the different soil layers are different convection cells, and CO2-rich air moves slowly from lower soil layers to upper soil layers due to thermal convection.

In summer (December), the soil temperature decreased with increasing soil depths, which does not favor TCV from the subsoil (Fig. 2a). This is likely the reason why soil CO2 concentrations were larger in December than in July (Fig. 3). Still, even in December, the soil was much warmer than the air at night in all measured depths (Fig. 4d-f). Given the temperature gradient in the soil (Fig. 2b), TCV across the whole soil profile is not likely, but TCV from the upper cm of the soil might still occur. Therefore, we calculated Ra for December only over the upper 15 cm of the soil (Fig. 5b). The small H requires high permeability ≥ 1 × 10–7 m2 for Ra to exceed Rac (Fig. 5b), hence TCV to occur. It is questionable if the studied soil has such a high permeability. While the soil texture and gravel content (Table 1) indicate that the permeability is likely lower than 1 × 10–7 m2 (Bear, 2013), it might still be that the many coarse stones in the topsoil (Fig. S1, Supplementary Information) facilitate gas transport along their surface. It might also be that small channels, created by insects or roots, increase the permeability of the soil and thus contribute to TCV of the soil. In addition, it can be speculated that numerous fractures in the upper m of the soil (Krone et al. 2021) which connect the subsoil with the atmosphere, might allow for nocturnal TCV to take place in December even against a temperature gradient in the topsoil, given that the subsoil was 10 °C warmer than the atmosphere at night (Fig. 4f). This has shown to be theoretically possible by Rose and Guo (1995). Finally, it needs to be taken into account that the soil profile is located on a hillslope (Fig. S1, Supplementary Information), where, the conventional Rayleigh criterion is not directly applicable, and TCV is much more likely to occur than on a plane, as mentioned above.

Furthermore, it needs to be considered that the oscillation of the soil CO2 concentration can potentially also be caused by the radiation-dependence of plant root respiration (Carbone and Trumbore 2007). However, it seems rather unlikely that variation in photosynthetic activity alone caused the strong diel oscillation of the soil CO2 concentration throughout the year, given the low precipitation (Fig. 1) and the sparse vegetation (see Sect. 2.1 and Fig. S1, Supplementary Information). Furthermore, radiation-dependent plant activity cannot explain why the CO2 concentration was higher in December than in July, particularly given that plant activity was likely elevated in July after the rainfall (Figs. 1).

Taken together, the diel oscillation of the CO2 concentration in the topsoil of this desert could be caused either by an oscillation of (1) the light- and soil temperature-dependent production of CO2 by biota or (2) the thermal convection of CO-rich air from the soil at night. Given the sparse vegetation, the low precipitation, and the large temperature differences between soil and air at night, the latter explanation seems more likely and it also explains the higher CO2 concentrations in December than in July.

Rainfall increases CO2 concentrations

We found the CO2 concentration increased after only a few mm of precipitation (Fig. 6). The reason for this is likely that microorganisms in the desert soil react very sensitively to water inputs. This is supported by a laboratory study with topsoil from the study site showing that even small inputs of water that increased the soil water content to only 20% of the soil’s water holding capacity, strongly increased the soil respiration (Seuss 2021). However, given that only 4 mm of rainfall caused an increase in the CO2 concentration in 50 cm depth (Fig. 6), it seems likely that the increase in the CO2 concentration was not only caused by microorganisms but at least partly also by root respiration, since the small amount of water likely did not reach a depth of 50 cm in the soil. Plants might have taken up water in the topsoil, which stimulated their activity leading to higher plant root respiration, also in a depth of 50 cm. This is supported by previous studies, showing that the soil CO2 concentration is affected by precipitation through plant root respiration in desert soils (Wang et al. 2014; Song et al. 2015).

The total amount of rainfall during the study period was very low, which is related to the mega-drought Chile is currently experiencing (Morales et al. 2020). Rainfall at the study site can usually be observed in June. Yet, in June 2019, the amount of total respiration was lower compared to the three previous years (Übernickel et al. 2020). It can be assumed that in years in which the precipitation is higher at the study site, the soil CO2 concentration is elevated due to increased biological activity.

Soil CO2 in desert soils

The CO2 concentrations measured here were relatively low throughout the year, and even at 185 cm depth did not exceed 1800 ppm. Besides low precipitation (see above), this can be attributed to the relatively low root density, low organic carbon content, and low microbial biomass (as indicated by the low DNA content), which cause low CO2 production. In addition, the coarse texture (Table 1) allows for fast emission of CO2, as discussed above. The low CO2 concentration found here is in agreement with Amundson and Davidson (1990) who reported based on a literature survey that the lowest soil CO2 concentrations across different ecosystems are typically observed in desert soils. Furthermore, Amundson et al. (1989) reported about a vegetation gradient in the eastern Mojave Desert that the lowest soil CO2 concentrations in this gradient were found at a site that had the lowest vegetation density and the lowest soil TOC content. This site had a MAP of 160 mm and the soil CO2 concentration in the upper 75 cm did not exceed 1000 ppm, which is similar to the site studied here. The low CO2 concentrations observed here are also in accordance with Terhune and Harden (1991) who studied seasonal variations of CO2 concentrations in stony, coarse-textured desert soils of southern Nevada, USA and observed that the CO2 concentration in the upper 100 cm exceeded 1000 ppm only in spring. The authors attributed the relatively small soil CO2 concentration to the coarse texture of the soils that allowed for fast transport of CO2 from the soil to the atmosphere. Furthermore, the low CO2 concentrations found here are in accordance with small CO2 concentrations in two soils located at the southern margin of the Great Basin Desert in Nevada, USA, where MAP was 151 and 180 mm, respectively, and the CO2 concentration measured over the whole year did not exceed 1000 ppm in the upper 100 cm of the soils (Oerter et al. 2018).

Throughout the year, the CO2 concentration increased with soil depth despite the strong decrease in roots, TOC, and DNA with depth (Table 1). An increase in the CO2 concentration with depth has been observed in many soils in different ecosystems (Solomon and Cerling 1987; Amundson et al. 1989; Castell and Galloway 1990; Richter and Markewitz 1995; Trumbore et al. 1995; Schwendenmann and Veldkamp 2006; Kim et al. 2017; Oerter et al. 2018). In the current study, the increase in the CO2 concentration can likely be attributed to the fact that the air density increases with increasing CO2 concentration, which leads to gravitational percolation of CO2-rich air through soil pores towards deeper soil horizons (Kowalski and Sánchez-Cañete 2010; Sánchez-Cañete et al. 2013).

Conclusions

The unique dataset of high frequency measurements of CO2 and temperature at different depths of a deep, coarse-textured desert soil presented here reveals diel oscillation of the CO2 concentration and a hysteretic relationship between the soil CO2 concentration and soil as well as air temperature, as hypothesized. The diel oscillation of the CO2 concentration and the hysteretic relationship between soil CO2 concentration and air temperature were very likely caused by thermal convection, leading to transport of CO2-rich air from the soil to the atmosphere at night. Our results have important implications as they indicate that the soil CO2 concentration can be strongly affected by the difference between soil and air temperature in soils. The results suggest that TCV might strongly contribute to gas exchange of soils, which is important for soil chemical reactions including weathering processes.