Constraining SWAT Calibration with Remotely Sensed Evapotranspiration Data

Historically, many watershed studies have been based on using the streamflow flux, typically from a single gauge at the basin's outlet, to support calibration. In this setting, there is great potential for equifinality of parameters during the optimization process, especially for parameters that are not directly related to streamflow. Therefore, some of the optimal parameter values achieved during the autocalibration process may be physically unrealistic. In recent decades a vast array of data from land surface models and remote sensing platforms can help to constrain hydrologic fluxes such as evapotranspiration (ET). While the spatial resolution of these ancillary datasets varies, the continuous spatial coverage of these gridded datasets provides flux measurements across the entire basin, in stark contrast to point‐based streamflow data. This study uses Global Land Evaporation: the Amsterdam Model data to constrain Soil and Water Assessment Tool parameter values associated with ET to a more physically realistic range. The study area is the Little Washita River Experimental Watershed, in southern Oklahoma. Traditional objective metrics such as the Nash‐Sutcliffe coefficients record no performance improvement after application of this method. However, there is a dramatic increase in the number of days with receding flow where simulations match observed streamflow.     

Does the Temporal Resolution of Precipitation Input Influence the Simulated Hydrological Components Employing the SWAT Model?

This study aimed to evaluate the influence of sub‐daily precipitation time steps on model performance and hydrological components by applying the Green and Ampt infiltration method using the Soil and Water Assessment Tool (SWAT). Precipitation was measured at a resolution of 0.1 mm and aggregated to 5‐, 15‐, 30‐, and 60‐min time steps. Daily discharge data over a 10‐year period were used to calibrate and validate the model. Following a global sensitivity analysis, relevant parameters were optimized through an automatic calibration procedure using SWAT‐CUP for each time step. Daily performance statistics were almost equal among all four time steps (NSE ≈ 0.47). Discharge mainly consisted of groundwater flow (55%) and tile flow (42%), in reasonable proportions for the investigated catchment. In conclusion, model outputs were almost identical, showing simulations responded nearly independently of the chosen precipitation time step. This held true for (1) the selection of sensitive parameters, (2) performance statistics, (3) the shape of the hydrographs, and (4) flow components. However, a scenario analysis revealed that the precipitation time step becomes important when saturated hydraulic conductivities are low and curve numbers are high. The study suggests that there is no need in using precipitation time steps <1 h for lowland catchments dominated by soils with a low surface runoff potential if daily flow values are being considered. 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.     

Watershed Regressions for Pesticides (warp) Models for Predicting Atrazine Concentrations in Corn Belt Streams1

Stone, Wesley W. and Robert J. Gilliom, 2012. Watershed Regressions for Pesticides (WARP) Models for Predicting Atrazine Concentrations in Corn Belt Streams. Journal of the American Water Resources Association (JAWRA) 48(5): 970‐986. DOI: 10.1111/j.1752‐1688.2012.00661.x Abstract: Watershed Regressions for Pesticides (WARP) models, previously developed for atrazine at the national scale, are improved for application to the United States (U.S.) Corn Belt region by developing region‐specific models that include watershed characteristics that are influential in predicting atrazine concentration statistics within the Corn Belt. WARP models for the Corn Belt (WARP‐CB) were developed for annual maximum moving‐average (14‐, 21‐, 30‐, 60‐, and 90‐day durations) and annual 95th‐percentile atrazine concentrations in streams of the Corn Belt region. The WARP‐CB models accounted for 53 to 62% of the variability in the various concentration statistics among the model‐development sites. Model predictions were within a factor of 5 of the observed concentration statistic for over 90% of the model‐development sites. The WARP‐CB residuals and uncertainty are lower than those of the National WARP model for the same sites. Although atrazine‐use intensity is the most important explanatory variable in the National WARP models, it is not a significant variable in the WARP‐CB models. The WARP‐CB models provide improved predictions for Corn Belt streams draining watersheds with atrazine‐use intensities of 17 kg/km² of watershed area or greater.     

Comparison of Stream Invertebrate Response Models for Bioassessment Metrics1

Waite, Ian R., Jonathan G. Kennen, Jason T. May, Larry R. Brown, Thomas F. Cuffney, Kimberly A. Jones, and James L. Orlando, 2012. Comparison of Stream Invertebrate Response Models for Bioassessment Metrics. Journal of the American Water Resources Association (JAWRA) 48(3): 570-583. DOI: 10.1111/j.1752-1688.2011.00632.x Abstract: We aggregated invertebrate data from various sources to assemble data for modeling in two ecoregions in Oregon and one in California. Our goal was to compare the performance of models developed using multiple linear regression (MLR) techniques with models developed using three relatively new techniques: classification and regression trees (CART), random forest (RF), and boosted regression trees (BRT). We used tolerance of taxa based on richness (RICHTOL) and ratio of observed to expected taxa (O/E) as response variables and land use/land cover as explanatory variables. Responses were generally linear; therefore, there was little improvement to the MLR models when compared to models using CART and RF. In general, the four modeling techniques (MLR, CART, RF, and BRT) consistently selected the same primary explanatory variables for each region. However, results from the BRT models showed significant improvement over the MLR models for each region; increases in R² from 0.09 to 0.20. The O/E metric that was derived from models specifically calibrated for Oregon consistently had lower R² values than RICHTOL for the two regions tested. Modeled O/E R² values were between 0.06 and 0.10 lower for each of the four modeling methods applied in the Willamette Valley and were between 0.19 and 0.36 points lower for the Blue Mountains. As a result, BRT models may indeed represent a good alternative to MLR for modeling species distribution relative to environmental variables.     

Modeling urban films using a dynamic multimedia fugacity model

A thin film coats impervious urban surfaces that can act as a source or sink of organic pollutants to the greater environment. We review recent developments in the understanding of film and film-associated pollutant behavior and incorporate them into an unsteady-state version of the fugacity based Multimedia Urban Model (MUM), focusing on detailed considerations of surface film dynamics. The model is used to explore the conditions under which these atmospherically-derived films act as a temporary source of chemicals to the air and/or storm water. Assuming film growth of 2.1 nm d⁻¹ (Wu et al., 2008a), PCB congeners 28 and 180 reach air-film equilibrium within hours and days, respectively. The model results suggest that the film acts as a temporary sink of chemicals from air during dry and cool weather, as a source to air in warmer weather, and as a source to storm water and soil during rain events. Using the downtown area of the City of Toronto Canada, as a case study, the model estimates that nearly 1 g d⁻¹ of ∑₅PCBs are transferred from air to film to storm water.     

Critical Steps for Continuing Advancement of Satellite Rainfall Applications for Surface Hydrology in the Nile River Basin1

Gebremichael, Mekonnen, Emmanouil N. Anagnostou, and Menberu M. Bitew, 2010. Critical Steps for Continuing Advancement of Satellite Rainfall Applications for Surface Hydrology in the Nile River Basin. Journal of the American Water Resources Association (JAWRA) 46(2):361-366. DOI: 10.1111/j.1752-1688.2010.00428.x. Abstract: Given the increasingly higher resolution and data accessibility, satellite precipitation products could be useful for hydrological application in the Nile River Basin, which is characterized by lack of reasonably dense hydrological in situ sensors and lack of access to the existing dataset. However, in the absence of both extreme caution and research results for the Nile basin, the satellite rainfall (SR) products may not be used, or may even be used erroneously. We identify two steps that are critical to enhance the value of SR products for hydrological applications in the Nile basin. The first step is to establish representative validation sites in the Nile basin. The validation site will help to quantify the errors in the different kinds of SR products, which will be used to select the best products for the Nile basin, include the errors in decision making, and design strategies to minimize the errors. Using rainfall measurements collected from the unprecedented high-density rain gauge network over a small region within the Nile basin, we indicate that SR estimates could be subject to significant errors, and quantification of estimation errors by way of establishing validation sites is critically important in order to use the SR products. The second step is to identify the degree of hydrologic model complexity required to obtain more accurate hydrologic simulation results for the Nile basin when using SR products as input. The level of model complexity may depend on basin size and SR algorithm, and further research is needed to spell out this dependence for the Nile basin.     

Multimedia fate modeling and comparative impact on freshwater ecosystems of pharmaceuticals from biosolids-amended soils

This study modeled the impact on freshwater ecosystems of pharmaceuticals detected in biosolids following application on agricultural soils. The detected sulfonamides and hydrochlorothiazide displayed comparatively moderate retention in solid matrices and, therefore, higher transfer fractions from biosolids to the freshwater compartment. However, the residence times of these pharmaceuticals in freshwater were estimated to be short due to abiotic degradation processes. The non-steroidal anti-inflammatory mefenamic acid had the highest environmental impact on aquatic ecosystems and warrants further investigation. The estimation of the solid-water partitioning coefficient was generally the most influential parameter of the probabilistic comparative impact assessment. These results and the modeling approach used in this study serve to prioritize pharmaceuticals in the research effort to assess the risks and the environmental impacts on aquatic biota of these emerging pollutants.     

宁夏地区光伏发电环境效益分析及建模

在阐述了宁夏地区光伏发电发展的优势及发展现状的基础上,分析了光伏发电开发利用对宁夏生态环境产生的积极和消极的影响,并为降低消极影响提出了相应的对策;为了定量直观的评估光伏发电积极影响和消极影响所产生的环境效益,分别建立了节能减排、防风固沙和土地增值、景观效应、植被恢复费用、噪声污染超标罚款、光污染费用等价值模型及光伏发电总的生态环境效益价值模型;最后以中电投宁夏能源铝中卫香山光伏电站为例进行计算验证,结果表明模型的合理性和准确性,同时光伏电站的节能减排效益显著。相对于传统的火力发电,表现了优越的环境效益。     

美国地震风险评估中灾害模型的探讨

介绍了以美国地震灾害模型为例的地震风险评估在金融保险业和土木工程两大领域的应用,回顾了地震灾害风险评估发展的历史背景,其中,金融保险业的模型更测重于对历史理赔数据的分析,土木工程业的模型则注重于对工程结构的分析。着重探讨了评估模型中地震灾害部分的建立,其中包括致灾因子模块及孕灾环境模块。针对建模中致灾因子、孕灾环境等的关键要素,如地震源及震级、地震衰减方式、场地土壤土质条件、地震复发率等,作了分析及归纳总结,并系统地介绍了所涉及要素的数据信息获取过程及对应的建模方法,进一步探讨了相应的风险量化的步骤和方法。     

Inhalation dose assessment of indoor radon progeny using biokinetic and dosimetric modeling and its application to Jordanian population

High indoor radon concentrations in Jordan result in internal exposures of the residents due to the inhalation of radon and its short-lived progeny. It is therefore important to quantify the annual effective dose and further the radiation risk to the radon exposure. This study describes the methodology and the biokinetic and dosimetric models used for calculation of the inhalation doses exposed to radon progeny. The regional depositions of aerosol particles in the human respiratory tract were firstly calculated. For the attached progeny, the activity median aerodynamic diameters of 50 nm, 230 nm and 2500 nm were chosen to represent the nucleation, accumulation and coarse modes of the aerosol particles, respectively. For the unattached progeny, the activity median thermodynamic diameter of 1 nm was chosen to represent the free progeny nuclide in the room air. The biokinetic models developed by the International Commission on Radiological Protection (ICRP) were used to calculate the nuclear transformations of radon progeny in the human body, and then the dosimetric model was applied to estimate the organ equivalent doses and the effective doses with the specific effective energies derived from the mathematical anthropomorphic phantoms. The dose conversion coefficient estimated in this study was 15 mSv WLM⁻¹ which was in the range of the values of 6-20 mSv WLM⁻¹ reported by other investigators. Implementing the average indoor radon concentration in Jordan, the annual effective doses were calculated to be 4.1 mSv y⁻¹ and 0.08 mSv y⁻¹ due to the inhalation of radon progeny and radon gas, respectively. The total annual effective dose estimated for Jordanian population was 4.2 mSv y⁻¹. This high annual effective dose calculated by the dosimetric approach using ICRP biokinetic and dosimetric models resulted in an increase of a factor of two in comparison to the value by epidemiological study. This phenomenon was presented by the ICRP in its new published statement on radon.