Next Article in Journal
Interactions of Polyamines and Phytohormones in Plant Response to Abiotic Stress
Previous Article in Journal
SNP Diversity and Genetic Structure of “Rogosija”, an Old Western Balkan Durum Wheat Collection
Previous Article in Special Issue
The Relationship and Influencing Factors between Endangered Plant Tetraena mongolica and Soil Microorganisms in West Ordos Desert Ecosystem, Northern China
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Plants
Volume 12
Issue 5
10.3390/plants12051158
plants-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Different Responses of Growing Season Ecosystem CO2 Fluxes to Rain Addition in a Desert Ecosystem

by
Xiaotian Xu
1,2,
Bo Wu
1,3,*,
Fang Bao
1,3,
Ying Gao
1,3,
Xinle Li
4,5,
Yanli Cao
1,3,
Qi Lu
1,3,
Junliang Gao
4,5,
Zhiming Xin
4,5 and
Minghu Liu
4,5
1
Institute of Desertification Studies, Chinese Academy of Forestry, Beijing 100091, China
2
Institute of Forestry and Pomology, Beijing Academy of Agriculture and Forestry Sciences, Beijing 100093, China
3
Key Laboratory of Desert Ecosystem and Global Change, State Administration of Forestry and Grassland, Beijing 100091, China
4
The Experimental Center of Desert Forestry of the Chinese Academy of Forestry, Bayannur 015200, China
5
Dengkou Desert Ecosystem Research Station of Inner Mongolia, Bayannur 015200, China
*
Author to whom correspondence should be addressed.
Plants 2023, 12(5), 1158; https://doi.org/10.3390/plants12051158
Submission received: 17 January 2023 / Revised: 21 February 2023 / Accepted: 22 February 2023 / Published: 3 March 2023
(This article belongs to the Special Issue Sand Vegetation and Restoration)
Browse Figures
Versions Notes

Abstract

:
Desert ecosystem CO2 exchange may play an important role in global carbon cycling. However, it is still not clear how the CO2 fluxes of shrub-dominated desert ecosystems respond to precipitation changes. We performed a 10-year long-term rain addition experiment in a Nitraria tangutorum desert ecosystem in northwestern China. In the growing seasons of 2016 and 2017, with three rain addition treatments (natural precipitation +0%, +50%, and +100% of annual average precipitation), gross ecosystem photosynthesis (GEP), ecosystem respiration (ER), and net ecosystem CO2 exchange (NEE) were measured. The GEP responded nonlinearly and the ER linearly to rain addition. The NEE presented a nonlinear response along the rain addition gradient, with a saturation threshold by rain addition between +50% and +100%. The growing season mean NEE ranged from −2.25 to −5.38 μmol CO2 m−2 s−1, showing net CO2 uptake effect, with significant enhancement (more negative) under the rain addition treatments. Although natural rainfall fluctuated greatly in the growing seasons of 2016 and 2017, reaching 134.8% and 44.0% of the historical average, the NEE values remained stable. Our findings highlight that growing season CO2 sequestration in desert ecosystems will increase against the background of increasing precipitation levels. The different responses of GEP and ER of desert ecosystems under changing precipitation regimes should be considered in global change models.

1. Introduction

Arid regions dominated by desert ecosystems occupy 35% of the Earth’s surface and store 15% of the organic carbon of the world [1,2]. However, the role of desert ecosystems in the global carbon cycle has long been underestimated due to long-term drought stress, sparse vegetation, low soil nutrient levels, and low CO2 fluxes. Recently, temperate shrublands have been found to dominate the interannual variability of the global terrestrial carbon cycle [3], and desert ecosystems have been reported to play a considerable role as a CO2 sink [4,5,6]. However, due to great uncertainties and a lack of research, whether desert ecosystems are CO2 sources or sinks remains controversial [7]. Global circulation models predict that the total amount of precipitation in mid and high latitudes will show a rising trend in the near future, with an increasing frequency of extreme precipitation events [8,9]. Future climate scenarios predict changes in the precipitation regimes in the desert regions of northwestern China, following an increasing trend [6,8,9,10]. For example, based on the RCP8.5 scenario, increases in mean annual precipitation of up to 25%, 50% [11], and even more than 100% [12] are expected in certain desert regions at the middle and end of the 21st century compared with the end of the 20th century.
In arid regions, soil water is the main factor limiting the diversity and vigor of perennial plant species and shaping the structure, function, and diversity of ecosystems [13]. Changes in precipitation amounts and patterns (e.g., magnitude, interval, and seasonality) will affect soil water availability and hence influence the physiological and ecological processes of plants growing in xeric ecosystems [14,15,16,17]. Increased precipitation stimulates ecosystem CO2 uptake through photosynthesis and CO2 release via respiration, subsequently influencing the net ecosystem CO2 exchange (NEE), which reflects the overall CO2 sink function of ecosystems. Previous studies have shown that increased precipitation can enhance CO2 fluxes and productivity more significantly than decreased precipitation across various terrestrial ecosystems [18]. However, different vegetation types may respond differently to similar precipitation changes [19,20]. Although numerous studies have been conducted to evaluate the effects of precipitation variations on ecosystem-level CO2 fluxes through changing water availability [17], most of them (more than 75%) [21] primarily focused on forest and grassland ecosystems [22,23,24].
Accordingly, we still lack data sets on how water controls carbon cycling, especially regarding the magnitudes of CO2 exchange fluxes in desert shrublands, and the results of existing studies remain largely various. For example, a study in the Chihuahuan Desert found that increased precipitation enhanced soil CO2 release and vegetation CO2 uptake, consequently enhancing NEE (becoming more negative in convention for the sign) [25]. Results from a desert grassland also showed that a large precipitation pulse enhanced NEE of both low- and high-cover areas [26]. However, in a study by Snyder et al., summertime irrigation shifted the soil water potential but did not increase either CO2 uptake or plant growth in shrubs in the Great Basin and the Mojave Desert [27]. The results from a rangeland experiment in Arizona indicated that NEE during precipitation pulses was strongly controlled by CO2 efflux and driven by species effects and soil microclimate [28,29]. In the Sonoran Desert, it has been found that the precipitation pulse effects on soil respiration were enlarged in fine-textured vegetated region, especially when the antecedent conditions were dry [30]. The results from an experiment in a temperate desert indicate that whether increased precipitation enhanced or dampened NEE depends on the responses of precipitation timing and ephemeral plant growing season coupling [31]. Given that many factors, including plant community type, soil type, seasonality of precipitation, and climate, may play a huge role, the explanation of the inconsistent conclusions of the above-mentioned studies were limited. As the review by Beier et al. mentioned, an important way to better clarify the responding mechanisms of desert CO2 exchange under changing precipitation regimes is to conduct multilevel treatments considering the current and future climate scenarios [21]. From the findings of the previous studies, it has been well accepted generally that ANPP and CO2 uptake increased with the increase in water availability, leveling off when extreme wet conditions are reached [32]. Both field experiments and model predictions found the response of CO2 fluxes rates showed similar concave down curves [33,34,35,36,37], with a saturation pattern. However, multiple treatment precipitation amounts involving comparisons between increases and decreases were only considered in a few studies [38,39], which is also the case for the responses of soil respiration to precipitation with gradients in grassland [14,40,41] and desert ecosystems [42]. In this sense, focusing on the integrates of CO2 fluxes along precipitation gradients, which are involved in extreme precipitation conditions that refer to future climate scenarios, will be helpful for clarifying the patterns of CO2 fluxes responding to precipitation change.
In this study, we conducted a manipulated experiment along a rain addition gradient to evaluate the CO2 fluxes responses of a Nitraria tangutorum-dominated desert ecosystem in Inner Mongolia, northwest China, investigating the responses of the ecosystem CO2 exchange processes. We attempted to answer the following questions: How do ecosystem CO2 fluxes (i.e., GEP, ER, and NEE) respond along a rain addition gradient? Does additional rain increase or decrease the net ecosystem CO2 uptake in a desert ecosystem? Given that desert ecosystems are often characterized with vegetation constrain and resource limitation, we tested the following hypotheses: (a) With increasing annual growing season rainfall, GEP, ER, and NEE will all increase nonlinearly, with different saturation thresholds. (b) The photosynthesis process (GEP) will benefit more than the respiration process (ER), leading to an enhanced net ecosystem CO2 exchange (NEE).

2. Results

2.1. Rain Addition Treatments and SWC Changes

The natural rainfall in the growing seasons (May to September) in 2016 and 2017 was 167.2 and 54.5 mm, respectively. Compared to the historical mean growing season value of 124 mm from 1961 to 2006 (Figure S1), they represented a wet year (2016) and a dry year (2017). On June 13, 2016, a hailstorm occurred in the study area, with a rainfall of 44.8 mm, causing a large number of leaves in the nabkhas to fall off (more than 90%; Figure 1a). The growing season rainfall was increased by the rain addition from 167.2 and 54.5 mm in 2016 and 2017 to 238.9 mm and 127 mm for A + 50% treatment, respectively, and for 311.4 and 199.5 mm for the A + 100% treatment.
As shown in Figure 1b, the SWC of the nabkhas in the growing season was low. If there was no large rainfall event, the SWC was mostly around 1.5%. Rain addition resulted in a rapid increase in SWC, reaching a peak within a short time, followed by a gradual decrease. In the case of heavy natural rainfall, the effect of rain addition and natural rainfall may be superposed. In addition, both A + 50% and A + 100% treatments resulted in an increase in SWC at a depth of 10 cm by 40.8% and 66.7% compared to the control plots, respectively. While at the depth of 20 cm, only A + 100% resulted in an increase in SWC. At the depth of 50 cm, the two treatments had little effect on SWC in the two growing seasons; this was also the case for natural rainfall events larger than 40 mm.

2.2. Response of Ecosystem CO2 Fluxes to Rain Addition

During the 9-day-a-month measurement period (from 1 day before to 7 days after the rain addition treatment), significant pulse responses of CO2 fluxes to rain addition were observed (Figure 2 and Figure S4). In most cases, GEP increased more persistently and steadily after the treatment, while ER increased rapidly after rain addition and then decreased to a level close to that prior to rain addition, showing a response peak (Figure S4), and NEE showed a similar pattern with GEP. In addition, natural rainfall events around the 9-day windows have certain influences on the response pulses of CO2 fluxes. As shown in Figure S4, in May and August in 2016 and May in 2017, there were natural rainfalls >2 mm within the measurement periods, where CO2 fluxes 7 days after the rain additions were enhanced. In June, July, and September in 2016, there were natural rainfalls >2 mm before rain additions, which may have enhanced the CO2 fluxes 1 day before the rain additions. Among them, only the two heaviest rainfall events in June and August in 2016 may have induced sound impacts (Figure 1 and Figure S4).
The pulse responses of CO2 fluxes showed seasonal variation, with different response amplitudes in different months (p < 0.05, Table S1). Correspondingly, all CO2 fluxes also changed greatly in different months (p < 0.001, Table 1). In 2016, the magnitude of all components of CO2 fluxes in July and August were about 60% higher (more negative for NEE) than those in other months, and the lowest CO2 fluxes were observed in May except June. In June 2016, the GEP was close to zero due to the loss of most of the leaves, caused by hail. In 2017, the peaks of the magnitude of CO2 fluxes were advanced to June and July, due to a 22 mm rainfall event in June. In September of these two growing seasons, when the leaves of N. tangutorum began to wither and fall, the magnitude of all CO2 fluxes rapidly decreased (Figure 2).
Overall, rain addition significantly affected all components of the ecosystem’s CO2 flux (p < 0.001, Table 1). The mean values of GEP, ER, and NEE across the two growing seasons for the control were 5.20, 2.95, and −2.25 μmol CO2 m−2 s−1, respectively. For the A + 50% treatment, they were 10.48, 5.01, and −5.38 μmol CO2 m−2 s−1, and for the A + 100% treatment, they were 10.27, 6.57, and −3.70 μmol CO2 m−2 s−1. Compared with the control, the magnitude of GEP, ER, and NEE values for the A + 50% treatment increased by 102%, 73% and 139%, respectively (more negative for NEE), and those for the A + 100% treatment increased by 98%, 123% and 64%, respectively (Figure 2). Across all treatments, the mean growing season values of GEP and ER significant decreased from 9.95 and 5.89 μmol CO2 m−2 s−1 in 2016 to 7.85 and 3.85 μmol CO2 m−2 s−1 in 2017, respectively, showing significant interannual variation (p < 0.05) (Figure 2). In contrast, NEE values in the growing season in 2016 (−3.55 μmol CO2 m−2 s−1) and 2017 (−4.00 μmol CO2 m−2 s−1) were not significantly different, without interannual variation (Table 2).

2.3. Response of CO2 Fluxes to Total Rainfall Amount (TRA)

In 2016 and 2017, the mean CO2 fluxes during the growing season were fitted with the total rainfall amount (TRA, the sum of added rainfall and natural rainfall) as the x-axis (Figure 3). For GEP (Figure 3a), the promotion effect of the A + 50% treatment was stronger than that of the A + 100% treatment and exhibited a quadratic nonlinear response pattern, while the relationship was significant only in 2017 (p < 0.01). For ER (Figure 3b), A + 100% had a stronger promotion effect than A + 50% and showed a significant linear response pattern (p < 0.05 in 2016 and p < 0.001 in 2017). The fitting with NEE was similar to GEP, exhibiting a nonlinear but not statistically significant pattern (p > 0.05) (Figure 3c).

2.4. Influences of Environmental Factors on CO2 Fluxes

Based on the results of the regression analysis, among the environmental factors which may affect CO2 fluxes, the air temperature in the chamber and soil moisture at a depth of 20 cm had better relationship with the CO2 fluxes than the soil temperature and soil moisture at 10 cm depth. Air temperature was positively correlated with ER and GEP and negatively correlated with NEE (p < 0.05). The relationship between ER and soil temperature at a depth of 10 cm was poor (p > 0.05), while the relationships of GEP and NEE with soil temperature at 10 cm were better (p < 0.05, except for GEP in 2016). Soil moisture at 10 cm had a weak influence on CO2 fluxes and was only significantly correlated with ER in 2017 (p < 0.001) and NEE in 2016 (p < 0.05), while soil moisture at 20 cm had stronger influence on CO2 fluxes (p < 0.05 except for GEP in 2016 and NEE in 2017) (Table 3). Further multiple regression analyses indicated that these changes in CO2 fluxes were primarily driven by the changes in soil moisture at 20 cm, especially in 2017. Soil moisture at 10 cm only exhibited stronger effects than 20 cm in 2017, while temperatures were generally ineffective except for ER. The effects of soil nutrient factors (total nitrogen, STN, and organic carbon, SOC) on CO2 fluxes were also evaluated, with generally no significant relationships (p > 0.05) except for STN and ER in 2017 (p < 0.05) (Figure 4). Further partial relationship analysis showed that, after controlling the meteorological drivers (air temperature, soil temperature, SWC), the general nonsignificant pattern did not change, with the largest partial correlation coefficients of 0.551 (p = 0.06) being between STN and ER.

3. Discussion

3.1. Differential Response Patterns of CO2 Fluxes

In our study, all components of the CO2 flux increased in magnitude with rain addition, but their response patterns were different. In line with our hypothesis, the enhancement of GEP showed a nonlinear response, and its possible saturation threshold was between A + 50% and A + 100%. However, ER showed different response patterns with GEP, with a linear response, which did not support our hypothesis. As a results, NEE showed a nonlinear response (Figure 3).
As a shrub of the genus Nitraria in the family Tribulus (Zygophyllaceae), N. tangutorum is endemic to China and mainly distributed in Inner Mongolia, Gansu, Shaanxi, Xinjiang, and other provinces in northwestern China; it also grows in arid Gobi deserts and at the edges of sandy deserts [43]. Because of its strong tolerance and adaptability to drought, salinity, cold, and wind, it has become the main dominant species of desert vegetation and a good sand-fixing plant in the arid region of northwestern China [44,45,46]. In the N. tangutorum distribution area, long-term annual average precipitation is about 80–300 mm, which means that the total precipitation corresponding to the A + 100% treatment in this study reached the upper limit of precipitation in the natural distribution area of N. tangutorum. N. tangutorum may adapt to low amounts of water available, with the photosynthesis capacity parameters saturated at a high precipitation level [47]. Therefore, it is not surprising that GEP and NEE tended to show saturation thresholds at a TRA value close to 300 mm (Figure 3). Compared to ER, the patterns of GEP and NEE against the total rainfall amounts were substantially scattered, and although the plants received more precipitation in 2016 than 2017, the GEP almost unchanged in each plot (Figure 3). This indicated that the extra water did not play a significant role. Besides the possible limited photosynthesis capacity, another key factor regulating the interannual variations of CO2 fluxes was the difference of meteorological conditions. The year of 2016 was much wetter than 2017, with two heavy natural rainfall events >40 mm, largely promoting the CO2 fluxes. Nevertheless, one of these two rainfall events was a hailstorm, causing great damage to plant leaves, which limited the overall seasonal GEP. The extra water from these heavy rainfall events may not be fully utilized by the plants, reducing the water use efficiency, and they weakened the relationship between GEP and environmental factors (Table 3). In addition, these two rainfall events were both around the 9-day windows of the CO2 flux measurement, leading to uncertainties of our results. The conclusion of our study may still need more validation in the future works.
The mechanism underlying the differences in photosynthesis (GEP) and respiration (ER) response patterns can be explained by their different response modes to rainfall. Generally, as an assimilation process, the response of photosynthesis to rainfall lags behind that of respiration because it requires water to infiltrate into the soil depth available to the plant and to last long enough to stimulate the plant to absorb water [48,49]. In our study region, the root system of N. tangutorum can distribute into deep soil of about 3 m below the aboveground shoots, with about 1–2 m in the sandy nabkhas and 1 m in the clay plain below the sand dunes [50]. However, respiration responds rapidly after rainfall, because the soil microbial community in the top layer will respond rapidly to even small rainfall events [48,51]. Therefore, large rainfall can directly translate to an increase in respiration, while the increase in photosynthesis can be steadier and more colimited by other factors. Previous studies showed that in semiarid grasslands, the nonlinear response of respiration is mainly due to substrate limitation [22,36,41], and increased soil water availability may translate water limitation into nutrient limitation [52,53,54]. However, our data in this study and our previous publications [55,56] suggested that on the N. tangutorum nabkhas with poor sandy soil, rain addition may not bring substrate and nutrient limitation. First, SOC and STN generally increased among the rain addition gradient (p < 0.05, Table 4), showing a strong positive relationship between SWC and STN (R2 = 0.577, p < 0.05). Second, rain addition promoted the growth of N. tangutorum [55] as well as of annual plants (Salsola collina and Agriophyllum squarrosum) in this field, which can be indicated by the increase in total vegetation cover (Table 4). Third, our previous study in this field had also demonstrated that rain addition increased microbial biomass carbon [56], indicating that the promotion of vegetation growth on the nabkhas may enhance microbially mediated litter decomposition, thus providing more nutrients for the soil. In addition, a former study in an N. tangutorum desert ecosystem confirmed that soil respiration increased linearly along a rain addition gradient [42]. Therefore, for ER, the linear response can be attributed to sufficient amounts of substrate and the rapid promotion of ER by the instantaneous reaction of soil microbes. These results indicate that under the TRA condition in this study, respiration did not reach the saturation threshold.

3.2. Rain Addition Enhanced NEE of the Desert Ecosystem

From the perspective of the growing season NEE of the N. tangutorum nabkhas when plants are leafed out and growing, the mean values ranged from −2.25 to −5.38 μmol CO2 m−2 s−1, showing net CO2 uptake effect. After rain addition, CO2 fluxes were enhanced, and NEE were promoted by up to more than 200% in magnitude (Figure 2). Under completely different natural rainfall conditions, based on the observation of the three CO2 flux components for 2 years, the growing season NEE was stable and robust.
The enhancement effect of NEE was the result of the response of GEP and ER to rain addition. In our study, the magnitude of all CO2 flux components were promoted by rain addition, and the response amplitude of GEP was greater (from 5.20 to 10.48 μmol CO2 m−2 s−1) than that of ER (2.95 from 6.57 μmol CO2 m−2 s−1). This is in agreement with studies of the N. tangutorum desert ecosystem in Minqin, Gansu Province, and with other studies conducted on the same rain addition experimental platform used in this study [42,47,55,57]. In these studies, leaf photosynthesis and soil respiration rates were enhanced, and plant shoot growth also increased significantly with rain addition. Our results indicate that NEE was more closely related to GEP (R2 = 0.76) than ER (R2 = 0.25), similar to previous studies in other regions [34]. Previous studies suggested that plant photosynthesis is predominately limited by available water, nutrients, and leaf photosynthesis capacity [36,58], with greater sensitivity than respiration, which is often determined by belowground activities of roots and microbes colimited by soil substrate and nutrient [59,60]. Therefore, under changing water availability, GEP can show greater variability and had greater impact on NEE. In arid ecosystems, such as grasslands and deserts, similar findings have been widely reported, where increased precipitation stimulated higher CO2 uptake than release [26,61,62,63]. These results highlight that an increase in water supply would favor CO2 sequestration by promoting aboveground plant activity (see also [58]), which could back up modeling studies reporting that more extreme rainfall events may benefit arid ecosystems by providing more available water [16].
Although natural precipitation varied greatly in 2016 and 2017, in both growing seasons, rain addition generally enhanced the NEE in magnitude. In June 2016, an unexpected hailstorm resulted in a large loss of leaves and, consequently, a sharp decline in photosynthesis. However, from the perspective of the entire growing season, NEE was not significantly affected. This is mainly because the hailstorm brought a large amount of rain, but also because the negative impact of hail on plants [64] and the positive impact of the sufficient water supply in the early growing season on the net CO2 uptake of the ecosystem [65,66] cancel each other out. The year 2017 was a dry year, but limited natural precipitation and low SWC did not significantly reduce the growing season NEE, even in the control plots without rain addition (Figure 1, Table 1). Generally, arid regions show a high interannual variability of precipitation. Based on the historical records from 1961–2006, the average variability of growing season precipitation in our study site was 31.87%, ranging from −69.68 to +112.82% (37.6 to 263.9 mm) (Figure S1). Therefore, the vegetation community in this region may show intrinsic insensitiveness to short-term precipitation fluctuations and other environmental variations.

3.3. Implications for the Carbon Cycling of Desert Ecosystems

Our findings based on the long-term rain addition experiment showed that the N. tangutorum desert ecosystem acts as net CO2 uptake during the growing season, with a significant increase under the A + 50% and A + 100% treatments. Although natural precipitation in 2016 and 2017 varied greatly, rainfall in the growing seasons of 2016 and 2017 was 134.8 and 44.0%, respectively, of the historical mean value of 124 mm (1961–2006). Although the hailstorm in June 2016 caused the loss of more than 90% of the leaves of N. tangutorum, the NEE remained stable. Therefore, against the background of increasing precipitation levels, CO2 sequestration in desert ecosystems will increase. However, our study primarily focused on the fluxes during the 5-month growing season, while the actual annual NEE can be codetermined by how rain addition treatments influence the respiration (both for the plant and soil) during the 7-month nongrowing season when most CO2 fluxes are toward the atmosphere, which means the whole picture of annual ecosystem CO2 exchange still need more investigations.
The NEE presented a nonlinear response along the rain addition gradient, and its possible saturation threshold was between A + 50% and A + 100%. Previous studies have indicated that both the GEP and ER of grassland ecosystems showed a nonlinear increase with increasing rainfall levels [22,36,41]. In contrast to the grassland ecosystem, GEP showed a nonlinear and ER a linear response along the rain addition gradient in the desert ecosystem (Figure 3). In this sense, the different response patterns of photosynthesis and respiration of desert ecosystems to changing precipitation regimes should be considered in a global change model.

4. Materials and Methods

4.1. Study Area

The experiment was carried out in a desert ecosystem on the northeast edge of the Ulan Buh Desert, Dengkou County, Inner Mongolia Autonomous Region, China (E106°43′, N40°24′, 1050 m a. s. l., Figure S2). The area is characterized by an arid continental climate, with a mean annual temperature of 7.8 °C and a mean annual precipitation of 145 mm, concentrated from July to September. Mean annual potential evaporation is 2327 mm, and the frost-free period is 136–205 days. The vegetation type is temperate desert shrubland with a vegetation cover of 20–35%. The community develops on nabkhas (oval sand dunes with deep-rooted shrubs) and is solely dominated by N. tangutorum. Each nabkhas was 0.8–1.8 m in height, with a short and long axis length of 5–9 m. The root system of N. tangutorum can distribute about 3 m below the aboveground shoots, with about 1–2 m in the sand dunes and 1 m in the clay plain below the sand dunes [50]. Generally, nabkhas scatter on a red clay ground to form a mosaic desert landscape and are 3–8 m away from each other. Occasionally, some species, e.g., Salsola collina and Agriophyllum squarrosum, also occur in the community, albeit at lower numbers. The soils covering the nabkhas consist of >90% sand and are thus categorized as sandy soils according to the China Soil Taxonomy System [67].

4.2. Experimental Design

According to historical records obtained from a local weather station (downloaded from the National Meteorological Information Center, http://data.cma.cn/ (accessed on 22 November, 2018)), the annual amount of precipitation has increased markedly from 1961–2006, which was more evident in the growing season (May to September, Figure S1). This trend was consistent with the future trend predicted by models for this region [8,68,69]. The rain addition experiment began in 2008 and is still ongoing based on a mean annual precipitation of 145 mm (1961–2006) at the study site. We used a completely randomized design with three treatments of rain addition by 0% (C), 50% (A + 50%), and 100% (A + 100%), and three replicates for each treatment (nine plots in total). In each plot of a circle area with a diameter of 12 m (113.04 m2), only one natural nabkha is located in the center of the plot. The nabkhas in all plots were similar regarding growing condition and size. The morphological characteristics of the nabkhas (length, height, area, and volume) as well as the plant cover are shown in Table 4. The distance between every two plots was more than 5 m to avoid possible interactions. During the growing season (May to September), rain addition was applied on the 15th of each month, and the rain amount added was 14.5 mm per month for A + 50% (72.5 mm per year) and 29 mm per month for A + 100% (145 mm per year). Water was pumped into a tank from a well near the plots and then irrigated into the plots via an irrigation system with a water-pump, water meters, and spraying arms (Figure S3a). The irrigation systems were installed on the top of the nabkha in the center of the plot. Two spraying arms can rotate freely about 0.3 m above the shrub crown on the top of nabkha to evenly spray the simulated rainwater over the treatment area. To reduce water evaporation at the time of rain addition, water addition was conducted in the early morning when the air temperature was relatively low and the air above the land surface was usually calm. The groundwater in the experimental site was below a depth of 5 m and did therefore not affect plant growth.

4.3. Ecosystem CO2 Flux Measurements

In each plot, the middle position of the southern slope of the nabkhas was selected. The vegetation at this position grows well and distributes evenly and can represent the general vegetation coverage on the nabkhas well, and it was used to measure the CO2 fluxes. In April 2016, a square stainless-steel frame (0.50 m × 0.50 m) was inserted into the soil at the selected measurement site on all nabkhas to a depth of 3 cm (Figure S3b). During installation, soil disturbance was largely avoided. In the growing seasons of 2016 and 2017, we used an infrared gas analyzer (IRGA, Li-8100, LI-COR Inc., Lincoln, NE, USA) attached to a transparent measuring chamber (0.50 m × 0.50 m × 0.90 m) to measure NEE and ecosystem respiration (ER) (Figure S3c). The frames provided a flat base between the soil surface and the chamber, and the box covered the entire vegetation within the frame. Two plastic tubes with a length of about 1 m connected the analyzer and the box to make sure that the air that was pulled through the tubing to the IRGA was recirculated back to the chamber during each measurement. During the measurements, two small electric fans installed within the box worked continuously to ensure that the CO2 was evenly mixed and the air temperature in the box did not increase too quickly. Each measurement of 2 min included 20 s of prepurge, a dead band of 30 s to achieve steady state, 10 s of postpurge, and a ventilation interval of 30 s before the next measurement. Therefore, the effective measurement time was a period of 1 min, with 60 consecutive recordings of CO2 concentrations at 1 s intervals under a relatively steady condition (change of air and soil temperature less than 0.2 °C, which were monitored immediately during the measurement using the JM 222 hand thermometer (Jinming, Tianjin, China) and the auxiliary temperature sensor of Li-8100, respectively). After measuring NEE, the chamber was lifted and vented to recover the CO2 concentration for 30 s, placed on the base frame, and covered with a cardboard box, creating a dark environment for the chamber by blocking the sunlight to measure ER. Since light (and hence photosynthesis) was eliminated, the value of the second CO2 exchange measurement represented ER.
In each month, we measured NEE and ER between 10:00 a.m. and 12:00 a.m. 1 day before, 1 day after, and 7 days after rain addition. Every three treatment plots were taken as one set, and the three sets of plots were measured in a consistent order.
After measuring NEE and ER, the gross ecosystem photosynthesis (GEP) was calculated as follows:
GEP = ER − NEE
Based on this equation, the values of both GEP and ER were positive, while for NEE, positive values represent CO2 release, and negative values represent CO2 absorption.

4.4. Soil Property Measurements

We used an EM-50 (decagon, Pullman, WA, USA) to measure the soil water content (SWC) of three typical plots (representing three treatments) at the depths of 10, 20, and 50 cm. Due to probe failure, data were lost for the depth of 20 cm of A + 100% after 10 August 2017. Since the EM-50 ECH2O EC-5 sensor is less accurate when used on sandy soils with low SWC, it was corrected based on previous studies under similar soil conditions [70].
For all nabkhas in all plots, soil nitrogen and organic carbon contents were determined, as shown in Table 4, using soil samples collected at depths of 0–20 cm. Five uniformly distributed replicates per plot were obtained using a soil core sampler (5 cm in diameter) after the end of rain addition in September 2016. All soil samples were air-dried and ground to pass a 100-mesh (mesh size of 150 μm) screen for analysis. Soil nitrogen content was determined by the Auto-Kjeldahl method (Kjektec System 1026 Distilling Unit, Hoganas, Sweden), while organic carbon was determined using the sulfuric acid and an aqueous potassium dichromate (K2Cr2O7) mixture with external heating. Measurements were conducted at the Research Center of Plant Ecology, Institute of Botany, Chinese Academy of Sciences, Beijing, China.

4.5. Statistic Analysis

Repeated measures analyses of variances (ANOVAs) were conducted to examine the interannual variability in the seasonal means of the ecosystem CO2 fluxes (NEE, ER, and GEP) when combined with the rain addition treatments. Repeated measures ANOVAs were used again to examine the treatment effects over the growing season each year. The between-subject effect was rain addition, and the within-subject effect was the month of the year and its interactions with rain addition. We also used repeated measures ANOVAs to investigate the effects of rain addition in each measurement month; the within-subject effects were the measurement dates and the interactions with the rain addition. Linear and exponential regressions, multiple regression, and partial correlation analysis were conducted between the CO2 fluxes and the environmental factors (air temperature, soil temperature, soil moisture at a depth of 10 cm and 20 cm, STN, and SOC). All analyses were performed using the statistical software SPSS 21.0 (IBM, Amonk, NY, USA).

5. Conclusions

In summary, our field experiment demonstrated that CO2 fluxes of the N. tangutorum desert ecosystem can be enhanced by the increased rainfall amounts, with the GEP responding nonlinearly and the ER linearly to rain additions. Correspondingly, the NEE presented a nonlinear response along the rain addition gradient, with enhanced CO2 uptake of −2.25 to −5.38 μmol CO2 m−2 s−1. These findings indicate that, in the forthcoming climate change scenarios, the growing season CO2 sequestration in desert ecosystems will be more active when the precipitation increases, but this trend may reach a saturation when the precipitation doubles due to the different response patterns of photosynthesis and respiration. We suggest that these influencing factors should be considered in the future global change models to improve the accuracy of prediction.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/plants12051158/s1, Table S1: Results (F values) of repeated measures ANOVA on the effects of the rain addition treatments, measurement days relative to the treatment day and their interactions on NEE (net ecosystem carbon exchange), ER (ecosystem respiration), and GEP (gross ecosystem photosynthesis) during each month in 2016 and 2017; Figure S1: Geographical location of the study area; Figure S2: Variation trend of annual precipitation and precipitation in growing season before the experiment (1960–2006); Figure S3: Photo of (a) the simulated rain addition equipment on the top of a nabkha, (b) measurement position, and (c) assimilation box and Li-8100 analyzer; Figure S4: ER (ecosystem respiration) and GEP (gross ecosystem photosynthesis) on 1 day before, 1 day after, and 7 days after the rain addition treatment in the N. tangutorum desert ecosystem in the two growing seasons.

Author Contributions

Conceptualization, B.W.; methodology, X.X. and B.W.; data curation, X.X., X.L., Y.C. and Z.X.; writing—original draft preparation, X.X. and B.W.; writing—review and editing, X.X., B.W., F.B., Y.G., J.G. and Z.X.; funding acquisition, X.X., B.W., X.L., Q.L. and M.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (grant number 31600394), the Science and Technology Foundation of the Chinese Academy of Forestry (grant number CAFYBB2007008), the State Key Research and Development Program of China (grant number 2016YFC0500806), the Funding of Basic Scientific Research Operations of the Chinese Academy of Forestry (grant number CAFYBB2020MB007), and the Central and local science and technology development fund project (Grant Nos. 2022ZY0177).

Data Availability Statement

The data that support the findings of this study are available on request from the corresponding author.

Acknowledgments

We thank Zhiqiang Wang for assisting with the data collection in the field.

Conflicts of Interest

The authors declare no competing interests.

References

  1. Lal, R. Carbon sequestration in dryland ecosystems. Environ. Manag. 2004, 33, 528–544. [Google Scholar] [CrossRef]
  2. Stone, R. Ecosystems—Have desert researchers discovered a hidden loop in the carbon cycle? Science 2008, 320, 1409–1410. [Google Scholar] [CrossRef] [PubMed]
  3. Ahlstrom, A.; Raupach, M.R.; Schurgers, G.; Smith, B.; Arneth, A.; Jung, M.; Reichstein, M.; Canadell, J.G.; Friedlingstein, P.; Jain, A.K.; et al. The dominant role of semi-arid ecosystems in the trend and variability of the land CO2 sink. Science 2015, 348, 895–899. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Jasoni, R.L.; Smith, S.D.; Arnone, J.A. Net ecosystem CO2 exchange in Mojave Desert shrublands during the eighth year of exposure to elevated CO2. Glob. Chang. Biol. 2005, 11, 749–756. [Google Scholar] [CrossRef]
  5. Wohlfahrt, G.; Fenstermaker, L.F.; Iii, J.A.A. Large annual net ecosystem CO2 uptake of a Mojave Desert ecosystem. Glob. Chang. Biol. 2008, 14, 1475–1487. [Google Scholar] [CrossRef]
  6. Xie, J.; Li, Y.; Zhai, C.; Li, C.; Lan, Z. CO2 absorption by alkaline soils and its implication to the global carbon cycle. Environ. Geol. 2009, 56, 953–961. [Google Scholar] [CrossRef]
  7. Schlesinger, W.H.; Belnap, J.; Marion, G. On carbon sequestration in desert ecosystems. Glob. Chang. Biol. 2009, 15, 1488–1490. [Google Scholar] [CrossRef]
  8. IPCC. Summary for Policymakers. In Climate Change 2013: The Physical Science Basis; Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Cambridge University Press: Cambridge, UK, 2013. [Google Scholar]
  9. Zhang, X.; Zwiers, F.W.; Hegerl, G.C.; Lambert, F.H.; Gillett, N.P.; Solomon, S.; Stott, P.A.; Nozawa, T. Detection of human influence on twentieth-century precipitation trends. Nature 2007, 448, 461–465. [Google Scholar] [CrossRef]
  10. Li, Y.; Huang, J.; Ji, M.; Ran, J. Dryland Expansion in Northern China from 1948 to 2008. Adv. Atmos. Sci. 2015, 32, 870–876. [Google Scholar] [CrossRef]
  11. Gao, Y.; Li, X.; Liu, L.; Jia, R.; Yang, H.; Li, G.; Wei, Y. Seasonal variation of carbon exchange from a revegetation area in a Chinese desert. Agric. For. Meteorol. 2012, 156, 134–142. [Google Scholar] [CrossRef]
  12. Wang, Y.; Zhou, B.; Qin, D.; Wu, J.; Gao, R.; Song, L. Changes in Mean and Extreme Temperature and Precipitation over the Arid Region of Northwestern China: Observation and Projection. Adv. Atmos. Sci. 2017, 34, 9–25. [Google Scholar] [CrossRef]
  13. Rodríguez-Iturbe, I.; Porporato, A. Ecohydrology of Water-Controlled Ecosystems: Soil Moisture and Plant Dynamics; Cambridge University Press: Cambridge, UK, 2005. [Google Scholar]
  14. Fay, P.A.; Kaufman, D.M.; Nippert, J.B.; Carlisle, J.D.; Harper, C.W. Changes in grassland ecosystem function due to extreme rainfall events: Implications for responses to climate change. Glob. Chang. Biol. 2008, 14, 1600–1608. [Google Scholar] [CrossRef]
  15. Felton, A.J.; Slette, I.J.; Smith, M.D.; Knapp, A.K. Precipitation amount and event size interact to reduce ecosystem functioning during dry years in a mesic grassland. Glob. Chang. Biol. 2020, 26, 658–668. [Google Scholar] [CrossRef] [PubMed]
  16. Knapp, A.K.; Beier, C.; Briske, D.D.; Classen, A.T.; Luo, Y.; Reichstein, M.; Smith, M.D.; Smith, S.D.; Bell, J.E.; Fay, P.A.; et al. Consequences of More Extreme Precipitation Regimes for Terrestrial Ecosystems. Bioscience 2008, 58, 811–821. [Google Scholar] [CrossRef]
  17. Weltzin, J.F.; Loik, M.E.; Schwinning, S.; Williams, D.G.; Fay, P.A.; Haddad, B.M.; Harte, J.; Huxman, T.E.; Knapp, A.K.; Lin, G.; et al. Assessing the response of terrestrial ecosystems to potential changes in precipitation. Bioscience 2003, 53, 941–952. [Google Scholar] [CrossRef]
  18. Wu, Z.; Dijkstra, P.; Koch, G.W. Responses of terrestrial ecosystems to temperature and precipitation change: A meta-analysis of experimental manipulation. Glob. Chang. Biol. 2011, 17, 927–942. [Google Scholar] [CrossRef] [Green Version]
  19. Wilcox, K.R.; von Fischer, J.C.; Muscha, J.M.; Petersen, M.K.; Knapp, A.K. Contrasting above- and belowground sensitivity of three Great Plains grasslands to altered rainfall regimes. Glob. Chang. Biol. 2015, 21, 335–344. [Google Scholar] [CrossRef] [PubMed]
  20. Zeppel, M.J.B.; Wilks, J.V.; Lewis, J.D. Impacts of extreme precipitation and seasonal changes in precipitation on plants. Biogeosciences 2014, 11, 3083–3093. [Google Scholar] [CrossRef] [Green Version]
  21. Beier, C.; Beierkuhnlein, C.; Wohlgemuth, T.; Penuelas, J.; Emmett, B.; Körner, C.; de Boeck, H.; Christensen, J.H.; Leuzinger, S.; Janssens, I.A.; et al. Precipitation manipulation experiments—Challenges and recommendations for the future. Ecol. Lett. 2012, 15, 899–911. [Google Scholar] [CrossRef]
  22. Chen, S.; Lin, G.; Huang, J.; Jenerette, G.D. Dependence of carbon sequestration on the differential responses of ecosystem photosynthesis and respiration to rain pulses in a semiarid steppe. Glob. Chang. Biol. 2009, 15, 2450–2461. [Google Scholar] [CrossRef]
  23. Luyssaert, S.; Inglima, I.; Jung M Richardson, A.D.; Papale, D.; Piao, S.; Schulze, E.D.; Wingate, L.; Matteucci, G. CO2 balance of boreal, temperate, and tropical forests derived from a global database. Glob. Chang. Biol. 2007, 13, 2509–2537. [Google Scholar] [CrossRef] [Green Version]
  24. Ru, J.; Zhou, Y.; Hui, D.; Zheng, M.; Wan, S. Shifts of growing-season precipitation peaks decrease soil respiration in a semiarid grassland. Glob. Chang. Biol. 2018, 24, 1001–1011. [Google Scholar] [CrossRef] [PubMed]
  25. Patrick, L.; Cable, J.; Potts, D.; Ignace, D.; Barron-Gafford, G.; Griffith, A.; Alpert, H.; Van Gestel, N.; Robertson, T.; Huxman, T.E.; et al. Effects of an increase in summer precipitation on leaf, soil, and ecosystem fluxes of CO2 and H2O in a sotol grassland in Big Bend National Park, Texas. Oecologia 2007, 151, 704–718. [Google Scholar] [CrossRef] [PubMed]
  26. Hamerlynck, E.P.; Scott, R.L.; Stone, J.J. Soil moisture and ecosystem function responses of desert grassland varying in vegetative cover to a saturating precipitation pulse. Ecohydrology 2012, 5, 297–305. [Google Scholar] [CrossRef]
  27. Snyder, K.A.; Donovan, L.A.; James, J.J.; Tiller, R.L.; Richards, J.H. Extensive summer water pulses do not necessarily lead to canopy growth of Great Basin and northern Mojave Desert shrubs. Oecologia 2004, 141, 325–334. [Google Scholar] [CrossRef]
  28. Huxman, T.E.; Cable, J.M.; Ignace, D.D.; Eilts, J.A.; English, N.B.; Weltzin, J.; Williams, D.G. Response of net ecosystem gas exchange to a simulated precipitation pulse in a semi-arid grassland: The role of native versus non-native grasses and soil texture. Oecologia 2004, 141, 295–305. [Google Scholar] [CrossRef]
  29. Potts, D.L.; Huxman, T.E.; Cable, J.M.; English, N.B.; Ignace, D.D.; Eilts, J.A.; Mason, M.J.; Weltzin, J.F.; Williams, D.G. Antecedent moisture and seasonal precipitation influence the response of canopy-scale carbon and water exchange to rainfall pulses in a semi-arid grassland. New Phytol. 2006, 170, 849–860. [Google Scholar] [CrossRef]
  30. Cable, J.M.; Ogle, K.; Williams, D.G.; Weltzin, J.F.; Huxman, T.E. Soil texture drives responses of soil respiration to precipitation pulses in the Sonoran Desert: Implications for climate change. Ecosystems 2008, 11, 961–979. [Google Scholar] [CrossRef]
  31. Huang, G.; Li, Y.; Padilla, F.M. Ephemeral plants mediate responses of ecosystem carbon exchange to increased precipitation in a temperate desert. Agric. For. Meteorol. 2015, 201, 141–152. [Google Scholar] [CrossRef]
  32. Knapp, A.K.; Ciais, P.; Smith, M.D. Reconciling inconsistencies in precipitation–productivity relationships: Implications for climate change. New Phytol. 2017, 214, 41–47. [Google Scholar] [CrossRef] [Green Version]
  33. Burkett, V.R.; Wilcox, D.A.; Stottlemyer, R.; Barrow, W.; Fagre, D.; Baron, J.; Price, J.; Nielsen, J.L.; Allen, C.D.; Peterson, D.L.; et al. Nonlinear dynamics in ecosystem response to climatic change: Case studies and policy implications. Ecol. Complex. 2005, 2, 357–394. [Google Scholar] [CrossRef] [Green Version]
  34. Flombaum, P.; Yahdjian, L.; Sala, O.E. Global-change drivers of ecosystem functioning modulated by natural variability and saturating responses. Glob. Chang. Biol. 2017, 23, 503–511. [Google Scholar] [CrossRef] [PubMed]
  35. Li, L.; Kang, X.; Biederman, J.A.; Wang, W.; Qian, R.; Zheng, Z.; Zhang, B.; Ran, Q.; Xu, C.; Liu, W.; et al. Nonlinear carbon cycling responses to precipitation variability in a semiarid grassland. Sci. Total Environ. 2021, 781, 147062. [Google Scholar] [CrossRef]
  36. Zhang, F.; Quan, Q.; Ma, F.; Tian, D.; Zhou, Q.; Niu, S. Differential responses of ecosystem carbon flux components to experimental precipitation gradient in an alpine meadow. Funct. Ecol. 2019, 33, 889–900. [Google Scholar] [CrossRef]
  37. Zhou, X.; Weng, E.; Luo, Y. Modeling patterns of nonlinearity in ecosystem responses to temperature, CO2, and precipitation changes. Ecol. Appl. 2008, 18, 453–466. [Google Scholar] [CrossRef] [Green Version]
  38. Grime, J.P.; Fridley, J.D.; Askew, A.P.; Thompson, K.; Hodgson, J.G.; Bennett, C.R. Long-term resistance to simulated climate change in an infertile grassland. Proc. Natl. Acad. Sci. USA 2008, 105, 10028–10032. [Google Scholar] [CrossRef] [Green Version]
  39. Yahdjian, L.; Sala, O.E. Vegetation structure constrains primary production response to water availability in the Patagonian steppe. Ecology 2006, 87, 952–962. [Google Scholar] [CrossRef]
  40. Harper, C.W.; Blair, J.M.; Fay, P.A.; Knapp, A.K.; Carlisle, J.D. Increased rainfall variability and reduced rainfall amount decreases soil CO2 flux in a grassland ecosystem. Glob. Chang. Biol. 2005, 11, 322–334. [Google Scholar] [CrossRef]
  41. Miao, Y.; Han, H.; Du, Y.; Zhang, Q.; Jiang, L.; Hui, D.; Wan, S. Nonlinear responses of soil respiration to precipitation changes in a semiarid temperate steppe. Sci. Rep. 2017, 7, 45782. [Google Scholar] [CrossRef] [Green Version]
  42. Song, W.; Chen, S.; Wu, B.; Zhu, Y.; Zhou, Y.; Li, Y.; Cao, Y.; Lu, Q.; Lin, G. Vegetation cover and rain timing co-regulate the responses of soil CO2 efflux to rain increase in an arid desert ecosystem. Soil Biol. Biochem. 2012, 49, 114–123. [Google Scholar] [CrossRef]
  43. Pan, X.; Shen, G.; Chen, P. A preliminary research on taxonomy and systematics of genus Nitraria. Acta Bot. Yunanica 1999, 21, 287–295. (In Chinese) [Google Scholar]
  44. Li, Q.; Xu, J.; Li, H.; Wang, S.; Yan, X.; Xin, Z.; Jiang, Z.; Wang, L.; Jia, Z. Effects of aspect on clonal reproduction and biomass allocation of layering modules of Nitraria tangutorum in nebkha dunes. PLoS ONE 2013, 8, e79927. [Google Scholar] [CrossRef] [PubMed]
  45. Lv, J. Brief on exploitation and utility of Nitraria in Chaidamu. Sci. Technol. Qinghai Agric. For. 2005, 4, 38–40. (In Chinese) [Google Scholar]
  46. Ren, J.; Tao, L. A numerical taxonomy of the genus Nitraria from Gansu provence, China. Acta Bot. Boreali-Occident. Sin. 2003, 23, 572–576. (In Chinese) [Google Scholar]
  47. He, J.; Bao, F.; Wu, B.; Zhang, J.; Jia, Z. Photosynthetic acclimation of a desert species Nitraria tangutorum to long term rain addition in north-western China. Int. J. Agric. Biol. 2019, 22, 1050–1058. [Google Scholar]
  48. Huxman, T.E.; Snyder, K.A.; Tissue, D.; Leffler, A.J.; Ogle, K.; Pockman, W.T.; Sandquist, D.R.; Potts, D.L.; Schwinning, S. Precipitation pulses and carbon fluxes in semiarid and arid ecosystems. Oecologia 2004, 141, 254–268. [Google Scholar] [CrossRef]
  49. Reynolds, J.F.; Kemp, P.R.; Ogle, K.; Fernandez, R.J. Modifying the “pulse reserve” paradigm for deserts of North America: Precipitation pulses, soil water, and plant responses. Oecologia 2004, 141, 194–210. [Google Scholar] [CrossRef]
  50. Zhang, J.; Gu, L.; Bao, F.; Cao, Y.; Hao, Y.; He, J.; Meng, P. Nitrogen control of 13C enrichment in heterotrophic organs relative to leaves in a landscape-building desert plant species. Biogeosciences 2015, 12, 15–27. [Google Scholar] [CrossRef] [Green Version]
  51. Austin, A.T.; Yahdjian, L.; Stark, J.M.; Belnap, J.; Porporato, A.; Norton, U.; Ravetta, D.A.; Schaeffer, S.M. Water pulses and biogeochemical cycles in arid and semiarid ecosystems. Oecologia 2004, 141, 221–235. [Google Scholar] [CrossRef]
  52. Ren, H.; Xu, Z.; Isbell, F.; Huang, J.; Han, X.; Wan, S.; Chen, S.; Wang, R.; Zeng, D.; Jiang, Y.; et al. Exacerbated nitrogen limitation ends transient stimulation of grassland productivity by increased precipitation. Ecol. Monogr. 2017, 87, 457–469. [Google Scholar] [CrossRef]
  53. Hooper, D.U.; Johnson, L. Nitrogen limitation in dryland ecosystems: Responses to geographical and temporal variation in precipitation. Biogeochemistry 1999, 46, 247–293. [Google Scholar] [CrossRef]
  54. Maestre, F.T.; Reynolds, J.F. Amount or pattern? Grassland responses to the heterogeneity and availability of two key resources. Ecology 2007, 88, 501–511. [Google Scholar] [CrossRef] [PubMed]
  55. Zhu, Y.; Jia, Z.; Wu, B.; Lu, Q.; Yao, B. The role of increased precipitation in promoting branch and leaf growth of Nitraria tangutorum. For. Res. 2012, 25, 626–631. (In Chinese) [Google Scholar]
  56. Gao, Y.; Xu, X.; Ding, J.; Bao, F.; De Costa, Y.G.; Zhuang, W.; Wu, B. The responses to long-term water addition of soil bacterial, archaeal and fungal communities in a desert ecosystem. Microorganisms 2021, 9, 981. [Google Scholar] [CrossRef]
  57. Zhang, J.; Lu, Q.; Wu, B.; Zhu, Y.; Li, Y. Responses of Nitraria tangutorum branch and leaf growth to simulated rainfall. For. Res. 2012, 25, 130–137. (In Chinese) [Google Scholar]
  58. Niu, S.; Wu, M.; Han, Y.; Xia, J.; Li, L.; Wan, S. Water-mediated responses of ecosystem carbon fluxes to climatic change in a temperate steppe. New Phytol. 2008, 177, 209–219. [Google Scholar] [CrossRef] [PubMed]
  59. Bell, C.W.; Acosta-Martinez, V.; McIntyre, N.E.; Cox, S.; Tissue, D.T.; Zak, J.C. Linking microbial community structure and function to seasonal differences in soil moisture and temperature in a Chihuahuan desert grassland. Microb. Ecol. 2009, 58, 827–842. [Google Scholar] [CrossRef]
  60. Xiao, C.; Janssens, I.A.; Liu, P.; Zhou, Z.; Sun, O.J. Irrigation and enhanced soil carbon input effects on below-ground carbon cycling in semiarid temperate grasslands. New Phytol. 2007, 174, 835–846. [Google Scholar] [CrossRef]
  61. Hastings, S.J.; Oechel, W.C.; Muhlia-Melo, A. Diurnal, seasonal and annual variation in the net ecosystem CO2 exchange of a desert shrub community (Sarcocaulescent) in Baja California, Mexico. Glob. Chang. Biol. 2005, 11, 927–939. [Google Scholar] [CrossRef]
  62. Schwinning, S.; Sala, O.E. Hierarchy of responses to resource pulses in and and semi-arid ecosystems. Oecologia 2004, 141, 211–220. [Google Scholar] [CrossRef]
  63. Xie, J.; Zha, T.; Jia, X.; Qian, D.; Wu, B.; Zhang, Y.; Bourque, C.P.A.; Chen, J.; Sun, G.; Peltola, H. Irregular precipitation events in control of seasonal variations in CO2 exchange in a cold desert-shrub ecosystem in northwest China. J. Arid Environ. 2015, 120, 33–41. [Google Scholar] [CrossRef]
  64. Turner, M.G.; Dale, V.H. Comparing large, infrequent disturbances: What have we learned? Ecosystems 1998, 1, 493–496. [Google Scholar] [CrossRef]
  65. Niu, S.; Yang, H.; Zhang, Z.; Wu, M.; Li, L.; Han, X.; Wan, S. Non-additive effects of water and nitrogen addition on ecosystem carbon exchange in a temperate steppe. Ecosystems 2009, 12, 915–926. [Google Scholar] [CrossRef]
  66. Kwon, H.; Pendall, E.; Ewers, B.E.; Cleary, M.; Naithani, K. Spring drought regulates summer net ecosystem CO2 exchange in a sagebrush-steppe ecosystem. Agric. For. Meteorol. 2008, 148, 381–391. [Google Scholar] [CrossRef]
  67. Chen, L.; Li, F.; Di, X.; Zhao, J. Aeolian Sandy Soils in China; Science Press: Beijing, China, 1998. [Google Scholar]
  68. Sun, Y.; Ding, Y. A projection of future changes in summer precipitation and monsoon in East Asia. Sci. China Earth Sci. 2010, 53, 284–300. [Google Scholar] [CrossRef]
  69. Zhu, Y.; Wu, B.; Lu, Q. Progress in the study on response of arid zones to precipitation change. For. Res. 2012, 25, 100–106. (In Chinese) [Google Scholar]
  70. Wu, B.; Han, H.; He, J.; Zhang, J.; Cui, L.; Jia, Z.; Yang, W. Field-specific calibration and evaluation of ECH2O EC-5 sensor for sandy soils. Soil Sci. Soc. Am. J. 2014, 78, 70–78. [Google Scholar] [CrossRef]
Plants 12 01158 g001 550
Figure 1. Daily precipitation (a) and soil water content (SWC) at depths of 10, 20, and 50 cm (b) in the two growing seasons in 2016 and 2017 on the N. tangutorum-dominated nabkhas. Black arrows represent rain addition treatments. C = control, A + 50% = rain addition by 50%, and A + 100% = rain addition by 100%. SWC at depth of 20 cm in August and September of 2017 was missing due to instrument failure.
Figure 1. Daily precipitation (a) and soil water content (SWC) at depths of 10, 20, and 50 cm (b) in the two growing seasons in 2016 and 2017 on the N. tangutorum-dominated nabkhas. Black arrows represent rain addition treatments. C = control, A + 50% = rain addition by 50%, and A + 100% = rain addition by 100%. SWC at depth of 20 cm in August and September of 2017 was missing due to instrument failure.
Plants 12 01158 g001
Plants 12 01158 g002 550
Figure 2. Changes in the CO2 fluxes ((a), gross ecosystem photosynthesis, GEP; (b), ecosystem respiration, ER; and (c), net ecosystem CO2 exchange, NEE) on the N. tangutorum-dominated nabkhas in two growing seasons under rain addition treatments. Error bars represent standard errors. C = control, A + 50% = rain addition by 50%, and A + 100% = rain addition by 100%. The monthly, interannual, and overall effects can be found to Table S1, Table 1 and Table 2, respectively.
Figure 2. Changes in the CO2 fluxes ((a), gross ecosystem photosynthesis, GEP; (b), ecosystem respiration, ER; and (c), net ecosystem CO2 exchange, NEE) on the N. tangutorum-dominated nabkhas in two growing seasons under rain addition treatments. Error bars represent standard errors. C = control, A + 50% = rain addition by 50%, and A + 100% = rain addition by 100%. The monthly, interannual, and overall effects can be found to Table S1, Table 1 and Table 2, respectively.
Plants 12 01158 g002
Plants 12 01158 g003 550
Figure 3. Regression relationships between total rainfall amounts (TRAs) and mean growing season CO2 fluxes ((a), gross ecosystem photosynthesis, GEP; (b), ecosystem respiration, ER; and (c), net ecosystem CO2 exchange, NEE) in 2016 and 2017. For GEP and NEE, nonlinear regressions are shown, while for ER, linear regressions are shown according to smaller p values during regression. The regression equations with statistical significance were GEP in 2017 (y = −1.339 + 0.119x − 3.03 × 10−4x2), ER in 2016 (y = 0.027x − 0.448) and 2017 (y = 0.024x + 0.853), and NEE in 2017 (y = 2.769 − 0.107x + 3.94 × 10−4x2).
Figure 3. Regression relationships between total rainfall amounts (TRAs) and mean growing season CO2 fluxes ((a), gross ecosystem photosynthesis, GEP; (b), ecosystem respiration, ER; and (c), net ecosystem CO2 exchange, NEE) in 2016 and 2017. For GEP and NEE, nonlinear regressions are shown, while for ER, linear regressions are shown according to smaller p values during regression. The regression equations with statistical significance were GEP in 2017 (y = −1.339 + 0.119x − 3.03 × 10−4x2), ER in 2016 (y = 0.027x − 0.448) and 2017 (y = 0.024x + 0.853), and NEE in 2017 (y = 2.769 − 0.107x + 3.94 × 10−4x2).
Plants 12 01158 g003
Plants 12 01158 g004 550
Figure 4. Relationships between the CO2 fluxes ((a,d), gross ecosystem photosynthesis, GEP; (b,e), ecosystem respiration, ER; and (c,f), net ecosystem CO2 exchange, NEE) and soil nutrient factors (STN = soil total nitrogen, SOC = soil organic carbon).
Figure 4. Relationships between the CO2 fluxes ((a,d), gross ecosystem photosynthesis, GEP; (b,e), ecosystem respiration, ER; and (c,f), net ecosystem CO2 exchange, NEE) and soil nutrient factors (STN = soil total nitrogen, SOC = soil organic carbon).
Plants 12 01158 g004
Table
Table 1. Results (F values) of repeated measures ANOVA on the effects of the measurement months, rain addition treatments, and their interactions on the CO2 fluxes (a, gross ecosystem photosynthesis, GEP; b, ecosystem respiration, ER; and c, net ecosystem CO2 exchange, NEE) during the growing seasons of 2016 and 2017.
Table 1. Results (F values) of repeated measures ANOVA on the effects of the measurement months, rain addition treatments, and their interactions on the CO2 fluxes (a, gross ecosystem photosynthesis, GEP; b, ecosystem respiration, ER; and c, net ecosystem CO2 exchange, NEE) during the growing seasons of 2016 and 2017.
CO2 Flux20162017
Month (M)GEP35.974 ***13.709 ***
ER24.564 ***17.487 ***
NEE31.505 ***10.864 ***
Treatment (T)GEP4.319 *6.219 **
ER8.756 **15.993 ***
NEE3.678 *3.106
Interaction (M × T)GEP1.8611.391
ER1.0001.160
NEE2.1181.169
Notes: *, **, and *** represent significant differences at p < 0.05, p < 0.01, and p < 0.001, respectively. Bold values indicate a significant difference at p = 0.05. Sample sizes (n) of each variable were 9 and 15 for month and treatment, respectively.
Table
Table 2. Results (F values) of repeated measures ANOVA on the effects of the measurement years, rain addition treatments, and their interactions on the seasonal means of the CO2 fluxes (a, gross ecosystem photosynthesis, GEP; b, ecosystem respiration, ER; and c, net ecosystem CO2 exchange, NEE).
Table 2. Results (F values) of repeated measures ANOVA on the effects of the measurement years, rain addition treatments, and their interactions on the seasonal means of the CO2 fluxes (a, gross ecosystem photosynthesis, GEP; b, ecosystem respiration, ER; and c, net ecosystem CO2 exchange, NEE).
VariableF
Year (Y)GEP5.236 *
ER55.496 ***
NEE0.798
Treatment (T)GEP24.496 ***
ER58.812 ***
NEE12.773 ***
Interaction (Y × T)GEP1.975
ER2.157
NEE2.633 ***
Notes: * and *** represent significant differences at p < 0.05 and p < 0.001, respectively. Bold values indicate a significant difference at p = 0.05. Sample sizes (n) of each variable were 9 and 6 for year and treatment, respectively.
Table
Table 3. Relationships between the CO2 fluxes (a, gross ecosystem photosynthesis, GEP; b, ecosystem respiration, ER; and c, net ecosystem CO2 exchange, NEE) and environmental factors.
Table 3. Relationships between the CO2 fluxes (a, gross ecosystem photosynthesis, GEP; b, ecosystem respiration, ER; and c, net ecosystem CO2 exchange, NEE) and environmental factors.
YearFactorGEPERNEE
R2p ValueR2p ValueR2p Value
2016Air temperature0.0510.0170.178<0.0010.0350.03
Soil temperature0.0270.0830.0050.4290.0550.006
Soil moisture at 10 cm0.0260.0900.0020.5760.0430.016
Soil moisture at 20 cm0.0000.8480.0330.0400.129<0.001
2017Air temperature0.1000.0010.0770.0040.0470.024
Soil temperature0.0550.0070.0040.5050.0330.038
Soil moisture at 10 cm0.0040.4860.103<0.0010.0090.275
Soil moisture at 20 cm0.143<0.0010.383<0.0010.0310.059
Note: The regressions between ER, GEP, and temperature are exponential, and the others are linear. Bold values indicate a significant difference at p = 0.05. Sample sizes (n) of each variable were 135 for each year.
Table
Table 4. Characteristics of the experimental field in the N. tangutorum desert ecosystem.
Table 4. Characteristics of the experimental field in the N. tangutorum desert ecosystem.
CA + 50%A + 100%
East–west length of the nabkhas (m)5.75 ± 0.766.30 ± 0.897.00 ± 1.34
South–north length of the nabkhas (m)6.45 ± 1.036.05 ± 0.698.83 ± 1.83
Height of the nabkhas (m)1.25 ± 0.181.18 ± 0.071.40 ± 0.23
Area of the nabkhas (m2)30.51 ± 7.6729.94 ± 5.0651.64 ± 15.38
Volume of the nabkhas (m3)13.98 ± 4.9511.94 ± 2.4526.11 ± 10.20
Plant cover (%)25.00 ± 3.54 a26.25 ± 5.54 a33.75 ± 4.73 a
Relative cover of N. t.0.74 ± 0.09 a0.83 ± 0.06 a0.63 ± 0.07 a
Canopy height (cm)51.88 ± 3.73 a60.69 ± 3.69 a61.81 ± 2.95 a
Soil organic carbon (%)0.104 ± 0.012 a0.099 ± 0.006 a0.157 ± 0.024 b
Soil total nitrogen (‰)0.143 ± 0.009 a0.140 ± 0.004 a0.178 ± 0.013 b
Notes: C = control, A + 50% = rain addition by 50%, and A + 100% = rain addition by 100%. N. t. means Nitraria tangutorum. Data represent mean values ± standard errors; different letters within rows represent significant differences in the Duncan test. These data were measured in 2016.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Xu, X.; Wu, B.; Bao, F.; Gao, Y.; Li, X.; Cao, Y.; Lu, Q.; Gao, J.; Xin, Z.; Liu, M. Different Responses of Growing Season Ecosystem CO2 Fluxes to Rain Addition in a Desert Ecosystem. Plants 2023, 12, 1158. https://doi.org/10.3390/plants12051158

AMA Style

Xu X, Wu B, Bao F, Gao Y, Li X, Cao Y, Lu Q, Gao J, Xin Z, Liu M. Different Responses of Growing Season Ecosystem CO2 Fluxes to Rain Addition in a Desert Ecosystem. Plants. 2023; 12(5):1158. https://doi.org/10.3390/plants12051158

Chicago/Turabian Style

Xu, Xiaotian, Bo Wu, Fang Bao, Ying Gao, Xinle Li, Yanli Cao, Qi Lu, Junliang Gao, Zhiming Xin, and Minghu Liu. 2023. "Different Responses of Growing Season Ecosystem CO2 Fluxes to Rain Addition in a Desert Ecosystem" Plants 12, no. 5: 1158. https://doi.org/10.3390/plants12051158

APA Style

Xu, X., Wu, B., Bao, F., Gao, Y., Li, X., Cao, Y., Lu, Q., Gao, J., Xin, Z., & Liu, M. (2023). Different Responses of Growing Season Ecosystem CO2 Fluxes to Rain Addition in a Desert Ecosystem. Plants, 12(5), 1158. https://doi.org/10.3390/plants12051158

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop