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.     

A STRIP MODEL APPROACH TO PARAMETERIZE A COUPLED GREEN‐AMPT KINEMATIC WAVE MODEL1

ABSTRACT: Infiltration processes at the plot scale are often described and modeled using a single effective hydraulic conductivity (Kg) value. This can lead to errors in runoff and erosion prediction. An integrated field measurement and modeling study was conducted to evaluate: (1) the relationship among rainfall intensity, spatially variable soil and vegetation characteristics, and infiltration processes; and (2) how this relationship could be modeled using Green and Ampt and a spatially distributed hydrologic model. Experiments were conducted using a newly developed variable intensity rainfall simulator on 2 m by 6 m plots in a rangeland watershed in southeastern Arizona. Rainfall application rates varied between 50 and 200 mm/hr. Results of the rainfall simulator experiments showed that the observed hydrologic response changed with changes in rainfall intensity and that the response varied with antecedent moisture condition. A distributed process based hydrologic simulation model was used to model the plots at different levels of hydrologic complexity. The measurement and simulation model results show that the rainfall runoff relationship cannot be accurately described or modeled using a single Kg value at the plot scale. Multi‐plane model configurations with infiltration parameters based on soil and plot characteristics resulted in a significant improvement over single‐plane configurations.     

A NOTE ON AREAL RAINFALL DEFINITION

A discussion is presented of the likely sources of error in defining areal rainfall on a storm basis. These include the instrumental error, sampling fluctuations over the area, and network density. The analysis of dense raingage data provides some perspective of the magnitude of the errors that might be encountered from the natural variability of rainfall. Except for one watershed in Arizona, the coefficient of variation, based on a sample of storm totals from the individual gages in various size areas, remains relatively constant with increasing area for a particular storm. The error due to rainfall variability over the area is probably the most important and must be considered in experiments which attempt to resolve small-area hydrologic problems.     

EVALUATING THE ACCURACY OF RAINFALL CATCH BY THREE DIFFERENT GAGES1

ABSTRACT: A U.S. standard gage, a weighing-type recording gage, a standard gage fitted with an Alter windshield, and a pit gage were installed to evaluate the accuracy and wind effects on rainfall catch. The study was conducted at the Stephen F. Austin Experimental Forest, about 20 km SW of Nacogdoches, Texas. A recording anemometer was also installed at a height corresponding to the standard gage orifice. Based on data from 67 storms collected over a one-year period (July 1995-August 1996), all three conventional gages consistently caught less rainfall than the reference pit gage with an average percent deficiency greater than 10 percent. However, the recording gage caught 2.7 percent less and the shielded gage caught 1 percent more than the standard gage—differences less than those reported elsewhere. The deficiencies were highly correlated with storm intensity, duration, or total rainfall. When the correction for wind effect on angle of raindrop inclination is included, the percent catch deficiency of the standard gage was reduced from 11 percent to 6 percent. The remaining errors may be attributed to wind effects (streamline vs. turbulent flow), nonrandom errors, or other unknown sources.     

Estimating Lag to Peak between Rainfall and Peak Streamflow with a Mixed‐Effects Model

We test the use of a mixed‐effects model for estimating lag to peak for small basins in Maine (drainage areas from 0.8 to 78 km²). Lag to peak is defined as the time between the center of volume of the excess rainfall during a storm event and the resulting peak streamflow. A mixed‐effects model allows for multiple observations at sites without violating model assumptions inherent in traditional ordinary least squares models, which assume each observation is independent. The mixed model includes basin drainage area and maximum 15‐min rainfall depth for individual storms as explanatory features. Based on a remove‐one‐site cross‐validation analysis, the prediction errors of this model ranged from ?42% to +73%. The mixed model substantially outperformed three published models for lag to peak and one published model for centroid lag for estimating lag to peak for small basins in Maine. Lag to peak estimates are a key input to rainfall–runoff models used to design hydraulic infrastructure. The improved accuracy and consistency with model assumptions indicates that mixed models may provide increased data utilization that could enhance models and estimates of lag to peak in other regions.     

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.     

ERRORS IN ESTIMATED STREAMFLOW PARAMETERS AND STORAGES FOR UNGAUGED CATCHMENTS1

ABSTRACT: For a set of 81 catchments in southeast Victoria, Australia, predictive equations were developed by least squares regression of the mean and coefficient of variation of annual Streamflow against a variety of rainfall and physiographic parameters. The data were also divided into subsets according to catchment size, subregion, or record length of investigate if the relationships differed significantly between subsets. Only the catchment area and some rainfall statistical parameters were found to be significant. Streamflow parameters predicted by the regression equations were used to estimate storage requirements in ungauged catchments. The influence of errors in the Streamflow parameters on the storage error was examined.     

Nonparametric Monte Carlo Simulation for Flood Frequency Curve Derivation: An Application to a Korean Watershed1

Abstract: A mix of causative mechanisms may be responsible for flood at a site. Floods may be caused because of extreme rainfall or rain on other rainfall events. The statistical attributes of these events differ according to the watershed characteristics and the causes. Traditional methods of flood frequency analysis are only adequate for specific situations. Also, to address the uncertainty of flood frequency estimates for hydraulic structures, a series of probabilistic analyses of rainfall‐runoff and flow routing models, and their associated inputs, are used. This is a complex problem in that the probability distributions of multiple independent and derived random variables need to be estimated to evaluate the probability of floods. Therefore, the objectives of this study were to develop a flood frequency curve derivation method driven by multiple random variables and to develop a tool that can consider the uncertainties of design floods. This study focuses on developing a flood frequency curve based on nonparametric statistical methods for the estimation of probabilities of rare floods that are more appropriate in Korea. To derive the frequency curve, rainfall generation using the nonparametric kernel density estimation approach is proposed. Many flood events are simulated by nonparametric Monte Carlo simulations coupled with the center Latin hypercube sampling method to estimate the associated uncertainty. This study applies the methods described to a Korean watershed. The results provide higher physical appropriateness and reasonable estimates of design flood.     

FLOOD EVENT MODELING - A STUDY OF TWO METHODS1

: Despite the advances in catchment modeling in recent years, engineers still face major problems in estimating flood flows. Application of unit hydrograph and runoff routing models to five United Kingdom catchments shows that either can be tuned to predict, on a test event, the routing effects of a catchment with equal accuracy. The larger remaining problem is the prediction of losses from rainfall and this study shows how alternative ways of describing the within event distribution of these losses can be an important factor controlling the success of the overall model. Other problems include the risks of extrapolation to larger events, baseflow separation methods, rainfall patterns, and inevitable errors in the data.     

Use of Predictive Weather Uncertainties in an Irrigation Scheduling Tool Part II: An Application of Metrics and Adjoints

We apply predictive weather metrics and land model sensitivities to improve the Colorado State University Water Irrigation Scheduler for Efficient Application (WISE). WISE is an irrigation decision aid that integrates environmental and user information for optimizing water use. Rainfall forecasts and verification performance metrics are used to estimate predictive rainfall probabilities that are used as input data within the irrigation decision aid. These input data errors are also used within a land model sensitivity study to diagnose important prognostic water movement behaviors for irrigation tool development purposes simultaneously performing the analysis in space and time. Thus, important questions such as “how long can a crop water application be delayed while maintaining crop yield production?” are addressed by evaluating crop growth stage interactions as a function of soil depth (i.e., space), rainfall events (i.e., time), and their probabilistic uncertainties. 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.     

Hydrologic Modeling Uncertainty Resulting From Land Cover Misclassification1

Abstract: A stochastic, spatially explicit method for assessing the impact of land cover classification error on distributed hydrologic modeling is presented. One‐hundred land cover realizations were created by systematically altering the North American Landscape Characterization land cover data according to the dataset’s misclassification matrix. The matrix indicates the probability of errors of omission in land cover classes and is used to assess the uncertainty in hydrologic runoff simulation resulting from parameter estimation based on land cover. These land cover realizations were used in the GIS‐based Automated Geospatial Watershed Assessment tool in conjunction with topography and soils data to generate input to the physically‐based Kinematic Runoff and Erosion model. Uncertainties in modeled runoff volumes resulting from these land cover realizations were evaluated in the Upper San Pedro River basin for 40 watersheds ranging in size from 10 to 100 km² under two rainfall events of differing magnitudes and intensities. Simulation results show that model sensitivity to classification error varies directly with respect to watershed scale, inversely to rainfall magnitude and are mitigated or magnified by landscape variability depending on landscape composition.     

EFFECT OF RAINFALL EXCESS CALCULATIONS ON MODELED HYDROGRAPH ACCURACY AND UNIT-HYDROGRAPH PARAMETERS1

ABSTRACT: Two methods of computing rainfall excess in the U.S. Army Corps of Engineers’flood hydrograph package (HEC-1), the Initial and Uniform method and the Exponential method, are compared to evaluate the effects on modeled hydrograph accuracy. Two computed unit-hydrograph parameters, time of concentration and storage coefficient, were also compared. Rainfall and runoff data from 209 storms in 32 gaged basins in Illinois were used to calibrate the HEC-1 model. Three hydrograph characteristics - sum of incremental flows, peak discharge, and time of peak discharge - were used to evaluate modeled hydrograph accuracy. Mean percent error for each basin and hydrograph characteristic was computed. An evaluation of the mean errors indicates that, although some bias in modeled hydrograph accuracy is evident, rainfall excess computed using either method results in a computed hydrograph accuracy that is within generally accepted limits. Application of a linear-regression model shows no significant differences in computed values of unit-hydrograph parameters.     

CALCULATION OF NONLINEAR CONFIDENCE AND PREDICTION INTERVALS FOR GROUND-WATER FLOW MODELS1

ABSTRACT: A method is derived to efficiently compute nonlinear confidence and prediction intervals on any function of parameters derived as output from a mathematical model of a physical system. The method is applied to the problem of obtaining confidence and prediction intervals for manually-calibrated ground-water flow models. To obtain confidence and prediction intervals resulting from uncertainties in parameters, the calibrated model and information on extreme ranges and ordering of the model parameters within one or more independent groups are required. If random errors in the dependent variable are present in addition to uncertainties in parameters, then calculation of prediction intervals also requires information on the extreme range of error expected. A simple Monte Carlo method is used to compute the quantiles necessary to establish probability levels for the confidence and prediction intervals. Application of the method to a hypothetical example showed that includsion of random errors in the dependent variable in addition to uncertainties in parameters can considerably widen the prediction intervals.     

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.     

ASSESSMENT OF SAMPLING ERROR ASSOCIATED WITH SOIL MOISTURE ESTIMATION DESIGNS1

A spectral formalism was developed and applied to quantify the sampling errors due to spatial and/or temporal gaps in soil moisture measurements. A design filter was developed to compute the sampling errors for discrete measurements in space and time. This filter has as its advantage a general form applicable to various types of sampling design. The lack of temporal measurements of the two‐dimensional soil moisture field made it difficult to compute the spectra directly from observed records. Therefore, the wave number frequency spectra of soil moisture data derived from stochastic models of rainfall and soil moisture were used. Parameters for both models were estimated using data from the Southern Great Plains Hydrology Experiment (SGP97) and the Oklahoma Mesonet. The estimated sampling error of the spatial average soil moisture measurement by airborne L‐band microwave remote sensing during the SGP97 hydrology experiment is estimated to be 2.4 percent. Under the same climate conditions and soil properties as the SGP97 experiment, equally spaced ground probe networks at intervals of 25 and 50 km are expected to have about 16 percent and 27 percent sampling error, respectively. Satellite designs with temporal gaps of two and three days are expected to have about 6 percent and 9 percent sampling errors, respectively.     

耕作措施对土壤入渗产流影响的实验研究

采用人工模拟降雨和室内分析相结合的方法,研究了黄土区不同耕作措施对降雨入渗的影响。结果表明：①不同耕作管理措施对降雨入渗的影响效用不同,在相同雨强和坡度下,降雨入渗速率表现为：耙耱地〉人工掏挖〉直线坡,在中小雨强和较短滞后情况下,这种情况表现更为显著;②不同耕作管理措施对降雨产流的影响效用不同,在相同雨强和坡度下,产流滞后表现为：耙耱地〉人工掏挖〉直线坡,在中小雨强和较短滞后情况下,这种情况表现更为显著;③根据水量平衡原理,得出了不同耕作管理措施不同坡度下入渗及产流滞后随雨强的变化关系式。上述结果为黄土高原坡耕地水土流失的治理和管理,提供了一定的理论依据。     

Evaluation of rainfall erosivity in Bheta Gad catchment, Kumaun Hills of Uttar Pradesh, central Himalayas

The importance of the size of raindrop in causing soil detachment and splash has long been recognized, although the total energy expended on erosion by splash may be small. The aggressiveness of rainfall or its capacity to cause detachment can be expressed in terms of drop size, rainfall intensity and kinetic energy or momentum. An attempt has been made to determine the rainfall erosivity (EI) of two gauged stations where continuous rainfall recorders were installed, on the basis of rainfall characteristics. Thus, the relationship between average storm EI₃₀ (rainfall erosivity for 30 minutes interval) values and average depths of rainfall could be developed for the Bheta Gad basin of the Gomati River in the Hindu-Kush Himalayas. The analysis has revealed that if factors other than rainfall remain constant, soil splash erosion from cultivated fields is directly proportional to the rainstorm parameter identified as EI.     

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.     

WATER QUALITY MODEL CALIBRATION: A COMPARISON OF INPUT AND OUTPUT ERROR CRITERIA1

ABSTRACT: A common framework for the analysis of water resources systems is the input-parameter-output representation. The system, described by its parameters, is driven by inputs and responds with outputs. To calibrate (estimate the parameters) models of these systems requires data on both inputs and outputs, both of which are subject to random errors. When one is uncertain as to whether the predominant source of error is associated with inputs or outputs, uncertainty also exists as to the correct specification of a calibration criterion. This paper develops and analyzes two alternative least squares criteria for calibrating a numerical water quality model. The first criterion assumes that errors are associated with inputs while the second assumes output errors. Statistical properties of the resulting estimators are examined under conditions of pure input or output error and mixed error conditions from a theoretical perspective and then using simulated results from a series of Monte Carlo experiments.     

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.