Assessing the Impact of Areal Precipitation Input on Streamflow Simulations Using the SWAT Model1

Masih Ilyas, Shreedhar Maskey, Stefan Uhlenbrook, and Vladimir Smakhtin, 2011. Assessing the Impact of Areal Precipitation Input on Streamflow Simulations Using the SWAT Model. Journal of the American Water Resources Association (JAWRA) 47(1):179‐195. DOI: 10.1111/j.1752‐1688.2010.00502.x Abstract: Reduction of input uncertainty is a challenge in hydrological modeling. The widely used model Soil Water Assessment Tool (SWAT) uses the data of a precipitation gauge nearest to the centroid of each subcatchment as an input for that subcatchment. This may not represent overall catchment precipitation conditions well. This paper suggests an alternative – using areal precipitation obtained through interpolation. The effectiveness of this alternative is evaluated by comparing its simulations with those based on the standard SWAT precipitation input procedure. The model is applied to mountainous semiarid catchments in the Karkheh River basin, Iran. The model performance is evaluated at daily, monthly, and annual scales by using a number of performance indicators at 15 streamflow gauging stations each draining an area in the range of 590‐42,620 km². The comparison suggests that the use of areal precipitation improves model performance particularly in small subcatchments in the range of 600‐1,600 km². The modified areal precipitation input results in increased reliability of simulated streamflows in the areas of low rain gauge density. Both precipitation input methods result in reasonably good simulations for larger catchments (over 5,000 km²). The use of areal precipitation input improves the accuracy of simulated streamflows with spatial resolution and density of rain gauges having significant impact on results.     

APPLICATION OF KRIGING TO ESTIMATING MEAN ANNUAL PRECIPITATION IN A REGION OF OROGRAPHIC INFLUENCE1

ABSTRACT: Estimates of mean annual precipitation (MAP) over areas are the starting point for all computations of water and chemical balances for drainage basins and surface water bodies. Any errors in the estimates of MAP are propagated through the balance computations. These errors can be due to: (1) failures of individual gages to collect the amount of precpitation that actually falls; (2) operator errors; and (3) failure of the raingage network to adequately sample the region of interest. This paper attempts to evaluate the last of these types of error by applying kriging in two different approaches to estimating MAP in New Hampshire and Vermont, USA. The data base is the 1951–1980 normal precipitation at 120 raingages in the two states and in adjacent portions of bordering states and provinces. In the first approach, kriging is applied directly to the MAP values, while in the second, kriging is applied to a “precipitation delivery factor” that represents the MAP with the orographic effect removed. The first approach gives slightly better kriged estimates of MAP at seven validation stations that were not included in the original analysis, but results in an error surface that is highly contorted and in larger maximum errors over most of the region. The second approach had a considerably smoother error surface and, thus, is generally preferable as a basis for point and areal estimates of MAP. MAP estimates in the region have 95 percent confidence intervals of about 20 cm/yr at low and moderate elevations, and up to 35 cm/yr at high elevations. These uncertainties amount to about 20 percent of estimated MAP values.     

VERIFICATION OF THE RHEA-OROGRAPHIC-PRECIPITATION MODEL1

ABSTRACT: Observed April 1 snowpack accumulations within and near the Gunnison River basin in southwestern Colorado are compared with simulations from the Rhea-orographic-precipitation model to determine if the model simulates reliable magnitudes and temporal and spatial variability in winter precipitation for the basin. Twenty simulations of the Rhea model were performed using‘optimal’parameter sets determined for 10-kilometer (km) grids (10-km by 10-km grid cells) through stochastic calibration. Comparisons of Rhea-model simulations of winter precipitation with April 1 snowpack accumulations at 32 snowcourse stations were performed for the years 1972–1990. For most stations and most years the Rhea model reliably simulates the temporal and spatial variability in April 1 snowpack accumulations. However, in general, the Rhea-model underestimates April 1 snowpack accumulations in the Gunnison River basin area, and the underestimation is greatest for locations that receive the largest amount of snow. A significant portion of the error in Rhea-model simulations is due to the calibration of the Rhea model using gauge-catch precipitation measurements which can be as much as 50 percent below actual snowfall accumulations. Additional error in the Rhea-model simulations is a result of the comparison of gridded precipitation values to observed values measured at points.     

GIS‐Based Spatial Precipitation Estimation: A Comparison of Geostatistical Approaches1

Abstract: As one of the primary inputs that drive watershed dynamics, the estimation of spatial variability of precipitation has been shown to be crucial for accurate distributed hydrologic modeling. In this study, a Geographic Information System program, which incorporates Nearest Neighborhood (NN), Inverse Distance Weighted (IDW), Simple Kriging (SK), Ordinary Kriging (OK), Simple Kriging with Local Means (SKlm), and Kriging with External Drift (KED), was developed to facilitate automatic spatial precipitation estimation. Elevation and spatial coordinate information were used as auxiliary variables in SKlm and KED methods. The above spatial interpolation methods were applied in the Luohe watershed with an area of 5,239 km², which is located downstream of the Yellow River basin, for estimating 10 years’ (1991‐2000) daily spatial precipitation using 41 rain gauges. The results obtained in this study show that the spatial precipitation maps estimated by different interpolation methods have similar areal mean precipitation depth, but significantly different values of maximum precipitation, minimum precipitation, and coefficient of variation. The accuracy of the spatial precipitation estimated by different interpolation methods was evaluated using a correlation coefficient, Nash‐Sutcliffe efficiency, and relative mean absolute error. Compared with NN and IDW methods that are widely used in distributed hydrologic modeling systems, the geostatistical methods incorporated in this GIS program can provide more accurate spatial precipitation estimation. Overall, the SKlm_EL_X and KED_EL_X, which incorporate both elevation and spatial coordinate as auxiliary into SKlm and KED, respectively, obtained higher correlation coefficient and Nash‐Sutcliffe efficiency, and lower relative mean absolute error than other methods tested. The GIS program developed in this study can serve as an effective and efficient tool to implement advanced geostatistics methods that incorporate auxiliary information to improve spatial precipitation estimation for hydrologic models.     

RAINGAGE NETWORK DESIGN USING NEXRAD PRECIPITATION ESTIMATES1

ABSTRACT: A general framework is proposed for using precipitation estimates from NEXRAD weather radars in raingage network design. NEXRAD precipitation products are used to represent space time rainfall fields, which can be sampled by hypothetical raingage networks. A stochastic model is used to simulate gage observations based on the areal average precipitation for radar grid cells. The stochastic model accounts for subgrid variability of precipitation within the cell and gage measurement errors. The approach is ideally suited to raingage network design in regions with strong climatic variations in rainfall where conventional methods are sometimes lacking. A case study example involving the estimation of areal average precipitation for catchments in the Catskill Mountains illustrates the approach. The case study shows how the simulation approach can be used to quantify the effects of gage density, basin size, spatial variation of precipitation, and gage measurement error, on network estimates of areal average precipitation. Although the quality of NEXRAD precipitation products imposes limitations on their use in network design, weather radars can provide valuable information for empirical assessment of rain‐gage network estimation errors. Still, the biggest challenge in quantifying estimation errors is understanding subgrid spatial variability. The results from the case study show that the spatial correlation of precipitation at subgrid scales (4 km and less) is difficult to quantify, especially for short sampling durations. Network estimation errors for hourly precipitation are extremely sensitive to the uncertainty in subgrid spatial variability, although for storm total accumulation, they are much less sensitive.     

DEVELOPMENT AND APPLICATION OF A SPATIAL HYDROLOGY MODEL OF OKEFENOKEE SWAMP,GEORGIA1

ABSTRACT: :The model described herein was used to assess effects of the Suwannee River sill (a low earthen dam constructed to impound the Suwannee River within the Okefenokee National Wildlife Refuge to eliminate wildfires) on the hydrologic environment of Okefenokee Swamp, Georgia. Developed with Arc/Info Macro Language routines in the GRID environment, the model distributes water in the swamp landscape using precipitation, inflow, evapotranspiration, outflow, and standing water. Water movement direction and rate are determined by the neighborhood topographic gradient, determined using survey grade Global Positioning Systems technology. Model data include flow rates from USGS monitored gauges, precipitation volumes and water levels measured within the swamp, and estimated evapotranspiration volumes spatially modified by vegetation type. Model output in semi‐monthly time steps includes water depth, water surface elevation above mean sea level, and movement direction and volume. Model simulations indicate the sill impoundment affects 18 percent of the swamp during high water conditions when wildfires are scarce and has minimal spatial effect (increasing hydroperiods in less than 5 percent of the swamp) during low water and drought conditions when fire occurrence is high but precipitation and inflow volumes are limited.     

PROBABLE MAXIMUM FLOOD ESTIMATION-EASTERN UNITED STATES1

In 1982, the National Weather Service (NWS) published criteria for developing the spatial and temporal precipitation distribution characteristics of Probable Maximum Storms. The criteria, which are intended for use in the United States east of the 105th meridian, involve four variables: (1) location of the storm center, (2) storm-area size, (3) storm orientation, and (4) temporal arrangement of precipitation amounts. A computer program has been developed which applies the NWS criteria to produce hyetographs of spatially-averaged precipitation for a basin, or for each subbasin if the basin is subdividided. The basis and operational characteristics of the program are described, and an application is illustrated in which the program is used in conjunction with a precipitation-runoff simulation program (HEC-1) to compute a Probable Maximum Flood.     

Generation of Monthly Precipitation Under Climate Change for the Upper Blue Nile River Basin,Ethiopia1

Abstract: This work develops a methodology to project the future precipitation in large river basins under limited data and climate change while preserving the historical temporal and spatial characteristics. The computationally simple and reliable conditional generation method (CGM) is presented and applied to generate reliable monthly precipitation data in the upper Blue Nile River Basin of Ethiopia where rain‐fed agriculture is prevalent. The results showed that the temporal analysis with the CGM performs better to reproduce the historical long‐term characteristics than other methods, and the spatial analysis with the CGM reproduced the historical spatial structure accurately. A 100‐year time series analysis using the outcomes of the six general circulation models showed that precipitation changes by the 2050s (2040 through 2069) can be ?7 to 28% with a mean increase of about 11%. The seasonal results showed increasing wet conditions in all seasons with changes of mean precipitation of 5, 47, and 6% for wet, dry, and mild seasons, respectively.     

Spatial extrapolation of soil characteristics using whole-soil particle size distributions

Soils support ecosystem functions such as plant growth and water quality because of certain physical, chemical, and biological properties. These properties have been studied at different spatial scales, including point scales to satisfy basic research needs, and regional scales to satisfy monitoring needs. Recently, soil property data for the entire USA have become available in the State Soil Geographic Data Base (STATSGO), which is appropriate for regional-scale research. We analyzed and created models of STATSGO data in this study to serve as a research tool, for example, for linking the soil to regional water quality monitoring data in our companion paper. Map units in STATSGO define geographic land areas by soil characteristics (SCs) of similar soil series. We selected 27 SCs that influenced water properties (in varying degrees), aggregated the layer and component SCs to map unit SCs, and used SCs to calculate relationships among map units. The relationships were defined by equations of conditional mean for the qth SC (SCq), while using the remaining 26 SCs as predictors. The relative standard errors for 22 of the 27 SCs were less than 10%, and less than 22% for the remaining five. We conclude that spatial extrapolation of SCs is feasible and the procedures are a first step toward extrapolating information across a region using SC-water property relationships. Although our procedure is for regional scale monitoring, it is also applicable to finer spatial scales commensurate with available soil data.     

EFFECTS OF IMPROVED PRECIPITATION ESTIMATES ON AUTOMATED RUNOFF MAPPING: EASTERN UNITED STATES1

ABSTRACT: We evaluated maps of runoff created by means of two automated procedures. We implemented each procedure using precipitation estimates of both 5-km and 10-km resolution from PRISM (Parameter-elevation Regressions on Independent Slopes Model). Our goal was to determine if using the 5-km PRISM estimates would improve map accuracy. Visual inspection showed good general agreement among our runoff maps, as well as between our maps and one produced using a manual method. A quantitative uncertainty analysis comparing runoff interpolated from our maps with gage data that had been withheld showed slightly smaller actual and percentage interpolation errors for the 5-km PRISM-based maps. Our analyses suggest a modest region-wide improvement in runoff map accuracy with the use of PRISM-based precipitation estimates of 5-km (compared to 10-km) resolution.     

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.     

A New Open‐Access HUC‐8 Based Downscaled CMIP‐5 Climate Model Forecast Dataset for the Conterminous United States

Watershed‐scale hydrologic simulation models generally require climate data inputs including precipitation and temperature. These climate inputs can be derived from downscaled global climate simulations which have the potential to drive runoff forecasts at the scale of local watersheds. While a simulation designed to drive a local watershed model would ideally be constructed at an appropriate scale, global climate simulations are, by definition, arbitrarily determined large rectangular spatial grids. This paper addresses the technical challenge of making climate simulation model results readily available in the form of downscaled datasets that can be used for watershed scale models. Specifically, we present the development and deployment of a new Coupled Model Intercomparison Project phase 5 (CMIP5) based database which has been prepared through a scaling and weighted averaging process for use at the level of U.S. Geological Survey (USGS) Hydrologic Unit Code (HUC)‐8 watersheds. The resulting dataset includes 2,106 virtual observation sites (watershed centroids) each with 698 associated time series datasets representing average monthly temperature and precipitation between 1950 and 2099 based on 234 unique climate model simulations. The new dataset is deployed on a HydroServer and distributed using WaterOneFlow web services in the WaterML format. These methods can be adapted for downscaled General Circulation Model (GCM) results for specific drainage areas smaller than HUC‐8. Two example use cases for the dataset also are presented.     

PRECIPITATION DISTRIBUTION IN COASTAL BRITISH COLUMBIA1

ABSTRACT: The spatial distribution and the temporal and spatial variation of the annual, seasonal, and monthly precipitation in two mountainous watersheds in southwestern British Columbia, Canada, have been analyzed using a detailed database for 1971–1990 in the Capilano and Seymour watersheds. The analysis showed that the precipitation increases up to the mid-position of the watersheds, and then either levels off or decreases. Precipitation on mountain slopes and in the valley at the same distance from the beginning of the slope is similar, and the barrier height is identified as the dominant parameter which influences the precipitation distribution. The temporal variation of the precipitation is the smallest at the mid-position of the watersheds. This variability is the least in the fall and winter and largest in the summer. Correlation between the precipitation accumulations at various stations is large, ranging from 0.80 for the wet period of October-March to 0.65 for the dry period of April-September for distances less than 32 km. Comparison with other studies and the analyses of precipitation and runoff data from coastal British Columbia showed that the results of this study are perhaps general and thus transferable to similar areas in the coastal Pacific Northwest.     

Chemometrical exploration of the wet precipitation chemistry from the Austrian Monitoring Network (1988-1999)

The present paper deals with the application of different chemometric methods to an environmental data set derived from the monitoring of wet precipitation in Austria (1988-1999). These methods are: principal component analysis (PCA); projection pursuit (PP); density-based spatial clustering of application with noise (DBSCAN); ordering points to identify the clustering structures (OPTICS); self-organizing maps (SOM), also called the Kohonen network; and the neural gas (NG) network. The aim of the study is to introduce some new approaches into environmetrics and to compare their usefulness with already existing techniques for the classification and interpretation of environmental data. The density-based approaches give information about the occurrence of natural clusters in the studied data set, which, however, do not occur in the case presented here; information about high-density zones (very similar samples) and extreme samples is also obtained. The partitioning techniques (clustering, but also neural gas and Kohonen networks) offer an opportunity to classify the objects of interest into several defined groups, the patterns of ionic concentration of which can be studied in detail. The visual aids, such as the color map and the Kohonen map, for each site are very helpful in understanding the relationships between samples and between samples and variables. All methods, and in particular projection pursuit, give information about samples with extreme characteristics.     

DEVELOPMENT OF WATERSIIED MODELS FOR TWO SIERRA NEVADA BASINS USING A GEOGRAPHIC INFORMATION SYSTEM1

ABSTRACT: Techniques were developed using vector and raster data in a geographic information system (GIS) to define the spatial variability of watershed characteristics in the north-central Sierra Nevada of California and Nevada and to assist in computing model input parameters. The U.S. Geological Survey's Precipitation-Runoff Modeling System, a physically based, distributed-parameter watershed model, simulates runoff for a basin by partitioning a watershed into areas that each have a homogeneous hydrologic response to precipitation or snowmelt. These land units, known as hydrologic-response units (HRU's), are characterized according to physical properties, such as altitude, slope, aspect, land cover, soils, and geology, and climate patterns. Digital data were used to develop a GIS data base and HRIJ classification for the American River and Carson River basins. The following criteria are used in delineating HRU's: (1) Data layers are hydrologically significant and have a resolution appropriate to the watershed's natural spatial variability, (2) the technique for delineating HRU's accommodates different classification criteria and is reproducible, and (3) HRU's are not limited by hydrographic-subbasin boundaries. HRU's so defined are spatially noncontiguous. The result is an objective, efficient methodology for characterizing a watershed and for delineating HRU's. Also, digital data can be analyzed and transformed to assist in defining parameters and in calibrating the model.     

Urbanization and Its Effect On Runoff in the Whiteoak Bayou Watershed,Texas1

Abstract: The capacity of a watershed to urbanize without changing its hydrologic response and the relationship between that response and the spatial configuration of the developed areas was studied. The study was conducted in the Whiteoak Bayou watershed (223 km²), located northwest of Houston, Texas, over an analysis period from 1949 to 2000. Annual development data were derived from parcel data collected by the Harris County Appraisal District. Using these data, measures of the spatial configuration of the watershed urban areas were calculated for each year. Based on regression models, it was determined that the annual runoff depths and annual peak flows depended on the annual precipitation depth, the developed area and the maximum 12‐h precipitation depth on the day and day before the peak flow took place. It was found that, since the early 1970s, when the watershed reached a 10% impervious area, annual runoff depths and peak flows have increased by 146% and 159%, respectively. However, urbanization is responsible for only 77% and 32% of the increase, respectively, while precipitation changes are responsible for the remaining 39% and 96%, respectively. Likewise, an analysis of the development data showed that, starting in the early 1970s, urbanization in the watershed consisted more of connecting already developed areas than of creating new ones, which increases the watershed’s conveyance capacity and explains the change in its response. Before generalizing conclusions, though, further research on other urban watersheds with different urbanization models appears to be necessary.     

Annual Estimates of Recharge,Quick‐Flow Runoff,and Evapotranspiration for the Contiguous U.S. Using Empirical Regression Equations

This study presents new data‐driven, annual estimates of the division of precipitation into the recharge, quick‐flow runoff, and evapotranspiration (ET) water budget components for 2000‐2013 for the contiguous United States (CONUS). The algorithms used to produce these maps ensure water budget consistency over this broad spatial scale, with contributions from precipitation influx attributed to each component at 800 m resolution. The quick‐flow runoff estimates for the contribution to the rapidly varying portion of the hydrograph are produced using data from 1,434 gaged watersheds, and depend on precipitation, soil saturated hydraulic conductivity, and surficial geology type. Evapotranspiration estimates are produced from a regression using water balance data from 679 gaged watersheds and depend on land cover, temperature, and precipitation. The quick‐flow and ET estimates are combined to calculate recharge as the remainder of precipitation. The ET and recharge estimates are checked against independent field data, and the results show good agreement. Comparisons of recharge estimates with groundwater extraction data show that in 15% of the country, groundwater is being extracted at rates higher than the local recharge. These maps of the internally consistent water budget components of recharge, quick‐flow runoff, and ET, being derived from and tested against data, are expected to provide reliable first‐order estimates of these quantities across the CONUS, even where field measurements are sparse.     

AN EVALUATION OF PRECIPITATION GAGE DENSITY IN A MOUNTAINOUS TERRAIN1

ABSTRACT: Inter-station analysis was employed to evaluate the adequacy of the precipitation network in topographically complex West Virginia. A 25-year period was determined as the minimum lingth of record needed for relatively stable and fairly accurate estimates of long-term (50-year) precipitation and in frequency analysis. Data from the 83 National Weather Service stations with 25-year records were adjusted for consistency and evaluated separately by zones east (31 stations) and west (52 stations) of the Appalachian divide. Correlation coefficients (r) and average standard errors of estimate were computed for all station pairs within 50 miles distance and 1000 feet elevation difference of each other. The third polynomial equation of inter-station distance eliminated using elevation and land slope as the criteria in network design in this mountainous terrain. A network with (r) = 0.9 estimates annual precipitation with accuracy as great as 5 percent, but requires about 250 additional gages (i.e., about 200 percent of the present density).     

宁夏近51年气候变化特征分析

利用宁夏10个气象站1960-2010年日气温和降水资料,对宁夏气温和降水的时空变化特征进行了分析.结果显示:①宁夏整体呈现暖干化趋势,但存在区域间差异.②近51年来宁夏年平均气温呈上升趋势,冬季升温幅度大,秋季升温幅度最小.③年平均降水量总体呈现减小趋势,春季微弱增加,冬季降水量显著增加,夏季降水微弱减少,秋季降水减少明显.④通过Kriging空间差值法分析表明,年平均气温总体呈现由西南向西北递增的趋势,中、北部增温较明显；平均年总降水量总体呈现由西南向西北递减的规律,中、北部降水量减小最为明显.