Using SWAT to Model Streamflow in Two River Basins With Ground and Satellite Precipitation Data1

Abstract: Both ground rain gauge and remotely sensed precipitation (Next Generation Weather Radar – NEXRAD Stage III) data have been used to support spatially distributed hydrological modeling. This study is unique in that it utilizes and compares the performance of National Weather Service (NWS) rain gauge, NEXRAD Stage III, and Tropical Rainfall Measurement Mission (TRMM) 3B42 (Version 6) data for the hydrological modeling of the Middle Nueces River Watershed in South Texas and Middle Rio Grande Watershed in South Texas and northern Mexico. The hydrologic model chosen for this study is the Soil and Water Assessment Tool (SWAT), which is a comprehensive, physical‐based tool that models watershed hydrology and water quality within stream reaches. Minor adjustments to selected model parameters were applied to make parameter values more realistic based on results from previous studies. In both watersheds, NEXRAD Stage III data yields results with low mass balance error between simulated and actual streamflow (±13%) and high monthly Nash‐Sutcliffe efficiency coefficients (NS > 0.60) for both calibration (July 1, 2003 to December 31, 2006) and validation (2007) periods. In the Middle Rio Grande Watershed NEXRAD Stage III data also yield robust daily results (time averaged over a three‐day period) with NS values of (0.60‐0.88). TRMM 3B42 data generate simulations for the Middle Rio Grande Watershed of variable qualtiy (MBE = +13 to ?16%; NS = 0.38‐0.94; RMSE = 0.07‐0.65), but greatly overestimates streamflow during the calibration period in the Middle Nueces Watershed. During the calibration period use of NWS rain gauge data does not generate acceptable simulations in both watersheds. Significantly, our study is the first to successfully demonstrate the utility of satellite‐estimated precipitation (TRMM 3B42) in supporting hydrologic modeling with SWAT; thereby, potentially extending the realm (between 50°N and 50°S) where remotely sensed precipitation data can support hydrologic modeling outside of regions that have modern, ground‐based radar networks (i.e., much of the third world).     

Long‐Term High‐Resolution Radar Rainfall Fields for Urban Hydrology

Accurate records of high‐resolution rainfall fields are essential in urban hydrology, and are lacking in many areas. We develop a high‐resolution (15 min, 1 km²) radar rainfall data set for Charlotte, North Carolina during the 2001‐2010 period using the Hydro‐NEXRAD system with radar reflectivity from the National Weather Service Weather Surveillance Radar 1988 Doppler weather radar located in Greer, South Carolina. A dense network of 71 rain gages is used for estimating and correcting radar rainfall biases. Radar rainfall estimates with daily mean field bias (MFB) correction accurately capture the spatial and temporal structure of extreme rainfall, but bias correction at finer timescales can improve cold‐season and tropical cyclone rainfall estimates. Approximately 25 rain gages are sufficient to estimate daily MFB over an area of at least 2,500 km², suggesting that robust bias correction is feasible in many urban areas. Conditional (rain‐rate dependent) bias can be removed, but at the expense of other performance criteria such as mean square error. Hydro‐NEXRAD radar rainfall estimates are also compared with the coarser resolution (hourly, 16 km²) Stage IV operational rainfall product. Stage IV is adequate for flood water balance studies but is insufficient for applications such as urban flood modeling, in which the temporal and spatial scales of relevant hydrologic processes are short. We recommend the increased use of high‐resolution radar rainfall fields in urban hydrology.     

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.     

Physically Based,Hydrologic Model Results Based on Three Precipitation Products1

Abstract: The main objective of the study is to examine the accuracy of and differences among simulated streamflows driven by rainfall estimates from a network of 22 rain gauges spread over a 2,170 km² watershed, NEXRAD Stage III radar data, and Tropical Rainfall Measuring Mission (TRMM) 3B42 satellite data. The Gridded Surface Subsurface Hydrologic Analysis (GSSHA), a physically based, distributed parameter, grid‐structured, hydrologic model, was used to simulate the June‐2002 flooding event in the Upper Guadalupe River watershed in south central Texas. There were significant differences between the rainfall fields estimated by the three types of measurement technologies. These differences resulted in even larger differences in the simulated hydrologic response of the watershed. In general, simulations driven by radar rainfall yielded better results than those driven by satellite or rain‐gauge estimates. This study also presents an overview of effects of land cover changes on runoff and stream discharge. The results demonstrate that, for major rainfall events similar to the 2002 event, the effect of urbanization on the watershed in the past two decades would not have made any significant effect on the hydrologic response. The effect of urbanization on the hydrologic response increases as the size of the rainfall event decreases.     

Evaluation of Gauge‐Radar Merging Methods Using a Semi‐Distributed Hydrological Model in the Upper Thames River Basin,Canada

Gauge‐radar merging methods combine rainfall estimates from rain gauges and radar to capitalize on the strengths of the individual instruments. The performance of four well‐known gauge‐radar merging methods, including mean field bias correction, Brandes spatial adjustment, local bias correction using kriging, and conditional merging, are examined using Environment Canada radar and the Upper Thames River Basin in southwestern Ontario, Canada, as a case study. The analysis assesses the effect of gauge‐radar merging methods on: (1) the accuracy of predicted rainfall accumulations; and (2) the accuracy of predicted streamflows using a semi‐distributed hydrological model. In addition, several influencing factors (i.e., gauge density, storm type, basin type, proximity to the radar tower, and time‐step of adjustment) are analyzed to determine their effect on the performance of the rainfall estimation techniques. Confirming results of previous studies, the merging methods provide an increase in the accuracy of both rainfall accumulation estimations and predicted streamflows. The results also indicate specific factors such as gauge density, rainfall intensity, and time‐step of adjustment can reduce the accuracy of merging methods and play a key role in the examination of its use for operational purposes. Results provide guidance for hydrologists and engineers assessing how best to apply corrected radar products to improve rainfall estimation and hydrological modeling accuracy.     

Satellite Rainfall Estimates Over South America – Possible Applicability to the Water Management of Large Watersheds1

Pereira Filho, Augusto J., Richard E. Carbone, John E. Janowiak, Phillip Arkin, Robert Joyce, Ricardo Hallak, and Camila G.M. Ramos, 2010. Satellite Rainfall Estimates Over South America – Possible Applicability to the Water Management of Large Watersheds. Journal of the American Water Resources Association (JAWRA) 46(2):344-360. DOI: 10.1111/j.1752-1688.2009.00406.x Abstract: This work analyzes high-resolution precipitation data from satellite-derived rainfall estimates over South America, especially over the Amazon Basin. The goal is to examine whether satellite-derived precipitation estimates can be used in hydrology and in the management of larger watersheds of South America. High spatial-temporal resolution precipitation estimates obtained with the CMORPH method serve this purpose while providing an additional hydrometeorological perspective on the convective regime over South America and its predictability. CMORPH rainfall estimates at 8-km spatial resolution for 2003 and 2004 were compared with available rain gauge measurements at daily, monthly, and yearly accumulation time scales. The results show the correlation between satellite-derived and gauge-measured precipitation increases with accumulation period from daily to monthly, especially during the rainy season. Time-longitude diagrams of CMORPH hourly rainfall show the genesis, strength, longevity, and phase speed of convective systems. Hourly rainfall analyses indicate that convection over the Amazon region is often more organized than previously thought, thus inferring that basin scale predictions of rainfall for hydrological and water management purposes have the potential to become more skillful. Flow estimates based on CMORPH and the rain gauge network are compared to long-term observed average flow. The results suggest this satellite-based rainfall estimation technique has considerable utility. Other statistics for monthly accumulations also suggest CMORPH can be an important source of rainfall information at smaller spatial scales where in situ observations are lacking.     

A PROPOSED METHOD FOR DROUGHT MONITORING1

ABSTRACT: Existing definitions of drought have focused on limited hydrologic indicators and are less effective for the purpose of drought monitoring. This study uses historical records of streamflow, precipitation, ground water, temperature, and lake elevation to define drought. Based on the method of truncation, drought durations and conditional probabilities of each indicator were estimated to define the drought severity levels, namely, 70 percent, 80 percent, 90 percent, and 95 percent. A drought monitoring method was developed by a combination of truncation level, duration, and conditional probabilities of five indicators. A six-month period of the 1988 drought in the central Ohio region was used to test the monitoring method. It was found that the developed method could effectively detect an occurrence of drought.     

Analysis of Meteorological Drought Pattern During Different Climatic and Cropping Seasons in Bangladesh

Drought is one of the most frequent natural disasters in Bangladesh which severely affect agro‐based economy and people's livelihood in almost every year. Characterization of droughts in a systematic way is therefore critical in order to take necessary actions toward drought mitigation and sustainable development. In this study, standardized precipitation index is used to understand the spatial distribution of meteorological droughts during various climatic seasons such as premonsoon, monsoon, and winter seasons as well as cropping seasons such as Pre‐Kharif (March‐May), Kharif (May‐October), and Rabi (December‐February). Rainfall data collected from 29 rainfall gauge stations located in different parts of the country were used for a period of 50 years (1961‐2010). The study reveals that the spatial characteristics of droughts vary widely according to season. Premonsoon droughts are more frequent in the northwest, monsoon droughts mainly occur in the west and northwest, winter droughts in the west, and the Rabi and Kharif droughts are more frequent in the north and northwest of Bangladesh. It is expected that the findings of the study will support drought monitoring and mitigation activities in Bangladesh.     

A COMPARISON OF METHODS FOR CALCULATION OF RADAR-RAINFALL HOURLY ACCUMULATIONS1

ABSTRACT: We compared two interpolation schemes for calculation of hourly accumulation of radar-rainfall. The schemes are: (1) the Advection Method, and (2) the Space-Time Kriging Method. The performance of the methods is investigated using numerical simulation experiments. Space-time evolution of rainfall fields is generated from a stochastic model. The generated fields are sampled following typical radar scanning strategies, and the investigated schemes are applied to obtain accumulated rainfall patterns. The statistical results and a visual analysis of the graphical images suggest that it is advisable to use an interpolation scheme for radar observations even when storm velocity is not high. The Space-Time Kriging Method provides the best results for low wind velocity. The Advection Method has the smallest standard deviation and mean absolute error, and preserves well the true rainfall pattern for high wind velocity.     

Regional scale flood modeling using NEXRAD rainfall, GIS, and HEC-HMS/RAS: a case study for the San Antonio River Basin Summer 2002 storm event

This paper develops a framework for regional scale flood modeling that integrates NEXRAD Level III rainfall, GIS, and a hydrological model (HEC-HMS/RAS). The San Antonio River Basin (about 4000 square miles, 10,000 km2) in Central Texas, USA, is the domain of the study because it is a region subject to frequent occurrences of severe flash flooding. A major flood in the summer of 2002 is chosen as a case to examine the modeling framework. The model consists of a rainfall-runoff model (HEC-HMS) that converts precipitation excess to overland flow and channel runoff, as well as a hydraulic model (HEC-RAS) that models unsteady state flow through the river channel network based on the HEC-HMS-derived hydrographs. HEC-HMS is run on a 4 x 4 km grid in the domain, a resolution consistent with the resolution of NEXRAD rainfall taken from the local river authority. Watershed parameters are calibrated manually to produce a good simulation of discharge at 12 subbasins. With the calibrated discharge, HEC-RAS is capable of producing floodplain polygons that are comparable to the satellite imagery. The modeling framework presented in this study incorporates a portion of the recently developed GIS tool named Map to Map that has been created on a local scale and extends it to a regional scale. The results of this research will benefit future modeling efforts by providing a tool for hydrological forecasts of flooding on a regional scale. While designed for the San Antonio River Basin, this regional scale model may be used as a prototype for model applications in other areas of the country.     

USE OF THE KRIGING METHOD FOR STUDYING CHARACTERISTICS OF GROUND WATER DROUGHTS1

round water drought events were derived by taking a truncation level through the time series of daily ground water depth that are recorded elevation differences between the water table and land surface at a well site. Droughts of various truncation levels at 70, 80, 90, and 95 percent, were obtained, where a 70 percent truncation level means that 70 percent of ground water depth data are less than or equal to the truncated value. The conditional probability that a drought occurring at a certain truncation level will prolong and advance to that of the next higher level was estimated. The regionalization analysis was conducted assuming that conditional probabilities estimated at selected wells are regionalized variables. Contour lines of conditional probabilities for each truncation level were constructed to express their spatial variability in the region. Estimation errors associated with the regionalization were reasonably small.     

Validation of Satellite Precipitation Adjustment Methodology From Seven Basins in the Continental United States1

Tobin, Kenneth J. and Marvin E. Bennett, 2012. Validation of Satellite Precipitation Adjustment Methodology From Seven Basins in the Continental United States. Journal of the American Water Resources Association (JAWRA) 48(2): 221‐234. DOI: 10.1111/j.1752‐1688.2011.00604.x Abstract: The precipitation science community has expressed concern regarding the ability of satellite‐based precipitation products to accurately capture rainfall values over land. There has been some work that has focused on addressing the deficiencies of satellite precipitation products, particularly on the adjustment of bias. This article outlines a methodology that adjusts satellite products utilizing ground‐based precipitation data. The approach is not a simple bias adjustment, but is a three‐step process that transforms a satellite product based on a ground‐based precipitation product (NEXRAD‐derived Multisensor Precipitation Estimator [MPE] product or rain‐gauge data). The developed methodology was successfully applied to seven moderate‐to‐large sized watersheds from continental United States (CONUS) and northern Mexico over a spectrum of climatic regimes ranging from dry to humid settings. Methodology validation is based on comparison of observed and simulated streamflow generated with SWAT (Soil and Water Assessment Tool) model using unadjusted and adjusted precipitation products as input. Streamflow comparison is based on mass balance error and Nash‐Sutcliffe efficiency coefficient. Finally, the contribution of how adjustment to correct misses, false alarms, and bias impacts adjusted datasets and the potential impact that the adjustment methodology can have on hydrological applications such as water resource monitoring and flood prediction are explored.     

COMBINATION OF RADAR AND GAUGE DATA IN A RAINFALL ARCHIVE SYSTEM1

ABSTRACT: Near real time daily rainfall estimates for the UK are available from three sources: a sparse network of gauges, radar data, or radar data adjusted by the sparse gauges. The PARAGON rainfall archive system, which has been developed by the UK Meteorological Office, is able to produce these estimates in near real time on a 5 km grid. The ability of these estimates to reproduce the 5 km grid point field derived later from a dense network of gauges is compared using case studies. Five techniques have been used to assess the relative quality of the various estimates. There is general agreement between the results of the various techniques. For the London radar there are examples of days when the rainfall estimate was improved by incorporating radar data; conversely, there are days when the radar data make it worse. Overall little evidence was found to suggest that adjusted radar data are consistently markedly better than gauge estimates. Discriminate use of radar data is recommended.     

MODIS Biophysical States and NEXRAD Precipitation in a Statistical Evaluation of Antecedent Moisture Condition and Streamflow1

Abstract: The potential of remotely sensed time series of biophysical states of landscape to characterize soil moisture condition antecedent to radar estimates of precipitation is assessed in a statistical prediction model of streamflow in a 1,420 km² watershed in south‐central Texas, Moderate Resolution Imaging Spectroradiometer (MODIS) time series biophysical products offer significant opportunities to characterize and quantify hydrologic state variables such as land surface temperature (LST) and vegetation state and status. Together with Next Generation Weather Radar (NEXRAD) precipitation estimates for the period 2002 through 2005, 16 raw and deseasoned time series of LST (day and night), vegetation indices, infrared reflectances, and water stress indices were linearly regressed against observed watershed streamflow on an eight‐day aggregated time period. Time offsets of 0 (synchronous with streamflow event), 8, and 16 days (leading streamflow event) were assessed for each of the 16 parameters to evaluate antecedent effects. The model results indicated a reasonable correlation (r² = 0.67) when precipitation, daytime LST advanced 16 days, and a deseasoned moisture stress index were regressed against log‐transformed streamflow. The estimation model was applied to a validation period from January 2006 through March 2007, a period of 12 months of regional drought and base‐flow conditions followed by three months of above normal rainfall and a flood event. The model resulted in a Nash‐Sutcliffe estimation efficiency (E) of 0.45 for flow series (in log‐space) for the full 15‐month period, ?0.03 for the 2006 drought condition period, and 0.87 for the 2007 wet condition period. The overall model had a relative volume error of ?32%. The contribution of parameter uncertainties to model discrepancy was evaluated.     

Hydrologic Model Framework for Water Resources Planning in the Santa Cruz River,Southern Arizona1

Abstract: The authors develop a model framework that includes a set of hydrologic modules as a water resources management and planning tool for the upper Santa Cruz River near the Mexican border, Southern Arizona. The modules consist of: (1) stochastic generation of hourly precipitation scenarios that maintain the characteristics and variability of a 45‐year hourly precipitation record from a nearby rain gauge; (2) conceptual transformation of generated precipitation into daily streamflow using varied infiltration rates and estimates of the basin antecedent moisture conditions; and (3) surface‐water to ground‐water interaction for four downstream microbasins that accounts for alluvial ground‐water recharge, and ET and pumping losses. To maintain the large inter‐annual variability of streamflow as prevails in Southern Arizona, the model framework is constructed to produce three types of seasonal winter and summer categories of streamflow (i.e., wet, medium, or dry). Long‐term (i.e., 100 years) realizations (ensembles) are generated by the above described model framework that reflects two different regimes of inter annual variability. The first regime is that of the historic streamflow gauge record. The second regime is that of the tree ring reconstructed precipitation, which was derived for the study location. Generated flow ensembles for these two regimes are used to evaluate the risk that the regional four ground‐water microbasins decline below a preset storage threshold under different operational water utilization scenarios.     

Drivers of future streamflow changes in watersheds across the Northeastern United States

Accurate projections of streamflow, which have implications for flooding, water resources, hydropower, and ecosystems, are critical to climate change adaptation and require an understanding of streamflow sensitivity to climate drivers. The northeastern United States has experienced a dramatic increase in extreme precipitation over the past 25 years; however, the effects of these changes, as well as changes in other drivers of streamflow, remain unclear. Here, we use a random forest model forced with a regional climate model to examine historical and future streamflow dynamics of four watersheds across the Northeast. We find that streamflow in the cold season (November–May) is primarily driven by 3-day rainfall and antecedent wetness (Antecedent Precipitation Index) in three rainfall-dominant watersheds, and 30-day rainfall, antecedent wetness, and 30-day snowmelt in the fourth, more snowmelt-dominated watershed. In the warm season (June–October), streamflow is driven by antecedent wetness and rainfall in all watersheds. By the end of the century (2070–2099), cold season streamflow depends on the importance placed on snow in the machine learning model, with changes ranging from −7% (with snow) to +40% (without snow) in a single watershed. Simulated future warm season streamflow increases in two watersheds (56% and 193%) due to increased precipitation and antecedent soil wetness, but decreases in the other two watersheds (−6% and −27%) due to reduced precipitation.     

LONG TERM TRENDS OF ANNUAL AND MONTHLY PRECIPITATION IN JAPAN1

ABSTRACT: Long term trends in Japan's annual and monthly precipitation are investigated in this study. The statistical significance of a trend at a study site is assessed by the Mann‐Kendall (MK) test, and field significance of trends in climatic Regions II, III, and IV is evaluated using the bootstrap test preserving cross correlation. The practical significance of a trend is judged by a percentage change of the sample mean over an observation period. The field significance assessment demonstrates that annual precipitation in Region II did not show any significant change, but regional precipitation shifts occurred in different months. Precipitation significantly increased by 12.2 percent in May, while it significantly decreased by 12.0, 10.5, 15.6, and 19.7 percent, respectively, in April, September, October, and December. In Region III, annual precipitation declined by 11.8 percent, and monthly precipitation significantly decreased from September through January and in April, with the greatest decrease (38.2 percent) in December. In Region IV, significant reductions occurred in both annual precipitation (by 15.6 percent) and monthly precipitation from September through February and in June and July, with the worst reduction (44.7 percent) in December.     

Assessment of the Extreme Rainfall Event at Nashville,TN and the Surrounding Region on May 1–3, 2010

This paper analyzes the May 1–3, 2010 rainfall event that affected the south‐central United States, including parts of Mississippi, Tennessee, and Kentucky. The storm is evaluated in terms of its synoptic setting, along with the temporal distributions, and spatial patterns of the rainfall. In addition, the recurrence interval of the storm is assessed and the implications for hydrologic structure designs are discussed. The event was associated with an upper‐level trough and stationary frontal boundary to the west of the rainfall region, which remained quasi‐stationary for a period of 48 h. Heavy rainfall was produced by two slow‐moving mesoscale convective complexes, combined with abundant atmospheric moisture. Storm totals exceeding 330 mm occurred within a large elongated area extending from Memphis to Nashville. Isolated rainfall totals over 480 mm were reported in some areas, with NEXRAD weather radar rainfall estimates up to 501 mm. An extreme value analysis was performed for one‐ and two‐day rainfall totals at Nashville and Brownsville, Tennessee, as well as for gridded rainfall estimates for the entire region using the Storm Precipitation Analysis System. Results suggest maximum rainfall totals for some durations during the May 1–3, 2010 event exceeded the 1,000‐year rainfall values from National Oceanic and Atmospheric Administration Atlas 14 for a large portion of the region and reached up to 80% of the probable maximum precipitation values for some area sizes and durations.     

Probabilistic Streamflow Generation Model for Data Sparse Arid Watersheds1

Abstract: The authors present a model that generates streamflow for ephemeral arid streams. The model consists of a stochastic hourly precipitation point process model and a conceptual model that transforms precipitation into flow. It was applied to the Santa Cruz River at the border crossing from Mexico into Southern Arizona. The model was constructed for four different seasons and three categories of inter‐annual variability for the wet seasons of summer and winter. The drainage area is ungauged and precipitation information was inferred from a precipitation gauge downstream. The precipitation gauge record was evaluated against simulated precipitation from a mesoscale numerical weather prediction model, and was found to be the representative of the regional precipitation variability. The flow generation was found to reproduce the variability in the observed record at the daily, seasonal and annual time scales, and it is suitable for use in planning studies for the study site.