Introduction
The High Plains of Texas have a semi-arid climate where plant growth and agricultural production are limited by high summer temperatures and low water availability. Increased climate extremes such as increased temperature, and more frequent and severe droughts driven by changes in precipitation patterns, have further exacerbated the agricultural soil degradation in this region. As a result, the physical, chemical, and biological aspects of soils have been severely affected, especially in dryland soils. Despite this soil degradation, the total cotton acreage planted in the Texas High Plains has increased over the past few years. However, this increased cotton acreage has not led to a concurrent increase in total cotton production [
1]. Furthermore, groundwater, which is the region’s main source of irrigation, is being depleted at an increasing rate [
2]. Thus, growers are being forced to reduce irrigation or switch to dryland (i.e., no irrigation) cultivation, a practice that could reduce yields by up to half [
3]. Regenerative agricultural practices, for example no tillage, crop rotation, residue retention, and cover crops, have been proposed as a potential mitigation approach as these practices are believed to stabilize the soil microenvironment, thereby maintaining better soil conditions for plant and microbial activity. Several cotton growers in the High Plains of Texas have already taken the initiative to incorporate these practices in their cropping systems [
4]. Nevertheless, little is known about how effective these practices would be in moisture-limited arid and semi-arid environments to improve soil health and crop productivity under future projected climate change.
Soil organic matter content in the soil is considered a key soil health indicator. It promotes plant growth by supplying nutrients, improves soil aggregate formation [
5], improves water retention [
6], and supports soil biological activity [
7]. Soil carbon in agricultural soils is declining because of unsustainable land management practices and increasing extreme climate events such as extreme temperature and precipitation. In total, agricultural soils have lost 25-75% of their soil organic carbon pool [
8]. In arable agricultural soil, depending on the extent of soil manipulation and farming technique used, cultivation can positively or negatively impact soil organic carbon storage and release [
9]. For instance, Luo et al. [
10] showed that soil organic carbon at 10 cm below the surface in cultivated land was 51 percent lower than in a natural ecosystem after five decades of farming, indicating that farming decreases organic matter in the soil. Further climate related changes exacerbate this soil carbon loss. Therefore, it becomes crucial to investigate the various factors that have the potential to influence soil organic matter stocks and carbon loss in a warmer world. By understanding these key driving forces, we can effectively address the challenges posed by climate change, ensuring agricultural sustainability and maintaining a balanced ecosystem.
Soil organic matter in the cultivated soil can be increased by adopting regenerative agricultural management practices such as no-tillage, cover cropping, mulching, residue retention, and crop rotation [
11,
12]. Residue mulching is an important agronomic practice that involves the covering of soil surfaces using organic or inorganic materials. Organic mulches include straw, husk, grasses, compost, and plant residue, while polyethylene plastic mulch is the most used inorganic mulch [
13]. The use of mulches helps to minimize water runoff, improve infiltration, and regulate soil temperature and moisture [
14]. Organic mulches such as grasses or plant residues not only regulate the soil environment but also contribute to soil health by providing carbon input and plant nutrients, enhancing biological activity in the soil [
14]. Residue also reduces evaporation and enhances the water retention rate of the soil, minimizing moisture loss due to excessive evaporation at higher temperatures [
15]. Several studies have shown that adding residue to the soil can increase soil organic matter content [
16,
17]. Conversely, it is important to note that the increased presence of active soil organic matter resulting from residue addition may accelerate the decomposition rate and lead to carbon loss to the atmosphere [
18]. Nevertheless, this process keeps the soil system active and dynamic by promoting microbial growth and activities, soil aggregate stability, and continuous recycling of nutrients in the soils.
Climate change can, directly and indirectly, affect soil health by changing the amount of carbon stored in soils and soil biological activity [
19,
20] although whether it causes net carbon loss or net increases soil carbon is still debated [
21,
22,
23]. Warming stimulates soil respiration, organic matter decomposition, and nutrient mineralization [
24,
25,
26,
27], thereby releasing more carbon from the soil as CO
2 [
23,
28]. Warming also increases plant carbon assimilation, which can enhance soil carbon inputs [
29,
30,
31]. The net change, expressed as the difference between increased carbon loss and increased net primary production in response to warming, determines whether carbon is stored or released from the soil in a warmer world [
30,
32].
Temperature is not the sole factor influencing soil organic matter decomposition and soil respiration; moisture, microbial abundance, availability and accessibility of soil microbes to the substrate, enzyme activity, and soil properties all play a crucial role [
33,
34]. Soil organic matter decomposition and soil respiration both increase with temperature if all other factors remain constant [
22]. However, in a natural environment, temperature interacts with various other factors within the soil system [
35] with the temperature–moisture interaction being the most important factor determining the soil carbon response to warming. Warming reduces soil water availability by increasing evapotranspiration and decreasing soil moisture [
36], thereby lowering the rate of organic matter decomposition, even to the point at which soil respiration no longer responds to warming [
27]. Soil moisture regulates the warming induced daily temperature fluctuations in the soil. Higher soil moisture increases the specific heat capacity of the soil, which increases the amount of heat needed to raise the soil temperature [
37,
38]. Dry soils show quicker and larger temperature variation than wet soils under climate extremes. Hence, the complexity of the temperature–moisture interaction in the soil, and its subsequent effects on soil organic matter decomposition need further attention to better understand the effects of future climate change on soil carbon dynamics.
Soil microorganisms and their activity are another important factor influencing soil carbon dynamics by driving plant litter decomposition and soil organic matter formation or directly contributing to the soil carbon pool through microbial biomass [
19]. Warming alters plant growth, litter production, and root-derived carbon via root exudates, stimulating soil microbial growth and activity [
39,
40]. Microbial biomass carbon is highly correlated with plant-derived carbon via root exudation and decomposition [
41]. Warming may increase [
42], decrease [
11,
43], or have no effects [
44] on microbial biomass carbon. The precipitation pattern, which may govern soil moisture regime and substrate availability, influences the response of microbial biomass carbon to warming. Microbial biomass is negatively correlated with warming when soil moisture is a limiting factor, but not under abundant moisture conditions [
45]. Therefore, the microbial contribution to soil organic matter is sensitive to temperature–moisture interactions and its resulting effects on microbial growth and activity.
While retaining crop residue from the cash crop is a common management practice in the Texas High Plains, some growers have recently shown interest in adding extra residue, such as planting perennial grasses between crop rows and terminating them shortly after crop germination. Here, we evaluated whether the physical aspect of adding dried grasses to the soil surface would show potential as a beneficial soil amendment practice in agriculture due to its role as a biodegradable cover and carbon source. A few previous studies have already shown that adding a layer of dried grasses on the soil surface is beneficial in improving irrigation efficiency and reducing the irrigation water demand in cotton farms in semiarid ecosystems [
46,
47]. In our study, we evaluated the effectiveness of using multispecies dried grass mulching (referred to as residue addition hereafter) as a viable strategy for reducing temperature and moisture fluctuations and increasing organic matter in the soil profile, thereby minimizing soil health degradation during climate extremes.
We examined the effects of summer warming and residue addition on soil carbon dynamics and cotton yield in both irrigated and dryland soils of the semi-arid Texas High Plains. Irrigation would not be a sustainable crop management strategy, because of dwindling water sources in semiarid regions, including the rapidly depleting Ogallala aquifer in our study area [
48]. Hence, we used the irrigated fields to examine whether residue addition and warming effects were similar under drier versus wetter conditions. We hypothesized that residue addition would lower daily fluctuations in soil temperature and soil moisture and reduce evaporation rates such that soil moisture levels would be higher in residue-added plots compared to plots without residue. We also hypothesized that organic matter content and soil carbon respiration would be greatest in irrigated, warmed plots with added residue due to the increased decomposition rate resulting from temperature-induced changes in microbial enzyme activity and the increased carbon substrate availability from the added residue. Overall, this study aimed to evaluate the importance of multispecies grass residue addition on buffering the negative impacts of soil temperature and moisture extremes on soil organic matter, cotton biomass, and yield. If so, mulching could be implemented as one of many practices aimed at mitigating the negative effects of climate change on cotton production.
Materials And Methods
Site Characteristics
The research was carried out during the growing season of 2021 at the Texas Tech Quaker Avenue Research Farm, Lubbock, Texas (33° 41’ 36.4596” N, -101° 54’ 18.612 “W, 992 m a.s.l.). The study site was in a semi-arid climate with a 30-year mean annual precipitation (MAP) of 466 mm and a mean annual temperature of 16.3 ⁰C [
49].The hottest month was July, with an average monthly temperature of 27.3 ⁰C, and the coolest month was January, with an average monthly temperature of 5.1 ⁰C [
49]. A weather station installed in the center of the research field was used to record field-level temperature, precipitation, relative humidity, and wind speed. During the growing season of 2021 the average temperature was 24.2 ⁰C (the hottest month was June, with an average monthly temperature of 27.3 ⁰C, and the coolest month was October with an average monthly temperature of 18.21 ⁰C). During our study period (June through October), the field received a total rainfall of 337 mm, i.e., 72.3% of MAP. The mean soil pH was 8.49. The soil had 1.042 + 0.10 % organic matter and a bulk density of 1.29 g/cm
3 at 0-10 cm depth. The soil had a sandy clay loam texture with 61.45 % sand, 15 % silt, and 23.55 % clay. The soil was classified as Amarillo-Acuff sandy clay loam (Fine-loamy, mixed, superactive, thermic Aridic Paleustalfs).
Experimental Design
The field was divided into two adjacent field sections, irrigated by a drip irrigation system and non-irrigated (i.e., dryland). Prior to this experiment, both sections had been operating under an irrigated cotton monocropping system. During the experiment period, the irrigated section received drip irrigation in addition to rainfall, while the dryland section had no additional irrigation (i.e., rainfall was the sole water source). During the growing season, a total of 218 mm of irrigation water was provided via drip lines to the irrigated section. There was a 4 m buffer zone between the irrigated and dryland sections. Each irrigation section was then divided into 3 blocks each (a total of 6 blocks) to capture the spatial gradient in soil properties. There were eight 1 meter × 1 meter plots in each block. The passive warming treatment was installed during the growing season. The warming treatments were implemented using 1m ×1m × 1m open-top chambers (OTC) made of aluminum rods and clear polycarbonate sheets. We set up the OTCs in the field immediately after sowing cotton seeds, using stakes and zip ties to secure them to the ground. In the plots with residue treatments, multispecies grass residue (Bermuda (Cynodon dactylon (L.) Pers.), blue grama (Bouteloua gracilis (Kunth) Lag. ex Griffiths), and fescue grasses (Festuca arundinacea Schreb.)) was added to the soil surface at a rate of 3 kg residue/m2. The added residue was covered by plastic garden netting to keep the mulches in place. Cotton (variety: Phytogen 394) was planted in early June continuously in a row and harvested in late October. Each plot had a single crop row containing 7-8 cotton plants spaced approximately 10 -12 cm apart.
Measurement of Environmental Variables
5TM sensors linked to EM50 data loggers (Meter Group, Inc., Pullman, Washington, USA) were used to record soil temperature and volumetric moisture content every 30 minutes at 10 cm soil depth in each plot. Ibuttons (Maxim Integrated, California, USA) were used to record air temperature and relative humidity every four hours at the canopy level of mature cotton plants (i.e., 50 cm above the ground). The ibuttons were inside radiation shields to prevent the heating of sensors from direct solar radiation. Ventilation of the radiation shields was achieved by overlapping two perforated plastic funnels in such a way that the holes in one funnel did not line up with the holes in the other. On top of the radiation shield, we placed HOBO Pendant Temperature/light data loggers (MX2022; Onset Computer Corp., Massachusetts, USA) to monitor the amount of light intercepted at the canopy level.
Soil Sample Collection and Laboratory Analysis
Soil samples were collected from each plot shortly after crop harvest in late October. We took samples from 0-15 cm deep with a soil corer (3 cm diameter). Two soil samples were collected per plot, one from each side of the crop row in the plot. For each soil sample, the soil was taken from three randomly selected sites within a side of the crop row in the plot and mixed to create one composite sample (i.e., one of two composite samples per plot). As a result, the field yielded a total of 96 soil samples from 48 plots. The soil samples were transferred to the laboratory in a refrigerated container. The samples were kept at 4 ⁰C after passing through a 2-mm sieve to remove bigger plant roots, debris, and stones and analyzed by Waters Agricultural Laboratories Inc. for soil macro- and micronutrients, soil organic matter, pH, and cation exchange capacity. Soil inorganic nitrogen availability (NH4+-N and NO3--N) were extracted with 2 M KCl and measured on a FIA analyzer (FIA Lab Instruments, Inc, Seattle WA). Other nutrients were extracted with a Mehlich III solution and analyzed on an ICP. Soil organic matter was measured using the loss-on-ignition method at 350oC for 2 hours.
Microbial biomass was measured using the chloroform fumigation extraction procedure [
50]. Four, 5 g dry weight equivalent soil samples were weighed in glass beakers, two of which were fumigated for 48 hours with 25 ml of chloroform and the other left unfumigated. Extractable carbon was extracted from fumigated and non-fumigated samples using 50 ml of 0.5 M K
2SO
4 and filtered through Whatman 43 filter paper. We measured the extracts at 280 nm wavelength using GENESYS 150 UV-Visible Spectrophotometer (ThermoFisher Scientific, Madison, USA). The difference in absorbance between the fumigated and unfumigated samples was used to calculate soil microbial biomass [
51], using a K
EC (extractable fraction of microbial biomass carbon) value of 0.45 for the calculation [
52].
Soil Respiration Measurement
We used the LI-8100A soil CO2 flux system (LI-COR Inc, Nebraska, USA) to measure soil respiration rates monthly starting in July. To enable a good seal with the soil surface for soil respiration measurements we installed soil collars in each plot at the beginning of the experiment; the 20-cm diameter soil collars were installed 2-3 cm deep into the soil in the middle of the plot, 5 cm apart from the crop row. Plant structures inside the soil collar were periodically removed to exclude aboveground plant tissue respiration. The measurement time was set to 2 minutes for each soil respiration measurement. Soil respiration data were taken during the same time window, from 8:30 AM to 11:30 AM, to eliminate measurement variability due to time-of-day.
Harvesting and Biomass Measurements
Cotton was harvested by hand in late October when most of the bolls were fully open. Cotton bolls were harvested from all plants within a plot and stored in separate plastic bags. The weight of the harvested seed cotton was recorded after it was air-dried for a week. The total number of plants per plot and the number of bolls in each plant were also recorded.
Plant biomass in each plot was recorded after a week of drying to remove all moisture. To quantify relative changes in belowground root biomass, we used soil cores with a diameter of 3 cm and a length of 10 cm. Three plants were chosen at random within each plot, and two soil cores of root samples (one sample from each side of the plant row) were obtained. Root samples were taken at 3-5 cm from the plant stem. As a result, six root samples were collected from each plot. We used a 2 mm sieve to separate roots from the soil. The roots were hand-picked from the sieved sample, washed, and dried before taking dry root weight.
Statistical Analysis
We evaluated the interaction between warming, residue, and irrigation treatments on soil temperature, air temperature, volumetric soil moisture content, and soil organic matter using linear mixed effects models in R [
53]. Microbial biomass and soil respiration rate were evaluated using generalized linear mixed effects models. For these variables, the residuals showed non-normal error distribution, and hence we chose generalized linear mixed effects models. The distribution for the generalized linear mixed model was selected based on AIC values. The best model for microbial biomass carbon and soil respiration had a log-linked gamma distribution and inverse-linked gamma distribution, respectively, which improved the behavior of residuals and had lower AIC values. Seed cotton yield, the number of bolls per plant, aboveground biomass, and belowground biomass was analyzed using linear mixed effect models. We used ‘lmer’ and ‘glmer’ function in ‘lme4’ package [
54] for linear mixed effects, and generalized linear mixed effects models, respectively. Since we took two soil samples from each plot, the data from the two samples were averaged to get plot level data before fitting the model. Blocks were included as a random intercept term in each of the models. For repeated time measurements (soil temperature, air temperature, volumetric water content, and soil respiration rate), we first calculated monthly averages for each plot, and then included the month and plots as an additional random intercept term. We also fit a separate model to evaluate the effects of climate data and soil variables for each response variable. First, we shortlisted a few predictors for each of our response variables based on the literature, then, we fit the liner mixed effects models in R. We used ‘car’ package [
55] to generate the ANOVA tables and p values for fixed effect predictors. Following that, post hoc analyses were performed using Tukey’s HSD with a 95% confidence interval to determine if there were significant differences between treatments. We used the ‘emmeans’ package [
56] for post hoc analysis. The ‘ggplot2’ package [
57] was used to visualize the data.
Results
We observed strong temporal variation in soil and air temperature throughout the growing season (
Figure 1). OTCs (
χ2 = 498.31,
P < 0.001) increased average air temperature by 2.2 ⁰C (
Figure 1b) but did not affect soil temperatures. Rather, soil temperatures were more affected by residue (
χ2 = 17.18,
P < 0.001) and irrigation (
χ2 = 18.51,
P < 0.001). Residue decreased soil temperature by 0.5 ⁰C, while irrigation decreased soil temperature by 0.7 ⁰C (
Figure 1a). Irrigation also decreased the air temperature by 0.6 °C (
χ2 = 12.64,
P = 0.003). Residue addition decreased (
χ2 = 80.15,
P < 0.001) the average daily temperature range in the soil system by 1.4 °C (
Figure 2).
Volumetric water content also showed a temporal fluctuation throughout the growing season (
Figure 3a). Irrigation and OTCs did not have a significant main effect on volumetric water content, but residue (
χ2 = 8.56,
P = 0.003) significantly changed volumetric water content. Our results also showed a significant three-way interaction between OTCs, residue, and irrigation (
χ2 = 4.96,
P = 0.025). In dryland, in the presence of OTCs, residue increased volumetric water content by 9.16 %, but decreased by 18.12 % when OTCs were not present (
Figure 3b). In irrigated fields residue decreased volumetric water content irrespective of OTCs treatment.
OTCs and residue did not affect soil organic matter content. However, there was a significant effect of irrigation on soil organic matter (
χ2 = 8.05,
P = 0.0046), with dryland soils having 36.8 % lower soil organic matter compared to irrigated soils (
Figure 4a). Soil organic matter was also negatively correlated with soil temperature (
χ2 = 4.02,
P = 0.044;
Figure 5a). We observed a significant interaction effect of OTC and residue (
χ2 = 7.37,
P = 0.0066) on microbial biomass: OTCs increased microbial biomass by 34.9 % under residue-added condition, but OTCs had no effect on microbial biomass when residue was not applied (
Figure 4b). Irrigation (
χ2 = 4.73,
P = 0.029) increased microbial biomass by 27.5 %. Additionally, a significant positive relation was observed between microbial biomass and organic matter (
χ2 = 8.39,
P = 0.003;
Figure 5b), as well as microbial biomass and available nitrates (
χ2 = 7.42,
P = 0.0064;
Figure 5c).
Residue addition (
χ2 = 72.84,
P < 0.001) significantly increased soil respiration; the residue-added plots had a 78.2 % higher soil CO
2 flux rate than plots without residue (
Figure 4c). We also observed a significant interaction effect between the irrigation and warming treatments on soil respiration (
χ2 = 5.64,
P = 0.017): OTCs increased soil respiration in dryland by 35.1% but had no effect in irrigated plots.
Seed cotton yield was not affected by OTCs, but irrigation (
χ2 = 6.87,
P = 0.0087) and residue (
χ2 = 4.83,
P < 0.027) had a significant impact on seed cotton yield. The residue addition and irrigation increased seed cotton yield by 15.2 % and 37.3 %, respectively (
Figure 6c). Moreover, both OTCs and residue addition did not change aboveground biomass and belowground biomass. Irrigation, however, increased aboveground biomass by 150.5 % (
χ2 = 22.61,
P < 0.0001;
Figure 6a) and belowground biomass by 129.7% (
χ2 = 9.01,
P < 0.0026
Figure 6b).
Author Contributions
Conceptualization, P.D. and N.V.G.; Methodology, P.D, R.K.S, N.G.S, L.C.S, N.V.G.; Formal Analysis, P.D.; Investigation, P.D, R.K.S, N.V.G.; Resources, P.D, R.K.S, N.G.S, L.C.S, N.V.G .; Data Curation, P.D., R.K.S; Writing – Original Draft Preparation, P.D.; Writing – Review & Editing, P.D, R.K.S, N.G.S, L.C.S, N.V.G .; Visualization, P.D., R.K.S; Supervision, N.G.S, L.C.S, N.V.G .; Project Administration, P.D, R.K.S, N.V.G.; Funding Acquisition, N.V.G. All authors have read and agreed to the published version of the manuscript.
Figure 1.
Weekly average of (a) soil temperature measured at 10 cm soil depth and (b) air temperature measured near leaf canopy, across the treatments. C: Ambient temperature, No Residue; OTC: Open Top Chamber, No Residue; R: Ambient temperature, Residue; and OTC + R: Open Top Chamber, Residue.
Figure 1.
Weekly average of (a) soil temperature measured at 10 cm soil depth and (b) air temperature measured near leaf canopy, across the treatments. C: Ambient temperature, No Residue; OTC: Open Top Chamber, No Residue; R: Ambient temperature, Residue; and OTC + R: Open Top Chamber, Residue.
Figure 2.
Average daily temperature range (DTR) measured at 10 cm soil depth across OTC, irrigation, and residue treatments. The larger square dots indicate the model-predicted average temperature for each treatment, while error bars represent 95% confidence intervals. Each smaller dot represents the temperature range (averaged by month) of an individual plot. The boxes show the lower quartiles, median, and upper quartiles, and the whiskers extend to the most extreme observations.
Figure 2.
Average daily temperature range (DTR) measured at 10 cm soil depth across OTC, irrigation, and residue treatments. The larger square dots indicate the model-predicted average temperature for each treatment, while error bars represent 95% confidence intervals. Each smaller dot represents the temperature range (averaged by month) of an individual plot. The boxes show the lower quartiles, median, and upper quartiles, and the whiskers extend to the most extreme observations.
Figure 3.
(a)Weekly average volumetric water contents (VWC) measured at 10 cm soil depth. The light blue, vertical bars show the weekly total rainfall. (b) Average VWC, across OTC, irrigation, and residue treatments. The larger square dots indicate the model-predicted average temperature, while error bars represent 95% confidence intervals. Each smaller dot represents the average temperature (averaged by months) of an individual plot. The boxes show the lower quartiles, median, and upper quartiles, and the whiskers extend to the most extreme observations. C: No Open Top Chamber, No Residue; OTC: Open Top Chamber, No Residue; R: No Open Top Chamber, Residue; and OTC + R: Open Top Chamber, Residue.
Figure 3.
(a)Weekly average volumetric water contents (VWC) measured at 10 cm soil depth. The light blue, vertical bars show the weekly total rainfall. (b) Average VWC, across OTC, irrigation, and residue treatments. The larger square dots indicate the model-predicted average temperature, while error bars represent 95% confidence intervals. Each smaller dot represents the average temperature (averaged by months) of an individual plot. The boxes show the lower quartiles, median, and upper quartiles, and the whiskers extend to the most extreme observations. C: No Open Top Chamber, No Residue; OTC: Open Top Chamber, No Residue; R: No Open Top Chamber, Residue; and OTC + R: Open Top Chamber, Residue.
Figure 4.
Average (a) soil organic matter, (b) microbial biomass carbon, and (c) soil respiration, across OTC, irrigation, and residue treatments. The larger square dots indicate the model-predicted average values, while error bars represent 95% confidence intervals. Each smaller dot represents the data for an individual plot. The boxes show the lower quartiles, median, and upper quartiles, and the whiskers extend to the most extreme observations.
Figure 4.
Average (a) soil organic matter, (b) microbial biomass carbon, and (c) soil respiration, across OTC, irrigation, and residue treatments. The larger square dots indicate the model-predicted average values, while error bars represent 95% confidence intervals. Each smaller dot represents the data for an individual plot. The boxes show the lower quartiles, median, and upper quartiles, and the whiskers extend to the most extreme observations.
Figure 5.
Regression plots showing a relationship between (a) soil organic matter and soil temperature, (b) microbial biomass carbon and soil organic matter, and (c) microbial biomass carbon and available nitrate nitrogen in soil. The solid line represents a regression line predicted from linear mixed effect model. The shaded region represents 95% confidence intervals. The light red dots show the dryland while blue dots show the irrigated fields. The square dots indicate the residue added plots whereas the round dots indicate plots without residue.
Figure 5.
Regression plots showing a relationship between (a) soil organic matter and soil temperature, (b) microbial biomass carbon and soil organic matter, and (c) microbial biomass carbon and available nitrate nitrogen in soil. The solid line represents a regression line predicted from linear mixed effect model. The shaded region represents 95% confidence intervals. The light red dots show the dryland while blue dots show the irrigated fields. The square dots indicate the residue added plots whereas the round dots indicate plots without residue.
Figure 6.
Average (a) above ground biomass, (b) belowground biomass and (c) seed cotton yield across OTC, irrigation, and residue treatments. The larger square dots indicate the model-predicted average values, while error bars represent 95% confidence intervals. Each smaller dot represents the data for an individual plot.
Figure 6.
Average (a) above ground biomass, (b) belowground biomass and (c) seed cotton yield across OTC, irrigation, and residue treatments. The larger square dots indicate the model-predicted average values, while error bars represent 95% confidence intervals. Each smaller dot represents the data for an individual plot.