Abstract
This study combines geospatial data and a classification scheme that uses landform elements to derive landform complexes that codify the collection of soils data at variable scale within a single field site. Our experiment was initiated in 2018 on three, 1.6 ha self-contained watersheds representing a southern tall grass prairie (STGP), a system of continuous winter wheat (Triticum aestivum) tilled via offset disking and chisel plow (WWCT), and a minimally disturbed winter wheat system that was periodically planted to a warm season forage, typically sorghum-sudangrass (Sorghum bicolor L.) (WWMT) from 1978 to 2018. A class I soil survey was conducted in 2018 by grid sampling the landscape of all watershed systems at the site. The survey indicated four distinct catena were present across all watersheds, which enabled us to utilize a split block design. This statistical approach allowed for testing of interactions among management practices, landscape position, and soil depth to obtain means and standard errors for different edaphic properties. Using a hydraulic probe, 144 random soil cores were collected to a 30 cm depth at each of the four 4.6 m by 3.8 m replicated blocks per landscape position (tread, riser, and toe) within the three watersheds. Cores were further divided into three depths (0 to 5, 5 to 15, and 15 to 30 cm). Baseline analyses included Mehlich-3, soil sulfate (SO4) and DTPA-sorbitol extractions, soil texture, bulk density, pH, total soil organic carbon (TSOC) and total soil nitrogen (TSN), particulate organic matter (POM), and non-hydrolysable C (RCAH), the resistant fraction of soil organic C. The majority of edaphic properties associated with soil classification varied with landscape position and depth. These included clay content, base saturation (calcium [Ca], magnesium [Mg], and potassium [K]), pH, and sulfur in the form of sulfate (SO4-S) and phosphorus (P). Carbon and N fractions varied with land use, conservation practices, and/or depth. The establishment of replicate sampling stations that account for and limit the spatial variability of edaphic properties within defined landform complexes enables researchers to more accurately quantify the effects of conservation practices and land management.
The creation of replicate sampling stations that represent landform complexes within uniform zones of management can aid researchers in addressing variations in scale and the inherent spatial variability associated with field research that limits comparisons among management practices. Both crop and animal agriculture in the Southern Great Plains are subject to high degrees of variability in both climate and natural resources (Steiner et al. 2018). Management practices that conserve soil and water can aid in mitigation of and adaptation to climate variability and change (Delgado et al. 2011; Lal et al. 2011). The National Institute of Food and Agriculture has defined climate change as “a trend of significant fluctuations in major weather patterns and conditions brought about by increased greenhouse gases. The trend has seen significant variations in temperature, rain, wind, and other factors over a period of years” (USDA NIFA 2020). Research is being conducted on eight 1.4 ha sized watersheds (the Water Resources Erosion [WRE] site) located at the USDA Agricultural Research Service (ARS) Grazinglands Research Laboratory in El Reno, Oklahoma, to address “agricultural sustainability, climate change, ecosystem services, and conservation of natural resources at the watershed or landscape scale”—a primary objective of the USDA ARS Conservation Effects Assessment Project (CEAP) (Mausbach and Dedrick 2004; Duriancik et al. 2008) and the Long-Term Agro-Ecosystem Research (LTAR) network (Spiegal et al. 2018; Kleinman et al. 2019; Nelson et al. 2020).
Globally and regionally, retention and distribution of organic material and nutrients in agroecosystems vary due to landscape-scale management and terrain attributes (elevation, slope, aspect, and texture), and anthropogenic management (Liu et al. 2006; Viaud et al. 2010; Zhang et al. 2011; Blanchet et al. 2017). These observations indicate that topographic and soil controls strongly influence nutrient distributions, losses, and uptake patterns within fields, affecting soil and water quality at microsite and landscape scales. Adoption of precision agriculture and grid sampling techniques to assess and adjust for the effects of terrain attributes on yield became feasible in the 1980s and 1990s, when computational capabilities of computers increased and global positioning system (GPS) applications became feasible (Stafford 2000). Once these tools were more widely available, proponents of soil and water conservation recommended their use to manage natural resources, and to adapt to climatic variation and change (Delgado et al. 2011).
Conservation management and terrain attributes partially determine mineralization-immobilization and turnover of nutrients and their potential loss. Additionally, mixtures of residues in polycultures, such as the southern tall grass prairies (STGP) included in this experiment (a USDA Natural Resources Conservation Service [NRCS] recommended conservation practice standard, code 340), may act as reserves for excess nutrients preventing losses (Matson et al. 1997) and aiding in the timing and release of plant available N with crop uptake (Lawson et al. 2013). However, reduced tillage and year-round (living) plant cover stratifies nutrients and increases field heterogeneity, reducing our ability to discern differences in land management (Bauer et al. 2002).
Nutrient concentrations and terrain attributes of soils are integral components of soil quality and health (used interchangeably in the current discussion). Soil quality is evaluated using inherent and dynamic soil properties that serve as indicators of soil function. These include sustaining biological diversity, activity, and productivity regulating water and solute flow; filtering, buffering, degrading organic and inorganic materials; storing and cycling nutrients and carbon (C); and providing physical stability and support (Larson and Pierce 1991). Inherent or use-invariant, soil properties such as aspect, slope, soil texture, type of clay, depth to bedrock, and drainage class change little or not at all with management forming over time (Doran et al. 1996). Management and natural disturbances that typically occur among growing seasons and decades affect dynamic soil properties. Indicators of soil quality are divided into physical, chemical, and biological categories relating to the soil function affected. Each category represents both inherent and dynamic indicators (Doran et al. 1996). Implementation of soil and water conservation practices in cropping and grassland systems were shown to reduce erosion while sustaining soil quality through the retention of residues and soil organic matter (Sanderson et al. 2013).
High-intensity soil survey methods are used when “detailed and very precise knowledge of soils and their variability’ are required (USDA NRCS 2019). Despite the effectiveness of such surveys, greater detail with respect to soil properties is required for farm and range management. During the last several decades geostatistics and precision agriculture have become tools used to assess and interpret data related to soil quality in the context of spatial variability and conservation practices (Allen et al. 2010; Killham and Staddon 2002; Lerch et al. 2005).
Seasonal variations in plant available nitrogen (N) and C turnover can be estimated via measurements of particulate organic matter C (POMC) and N that contain significant proportions of plant material (Cambardella and Elliott 1993). Estimates of nonhydrolysable C (RACH) turnover rates are in the 100s to 1,000s of years’ timeframe (Paul et al. 2001a). Despite lengthy turnover rates, the amount of C contained in RCAH can increase in years or decades with the addition of compost, animal wastes, or other materials containing recalcitrant compounds (Fortuna et al. 2003b). The POM fraction typically contains material with biologically active C and N and other materials with intermediate turnover rates (Fortuna et al. 2018). The pool of resistant C, defined in this paper as RCAH, can serve as a stable measure or marker that remains constant over time. In contrast, soil organic C (SOC) contains both biologically active and inactive components that temporarily mask the effects of conservation management on C storage and cycling (Stockmann et al. 2013). Particulate organic matter and RCAH provide information with respect to the short- and long-term effects of land use and conservation management on soil C and N processes.
Our research focuses on identifying and testing diversified, adaptive systems of crop and livestock production that reduce disturbance and improve soil and water quality while enhancing productivity. Previous research on the WRE determined that landscape position and applied management had significant effects on the concentration and distribution of soil organic matter (SOM) (Northup and Daniel 2010). Other researchers in the Plains region have also determined that terrain and land management affect the distribution and cycling of nutrients. Specifically, toe slopes had higher total soil C, N, phosphorus (P), and mineralization rates across a range of landscapes in Colorado (Burke et al. 1995). The same study found that treatments involving no-till or perennial grasses immobilized N and resulted in higher rates of C mineralization than tilled treatments on sites and landscape positions with high productivity and organic matter concentrations after five years of management.
An objective of our study is to combine geospatial data and a landform element classification scheme to derive landform complexes that codify the collection of soil and water data at variable scale within a single field experiment. An additional objective is to establish replicate sampling stations within defined landform complexes to explain and limit spatial variability that confounds researchers’ interpretation of the effects of soil conservation practices and land use on edaphic properties that also serve as soil quality measurements. Indicators in this experiment include terrain attributes (elevation, slope, aspect, and texture) and measurements obtained from routine soil testing analysis. We also contrast the utilization of current soil quality measures, total soil organic C (TSOC) and total soil N (TSN), with alternate measures of C fractions used to assess conservation management and land use. Our long-term goal is to derive a set of flexible real-time guidelines that growers can use in on site management of a field or portion of a field, and to determine the optimal time required to improve soil quality in agroecosystems under conservation management practices (reduced tillage, forage-based cover crops, and rotational grazing at low stocking rates).
Materials and Methods
Site Description. This study was undertaken within the WRE watersheds located on the USDA ARS Grazinglands Research Laboratory, near El Reno, Oklahoma (35°33′29″ N, 98°1′50″ W). The watersheds were constructed in 1976 and included eight experimental paddocks, with four managed as STGP and four managed as cropland. The site was STGP prior to the WRE initiative, with the larger land unit containing the paddocks grazed by groups of cow/calf pairs. Manmade berms, ridges, and natural boundaries that create self-contained watersheds surround each unit. Each unit is 80 m in width, 200 m in length, and encloses 1.6 ha. Slopes of the watersheds are oriented east to west and range between 3% and 4%. The plant community of all units in STGP for the site are dominated by big bluestem (Andropogon gerardii Vitman), indiangrass (Sorghastrum nutans [L.] Nash), and little bluestem (Schizachyrium scoparium [Michx.] Nash), with a mixture of other tall, short, and mid-grass species combined with a range of broadleaf species (forbs). Depending on growing conditions and management, the three dominant species generate 60% (±10%) of total biomass produced by native paddocks. All watersheds maintained in STGP were remnant prairies that were never tilled or planted to agronomic crops.
Three of the eight paired watersheds representing each land use and management were selected for this experiment. Watersheds 2 and 4 were to be converted from grassland to conservation management in the fall of 2018. Therefore, they were not included in the current experiment. Watersheds 1 and 5 were chosen in part because the spatial variability of SOM had been previously mapped by Northup and Daniel (2010) in 2004. The third watershed unit 8 was chosen from the two remaining watersheds maintained as continuously tilled winter wheat (Triticum aestivum). Watershed unit 8 was selected because it contained the least soil variability and more consistent management relative to unit 6. In the current study, watershed unit 1 is maintained as STGP. Watershed unit 5 (WWMT) is under low-disturbance management using no or minimal tillage and planted to winter wheat since 1976. Watershed 5 was disturbed in 1986 only by means of offset disking twice and harrowing. In years prior to 2000, watershed 5 was also planted to summer forage crops that included forage sorghum (Sorghum bicolor L.) and Sudan grass (Sorghum sudanense) in rotation with winter wheat. Lespedeza (Lespedeza cuneate) was planted as a summer forage in 1999 only. Watershed unit 8 (WWCT) has been managed using a combination of offset disking and chisel plowing and planted to winter wheat since 1976 (figure 1). The above truncated description of management applied to the watersheds was derived from Nelson et al. (2020) and Vogel et al. (2000). A detailed description of all management practices including fertilizer applications is available at https://doi.org/10.15482/USDA.ADC/1518421. Fertilization of the WWMT treatment (unit 5) was the same as WWCT (unit 8) apart from the equivalent of a 38 kg ha–1 N application to the summer forage crops.
The eight Water Resources and Erosion (WRE) watersheds located at the USDA Agricultural Research Service (ARS) Grazinglands Research Lab in El Reno, Oklahoma, United States. This paired watershed experiment has been in continuous management since 1976.
Average annual maximum and minimum air temperature is 22°C and 9°C with ranges as high as 46°C and low as –26°C. Average annual precipitation is 838 mm with highs of 1,372 mm and lows of 381 mm. The average number of days with precipitation per year is 68. Average monthly precipitation ranges from 25.4 to 127 mm and has a bimodal distribution falling primarily in May to June and again in September and October (Mesonet, El Reno, Oklahoma; https://www.mesonet.org).
Grid and Zone Sampling of Soil. Our sampling scheme was based on a class I soil survey from 2018 conducted in partnership with NRCS. The survey included all eight WRE paired watersheds. Variations in SOM measurements taken by Northup and Daniel (2010) in 2004 from watershed unit 1 in STGP management and watershed unit 5 in low-disturbance management cropped to winter wheat and periodically planted to summer forages were examined prior to implementing a sampling design. This previous research on the WRE confirmed that responses of SOM to long-term management were confounded by terrain attributes within individual paddocks and across paddocks (Northup and Daniel 2010). Their findings indicated that SOM remained constant in ~4.5 m area patches that varied due to the effects of soil catenas located at the tread, riser, and toe. Additionally, some variability could be removed by excluding the ~15 m area boarding the berms surrounding each watershed. Watershed treatments on the WRE were grid sampled in 2018 outside of the 15 m boarder described above at four 4.6 m by 3.8 m replicated blocks per landscape position—a sampling density that reflected the spatial variability measured by Northup and Daniel (2010). Baseline samples of soils were taken in 2018 and analyzed for TSOC, TSN, pH, and plant nutrients. Once we determined that baseline measurements varied within the same sampling density, each watershed treatment was divided into replicate sampling stations that represent landform complexes within uniform zones of management.
The majority of soils located within the WRE are silt loams, derived from undifferentiated Dog Creek Shale and Blaine Formations of the El Reno Group (Daniel 2001). For this study, soil samples were obtained from watershed units 1 (STGP), 5 (WWMT), and 8 (WWCT). Each watershed was grid sampled at the tread, riser, and toe slopes to create four replicate sampling zones within each landscape position by watershed, capturing the catena effects, slopes, and aspects. Using a hydraulic probe, 144 soil cores were collected from three of the watersheds at tread, riser, and toe slopes. Each watershed by slope contained four, 4.6 m by 3.8 m blocks from which four random soil cores 0 to 30 cm were taken. A Garmin ETrex 10 Outdoor Handheld GPS Navigation Unit-AW16 (Garmin Ltd., Olathe, Kansas) was used to record the locations of the corners of each of the sampling zones at each of the three landscape positions on the three watersheds (Fortuna et al. 2019). Creation of the four replicated blocks within a landscape position (tread, riser, and toe) and watershed unit enabled us to define the landform complexes needed to create uniform zones within managed paddocks that represented the variable characteristics of soils that affect production within and across watersheds. A Giddings (Windsor, Colorado) Probe, Model HDGSRPS, fitted with poly-ethylene terephthalate glycol (PTEG) soil tube liners (inert to chemicals and contaminants) was used to collect soil samples. Each soil core was capped at both ends and immediately frozen at −18°C.
The volume of the core was recorded. Whole soil cores were weighed prior to freezing. Weights corrected for gravimetric moisture and core volume were used to calculate bulk density to 30 cm depth via the core method. Frozen soil cores were cut into three depths (0 to 5, 5 to 15, and 15 to 30 cm) using a miter saw. The volume of each depth increment was recorded as was the field moist weight. Weights corrected for gravimetric moisture and core volume per depth were used to calculate a bulk density at 0 to 5, 5 to 15, and 15 to 30 cm depths (Blake and Hartge 1986). Soil samples were thawed, thoroughly mixed, and passed through a 2 mm sieve. Sieving removed plant residue to allow for further mixing and homogenized soil samples. Each watershed contained four replicate soil samples per landscape position by depth from which edaphic properties were measured.
Soil Analysis. Edaphic indicators chosen for this experiment varied within the same spatial scale of a given watershed management and landscape position. Gravimetric water content was determined for each of the original 144 soil samples (Gee and Bauder 1986). The texture of each soil sample was determined by the soil hydrometer method (Huluka and Miller 2014). Subsamples of each soil were dried to a constant weight at 60°C in a Fisherbrand Isotemp heating and drying oven (Fisher Scientific, Waltham, Massachusetts). After oven drying, each of the 144 soil samples were further subsampled for POM and nonhydrolysable C, i.e., residual C by acid hydrolysis (RCAH) (Paul et al. 2001b). Several steps are required to extract POM, defined as the organic C and N associated with sand-sized soil constituents remaining on a 53 μm sieve after soil is passed through a 2 mm sieve and dispersed with hexametaphosphate (Gregorich et al. 2008). Rinsed POM materials were transferred into aluminum weighing boats and oven dried at 60°C for 24 hours. Dried POM samples were ground with a mortar and pestle to pass a 250 mm sieve.
The resistant C fraction defined as the RCAH fraction constituting ~50% to 70% of TSOC is not biologically active but contributes to the capacity of soil to hold water, soil structure, aggregation, and other physical properties of soil (Paul et al. 2001b). RCAH in soil was isolated using oven-dried subsamples weighing 1 g. Each 1 g soil sample was added to a 50 mL digestion tube in which 25 mL of 6 mol L−1 HCl was added. Tubes were then vortexed briefly for 30 seconds. Tubes were added to a DigiPREP HT Digestion block (SCP SCIENCE, Champlain, New York), a manifold was placed on top of the tubes, and the block was maintained at 115°C for 18 hours. Thereafter, the tubes were cooled to room temperature and centrifuged for 5 minutes. The majority of the acid was decanted off without loss of the remaining soil material, and 50 mL of deionized water (diH2O) was added to each tube. Each digestion tube was then vortexed, allowed to settle for 24 hours, centrifuged, and the liquid decanted. This process was repeated three times. Remaining soils in samples were then transferred to aluminum weigh tins and oven dried at 60°C for 24 hours, cooled, and weighed immediately.
Oven-dried whole soil, POM, and RCAH were weighed and ground using a mortar and pestle (<0.25 mm). Half-gram samples of oven-dried soil, POM, and RCAH were then weighed for TSN, SOC, POMN, POMC, and RCAH carbon analysis on a 928 Series LECO C/N Analyzer (Leco Corporation, St. Joseph, Michigan) via dry combustion.
Additional baseline measurements related to soil fertility were analyzed at The Soil, Water, and Forage Analytical Laboratory located at Oklahoma State University in Stillwater, Oklahoma. Analyses included buffered and unbuffered pH (Sims 1996; Sikora 2006), Mehlich-3 extraction (Mehlich 1984), DTPA-sorbitol extraction of micronutrients (Gavlak et al. 2003), and soil sulfate (SO4-S) extraction with 0.008 M calcium phosphate [Ca3(PO4)2]. Mehlich 3 extracted P, calcium (Ca), magnesium (Mg), and potassium (K) and SO4-S were analyzed on an inductively coupled plasma atomic emission spectrometer (ICP) (Soltanpour et al. 1996). Micronutrients (zinc [Zn], copper [Cu], boron [B], and iron [Fe]) were analyzed via ICP using the method of Gavlak et al. (2003).
Statistical Analyses. The use of replicate zones within the landscape of the larger managed watershed treatments with variable soil characteristics allowed for the creation of uniform zones where responses to management could be more accurately defined. Means, standard deviations, and coefficients of variation (CV) were calculated separately for all measured variables from each of the watersheds (three) by landscape positions (three) by replicate blocks (four).
In addition, the split block design applied in this study allowed us to test for interactions among watershed treatments (three) by landscape positions (three) by replicate blocks (four) by soil depth (three). Data were analyzed as repeated measures in mixed models analyses, with three watershed treatments (management) as main effects, and three landscape positions within paddocks as split-blocks. Soil depth (n = 3) was included in the model as a repeated element and was tested with a compound symmetry (CS) covariance matrix to account for autocorrelation and covariance among depths (Littell et al. 1998). Replicate by landscape position interactions were included in the repeated statement as local subjects to order analyses. Fishers least significant difference (LSD) was applied to test significance of main and interaction effects; level of significance for tests was p = 0.05.
Nonmetric multidimensional scaling (NMDS) was also employed to determine the similarities and differences in a suite of chemical, biological, and physical measures of edaphic properties among treatments and landscape positions at a given depth. “NMDS represents the original position of data in multidimensional space using a reduced number of dimensions that can be plotted and visualized” (Lefcheck 2020). All soil variables were log transformed prior to NMDS analyses.
Results and Discussion
Soil Texture. Variations in scale between soil properties and inherent spatial variability were addressed by the derivation of landform complexes located within uniform zones of management using soil survey and grid-based sampling techniques. The class I soil survey conducted on the WRE confirmed that Bethany (fine, mixed, super active, thermic Pachic Paleustolls) and Renfrow (fine, mixed, super active, thermic Udertic Paleustolls) were the two primary soil series mapped on the watersheds, with smaller interfingering areas of Norge (fine-silty, mixed, active, thermic Udic Paleustolls) present (figure 2). There were no carbonates in any of the soil samples collected to a 30 cm depth. A pH of 6.0 across all three watershed treatments and landscape positions verified the absence of pedogenic C. Silt loam was the predominant soil texture in all three watersheds (STGP, WWMT, and WWCT) and landscape positions (tread, riser, and toe slopes) at the 0 to 5 cm and 5 to 15 cm depths (supplementary tables S1, S2, and S3). Units 1 (STGP) and 8 (WWCT) contained more clay, particularly at the 15 to 30 cm depth.
Map derived from the class I soil survey conducted on the Water Resources Erosion Units (WRE) in November of 2018.
The greatest variation in soil series occurred from north to south, and was more prominent in the toe slope position (supplementary tables S1, S2, and S3 and figure 2). However, there was a change in soil series running east to west between the toe and tread slopes, which was noted in watershed units 6 through 8 (WWCT) (figure 2). Despite variations in soil series, each soil type contained comparable mineralogy (illite, montmorillonite, kaolinite, and quartz), and was typically found in similar locations within the landscape of the region (Voss et al. 1978). Additional variability in soil texture was likely due to the interactions among soil series, occurrence of small inclusions of different soils (Williams et al. 1987), land management (grassland versus agronomic), conservation management (vegetative cover, level of disturbance), and erosion. Previous studies on the WRE by Ahuja et al. (1988) and Garbrecht et al. (2015) support the above findings. Other researchers have noted similar patterns where clay content varied due in part to erosion associated with land use (Moges et al. 2013). Golchin and Asgari (2008) found a relationship between increased clay content in cultivated versus uncultivated soils.
Spatial Variability of Edaphic Properties. The spatial variability of all edaphic measurements are presented as means, standard deviations, and coefficients of variation (CV). The spatial variability of all nutrients by depth, averaged across watershed and landscape position, are available in table 1. Coefficients of variation for all nutrients except P, Zn, Fe, and Cu ranged between 5% and 28% across the three watershed treatments and landscape positions that constituted a 4.8 ha area. Thus, the choice of areas to block by landscape position and watershed reduced the variability of most edaphic properties. The interactive effects of watershed treatment by depth on the spatial variability of TSOC, TSN, P, SO4, Ca, and Mg are contained in table 2. The spatial variability of TSOC, TSN, and Ca concentrations were reduced by approximately half when watershed treatment by depth was taken into account. The interactive effects of watershed treatment by depth on the spatial variability of soil K, Fe, Cu, Zn, and B averaged across treatment by depth are found in table 3. The spatial variability of nutrients contained in table 3 were similar to the variability in table 1 in which nutrients were averaged across watershed treatment and landscape position. Despite variable landscape position, management, and soil series within the WRE watersheds, Ca, Mg, TSOC, and TSN (tables 1, 2, and 3) had CV values between 5.11% and 17.4% within the 0 to 5, 5 to 15, and 15 to 30 cm depths. Therefore, further separation of nutrients by landscape position did little to reduce the variability of nutrients (supplementary tables S4, S5, S6, S7, S8, and S9).
Spatial variability of all nutrients by depth, averaged across watershed and landscape position as measured by coefficients of variation (CV), means (Average), and standard deviations (SD).
The interactive effects of watershed by depth averaged across landscape position on the spatial variability of total soil organic carbon (TSOC), total soil nitrogen (TSN), soil phosphorus (P), sulfate (SO4), calcium (Ca), and magnesium (Mg) as measured by coefficients of variation (CV), means (Average), and standard deviations (SD).
The interactive effects of watershed treatment by depth averaged across landscape position on the spatial variability of soil potassium (K), soil iron (Fe), soil copper (Cu), zinc (Zn), and boron (B) across treatment by depth as measured by coefficients of variation (CV), means (Average), and standard deviations (SD).
Macro and Secondary Nutrients. The majority of measured variables associated with soil classification and terrain attributes on the WRE varied with depth and landscape position. These variables included clay content, macro (P and K), secondary (cations included in base saturation [Ca, Mg, K] and SO4-S), and micronutrients (Fe, Zn, Cu, and B). Phosphorus, K, and SO4-S concentrations were highest in the WWCT (Unit 8) at 0 to 5 cm depth (table 4). The WRE units were not fertilized with P or K for a decade prior to soil sampling in 2018, and sulfur in the form of fertilizer was never applied to the WRE.
The interactive effects of watershed treatment by depth averaged across landscape position on concentrations of macro and secondary nutrients in soil.
Soil texture is an integral component of soil fertility and quality assessment that has a direct effect on concentrations of soil macro, secondary, and micronutrients, physical parameters, and water availability (Fine et al. 2017; Hurisso et al. 2018). Land grant universities have conducted nutrient testing for recommendations related to soil fertility on a regional scale for ~80 years, but only recently at a national level (Lyons et al. 2020). Regional and local variations in soil types, crops, and climate often require different fertility recommendations, and in some instances protocols to extract soil macro, secondary, and micronutrients. Many of these nutrients are also indicators of soil quality reflecting soil condition/productivity, and measured with the same protocols and exhibiting the same spatio-temporal variability (Harisso et al. 2018). Recently, long-term data from different verified sources were aggregated into a national database, the Fertilizer Recommendation Support Tool (FRST). This tool was designed to improve recommendations related to soil fertility, but can also serve to aid in the interpretation of indicators of soil quality included in the Comprehensive Assessment of Soil Health (CASH) (Moebius-Clune et al. 2016).
Increases in K concentrations at the 0 to 5 cm depth of the WWCT watershed were likely due to the release of K from minerals. The textural designations within unit 8 (WWCT) are more variable than units 1 (STGP) and 5 (WWMT) at the 0 to 5 cm depth (tables S1, S2, and S3). In addition, erosion and disturbance within areas of the WRE have resulted in the removal of part of the organic-rich A horizon, the result of which has been the mixing of the upper portions of the B horizon into the 0 to 5 cm depth (Ahuja et al. 1988), likely causing an increase in mineral derived K. Previous research on the WRE has verified that tilled winter wheat management increased soil erosion while conservation tillage management coupled with winter wheat and a summer forage mitigated the negative effects of agronomic management on soil erosion (Garbrecht et al. 2015). Sources of K other than fertilizer are primarily clay minerals (~90%) and plant residues (Brady and Weil 2008). Greater K is removed in grain and straw than stored in residues particularly in a tilled monoculture that is grazed and/or hayed. The distribution of K across a 1,670 km2 area of Switzerland was due to parent material, soil type, and agricultural practices (Blanchet et al. 2017).
Similarly, differences in P and SO4 concentrations were more likely related to soil mineralogy than the turnover of organic material in SOM and residues. Mollisols in the Great Plains contain sulfur (S) bearing minerals associated with genetic horizons that supply SO4 (Ellert and Bettany 1995). Globally, mineral derived S is less common than organically bound forms of S and S associated from industrial deposition (Brady and Weil 2008).
Concentrations of secondary nutrients (Ca and Mg) were significantly higher in unit 1 (STGP) (table 4). Calcium concentrations were greater at all depths of the STGP relative to the WWMT and WWCT treatments (table 4). In contrast, Mg concentrations associated with Mg bearing minerals were significantly greater at the 15 to 20 cm depth of the STGP. A similar trend with depth was found in the WWMT and WWCT treatments. The depth and concentrations of Ca and Mg were previously shown to vary with mineralogy, SOM, and inputs of organic material (Takoutsing et al. 2017). As previously mentioned, other researchers have measured higher K concentrations in no-till systems across a range of soil textures (Franzluebbers and Hons 1996; Karathanasis and Wells 1990). The association of higher Ca and Mg concentrations in grasslands has also previously been documented (Frau et al. 2020).
Phosphorus had the greatest variation resulting from the statistically significant interaction between land management and landscape position (table 4). Specifically, unit 8 (WWCT) contained the highest P concentrations across landscape positions (42 to 30 mg kg–1), unit 5 (WWMT) contained around half as much P (18 to 11 mg kg–1), and unit 1 (STGP) around a third the concentration (9.19 to 7.75 mg kg–1) (table 4). Erosion of the Ap horizon is the likely cause of significant increases in P concentrations in the WWCT management at all landscape positions, Ca in the toe slope and SO4-S at the riser position. Miheretu and Yimer (2018) found a similar pattern in Ca and P concentrations between landscape position and land use. Clay minerals that contain and release SO4-S and P are associated with the B horizon where clay accumulates in the soil solum.
Micronutrients. Terrain attributes had a significant effect on the concentration and location of micronutrients. Boron concentrations were highest in the STGP treatment in the 0 to 5 cm depth. Boron and Zn varied significantly due to land management by depth (table 5). Other researchers have confirmed that Zn concentrations varied by landscape position and depth in four soils ranging from silt loam to silty clay textures (Brubaker et al. 1993). Zinc concentrations were highest in the STGP and WWMT treatments at the 0 to 5 cm depth. This finding is supported by other no-till research conducted for 8.5 years in a silty clay loam resulting in higher Zn concentrations than under tilled management (Franzluebbers and Hons 1996) and in an Orthic Humic Gleysol, Brookston clay loam (Van Eerd et al. 2014). Copper and B varied by watershed treatment and landscape position (table 5). Distribution of Cu with depth followed a similar pattern in a coarse-silty, mixed, superactive, mesic Typic Haploxeroll (Shiwakoti et al. 2019). Copper was highest in the WWCT treatment on the riser and toe positions. Iron varied with watershed treatment by depth and by landscape position (table 6). Iron was significantly higher or was of a higher concentration in the 0 to 5 cm depth. The STGP treatment had significantly lower concentrations of Fe relative to the WWCT treatment at all depths and landscape positions (table 6). All of these micronutrients are associated with the constituents’ mineralogical micas and feldspars (Singh and Schulze 2015).
The interactive effects of watershed treatment by depth averaged across landscape position on concentrations of micronutrients in soil.
The interactive effects of watershed treatment by landscape position and by depth on soil iron concentrations (mg kg–1).
Carbon and Nitrogen Indicators Associated with Conservation Management. The repeated measures in mixed models analyses revealed that TSOC, TSN, and POMC varied with management by depth (table 7). Tillage redistributes plant materials, no-till systems stratify nutrients, and prairies concentrate nutrients in the rooting zone (Bruulsema and Ketterings 2008). Specifically, C and N are redistributed within the tilled winter wheat system (WWCT) or stratified within the minimal till and grassland systems (WWMT and STGP) resulting in a management associated artifact. Therefore, soil C and N concentrations (TSOC, TSN, and POMC) across managements (STGP, WWMT, and WWCT) were measured and presented on a 0 to 30 cm depth after concentrations were adjusted for bulk density (table 7).
Integrators of land use and inherent soil properties at 0 to 30 cm depth: total soil organic carbon (TSOC), resistant fraction of soil organic C measured via the residual C of acid hydrolysis (RCAH), RCAH as a percentage of TSOC, particulate organic matter C (POMC), and percentage POMC as a percentage of TSOC.
The lack of a difference in TSOC, TSN, and POMC with landscape position is partly due to the creation of replicate sampling stations within the delineated landform complexes that successfully explained and limited the spatial variability that confounded the interpretation of the effects of soil conservation practices and land use on edaphic measurements. The uniform zones of management sufficiently reduced the variability associated with terrain attributes enabling us to assign the remaining variability to land management and specific conservation practices. Additionally, fertilizer management in the agronomic systems has reduced much of the spatial variability in plant biomass that can be associated with landscape position. Biomass measurements taken in the fall and spring between 2018 and 2020 on WWMT and WWCT varied among managements but did not vary statistically within a treatment (WWMT and WWCT) due to landscape position (data not shown). Other researchers have successfully coupled management zones, precision agriculture, and conservation practices (Kitchen et al. 2005; Lerch et al. 2005). Increases and decreases of soil C and N fractions are the result of the amount of C and N returned to the soil as organic material, in this instance plant residues and the rate at which these organic C and N fractions decompose (Larson et al. 1972). Disturbance redistributes residues and aids in the fragmentation and oxidation of organic materials, which increase their rate of decomposition (Krull et al. 2003). This is likely the cause of decreasing stocks of C and N in the WWCT management of the WRE.
Both the STGP and WWMT managements contained significantly higher concentrations of TSOC, TSN, and POMC than the WWCT management (table 7). Previous research has verified that TSOC, POMC, and RCAH are correlated with one another (Dou et al. 2008). Particulate organic matter C concentrations on the WRE contained between 17% and 18% of TSOC. Concentrations of POMC measured across a variety of crop and pasture land ranged from 14% to 31% of TSOC (Fortuna et al. 2003b; Dou et al. 2008; Li et al. 2018). Other researchers have measured greater TSOC and POMC after 17 years of no-till versus tilled management in a silt loam (Liu et al. 2014). Badgery et al. (2014) measured higher concentrations of TSOC and POMC in permanent pasture relative to cropland. Although the percentage of POMC in TSOC was not significantly different between managements, the concentrations of POMC varied significantly between managements. The STGP and WWMT managements contained ~15% more POMC and ~20% more TSOC and TSN then that of WWCT management. The RCAH fraction in the 0 to 30 cm depth of crop and grasslands has been previously reported to range from 44% to 71% of TSOC (Dou et al. 2008; Fortuna et al. 2003a; Bhowmik et al. 2017). There was a trend, although not statistically significant, indicating that the STGP management contained ~8% more RCAH than the agronomic managements (WWMT and WWCT).
Although the WWCT contained less POMC than the STGP and WWMT treatments, slightly more of the total C was contained in the biologically active fraction (POMC) relative to the resistant, nonbiologically active fraction (RCAH). Concentrations of C contained in the RCAH fraction in the two agronomic treatments (WWMT and WWCT) were equivalent, although TSOC values were significantly greater in the WWMT treatment. Differences in the quantity of biologically and nonbiologically active fractions among managements relates to how C and N are cycled and retained in each system. Tilled agronomic systems tend to retain more C in active C fractions with mean turnover rates of days and slow pool turnover rates of years to decades (Paul et al. 1999). In contrast, grasslands and minimally disturbed systems retain more C in resistant fractions (Paul 1999). Cultivation has been shown to increase N turnover but not C turnover in agronomic systems relative to grasslands (Schimel 1986).
Nonmetric Multidimensional Scaling Analyses. The NMDS analyses contain data from all edaphic indicators, not just C and N. Therefore, TSOC, TSN, POMC, and N are presented as observations taken from each of the separate soil depths (tables 8, 9, and 10). The NMDS analysis corroborates the significant differences in edaphic properties resulting from interactions among land management, landscape position, and depth represented by the univariate three-way ANOVA. The percentage variation explained by the two NMDS axes are represented by the proximity of the measured variables to the axes in each table. Tables represent the correlation coeficients (r) between the individual variables (edaphic indicators) for a given sampling depth and the two NMDS axes. Significant variables differ from random (p < 0.05), reduced stress, and improved the goodness of fit by ≥0.050, which is in the range of good to excellent (Šmilauer and Lepš 2014). Edaphic properties that correlated with one or both axes displayed in a table had absolute r values greater than 0.50 and are bolded.
Nonmetric multidimensional scaling (NMDS) analysis of edaphic indicators at the 0 to 5 cm depth measuring the effects of land management and landscape position. Correlation coefficients (r) resulting from correlations between measured variables and each NMDS axis.
Nonmetric multidimensional scaling (NMDS) analysis of edaphic indicators at the 5 to 15 cm depth measuring the effects of land management and landscape position. Correlation coefficients (r) resulting from correlations between measured variables and each NMDS axis.
Nonmetric multidimensional scaling (NMDS) analysis of edaphic indicators at the 15 to 30 cm depth measuring the effects of land management and landscape position. Correlation coefficients (r) resulting from correlations between edaphic properties and each NMDS axis.
The NMDS analysis verified that texture, chemical (P, K, Ca, SO4, Zn, and B), and biological (TSN and TSOC) edaphic indicators contributed to the variability of axis 1 within the 0 to 5 cm depth (table 8). In contrast, the variability of axis 2 related to Fe concentrations. Other researchers have observed an increase in total P in cropland relative to grassland (Gosling et al. 2017), and that land use is a primary driver of P variability at the field and watershed scales (Page et al. 2005). Cade-Menun et al. (2017) measured higher Zn concentrations in grassland relative to cropland that also varied with depth. Iron could be higher in agronomic managements due to crop residues, trace amounts of Fe in fertilizer applications, and erosion of the plow layer resulting in increased clay minerals that contain Fe. Total SOC and TSN explained a similar proportion of the variation in axis 1 at the 0 to 5 cm depth (table 8). Given the same climate and soil type, prairie and other uncultivated soils are typically higher in TSOC and TSN. Cultivation particularly when coupled with disturbance such as tillage reduces C and N fractions >20 μm, a fraction that includes POM (Hassink 1997). Therefore, POM fractions provide a more rapid (one to four year) assessment of a management practice on soil C and N that can be used as a positive, neutral, or negative indicator by land managers.
The majority of the variability in axis 1 of the 5 to 15 cm depth was correlated with edaphic parameters associated with inorganic soil properties, while axis 2 correlated strongly with C and N stocks (table 9). Phosphorus, SO4-S, Ca, and Mg partially explained the variation in axis 1, and K partially explained the variation in axis 2. Micronutrients explained a portion of the variation in axis 1 (B) and axis 2 (Zn). TSOC, TSN, POMC, and N explained the majority of the variability in axis 2. As previously noted, tillage disturbance, nutrient stratification in no-till, and the distribution of roots in the STGPs were the cause of the overlapping treatment observations. The 5 to 15 cm depth is subject to mixing and inconsistent redistribution of residues in the WWCT treatment and contains the majority of plant roots in all treatments.
Inorganic edaphic indicators were correlated with the variability of axis 1 and 2 at the 15 to 30 cm depth. Phosphorus, Mg, Fe, and Zn partially explained the variation in axis 1, and texture, K, SO4-S, Ca, and Mg partially explained the variation in axis 2 (table 10).
Summary and Conclusions
A number of soil fertility measurements analyzed for decades are also used as measures of soil quality, varying on the same spatial scale. These include macro, secondary, and micronutrients that differ due to taxonomic properties, terrain attributes, and conservation practices. The majority of variables associated with soil classification and terrain attributes on the WRE varied with depth and landscape position. Measurements included clay content, macro (P and K), secondary (cations included in base saturation [Ca, Mg, and K] and SO4-S), and micronutrients [Fe, Zn, B, and Cu]).
The uniform zones of management delineated using a class I soil survey, grid, and zone sampling sufficiently reduced the variability associated with terrain attributes, enabling us to assign the remaining variability in edaphic properties to land management and specific conservation practices. As a result, TSOC, TSN, and POMC (the sand-sized soil separates remaining on a 53 μm sieve) did not vary with landscape position. Cultivation particularly in the tilled treatment reduced soil C and N fractions. The sand associated C fraction that includes POMC provided a rapid (one to four year) assessment of conservation management practices effects on edaphic properties used to measure soil quality. Use of a suite of edaphic indicators that vary significantly and others that remain constant with soil classification and terrain attributes, land management, and specific conservation practices can aid land managers to quickly assess whether transitions in land use and management have a positive, neutral, or negative effect.
Our efforts to parameterize the inherent spatial variability of these watersheds in order to determine the effects of land management on soil quality were successful. We established replicate sampling stations within defined landform complexes to explain and limit spatial variability that confounds researchers’ interpretation of the effects of soil conservation practices and land use on soil quality measurements. Current measurements provide information with respect to soil properties that affect water (H2O) movement, soil fertility, and plant nutrient efficiency, which can directly influence H2O runoff and quality. Future work will include use of in-situ near infrared (NIRS) technology coupled with a set of soil quality measures calibrated against NIRS scans to verify that measures of soil quality can be predicted using NIRS and a minimal data set (Fortuna et al. 2019). Measurements collected from the 4.6 m by 3.8 m replicated sized spatial zones created in this study will be used to validate agricultural models like the Agricultural Policy/Environmental eXtender Model (APEX) to predict the effects of alternative land management on water quality.
Supplemental Material
The supplementary material for this article is available in the online journal at https://doi.org/10.2489/jswc.2022.00160.
Disclaimers
Mention of trade names or commercial products in this publication is solely for the purpose of providing specific information and does not imply recommendation or endorsement by the US Department of Agriculture (USDA). The USDA prohibits discrimination in all its programs and activities on the basis of race, color, national origin, age, disability, and where applicable, sex, marital status, familial status, parental status, religion, sexual orientation, genetic information, political beliefs, reprisal, or because all or part of an individual’s income is derived from any public assistance program. USDA is an equal opportunity provider and employer.
Acknowledgements
This research contributed to the Conservation Effects Assessment Project and is a contribution of the USDA Agricultural Research Service (ARS) Southern Plains Long-Term Agroecosystem Research Site. We also thank Jeremy Dennis and Jarred Goedeke (soil scientists, Stillwater Soil Survey Office, Stillwater, Oklahoma) for providing us the Class I Soil Survey. Research was funded by the United States Department of Agriculture.
- Received September 18, 2020.
- Revision received May 11, 2021.
- Accepted May 11, 2021.
- © 2022 by the Soil and Water Conservation Society