Recharge is a critical issue for water management. Recharge assessment and the factors affecting recharge are of scientific and practical importance. The purpose of this study was to develop a daily recharge assessment model (DREAM) on the basis of a water balance principle with input from conventional and generally available precipitation and evaporation data and demonstrate the application of this model to recharge estimation in the Western Mountain Aquifer (WMA) in Israel. The WMA (area 13,000 km2)isa karst aquifer that supplies 360–400 Mm3 yr−1 of freshwater, which constitutes 20% of Israel's freshwater and is highly vulnerable to climate variability and change. DREAM was linked to a groundwater flow model (FEFLOW) to simulate monthly hydraulic heads and spring flows. The models were calibrated for 1987–2002 and validated for 2003– 2007, yielding high agreement between calculated and measured values (R2 = 0.95; relative root‐mean‐square error = 4.8%; relative bias = 1.04). DREAM allows insights into the effect of intra‐annual precipitation distribution factors on recharge. Although annual precipitation amount explains ∼70% of the variability in simulated recharge, analyses with DREAM indicate that the rainy season length is an important factor controlling recharge. Years with similar annual precipitation produce different recharge values as a result of temporal distribution throughout the rainy season. An experiment with a synthetic data set exhibits similar results, explaining ∼90% of the recharge variability. DREAM represents significant improvement over previous recharge estimation techniques in this region by providing near‐real‐time recharge estimates that can be used to predict the impact of climate variability on groundwater resources at high temporal and spatial resolution.
This paper summarises innovative research into the assessment of long-term groundwater recharge from flood events in dryland environments of the Kuiseb (Namibia) and the Buffels (South Africa) rivers. The integrated water resource management (IWRM) policies and institutions affecting the exploitation of groundwater resources in each of these developing countries are compared. The relatively large alluvial aquifer of the Kuiseb River (similar to 240 Mm(3)) is recharged from irregular floods originating in the upper catchment. Reported abstraction of 4.6 Mm(3) per year is primarily consumed in the town of Walvis Bay, although the groundwater decay (pumping and natural losses along the period 1983-2005) was estimated in 14.8 Mm(3) per year. Recharge is variable, occurring in 11 out of 13 years in the middle Kuiseb River, but only in 11 out of 28 years in the middle-lower reaches. In contrast, the Buffels River has relatively minor alluvial aquifers (similar to 11 Mm(3)) and recharge sources derive from both lateral subsurface flow and floodwater infiltration, the latter limited to a recharge maximum of 1.3 Mm(3) during floods occurring once every four years. Current abstractions to supply the adjacent rural population and a few small-scale, irrigated commercial farms are 0.15 Mm(3) yr (-aEuro parts per thousand 1), well within the long-term sustainable yield estimated to be 0.7 Mm(3) yr (-aEuro parts per thousand 1). Since independence in 1990, Namibia`s water resource management approach has focussed on ephemeral river basin management of which the Kuiseb Basin Management Committee (KBMC) is a model. Here, some water points are managed independently by rural communities through committees while the national bulk water supplier provides for Walvis Bay Municipality from the lower aquifers. This provides a sense of local ownership through local participation between government, NGOs and CBOs (community-based organisations) in the planning and implementation of IWRM. Despite the potential for water resource development in the lower Buffels River, the scope for implementing IWRM is limited not only by the small aquifer size, but also because basin management in South Africa is considered only in the context of perennial rivers. Since 2001, water service delivery in the Buffels River catchment has become the responsibility of two newly created local municipalities. As municipal government gains experience, skills and capacity, its ability to respond to local needs related to water service delivery will be accomplished through local participation in the design and implementation of annual `integrated development plans`. These two case studies demonstrate that a variety of IWRM strategies in the drylands of developing countries are appropriate depending on scales of governance, evolving policy frameworks, scales of need and limitations inherent in the hydrological processes of groundwater resources.
Flood water infiltrates ephemeral channels, recharging local and regional aquifers, and it is the main water source in hyperarid regions. Quantitative estimations of these resources are limited by the scarcity of data from such regions. The floods of the Kuiseb River in the Namib Desert have been monitored for 46 years, providing a unique data set of flow hydrographs from one of the world's hyperarid regions. The study objectives were to: (1) subject the records to quality control; (2) model flood routing and transmission losses; and (3) study the relationships between flood characteristics, river characteristics and recharge into the aquifers. After rigorous quality-testing of the original gauge-station data, a flood-routing model based on kinematic flow with components accounting for channel-bed infiltration was constructed and applied to the data. A simplified module added to this routing model estimates aquifer recharge from the infiltrating flood water. Most of the model parameters were obtained from field surveys and GIS analyses. Two of the model parameters-Manning's roughness coefficient and the constant infiltration rate-were calibrated based on the high-quality measured flow data set, providing values of 0.025 and 8.5 mm/h, respectively. This infiltration rate is in agreement with that estimated from extensive direct TDR-based moisture measurements in the vadose zone under the Kuiseb River channel, and is low relative to those reported for other sites. The model was later verified with additional flood data and observed groundwater levels in boreholes. Sensitivity analysis showed the important role of large and medium floods in aquifer recharge. To generalize from the studied river to other streams with diverse conditions, we demonstrate that with increasing in infiltration rate, channel length or active channel width, the relative contribution of high-magnitude floods to recharge also increases, whereas medium and small floods contribute less, often not reaching the downstream parts of the arid ephemeral river at all. For example, more than three-quarters of the floods reaching the downstream Kuiseb River (with an infiltration rate of 8.5 mm/h) would not have reached similar distances in rivers with all other properties similar but with infiltration rates of 50 mm/h. The recharge volume in the downstream segment in the case of higher infiltration is mainly contributed by floods with magnitude ???93rd percentile, compared to floods in the 63rd percentile at an infiltration rate of 8.5 mm/h. ?? 2009 Elsevier B.V. All rights reserved.
Detailed hydrologic models require high-resolution spatial and temporal data. This study aims at improving the spatial interpolation of daily precipitation for hydrologic models. Different parameterizations of (1) inverse distance weighted (IDW) interpolation and (2) A local weighted regression (LWR) method in which elevation is the explanatory variable and distance, elevation difference and aspect difference are weighting factors, were tested at a hilly setting in the eastern Mediterranean, using 16 years of daily data. The preferred IDW interpolation was better than the preferred LWR scheme in 27 out of 31 validation gauges (VGs) according to a criteria aimed at minimizing the absolute bias and the mean absolute error (MAE) of estimations. The choice of the IDW exponent was found to be more important than the choice of whether or not to use elevation as explanatory data in most cases. The rank of preferred interpolators in a specific VG was found to be a stable local characteristic if a sufficient number of rainy days are averaged. A spatial pattern of the preferred IDW exponents was revealed. Large exponents (3) were more effective closer to the coast line whereas small exponents (1) were more effective closer to the mountain crest. This spatial variability is consistent with previous studies that showed smaller correlation distances of daily precipitation closer to the Mediterranean coast than at the hills, attributed mainly to relatively warm sea-surface temperature resulting in more cellular convection coastward. These results suggest that spatially variable, physically based parameterization of the distance weighting function can improve the spatial interpolation of daily precipitation
Quantitatively estimating rainfall-runoff relations in extremely arid regions is a challenging task, mainly because of lack of in situ data. For the past 40 years, rain and floods have been monitored in the Nahal Yael catchment (0.5 km2) in southern Israel, providing a unique data set of runoff hydrographs and rainfall in a hyper-arid region. Here we present an exploratory study focusing on rainfall-runoff modeling issues for a small (0.05 km2) sub-catchment of Nahal Yael. The event-based model includes the computation of rainfall excess, hillslope and channel routing. Two model parameters of the infiltration process were found by calibration. A resampling methodology of calibration group composition is suggested to derive optimal model parameters and their uncertainty range. Log-based objective functions were found to be more robust and less sensitive than non-log functions to calibration group composition. The fit achieved between observed and computed runoff hydrographs for the calibration and validation events is considered good relative to other modeling studies in arid and semi-arid regions. The study indicates that, under the calibration scheme used, a lumped model performs better than a model representing the catchment division into three sub-catchments. In addition, the use of rain data from several gauges improves runoff prediction as compared to input from a single gauge. It was found that rainfall uncertainty dominates uncertainties in runoff prediction while parameter uncertainties have only a minor effect. ?? 2009 Elsevier B.V. All rights reserved.
Flash-flood warning models can save lives and protect various kinds of infrastructure. In dry climate regions, rainfall is highly variable and can be of high-intensity. Since rain gauge networks in such areas are sparse, rainfall information derived from weather radar systems can provide useful input for flash- flood models. This paper presents a flash-flood warning model which utilizes radar rainfall data and applies it to two catchments that drain into the dry Dead Sea region. Radar-based quantitative precipita- tion estimates (QPEs) were derived using a rain gauge adjustment approach, either on a daily basis (allowing the adjustment factor to change over time, assuming available real-time gauge data) or using a constant factor value (derived from rain gauge data) over the entire period of the analysis. The QPEs served as input for a continuous hydrological model that represents the main hydrological processes in the region, namely infiltration, flow routing and transmission losses. The infiltration function is applied in a distributed mode while the routing and transmission loss functions are applied in a lumped mode. Model parameters were found by calibration based on the 5 years of data for one of the catchments. Val- idation was performed for a subsequent 5-year period for the same catchment and then for an entire 10- year record for the second catchment. The probability of detection and false alarm rates for the validation cases were reasonable. Probabilistic flash-flood prediction is presented applying Monte Carlo simulations with an uncertainty range for the QPEs and model parameters. With low probability thresholds, one can maintain more than 70% detection with no more than 30% false alarms. The study demonstrates that a flash-flood warning model is feasible for catchments in the area studied.
viously applied in the Alps of Europe. Adjustment coefficients have been derived for 28 rainfall periods using 59 independent gauges of a quality-checked training data set. The validation was based on an independent data set composed of gauges located in eleven 20 ? 20 km2 validation areas, which are representative of different climate, topography and radar distance conditions. The WR and WMR methods were found preferable with a slight better performance of the latter. Furthermore, a novel approach has been adopted in this study, whereby radar estimates are considered useable if they provide information that is better than gauge-only estimates. The latter was derived by spatial interpolation of the gauges belonging to the training data set. Note that these gauges are outside the validation areas. As for the radar-adjusted estimates, gauge-derived estimates were assessed against gauge data in the validation areas. It was found that radar-based estimates are better for the validation areas at the dry climate regime. At distances larger than 100 km, the radar underestimation becomes too large in the two northern validation areas, while in the southern one radar data are still better than gauge interpolation. It is concluded that in ungauged areas of Israel it is preferable to use WMR-adjusted (or alternatively, simply WR-adjusted) radar echoes rather than the standard bulk adjustment method and for dry ungauged areas it is preferable over the conventional gauge-interpolated values derived from point measurements, which are outside the areas themselves. The WR and WMR adjustment methods provide useful rain depth estimates for rainfall periods for the examined areas but within the limitation stated above.
Analysis of extreme hydrometeorological events is important for characterizing and better understanding the meteorological conditions that can generate severe rainstorms and the consequent catastrophic flooding. According to several studies (e.g., Alpert et al., 2004; Wittenberg et al., 2007), the occurrence of such extreme events is increasing over the eastern Mediterranean although total rain amounts are generally decreasing. The current study presents an analysis of an extreme event utilizing different methodologies: (a) synoptic maps and high resolution satellite imagery for atmospheric condition analysis; (b) rainfall analysis by rain gauges data; (c) meteorological radar rainfall calibration and analysis; (d) field measurements for estimating maximum peak discharges; and, (e) high resolution aerial photographs together with field surveying for quantifying the geomorphic impacts. The unusual storm occurred over Israel between 30 March and 2 April, 2006. Heavy rainfall produced more than 100mm in some locations in only few hours and more than 200mm in the major core area. Extreme rain intensities with recurrence intervals of more than 100 years were found for durations of 1 h and more as well as for the daily rain depth values. In the most severely affected area,Wadi Ara, extreme flash floods caused damages and casualties. Specific peak discharges were as high as 10–30m3/s/km2 for catchments of the size of 1–10 km2, values larger than any recorded floods in similar climatic regions in Israel.
Weather radar data contain detailed information about the spatial structures of rain fields previously unavail- able from conventional rain gauge networks. This information is of major importance for enhancing our understanding of precipitation and hydrometeorological systems. This study focuses on spatial features of convective rain cells in southern Israel where the climate ranges fromMediterranean to hyper-arid. Extensive data bases from two study areas covered by radar systems were analyzed. Rain cell features were extracted such as center location, area, maximal rain intensity, spatial integral of rain intensity, major radius length, minor radius length, ellipticity, and orientation. Rain cells in the two study areas were compared in terms of feature distributions and the functional relationships between cell area and cell magnitude, represented by maximal rain intensity and spatial integral of rain intensity. Analytical distribution functions were fitted to the empirical distributions and the log-normal function was found to fit well the distributions of cell area, maximal rain intensity and major and minor radius lengths. The normal distribution fits well ellipticity em- pirical distribution, and orientation distribution was well-represented by the normal or uniform distribution functions. The effect of distance fromtheMediterranean coastline on cell features was assessed. Amaximum of cell rain intensity at the coastline and maximum cell density 15 km inland from the coastline were found. In addition, a gradual change of cell orientation was observed with a northwest-southeast orientation 30 km from the coastline at the Mediterranean Sea and to almost a west-east orientation 30 km from the coastline inland
Weather radar systems provide detailed information on spatial rainfall patterns known to play a significant role in runoff generation processes. In the current study, we present an innovative approach to exploit spatial rainfall information of air mass thunderstorms and link it with a watershed hydrological model. Observed radar data are decomposed into sets of rain cells conceptualized as circular Gaussian elements and the associated rain cell parameters, namely, location, maximal intensity and decay factor, are input into a hydrological model. Rain cells were retrieved from radar data for several thunderstorms over southern Arizona. Spatial characteristics of the resulting rain fields were evaluated using data from a dense rain gauge network. For an extreme case study in a semi-arid watershed, rain cells were derived and fed as input into a hydrological model to compute runoff response. A major factor in this event was found to be a single intense rain cell (out of the five cells decomposed from the storm). The path of this cell near watershed tributaries and toward the outlet enhanced generation of high flow. Furthermore, sensitivity analysis to cell characteristics indicated that peak discharge could be a factor of two higher if the cell was initiated just a few kilometers aside.
A spatial rainfall model was applied to radar data of air mass thunderstorms to yield a rainstorm representation as a set of convective rain cells. The modeled rainfall was used as input into hydrological model, instead of the standard radar-grid data. This approach allows a comprehensive linkage between runoff responses and rainfall structures
Radar-based estimates of rainfall rates and accumulations are one of the principal tools used by the National Weather Service (NWS) to identify areas of extreme precipitation that could lead to flooding. Radar-based rainfall estimates have been compared to gauge observations for 13 convective storm events over a densely instrumented, experimental watershed to derive an accurate reflectivity–rainfall rate (i.e., Z–R) relationship for these events. The resultant Z–R relationship, which is much different than the NWS operational Z–R, has been examined for a separate, independent event that occurred over a different location. For all events studied, the NWS operational Z–R significantly overestimates rainfall compared to gauge measurements. The gauge data from the experimental network, the NWS operational rain estimates, and the improved estimates resulting from this study have been input into a hydrologic model to “predict” watershed runoff for an intense event. Rainfall data from the gauges and from the derived Z–R relation produce predictions in relatively good agreement with observed streamflows. The NWS Z–R estimates lead to predicted peak discharge rates that are more than twice as large as the observed discharges. These results were consistent over a relatively wide range of subwatershed areas (4–148 km2). The experimentally derived Z–R relationship may provide more accurate radar estimates for convective storms over the southwest United States than does the operational convective Z–R used by the NWS. These initial results suggest that the generic NWS Z–R relation, used nationally for convective storms, might be substantially improved for regional application.
Intense rainstorms cause debris flows on escarpments in hyperarid environments. In contrast with more temperate environments, there have been no direct observations on rainfall intensities and durations required for initiating debris flows in hyperarid environments. Here, we report rainfall volume and intensities, acquired by gauge and radar measurements, for two successive storms along the hyperarid (\textless50 mm/yr) western escarpment of the Dead Sea basin. These rainfall data were analyzed in conjunction with detailed mapping of debris flows that occurred during these storms to determine values of rainfall intensity and duration required to generate debris flows on the Dead Sea western escarpment. The first of the two analyzed storms occurred on 2 November 1995. During this storm, two convective cells rained sequentially within a 5 It period at the lower reaches of the Nahal David and the Nahal 'Arugot that dissects the western escarpment of the Dead Sea, Israel. This storm triggered debris flows in 38 small (\textless3 km(2)) and high-gradient drainage basins along the escarpment. Total rainfall volume and spatial distribution were determined by 10 cumulative rain gauges that were also used to calibrate rainfall-intensity distributions from radar data. For this storm, region, and landscape, rainfall intensities exceeding 30 mm/h for a duration of I h were required to initiate debris flows. A second storm in the same area on 1718 October 1997 allowed the evaluation of the results determined from the 1995 storm. In this second, more regional storm, maximum rainfall intensities were 19-27 mm/h for a duration of 45 min. These values, lower than the 30 mm/h minimal threshold defined in the previous storm, are consistent with the occurrence of only three debris flows. The small number of debris flows resulted from the concentration of the highest intensities of rainfall on the desert plateau and not directly on top of the canyon walls. Most first- to third-order basins draining the Dead Sea escarpment contain evidence of zero to three late Holocene (\textless3000 yr) debris flows. From analysis of the two storms, we propose that most of these debris flows were triggered by storms similar to the 2 November 1995 event in which localized convective cells had rainfall intensities of \textgreater30 mm/h and durations of at least 1 h. The small number of debris flows that has occurred during the late Holocene indicates that such events are rare at the scale of individual drainage basins.
Global precipitation is monitored from a variety of platforms including spaceborne, ground-, and ocean-based platforms. Intercomparisons of these observations are crucial to validating the measurements and providing confidence for each measurement technique. Probability distribution functions of rain rates are used to compare satellite and ground-based radar observations. A preferred adjustment technique for improving rain rate distribution estimates is identified using measurements from ground-based radar and rain gauges within the coverage area of the radar. The underwater measurement of rainfall shows similarities to radar measurements, but with intermediate spatial resolution and high temporal resolution. Reconciling these different measurement techniques provides understanding and confidence for all of the methods.
Meteorological radar is a remote sensing system that provides rainfall estimations at high spatial and temporal resolutions. The radar-based rainfall intensities (R) are calculated from the observed radar reflectivities (Z). Often, rain gauge rainfall observations are used in combination with the radar data to find the optimal parameters in the Z–R transformation equation. The scale dependency of the power-law Z–R parameters when estimated from radar reflectivity and rain gauge intensity data is explored herein. The multiplicative (a) and exponent (b) parameters are said to be “scale dependent” if applying the observed and calculated rainfall intensities to objective function at different scale results in different “optimal” parameters. Radar and gauge data were analyzed from convective storms over a midsize, semiarid, and well-equipped watershed. Using the root-mean-square difference (rmsd) objective function, a significant scale dependency was observed. Increased time- and space scales resulted in a considerable increase of the a parameter and decrease of the b parameter. Two sources of uncertainties related to scale dependency were examined: 1) observational uncertainties, which were studied both experimentally and with simplified models that allow representation of observation errors; and 2) model uncertainties. It was found that observational errors are mainly (but not only) associated with positive bias of the b parameter that is reduced with integration, at least for small scales. Model errors also result in scale dependency, but the trend is less systematic, as in the case of observational errors. It is concluded that identification of optimal scale for Z–R relationship determination requires further knowledge of reflectivity and rain-intensity error structure.
A new characteristic timescale of a catchment is presented, the response timescale (RTS). It is a range of averaging time intervals which, when applied to catchment rainfall, yield smoothed time series that best approximate that of the resultant runoff. In determining the RTS, nothing is assumed about the nature of the rainfall-runoff transformation. In addition, this new measure is shown to be robust against measurement errors. An objective, automatic, observations-based algorithm is described that introduces the concept of peaks density for the estimation of RTS. Estimation is exemplified for single and multiple rainfall-runoff events through application to small catchments in Panama and Israel. In all cases, relatively stable values of response timescale are obtained. It is concluded that at least for the case studies, the response timescale is an intrinsic characteristic of the catchment and it is generally expected to be different from the catchment lag time and time of concentration. INDEX
The transformation of rainfall into runoff at a basin outlet is the combined effect of many hydrological processes, which occur at a wide range of spatial and temporal scales. However, determining the scale of the combined hydrological response of the basin is still problematic and concepts for its definition are yet to be identified. In this paper high-resolution meteorological radar data are used for the determination of a characteristic temporal scale for the hydrological response of the basin - the 'response time scale' (Ts*). Ts* is defined as the time scale at which the pattern of the time-averaged radar rainfall hietograph is most similar to the pattern of the measured outlet runoff hydrograph. The existence of such similarity at a relatively stable time scale for a specific basin indicates that it is an intrinsic property of the basin and is related to its hydrological response. The identification of the response time scale is carried out by analysis of observations only, without assuming a specific rainfall-runoff model. Ts* is examined in four small basins (10-100 km2) in Israel. The spatial scale is assumed as the entire basin. For all analyzed basins a stable response time scale is identified. Relatively short time scales are found for the urban and arid basins (15-30 min), while for the rural basins longer time scale are identified (90-180 min). The issues of relationship between the response time scale and basin properties and modeling at the response time scale have yet to be determined. ?? 2001 Elsevier Science B.V. All rights reserved.