Scale effects of hydrostratigraphy and recharge zonation on base flow

Uncertainty regarding spatial variations of model parameters often results in the simplifying assumption that parameters are spatially uniform. However, spatial variability may be important in resource assessment and model calibration. In this paper, a methodology is presented for estimating a critical basin size, above which base flows appear to be relatively less sensitive to the spatial distribution of recharge and hydraulic conductivity, and below which base flows are relatively more sensitive to this spatial variability. Application of the method is illustrated for a watershed that exhibits distinct infiltration patterns and hydrostratigraphic layering. A ground water flow model (MODFLOW) and a parameter estimation code (UCODE) were used to evaluate the influence of recharge zonation and hydrostratigraphic layering on base flow distribution. Optimization after removing spatial recharge variability from the calibrated model altered base flow simulations up to 53% in watersheds smaller than 40 km(2). Merging six hydrostratigraphic units into one unit with average properties increased base flow residuals up to 83% in basins smaller than 50 km(2). Base flow residuals changed <5% in watersheds larger than 40 and 50 km(2) when recharge and hydrostratigraphy were simplified, respectively; thus, the critical basin size for the example area is approximately 40 to 50 km(2). Once identified for an area, a critical basin size could be used to guide the scale of future investigations. By ensuring that parameter discretization needed to capture base flow distribution is commensurate with the scope of the investigation, uncertainty caused by overextending uniform parameterization or by estimating extra parameter values is reduced.     

The spatial and temporal variability of groundwater recharge in a forested basin in northern Wisconsin

Recharge varies spatially and temporally as it depends on a wide variety of factors (e.g. vegetation, precipitation, climate, topography, geology, and soil type), making it one of the most difficult, complex, and uncertain hydrologic parameters to quantify. Despite its inherent variability, groundwater modellers, planners, and policy makers often ignore recharge variability and assume a single average recharge value for an entire watershed. Relatively few attempts have been made to quantify or incorporate spatial and temporal recharge variability into water resource planning or groundwater modelling efforts. In this study, a simple, daily soil–water balance model was developed and used to estimate the spatial and temporal distribution of groundwater recharge of the Trout Lake basin of northern Wisconsin for 1996–2000 as a means to quantify recharge variability. For the 5 years of study, annual recharge varied spatially by as much as 18 cm across the basin; vegetation was the predominant control on this variability. Recharge also varied temporally with a threefold annual difference over the 5‐year period. Intra‐annually, recharge was limited to a few isolated events each year and exhibited a distinct seasonal pattern. The results suggest that ignoring recharge variability may not only be inappropriate, but also, depending on the application, may invalidate model results and predictions for regional and local water budget calculations, water resource management, nutrient cycling, and contaminant transport studies. Recharge is spatially and temporally variable, and should be modelled as such. Copyright © 2009 John Wiley & Sons, Ltd.     

Groundwater recharge: A hydrogeologic thought

The ensuing paper imparts vital information on an important component of the hydrologic cycle: recharge. Water flows through the porous media and forms a definite flow pattern that can be depicted in an elaborate manner through the micro level studies at the small watershed level. The estimation of recharge is indispensable for the groundwater budgeting studies. The advantages and disadvantages of some of the techniques have been reviewed. In the present paper an attempt is made to develop a fundamental understanding of the spatial and temporal distribution of the recharge component and to attract excellent research in the same field.     

Effects of Linking a Soil‐Water‐Balance Model with a Groundwater‐Flow Model

A previously published regional groundwater‐flow model in north‐central Nebraska was sequentially linked with the recently developed soil‐water‐balance (SWB) model to analyze effects to groundwater‐flow model parameters and calibration results. The linked models provided a more detailed spatial and temporal distribution of simulated recharge based on hydrologic processes, improvement of simulated groundwater‐level changes and base flows at specific sites in agricultural areas, and a physically based assessment of the relative magnitude of recharge for grassland, nonirrigated cropland, and irrigated cropland areas. Root‐mean‐squared (RMS) differences between the simulated and estimated or measured target values for the previously published model and linked models were relatively similar and did not improve for all types of calibration targets. However, without any adjustment to the SWB‐generated recharge, the RMS difference between simulated and estimated base‐flow target values for the groundwater‐flow model was slightly smaller than for the previously published model, possibly indicating that the volume of recharge simulated by the SWB code was closer to actual hydrogeologic conditions than the previously published model provided. Groundwater‐level and base‐flow hydrographs showed that temporal patterns of simulated groundwater levels and base flows were more accurate for the linked models than for the previously published model at several sites, particularly in agricultural areas.     

Automating Recession Curve Displacement Recharge Estimation

Recharge estimation is an important and challenging element of groundwater management and resource sustainability. Many recharge estimation methods have been developed with varying data requirements, applicable to different spatial and temporal scales. The variability and inherent uncertainty in recharge estimation motivates the recommended use of multiple methods to estimate and bound regional recharge estimates. Despite the inherent limitations of using daily gauged streamflow, recession curve displacement methods provide a convenient first‐order estimate as part of a multimethod hierarchical approach to estimate watershed‐scale annual recharge. The implementation of recession curve displacement recharge estimation in the United States Geologic Survey (USGS) RORA program relies on the subjective, operator‐specific selection of baseflow recession events to estimate a gauge‐specific recession index. This paper presents a parametric algorithm that objectively automates this tedious, subjective process, parameterizing and automating the implementation of recession curve displacement. Results using the algorithm reproduce regional estimates of groundwater recharge from the USGS Appalachian Valley and Piedmont Regional Aquifer‐System Analysis, with an average absolute error of less than 2%. The algorithm facilitates consistent, completely automated estimation of annual recharge that complements more rigorous data‐intensive techniques for recharge estimation.     

MODIS‐Aided Statewide Net Groundwater‐Recharge Estimation in Nebraska

Monthly evapotranspiration (ET) rates (2000 to 2009) across Nebraska at about 1‐km resolution were obtained by linear transformations of the MODIS (MODerate resolution Imaging Spectroradiometer) daytime surface temperature values with the help of the Priestley–Taylor equation and the complementary relationship of evaporation. For positive values of the mean annual precipitation and ET differences, the mean annual net recharge was found by an additional multiplication of the power‐function‐transformed groundwater vulnerability DRASTIC‐code values. Statewide mean annual net recharge became about 29 mm (i.e., 5% of mean annual precipitation) with the largest recharge rates (in excess of 100 mm/year) found in the eastern Sand Hills and eastern Nebraska. Areas with the largest negative net recharge rates caused by declining groundwater levels due to large‐scale irrigation are found in the south‐western region of the state. Error bounds of the estimated values are within 10% to 15% of the corresponding precipitation rates and the estimated net recharge rates are sensitive to errors in the precipitation and ET values. This study largely confirms earlier base‐flow analysis‐based statewide groundwater recharge estimates when considerations are made for differences in the recharge definitions. The current approach not only provides better spatial resolution than available earlier studies for the region but also quantifies negative net recharge rates that become especially important in numerical modeling of shallow groundwater systems.     

Comparative streamflow characteristics in urbanizing basins in the Portland Metropolitan Area,Oregon, USA

This study investigates changes in streamflow characteristics for urbanizing watersheds in the Portland Metropolitan Area of Oregon for the period from 1951 to 2000. The objective of this study was to assess how mean annual runoff ratio, mean seasonal runoff ratio, annual peak runoff ratio, changes in streamflow in response to storm amount, the fraction of time that the daily mean flow exceeds the annual mean flow, 3‐day recession constants, and dry/wet flow ratio vary among watersheds with different degrees of urban development. There were no statistically significant changes in annual runoff ratio and annual peak runoff ratio for the mixed land‐use watershed (Tualatin River watershed) and the urban watershed (Johnson Creek watershed) during the entire study period. The Tualatin River watershed, where most of the urban development occurred in a lower part of the watershed, showed a statistically significant increase in annual peak runoff ratio during the 1976 and 2000 period. The Upper Tualatin River watershed illustrated a significant decrease in annual peak runoff ratio for the entire study period. With significant differences in seasonal runoff ratio, only Johnson Creek exhibited a significant increase in both wet and dry season runoff ratios. Streamflow during storm events declined rapidly in the urban watershed, with a high 3‐day recession constant. At an event storm scale, streamflow in Fanno Creek, which is the most urbanized watershed, responded quickly to precipitation input. The fraction of time that the daily mean flow exceeded the annual mean flow and dry/wet flow ratio are all lower in Johnson Creek. This suggests a shorter duration of storm runoff and lower baseflow in the urbanized watershed when compared to the mixed land use watershed. The findings of this study demonstrate the importance of spatial and temporal scale, climate variability, and basin physiographic characteristics in detecting the hydrologic effects of urbanization in the Pacific Northwest of the USA. Copyright © 2006 John Wiley & Sons, Ltd.     

Spatial and temporal analysis of groundwater recharge with application to sampling design

A numerical experiment of flow in variably saturated porous media was performed in order to evaluate the spatial and temporal distribution of the groundwater recharge at the phreatic surface for a shallow aquifer as a function of the input rainfall process and soil heterogeneity. The study focused on the groundwater recharge which resulted from the percolation of the excess rainfall for a 90-days period of an actual precipitation record. Groundwater recharge was defined as the water flux across the moving phreatic surface. The observed spatial non-uniformity of the groundwater recharge was caused by soil heterogeneity and is particularly pronounced during the stage of recharge peak (substantial percolation stage). During that stage the recharge is associated with preferential flow paths defined as soil zones of locally higher hydraulic conductivity. For the periods of low percolation intensity the groundwater recharge was exhibiting more uniform spatial characteristics. The temporal distribution of the recharge was found to be a function of the frequency and intensity of the rainfall events. Application of sampling design demonstrates the joint influence of the spatial and temporal recharge variability on the cost-effective monitoring of groundwater potentiometric surfaces.     

Regional estimation of base recharge to ground water using water balance and a base-flow index

Naturally occurring long-term mean annual base recharge to ground water in Nebraska was estimated with the help of a water-balance approach and an objective automated technique for base-flow separation involving minimal parameter-optimization requirements. Base recharge is equal to total recharge minus the amount of evapotranspiration coming directly from ground water. The estimation of evapotranspiration in the water-balance equation avoids the need to specify a contributing drainage area for ground water, which in certain cases may be considerably different from the drainage area for surface runoff. Evapotranspiration was calculated by the WREVAP model at the Solar and Meteorological Surface Observation Network (SAMSON) sites. Long-term mean annual base recharge was derived by determining the product of estimated long-term mean annual runoff (the difference between precipitation and evapotranspiration) and the base-flow index (BFI). The BFI was calculated from discharge data obtained from the U.S. Geological Survey's gauging stations in Nebraska. Mapping was achieved by using geographic information systems (GIS) and geostatistics. This approach is best suited for regional-scale applications. It does not require complex hydrogeologic modeling nor detailed knowledge of soil characteristics, vegetation cover, or land-use practices. Long-term mean annual base recharge rates in excess of 110 mm/year resulted in the extreme eastern part of Nebraska. The western portion of the state expressed rates of only 15 to 20 mm annually, while the Sandhills region of north-central Nebraska was estimated to receive twice as much base recharge (40 to 50 mm/year) as areas south of it.     

Refined chloride mass‐balance method and its application in Saudi Arabia

The rainfall and infiltration elements of the hydrological cycle in arid regions are characterized by temporal and spatial variations that are random and sporadic. Consequently, the chloride concentration in rainfall has a similar behaviour. Despite this, the classical chloride mass balance (CMB) approach only employs arithmetic and weighted averages for recharge estimation. In this paper, the classical CMB method is modified by taking into account some perceived deficiencies in the methodology. The modified CMB method takes into consideration additional statistical parameters, namely variances and the correlation coefficient between variables concerned based on the application of the perturbation method. Strategic aquifer planning in the Kingdom of Saudi Arabia requires a quick method for estimating groundwater recharge in order to determine the temporal management of available water resources. To demonstrate the difference between the classical and the refined CMB methods, both were applied to a representative basin, i.e. Wadi Yalamlam, in the western part of Saudi Arabia. Based on the refined calculations, recharge to groundwater is found to be 11% of the effective annual rainfall. This refined method provides higher recharge rates because it takes into account the actual variability in the variables concerned and can, thus, improve the accuracy of future groundwater recharge estimation studies. Copyright © 2006 John Wiley & Sons, Ltd.     

Estimation of potential evapotranspiration in the mountainous Panama Canal watershed

Spatially distributed hydrometeorological and plant information within the mountainous tropical Panama Canal watershed is used to estimate parameters of the Penman–Monteith evapotranspiration formulation. Hydrometeorological data from a few surface climate stations located at low elevations in the watershed are complemented by (a) typical wet‐ and dry‐season fields of temperature, wind, water vapour and pressure produced by a mesoscale atmospheric model with a 3 × 3 km² spatial and hourly temporal resolution, and (b) leaf area index fields estimated over the watershed during a few years using satellite data with two different spatial and temporal resolutions. The mesoscale model estimates of spatially distributed surface hydrometeorological variables provide the basis for the extrapolation of the surface climate station data to produce input for the Penman–Monteith equation. The satellite information and existing digital spatial databases of land use and land cover form the basis for the estimation of Penman–Monteith spatially distributed parameter values. Spatially distributed 3 × 3 km² potential evapotranspiration estimates are obtained for the 3300 km² Panama Canal watershed. Estimates for Gatun Lake within the watershed are found to reproduce well the monthly and annual lake evaporation obtained from submerged pans. Sensitivity analysis results of potential evapotranspiration estimates with respect to cloud cover, dew formation, leaf area index distribution and mesoscale model estimates of surface climate are presented and discussed. The main conclusion is that even the limited spatially distributed hydrometeorological and plant information used in this study contributes significantly toward explaining the substantial spatial variability of potential evapotranspiration in the watershed. These results also allow the determination of key locations within the watershed where additional surface stations may be profitably placed. Copyright © 2006 John Wiley & Sons, Ltd.     

Transient Recharge Estimability Through Field‐Scale Groundwater Model Calibration

The estimation of recharge through groundwater model calibration is hampered by the nonuniqueness of recharge and aquifer parameter values. It has been shown recently that the estimability of spatially distributed recharge through calibration of steady‐state models for practical situations (i.e., real‐world, field‐scale aquifer settings) is limited by the need for excessive amounts of hydraulic‐parameter and groundwater‐level data. However, the extent to which temporal recharge variability can be informed through transient model calibration, which involves larger water‐level datasets, but requires the additional consideration of storage parameters, is presently unknown for practical situations. In this study, time‐varying recharge estimates, inferred through calibration of a field‐scale highly parameterized groundwater model, are systematically investigated subject to changes in (1) the degree to which hydraulic parameters including hydraulic conductivity (K) and specific yield (S_y) are constrained, (2) the number of water‐level calibration targets, and (3) the temporal resolution (up to monthly time steps) at which recharge is estimated. The analysis involves the use of a synthetic reality (a reference model) based on a groundwater model of Uley South Basin, South Australia. Identifiability statistics are used to evaluate the ability of recharge and hydraulic parameters to be estimated uniquely. Results show that reasonable estimates of monthly recharge (<30% recharge root‐mean‐squared error) require a considerable amount of transient water‐level data, and that the spatial distribution of K is known. Joint estimation of recharge, S_y and K, however, precludes reasonable inference of recharge and hydraulic parameter values. We conclude that the estimation of temporal recharge variability through calibration may be impractical for real‐world settings.     

The effect of random recharge on uniform steady free-surface flow in heterogeneous porous formations

The effect of parametric uncertainty in recharge rate and spatial variability of hydraulic conductivity upon free-surface flow is investigated in a stochastic framework. We examine the three-dimensional free-surface gravitational flow problem for sloped mean uniform flow in a randomly heterogeneous porous medium under the influence of random recharge. We develop analytic solutions for the variance of free-surface position, head, and specific discharge on the free surface. Additionally, we obtain semi-analytic solutions for the statistical moments of head and specific discharge beneath the free-surface. Statistical moments are derived using a first-order approximation and then compared with their parallel in an unbounded medium. The effect of recharge mean and variability on the statistical moments is analyzed. Results can be applied to more complex flows, slowly varying in the mean.     

Assessment of transmission losses and groundwater recharge from runoff events in a wadi under shortage of data on lateral inflow,Negev, Israel

A hydrological–lithostratigraphical model was developed for assessment of transmission losses and groundwater recharge from runoff events in arid water courses where hydrological and meteorological records are incomplete. Water balance equations were established for reaches between hydrometric stations. Because rainfall and tributary flow data are scarce, lateral inflow, which is an essential component of the water balance equation, could not be estimated directly. The solution was obtained by developing a method which includes a hydrological–lithostratigraphical analogy. This is based on the following assumptions: (a) runoff resulting from a given rainfall event is related to the watershed surface lithology; (b) for a given event, the spatial distribution of runoff reflects the distribution of rainfall: and (c) transmission losses are uniquely related to the total inflow to the reach. The latter relationship, called the loss function, and the water balance equation comprise a model which simultaneously assesses lateral inflow and transmission losses for runoff events recorded at the terminal stations. The model was applied to three reaches of the arid Nahal Tsin in Israel. In this case study, the transmission losses were of the same order of magnitude as the flow at the major hydrometric stations. The losses were subdivided into channel moistening, which subsequently evaporates, and deep percolation, which recharges groundwater. For large runoff events, evaporation was substantially smaller than the losses. The mean annual recharge of groundwater from runoff events in the Tsin watershed was 4·1×10⁶ m³, while the mean annual flow volume at the major stations ranged from 0·6 to 1·5×10⁶ m³. Once in 100 years, the annual recharge may be seven times higher than the mean annual value, but the recharge during most years is very small. Copyright © 1999 John Wiley & Sons, Ltd.     

Ground water flow analysis of a mid-Atlantic outer coastal plain watershed, Virginia, U.S.A

Models for ground water flow (MODFLOW) and particle tracking (MODPATH) were used to determine ground water flow patterns, principal ground water discharge and recharge zones, and estimates of ground water travel times in an unconfined ground water system of an outer coastal plain watershed on the Delmarva Peninsula, Virginia. By coupling recharge and discharge zones within the watershed, flowpath analysis can provide a method to locate and implement specific management strategies within a watershed to reduce ground water nitrogen loading to surface water. A monitoring well network was installed in Eyreville Creek watershed, a first-order creek, to determine hydraulic conductivities and spatial and temporal variations in hydraulic heads for use in model calibration. Ground water flow patterns indicated the convergence of flow along the four surface water features of the watershed; primary discharge areas were in the nontidal portions of the watershed. Ground water recharge zones corresponded to the surface water features with minimal development of a regional ground water system. Predicted ground water velocities varied between < 0.01 to 0.24 m/day, with elevated values associated with discharge areas and areas of convergence along surface water features. Some ground water residence times exceeded 100 years, although average residence times ranged between 16 and 21 years; approximately 95% of the ground water resource would reflect land use activities within the last 50 years.     

Throughflow variability on a forested hillslope underlain by compacted glacial till

This study examined the spatial and temporal variability of throughflow in a shallow forest oil in terms of the controls exerted by topography, soil properties and the shape of the water table. Hillslope outflow was measured at nine troughs installed at a road cut, and hydraulic head was measured at 37 piezometers installed to the confining till surface at the base of the soil. At the lowest flows, the distribution of throughflow across the hillslope conformed reasonably to the distribution of upslope contributing area computed from the topography of the confining basal till layer. At higher flows, however, the distribution of throughflow shifted, partly as a result of changes in the orientation in flow tubes as the piezometric surface changed. At high flows, surface topography appeared to provide a reasonable proxy for the shape of the water table. Shunting of water through discrete macropores, however, can apparently overwhelm topographic controls on throughflow, at least for slope widths less than about 10 m. Effective hydraulic conductivities computed at the scale of the troughs varied over two orders of magnitude and displayed no consistent relationship with the thickness of the saturated layer. Estimates at the plot scale (c. 10 m hillslope width) were far less variable, and for higher flows displayed an approximately linear relationship with saturated zone thickness. Copyright © 2000 John Wiley & Sons, Ltd.     

Effects of afforestation on groundwater recharge and water budgets in the western region of Uruguay

Uruguay has stimulated the development of its forest sector since the promulgation of Forest Law N° 15 939 in December of 1987. Nevertheless, the substitution of natural grasslands with forest plantations for industrial use has raised concerns regarding hydrological processes of groundwater recharge and water consumption involving evapotranspiration. The purpose of this study is to assess the effects of this substitution approach on water resources. Input data were collected from two small experimental watersheds of roughly 100–200 hectares located in western Uruguay. The watersheds are characterized by Eucalyptus Globulus ssp. Maidenni and natural grasslands for cattle use. Total rainfall, stream discharge, rainfall redistribution, soil water content and groundwater level data were collected. Groundwater recharge was estimated from water table fluctuations and from groundwater contributions to base flows. Seasonal and annual water budgets were computed from October of 2006 to September of 2014 to evaluate changes in the hydrological processes. The data show a decrease in annual specific discharge of roughly 17% for mean hydrological years and no conclusive effects on annual groundwater recharge in the forested watershed relative to the reference pasture watershed. Reduced annual specific discharge is equivalent to the mean annual interception. The computed actual annual evapotranspiration is consistent with international catchment measurements. Reduction rates vary seasonally and according to accumulated rainfall and its temporary distribution. The degree of specific discharge decline is particularly high for drier autumns and winters (32 to 28%) when the corresponding rainfall varies from 275 to 400 mm. These results are of relevance for water resources management efforts, as water uses downstream can be affected. These findings, based on a study period dominated by anomalous wet springs and summers and by dry autumns and winters, oppose earlier results based on 34 years of rainfall and discharge data drawn from Uruguayan large basins. Copyright © 2016 John Wiley & Sons, Ltd.     

EFFECT OF DIGITAL ELEVATION MODEL RESOLUTION ON EMPIRICAL ESTIMATION OF SOIL LOSS AND SEDIMENT TRANSPORT WITH GIS

The horizontal accuracy of topographic data represented by digital elevation model (DEM) resolution brings about uncertainties in landscape process modeling with raster GIS. This paper presents a study on the effect of topographic variability on cell-based empirical estimation of soil loss and sediment transport. An original DEM of 10m resolution for a case watershed was re-sampled to three realizations of higher grid sizes for a comparative examination. Equations based on the USLE are applied to the watershed to calculate soil loss from each cell and total sediment transport to streams. The study found that the calculated total soil loss from the watershed decreases with the increasing DEM resolution with a linear correlation as spatial variability is reduced by cell aggregation. The USLE topographic factors (LS) extracted from applied DEMs represent spatial variability, and determine the estimations as shown in the modeling results. The commonly used USGS 30m DEM appears to be able to reflect essential spatial variability and suitable for the empirical estimation. The appropriateness of a DEM resolution is dependent upon specific landscape characteristics, applied model and its parameterization. This work attempts to provide a general framework for the research in the DEM-based empirical modeling.     

Decoupling the effects of deforestation and climate variability in the Tapajós river basin in the Brazilian Amazon

Tropical river basins are experiencing major hydrological alterations as a result of climate variability and deforestation. These drivers of flow changes are often difficult to isolate in large basins based on either observations or experiments; however, combining these methods with numerical models can help identify the contribution of climate and deforestation to hydrological alterations. This paper presents a study carried out in the Tapaj?s River (Brazil), a 477,000 km² basin in South‐eastern Amazonia, in which we analysed the role of annual land cover change on daily river flows. Analysis of observed spatial and temporal trends in rainfall, forest cover, and river flow metrics for 1976 to 2008 indicates a significant shortening of the wet season and reduction in river flows through most of the basin despite no significant trend in annual precipitation. Coincident with seasonal trends over the past 4 decades, over 35% of the original forest (140,000 out of 400,000 km²) was cleared. In order to determine the effects of land clearing and rainfall variability to trends in river flows, we conducted hindcast simulations with ED2 + R, a terrestrial biosphere model incorporating fine scale ecosystem heterogeneity arising from annual land‐use change and linked to a flow routing scheme. The simulations indicated basin‐wide increases in dry season flows caused by land cover transitions beginning in the early 1990s when forest cover dropped to 80% of its original extent. Simulations of historical potential vegetation in the absence of land cover transitions indicate that reduction in rainfall during the dry season (mean of ?9 mm per month) would have had an opposite and larger magnitude effect than deforestation (maximum of +4 mm/month), leading to the overall net negative trend in river flows. In light of the expected increase in future climate variability and water infrastructure development in the Amazon and other tropical basins, this study presents an approach for analysing how multiple drivers of change are altering regional hydrology and water resources management.     

Estimating groundwater recharge in lowland watersheds

Little attention has been given to the role of groundwater in the hydrological cycle of lowland watersheds. Our objective in this study was to estimate total recharge to groundwater by analysing water table response to storm events and the rate at which water was transferred into the shallow aquifer. This was conducted at three sites in a rural watershed in the lower Atlantic coastal plain near Charleston, South Carolina, USA. A novel version of the water table fluctuation method was used to estimate total recharge to the shallow aquifer by comparing hourly data of water table position following storm events and measuring water table recession behavior, rather than subjective graphical analysis methods. Also, shallow aquifer recharge rates (vertical fluxes) were estimated using Darcy's Law by comparing static water levels in a water table well and in a shallow piezometer during dry periods. The total annual recharge estimated ranged from 107 ± 39 mm·yr^–1 (5–10% of annual precipitation) at a poorly drained topographic low area to 1140 ± 230 mm·yr^–1 (62–94% of annual precipitation) for a moderately well‐drained upland site. The average aquifer recharge rate was 114 ± 60 mm·yr^–1, which is similar to previous estimations of base flow for the ephemeral third‐order streams in this watershed. The difference in the two methods may have been caused by processes not accounted for in the Darcy flux method, soil moisture deficits, and average evapotranspiration demand, which is about 1100 mm·yr^–1 for this region. Although other factors also can affect partitioning of recharge, an integrated approach to inspecting easily gathered groundwater data can provide information on an often neglected aspect of water budget estimation. We also discuss the effects of land use change on recharge reduction, given a typical development scenario for the region. Copyright © 2011 John Wiley & Sons, Ltd.