Cover Crops for Resilience of a Limited-Irrigation Winter Wheat–Sorghum–Fallow Rotation: Soil Carbon, Nitrogen, and Sorghum Yield Responses

Cover crops can improve soil health by maintaining soil organic carbon (SOC) and nitrogen (N) contents, yet their dynamics in relation to crop yield in a semi-arid cropping system are poorly understood. The main objective of this study was to evaluate the response of diverse winter cover crop species and their mixture on SOC and N fractions and their relationship with sorghum (Sorghum bicolor L. Moench) yield in a winter wheat (Triticum aestivum L.)–sorghum–fallow rotation with limited irrigation management. Cover cropping treatments included pea (Pisum sativum L.), oat (Avena sativa L.), canola (Brassica napus L.), and mixtures of pea+oat (POM), pea+canola (PCM), peat+oat+canola (POCM), and a six-species mixture (SSM) of pea+oat+canola+hairy vetch (Vicia villosa Roth)+forage radish (Raphanussativus L.)+barley (Hordeum vulgare L.) as cover crops and a fallow. Soil samples were analyzed for residual inorganic N, potentially mineralizable carbon (PMC) and nitrogen (PMN), SOC, and total N. Response of labile inorganic N, PMC, and PMN varied with cover crop treatments. The SOC and total N contents did not differ among treatments but were 20% and 35% higher in 2020 than in 2019, respectively. Sorghum grain yield was 25% and 40% greater with oats than with PCM and canola cover crops in 2019, while it was 33–97% greater with fallow and oats than other treatments in 2020. Oat as a cover crop could improve the resilience of limited-irrigation cropping systems by increasing SOC, soil N, and crop yield in semi-arid regions.


Introduction
Sorghum represents the third-largest cereal grain in the United States, with 9,473,184 metric tons yield and 2,061,873 ha area in 2020 [1]. Because sorghum is a heat-and drought-tolerant crop, its acreage is increasing under low or no irrigation (dryland) conditions in semi-arid Kansas, Oklahoma, Texas, and New Mexico. However, the decline in biomass carbon (C) and nutrient inputs to regulate nutrient cycling during the transition from irrigated to limited irrigation or dryland may threaten sorghum production in the long-term and affect the sustainability of sorghum-based cropping systems [2]. Farmers are looking for crop management options that increase SOC and N storage and improve soil health while maintaining crop yield in water-limited environments.
In recent years, cover cropping has been increasingly adopted to increase SOC and N sequestration and improve soil health in semi-arid environments [3]. Additional aboveand below-ground biomass inputs through cover cropping increase both active and passive pools of SOC and N [4]. Grass cover crops with dense rooting systems produce a high C:N ratio but greater biomass, increasing soil C levels [5,6]. Legume cover crops symbiotically fix atmospheric N in their root nodules and increase soil N content [7]. They produce highquality residues and favor early mineralization and rapid recycling of soil nutrients [5,6]. Table 1. Maximum and minimum monthly air temperature, total monthly precipitation, and irrigation water applied during sorghum and cover crop growing periods from 2016-2020.

Month
Max and Min Monthly Air Temperature ( • C) Agronomy 2021, 11, x FOR PEER REVIEW 3 of 11 approximately 70% of which occurs from May through September (Table 1). The study area experiences a high seasonal and inter-annual variability in precipitation, with shortterm drought periods often occurring within crop-growing seasons. Soils are characterized as Olton clay loam (Fine, mixed, superactive, thermic Aridic Paleustolls) under the United States Department of Agriculture (USDA) soil classification system, with 43.7% sand, 21.5% silt, and 34.8% clay. The experiment was established under no-tillage management in 2015. The experimental field was previously under conventional management of irrigated corn and sorghum production for several years. Soil bulk density measured at the time of plot establishment ranged from 1.1-1.3 g cm −3 , gravimetric soil moisture from 14.2-20.7%, soil pH from 7.9-8.1, and electrical conductivity from 0.28-0.51 dS m −1 . The SOM varied from 13-16 g kg −1 , inorganic N from 2. 24-8.97 kg N ha −1 , and available phosphorus (P) from 24.9-34.5 P kg ha −1 at 0-15 cm depth [18]. Table 1. Maximum and minimum monthly air temperature, total monthly precipitation, and irrigation water applied during sorghum and cover crop growing periods from 2016-2020.  , and a six-species mixture [SSM] of pea+oat+canola+hairy vetch+forage radish+barley. Each year, experimental plots were treated with N-phosphonomethyl glycine 53.8% (glyphosate) at the rate of 0.38 L ha −1 and 2,4-D ester (6 lb gal −1 , LV-6) at the rate of 0.87 L ha −1 with ammonium sulfate and nonionic surfactant at rates of 20 g L −1 and 5 mL L −1 , respectively, two weeks before cover crop planting. Cover crops were planted in the last week of February in a fallow field using a plot drill (Great Plains 3P600, Salina, KS, USA). The seeding rates for sole cover crops were , gravimetric soil moisture , soil pH from 7.9-8.1, and electrical conductivity from 0.28-0.51 dS m −1 . from 13-16 g kg −1 , inorganic N from 2.24-8.97 kg N ha −1 , and available rom 24.9-34.5 P kg ha −1 at 0-15 cm depth [18].

Month
air temperature, total monthly precipitation, and irrigation water applied ods from 2016-2020.  Weather data recorded at the New Mexico State University Agricultural Science Center at Clovis, NM. The growing season for sorghum spanned from June-Oct in each crop year. Winter cover crops were grown from the last week of February through the last week of May (85-90 days). † Irrigation water applied to winter cover crops. ‡ Annual average air temperature, total annual precipitation, and irrigation water.
The cover crop plots (18.3 m × 12.2 m) were established in a randomized complete block design with eight treatments and three replications. Cover crops were planted in fallow periods before each winter wheat and sorghum crop in a winter wheat-sorghumfallow/cover crop rotation. Cover crop treatments included fallow (no cover crop); three sole cover crops: pea, oat, and canola; and four cover crop mixtures: pea+oat [POM], pea+canola [PCM], pea+oat+canola [POCM], and a six-species mixture [SSM] of pea+oat+canola+hairy vetch+forage radish+barley. Each year, experimental plots were treated with N-phosphonomethyl glycine 53.8% (glyphosate) at the rate of 0.38 L ha −1 and 2,4-D ester (6 lb gal −1 , LV-6) at the rate of 0.87 L ha −1 with ammonium sulfate and nonionic surfactant at rates of 20 g L −1 and 5 mL L −1 , respectively, two weeks before cover crop planting. Cover crops were planted in the last week of February in a fallow field using a plot drill (Great Plains 3P600, Salina, KS, USA). The seeding rates for sole cover crops were 22.4, 44.8, and 4.5 kg ha −1 for pea, oat, and canola, respectively. The seeding rates for hairy vetch, forage radish, and barley were 11.2, 4.48, and 44.8 kg ha −1 , respectively. Cover crop species used in two-, three-, and six-species mixtures used 50, 33, and 16.6% of the sole seeding rates. Irrigation water was applied to cover crops only for seed germination in each year (Table 1), after which no additional irrigation was applied. All cover crops were maintained in plots for three months before being chemically terminated at the flowering stage of oat (85-90 d). After termination, the cover crop residues were left on the soil surface. Fallow plots were sprayed once in March and again at the time of cover crop termination in May to control weeds.

Winter Wheat and Sorghum Management
Winter wheat (var. TAM113 2016-2018 and TAM114 2019-2020) was planted in the second week of October using a plot drill (Great Plains 3P600, Salina, KS, USA) at a seeding Agronomy 2021, 11, 762 4 of 11 rate of 62 kg ha −1 with the drill spacing maintained at 0.25 m. Details of crop management are also explained in Thapa et al. [10]. All experimental plots received 67 kg N ha −1 in 2016 and 2017, 70 kg N ha −1 in 2018 and 2019, and 12 kg sulfur (S) ha −1 each year through fertigation during the spring season. Sorghum (cultivar NK 5418) was planted in the first week of June using a no-till drill (John Deere, Moline, IL, USA) at a seeding rate of 123,553 seeds ha −1 with the row spacing maintained at 0.76 m. All sorghum plots received 97 kg N ha −1 and 15 kg S ha −1 from a mixture of urea, ammonium nitrate, and ammonium thiosulfate in liquid form at the time of planting each year. The experiment was maintained under limited-irrigation conditions. The amount of irrigation water applied to sorghum in each year is presented in Table 1. Soil fertility management was based on soil test recommendations for both wheat and sorghum. Weeds and pests were controlled by applying thifensulfuron-methyl 25% + tribenuron-methyl 25% (Affinity BS) at the labeled rate of 0.05 L ha −1 , 2,4-D ester 0.72 kg L −1 (LV-6) at the rate of 0.87 L ha −1 , insecticide chlorpyrifos 44.9% (Govern) at the rate of 0.2 L ha −1 , and Prowl H 2 O at the rate of 0.6 L ha −1 during crop growth as needed.

Soil Sampling and Laboratory Analysis
Soil samples were collected from 0-15 cm depth at the time of sorghum harvest (October in all years). The fallow plots were considered a control to compare soil N and C changes due to cover cropping. Three soil cores were collected diagonally from each plot using a core sampler (2 cm diam.), composited, and thoroughly homogenized, and all visible plant materials (roots, stems, and leaves) and crop residues were removed by hand. Soil samples were transported to the laboratory, and approximately 20 g subsamples were used for soil moisture estimation. The rest of the samples were stored in a refrigerator at 4 • C for inorganic N, PMN, and PMC estimation. About 20 g subsamples were air-dried and ground to estimate SOC and total N in 2019 and 2020. In the laboratory, gravimetric soil moisture was estimated by oven drying 20 g soil samples at field moisture at 105 • C for 24 h. Soil inorganic N concentration was determined by measuring nitrate (NO 3 − ) and ammonium (NH 4 + ) ions in an automated flow injection N analyzer (Timberline Instruments, LLC, Boulder, CO, USA). For this, 5 g of soil were extracted with 25 mL of 1M potassium chloride (KCL). The PMC concentration was estimated by aerobic incubation of approximately 20 g soil samples in 1-L glass jars for two weeks and measuring CO 2 released from incubation jars using an infrared gas analyzer (LI-COR Inc., Lincoln, NE, USA) as described in Ghimire et al. [23]. The PMN concentration was analyzed by extracting the incubated samples with KCL and determining inorganic N concentration as above without subtracting the unincubated control. The SOC and total N were analyzed by dry combustion method in a C:N analyzer (828 series, LECO corporation, St. Joseph, MO, USA). Soil bulk density was determined by collecting four undisturbed soil cores (2.3 cm diam. × 15 cm depth) using a hand probe within each plot, oven drying soil samples at 105 • C for 24 h, and dividing the weight of the oven-dried sample by the volume of the core. The contents of inorganic N, PMN, and PMC (kg ha −1 ) were determined by multiplying their concentrations (ppm or mg kg −1 ) by soil bulk density (g cm −3 ) and the thickness of the soil layer (i.e., 0.15 m) multiplied by 10.

Sorghum Grain and Total Yield
Sorghum was harvested at physiological maturity in the last week of October in all years by hand-harvesting a bundle grain sample from the 6th and 7th row at 6 m length, whereas stalks were harvested from the same rows at 1.5 m length in each plot. Sorghum aboveground biomass (head and stalk) was collected in plastic bags, brought to the laboratory, and heads were thrashed using a plot combine thresher (Wintersteiger, Ried im Innkreis, Austria) to separate the grain. The moisture percentage of sorghum grain was determined with a moisture meter (GAC 2100b, DICKEY-john Corporation, Auburn, IL, USA), whereas stalks were oven-dried at 65 • C for 72 h to determine dry weight. Sorghum grain yield was adjusted to 12% moisture. Harvest index was calculated by dividing sorghum grain yield by total yield (head + stalk).

Statistical Analysis
Data for soil N and C fractions and crop yield parameters were analyzed using a MIXED model procedure in SAS (ver. 9.4, SAS Institute Inc., Cary, NC, USA, 2013) for a completely randomized block experiment. All the data were tested and met the assumptions for normality of residuals and equality of variance. During analysis, treatment and year were considered fixed effects, whereas replication, replication × treatment, and replication × year interactions were considered random effects in the statistical model. The significance of fixed effects was tested using an F-test that utilized Type III sums of squares for balanced cases. Data were evaluated by year and treatment when there was significant treatment × year interaction. A post hoc Fisher's protected least significant difference (LSD) was used to separate means when treatment, year, and interaction effects were significant at p = 0.05, unless otherwise stated. The log transformation was applied for soil inorganic N and PMN because of the non-normal distribution of the data, and the back-transformed means were reported. Simple linear regression analysis was performed to explain the relationship of inorganic N and PMN with PMC and to predict sorghum grain yield as a function of PMN.

Results
Soil inorganic N and PMN contents were significantly different among years but not among treatments and treatment × year interaction ( Table 2). Averaged across treatments, the residual soil inorganic N content decreased over years of cover cropping. Inorganic N was 2-10 times higher in 2016 than in other years (Table 3). Soil PMN content in 2020 was 12.1 kg ha −1 , which was statistically similar to 11.3 kg ha −1 in 2018 but greater than in 2016 (8.21 kg ha −1 ), 2017 (4.88 kg ha −1 ), and 2019 (0.89 kg ha −1 ) ( Table 3). Soil PMN contents in 2016, 2017, and 2019 were also different from one another. Table 2. Analysis of variance (ANOVA) for soil inorganic nitrogen (IN), potentially mineralizable nitrogen (PMN) and carbon (PMC), soil organic carbon (SOC), total N (TN), sorghum grain yield (SGY), total yield (TY), harvest index (HI), and thousand-grain weight (1000-Gwt). Unlike soil N fractions, soil PMC content, the easily decomposable fraction of SOC, was different among treatments at p = 0.079 and years at p = 0.095, and not among treatment × year interaction ( Table 2). The SOC and total N evaluated during the fourth and fifth years of the study were different between years, but not among treatments and treatment × year interaction ( Table 2). Soil PMC was 254 kg ha −1 under SSM, followed by oat (239 kg ha −1 ), which were significantly greater than under pea (189 kg ha −1 ). Fallow (average 232 kg ha −1 ) and canola (216 kg ha −1 ) remained intermediate of SSM, oat, and pea ( Table 3). Irrespective of treatments, PMC was 46% and 36% higher in 2016 and 2018, respectively, than in 2017. Soil PMC contents in 2019 and 2020 remained in between 2016, 2017, and 2018, and were not statistically different from each other. The SOC and total N contents were 20% and 35% higher in 2020 than in 2019, respectively. Table 3. Means of treatment and year on soil inorganic nitrogen (IN), potentially mineralizable nitrogen (PMN) and carbon (PMC), soil organic carbon (SOC), total N (TN), sorghum grain yield (SGY), total yield (TY), harvest index (HI), and thousand-grain weight (1000-Gwt).

Treatment IN
d, 21.5% silt, and 34.8% clay. The experiment was established under no-tillage manement in 2015. The experimental field was previously under conventional management irrigated corn and sorghum production for several years. Soil bulk density measured the time of plot establishment ranged from 1.1-1.3 g cm −3 , gravimetric soil moisture m 14.2-20.7%, soil pH from 7.9-8.1, and electrical conductivity from 0.28-0.51 dS m −1 . e SOM varied from 13-16 g kg −1 , inorganic N from 2.24-8.97 kg N ha −1 , and available osphorus (P) from 24.9-34.5 P kg ha −1 at 0-15 cm depth [18]. mum monthly air temperature, total monthly precipitation, and irrigation water applied p growing periods from 2016-2020. Each year, experimental plots were ated with N-phosphonomethyl glycine 53.8% (glyphosate) at the rate of 0.38 L ha −1 and -D ester (6 lb gal −1 , LV-6) at the rate of 0.87 L ha −1 with ammonium sulfate and nonionic rfactant at rates of 20 g L −1 and 5 mL L −1 , respectively, two weeks before cover crop nting. Cover crops were planted in the last week of February in a fallow field using a t drill (Great Plains 3P600, Salina, KS, USA). The seeding rates for sole cover crops were PMN sand, 21.5% silt, and 34.8% clay. The experiment was established under no-tillage management in 2015. The experimental field was previously under conventional management of irrigated corn and sorghum production for several years. Soil bulk density measured at the time of plot establishment ranged from 1.1-1.3 g cm −3 , gravimetric soil moisture from 14.2-20.7%, soil pH from 7.9-8.1, and electrical conductivity from 0.28-0.51 dS m −1 . The SOM varied from 13-16 g kg −1 , inorganic N from 2.24-8.97 kg N ha −1 , and available phosphorus (P) from 24.9-34.5 P kg ha −1 at 0-15 cm depth [18]. um and minimum monthly air temperature, total monthly precipitation, and irrigation water applied and cover crop growing periods from 2016-2020.    , gravimetric soil moisture , soil pH from 7.9-8.1, and electrical conductivity from 0.28-0.51 dS m −1 . from 13-16 g kg −1 , inorganic N from 2.24-8.97 kg N ha −1 , and available rom 24.9-34.5 P kg ha −1 at 0-15 cm depth [18].
air temperature, total monthly precipitation, and irrigation water applied ods from 2016-2020.   The crop yield parameters, such as sorghum grain yield, total yield, harvest index (HI), and thousand-grain weight (1000-Gwt), did not significantly vary among treatments, while they did vary among sampling years (Table 2). Sorghum grain yield varied with treatment × year interaction, whereas treatment × year effect was not observed for total yield (Table 4). When averaged across treatments, sorghum grain yield was 47% and 27% higher in 2019 than in 2018 and 2020, respectively (Table 3). It was also 31% higher in 2016 compared to 2018. Sorghum grain yield in 2017 was not different from all other sampling years. Sorghum grain yield was consistent among treatments in 2016, 2017, and 2018 (Table 4). In 2019, sorghum grain yield was the greatest under oats (9708 kg ha −1 ), which was statistically similar to fallow (8867 kg ha −1 ) and pea (8970 kg ha −1 ), but greater than under canola (average 6939 kg ha −1 ) ( Table 4). Treatments POM, POCM, and SSM remained in between oat and canola and were not different from each other. In 2020, sorghum grain yield was the greatest under fallow (9853 kg ha −1 ), which was statistically similar to oat (8340 kg ha −1 ) but greater than all other treatments. Treatments such as pea, canola, POM, PCM, POCM, and SSM were not significantly different from each other. Within a treatment, sorghum grain yield was not different among sampling years for canola and PCM. Sorghum grain yield was lower in 2018 than in 2019 for fallow, pea, oat, POM, POCM, and SSM.

Total Monthly Precipitation (mm) Irrigation Water Applied to Sorghum and Cover
The total yield was 54-72% higher in 2020 than all other sampling years (Table 3). No significant differences in total yield were observed between 2016, 2017, 2018, and 2019. The HI was 0.49 in 2019, followed by 2016 and 2017 (0.45 each), but was greater than in 2018 (0.37) and 2020 (0.25). The 1000-Gwt was also highest in 2019 (25.8 g), which was similar to 2020 (25.2 g) but greater than all other sampling years, leaving 2018 the lowest (23.1 g). There was no difference in 1000-Gwt between 2016 (24.0 g) and 2017 (24.1 g).
Simple linear regression analysis showed that inorganic N and PMN increased with PMC content (Figure 1). Also, the linear relationship between residual PMN and sorghum grain yield was negative (Figure 2). It appears that greater amounts of residual PMN retained in the soil influence the availability of N to plants, which in turn reduces sorghum Agronomy 2021, 11, 762 7 of 11 grain yield. The relationship between sorghum grain and total yield with inorganic N and PMC was not significant.

Discussion
Soils with high SOC storage are often considered resilient because they can hold more water, resist erosion, and recycle nutrients. Diversifying cropping systems by integrating cover crops can improve soil resilience by increasing SOC accumulation and sustaining crop production. Our observation of higher PMC under SSM followed by oat and the lowest under pea across all study years suggests a continuous food supply (providing a C source) to microbes under diverse cover crop mixes, supporting higher microbial activity.

Discussion
Soils with high SOC storage are often considered resilient because they can hold more water, resist erosion, and recycle nutrients. Diversifying cropping systems by integrating cover crops can improve soil resilience by increasing SOC accumulation and sustaining crop production. Our observation of higher PMC under SSM followed by oat and the lowest under pea across all study years suggests a continuous food supply (providing a C source) to microbes under diverse cover crop mixes, supporting higher microbial activity. The highest biomass produced by oat may have also increased SOC accrual. Our previous study at the same site showed higher enzyme activity associated with C and nutrient cycling under SSM and larger microbial community size and fungal community under oats and SSM [10]. From the fourth to fifth years of cover cropping, the SOC and total N contents increased by 20% and 35%, respectively, suggesting gradual improvement in soil health over the years. Adoption of no-tillage could have complementary effects in increasing SOC. Although there were no significant differences in SOC and total N among treatments, we observed numerically greater SOC and total N under cover crops than under fallow, indicating that a small but positive change in SOC accumulation is underway with cover crops.
The inorganic N and PMN increased linearly with an increase in PMC (Figure 1). The C and N mineralization occurs simultaneously during microbial decomposition of organic residues. Increased microbial activity, as indicated by higher PMC, supported N mineralization. However, we did not observe significant differences among treatments on residual inorganic N and PMN at the time of sorghum harvest. Cover crops were terminated in the last week of May, and residues were left on the soil surface. During the summer, the presence of cover crop residues likely improved precipitation storage efficiency, increased infiltration, reduced evaporative loss, and helped conserve more water in the soil profile. Summer soil temperatures can be very high (up to 44 • C in July) in eastern New Mexico. Increased soil moisture and temperature may have accelerated the mineralization of cover crop biomass N, which may have been utilized by sorghum and supported yield formation. The residual N was measured five months after cover crop termination at the time of sorghum harvest. A previous study at the same site reported 41-49% less inorganic N under cover crops than under fallow at cover crop termination, while no difference among treatments was observed at the time of wheat harvest [18]. A significant negative linear relationship between residual PMN and sorghum grain yield suggests that yield depends on crop N uptake. Sorghum yield was higher in the years where residual inorganic N and PMN were lower, and vice-versa (Table 3). Residual inorganic N content was significantly higher in 2016 than all other years, irrespective of cover crop treatments. The field was previously under conventional management of irrigated corn and sorghum for several years before the experiment was established in October 2015 with no-tillage management, cover cropping, and about 50% of typical irrigated crop production. Reduced water and nutrient inputs from plant biomass recycling gradually depleted N in the soil profile, ultimately reducing yield and SOC storage.
Water and nutrient availability are the major limiting factors for crop production in semi-arid regions. When averaged across treatments, sorghum grain yield varied among years (Table 3). This could be attributed to a variable amount of water available from irrigation and precipitation during the crop-growing season. The amounts of precipitation plus irrigation received during sorghum growth periods for 2016 through 2020 were 411, 563, 510, 590, and 411 mm, respectively (Table 1). Stone and Schlegel [24] reported that the growing season irrigation explained 63% of variations in sorghum yield. The effect of cover crops on sorghum yield was not observed until the fourth year of the study. Cover crop residue after termination, and main crop residue after harvest, was left on the soil surface. Annual cover crop biomass input was <1 Mg ha −1 for these study plots [10], which was not enough to maintain SOC in our study. At least 5 Mg ha −1 of cover crop residue [23] or more (6.8 to 7.6 Mg ha −1 yr −1 ) [25] is needed to maintain SOC stocks under limited irrigation in the semi-arid Southern and Central High Plains of USA. However, leaving sorghum residue on the ground and years of no-tillage helped to maintain SOC across all treatments. No-tillage was in place across all treatments since the establishment of this study in a field previously under conventional tillage for several years.
Observations from interaction effects showed that sorghum yield was highest under oat in 2019. This could be associated with greater biomass input and surface residue cover with oat. The increased enmeshing action of oat roots may have also improved soil aggregation and infiltration, leading to more soil water available for the subsequent sorghum crop. Oat produced greater biomass (3343 kg ha −1 ) than PCM (2010 kg ha −1 ) and canola (2407 kg ha −1 ) in 2019. The oat residues, with higher C:N ratios than pea and canola residues [23], might have slowed down the decomposition process and provided soil cover for a longer period of time, thus improving soil moisture retention, sustaining nutrient cycling, and improving sorghum yield. Sorghum yield in 2020 was highest under fallow, which was similar to oat but significantly greater than all other treatments (Table 4). Fallow treatment could have preserved soil water and nutrients and consequently increased yield, while improved SOC and total N under oats with approximately half of the sorghum water requirement applied [18] was able to maintain yield as under fallow.
Total yield (grain + biomass) also did not vary among treatments, but it was greater in 2020 than other sampling years. The amount of irrigation water applied during the sorghum period was relatively higher in 2020 than in previous years because of extreme dry conditions (Table 1). Despite higher biomass production, grain yield did not increase in 2020, possibly due to water stress during the grain filling stage. Most of the irrigations were applied early in the growing season. However, higher 1000-grain weight for sorghum in 2019 and 2020 also supported the gradual improvements in soil quality through cover cropping. The gradual release of soil nutrients might have improved nutrient use efficiency and led to greater grain dry matter, resulting in greater 1000-grain weight.

Conclusions
Increased SOC storage is a key strategy to improve soil health and resilience in semiarid cropping systems. The results of this study suggested that short-season, spring-planted cover crops could increase SOC storage with no or low yield penalty to grain sorghum in limited-irrigation conditions. The numerically greater SOC and total N contents with cover crops compared to fallow suggest positive change on SOC sequestration is underway with cover crops. The SSM had the greatest PMC, indicating that a mixture of grasses, legumes, and brassicas as cover crops supports higher microbial activity. Oat as a cover crop increased sorghum yield by 8-40% in a good rainfall year. Even in a drier year, oat was able to maintain yield as under fallow. Among cover crops, oat increased sorghum yield by 33-67% in the drier year. However, the effect of cover crops on sorghum yield was not observed until the fourth year of the study. Oat as a cover crop increased sorghum yield by 8-40% when precipitation plus irrigation was 590 mm, and maintained yield as under fallow when precipitation plus irrigation was 411 mm. Oat and SSM as cover crops could be used to improve soil health and resilience in limited-irrigation cropping systems in the semi-arid SHP and similar agroecosystems worldwide.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.