Abstract
Conventional cultivation of claypan soils leads to soil and water quality degradation because of high runoff and associated soil erosion. The Goodwater Creek Experimental Watershed, which is part of the USDA Agricultural Research Service Benchmark Conservation Effects Assessment Project, Watershed Assessment Studies, was established to address these issues. Plot studies have highlighted trade-offs between erosion control and herbicide or nutrient runoff. There is a need for long-term field-scale evaluation of combined practices that reduce sediment, nutrient, and herbicide losses by runoff. A 36 ha field located in Missouri was under a conventional corn (Zea mays L.)–soybean (Glycine max L.) system from 1993 to 2003 with fertilizer application and tillage prior to planting in the spring. A precision agriculture system defined by two main management zones was implemented from 2004 to 2014: wheat (Triticum aestivum L.) and soybean in 60% of the field, and corn and soybean in the remaining 40%. The system included no-till, cover crops, atrazine split-applications based on weed pressure, variable rates of nitrogen (N), and variable rates of fall-applied phosphorus (P). The objective of this study was to compare runoff water quality from the two management systems, based on flow and load duration curves, cumulative distribution functions, and conclusions from replicated plot studies. The precision agriculture system did not affect annual runoff, but it did increase the frequency of low flows. Sediment losses were reduced by 87% as a result of no-till and cover crops. Atrazine and P losses were lower than expected, despite the lack of incorporation into the soil. Atrazine losses were possibly lower because of the wheat area acting as a buffer, greater atrazine adsorption and retention in the field, and faster decay. Dissolved P losses as a fraction of applied remained the same, likely because of greater adsorption and lower runoff risk when applying P. Finally, nitrate-N (NO3-N) losses decreased and resulted in an overall decrease of N losses despite a slight increase of ammonium-N (NH4-N) losses. Explanations included a greater soil water content, a different timing of N applications, and N uptake by cover crops. Building on these successes, an aspirational management system is proposed to further improve on the performance and practicality of the precision agriculture system.
- Key words
- conservation practices
- cover crops
- no-till
- nutrient management
- precision agriculture
- water quality
Introduction
The Central Claypan Area (Major Land Resource Area 113) encompasses 3 million ha in northeastern Missouri and south-central Illinois (USDA NRCS 2006) (figure 1). The area is dominated by high runoff potential claypan soils, and soil and water quality degradation within the northeastern Missouri portion of the region has been well documented over the last 20 years (Sadler et al. 2015; Veum et al. 2015). The Salt River Basin encompasses the majority of the Missouri portion of the Central Claypan Area, and numerous studies have documented significant stream and groundwater contamination by sediment (Baffaut et al. 2013), nutrients (Kitchen et al. 2015; Lerch et al. 2015c), and herbicides (Blanchard and Lerch 2000; Lerch et al. 2011b) within the basin. Long-term studies within the region have documented soil degradation in systems with intensive tillage, two-year crop rotations, and high chemical inputs (Veum et al. 2015). The Central Claypan Area has the highest estimated soil erosion rates for cropland within the state (Willett et al. 2012). Because of the soil and water quality degradation, this region was one of the sites chosen for both the Conservation Effects Assessment Project (CEAP) (Lerch et al. 2008) and the Long-Term Agroecosystem Research (LTAR) project (Sadler et al. 2015). One focus of CEAP was to conduct watershed-scale assessments to determine the collective effects of existing conservation practices on trends in water quality (Lerch et al. 2011a, 2011b, 2015a) as well as to identify sources of contaminants (Willett et al. 2012; Peacher et al. 2018) to support future targeting efforts.
The 72 km2 Goodwater Creek Experimental Watershed is located in the headwaters of the Salt River Basin. Overall, long-term changes in water quality at the watershed scale were modest and attributed to changes in herbicide usage, fertilizer management, and land use. For example, trends in herbicide concentrations were only observed when major changes in usage occurred in the watershed (Lerch et al. 2011a). Similarly, decrease in nutrient transport were attributed to land use changes and improved nitrogen (N) fertilizer management in the watershed, which shifted the amount, timing, and placement of fertilizer applications (Lerch et al. 2015a).
Two factors help explain the limited improvements in water quality observed in Goodwater Creek. One is the importance of streambank erosion as sources of sediment and nutrient contamination. A study in two watersheds within the Salt River Basin showed that streambanks contributed an average of 83% of sediment, 67% of total phosphorus (P), and 23% of the total N transported from these watersheds underlain by claypan soils (Willett et al. 2012; Peacher et al. 2018). Thus, cropped fields may not currently be the major source of sediment and P within these watersheds, and practices that stabilize streambanks will be needed to see improvements in stream water quality.
Another factor is the extent of conservation practice implementation. In 2006, 14.7% of the Goodwater Creek Experimental Watershed was protected by at least one conservation practice (Lerch et al. 2011a). Assuming the same implementation rate as occurred between 1990 and 2006, there may currently be 23% of the watershed protected by the time of this writing. Moreover, practices heavily focused on grass waterways and terraces, which are very effective for field erosion control but have limited effect on dissolved-phase herbicide or nutrient transport. They are also among the most costly conservation practices to implement. There is a need for practices that producers are willing to implement and that control erosion as well as reduce nutrient and herbicide transport. Plot studies have highlighted trade-offs between practices designed to reduce edge-of-field transport of specific constituents. One major trade-off concerns no-till, which is effective at reducing soil erosion, but increases the transport of surface applied herbicides (Ghidey et al. 2005; Lerch et al. 2013) and dissolved P (Smith et al. 2015; Verbree et al. 2010). Thus, there is a need to evaluate the combined effect of multiple practices on water quality with the goal of addressing multiple contaminants and transport modes (i.e., dissolved-phase versus sediment-bound).
One approach to an integrated system of practices that could sustain production yet benefit water quality is a precision agriculture system (PAS). Using profitability data analysis, Kitchen et al. (2005) developed a comprehensive PAS for a 36 ha field located in the Goodwater Creek Experimental Watershed. From 1993 to 2003, the field was managed under a conventional system that represented the dominant cropping practices in the region: a corn (Zea mays L.)–soybean (Glycine max L.) rotation with high chemical inputs and intensive tillage. The depth of topsoil above the claypan, a dense and high clay content layer with very low permeability (Kitchen et al. 1999), has been linked to runoff generation and associated transport of pollutants (Lerch et al. 2005) as well as to productivity (Kitchen et al. 2005) on this field. The spatial variability in topsoil depth is related to soil loss over the last 100 to 150 years within the field, and controls the soil quality and crop productivity patterns currently observed. It may also control water quality, as Mudgal et al. (2012) demonstrated using a simulation model.
The PAS design included three zones with one management system for zones A and B (figure 2), and one for zone C, based on the expectation that these systems would improve profitability and conserve soil and water resources. The PAS was targeted to these areas with the following objectives: (1) reduce production costs while stabilizing crop yields and improving crop water use efficiency; (2) improve surface runoff quality by reducing sediment, herbicide, and nutrient loss; (3) reduce the leaching of nitrates (NO3) to groundwater; and (4) improve soil quality by reducing topsoil loss, improving soil structure and infiltration, and building organic matter. After 11 years of PAS on this field, Yost et al. (2017) concluded that the greatest production benefit of PAS was reduced temporal variability of yields. They did not observe improvement in yields. Yet, profitability was sustained (Yost et al. 2019) in spite of the additional costs incurred under the PAS. Producers and farm land owners would most likely be motivated to implement these practices because of reduced variability (and therefore uncertainty) in profitability and income. Enhancements of soil and water quality by the PAS are desirable environmental outcomes that need to be evaluated.
One challenge of field experiments is to have access to replicated fields for statistical comparison of the treatments. Helmers et al. (2012) used a balanced incomplete block design conducted at three locations across 12 field-size row crop watersheds in Iowa to study the effects of varying filter strip size and placement on contaminant transport in runoff. This highlights the difficulty in designing a replicated catchment experimental study with a randomized complete block design. An alternative is to define blocks and drainage areas that are sufficiently small to be replicated within a larger area and randomly assigned a treatment (Bosch et al. 2012). Such designs can result in drainage areas too small to address the problem of spatial variability within a drainage area. Another is to have paired watersheds (Udawatta et al. 2011), in which a baseline period is used to derive the relationship between the outcome in two or more watersheds, and a treatment period is used to investigate the effect of a specific treatment in one of the watersheds.
Implementing treatments successively in the same field or watershed is another option: these designs are sometimes called pre-post, before-after, or paired in-time. Analysis of variance techniques cannot be applied because the data points for each treatment are not independent of each other nor are treatments randomly assigned to experimental units. However, the data can be analyzed on the basis of cumulative distribution functions (Spooner et al. 1985).
This is the approach used in this study, in which we use statistical methods based on daily runoff and export of sediment and nutrients to compare the impacts of management over the two periods. The goal of this study was to compare the environmental impacts of the conventional and PAS management systems on the same field over separate 11-year periods (1993 to 2003 for conventional and 2004 to 2014 for PAS). Specific objectives of this study were to (1) compare runoff discharge and water quality over the two periods, and (2) use adjacent replicated plot studies on soil and water quality to inform the results.
Materials and Methods
Study Field. The area under study is a 36 ha field located in Boone County in Missouri (39°13’48” N, 92°07’12” W). The field soils and topography are characteristic of those found on a dominant fraction of the cropland in the Goodwater Creek Experimental Watershed: low slopes (<3%) and claypan soils that belong to the Adco-Mexico-Putnam association. In this field, the top of the claypan was found at depths ranging from the soil surface on the steepest parts of the field to as much as 100 cm around the channel at the bottom of the field. Average annual precipitation over the study period (1993 to 2014) was 1,022 mm, and average annual minimum and maximum daily temperatures were 6.8°C and 17.6°C, respectively.
During the conventional period, the field was managed under a corn–soybean rotation with fertilizer application and tillage prior to planting in the spring. In 1995, heavy spring rains prevented timely planting of corn, and a sorghum (Sorghum bicolor L.) crop was planted in late June. In addition, a tillage operation took place in December of 1995 to address the soil compaction caused by the heavy spring rains. Otherwise, tillage occurred in the spring: one or two passes (disk or cultivator) prior to planting to loosen the ground, kill weeds, and prepare the seedbed, and one pass just before planting to incorporate fertilizers and herbicides. Fertilizer applications were uniform throughout the field and made during the corn year of the rotation (table 1). Fertilization usually consisted of preplant broadcast liquid urea-ammonium nitrate (UAN) solution, inorganic P (either diammonium phosphate [DAP] or monoammonium phosphate [MAP]), potassium (K), and lime, which were then incorporated by tillage. Herbicides included atrazine and alachlor or metolachlor on corn, and alachlor or metolachlor and imazaquin on soybeans. Table 1 shows the application rates of agricultural inputs considered in this study: N, P, and atrazine.
During the PAS period, the crop area was divided into three zones (figure 2): a wheat (Triticum aestivum L.)–soybean rotation in zones A and B (21 ha or 60% of the field) where corn had been found unprofitable for the most part, and a corn–soybean rotation in zone C (15 ha or 40% of the field). Both cropping systems included winter or summer cover crops as needed, and there was no use of tillage for any of the crops. One exception to that protocol occurred in August of 2007 when sediment was leveled to mitigate pooling of water in the area immediately above the weir, resulting in soil disturbance. Zone B included the land around the drainage channel of the field. For practicality reasons, zone B was initially managed with the same system as in zone A (soybean–wheat) but included permanent stiff-stem grasses planted at critical points along the drainage channel to control erosion. Herbicide applications for the grain crops caused the stiff-stem grasses to practically disappear by 2008, a time by which benefits from no-till and cover crops had become sufficient to not require specific management in zone B.
Both N and herbicide applications were split between planting and in-season (table 2). Nitrogen fertilization for corn and wheat included surface-applied UAN at planting and in-season top-dress applications with variable rates as a function of need determined using commercial ground-based crop canopy light reflectance sensors (USDA NRCS 2009; Kitchen et al. 2010). A small amount of starter fertilizer was placed with the seed at planting. In 2009, subsurface application of N fertilizer at planting was implemented to decrease volatilization losses. This was not repeated in 2011 and 2013 because of practical issues with the subsurface application, which resulted in uneven application rate. Instead, a volatilization inhibitor was added each time that N fertilizer was surface applied. In-season side-dressed N was dribbled on the ground surface near the plant instead of being knifed-in because the size of the plants at the time of the operation required a high-clearance applicator, and the one available could not support soil-engaging tools. Although sensor-controlled N application is possible earlier in the season when a tractor system can be used for subsurface fertilizer placement, later application is advantageous both to better synchronize N application with corn N uptake and to improve the sensitivity of sensor measurements to crop N differences (Kitchen et al. 2010). Table 2 includes rates for the compounds considered in this study: N, P, and atrazine. A full description with information for K and lime applications is given by Yost et al. (2017).
Runoff Quantity and Quality Data. Measurement of runoff stage and discharge over a 3:1 broad-crested v-notch weir constructed at the outlet of the field has been described in detail by Baffaut et al. (2015). The vertical wall of the weir extends to 1 m belowground such that the structure captured all surface runoff and the shallow subsurface lateral flow from the field. While the field total area was 36 ha, the drainage area at the weir was limited to 31.6 ha (Mudgal et al. 2012) and included 12.5 ha of zone C (41% of drainage area) and 19.1 ha of zones A and B (59% of drainage area). Baffaut et al. (2013), Lerch et al. (2015b, 2015c), and Mudgal et al. (2012) have described the equipment used to collect automated flow proportioned runoff samples; the laboratory analysis methods used to determine sediment, nutrient, and herbicide concentrations; and the data processing methods to calculate daily loads. While there were some changes in equipment over the 22 years, the sampling regime (discharge to trigger sampling and sampling flow interval) remained the same over the whole period.
Data Analysis. The comparison of environmental outcomes during the conventional and PAS periods relied on statistical methods based on daily edge-of-field runoff and export of sediment, atrazine, and nutrients (NO3-N, ammonium-N [NH4-N], and dissolved P). Daily discharge (expressed as a depth over the drainage area) and constituent loads obtained during the conventional and PAS periods were compared based on (1) average annual totals expressed as a fraction of applied amounts, (2) flow and load duration curves, and (3) cumulative distribution functions. Annual mean, median, and standard deviation of daily precipitation, flow, and loads were calculated for each management system. For atrazine, the analysis used years when atrazine was applied and measured: odd years from 1993 to 2001 and 2005 to 2013 for the conventional and PAS systems, respectively. The year 2003 was excluded because atrazine was applied but inadvertently not analyzed in the water samples during that transition year. For other variables, all years were used. Losses were also calculated as a fraction of applied amounts. If the application was not annual, losses and applied amounts were aggregated over the relevant number of years. For N, we combined the NO3-N and NH4-N losses in order to make the relative comparison with total N inputs. No statistical test was applied since annual values do not represent a true replication of the system. Average annual values determined the direction and magnitude of the change.
Flow duration curves relate daily flow value to their frequency of exceedance, also known as the flow duration interval. After ranking daily flow values in decreasing order, each value was assigned a frequency calculated with the Cunnane (1978) formula (equation 1): 1 where pk is the frequency of exceedance (i.e., the frequency at which a daily value exceeds this value), k is the rank of the daily value, and N is the total number of days in the data set. With k varying from 1 to N, pk varies between 0 and 1. The load duration curve is the plot of daily constituent loads as a function of the flow duration interval of the corresponding daily flow value. To determine if a water quality change had occurred as a result of the management, we applied a variant of the quantile-quantile plot method proposed by Spooner et al. (1985) where concentration data are ranked, and mean concentrations are calculated for each percentile for both periods and plotted against each other. The regression slope is then compared to one. In this study, availability of continuous discharge and flow weighted concentration data allowed a similar analysis based on daily loads. Median loads calculated for each percentile of flow duration interval for the two periods were plotted against each other, and the regression slope statistically compared relative to one. The method requires a similar range of flows during the before and after treatment periods, and this condition was met as well (see below).
Cumulative distribution functions (CDF) relate daily values to the frequency at which the variable (flow or load) was less or equal to that value. Although similar to the flow or load duration curve, the CDF considers the frequency of the daily load value directly, not the frequency of the corresponding flow value. Cumulative distribution functions were statistically compared using the non-parametric Kolmogorov-Smirnov test using the NPAR1WAY Procedure of SAS (SAS Institute 2011) at p ≤ 0.05 for daily precipitation, daily flow, and daily NO3-N, NH4-N, dissolved P, and sediment loads using 1993 to 2003 data for the conventional period, and 2004 to 2014 for the PAS period. For atrazine, the test was based on only the years where atrazine was applied and analyzed in the runoff samples: odd years from 1993 to 2001 for the conventional period, and odd years from 2005 to 2013 for the PAS period.
Adjacent Replicated Plots. Thirty large plots (20 × 200 m) adjacent to the field under study were established in 1991. Replicated management systems including the conventional system in this study, a no-till corn–soybean system, and a no-till corn–soybean–wheat with cover crops were consistently implemented from 1996 to 2010, with each phase of the rotation replicated three times each year. Runoff and water quality were monitored in the corn phase of the rotation from 1997 to 2002 on two of the three replicates. No management replicated the PAS because it was specifically designed for the field and included multiple zones, which was not feasible on plots. Data from these plots were analyzed for runoff (Baffaut et al. 2015) and herbicide losses (Ghidey et al. 2005). In 2008, a soil quality assessment was conducted (Veum et al. 2015) using soil samples collected from 0 to 5 and 5 to 15 cm depth layers across three landscape positions (summit, backslope, and toeslope) to contrast the conventional system and the no-till corn–soybean–wheat rotation with cover crops. In addition, long-term effects of these cropping systems were estimated using model simulation (Baffaut et al. 2019). The conclusions of these analyses were used as needed to explain the differences in water quality from the field under study during the conventional and PAS periods.
Results and Discussion
Average Annual Values. Precipitation did not change appreciably over the two periods (table 3). The number of years under each management (11) was sufficiently long that each period included normal, wet, and dry years. Average annual precipitation varied by only 4% between the two periods. The driest summer during the conventional period was 1999 (74 mm) and the wettest was 1993 (654 mm). During PAS, the driest summer was 2001 (113 mm); the wettest was 2008 (785 mm). The wettest springs occurred in 1995 (568 mm) and 2010 (469 mm). The most important differences in seasonal precipitation between the two periods were obtained for summer (15% difference) and fall (14%) with the PAS period experiencing more precipitation than the conventional period during these two seasons. Similar to precipitation, annual runoff volumes were similar over the two periods.
Sediment losses under the PAS were reduced by 87%, despite including the high sediment losses caused by waterway grading in August of 2007. In that year, losses (3,830 kg ha−1) were seven times the next highest annual sediment loss (546 kg ha−1) during the PAS period. No-till combined with cover crops at all times reduced soil movement to very small amounts. If we exclude 2007, the average annual soil loss during the PAS period was 285 kg ha−1, or only 6% of the average annual loss during the conventional period.
There was ~60% less atrazine applied to the field during the PAS because corn was only grown in 40% of the field. These lower quantities translated to 40% less atrazine transported with runoff. As a result, loss as a function of applied was greater during the PAS period (2.7% compared to 1.7%; table 3).
Applied P increased under the PAS. Fertilization was informed by soil P tests, which showed excess at the boundary between zones C and A, and deficiencies in the area just upstream of the weir. During most years of the PAS, fertilizer P applications were designed to increase soil P to support the growth of two crops and a cover crop annually. As a result, P applications almost tripled (2.8 ratio) during the PAS period and were especially large in the area closest to the weir. The P fertilization timing was quite different with some P applied each year under the PAS period compared to every other year under the conventional system. The increase of applied P led to a tripling of P (3.6 ratio) transported by runoff. However, the ratios of lost P to applied P remained similar between the two systems, at 5% and 6% of applied during the conventional and PAS periods, respectively.
Under PAS, N was applied on corn and wheat at planting and in-season with variable rates. In-season variable N applications resulted in 50 to 110 kg N ha−1 applied on wheat and 0 to 235 kg N ha−1 on corn. Average total N applied over the whole areas ranged from 160 to 200 kg N ha−1 for corn and 107 to 129 kg N ha−1 for wheat. In addition, uniform rates of 18 to 70 kg N ha−1 were applied to cover crops because it was included in the DAP or MAP used for P fertilization. Overall, total N applied was almost the same (1% difference) during the two periods. However, the application timing was more spread out during the growing seasons of the various crops grown under PAS. In spite of equal amounts of applied N, average annual N losses under PAS were 39% lower than those under the conventional system. The N lost in runoff as a fraction of applied N was reduced by close to 44% (table 3).
Flow Duration Curves and Cumulative Distribution Functions. Figure 3 shows the flow duration curves and load duration curves for sediment, atrazine, dissolved P, and N. These curves and the regression slopes between median loads for each flow duration interval (figure 4) highlight the impact of PAS on a daily basis. Precipitation duration curves are not shown but were identical.
The range of daily flow values can be divided into two parts. For flow volumes >1 mm d−1 over the drainage area, the two flow duration curves appear identical. However, the quantile-quantile plot showed greater flow during PAS than during the conventional period for high runoff events, and the regression slope of the quantile-quantile plot was significantly greater than 1 (table 4). This corresponded to when surface runoff on nearly saturated soil was the dominant process (high runoff events). Low daily runoff volumes (<1 mm d−1) were more frequent under PAS than during the conventional period. The accumulated flow volume when daily flow <1 mm d−1 represented only 3% of the annual flow volume during either period. However, this happened during 78% of the non-zero flow days, i.e., very frequently. In spite of this, the Kolmogorov-Smirnov test showed no significant difference between the two flow CDF (table 4).
Comparing loads for the same flow interval between the two systems clearly shows that daily NO3 and sediment loads decreased (slope significantly less than 1) as a result of PAS, which was confirmed by the Kolmogorov-Smirnov tests. Given the similarity of flows during the two periods, decreasing sediment and NO3-N loads were clear benefits of the PAS for these two constituents, decreasing concentrations and loads. On the other hand, results for NH4-N were inconclusive. The regression slope of the quantile-quantile plot indicates an increase of NH4-N transport, but the Kolmogorov-Smirnov test did not indicate any difference. Results were similarly inconclusive for atrazine, with the slope indicating a decrease of daily loads while no significant difference was detected between the two CDF. As we have seen previously, dissolved P loads increased but the increases were explained by increased inputs.
Effect of Management on Discharge. This analysis confirmed previous analyses of runoff volumes under tilled and no-till systems. Using data from the adjacent plots, Ghidey et al. (2005) showed no difference in growing season runoff volume from the conventional system and a no-till corn–soybean system. Similarly, analysis of event data from these large plots confirmed no difference between tilled and no-till cropping systems, with or without cover crops (Baffaut et al. 2015). Field-scale data did show some differences in discharge between this field during the conventional period and other fields with no-till management systems. However, different bio-physical characteristics of the fields (i.e., slope, slope length, and depth to clay) and differences in management prior to the experiment may have affected the results (Baffaut et al. 2015). In this study, the field was the same over the two periods, and the management prior to the conventional management period, which had been in place since 1981, was similar to the conventional management: soybean, grain sorghum, and wheat production under disc and cultivator tillage with inorganic fertilizers (Lerch et al. 2005). Thus, initial conditions at the beginning of the conventional and PAS periods were likely to be similar. While full soil characterization was not repeated in 2003, 1993 and 2003 field mean levels of P (29 and 31 kg P pentoxide [P2O5] ha−1, respectively), K (189 and 158 kg K2O ha−1, respectively), and pH (6.1 and 6.3, respectively) confirm this assumption.
The increased frequency of daily runoff volumes <1 mm d−1 suggests that soil water content was greater during that period and led to excess water moving as subsurface flow resurfacing before getting to the field outlet, or captured by the foundation of the weir. This result could indicate greater water holding capacity in the topsoil layer under PAS. However, studies on the plots adjacent to this field did not detect any significant differences in hydraulic conductivity in tilled and no-till systems (Jiang et al. 2007; Jung et al. 2007), or macroaggregate stability and bulk density (Veum et al. 2015). Similarly, Haruna et al. (2018) have detected no significant differences in bulk density or water holding capacity of the topsoil layer (20 cm) between tilled and no-till systems with and without cover crops on a somewhat poorly drained soil. Thus, the increase of interflow remains unexplained by supporting experiments and analyses. One possibility is that PAS did increase the duration of runoff events, suggesting that runoff velocity was reduced. Analysis of hydrograph dynamics has not been completed at this time.
Effect of Management on Sediment Loss. The quantitatively largest impact of PAS on surface runoff water quality was the sharp reduction in sediment transport. Average annual losses were reduced by 87% (95% if we exclude 2007; table 3), and daily loads for similar discharge were reduced by 50% to 90%. These data demonstrated that no-till in combination with cover crops was very effective at reducing soil erosion on these soils. While no-till alone is very effective at reducing sediment losses (Verbree et al. 2010; Zeimen et al. 2006), our own observations on adjacent no-till corn–soybean plots over the 1991 to 2009 period showed that when no-till was implemented without cover crops, rills formed in the area of concentrated flow. These rills then increased in size with subsequent years unless they were flattened by localized tillage. Cover crops provided greater stability to the soil and prevented these rills from occurring. This is consistent with results from replicated plots on claypan soils by Zhu et al. (1989), who showed that cover crops were a critical aspect for a successful no-till soybean crop when implemented on claypan soils: cover crops reduced soil loss from plots by 85% to 90% compared to the same no-till system without cover crops.
Effect of Management on Herbicide Loss. Water quality data from the adjacent replicated plots showed that atrazine losses were two to three times greater under the no-till corn–soybean system than when the chemical was incorporated into the soil by tillage (Ghidey et al. 2005). A modeling study based on the same plots corroborated the finding over a 14-year long simulation period (Baffaut et al. 2019). The result was confirmed at field scale by contrasting losses from this field under conventional management to a field under no-till corn–soybean during a subset of the conventional period (1997 to 2001). Lerch et al. (2013) also reported significantly greater atrazine losses in runoff from a claypan soil when comparing a no-till system to two tilled treatments in which atrazine was incorporated. Thus we expected treated area atrazine losses to be two to three times as high as during the conventional period. During PAS, the treated area was 40% of the field. The overall total applied amount was 37% of what had been applied during the conventional period because of varying application rates from year to year as a result of the split application. Based on these inputs and the conclusions derived from these previous studies, we expected losses to be 74% to 111% of what was lost during the conventional period. Yet, average annual losses were only ~60% and daily losses were 56% to 79% of those under the conventional system (table 4). Expressed as a fraction of applied, losses were only 60% greater (2.7% versus 1.7%; table 3 and figures 3 and 4) as opposed to 100% to 200% as expected from previous studies. Overall, the no-till PAS did not increase atrazine losses in runoff as much as expected.
One reason was the placement of atrazine on the far end of the field (zone C), relative to the monitoring weir. Under the zone management approach used, the remainder of the field was in wheat (zones A and B) during the vulnerable atrazine transport period from April through June (Lerch et al. 2011b). Thus, atrazine had to move from the south end of the field, by either surface runoff or interflow, through the growing wheat, which functioned as an effective vegetative buffer strip that retained atrazine in the field. Further, the wheat area was greater than that of the corn, creating a source-to-buffer-area ratio of <1:1. On claypan soils, grass buffers are estimated to reduce atrazine transport by 15% to 21% at source-to-buffer-area ratios up to 20:1 (Lerch et al. 2017) through a combination of runoff infiltration and herbicide adsorption (to plant residues and surface soil). In that study, the role of adsorption was greater than infiltration, demonstrating that vegetative buffers are effective for poorly drained soils with limited infiltration or when saturated. Since PAS did not decrease runoff on an annual basis compared to the conventional system, increased infiltration may have been limited as a retention mechanism. However, if the PAS did indeed reduce runoff velocity (figure 3; runoff flow duration curves), this allowed more time for atrazine to interact with crop residues and surface soil, presumably leading to greater atrazine adsorption and retention in the field compared to that previously observed for conventional no-till systems (i.e., corn grown on whole field or plot).
Once retained in the wheat, the atrazine may have also degraded faster in the soil. Grasses have been shown to increase atrazine degradation in the rhizosphere by a combination of stimulated microbial activity and production of atrazine degrading phytochemicals (Willett et al. 2016). Members of the Poacea family, which includes wheat, have been documented to produce benzoxazinone compounds that can nonenzymatically degrade atrazine (Niemeyer 2009; Willett et al. 2016). Long-term atrazine usage on this field may have also resulted in microbial adaptation and accelerated atrazine degradation in soil over time as has been widely documented (Krutz et al. 2010; Lerch et al. 2017). The 2008 soil assessment (Veum et al. 2015) showed that no-till and cover crops increased microbial biomass carbon (C) by 64%. Microbial biomass C represents an estimate of the microbial population that is available to decompose residues and organic compounds (Vance et al. 1987). The process of adaptation occurs quickly (Zablotowicz et al. 2007), and atrazine half-life in soil can decrease from 60 to 120 days in nonadapted soil to 1 to 12 days in adapted soil (Krutz et al. 2010; Lerch et al. 2017). Atrazine has been applied to this field 11 times since 1993, and microbial adaptation and accelerated degradation in this field are probable.
Finally, there also might have been greater adsorption of atrazine to soil caused by the increased soil organic C content of the top-soil layer under PAS. The 2008 soil assessment (Veum et al. 2015) showed that no-till and cover crops increased organic C in the top 5 cm of soil by 60%. Soil organic matter at the soil surface does increase adsorption of atrazine to soil particles (Ben-Hur et al. 2003), and decrease its ability to move with surface runoff.
Effect of Management on Phosphorus Loss. Many studies report greater dissolved P losses from no-till grain systems compared to tilled ones because the fertilizer remains available on the soil surface instead of being incorporated by tillage: Smith et al. (2015) reported growing season losses 2 to 4.5 times those from tilled systems, and Verbree et al. (2010) reported losses 2 to 9 times greater for no-till. With little reduction of runoff on claypan soils, not incorporating P only increases the amount picked up by runoff (Zeimen et al. 2006). Modeling results on the plots adjacent to this field showed that no-till would double dissolved P losses (Baffaut et al. 2019). As we have seen before, PAS resulted in 2.8 times as much P applied than during the conventional period. Based on published studies and with almost 3 times as much P applied, we expected the losses to increase by as much as a factor of 6. Yet, annual dissolved P losses increased by a ratio of only 3.6 (an increase comparable to the increase in inputs), and daily losses for similar discharge increased by 38% to 131%.
Phosphorus applications during the PAS period were spread out during the year instead of being coordinated with planting operations in April or May, just ahead of the period with the highest risk for runoff. While some P was still applied in April or May from 2004 to 2014 with the starter fertilizer, buildup applications were scheduled in October after corn or soybean harvest, when the ground was clear and runoff less likely. In 2005, an October application was not feasible and it occurred in February of 2006. Risk of P runoff losses is greatest immediately after application (Smith et al. 2018). For three of the four conventional period P applications, runoff occurred within one week of P application. The fourth one was followed by a 90 mm event nine days after application. In contrast, no runoff event occurred within one week of any of the fall or winter P applications during the PAS period. Some did follow the spring applications, but the applied amounts were much smaller. Fall and winter P applications are problematic in regions where soils may freeze, which includes Missouri. Frozen conditions were considered in scheduling the February application, and soil temperatures were always >0°C. In Missouri, the long, warm, and dry autumns provide opportunities for P buildup at a time of minimal risk of runoff. In addition, the presence of growing wheat and cover crops increased evapotranspiration and P uptake and further decreased the likelihood of P loss with runoff. This helped explain identical P loss as a fraction of applied in spite of not incorporating the fertilizer.
The increased P sorption to soil may also explain the lower than expected P losses. In the 2008 soil quality assessment (Veum et al. 2015), the no-till system with cover crops resulted in increased soil organic C, which increases the affinity of P to soil particles similarly to herbicides. Indeed, the same soil quality assessment found 33% less extractable P under the no-till system with cover crops than under the conventional system.
Effect of Management on Nitrogen Loss. Applied UAN consists of 50% urea, 25% NH4-N, and 25% NO3-N. Ammonium-N typically does not remain in the soil very long because of the nitrification and volatilization processes. Additionally, when soil water content is near saturation, denitrification can occur. During the PAS period, NO3-N daily losses in runoff were significantly smaller (about half) than the losses observed during the conventional period (table 4). Ammonium-N losses in runoff were greater during the PAS period than during the conventional period according to the quantile-quantile plots. The Kolmogorov-Smirnov test showed no significant change. Zhu et al. (1989) also found less NO3-N and more NH4-N in surface runoff when comparing losses from a soybean field with and without cover crops. As in this study, the field was underlain by a soil with a claypan within its profile.
All N transformation processes (nitrification, volatilization, denitrification, and immobilization) were likely taking place, but in varying degrees during both management periods. These processes largely are governed by timing, amount, placement, and source of N fertilization; C:N ratio of plant residues and the degree of incorporation of these residues into the soil; and soil water and its impact on soil oxygen (O). Aspects of all these processes were different between the two management systems contrasted in this investigation. As seen previously, the higher frequency of low flow days (<1 mm d−1) under PAS suggests that soil water content was greater during that period. Anaerobic conditions may have been more frequent, which, among other things, prevent nitrification and promote denitrification. With less nitrification, more NH4-N and less NO3-N was available for transport by surface runoff and subsurface flow.
In addition, during the conventional period, N fertilizer was applied in April at planting, a period that is near the beginning of the high runoff period (April to June). Thus, there is high potential for loss of N with runoff and subsurface flow. During PAS, some N was applied in April, but less. The side-dress application occurred in June, when plants are waist high and dewater the soil profile more effectively, and precipitation events start to decrease in magnitude and frequency. The timing of the N application may have contributed to overall lower losses of NO3-N during the PAS period.
Cover crops may also have had an effect on the availability of N in the soil profile. Farm managers use cover crops in part because these crops can capture excess N remaining in the soil when grown following a harvested cash crop, and immobilize it in their biomass when it would otherwise leach through the soil. That N is then released when the cover crops are killed. Thus cover crops reduce the availability of NO3 in the soil profile and the likelihood of transport by surface and subsurface flow.
Toward an Improved Management System. The PAS showed significant benefits in terms of reduced sediment and NO3 losses in runoff. Anticipated trade-offs between sediment reductions and increases in atrazine and dissolved P were not observed. One challenge with PAS was the increasingly difficult management of weeds where soybean and wheat were grown, because of herbicide resistance of weeds like amaranth (Amaranthus palmeri). The problem did not occur where corn was grown because atrazine was applied every other year. The zonal management also complicated the organization and timing of field operations. A solution may be to switch to a three-year no-till corn–soybean–wheat with cover crops, variable rates of fertilizers, and split applications of atrazine. This aspirational system would combine the advantages of single crops across the whole field with the ability to apply atrazine once every three years to manage weed resistance.
Without the presence of a large area in wheat to trap and decay atrazine and its derivatives, atrazine losses may increase. To avoid this possibility, a rotary harrow could be used to incorporate the chemical into the first 5 cm of soil without significantly disturbing the residue cover. Lerch et al. (2013) showed that such incorporation can significantly reduce herbicide losses compared to a no-till system, without significantly increasing soil losses. Such a compromise could be beneficial for other chemicals as well: other herbicides and P fertilizer in particular. In addition, a vegetated buffer placed at the outlet of the field could further reduce sediment-bound and dissolved herbicides by trapping sediment, increasing infiltration, and increasing adsorption of the chemicals. A source-to-buffer-area ratio of 20:1 would amount to 1.6 ha of perennial grasses at the bottom of the field.
In this study, all but the starter N fertilizer was surface applied, and a stabilizer was added to control volatilization of urea. Plant uptake could only occur after precipitation and infiltration moved N into the root zone. However, precipitation also increased the risk of loss by runoff. Subsurface placement of N at planting and in season would decrease the potential for volatilization and transport by runoff, and make the fertilizer readily available to the plants. This would require appropriate equipment for doing so, especially in-season when plants are knee high or taller. Costs will certainly affect this decision.
Summary and Conclusions
After 11 years of conventional management in a 36 ha field, productivity data led to the definition of two main management zones: wheat and soybean in 60% of the field, and corn and soybean in the remaining 40%. Additionally the system included no-till, cover crops, atrazine split-applications based on weed pressure, variable N rates, and variable rates of fall-applied P. As expected, the system resulted in greatly reduced sediment losses (87% less). While plot-scale experimental and modeling studies suggested that no-till would result in increased herbicide and dissolved P losses, the assessment showed lower losses than expected. Overall N losses decreased by 40% but NO3-N and NH4-N behaved differently: NO3-N losses decreased by ~57% while NH4-N losses increased by 47%. Conclusions from plot-scale experiments were consistent with these results: increased adsorption caused by gains in soil organic C at the soil surface resulting from no-till and cover crops, adaptations of microorganisms to atrazine, and the wheat area functioning as a buffer for herbicides applied at the top end of the field. Other not-specifically tested explanations included fertilizer applications when runoff risk was lower, and uptake of nutrients by the cover crops. These results demonstrate the benefits of a comprehensive conservation plan that includes multiple practices, some counteracting the negative aspects of others.
The results obtained during the 11 years of the PAS and the farming experience gained led to the proposal for an aspirational system: a three-year no-till corn–soybean–wheat with cover crops, variable rates of fertilizers, split applications of atrazine, incorporation of agrichemical inputs with a rotary harrow, 1.6 ha of perennial grasses at the bottom of the field functioning as a buffer, and subsurface N application. Some of these aspects have been implemented and future analyses will determine their benefits.
Acknowledgements
The authors sincerely thank numerous USDA and University of Missouri personnel for assisting in the establishment, maintenance, and management of the field used in this study, in the associated monitoring infrastructure, and in the collection and chemical analysis of water samples. We gratefully acknowledge that partial funding for this research was provided by the USDA Natural Resources Conservation Service Conservation Effects Assessment Project Watershed Assessment Studies and the USDA–CSREES (now National Institute for Food and Agriculture) National Integrated Water Quality Program–Conservation Effects Assessment Project (CEAP) (Grant no. 2005-51130-02380). This research was also supported by the USDA Agricultural Research Service (ARS) Long-Term Agroecosystem Research network, and ARS National Program 211.
Footnotes
Disclaimer
Mention of trade names or commercial products in this article is solely for the purpose of providing specific information and does not imply recommendation or endorsement by the USDA.
- © 2020 by the Soil and Water Conservation Society