Modeling Discharge Rates Using a Coupled Modeled Approach for the Merced River in Yosemite National Park

This study describes the application of the NASA version of the Carnegie‐Ames‐Stanford Approach (CASA) ecosystem model coupled with a surface hydrologic routing scheme previously called the Hydrological Routing Algorithm (HYDRA) to model monthly discharge rates from 2000 to 2007 on the Merced River drainage in Yosemite National Park, California. To assess CASA‐HYDRA's capability to estimate actual water flows in extreme precipitation years, the focus of this study is the 2007 water year, which was very dry, and the 2005 water year, which was a moderately wet year in the historical record. Prior to comparisons to gauge records, CASA‐HYDRA snowmelt algorithms were modified with equations from the U.S. Department of Agriculture Snowmelt‐Runoff Model (SRM), which has been designed to predict daily streamflow in mountain basins where snowmelt is a major runoff factor. Results show that model predictions closely matched monthly flow rates at the Pohono Bridge gauge station (USGS#11266500), with R² = 0.67 and Nash‐Sutcliffe (E) = 0.65. By subdividing the upper Merced River basin into subbasins with high spatial resolution in the gridded modeling approach, we were able to determine which biophysical characteristics in the Sierra differed to the largest degree in extreme low‐flow and high‐flow years. Average elevation and snowpack accumulation were found to be the most important explanatory variables to understand subbasin contributions to monthly discharge rates.     

HYDROLOGICAL MODELING OF THE IROQUOIS RIVER WATERSHED USING HSPF AND SWAT1

ABSTRACT: The performance of two popular watershed scale simulation models — HSPF and SWAT — were evaluated for simulating the hydrology of the 5,568 km² Iroquois River watershed in Illinois and Indiana. This large, tile drained agricultural watershed provides distinctly different conditions for model comparison in contrast to previous studies. Both models were calibrated for a nine‐year period (1987 through 1995) and verified using an independent 15‐year period (1972 through 1986) by comparing simulated and observed daily, monthly, and annual streamflow. The characteristics of simulated flows from both models are mostly similar to each other and to observed flows, particularly for the calibration results. SWAT predicts flows slightly better than HSPF for the verification period, with the primary advantage being better simulation of low flows. A noticeable difference in the models' hydrologic simulation relates to the estimation of potential evapotranspiration (PET). Comparatively low PET values provided as input to HSPF from the BASINS 3.0 database may be a factor in HSPF's overestimation of low flows. Another factor affecting baseflow simulation is the presence of tile drains in the watershed. HSPF parameters can be adjusted to indirectly account for the faster subsurface flow associated with tile drains, but there is no specific tile drainage component in HSPF as there is in SWAT. Continued comparative studies such as this, under a variety of hydrologic conditions and watershed scales, provide needed guidance to potential users in model selection and application.     

Spatial Calibration and Temporal Validation of Flow for Regional Scale Hydrologic Modeling1

Abstract: Physically based regional scale hydrologic modeling is gaining importance for planning and management of water resources. Calibration and validation of such regional scale model is necessary before applying it for scenario assessment. However, in most regional scale hydrologic modeling, flow validation is performed at the river basin outlet without accounting for spatial variations in hydrological parameters within the subunits. In this study, we calibrated the model to capture the spatial variations in runoff at subwatershed level to assure local water balance, and validated the streamflow at key gaging stations along the river to assure temporal variability. Ohio and Arkansas‐White‐Red River Basins of the United States were modeled using Soil and Water Assessment Tool (SWAT) for the period from 1961 to 1990. R² values of average annual runoff at subwatersheds were 0.78 and 0.99 for the Ohio and Arkansas Basins. Observed and simulated annual and monthly streamflow from 1961 to 1990 is used for temporal validation at the gages. R² values estimated were greater than 0.6. In summary, spatially distributed calibration at subwatersheds and temporal validation at the stream gages accounted for the spatial and temporal hydrological patterns reasonably well in the two river basins. This study highlights the importance of spatially distributed calibration and validation in large river basins.     

Nutrient retention in tile-fed and non-tile reconstructed oxbows in north central Iowa

Nutrient export from the agricultural Midwest threatens the Gulf of Mexico and new conservation practices are needed to reduce the loss of nutrient from subsurface tile drainage systems. Oxbows are natural waterbodies formed when a river cuts off a meander loop and water quality benefits of reconstructed oxbows are being increasingly recognized. In this study, we monitored four reconstructed oxbow sites (two tile-fed, two non-tile) over a 2-year period in north-central Iowa and assessed their capacity for NO₃-N and dissolved reactive phosphorus (DRP) reductions. Water flow and quality monitoring of tiles, shallow groundwater, oxbow and receiving streams documented that the oxbows were dominated by tile drainage inputs. NO₃-N concentrations were highest in the drainage tiles flowing into the tile-fed oxbows (mean 8–10 mg/L) and much lower in floodplain groundwater (<1–2 mg/L). Annual NO₃-N loads into the tile-fed oxbows were substantially larger than input loads into the non-tiled oxbows. For the two tile-fed oxbows, the 2-year NO₃-N retention efficiencies were very similar (0.76–0.77) and on a monthly basis, greater retention efficiencies were measured in summer and fall. DRP concentrations and loads into the tile-fed oxbows were too low to allow for meaningful estimates of retention. Reconstructing oxbows to receive tile drainage water should be considered a sustainable conservation practice for tile drainage treatment in agricultural areas.     

Modeling of Sediment Yield From Anjeni-Gauged Watershed,Ethiopia Using SWAT Model1

Setegn, Shimelis G., Bijan Dargahi, Ragahavan Srinivasan, and Assefa M. Melesse, 2010. Modeling of Sediment Yield From Anjeni-Gauged Watershed, Ethiopia Using SWAT Model. Journal of the American Water Resources Association (JAWRA) 46(3):514-526. DOI: 10.1111/j.1752-1688.2010.00431.x Abstract: The Soil and Water Assessment Tool (SWAT) was tested for prediction of sediment yield in Anjeni-gauged watershed, Ethiopia. Soil erosion and land degradation is a major problem on the Ethiopian highlands. The objectives of this study were to evaluate the performance and applicability of SWAT model in predicting monthly sediment yield and assess the impacts of subbasin delineation and slope discretization on the prediction of sediment yield. Ten years monthly meteorological, flow and sediment data were used for model calibration and validation. The annual average measured sediment yield was 24.6 tonnes/ha. The annual average simulated sediment yield was 27.8 and 29.5 tones/ha for calibration and validation periods, respectively. The study found that the observed values showed good agreement with the simulated sediment yield with Nash-Sutcliffe efficiency (NSE) = 0.81, percent bias (PBIAS) = 28%, RMSE-observations standard deviation ratio (RSR) = 0.23, and coefficient of determination (R²) = 0.86 for calibration and NSE = 0.79, PBIAS = 30%, RSR = 0.29, and R² = 0.84 for validation periods. The model can be used for further analysis of different management scenarios that could help different stakeholders to plan and implement appropriate soil and water conservation strategies.     

Simulation of Daily Flow Pathways,Tile‐Drain Nitrate Concentrations,and Soil‐Nitrogen Dynamics Using SWAT

Tile drainage significantly alters flow and nutrient pathways and reliable simulation at this scale is needed for effective planning of nutrient reduction strategies. The Soil and Water Assessment Tool (SWAT) has been widely utilized for prediction of flow and nutrient loads, but few applications have evaluated the model's ability to simulate pathway‐specific flow components or nitrate‐nitrogen (NO₃‐N) concentrations in tile‐drained watersheds at the daily time step. The objectives of this study were to develop and calibrate SWAT models for small, tile‐drained watersheds, evaluate model performance for simulation of flow components and NO₃‐N concentration at daily intervals, and evaluate simulated soil‐nitrogen dynamics. Model evaluation revealed that it is possible to meet accepted performance criteria for simulation of monthly total flow, subsurface flow (SSF), and NO₃‐N loads while obtaining daily surface runoff (SURQ), SSF, and NO₃‐N concentrations that are not satisfactory. This limits model utility for simulating best management practices (BMPs) and compliance with water quality standards. Although SWAT simulates the soil N‐cycle and most predicted fluxes were within ranges reported in agronomic studies, improvements to algorithms for soil‐N processes are needed. Variability in N fluxes is extreme and better parameterization and constraint, through use of more detailed agronomic data, would also improve NO₃‐N simulation in SWAT. Editor's note : This paper is part of the featured series on SWAT Applications for Emerging Hydrologic and Water Quality Challenges. See the February 2017 issue for the introduction and background to the series.     

Estimation of Nonpoint Source Nitrate Concentrations in Indiana Rivers Based on Agricultural Drainage in the Watershed

Subsurface tile‐drained agricultural fields are known to be important contributors to nitrate in surface water in the Midwest, but the effect of these fields on nitrate at the watershed scale is difficult to quantify. Data for 25 watersheds monitored by the Indiana Department of Environmental Management and located near a U.S. Geological Survey stream gage were used to investigate the relationship between flow‐weighted mean concentration (FWMC) of nitrate‐N and the subsurface tile‐drained area (DA) of the watershed. The tile DA was estimated from soil drainage class, land use, and slope. Nitrate loads from point sources were estimated based on reported flows of major permitted facilities with mean nitrate‐N concentrations from published sources. Linear regression models exhibited a statistically significant relationship between annual/monthly nonpoint source (NPS) nitrate‐N and DA percentage. The annual model explained 71% of the variation in FWMC of nitrate‐N. The annual and monthly models were tested in 10 additional watersheds, most with absolute errors within 1 mg/l in the predicted FWMC. These models can be used to estimate NPS nitrate for unmonitored watersheds in similar areas, especially for drained agricultural areas where model performance was strongest, and to predict the nitrate reduction when various tile drainage management techniques are employed.     

Model Based Nitrate TMDLs for Two Agricultural Watersheds of Southeastern Minnesota1

Abstract: In this study, a set of nitrogen reduction strategies were modeled to evaluate the feasibility of improving water quality to meet total maximum daily loads (TMDLs) in two agricultural watersheds. For this purpose, a spatial‐process model was calibrated and used to predict monthly nitrate losses (1994‐96) from Sand and Bevens Creek watersheds located in south‐central Minnesota. Statistical comparison of predicted and observed flow and nitrate losses gave r² coefficients of 0.75 and 0.70 for Sand Creek watershed and 0.72 and 0.67 for Bevens Creek watershed, respectively. Modeled alternative agricultural management scenarios included: six different N application rates over three application timings and three different percentages of crop land with subsurface drainage. Predicted annual nitrate losses were then compared with nitrate TMDLs assuming a 30% reduction in observed nitrate losses is required. Reductions of about 33 (8.6 to 5.8 kg/ha) and 35% (23 to 15 kg/ha) in existing annual nitrate losses are possible for Sand and Bevens Creek watersheds, respectively, by switching the timing of fertilizer application from fall to spring. Trends towards increases in tile‐drained crop land imply that attaining nitrate TMDLs in future may require other alternative management practices in addition to fertilizer management such as partial conversion of crop land to pasture.     

Integrating Temporal Inequality into Conservation Planning to Improve Practice Design and Efficacy

In contrast to spatial inequality, there are currently no methods for leveraging information on temporal inequality to improve conservation efficacy. The objective of this study was to use Lorenz curves to quantify temporal inequality in surface runoff and tile drainage, identify controls on nutrient loading in these flowpaths, and develop design flows for structural conservation practices. Surface runoff (n = 94 site‐years) and tile drainage (n = 90 site‐years) were monitored on 40 fields in Ohio. Results showed, on average, 80% of nitrate‐nitrogen, soluble reactive phosphorus (P), and total P loads occurred between 7 and 12 days per year in surface runoff and between 32 and 58 days per year in tile drainage. Similar temporal inequality between discharge and load provided evidence that loading was transport‐limited and highlighted the critical role hydrologic connectivity plays in nutrient delivery from tile‐drained fields. Design flow criterion for sizing structural practices based on load reduction goals was developed by combining Lorenz curves and flow duration curves. Comparing temporal inequality between fields and the Maumee River, the largest tributary to the western Lake Erie Basin, revealed challenges associated with achieving watershed load reduction goals with field‐scale conservation. In‐field (i.e., improved nutrient and water management), edge‐of‐field (i.e., structural practices), and instream practices will all be required to meet nutrient reduction goals from tile‐drained watersheds.     

Calibration and Validation of ADAPT and SWAT for Field‐Scale Runoff Prediction1

Abstract: The pollutant reduction possible with a given agricultural best‐management practice (BMP) is complex and site‐specific. Water‐quality models can evaluate BMPs, but model results are often limited by the lack of calibrated parameters for a given BMP. This study calibrated runoff prediction of two models (ADAPT and SWAT) for individual field plots having one till and two no‐till management practices. The factors used for runoff calibration were curve number II (CN_II) and saturated hydraulic conductivity (Ksat) for ADAPT, and CN_II, Ksat, and available water capacity for SWAT. Results were evaluated using coefficient of determination (R²), Nash‐Sutcliffe efficiency (E_f), root‐mean square error, median‐based E_f, and sign tests. Results indicated that for ADAPT, the best‐fit CN_II was 66 for the NT/SB (no‐till plot with surface‐broadcast fertilizer) treatment, 68 for the NT/DB (no‐till with deep‐banded fertilizer) treatment, and 70 for the tilled plot, whereas for SWAT the best‐fit CN_II was much higher, 86, for all treatments. Neither agreed with the textbook CN_II, 78, for sorghum in silty clay loam soil. The best‐fit model parameters for both runoff calibration phases had excellent correlation to monthly totals and moderate correlation to individual events.     

Modeling Parent and Metabolite Fate and Transport in Subsurface Drained Fields With Directly Connected Macropores1

Abstract: Few studies exist that evaluate or apply pesticide transport models based on measured parent and metabolite concentrations in fields with subsurface drainage. Furthermore, recent research suggests pesticide transport through exceedingly efficient direct connections, which occur when macropores are hydrologically connected to subsurface drains, but this connectivity has been simulated at only one field site in Allen County, Indiana. This research evaluates the Root Zone Water Quality Model (RZWQM) in simulating the transport of a parent compound and its metabolite at two subsurface drained field sites. Previous research used one of the field sites to test the original modification of the RZWQM to simulate directly connected macropores for bromide and the parent compound, but not for the metabolite. This research will evaluate RZWQM for parent/metabolite transformation and transport at this first field site, along with evaluating the model at an additional field site to evaluate whether the parameters for direct connectivity are transferable and whether model performance is consistent for the two field sites with unique soil, hydrologic, and environmental conditions. Isoxaflutole, the active ingredient in BALANCE^® herbicide, was applied to both fields. Isoxaflutole rapidly degrades into a metabolite (RPA 202248). This research used calibrated RZWQM models for each field based on observed subsurface drain flow and/or edge of field conservative tracer concentrations in subsurface flow. The calibrated models for both field sites required a portion (approximately 2% but this fraction may require calibration) of the available water and chemical in macropore flow to be routed directly into the subsurface drains to simulate peak concentrations in edge of field subsurface drain flow shortly after chemical applications. Confirming the results from the first field site, the existing modification for directly connected macropores continually failed to predict pesticide concentrations on the recession limbs of drainage hydrographs, suggesting that the current strategy only partially accounts for direct connectivity. Thirty‐year distributions of annual mass (drainage) loss of parent and metabolite in terms of percent of isoxaflutole applied suggested annual simulated percent losses of parent and metabolite (3.04 and 1.31%) no greater in drainage than losses in runoff on nondrained fields as reported in the literature.     

Spatio-Temporal Variability of Soil Respiration of Forest Ecosystems in China: Influencing Factors and Evaluation Model

Understanding the influencing factors of the spatio-temporal variability of soil respiration (R _s) across different ecosystems as well as the evaluation model of R _s is critical to the accurate prediction of future changes in carbon exchange between ecosystems and the atmosphere. R _s data from 50 different forest ecosystems in China were summarized and the influences of environmental variables on the spatio-temporal variability of R _s were analyzed. The results showed that both the mean annual air temperature and precipitation were weakly correlated with annual R _s, but strongly with soil carbon turnover rate. R _s at a reference temperature of 0°C was only significantly and positively correlated with soil organic carbon (SOC) density at a depth of 20 cm. We tested a global-scale R _s model which predicted monthly mean R _s (R _s,monthly) from air temperature and precipitation. Both the original model and the reparameterized model poorly explained the monthly variability of R _s and failed to capture the inter-site variability of R _s. However, the residual of R _s,monthly was strongly correlated with SOC density. Thus, a modified empirical model (TPS model) was proposed, which included SOC density as an additional predictor of R _s. The TPS model explained monthly and inter-site variability of R _s for 56% and 25%, respectively. Moreover, the simulated annual R _s of TPS model was significantly correlated with the measured value. The TPS model driven by three variables easy to be obtained provides a new tool for R _s prediction, although a site-specific calibration is needed for using at a different region.     

CALIBRATION OF A WATER-BALANCE MODEL FOR SMALL WATERSHEDS IN EASTERN OREGON1

ABSTRACT: The BURP water-balance model was calibrated for 13 small (0.46 to 7.00 mi²), forested watersheds in the Blue Mountains of eastern Oregon where snowmelt is the dominant source of runoff. BURP is the model name and is not an acronym. Six of the 16 parameters in BURP were calibrated. The subsurface recession coefficient and three subsurface water-storage parameters were most sensitive for simulating monthly flow. Calibrated subsurface recession coefficients ranged from 0.988 to 0.998. The subsurface-water storage parameters were calibrated at between 20 to 120 percent of their initial values obtained from a category III soil survey. That reconnaissance-level survey was apparently too broad to accurately reflect subsurface-water storage in small watersheds. Tests of model performance showed BURP is capable of producing accurate simulations of monthly flow for mountainous, snow-dominated watersheds with shallow (< 4 ft) soils when calibrated with 2 to 4 years of streamflow data. A regression of observed versus simulated monthly flows with data from all watersheds combined showed that BURP accounted for 85 percent of the variability in observed flows, which ranged from 0.01 to 20.8 inches, but underpredicted high flow months, with a slope of 1.15 that is significantly different from 1.0 (p = 0.05). Without prior calibration, subsurface-water storage parameters appeared to be the greatest source of potential error.     

Water quality modeling of fertilizer management impacts on nitrate losses in tile drains at the field scale

Nitrate losses from subsurface tile drained row cropland in the Upper Midwest U.S. contribute to hypoxia in the Gulf of Mexico. Strategies are needed to reduce nitrate losses to the Mississippi River. This paper evaluates the effect of fertilizer rate and timing on nitrate losses in two (East and West) commercial row crop fields located in south-central Minnesota. The Agricultural Drainage and Pesticide Transport (ADAPT) model was calibrated and validated for monthly subsurface tile drain flow and nitrate losses for a period of 1999-2003. Good agreement was found between observed and predicted tile drain flow and nitrate losses during the calibration period, with Nash-Sutcliffe modeling efficiencies of 0.75 and 0.56, respectively. Better agreements were observed for the validation period. The calibrated model was then used to evaluate the effects of rate and timing of fertilizer application on nitrate losses with a 50-yr climatic record (1954-2003). Significant reductions in nitrate losses were predicted by reducing fertilizer application rates and changing timing. A 13% reduction in nitrate losses was predicted when fall fertilizer application rate was reduced from 180 to 123 kg/ha. A further 9% reduction in nitrate losses can be achieved when switching from fall to spring application. Larger reductions in nitrate losses would require changes in fertilizer rate and timing, as well as other practices such as changing tile drain spacings and/or depths, fall cover cropping, or conversion of crop land to pasture.     

SIMULATING TILE DRAINAGE AND NITRATE LEACHING UNDER A POTATO CROP1

ABSTRACT: A 2.2-hectare potato (Solanum tuberosum L. cv Chieftain) field at Saint Leonard d'Aston, Quebec (lat. 72° 24′ 30″ long. 46° 5′ 30″) was instrumented to measure tile drain flow over two growing seasons, 1989 and 1990. The soil was a Sainte Jude sandy loam. Soil properties and nitrate concentrations in the drain flow were measured. The CREAMS (Chemicals, Runoff and rosion from agricultural Management systems) computer simulation model was validated for the study site. CREAMS underpredicted event percolation depths. However, total monthly percolation depths were close to observed values. CREAMS overpredicted event nitrate concentrations leached to tile drainage. There was a poor match between predicted and observed event nitrate concentrations in drain flow (coefficient of predictability, CP_A= 104.95). Based on a sensitivity analysis, input parameters, representative of local conditions, were determined for the CREAMS hydrology and nutrient submodels.     

Integrating Seasonal Information on Nutrients and Benthic Algal Biomass into Stream Water Quality Monitoring

Benthic chlorophyll a (BChl a) and environmental factors that influence algal biomass were measured monthly from February through October in 22 streams from three agricultural regions of the United States. At‐site maximum BChl a ranged from 14 to 406 mg/m² and generally varied with dissolved inorganic nitrogen (DIN): 8 out of 9 sites with at‐site median DIN >0.5 mg/L had maximum BChl a >100 mg/m². BChl a accrued and persisted at levels within 50% of at‐site maximum for only one to three months. No dominant seasonal pattern for algal biomass accrual was observed in any region. A linear model with DIN, water surface gradient, and velocity accounted for most of the cross‐site variation in maximum chlorophyll a (adjusted R² = 0.7), but was no better than a single value of DIN = 0.5 mg/L for distinguishing between low and high‐biomass sites. Studies of nutrient enrichment require multiple samples to estimate algal biomass with sufficient precision given the magnitude of temporal variability of algal biomass. An effective strategy for regional stream assessment of nutrient enrichment could be based on a relation between maximum BChl a and DIN based on repeat sampling at sites selected to represent a gradient in nutrients and application of the relation to a larger number of sites with synoptic nutrient information.     

Tillage and nutrient source effects on surface and subsurface water quality at corn planting

This study quantified the effects of tillage (moldboard plowing [MP], ridge tillage [RT]) and nutrient source (manure and commercial fertilizer [urea and triple superphosphate]) on sediment, NH4+ -N, NO3- -N, total P, particulate P, and soluble P losses in surface runoff and subsurface tile drainage from a clay loam soil. Treatment effects were evaluated using simulated rainfall immediately after corn (Zea mays L.) planting, the most vulnerable period for soil erosion and water quality degradation. Sediment, total P, soluble P, and NH4+ -N losses mainly occurred in surface runoff. The NO3- -N losses primarily occurred in subsurface tile drainage. In combined (surface and subsurface) flow, the MP treatment resulted in nearly two times greater sediment loss than RT (P < 0.01). Ridge tillage with urea lost at least 11 times more NH4+ -N than any other treatment (P < 0.01). Ridge tillage with manure also had the most total and soluble P losses of all treatments (P < 0.01). If all water quality parameters were equally important, then moldboard plow with manure would result in least water quality degradation of the combined flow followed by moldboard plow with urea or ridge tillage with urea (equivalent losses) and ridge tillage with manure. Tillage systems that do not incorporate surface residue and amendments appear to be more vulnerable to soluble nutrient losses mainly in surface runoff but also in subsurface drainage (due to macropore flow). Tillage systems that thoroughly mix residue and amendments in surface soil appear to be more prone to sediment and sediment-associated nutrient (particulate P) losses via surface runoff.     

CLIMATE CHHANGE SENSITIVITY ASSESSMENT ON UPPER MISSISSIPPI RIVER BASIN STREAMFLOWS USING SWAT1

ABSTRACT: The Soil and Water Assessment Tool (SWAT) model was used to assess the effects of potential future climate change on the hydrology of the Upper Mississippi River Basin (UMRB). Calibration and validation of SWAT were performed using monthly stream flows for 1968–1987 and 1988–1997, respectively. The R² and Nash‐Sutcliffe simulation efficiency values computed for the monthly comparisons were 0.74 and 0.69 for the calibration period and 0.82 and 0.81 for the validation period. The effects of nine 30‐year (1968 to 1997) sensitivity runs and six climate change scenarios were then analyzed, relative to a scenario baseline. A doubling of atmospheric CO₂ to 660 ppmv (while holding other climate variables constant) resulted in a 36 percent increase in average annual streamflow while average annual flow changes of ?49, ?26, 28, and 58 percent were predicted for precipitation change scenarios of ?20, ?10, 10, and 20 percent, respectively. Mean annual streamflow changes of 51,10, 2, ?6, 38, and 27 percent were predicted by SWAT in response to climate change projections generated from the CISRO‐RegCM2, CCC, CCSR, CISRO‐Mk2, GFDL, and HadCMS general circulation model scenarios. High seasonal variability was also predicted within individual climate change scenarios and large variability was indicated between scenarios within specific months. Overall, the climate change scenarios reveal a large degree of uncertainty in current climate change forecasts for the region. The results also indicate that the simulated UMRB hydrology is very sensitive to current forecasted future climate changes.     

MODELING STREAMFLOW FROM ARTIFICIALLY DRAINED AGRICULTURAL WATERSHEDS IN ILLINOIS1

ABSTRACT: Evaluation of the applicability and validity of hydrologic simulation models for various cropping systems in different hydrogeologic and soil conditions is needed for a range of spatial scales. We calibrated and tested the ADAPT model for simulating streamflow from 552 to 1,985 km² watersheds in central Illinois, where more than 79 percent of the land is used for maize‐soybean production and tile drainage is common. Model calibration was performed with a seven year period (1987 through1993) of measured streamflow from one of the watersheds, and model testing was done using independent weather and measured streamflow data from the two neighboring watersheds for the same seven year period. Simulations of annual streamflow were accurate with a coefficient of determination and Willmott's index of agreement of 0.98 and 0.99, respectively. For simulation of monthly streamflow, Willmott's index of agreement ranged from 0.93 to 0.95. For simulation of daily streamflow, Willmott's index of agreement ranged from 0.84 to 0.85. The daily simulations challenged the temporal and spatial resolution of our measured precipitation data. Discrepancies between simulated and measured data may result from the model's inability to effectively address frozen soils and snowmelt runoff processes and in accurately representing evapotranspiration.     

Model Calculations of Total Maximum Daily Loads of Mercury for Drainage Lakes1

Abstract: The Watershed Analysis Risk Management Framework watershed model was enhanced to simulate the transport and fate of mercury and to calculate the fish mercury concentrations (FMC) attained by fish through the food web. The model was applied to Western Lake Superior Basin of Minnesota, which has many peat lands and lakes. Topographic, land use, and soil data were used to set up the model. Meteorology and precipitation chemistry data from nearby monitoring stations were compiled to drive the model. Simulated flow and mercury concentrations for several stream stations were comparable to available data. The model was used to perform mercury total maximum daily load calculations for two contrasting drainage lakes (Wild Rice Lake and Whiteface Reservoir). The model results for wet deposition, dry deposition, evasion, watershed yield, and soil sequestration of mercury were comparable with available actual data. The model predicted lake ice cover from November to April and weak stratification in summer, typical of shallow lakes in cold regions. The simulated sulfate decrease and methylmercury increase near the lake bottom in late summer are caused by sulfate reduction and mercury methylation that occur in the surficial sediment. Simulated FMC were within the range of observed values and the R² of correlation between the simulated and observed FMC was 0.77. Under the 1989‐2004 base condition, the average simulated FMC of four‐year‐old walleye was 0.31 μg/g for Whiteface Reservoir and 0.15 μg/g for Wild Rice Lake. The FMC criterion in Minnesota is 0.2 μg/g. Wild Rice Lake already meets this criterion without any load reduction. The model showed that a 65% reduction in atmospheric mercury deposition will not, by itself, allow Whiteface Reservoir to meet the criterion in 15 years. Additional best management practices will be needed to reduce 50% of the watershed input.