Methods of groundwater recharge estimation in eastern Mediterranean—a water balance model application in Greece,Cyprus and Jordan

Groundwater recharge studies in semi‐arid areas are fundamental because groundwater is often the only water resource of importance. This paper describes the water balance method of groundwater recharge estimation in three different hydro‐climatic environments in eastern Mediterranean, in northwest Greece (Aliakmonas basin/Koromilia basin), in Cyprus (Kouris basin and Larnaka area) and in Jordan (northern part of Jordan). For the Aliakmonas basin, groundwater recharge was calculated for different sub‐catchments. For the Upper Aliakmonas basin (Koromilia basin), a watershed‐distributed model was developed and recharge maps were generated on a daily basis. The mean annual recharge varied between 50 and 75 mm/year (mean annual rainfall 800 mm/year). In Cyprus, the mean groundwater recharge estimates yielded 70 mm/year in the Kouris basin. In the Larnaka area, groundwater recharge ranged from 30 mm/year (lowland) to 200 mm/year (mountains). In Jordan, the results indicated recharge rates ranging from 80 mm/year for very permeable karstified surfaces in the upper part of the Salt basin, where rainfall reaches 500 mm/year to less than 10 mm/year and to only about 1 mm/year in the southernmost part of the basin. For the north part of Jordan, a watershed‐distributed model was developed and recharge maps were generated. This water balance model was used for groundwater recharge estimations in many regions with different climatic conditions and has provided reliable results. It has turned out to be an important tool for the management of the limited natural water resources, which require a detailed understanding of regional hydro(geo)logical processes. Copyright © 2007 John Wiley & Sons, Ltd.     

Rate and historical change of direct recharge from precipitation constrained by unsaturated zone profiles of chloride and oxygen‐18 in dry river bed of North China Plain

Quantifying of direct recharge derived from precipitation is crucial for assessing sustainability of well‐irrigated agriculture. In the North China Plain, the land use is dominated by groundwater‐irrigated farmland where the direct recharge derived from precipitation and irrigation. To characterize the mean rate and historical variance of direct recharge derived from precipitation, unsaturated zone profiles of chloride and δ¹⁸O in the dry river bed of the Beiyishui River were employed. The results show that archival time scale of the profile covers the duration from 1980 to 2002 (corresponding to depths from 5 to 2 m) which is indicated by matching the δ¹⁸O peaks in the isotope profile with the aridity indexes gained by instrumental records of annual precipitation and annual potential evaporation. Using the chloride mass balance method, the mean rate of the direct recharge corresponding to the archival time scale is estimated to be 3·8 ± 0·8 mm year^?1, which accounts for about 0·7% of the long‐term average annual precipitation. Further, the direct recharge rates vary from 2·1 to 6·8 mm year^?1 since 1980. Despite the subhumid climate, the estimate of recharge rates is in line with other findings in semiarid regions. The low rate of direct recharge is considered as a result of the relative dry climate in recent decades. In dry river bed, unsaturated zone profiles of chloride and δ¹⁸O combined with instrumental records could offer valuable information about the direct recharge derived from precipitation during droughts. Copyright © 2011 John Wiley & Sons, Ltd.     

Groundwater recharge assessment in a rural,arid, mid-mountain basin in North-Central Chile

ABSTRACT

A basic component of any hydrogeological study is the magnitude and temporal variation of groundwater recharge. This can be difficult to assess accurately, particularly in arid and semi-arid rainfed mid-mountain zones, as is the situation in the rural, low population density zones of North-Central Chile. In this study, recharge in the Punitaqui Basin, North-Central Chile, was characterized, contrasting the results of two methods: a modified Thornthwaite-Mather (MTM) and discharge recession analysis (DRA). We found a recharge rate of between 1 and 4% of average annual precipitation. Average recharge estimated by the MTM method is consistently higher than that estimated by DRA. Also, DRA tends to smooth the recharge values, resulting in a lower inter-annual variation coefficient. Both methods identified a threshold value of total annual precipitation, above which recharge can be expected to occur, of the order of 180 mm year^?1, consistent with values reported in similar areas.     

Soil profile evolution following land‐use change: implications for groundwater quantity and quality

Soil and vadose zone profiles are used as an archive of changes in groundwater recharge and water quality following changes in land use in an area of the Loess Plateau of China. A typical rain‐fed loess‐terrace agriculture region in Hequan, Guyuan, is taken as an example, and multiple tracers (chloride mass balance, stable isotopes, tritium and water chemistry) are used to examine groundwater recharge mechanisms and to evaluate soil water chloride as an archive for recharge rate and water quality. Results show that groundwater recharge beneath natural uncultivated grassland, used as a baseline, is about 94–100 mm year^?1 and that the time it takes for annual precipitation to reach water table through the thick unsaturated zone is from decades to hundreds of years (tritium free). This recharge rate is 2–3 orders of magnitude more than in the other semiarid areas with similar annual rainfall but with deep‐rooted vegetation and relatively high temperature. Most of the water that eventually becomes recharge originally infiltrated in the summer months. The conversion from native grassland to winter wheat has reduced groundwater recharge by 42–50% (50–55 mm year^?1 for recharge), and the conversion from winter wheat to alfalfa resulted in a significant chloride accumulation in the upper soil zone, which terminated deep drainage. The paper also evaluates the time lag between potential recharge and actual recharge to aquifer and between increase in solute concentration in soil moisture and that in the aquifer following land‐use change due to the deep unsaturated zone. Copyright © 2012 John Wiley & Sons, Ltd.     

Application of chloride profile and water balance methods in estimating groundwater recharge in Luanjing Irrigation Area,Inner Mongolia / Application des méthodes du profil de chlorure et du bilan hydrique pour l'estimation de la recharge hydrogéologique dans le Périmètre Irrigué de Luanjing,Mongolie Intérieure

Abstract

A study was carried out to investigate the use of the chloride profile method in conjunction with the water balance method to estimate the annual groundwater recharge in both natural and irrigation sites in Luanjing Irrigation Area, Inner Mongolia. Groundwater recharge from precipitation, estimated by the chloride profile method, is less than 0.1 mm year^?1 which accounts for just 0.06% of the long-term average annual rainfall, indicating that rainfall presently plays a minor role in the groundwater recharge. It appears that recharge events only occurred after heavy rain or sustained rainfall events. In the cropped area, the chloride profile method indicated that the average annual recharge is 268 mm year^?1 with an infiltration rate of 32.5%, which is reasonably consistent with the 33.1% obtained by the water balance method in 2007. The study shows that about one third of that water is discharged back to the groundwater.     

A combined methodology for estimating the potential natural aquifer recharge in an arid environment

ABSTRACT

An innovative methodology that combines an indirect physiography-based method for determining the runoff coefficient at a sub-basin scale and a water balance model applied on a daily time scale was developed to calculate the natural groundwater recharge in three watersheds within the Oum Zessar arid area, Tunisia. The effective infiltration was calculated as part of the water surplus by considering the average available water content (AWC) of soil and an average runoff coefficient for each sub-basin. The model indicates that the sub-basins covered mainly by the “artificial” soils of tabias and jessour, characterized by average AWC values greater than 150 mm, did not contribute to natural groundwater recharge over the 10-year period (2003–2012) considered. The estimated volume for the Triassic aquifer amounted to about 4.5 hm³ year^?1, which is consistent with previous studies. For the Jurassic and Cretaceous aquifers, the estimated volumes amounted to about 200 dm³ year^?1.     

Hydrological and erosional response to natural rainfall in a semi‐arid area of south‐east Spain

A better knowledge of soil erosion by water is essential for planning effective soil and water conservation practices in semi‐arid Mediterranean environments. The special climatic and hydrological characteristics of these areas, however, make accurate soil loss predictions difficult, particularly in the absence of minimal data. Two zero‐order experimental microcatchments (328–759 m²), representative of an extensive semi‐arid watershed with a high potential erosion risk in the south‐east of Spain, were selected and monitored for 3 years (1991–93) in order to provide information on the hydrological and erosional response. A pluviogram and hydrograph recorded data at 1‐min intervals during each storm, after which the soil loss was collected and the particle size of the sediment was analysed. Runoff coefficients of about 9% and soil losses of between 84·83 and 298·9 g m^?2 year^?1 were observed in the area. Rapid response times (geometric mean values lower than 2 h) and low runoff thresholds (mean values between 3·5 to 5·9 mm) were the norm in the experimental areas. A rain intensity of over 15 mm h^?1 was considered as ‘erosive rainfall’ in these areas because of the total soil loss and the transport capacity of the overland flow. Differences in pore‐size distribution explained the different hydrological responses observed between areas. The erosional response was more complex and basically seemed to be determined by soil aggregate stability and topographical properties. A greater proportion of finer particles in the eroded material than in the soil matrix indicated selective erosion and the transport of finer material. Copyright © 2001 John Wiley & Sons, Ltd.     

A recharge model for high altitude,arid, Andean aquifers

Evidence for groundwater recharge in arid zones is mounting, despite early ideas that recharge was unlikely where evaporation greatly exceeded precipitation. The mechanisms and magnitude of groundwater recharge in the Andes and Atacama Desert are not well known but the subject of current research. Diffuse recharge is expected to be limited to high altitude areas with coarse‐grained soils devoid of vegetation. A recharge model for this environment is developed based on a simple soil moisture budgeting technique and the calculation of actual evaporation based on empirical studies. The model is run with data for the Linzor basins, over 4000 m elevation at 22·2°S on the west slope of the Andes. It is checked against independent estimates based on the chloride mass balance (CMB) method and flood events measured downstream in the Río Salado and found to provide robust and reliable results. The results indicate that irregular and volumetrically limited amounts of diffuse recharge occur at high elevations in half of all years, with a tendency to cluster during La Niña episodes. For the Linzor Basins, mean annual recharge is found to be equivalent to 28 mm a^?1, although no recharge occurs in years with precipitation less than 120 mm, and increases proportionately with annual rainfall amounts above this limit. Copyright © 2009 John Wiley & Sons, Ltd.     

Impacts of human activities on recharge in a multilayered semiarid aquifer (Campo de Cartagena,SE Spain)

The development of intense agriculture in semiarid areas modifies intensity and spatial distribution of groundwater recharge by summing irrigation return flow to limited rainfall infiltration. Environmental tracers provide key information, but their interpretation is complicated by more complex groundwater flow patterns. In multilayered aquifers, the real origin of the groundwater samples is hard to assess because of local mixing processes occurring inside long‐screened boreholes. We use environmental tracers (¹⁴C, ¹³C, ²H, ¹⁸O, ³H) to investigate the long‐term evolution of recharge in the five‐layer Campo de Cartagena aquifer in South‐Eastern Spain, in addition to high‐resolution temperature loggings to identify the depth of origin of groundwater. Despite the complex background, this methodology allowed a reliable interpretation of the geochemistry and provided a better understanding of the groundwater flow patterns. The tritium method did not give good quantitative results because of the high variability of the recharge signal but remained an excellent indicator of recent recharge. Nonetheless, both pre‐anthropization and post‐anthropization recharge regime could be identified and quantified by radiocarbon. Before the development of agriculture, recharge varied from 17 mm.year^‐1 at the mountain ranges to 6 mm.year^‐1 in the plain, whereas the mean annual rainfall is about 300 mm. In response to the increase of agricultural activity, recharge fluxes to the plain were amplified and nowadays reach up to 210 mm.year^‐1 in irrigated areas. These values are strengthened by global water budget and local unsaturated zone studies. Copyright © 2013 John Wiley & Sons, Ltd.     

Hydro‐climatic variability and trends in Washington State for the last 50 years

Historical records of monthly streamflow and precipitation coupled with mean, minimum, and maximum air temperatures for Washington State were used to study the variation and the trend characteristics that occurred over the last 50 years (1952–2002). Results indicate that the 1967 statewide water resource assessment needs to be updated because all of the stations used in that study exhibited a decreasing trend in annual streamflow ranging from ?0·9% to ?49·3%, with an arithmetic mean of ?11·7% and a median value of ?9·8%. Furthermore, a slightly decreasing trend in annual streamflow, although not statistically significant, was detected. The decreasing streamflow magnitude was about ?1·178 mm year^?2, or 4·88 m³ s^?1 year^?1, which caused a decrease in annual streamflow in the state of about 58·9 mm, or 244 m³ s^?1. This magnitude was about 9·6% of the average annual streamflow for the entire state from 1952 to 2002. Contrastingly, the overall annual precipitation in the entire state increased 1·375 mm year^?2. Overall the annual means of daily mean, maximum, and minimum temperature increased by 0·122, 0·048, and 0·185 °C/10 years, respectively, during the study period. Thus the corresponding annual means of daily mean, maximum, and minimum temperatures increased by 0·61, 0·24, and 0·93 °C, respectively. All of these trends and magnitudes were found to vary considerably from station to station and month to month. The possible reasons resulting in these detected trends include, but are not limited to, human activities, climate variability and changes, and land use and land cover changes. Copyright © 2009 John Wiley & Sons, Ltd.     

Atmospheric chloride deposition in continental Spain

The atmospheric bulk deposition rate of chloride in continental Spain was studied to get basic information in order to help in the evaluation of diffuse recharge to aquifers through an environmental chemical balance. Both new, recent data and bibliographic data have been used. Most sampling records are less than 5 years long and often only 1 year long. This means that the calculated mean yearly bulk deposition rate of chloride is quite uncertain by 30% on average, and larger than the values derived form records up to 15 years long. A map of atmospheric bulk deposition of chloride has been drawn using ordinary kriging. The mean bulk deposition rate of chloride varies from 1 to 30 g m^?2 year^?1 in coastal areas, with strong negative landward gradients between 0·1 and 1 g m^?2 year^?1 km^?1. In the centre of the Iberian Peninsula, chloride deposition rates vary from 0·2 to 0·5 g m^?2 year^?1, with gradients around or less than 5 × 10^?3 g m^?2 year^?1 km^?1. The coefficient of variation of the mean bulk atmospheric deposition rate of chloride, for any place, ranges from 0·1 to 1. Values larger than ～0·5 are not a good indicator of natural uncertainty for this series of data that has a skewed distribution. The map of bulk deposition rate and its error is one of the terms needed for aquifer recharge estimation by means of the chloride ion balance. Copyright © 2008 John Wiley & Sons, Ltd.     

Modelling the effects of land‐use change on water and salt delivery from a catchment affected by dryland salinity in south‐east Australia

A comprehensive framework for the assessment of water and salt balance for large catchments affected by dryland salinity is applied to the Boorowa River catchment (1550 km²), located in south‐eastern Australia. The framework comprised two models, each focusing on a different aspect and operating on a different scale. A quasi‐physical semi‐distributed model CATSALT was used to estimate runoff and salt fluxes from different source areas within the catchment. The effects of land use, climate, topography, soils and geology are included. A groundwater model FLOWTUBE was used to estimate the long‐term effects of land‐use change on groundwater discharge. Unlike conventional salinity studies that focus on groundwater alone, this study makes use of a new approach to explore surface and groundwater interactions with salt stores and the stream. Land‐use change scenarios based on increased perennial pasture and tree‐cover content of the vegetation, aimed at high leakage and saline discharge areas, are investigated. Likely downstream impacts of the reduction in flow and salt export are estimated. The water balance model was able to simulate both the daily observed stream flow and salt load at the catchment outlet for high and low flow conditions satisfactorily. Mean leakage rate of about 23·2 mm year^?1 under current land use for the Boorowa catchment was estimated. The corresponding mean runoff and salt export from the catchment were 89 382 ML year^?1 and 38 938 t year^?1, respectively. Investigation of various land‐use change scenarios indicates that changing annual pastures and cropping areas to perennial pastures is not likely to result in substantial improvement of water quality in the Boorowa River. A land‐use change of about 20% tree‐cover, specifically targeting high recharge and the saline discharge areas, would be needed to decrease stream salinity by 150 µS cm^?1 from its current level. Stream salinity reductions of about 20 µS cm^?1 in the main Lachlan River downstream of the confluence of the Boorowa River is predicted. The FLOWTUBE modelling within the Boorowa River catchment indicated that discharge areas under increased recharge conditions could re‐equilibrate in around 20 years for the catchment, and around 15 years for individual hillslopes. Copyright © 2004 John Wiley & Sons, Ltd.     

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.     

Likely effects of climate change on groundwater availability in a Mediterranean region of Southeastern Spain

Groundwater resources are typically the main fresh water source in arid and semi‐arid regions. Natural recharge of aquifers is mainly based on precipitation; however, only heavy precipitation events (HPEs) are expected to produce appreciable aquifer recharge in these environments. In this work, we used daily precipitation and monthly water level time series from different locations over a Mediterranean region of Southeastern Spain to identify the critical threshold value to define HPEs that lead to appreciable aquifer recharge in this region. Wavelet and trend analyses were used to study the changes in the temporal distribution of the chosen HPEs (≥20 mm day^?1) over the observed period 1953–2012 and its projected evolution by using 18 downscaled climate projections over the projected period 2040–2099. The used precipitation time series were grouped in 10 clusters according to similarities between them assessed by using Pearson correlations. Results showed that the critical HPE threshold for the study area is 20 mm day^?1. Wavelet analysis showed that observed significant seasonal and annual peaks in global wavelet spectrum in the first sub‐period (1953–1982) are no longer significant in the second sub‐period (1983–2012) in the major part of the ten clusters. This change is because of the reduction of the mean HPEs number, which showed a negative trend over the observed period in nine clusters and was significant in five of them. However, the mean size of HPEs showed a positive trend in six clusters. A similar tendency of change is expected over the projected period. The expected reduction of the mean HPEs number is two times higher under the high climate scenario (RCP8.5) than under the moderate scenario (RCP4.5). The mean size of these events is expected to increase under the two scenarios. The groundwater availability will be affected by the reduction of HPE number which will increase the length of no aquifer recharge periods (NARP) accentuating the groundwater drought in the region. Copyright © 2016 John Wiley & Sons, Ltd.     

Hydrological processes controlling ground and surface water flow from a hypermaritime forest–peatland complex,Diana Lake Provincial Park,British Columbia,Canada

The proposed harvesting of previously undeveloped forests in north coastal British Columbia requires an understanding of hydrological responses. Hydrometric and isotopic techniques were used to examine the hydrological linkages between meteoric inputs to the surface‐groundwater system and runoff response patterns of a forest‐peatland complex. Quickflow accounted for 72–91% of peak storm discharge. The runoff ratio was lowest for open peatland areas with thick organic horizons (0·02–0·05) due to low topographic gradients and many surface depressions capable of retaining surface water. Runoff ratio increased comparatively for ephemeral surface seep flows (0·06–0·40) and was greatest in steeply sloping forest communities with more permeable soils (0·33–0·69). The dominant mechanism for runoff generation was saturated shallow subsurface flow. Groundwater fluxes from the organic horizon of seeps (1·70–1·72 m³ day^?1 m^?1) were an important component of quickflow. The homogeneous δ²H? δ¹⁸O composition of groundwater indicated attenuation of the seasonal rainfall signal by mixing during recharge. The positive correlation (r² = 0·64 and 0·38, α = 0·05) between slope index and δ¹⁸O values in groundwater suggests that the spatial pattern in the δ¹⁸O composition along the forest‐peatland complex is influenced by topography and provides evidence that topographic indices may be used to predict groundwater residence time. Copyright © 2006 John Wiley & Sons, Ltd.     

Climatic factors controlling chemical and isotopic characteristics of precipitation in Syria

The chemical and isotopic composition of monthly composite rain samples collected at 13 meteoric stations in Syria during two hydrological cycles from 1991 to 1993 have been measured. The chemical analysis of the samples revealed at a number of stations pollution due to industry and sand storms. The temporal and spatial variation of the isotopic composition has been found to be comparable with one of the neighbouring countries such as Jordan and others. The average weighted oxygen‐18 and deuterium values are − 7·5 and − 39·11‰, respectively, and the deuterium excess is 21‰ . The individual isotope values can be divided into two groups. One group is represented by winter precipitation and fits closely the Mediterranean Meteorological Water Line (MMWL). Thus, for winter precipitation, condensation of Mediterranean water vapour appears to be the dominating isotope fractionation process. The other group represents spring precipitation and is spread along an evaporation line below the MMWL, thus indicating the influence of sub‐cloud evaporation. The d‐excess has been found to be lower in the north of Syria (19·9‰ at Tartous, 18·1‰ at Jarablous) than in the south of the country (23·4‰ at Sweida, 24·1‰ at Izra) where Mediterranean air mass dominates. The d‐excess of precipitation in neighbouring countries is also close to the average value of the eastern Mediterranean basin of 22‰ , e.g. for Jordan the value is 23‰ , which signifies that Mediterranean water vapour is, for all these countries, the dominant source of precipitation. The tritium content of precipitation was found to increase with distance from the coast (5·3 TU at coastal station Tartous, 8·8 TU at continental station Palmyra). Low tritium content and high d‐excess at coastal stations clearly indicate a Mediterranean air moisture source. Copyright © 2010 John Wiley & Sons, Ltd.     

A water balance model to estimate climate change impact on groundwater recharge in Yucatan Peninsula,Mexico

ABSTRACT

The aim of this paper is to estimate the effect that climate change will have on groundwater recharge at the Yucatan Peninsula, Mexico. The groundwater recharge is calculated from a monthly water balance model considering eight methods of potential and actual evapotranspiration. Historical data from 1961–2000 and climate model outputs from five downscaled general circulation models in the near horizon (2015–2039), with representative concentration pathway (RCP) 4.5 and 8.5 are used. The results estimate a recharge of 118 ± 33 mm·year^–1 (around 10% of precipitation) in the historical period. Considering the uncertainty from GCMs under different RCP and evapotranspiration scenarios, our monthly water balance model estimates a groundwater recharge of 92 ± 40 mm·year^–1 (RCP4.5) and 94 ± 38 mm·year^–1 (RCP8.5) which represent a reduction of 23% and 20%, respectively, a result that threatens the socio-ecological balance of the region.     

A catchment framework for one‐dimensional models: introducing FLUSH and its application

A new modelling framework capable of incorporating detailed one‐dimensional models in a catchment context is presented which can be used to asses the hydrological implications (recharge, discharge, salt movement) of different land uses on different parts of the catchment. The modelling framework incorporates farming systems models and, thus, simulates crop and pasture production, whilst also accounting for lateral fluxes of water (surface and subsurface) and groundwater recharge and discharge. The framework was applied to Simmons Creek catchment, a subcatchment of the Billabong Creek in southern New South Wales, comprising gentle uplands and substantial low‐relief areas containing swamps. An integrated approach incorporating soil, hydrology, hydrogeology, and terrain analysis resulted in interpretation of landscape function and the necessary parameterization of the modelling framework. Current land use (crop rotation and pasture) and an alternative land use (10% trees on uphill units and pasture in the lower lying lands) were simulated to compare the relative contribution of parts of the catchment with total recharge. Comparison between current and alternative land use over 44 years of simulations indicated a decrease of mean annual drainage from 39 to 29 mm year^?1 and an average reduction of the groundwater level of about 0·4 m. A more substantial decrease in water‐table depth would require targeted tree planting over larger areas. This can be investigated further with the spatial framework. Copyright © 2007 John Wiley & Sons, Ltd.