Disentangling Drought and Grazing Effects on Soil Carbon Stocks and CO2 Fluxes in a Semi-Arid African Savanna

Grasslands cover ca. 30% of the global land surface and provide critical ecosystem services. Among them, carbon storage is one of the most important. However, grasslands are increasingly threatened by drought and overgrazing which might negatively affect soil carbon stocks. Despite this threat, there is a dearth of information on how drought and grazing jointly impact soil carbon stocks and CO2 fluxes in dryland grasslands. With the aid of a large field experiment, we studied the combined effects of a 5-year extreme drought and moderate grazing on soil carbon stocks, CO2 fluxes and soil chemical properties. Extreme drought was induced by reducing ambient rainfall by 66% using large rainout shelters. We found CO2 fluxes to strongly respond to the 5-year experimental drought. Extreme drought reduced CO2 emission rates by 32% compared to ambient conditions. CO2 fluxes averaged 5.7 mg m–2min–1 under drought compared to 8.3 mg m–2 min–1 under ambient conditions. CO2 fluxes were, however, not influenced by grazing. At the end of the growth period, grazed plots under ambient rainfall had released 16.3 tons of CO2 ha–1 which was 58% higher than observed on grazed plots subjected to severe drought. Soil carbon stocks were higher under drought conditions due to slower decomposition rates. Drought resulted in increased concentrations of primary macronutrients (N, P, and K), micronutrients (Zn and Mn) and pH in the top 30 cm of the soil relative to ambient conditions. The results also showed that grazing reduced the concentration of N and P in the topsoil compared to the ungrazed plots. This study provided insights on the soil carbon storage, CO2 emission rates and nutrient dynamics in a semi-arid dryland grassland as influenced by both drought and grazing. Our study also revealed that long-term extreme drought may be favorable in terms of preserving the existing soil carbon stocks through reduced CO2 release. This finding is critical for understanding future soil carbon dynamics in dryland grasslands in the face of climate change.

Grasslands cover ca. 30% of the global land surface and provide critical ecosystem services. Among them, carbon storage is one of the most important. However, grasslands are increasingly threatened by drought and overgrazing which might negatively affect soil carbon stocks. Despite this threat, there is a dearth of information on how drought and grazing jointly impact soil carbon stocks and CO 2 fluxes in dryland grasslands. With the aid of a large field experiment, we studied the combined effects of a 5-year extreme drought and moderate grazing on soil carbon stocks, CO 2 fluxes and soil chemical properties. Extreme drought was induced by reducing ambient rainfall by 66% using large rainout shelters. We found CO 2 fluxes to strongly respond to the 5-year experimental drought. Extreme drought reduced CO 2 emission rates by 32% compared to ambient conditions. CO 2 fluxes averaged 5.7 mg m −2 min −1 under drought compared to 8.3 mg m −2 min −1 under ambient conditions. CO 2 fluxes were, however, not influenced by grazing. At the end of the growth period, grazed plots under ambient rainfall had released 16.3 tons of CO 2 ha −1 which was 58% higher than observed on grazed plots subjected to severe drought. Soil carbon stocks were higher under drought conditions due to slower decomposition rates. Drought resulted in increased concentrations of primary macronutrients (N, P, and K), micronutrients (Zn and Mn) and pH in the top 30 cm of the soil relative to ambient conditions. The results also showed that grazing reduced the concentration of N and P in the topsoil compared to the ungrazed plots. This study provided insights on the soil carbon storage, CO 2 emission rates and nutrient dynamics in a semi-arid dryland grassland as influenced by both drought and grazing. Our study also revealed that long-term extreme drought may be favorable in terms of preserving the existing soil carbon stocks through reduced CO 2 release. This finding is critical for understanding future soil carbon dynamics in dryland grasslands in the face of climate change.
Keywords: C stocks, CO 2 flux, drought, dryland grasslands, grazing INTRODUCTION Grasslands are the largest terrestrial biome, covering up to 30% of the global land surface (Suttie et al., 2005). They are also a large reservoir of organic carbon, storing 10-30% of the global soil organic carbon (Zhou et al., 2019a). Despite their low productivity, dryland grasslands play an important role in global carbon sequestration as decomposition rates are generally low due to lack of water. It is estimated that about 32% of the world's soil organic carbon is stored in dryland soils (Plaza et al., 2018;Gaitán et al., 2019). In dryland grasslands, grazing and variable rainfall are among the major determinants of grassland productivity (Bat-Oyun et al., 2016;Guuroh et al., 2018). In recent years, the frequency and intensity of drought have increased considerably in global drylands, while overgrazing continues to be a problem (IPCC, 2019). This is also the case in drylands situated in Sub-Saharan Africa, SSA. Combined effects of grazing and drought can have particularly strongly impacts on the ecosystem functions and services delivered by these ecosystems (Ruppert et al., 2015;Pfeiffer et al., 2019).
The decomposition of the large soil organic matter pools could be a major source of CO 2 emission into the atmosphere, thus providing a potentially large climate change feedback and a major source of uncertainty in climate projections (Pries et al., 2017). Due to the large size of soil carbon pools in dryland grasslands (Plaza et al., 2018) and due to the uncertainty of emission capacity of these ecosystems, it is critical to carry out precise quantifications of carbon stocks and emissions. This would not only help to obtain reliable global carbon budgets but to also develop mitigation options for low-emission through proper land use management practices (Berger et al., 2019). The quantification of grassland carbon (C) fluxes is also required to identify major control factors and to clarify whether specific grassland ecosystems sequester more atmospheric CO 2 through photosynthesis than they emit to the atmosphere through respiration (Lefèvre et al., 2017). Currently, there is limited empirical carbon storage and emission data available from natural ecosystems, particularly in SSA. This is despite reports that CO 2 emissions are the largest contributor to greenhouse gas (GHG) emissions and global warming potential in SSA natural terrestrial systems (Kim et al., 2016).
Grazing and drought both have strong effects on soil carbon storage and GHG emission in dryland grasslands, mainly indirectly through their effect on biomass production (Jones et al., 2005;Ferner et al., 2018;Gaitán et al., 2019). Among these two drivers, drought has been reported to have a greater influence on biomass production than grazing (Bat-Oyun et al., 2016). It is well documented that precipitation promotes plant growth which in turn leads to increased CO 2 assimilation (Berger et al., 2019). However, that same precipitation that promotes CO 2 uptake also stimulates its emission through increased respiration (Xu et al., 2004). Accordingly, the lack of precipitation during a drought may suppress both CO 2 assimilation and emission due to constrained plant growth and reduced soil respiration. Grazing on the other hand, can impact carbon storage and emission by altering the species composition of the vegetation, and via a reduction of aboveground biomass from livestock feeding (Chen et al., 2017). In overgrazed dryland grasslands, annual grass species tend to dominate relative to perennial species Mudongo et al., 2016). Annual grasses are known to invest more on aboveground biomass and less on belowground biomass, leading to low root respiration (Bat-Oyun et al., 2016). Chen et al. (2017) reported a reduction of between 25 and 47% in CO 2 emission due to grazing in different rangeland types.
Although grazing mainly affects soil carbon storage emission via its imprint on vegetation structure and composition, grazing may also influence these variables through its impact on soil properties (Sandhage-Hofmann et al., 2015;Linstädter et al., 2016). Though not well understood, overgrazing has the potential to directly alter soil properties and processes. For example, trampling by livestock can cause soil compaction leading to reduced soil water infiltration and, consequently, reduced moisture availability to plants and microbes (Bilotta et al., 2007;Savadogo et al., 2007;Batey, 2009). Such changes in soil properties can also affect soil organic matter build-up, because processes such as decomposition and root growth are hindered (Okach et al., 2019a). Nonetheless, the effect of grazing on carbon content has not been conclusive. For instance, two studies conducted in a typical steppe Mongolia showed opposite responses: While Liu et al. (2012) reported higher soil organic carbon  found decreases in soil organic carbon due to grazing. On the other hand, grazing may actually promote soil aeration through the breaking of soil caps. Livestock also adds dung to the soil, thus increasing soil organic carbon content (Sandhage-Hofmann et al., 2015) and promoting plant growth. In addition, grazing and drought both have the potential to alter soil chemical properties, which in turn might also influence carbon storage and emission. However, we still know very little how long-term drought and grazing jointly affect soil chemical properties (see Lucci, 2019). The only experimental study we are aware of Okach et al. (2019a) implemented a short-term drought, and mainly focused on soil physical properties such as bulk density and hydraulic properties, and on the commonly reported C and N.
In summary, there is a need not only for empirically measured C storage and emissions data in dryland grasslands, but also for a better understanding of how grazing and drought jointly affect soil organic carbon (Lucci, 2019). Accordingly, this study aims to disentangle joint effects of long-term drought and grazing on soil carbon stocks, CO 2 fluxes, and soil chemical properties, focusing on a semi-arid African savanna grassland. We hypothesized that (1) both long-term drought and moderate grazing reduce C storage leading to slower CO 2 emission rates, and that (2) soil chemical properties are interactively influenced by drought and grazing.

Study Site
The study was conducted at the University of Limpopo's Syferkuil experimental farm, Limpopo Province, South Africa (Figure 1). The climate is classified as arid steppe hot climate (Kottek et al., 2006) which receives annual average rainfall ranging from 400 to 600 mm (Benhin, 2006). Average minimum and maximum temperatures are 4-20 • C in winter and 17-27 • C in summer. In the 2018/19 rainy season, a total of 352 mm of rain fall was received from September 2018 to April 2019. Daily temperature, rainfall and other weather date were recorded from an Agricultural Research Council automatic weather station installed at the experimental farm and was about 2 km from the drought act experiment.
The camp (i.e., paddock) on which the experiment was implemented has a size of about 40 ha. It is part of a rotational camp system that is moderately grazed at a stocking density of 9 ha LSU −1 where a livestock unit (LSU) is equivalent to a mature cow with a bodyweight of 450 kg (Meissner, 1983). Grazing periods did not exceed 30 days, with intermittent recovery of not less than 6 weeks during the growing season (8 weeks during the dry season, respectively). The rotational management was established 6 years prior to the set-up of the experiment, resulting in good rangeland condition at the onset of treatments.

Experimental Design
The experiment "DroughtAct" combines a severe drought treatment with a resting treatment. A full factorial experimental design was applied, with crossed drought treatments (two levels: severe drought D+, or ambient rainfall D−) and grazing treatments (two levels: moderate grazing G+ or resting, G−) replicated over four blocks. The resulting 24 plots (10 × 10 m) were separated by corridors with a width of 5 m. Cattle grazing on G-plots was excluded with the aid of permanent wire fences. On D + plots, rainfall was reduced by 66% through the use of large rainout shelters with a size of 36 m 2 (Figure 2A), established at the center of the plots. The rainout shelters allowed for cattle to roam and graze under them. Because the shelters were open on all sides and had a minimum height of 2 m, they had negligible effects on the air temperature beneath them (data not shown). We also oriented the lower sides of the shelters to intercept rain from the dominant wind direction (Carlyle et al., 2014). To prevent interference of lateral soil water movements, we trenched around the perimeter of the rainout shelters to maximum soil depth (≤70 cm) and inserted an impermeable plastic sheet. Sampling was restricted to a central subplot area in each plot (4.8 m × 4.8 m) to reduce marginal effects. Grazing and drought treatments were established in October 2014 and maintained permanently.

Installation of CO 2 Chambers and Measurement of CO 2
In December 2018, two CO 2 chambers were installed in the central subplot of each plot (32 chambers in total). The DroughtAct plots had been under the grazing and drought treatments for four growth periods prior to the installation of the chambers, and treatment effects on vegetation were clearly visible (Figures 2C-F). The CO 2 chambers were constructed and installed in the field according to the USDA-ARS GRACEnet Project Protocols (Parkin and Venterea, 2010). The gas chambers were, however, modified and interfaced with a CO 2 probe GMP343 along with MI70 data logger (Viasala, Vantaa, Finland) to allow for in situ CO 2 measurements ( Figure 2B). The chambers consisted of two separate PVC rings, i.e., a PVC collar (0.20 m diameter, 0.15 m height) and another PVC ring (0.20 m diameter, 0.10 m height) sealed on one side with a PVC circle to make a chamber lid. The chamber lid had a small vent on it to prevent pressure build-up in the chambers during measurements.
The chamber lid was also perforated on top to fix the CO 2 probe. The collars were hammered into the soil to a depth of 0.05 m leaving 0.1 m above the soil. These collars were left in the soil for a few days to settle before measurements were taken.
When taking measurements, the chamber lid with the fixed CO 2 probe on it was attached to the collar and secured with a tube strip to make it airtight. The CO 2 probe GMP343 was set to record measurements every 30 s for 5 min. All measurements were taken between 10:00 h and 14:00 h. CO 2 fluxes were measured over a period of 5 months from January to May 2019 at 2-week intervals, resulting in eight measurement events.

Calculation of CO 2 Fluxes
The CO 2 probe GMP343 gives measurements of CO 2 in parts per million (ppm). The measurements were first converted to mg m −3 using the ideal gas law (Eq. 1): (Collier et al., 2014) (1) Where P = pressure, V = volume, n = moles of gas, R = gas law constant (8.3145 J mol −1 K −1 ), and T = temperature in Kelvin.
The molar volume of an ideal gas at 1atm pressure (101.325 kPa) and 25 • C is 22.4 L mol −1 . Thus, at different pressures, the molar volume of the gas can be calculated as follows: At different temperature and pressure the concentration of CO 2 in mg/m 3 was calculated as follows: Where CO 2 ppm is the measured concentration of CO 2 at any given time, T is the chamber temperature (Temperature in • C + 273.15 K) and P is the ambient pressure. CO 2 concentration in mg m −3 was then plotted against time (min) giving a slope in mg m −3 min −1 . The slope of the resulting regression lines was then determined for each installed chamber. The slope was then multiplied by the volume of the chamber (0.00628 m 3 ) and divided by the area covered by the chamber (0.0314 m 2 ), giving the resultant flux in mg m −2 min −1 , which describes the CO 2 flux out of the soil. Cumulative CO 2 was calculated by assuming that the rate of CO 2 release was constant between two measurement points.

Determination of Bulk Density, Infiltration Rate, and Particle Size
Bulk density was measured using the core ring method. Cores with a diameter of 5 cm and a height of 5 cm were used to collect an undisturbed soil core sample in the topsoil of each treatment plot. The soils were then oven-dried at 105 • C for 24 h before being weighed. Infiltration rate was measured using mini-disk infiltrometers (Decagon Devices, United States). The infiltration measurements were taken once per plot. The particle size was determined by the pipette method (Day, 1965).

Soil Chemical Analyses
Soil samples were collected in plots' central subplot area for two depths i.e., 0-15 and 15-30 cm using augers. To minimize destructive sampling on monitoring plots, we only sampled once per plot. The samples were then air-dried and sieved to pass through a 2 mm sieve. Soil pH was determined in both water and KCl. De-ionized water and 1M KCl solution were used in a ratio of 1:2.5 and measured using a glass electrode. Calcium and Mg were determined by atomic absorption after being extracted using 1M KCl, filtered and diluted with 20 ml of 0.0356M SrCl 2 . Phosphorus, K, Zn, Cu, and Mn were extracted using Ambic-2 solution and determined using atomic absorption. Total C and N were analyzed by the Automated Dumas dry combustion method using a LECO CNS 2000 (LECO Corporation, Michigan).

Determination of Soil Carbon Stocks
Soil carbon stocks are the amount of carbon in a soil layer of known bulk density. Soil carbon stocks are normally given as weight per unit area and are commonly limited to the soil fraction of less 2 mm in size. Soil C stocks are determined by the following formula (Batjes, 1996): Here, soil C stock is soil carbon stock (kg C m −2 ), C tot is the total carbon content (g C g −1 soil), BD is bulk density (kg m −3 ); and d is soil depth (m). Carbon stocks were calculated from the top 5 cm since bulk density was only determined in that depth. This was mainly because the cores were of 5 cm height and there was need to minimize destructive sampling.

Aboveground Net Primary Production
Aboveground net primary production (ANPP) was assessed via destructive biomass sampling. Grass biomass was collected at the time of peak standing biomass (i.e., during May and June) of the years 2018 and 2019. During this time, most plant species present on the plots had fully grown and flowered. On grazed plots, standing biomass was harvested within three 1 m 2 quadrats which were protected from grazing offtake during the growth period with the aid of moveable cages (Figures 3D,F). We created a uniform starting conditions for plant growth at the onset of each growth period by removing residual biomass and carryover material from the previous growth period. This was done by clipping plants to stubble height (∼5 cm). On ungrazed plots, sampling quadrats (3 per plot) only had a size of 50 × 50 cm to minimize the impact of destructive sampling on monitoring plots (see Linstädter and Baumann, 2013). We discarded all moribund plant material, and only collected green and senescent material. Plant material was oven-dried (60 • C, 48 h), and weighed to the nearest gram. ANPP was then estimated from aboveground biomass data via the peak standing crop method (Scurlock et al., 2002;Ruppert and Linstädter, 2014).

Data Analysis
Two-way ANOVAs were performed to determine the effect of grazing, drought and their interaction on C stocks, CO 2 fluxes (individual points, cumulative and average rates), and chemical soil properties. Where effects of factors or their interaction were significant, Tukey HSD was used for post hoc comparisons. Linear regression analysis was done to determine the impact of ANPP, bulk density, and C stocks on CO 2 flux rates. We also performed linear regression between C stocks and CO 2 flux rates to search for dependencies. All data were analyzed using SPSS 25 (SPSS, United States).  (Figure 3). Daily rainfall and mean daily temperature data indicate that CO 2 fluxes were higher when rainfall was higher (Figure 3). This is clearly shown by the sharp decrease in CO 2 fluxes in mid-March when there was a month-long period of sporadic and ineffective daily rainfall (<5 mm). The CO 2 flux rate then sharply rose after 35 mm of rainfall were received. In contrast, air temperature did not influence the variation in CO 2 fluxes. A clear trend toward higher CO 2 emissions from grazed plots was discernible, the differences were, however, not significant (Table 1 and Figure 4). An exception was one sampling event in February, where grazed plots released significantly more CO 2 than rested plots (12.17 ± 4.86 mg m −2 min −1 compared to 9.95 ± 3.55 mg m −2 min −1 ; Figure 3A, marked with a star). Contrastingly, drought significantly reduced CO 2 flux almost for all measurement events (Figure 4) until the beginning of the dry season in May, when the measured CO 2 fluxes were almost equal on the last sampling event (Figure 3B, marked by a star). Average CO 2 emission from drought plots were 32% lower than those under ambient rainfall.

CO 2 Emission as Influenced by Grazing and Drought
The interaction effect between grazing and drought was not significant for CO 2 fluxes (Table 1). However, Figure 4B clearly  showed that the impact of drought overrode the impact of grazing as treatments that included drought (D + G− and D + G+) had lower CO 2 flux rates compared to those that excluded drought (D−G− and D−G+). Looking at each treatment individually, cumulative CO 2 emission were higher under ambient rainfall with moderate grazing followed by ambient conditions with grazing exclusion while the least emissions were observed for the treatments that were under drought ( Figure 3C). By the end of the growth period, grazed plots under ambient rainfall (D−G+) had released 16.3 tons of CO 2 ha −1 which was 58% higher than observed on grazed plots subjected to severe drought (D+ G+).
The results also showed that the average CO 2 flux rate under ambient conditions was 8.3 mg m −2 min −1 equating to 43.6 tons CO 2 ha −1 year −1 (assuming the rate remained constant throughout the year), while under drought, the CO 2 flux rate was 5.7 mg m −2 min −1 , resulting in 30 tons CO 2 ha −1 year −1 .

Effect of Drought and Grazing on Infiltration Rate, Bulk Density, and Carbon Stocks
Drought did not affect bulk density of the top 5 cm of the soil but grazing did (Table 1 and Figure 5A). Bulk density was higher under grazed plots (1.53 ± 0.07 g cm −3 ) compared to ungrazed plots (1.39 ± 0.08 g cm −3 ) ( Figure 5A). Infiltration rates also had a similar result where infiltration was also not influenced by drought but only with grazing, whereby the infiltration rate was higher under grazed plots compared to rested plots (results not shown). Soil carbon stocks, on the other hand, were significantly affected by drought (P < 0.05) but not by grazing (P > 0.05) ( Figure 5C). Soil carbon stocks were higher under drought (1.118 ± 0.130 kg C m −2 ) compared to ambient conditions (0.999 ± 0.115 kg C m −2 ). Although no variations were observed, soil carbon stocks tended to be relatively higher under grazed plots compared to ungrazed ones. Only bulk density was affected by the combined effect of grazing and drought ( Figure 5). Bulk density was higher under drought and grazing with a bulk density of 1.59 g cm −3 . Soil carbon stocks did not respond to the interaction effect of grazing and drought but was also relatively higher under drought and grazing (D + G +).

Relationships Between C Stocks, ANPP, Bulk Density, and CO 2 Flux
Regression analyses (Figure 6) were performed to enhance the understanding of the factors that affect carbon storage and CO 2 flux in the soil. The relationship between bulk density and soil C stocks was positive but not significant. Soils with high bulk density tended to have a higher carbon content. The relationship between soil C stocks and CO 2 flux revealed a negative association. ANPP from the previous growth period (2017/18) did not show any relationship with soil C stocks, while ANPP from the same growth period (2018/19) explained about 20% of the variation observed in soil C stocks ( Figure 6C). While soil C stocks slightly decreased with an increase in biomass production, CO 2 fluxes increased with an increase in production. ANPP of the previous growth period and of the current growth period explained 32% and 48% of the variation observed in CO 2 fluxes, respectively.

Soil Chemical Characteristics as Affected by Drought Treatments
Chemical properties in the topsoil (0-15 cm) did not significantly differ from properties in 15-30 cm depth (results not shown). However, visual observations revealed that P, K, C, Cu, and Zn tended to be relatively higher in the top 15 cm compared to the 15-30 cm, while Ca and Mg were relatively higher in the 15-30 cm layer compared to the 0-15 cm layer. Due to the lack of significant differences between the two depths, data were combined and jointly analyzed for the effect of drought and grazing ( Table 2) and their combined effect ( Table 3).
The results showed that N and P were lower in the grazed plots compared to the ungrazed, while Mg and Mn were higher in grazed plots. The concentration of potassium, Ca, Zn, Cu, and the soil pH were not influenced by grazing. On the other hand, drought significantly increased N, P, K, Zn, and Mn concentrations in the topsoil. Soil pH was also 0.16 units higher in the plots that were under drought compared to those that were under ambient rainfall conditions. The combined effect of drought and grazing had a significant influence on P, K, Zn, N, and on C:N ratios (Table 3). Phosphorus, Zn and N were all higher for plots under drought and ungrazed conditions (D + G−) compared to the other three treatments. Similarly, the C:N ratios were lowest under D + G−. Potassium, however, was highest under drought and grazed (D + G +) with a concentration of 237 mg kg −1 and lowest in ambient and grazed (D−G+) with a concentration of 140 mg kg −1 . Even though the combined effect of drought and grazing did not influence carbon, D + G− plots tended to have higher total carbon in soil compared to the plots. Soil pH was slightly acidic, ranging from 4.5 to 4.8 and did not respond to the impact of drought and grazing.

Effect of Drought and Grazing on Soil Carbon Storage and CO 2 Fluxes
This study used a large field experiment in a dryland grassland to provide novel insights into soil carbon storage and release under long term exposure to extreme drought and moderate grazing. While grazing had no clear effect, drought reduced CO 2 fluxes and increased soil C stocks. Intuitively, ambient conditions would be expected to have higher carbon storage due to higher biomass accumulation (see Figure 2) and higher photosynthetic rates (Li et al., 2017), which should result in more carbon input into the soil than is lost through respiration. There are indeed studies that have shown a higher soil organic carbon accumulation due to higher plant biomass production (Zhao et al., 2018;Prommer et al., 2020). The results of the current study are consistent with other studies comparing decomposition rates under drought and ambient conditions. Sanaullah et al. (2012) reported more than 50% decrease in grass litter decomposition due to drought when compared to ambient conditions. The lower CO 2 flux Significance levels: *P < 0.05, **P < 0.01, ***P < 0.001, ns means not significant. D− G− means ambient rainfall and ungrazed (fenced); D− G+ represents ambient rainfall and grazed; D + G− is drought and ungrazed (fenced) and D + G + is drought and grazed. Significance levels: *P < 0.05, **P < 0.01, ***P < 0.001, ns means not significant. Different letters in the same column indicate significant differences.
rate observed under drought conditions could possibly mean that less soil carbon was lost thus maintaining higher carbon storage ( Figure 5C). Decomposition rates are known to be slower under dry conditions (Sanaullah et al., 2012;Lefèvre et al., 2017) mainly because microbial activity decreases under water stress conditions (Manzoni et al., 2012;Schimel, 2018). Therefore, the observed higher soil C stocks under extreme drought could probably be attributed to the low CO 2 emission rates under drought. This means that less C was lost when compared to ambient conditions. In addition, the higher soil C stocks observed under drought could also be partially explained by the higher concentration of multivalent cations such as Ca, Zn, and Mn. These multivalent cations have been reported to contribute to the stabilization of soil organic carbon in the same study area (Dlamini et al., 2019). Multivalent cations, such as Ca, Al, Fe, and Mn offer chemical protection of soil organic carbon through the formation of organo-mineral complexes (Jastrow et al., 2007).
Studies have shown that most of the CO 2 released from the soil is due to respiration by roots (autotrophic) and/or soil microorganisms and fauna (heterotrophic) (Hanson et al., 2000;Bond-Lamberty et al., 2004), both of which are sensitive to soil moisture (Li et al., 2018;. Low soil moisture content reduces microbial biomass (Liu et al., 2009). It also hinders root growth due to inhibition of shoot growth and the subsequent decrease in translocation of assimilates to the roots and the rhizosphere (Li et al., 2018). Furthermore, drought may massively reduce root exudates which promote rhizospheric respiration. The current results point to a reduction in both autotrophic and heterotrophic respiration due to drought. However, from the field observations, plots under drought conditions had very sparse vegetation if not completely bare (Figures 2E,F), suggesting a higher contribution of heterotrophic respiration to the measured CO 2 flux than autotrophic respiration. In support, Balogh et al. (2016) found that rhizospheric respiration, a component of autotrophic respiration was highly sensitive to drought. Moreover, Hanson et al. (2000), asserts that the root contribution to total soil respiration is higher during the growing season and lower during the dormant periods of the year while  allude to the stronger relationship between plant growth and root respiration.
Grazing did not significantly influence both soil carbon storage ( Figure 5C) and CO 2 flux ( Figure 4B). These findings are contrary to the findings of Li et al. (2017). However, soil carbon storage was relatively higher in grazed plots compared to ungrazed plots which is in agreement with the findings of Sun et al. (2019). Increased belowground biomass under moderate grazing is reported to be due to the translocation of photosynthates to the roots (Shi et al., 2013;Sun et al., 2019). It is also believed that moderate grazing promotes root-associated microbes and fauna that leads to an increase in below-ground biomass. Li et al. (2017) on the other hand found that grazing reduced root biomass compared to ungrazed plots.
In this study we did not observe a strong confounding effects of drought and grazing on CO 2 flux and carbon stocks after a 5year long extreme drought. As observed in Figure 4B, the impact of the long extreme drought seemed to override any impact that grazing might have had. This is also evident on Figures 2E,F where there is rarely any vegetation on plots exposed to the extreme drought. It can be postulated that the interaction might have been significant and additive in the earlier stages of the experiment but got progressively weaker to non-existence with time due to the deleterious effects of extended extreme drought. In support of our postulation, Okach et al. (2019b) reported a minimal impact of grazing on plots with less rainfall compare to those under higher rainfall. Furthermore, Zhou et al. (2019b) also reported higher impact of drought on CO 2 flux compared to that of grazing. However, grazing and drought can have confounding impact on carbon storage and soil CO 2 flux (Okach et al., 2019b;Zhou et al., 2019a,b).
With reference to the challenges of global climate change and the need to reduce greenhouse gas emissions, this study has provided inventory data on the CO 2 emission levels in African dryland grasslands. The results show that by early May, ambient rainfall with moderate grazing had resulted in cumulative CO 2 emissions of about 16.3 tons ha −1 . These values are comparable to those reported in other parts of SSA (Kim et al., 2016;Berger et al., 2019). Kim et al. (2016) reported emission rates of 3.3-57 tons CO 2 ha −1 year −1 in natural terrestrial systems of the SSA of which grasslands emitted an average of 15.5 ± 3.8 tons CO 2 ha −1 year −1 . On the other hand, we found that ungrazed dryland grasslands to store about 0.996 kg C m −2 (i.e., 9.7 tons C ha −1 ), while grazed ones stored about 1.122 kg C m −2 (i.e., 11.2 tons C ha −1 ). These quantities are in line in with those obtained by Dlamini et al. (2019) in the same area which ranged from 0.40 to 1.87 kg C m −2 .

Relationships Between C Stocks, ANPP, Bulk Density, and CO 2 Flux
Our study showed that soils with higher bulk density had higher soil C stocks (R 2 = 0.19, P = 0.10). This is because soils with higher bulk density are compacted and have a curtailed air filled porosity (van Verseveld and Gebert, 2020). The curtailment of air filled porosity reduces oxygen availability to microbes leading to reduced carbon oxidation and thus preserves carbon in the soil. Neilson and Pepper (1990) found reduced CO 2 evolution on soils with higher bulk density. Consequently, a negative relationship between soil C stocks and CO 2 emission (R 2 = 0.25; P = 0.04) was observed. ANPP on the other hand was positively related to CO 2 emission rate with a third to almost half of the variation observed in CO 2 emission rate being explained by it. Higher ANPP suggests an increase in root biomass and thus also an increased root respiration (Luo and Zhou, 2006). According to the study by Luo and Zhou (2006), root respiration contributes about half of the total soil respiration. Finally, we found ANPP to be negatively related to carbon stocks. This is at the first glance surprising, but is in line with recent studies stating that ANPP may either increase or decrease soil organic carbon in grasslands, depending on the net result of plant-C inputs and C decomposition (Peterson and Lajtha, 2013;Chen et al., 2018).

Effect of Drought and Grazing on Soil Properties
This study provides an indication of combined effects of drought and grazing on soil chemical properties other than the commonly reported C and N (e.g., Okach et al., 2019a). The results revealed that grazing decreased P and N in the topsoil but did not influence other chemical properties such as pH, K, and the C:N ratio. The decrease in nitrogen due to grazing is a common phenomenon in heavily grazed grasslands (Han et al., 2008;He et al., 2019). Under low grazing intensities, the concentration of N in soils is usually unaffected due to the balance between N removal and N returned. It is believed that livestock returns the majority of the ingested N to the soil in the form of urine and feces (Han et al., 2008). In addition, the trampling by livestock increases the incorporation and the microbial mineralization of the litter thus returning C and N to the soil (Li et al., 2017). On the other hand, the effect of grazing on P is not much reported. The only study we found that reported on P showed the P response to grazing is dependent on mean annual precipitation. It was shown that grazing reduces soil P when precipitation is lower than 400 mm but increases when it higher (He et al., 2019).
While grazing reduced soil P, drought, on the other hand, resulted in increased soil P. Drought also increased other soil chemical properties such as pH, K, Zn, Mn, N as well as lowering the C:N ratio ( Table 2). The higher concentrations observed under drought could be attributed to less leaching of the nutrients as the plots under drought received less than 40% of the ambient rainfall. The generally low pH observed in this study is a testament to the prevalence of basic ions leaching from the topsoil. Leaching is a problem particularly in light-textured soils (Matichenkov et al., 2020). Also, the higher concentrations of some of the determined soil chemical properties could just be a result of reduced nutrient mining by plants as drought impeded growth and nutrient uptake. Nutrient mining which is the removal of nutrients from the soil by plants (Vanlauwe et al., 2015;Kodzwa et al., 2020) and it can be an issue in poorly managed dryland grasslands due to inefficient nutrient cycling (Bilotta et al., 2007;Han et al., 2008).
In addition, drought and grazing acted together to influence the concentration of macronutrients (NPK) as well as Zn and C:N ratio in the top soil. Exposure to drought with grazing (D + G−) resulted in increased concentration of P and N and subsequently lower C:N ratios. These findings of increased N and relatively lower C:N ratios under D+ G− are contrary to the findings of Okach et al. (2019a). It is not clear why N, P, and Zn were higher under D + G− but this could be the attributed to reduced leaching and no removal by feeding animals.
Grazing increased the bulk density of the surface soil. The average bulk density of grazed plots was 1.54 g cm −3 while that of ungrazed plots was 1.39 g cm −3 . These results are in agreement with findings from previous studies (Chen et al., 2017;Li et al., 2017;Okach et al., 2019a;Sun et al., 2019;Zhou et al., 2019a). Soils may become compacted because of trampling by livestock, especially when grazing takes place in damp areas. According to Bilotta et al. (2007), hooves of an adult cow weighing between 300 and 600 kg can exert static pressures of approximately 200 kPa on the soil, which is a high enough force to compact soil. Thus, if the stocking rate is high and the soil is wet, trampling could be highly detrimental for soil physical properties. The effect of drought on bulk density has rarely been reported and in this study it was observed that drought tended to increase the bulk density of the soil compared to under ambient conditions. Slightly higher bulk densities under drought could be due to the lack of vegetation which exposes the soil to compaction from external forces such as raindrop impact (Neave and Rayburg, 2007) and livestock trampling (Bilotta et al., 2007). Also, plant roots enhance soil aggregation, improving soil structure and hence reducing bulk density (Rillig et al., 2015).

CONCLUSION
This study provides an indication of carbon storage and emission rates in a dryland grassland. Both grazing and severe drought had a positive effect on soil carbon storage but had dissimilar impacts on CO 2 fluxes. Drought drastically reduced the amount of CO 2 released from the soil, while moderate grazing did not influence it. At the end of the growth period, grazing under ambient rainfall had released 16.3 tons of CO 2 ha −1 which was 58% higher than observed on grazed plots subjected to severe drought. Moderate grazing reduced P in the soil compared to the ungrazed plots. On the other hand, drought increased P concentration in the soil when compared to under ambient conditions. Our results imply that long-term extreme drought in dryland grasslands may be favorable in terms of preserving the existing soil carbon stocks through reduced CO 2 release to the atmosphere. This finding is critical in understanding future soil carbon dynamics in dryland grasslands in the face of climate change.

DATA AVAILABILITY STATEMENT
The datasets presented in this article are not readily available because some data were used as predictor variables and that data still needs to be published. Requests to access the datasets should be directed to AL, anja.linstaedter@uni-koeln.de.

AUTHOR CONTRIBUTIONS
AL and KKA conceived and designed the experiments. LM and TM redesigned some parts of the experiments. LM, KB, EM, and MM carried out the field and lab-work. LM analyzed the data. LM wrote the manuscript. All authors commented on the analysis, interpretation, presentation of the data, and critically revised the manuscript.