Applications of the SWAT Model Special Section: Overview and Insights.

The Soil and Water Assessment Tool (SWAT) model has emerged as one of the most widely used water quality watershed- and river basin-scale models worldwide, applied extensively for a broad range of hydrologic and/or environmental problems. The international use of SWAT can be attributed to its flexibility in addressing water resource problems, extensive networking via dozens of training workshops and the several international conferences that have been held during the past decade, comprehensive online documentation and supporting software, and an open source code that can be adapted by model users for specific application needs. The catalyst for this special collection of papers was the 2011 International SWAT Conference & Workshops held in Toledo, Spain, which featured over 160 scientific presentations representing SWAT applications in 37 countries. This special collection presents 22 specific SWAT-related studies, most of which were presented at the 2011 SWAT Conference; it represents SWAT applications on five different continents, with the majority of studies being conducted in Europe and North America. The papers cover a variety of topics, including hydrologic testing at a wide range of watershed scales, transport of pollutants in northern European lowland watersheds, data input and routing method effects on sediment transport, development and testing of potential new model algorithms, and description and testing of supporting software. In this introduction to the special section, we provide a synthesis of these studies within four main categories: (i) hydrologic foundations, (ii) sediment transport and routing analyses, (iii) nutrient and pesticide transport, and (iv) scenario analyses. We conclude with a brief summary of key SWAT research and development needs.

Sediment delivery estimates in water quality models altered by resolution and source of topographic data.

Moderate-resolution (30-m) digital elevation models (DEMs) are normally used to estimate slope for the parameterization of non-point source, process-based water quality models. These models, such as the Soil and Water Assessment Tool (SWAT), use the Universal Soil Loss Equation (USLE) and Modified USLE to estimate sediment loss. The slope length and steepness factor, a critical parameter in USLE, significantly affects sediment loss estimates. Depending on slope range, a twofold difference in slope estimation potentially results in as little as 50% change or as much as 250% change in the LS factor and subsequent sediment estimation. Recently, the availability of much finer-resolution (∼3 m) DEMs derived from Light Detection and Ranging (LiDAR) data has increased. However, the use of these data may not always be appropriate because slope values derived from fine spatial resolution DEMs are usually significantly higher than slopes derived from coarser DEMs. This increased slope results in considerable variability in modeled sediment output. This paper addresses the implications of parameterizing models using slope values calculated from DEMs with different spatial resolutions (90, 30, 10, and 3 m) and sources. Overall, we observed over a 2.5-fold increase in slope when using a 3-m instead of a 90-m DEM, which increased modeled soil loss using the USLE calculation by 130%. Care should be taken when using LiDAR-derived DEMs to parameterize water quality models because doing so can result in significantly higher slopes, which considerably alter modeled sediment loss.

Mapping landscape-level hydrological connectivity of headwater wetlands to downstream waters: A geospatial modeling approach - Part 1.

Headwater wetlands affect ecosystem integrity of downstream waters; however, many wetlands - particularly geographically isolated wetlands (GIWs) - continue to be at risk. A significant portion of US federal policy is based on the jurisdictional status of wetlands, which is partly determined by the relationship between wetlands and downstream waters, including the cumulative impact of wetlands on those waters. We present a novel multi-phase geospatial modeling method to help elucidate hydrological relationship between GIWs and downstream waters at the landscape scale. The presented approach in this study used inundation maps derived from time series remotely sensed data between 1985 and 2010, weather and hydrological records, and ancillary geospatial data including information from the US Fish and Wildlife Service National Wetlands Inventory (NWI). The study site was a headwater catchment (292 km2) of the Choptank River Basin, located in the Mid-Atlantic region of USA, which contained a large number of Delmarva bays. The results showed inundation extent within GIWs varied, in aggregate, in response to weather variability (r = 0.58; p-value = 0.05), and was well correlated with streamflow (r = 0.81; p-value < 0.01) and base flow (r = 0.57; p-value < 0.1) conditions. The relationship between inundation patterns and stream discharge also varied with NWI hydrologic modifiers. The GIWs with water regime characterized by longer durations of flooding exhibited stronger correlations with stream discharge, but those GIWs with shorter durations of flooding were less correlated with stream discharge. This analysis suggests the mutual reliance (i.e., connection) of wetlands and streams on groundwater. GIWs appeared to function in aggregate, and it is likely that the combined effect of these wetlands significantly influenced the functioning of downstream waters.

Mapping landscape-level hydrological connectivity of headwater wetlands to downstream waters: A catchment modeling approach - Part 2.

In Part 1 of this two-part manuscript series, we presented an effective assessment method for mapping inundation of geographically isolated wetlands (GIWs) and quantifying their cumulative landscape-scale hydrological connectivity with downstream waters using time series remotely sensed data (Yeo et al., 2018). This study suggested strong hydrological coupling between GIWs and downstream waters at the seasonal timescale via groundwater. This follow-on paper investigates the hydrological connectivity of GIWs with downstream waters and cumulative watershed-scale hydrological impacts over multiple time scales. Modifications were made to the representation of wetland processes within the Soil and Water Assessment Tool (SWAT). A version of SWAT with improved wetland function, SWAT-WET, was applied to Greensboro Watershed, which is located in the Mid-Atlantic Region of USA, to simulate hydrological processes over 1985-2015 under two contrasting land use scenarios (i.e., presence and absence of GIWs). Comparative analysis of simulation outputs elucidated how GIWs could influence partitioning of precipitation between evapotranspiration (ET) and terrestrial water storage, and affect water transport mechanisms and routing processes that generate streamflow. Model results showed that GIWs influenced the watershed water budget and stream flow generation processes over the long-term (30 year), inter-annual, and monthly time scales. GIWs in the study watershed increased terrestrial water storage during the wet season, and buffered the dynamics of shallow groundwater during the dry season. The inter-annual modeling analysis illustrated that densely distributed GIWs can exert strong hydrological influence on downstream waters by regulating surface water runoff, while maintaining groundwater recharge and ET under changing (wetter) climate conditions. The study findings highlight the hydrological connectivity of GIWs with downstream waters and the cumulative hydrological influence of GIWs as hydrologic sources to downstream ecosystems through different runoff processes over multiple time scales.

Reducing insecticide and fungicide loads in runoff from plastic mulch with vegetative-covered furrows.

A common management practice for the production of fresh-market vegetables utilizes polyethylene (plastic) mulch because it increases soil temperature, decreases weed pressure, maintains soil moisture, and minimizes soil contact with the product. However, rain events afford much more erosion and runoff because 50-75% of the field is covered with an impervious surface. A plot study was conducted to compare and to quantify the off-site movement of soil, insecticides, and fungicides associated with runoff from plots planted with Sunbeam tomatoes (Lycopersicon esculentum Mill) using the conventional polyethylene mulch management practice vs an alternative management practice-polyethylene mulch-covered beds with cereal rye (Secale cereale) planted in the furrows between the beds. The use of cereal rye-covered furrows with the conventional polyethylene system decreased runoff volume by more than 40%, soil erosion by more than 80%, and pesticide loads by 48-74%. Results indicate that vegetative furrows are critical to minimizing the negative aspects of this management practice.

Predicting perchlorate exposure in milk from concentrations in dairy feed.

Perchlorate has been detected in U.S. milk samples from many different states. Applying data from a recently reported 9-week experiment in which 16 Holstein dairy cows were administered perchlorate allowed us to derive an equation for the dose-response relationship between perchlorate concentrations in feed/drinking water and its appearance in milk. Examination of background concentrations of perchlorate in the total mixed ration (TMR) fed in addition to the variable dose supplied to treated cows as a ruminal infusate revealed that cows receive significant and variable exposure to perchlorate from the TMR. Weekly examination of the TMR disclosed that a change in ingredients midway through the experiment caused a significant (78%) change in TMR perchlorate concentration. Analyses of the ingredients comprising the TMR revealed that 41.9% of the perchlorate came from corn silage, 22.9% came from alfalfa hay and 11.7% was supplied by sudan grass. Finally, USDA Food and Nutrition Survey data on fluid milk consumption were used to predict potential human exposure from milk that contained concentrations of perchlorate observed in our previous dosing study. The study suggests that reducing perchlorate concentration in dairy feed may reduce perchlorate concentrations in milk as well as the potential to reduce human exposure to perchlorate in milk.

Evaluation of vegetable production management practices to reduce the ecological risk of pesticides.

The ability of agricultural management practices to reduce the ecological risks of pesticides was evaluated. Risk quotients, a mathematical description of the relationship between exposure and toxicity, and hazard ratings, a rank of the potential risk of pesticides to aquatic environments, were calculated for conventional and alternative cultivation practices for tomatoes: Poly-Bare, raised beds covered with polyethylene mulch with bare-soil furrows; Poly-Rye, raised beds covered with polyethylene mulch with cereal rye (Secale cereale) grown in the furrows; and Vetch, raised beds and furrows planted with hairy vetch seed (Vicia villosa). Evaluations were conducted using measured pesticide concentrations in runoff at the edge-of-field and estimated environmental concentrations in an adjacent creek and a theoretical pond receiving the runoff. Runoff from Poly-Bare presented the greatest risk to ecosystem health and to sensitive organisms, whereas the use of Vetch minimized these risks. Previous studies have shown that harvest yields were maintained and that runoff volume, soil loss, and off-site transport of pesticides measured in runoff were reduced using the alternative management practices (Poly-Rye and Vetch). Together, these results indicate that the alternative management practices (Poly-Rye and Vetch) have a less adverse impact on the environment than the conventional management practice (Poly-Bare) while providing growers with an acceptable economic return. In addition, the present study demonstrates the need to consider the management practice when assessing the potential risks and hazards for certain pesticides.

Impacts of Watershed Characteristics and Crop Rotations on Winter Cover Crop Nitrate-Nitrogen Uptake Capacity within Agricultural Watersheds in the Chesapeake Bay Region.

The adoption rate of winter cover crops (WCCs) as an effective conservation management practice to help reduce agricultural nutrient loads in the Chesapeake Bay (CB) is increasing. However, the WCC potential for water quality improvement has not been fully realized at the watershed scale. This study was conducted to evaluate the long-term impact of WCCs on hydrology and NO3-N loads in two adjacent watersheds and to identify key management factors that affect the effectiveness of WCCs using the Soil and Water Assessment Tool (SWAT) and statistical methods. Simulation results indicated that WCCs are effective for reducing NO3-N loads and their performance varied based on planting date, species, soil characteristics, and crop rotations. Early-planted WCCs outperformed late-planted WCCs on the reduction of NO3-N loads and early-planted rye (RE) reduced NO3-N loads by ~49.3% compared to the baseline (no WCC). The WCCs were more effective in a watershed dominated by well-drained soils with increased reductions in NO3-N fluxes of ~2.5 kg N·ha-1 delivered to streams and ~10.1 kg N·ha-1 leached into groundwater compared to poorly-drained soils. Well-drained agricultural lands had higher transport of NO3-N in the soil profile and groundwater due to increased N leaching. Poorly-drained agricultural lands had lower NO3-N due to extensive drainage ditches and anaerobic soil conditions promoting denitrification. The performance of WCCs varied by crop rotations (i.e., continuous corn and corn-soybean), with increased N uptake following soybean crops due to the increased soil mineral N availability by mineralization of soybean residue compared to corn residue. The WCCs can reduce N leaching where baseline NO3-N loads are high in well-drained soils and/or when residual and mineralized N availability is high due to the cropping practices. The findings suggested that WCC implementation plans should be established in watersheds according to local edaphic and agronomic characteristics for reducing N leaching.

Pesticide occurrence in selected South Florida canals and Biscayne Bay during high agricultural activity.

Climate and soil conditions in South Florida along with an extensive canal system facilitate movement of agricultural pesticides into surface waters. In a two-year study (2002-2004) of the currently used pesticides in South Florida, atrazine, endosulfan, metolachlor, chlorpyrifos, and chlorothalonil were the most frequently detected in the canals and in Biscayne Bay, with average concentrations of 16, 11, 9.0, 2.6, and 6.0 ng/L, respectively. Concentrations of atrazine and chlorpyrifos were highest near corn production. Chlorothalonil and endosulfan concentrations were highest near vegetable production, with no clear trend for metolachlor, which is used on multiple crops. Concentration data were used to calculate an aquatic life hazard potential for the planting period (November) versus the harvest period (March). This analysis indicated that a higher hazard potential occurs during harvest, primarily from the use of endosulfan. These data will also serve to document canal conditions prior to implementation of the Comprehensive Everglades Restoration Plan (CERP).

Atrazine distribution measured in soil and leachate following infiltration conditions.

Atrazine transport through packed 10 cm soil columns representative of the 0-10 cm soil horizon was observed by measuring the atrazine recovery in the total leachate volume, and upper and lower soil layers following infiltration of 7.5 cm water using a mechanical vacuum extractor (MVE). Measured recoveries were analyzed to understand the influence of infiltration rate and delay time on atrazine transport and distribution in the column. Four time periods (0.28, 0.8, 1.8, and 5.5 h) representing very high to moderate infiltration rates (26.8, 9.4, 4.2, and 1.4 cm/h) were used. Replicate soil columns were tested immediately and following a 2-d delay after atrazine application. Results indicate atrazine recovery in leachate was independent of infiltration rate, but significantly lower for infiltration following a 2-d delay. Atrazine distribution in the 0-1 and 9-10 cm soil layers was affected by both infiltration rate and delay. These results are in contrast with previous field and laboratory studies that suggest that atrazine recovery in the leachate increases with increasing infiltration rate. It appears that the difference in atrazine recovery measured using the MVE and other leaching experiments using intact soil cores from this field site and the rain simulation equipment probably illustrates the effect of infiltrating water interacting with the atrazine present on the soil surface. This work suggests that atrazine mobilization from the soil surface is also dependent on interactions of the infiltrating water with the soil surface, in addition to the rate of infiltration through the surface soil.

Comparison of copper levels in runoff from fresh-market vegetable production using polyethylene mulch or a vegetative mulch.

Runoff from tomato (Lycopersicon esculentum Mill.) production with polyethylene mulch has been implicated in the failure of commercial shellfish farms in the Mid-Atlantic Region of the United States. Copper, applied in the form of copper hydroxide, is the most widely used fungicide-bactericide for control of tomato diseases and recently has been detected in the Chesapeake Bay (USA) watershed. Elevated levels of copper have been shown to have adverse effects on shellfish, finfish, and other aquatic organisms. This research evaluates the off-site movement of copper with the dissolved phase and the particulate phase of runoff from controlled field plots containing tomato plants grown in either polyethylene mulch or a vegetative mulch, hairy vetch (Vicia villosa Roth.). Overall, runoff collected from polyethylene mulch plots contained significantly (p < or = 0.05) greater loads of dissolved- and particulate-phase copper than runoff from hairy vetch mulch plots. However, the loss of copper associated with the particulate phase was significantly greater (p < or = 0.05) than that associated with the dissolved phase of runoff from both mulch treatments, with the particulate phase accounting for more than 80% of the copper loads. The reported toxicity of copper to aquatic organisms and the greater runoff volume, soil loss, and off-site loading of copper measured in runoff from the polyethylene mulch suggests that this management practice is less sustainable and may have a more harmful impact on aquatic ecosystems.

Use of vegetative furrows to mitigate copper loads and soil loss in runoff from polyethylene (plastic) mulch vegetable production systems.

The transport of runoff with high copper concentrations and sediment loads into adjacent surface waters can have adverse effects on nontarget organisms as a result of increased turbidity and degraded water quality. Runoff from vegetable production utilizing polyethylene mulch can contain up to 35% of applied copper, a widely used fungicide/bactericide that has adverse effects on aquatic organisms. Copper is primarily transported in runoff with suspended particulates; therefore, implementation of management practices that minimize soil erosion will reduce copper loads. Replacing bare-soil furrows with furrows planted in rye (Secale cereale) significantly improved the sustainability of vegetable production with polyethylene mulch and reduced the potential environmental impact of this management practice. Vegetative furrows decreased runoff volume by >40% and soil erosion by >80%. Copper loads with runoff were reduced by 72% in 2001, primarily as a result of reduced soil erosion since more than 88% of the total copper loads were transported in runoff with suspended soil particulates. Tomato yields in both years were similar between the polyethylene mulch plots containing either bare-soil or vegetative furrows. Replacing bare-soil furrows with vegetative furrows greatly reduces the effects of sediments and agrochemicals on sensitive ecosystems while maintaining crop yields.

Current United States Department of Agriculture-Agricultural Research Service research on understanding agrochemical fate and transport to prevent and mitigate adverse environmental impacts.

Environmentally and economically viable agriculture requires a variety of cultivation practices and pest management options as no one system will be appropriate for every situation. Agrochemicals are some of the many pest control tools used in an integrated approach to pest management. They are applied with the intent of maximizing efficacy while minimizing off-site movement; however, their judicious use demands a practical knowledge of their fate and effects in agricultural and natural ecosystems. Agrochemical distribution into environmental compartments is influenced by the physical and chemical properties of the agrochemical and environmental conditions, ie soil type and structure, and meteorological conditions. Agricultural Research Service (ARS) researchers working in the area of agrochemical fate have focused on accurately describing those processes that govern the transport, degradation and bioavailability of these chemicals under conditions reflecting actual agronomic practices. Results from ARS research concerning the environmental fate and effects of agrochemicals have led to the development of science-based management practices that will protect vulnerable areas of the ecosystem. The new challenge is to identify these vulnerable areas and the temporal and spatial variations prior to use of the chemical by predicting how it will behave in environmental matrices, and using that information, predict its transport and transformation within an air- or watershed. With the development of better predictive tools and GIS (Geographic Information System)-based modeling, the risks of agricultural management systems can be assessed at the watershed and basin levels, and management strategies can be identified that minimize negative environmental impacts.

Metolachlor metabolite (MESA) reveals agricultural nitrate-N fate and transport in Choptank River watershed.

Over 50% of streams in the Chesapeake Bay watershed have been rated as poor or very poor based on the index of biological integrity. The Choptank River estuary, a Bay tributary on the eastern shore, is one such waterway, where corn and soybean production in upland areas of the watershed contribute significant loads of nutrients and sediment to streams. We adopted a novel approach utilizing the relationship between the concentration of nitrate-N and the stable, water-soluble herbicide degradation product MESA {2-[2-ethyl-N-(1-methoxypropan-2-yl)-6-methylanilino]-2-oxoethanesulfonic acid} to distinguish between dilution and denitrification effects on the stream concentration of nitrate-N in agricultural subwatersheds. The ratio of mean nitrate-N concentration/(mean MESA concentration * 1000) for 15 subwatersheds was examined as a function of percent cropland on hydric soil. This inverse relationship (R(2)=0.65, p<0.001) takes into consideration not only dilution and denitrification of nitrate-N, but also the stream sampling bias of the croplands caused by extensive drainage ditch networks. MESA was also used to track nitrate-N concentrations within the estuary of the Choptank River. The relationship between nitrate-N and MESA concentrations in samples collected over three years was linear (0.95 ≤ R(2) ≤ 0.99) for all eight sampling dates except one where R(2)=0.90. This very strong correlation indicates that nitrate-N was conserved in much of the Choptank River estuary, that dilution alone is responsible for the changes in nitrate-N and MESA concentrations, and more importantly nitrate-N loads are not reduced in the estuary prior to entering the Chesapeake Bay. Thus, a critical need exists to minimize nutrient export from agricultural production fields and to identify specific conservation practices to address the hydrologic conditions within each subwatershed. In well drained areas, removal of residual N within the cropland is most critical, and practices such as cover crops which sequester the residual N should be strongly encouraged. In poorly drained areas where denitrification can occur, wetland restoration and controlled drained structures that minimize ditch flow should be used to maximize denitrification.

Relating nutrient and herbicide fate with landscape features and characteristics of 15 subwatersheds in the Choptank River watershed.

Excess nutrients and agrochemicals from non-point sources contribute to water quality impairment in the Chesapeake Bay watershed and their loading rates are related to land use, agricultural practices, hydrology, and pollutant fate and transport processes. In this study, monthly baseflow stream samples from 15 agricultural subwatersheds of the Choptank River in Maryland USA (2005 to 2007) were characterized for nutrients, herbicides, and herbicide transformation products. High-resolution digital maps of land use and forested wetlands were derived from remote sensing imagery. Examination of landscape metrics and water quality data, partitioned according to hydrogeomorphic class, provided insight into the fate, delivery, and transport mechanisms associated with agricultural pollutants. Mean Nitrate-N concentrations (4.9 mg/L) were correlated positively with percent agriculture (R(2)=0.56) and negatively with percent forest (R(2)=0.60). Concentrations were greater (p=0.0001) in the well-drained upland (WDU) hydrogeomorphic region than in poorly drained upland (PDU), reflecting increased denitrification and reduced agricultural land use intensity in the PDU landscape due to the prevalence of hydric soils. Atrazine and metolachlor concentrations (mean 0.29 µg/L and 0.19 µg/L) were also greater (p=0.0001) in WDU subwatersheds than in PDU subwatersheds. Springtime herbicide concentrations exhibited a strong, positive correlation (R(2)=0.90) with percent forest in the WDU subwatersheds but not in the PDU subwatersheds. In addition, forested riparian stream buffers in the WDU were more prevalent than in the PDU where forested patches are typically not located near streams, suggesting an alternative delivery mechanism whereby volatilized herbicides are captured by the riparian forest canopy and subsequently washed off during rainfall. Orthophosphate, CIAT (6-chloro-N-(1-methylethyl)-1,3,5-triazine-2,4-diamine), CEAT (6-chloro-N-ethyl-1,3,5-triazine-2,4-diamine), and MESA (2-[(2-ethyl-6-methylphenyl) (2-methoxy-1-methylethyl)amino]-2-oxoethanesulfonic acid) were also analyzed. These findings will assist efforts in targeting implementation of conservation practices to the most environmentally-critical areas within watersheds to achieve water quality improvements in a cost-effective manner.

Pollutant fate and spatio-temporal variability in the choptank river estuary: factors influencing water quality.

Restoration of the Chesapeake Bay, the largest estuary in the United States, is a national priority. Documentation of progress of this restoration effort is needed. A study was conducted to examine water quality in the Choptank River estuary, a tributary of the Chesapeake Bay that since 1998 has been classified as impaired waters under the Federal Clean Water Act. Multiple water quality parameters (salinity, temperature, dissolved oxygen, chlorophyll a) and analyte concentrations (nutrients, herbicide and herbicide degradation products, arsenic, and copper) were measured at seven sampling stations in the Choptank River estuary. Samples were collected under base flow conditions in the basin on thirteen dates between March 2005 and April 2008. As commonly observed, results indicate that agriculture is a primary source of nitrate in the estuary and that both agriculture and wastewater treatment plants are important sources of phosphorus. Concentrations of copper in the lower estuary consistently exceeded both chronic and acute water quality criteria, possibly due to use of copper in antifouling boat paint. Concentrations of copper in the upstream watersheds were low, indicating that agriculture is not a significant source of copper loading to the estuary. Concentrations of herbicides (atrazine, simazine, and metolachlor) peaked during early-summer, indicating a rapid surface-transport delivery pathway from agricultural areas, while their degradation products (CIAT, CEAT, MESA, and MOA) appeared to be delivered via groundwater transport. Some in-river processing of CEAT occurred, whereas MESA was conservative. Observed concentrations of herbicide residues did not approach established levels of concern for aquatic organisms. Results of this study highlight the importance of continued implementation of best management practices to improve water quality in the estuary. This work provides a baseline against which to compare future changes in water quality and may be used to design future monitoring programs needed to assess restoration strategy efficacy.