Precipitation changes impact stream discharge, nitrate-nitrogen load more than agricultural management changes

Nitrate-N losses to surface waters in the Upper Midwest of the Untied States have increased in recent decades, contributing to hypoxia in the Gulf of Mexico. This paper investigates whether increasing nitrate-N export from cropland in the Upper Midwest since the late 1960s results from changes in land use or climate. The Agricultural Drainage and Pesticide Transport (ADAPT) Model simulated current and historical agricultural systems under past and recent wet climate for Seven Mile Creek in Minnesota. Simulations were run with management and climate for three distinctly different periods--namely, 1965 to 1969, 1976 to 1980, and 1999 to 2003 (wettest period). Results showed discharge and nitrate-N losses responded more to changes in climate than management. The wetter period (1999-2003) caused a simulated 70% increase in discharge under 1960s-era management compared with that period's observed climate and a simulated 51% increase in discharge under 1970s-era management compared with the 1976 to 1980 climate. The recent, wetter climate also produced a 62% increase in nitrate-N losses for 1960s-era management compared with the actual climate and a 137% increase in nitrate-N losses for 1978 management conditions compared with actual 1970s climate. Had recent climate been in place and stable since 1965, agricultural changes would have decreased discharge by 6.4% through the late 1970s and then by another 21.1% under modern management but would have increased nitrate-N losses by 184% through the late 1970s and then decreased nitrate-N losses by 13.5% between 1978 and 2001. Management changes that were important drivers included increasing N-fertilizer rates, increases in corn acreage, and increases in crop yield. But the most important factor driving increased nitrate-N losses from agriculture since the 1970s was an increasingly wetter climate.     

Operational Evapotranspiration Mapping Using Remote Sensing and Weather Datasets: A New Parameterization for the SSEB Approach

The increasing availability of multi‐scale remotely sensed data and global weather datasets is allowing the estimation of evapotranspiration (ET) at multiple scales. We present a simple but robust method that uses remotely sensed thermal data and model‐assimilated weather fields to produce ET for the contiguous United States (CONUS) at monthly and seasonal time scales. The method is based on the Simplified Surface Energy Balance (SSEB) model, which is now parameterized for operational applications, renamed as SSEBop. The innovative aspect of the SSEBop is that it uses predefined boundary conditions that are unique to each pixel for the “hot” and “cold” reference conditions. The SSEBop model was used for computing ET for 12 years (2000‐2011) using the MODIS and Global Data Assimilation System (GDAS) data streams. SSEBop ET results compared reasonably well with monthly eddy covariance ET data explaining 64% of the observed variability across diverse ecosystems in the CONUS during 2005. Twelve annual ET anomalies (2000‐2011) depicted the spatial extent and severity of the commonly known drought years in the CONUS. More research is required to improve the representation of the predefined boundary conditions in complex terrain at small spatial scales. SSEBop model was found to be a promising approach to conduct water use studies in the CONUS, with a similar opportunity in other parts of the world. The approach can also be applied with other thermal sensors such as Landsat.     

Estimating Evapotranspiration for Dryland Cropping Systems in the Semiarid Texas High Plains Using SWAT

The Soil and Water Assessment Tool (SWAT) is one of the most widely used watershed models for simulating hydrology in response to agricultural management practices. However, limited studies have been performed to evaluate the SWAT model's ability to estimate daily and monthly evapotranspiration (ET) in semiarid regions. ET values were simulated using ArcSWAT 2012 for a lysimeter field managed under dryland conditions at the USDA‐ARS Conservation and Production Research Laboratory at Bushland, Texas, and compared with measured lysimeter values from 2000 to 2010. Two scenarios were performed to compare SWAT's performance: (1) use of default plant leaf area index (LAI) values in the embedded plant database and (2) adjusted LAI values. Scenario 1 resulted in an “unsatisfactory” Nash‐Sutcliffe efficiency (NSE) of 0.42 and 0.38 for the calibration and validation periods, respectively. Scenario 2 resulted in a “satisfactory” NSE value for the calibration period while achieving a “good” NSE of 0.70 for the validation period. SWAT generally underestimated ET at both the daily and monthly levels. Overestimation during fallow years may be due to the limitations of the pothole function used to simulate furrow diking. Users should be aware of potential errors associated with using default LAI parameters. Inaccuracies in ET estimation may also stem from errors in the plant stress functions, particularly when evaluating water management practices for dryland watersheds.     

A Parallel Computation Tool to Enable Dynamic Sensitivity and Model Performance Analysis of APEX: Evapotranspiration Modeling

Global sensitivity analysis can be used for assessing the relative importance of model parameters on model outputs. The sensitivity of parameters usually indicates a temporal variation due to variation in the environmental conditions (e.g., variation in weather or plant growth). In addition, the size of averaging window by which the outputs of a model are aggregated or averaged may impact parameter sensitivities. In this study, temporal variation of parameters sensitives, model performance, as well as the impact of the size of time‐averaging window on evapotranspiration (ET) prediction using the Agricultural Policy/Environmental eXtender (APEX) model are investigated. To achieve these goals, an open‐source package named PARAPEX was developed in R and used to perform dynamic sensitivity and model performance analysis of APEX using parallel computation. PARAPEX reduced the computation time from 5,939 to 379 s (using 20 and 1 computation nodes, respectively). The sensitivity analysis results indicated the parameters accounting for the reducing effect of plant cover on evaporation from the soil surface, the effect of soil on the plant root growth, and the effect of cycling and transformation dynamics of organic matter at the top soil layer as the top sensitive parameters based on the mean daily simulated ET and the Nash–Sutcliffe model performance measure. The dynamic performance analysis indicated poor ET predictions by APEX during the growing seasons. Editor's note : This paper is part of the featured series on Optimizing Ogallala Aquifer Water Use to Sustain Food Systems. See the February 2019 issue for the introduction and background to the series.     

Transition Pathways to Sustainable Agricultural Water Management: A Review of Integrated Modeling Approaches

Agricultural water management (AWM) is an interdisciplinary concern, cutting across traditional domains such as agronomy, climatology, geology, economics, and sociology. Each of these disciplines has developed numerous process‐based and empirical models for AWM. However, models that simulate all major hydrologic, water quality, and crop growth processes in agricultural systems are still lacking. As computers become more powerful, more researchers are choosing to integrate existing models to account for these major processes rather than building new cross‐disciplinary models. Model integration carries the hope that, as in a real system, the sum of the model will be greater than the parts. However, models based upon simplified and unrealistic assumptions of physical or empirical processes can generate misleading results which are not useful for informing policy. In this article, we use literature and case studies from the High Plains Aquifer and Southeastern United States regions to elucidate the challenges and opportunities associated with integrated modeling for AWM and recommend conditions in which to use integrated models. Additionally, we examine the potential contributions of integrated modeling to AWM — the actual practice of conserving water while maximizing productivity. Editor's note : This paper is part of the featured series on Optimizing Ogallala Aquifer Water Use to Sustain Food Systems. See the February 2019 issue for the introduction and background to the series.     

Evaluation of Satellite‐Derived Rainfall Data for Multiple Physio‐Climatic Regions in the Santiago River Basin,Mexico

Assessment of water resources requires reliable rainfall data, and rain gauge networks may not provide adequate spatial representation due to limited point measurements. The Tropical Rainfall Measuring Mission (TRMM) provides rainfall data at global scale, and has been used with good results. However, TRMM data are an indirect measurement of rainfall, and therefore must be validated for its proper use. In this work, a validation scheme was designed and implemented to compare the TRMM Version 7 (V7) monthly rainfall product at different time frames with data measured in two hydrologic subregions of the Santiago River Basin (SRB) in Mexico: Río Alto Santiago and Río Bajo Santiago (RBS). Additionally, three physio‐climatic regions provide an assessment of the interplay of topography, distance from coastal regions, and seasonal weather patterns on the correspondence between both datasets. The TRMM V7 rainfall product exhibited good agreement with the rain gauge data particularly for the RBS and for the whole SRB during wettest summer and autumn seasons. However, strong regional dependence was observed due to differences in climate and topography. Overall, in spite of some noted underestimations, the monthly TRMM V7 rainfall product was found to provide useful information that can be used to complement limited monitoring as is the case of RBS. An improved combined rainfall product could be generated and thus gaining the most benefits from both data sources.     

Ozonation of the 5-fluorouracil anticancer drug and its prodrug capecitabine: Reaction kinetics,oxidation mechanisms,and residual toxicity

Specific second-order rate constants were determined for 5-FU and CAP with ozone. Reaction sites were confirmed by kinetics, Fukui analysis, and products. The olefin moiety was the main ozone reaction site for 5-FU and CAP. Carboxylic acids comprised most of the residual TOC for 5-FU. Ozonation removed the toxicity associated with 5-FU and products but not CAP. Anticancer drugs (ADs) have been detected in the environment and represent a risk to aquatic organisms, necessitating AD removal in drinking water and wastewater treatment. In this study, ozonation of the most commonly used antimetabolite ADs, namely 5-fluorouracil (5-FU) and its prodrug capecitabine (CAP), was investigated to determine reaction kinetics, oxidation mechanisms, and residual toxicity. The specific second-order rate constants between aqueous ozone and 5-FU, 5-FU⁻, 5-FU²⁻, CAP, and CAP⁻ were determined to be 7.07(±0.11)×10⁴ M⁻¹·s⁻¹, 1.36(±0.06)×10⁶ M⁻¹·s⁻¹, 2.62(±0.17)×10⁷ M⁻¹·s⁻¹, 9.69(±0.08)×10³ M⁻¹·s⁻¹, and 4.28(±0.07)×10⁵ M⁻¹·s⁻¹, respectively; furthermore, the second-order rate constants for ^•OH reaction with 5-FU and CAP at pH 7 were determined to be 1.85(±0.20)×10⁹ M⁻¹·s⁻¹ and 9.95(±0.26)×10⁹ M⁻¹·s⁻¹, respectively. Density functional theory was used to predict the main ozone reaction sites of 5-FU (olefin) and CAP (olefin and deprotonated secondary amine), and these mechanisms were supported by the identified transformation products. Carboxylic acids constituted a majority of the residual organic matter for 5-FU ozonation; however, carboxylic acids and aldehydes were important components of the residual organic matter generated by CAP. Ozone removed the toxicity of 5-FU to Vibrio fischeri, but the residual toxicity of ozonated CAP solutions exhibited an initial increase before subsequent removal. Ultimately, these results suggest that ozone is a suitable technology for treatment of 5-FU and CAP, although the residual toxicity of transformation products must be carefully considered.     

Impact of Length of Dataset on Streamflow Calibration Parameters and Performance of APEX Model

Due to resource constraints, long‐term monitoring data for calibration and validation of hydrologic and water quality models are rare. As a result, most models are calibrated and, if possible, validated using limited measured data. However, little research has been done to determine the impact of length of available calibration data on model parameterization and performance. The main objective of this study was to evaluate the impact of length of calibration data (LCD) on parameterization and performance of the Agricultural Policy Environmental eXtender model for predicting daily, monthly, and annual streamflow. Long‐term (1984‐2015) measured daily streamflow data from Rock Creek watershed, an agricultural watershed in northern Ohio, were used for this study. Data were divided into five Short (5‐year), two Medium (15‐year), and one Long (25‐year) streamflow calibration data scenarios. All LCD scenarios were calibrated and validated at three time steps: daily, monthly, and annual. Results showed LCD affected the ability of the model to accurately capture temporal variability in simulated streamflow. However, overall average streamflow, water budgets, and crop yields were simulated reasonably well for all LCD scenarios.     

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.