Characterization of the shallow groundwater system in an alpine watershed: Handcart Gulch,Colorado, USA

Water-table elevation measurements and aquifer parameter estimates are rare in alpine settings because few wells exist in these environments. Alpine groundwater systems may be a primary source of recharge to regional groundwater flow systems. Handcart Gulch is an alpine watershed in Colorado, USA comprised of highly fractured Proterozoic metamorphic and igneous rocks with wells completed to various depths. Primary study objectives include determining hydrologic properties of shallow bedrock and surficial materials, developing a watershed water budget, and testing the consistency of measured hydrologic properties and water budget by constructing a simple model incorporating groundwater and surface water for water year 2005. Water enters the study area as precipitation and exits as discharge in the trunk stream or potential recharge for the deeper aquifer. Surficial infiltration rates ranged from 0.1–6.2×10^?5 m/s. Discharge was estimated at 1.28×10^?3 km³. Numerical modeling analysis of single-well aquifer tests predicted lower specific storage in crystalline bedrock than in ferricrete and colluvial material (6.7×10^?5–2.0×10^?3 l/m). Hydraulic conductivity in crystalline bedrock was significantly lower than in colluvial and alluvial material (4.3×10^?9–2.0×10^?4 m/s). Water budget results suggest that during normal precipitation and temperatures water is available to recharge the deeper groundwater flow system.     

Groundwater balance in the Khor Arbaat basin, Red Sea State, eastern Sudan

The Khor Arbaat basin is the main source of potable water supply for the more than 750,000 inhabitants of Port Sudan, eastern Sudan. The variation in hydraulic conductivity and storage capacity is due to the heterogeneity of the sediments, which range from clay and silt to gravely sand and boulders. The water table rises during the summer and winter rainy seasons; it reaches its lowest level in the dry season. The storage capacity of the Khor Arbaat aquifer is estimated to be 21.75?×?10⁶ m³. The annual recharge through the infiltration of flood water is about 1.93?×?10⁶ m³. The groundwater recharge, calculated as underground inflow at the ‘upper gate’, is 1.33?×?10⁵ m³/year. The total annual groundwater recharge is 2.06?×?10⁶ m³. The annual discharge through underground outflow at the ‘lower gate’ (through which groundwater flows onto the coastal plain) is 3.29?×?10⁵ m³/year. Groundwater discharge due to pumping from Khor Arbaat basin is 4.38?×?10⁶ m³/year on average. The total annual groundwater discharge is about 4.7?×?10⁶ m³. A deficit of 2.6?×?10⁶ m³/year is calculated. Although the total annual discharge is twice the estimated annual recharge, additional groundwater flow from the fractured basement probably balances the annual groundwater budget since no decline is observed in the piezometric levels.     

Spatial characteristics of groundwater temperature in the Ishikari Lowland, Hokkaido, northern Japan: analytical and numerical applications

In porous sediments of the Ishikari Lowland, there is a gradual increase in the background geothermal gradient from the Ishikari River (3–4 °C 100 m^–1) to the southwest highland area (10 °C 100 m^–1). However, the geothermal gradient at shallow depths differs in detail from the background distribution. In spite of convective heat-flow loss generally associated with groundwater flow, heat flow remains high (100 mW m^–2) in the recharge area in the southwestern part of the Ishikari basin, which is part of an active geothermal field. In the northeastern part of the lowland, heat flow locally reaches 140 mW m^–2, probably due to upward water flow from the deep geothermal field. Between the two areas the heat flow is much lower. To examine the role of hydraulic flow in the distortion of the isotherms in this area, thermal gradient vs. temperature analyses were made, and they helped to define the major components of the groundwater-flow system of the region. Two-dimensional simulation modeling aided in understanding not only the cause of horizontal heat-flow variations in this field but also the contrast between thermal properties of shallow and deep groundwater reservoirs. Electronic Publication     

Groundwater sustainability and groundwater/surface-water interaction in arid Dunhuang Basin,northwest China

The Dunhuang Basin, a typical inland basin in northwestern China, suffers a net loss of groundwater and the occasional disappearance of the Crescent Lake. Within this region, the groundwater/surface-water interactions are important for the sustainability of the groundwater resources. A three-dimensional transient groundwater flow model was established and calibrated using MODFLOW 2000, which was used to predict changes to these interactions once a water diversion project is completed. The simulated results indicate that introducing water from outside of the basin into the Shule and Danghe rivers could reverse the negative groundwater balance in the Basin. River-water/groundwater interactions control the groundwater hydrology, where river leakage to the groundwater in the Basin will increase from 3,114?×?10⁴ m³/year in 2017 to 11,875?×?10⁴ m³/year in 2021, and to 17,039?×?10⁴ m³/year in 2036. In comparison, groundwater discharge to the rivers will decrease from 3277?×?10⁴ m³/year in 2017 to 1857?×?10⁴ m³/year in 2021, and to 510?×?10⁴ m³/year by 2036; thus, the hydrology will switch from groundwater discharge to groundwater recharge after implementing the water diversion project. The simulation indicates that the increased net river infiltration due to the water diversion project will raise the water table and then effectively increasing the water level of the Crescent Lake, as the lake level is contiguous with the water table. However, the regional phreatic evaporation will be enhanced, which may intensify soil salinization in the Dunhuang Basin. These results can guide the water allocation scheme for the water diversion project to alleviate groundwater depletion and mitigate geo-environmental problem.     

Effect of thermal-hydrogeological and borehole heat exchanger properties on performance and impact of vertical closed-loop geothermal heat pump systems

Ground-source geothermal systems are drawing increasing attention and popularity due to their efficiency, sustainability and being implementable worldwide. Consequently, design software and regulatory guidelines have been developed. Interaction with the subsurface significantly affects the thermal performance, sustainability, and impacts of such systems. Reviewing the related guidelines and the design software, room for improvement is evident, especially in regards to interaction with groundwater movement. In order to accurately evaluate the thermal effect of system and hydrogeological properties on a borehole heat exchanger, a fully discretized finite-element model is used. Sensitivity of the loop outlet temperatures and heat exchange rates to hydrogeological, system and meteorological factors (i.e. groundwater flux, thermal conductivity and volumetric heat capacity of solids, porosity, thermal dispersivity, grout thermal conductivity, background and inlet temperatures) are analyzed over 6-month and 25-year operation periods. Furthermore, thermal recovery during 25  years after system decommissioning has been modeled. The thermal plume development, transport and dissipation are also assessed. This study shows the importance of subsurface thermal conductivity, groundwater flow (flux > 10^?7 m/s), and background and inlet temperature on system performance and impact. It also shows the importance of groundwater flow (flux > 10^?8 m/s) on thermal recovery of the ground over other factors.     

Assessment of surface water potential and groundwater recharge in ungauged watersheds: a case study in Tamil Nadu,India

Many of the states in India have been facing water scarcity for more than 2 decades due to increased demand, because of the increase in population and higher living standards. Consequently, many states have almost fully utilized the available surface water resources and are exploiting groundwater to augment water supplies. Investigations were carried out in the upper Thurinjalar watershed of Ponnaiyar basin in Tamil Nadu to determine the availability of surface water and to investigate the potential for enhancing groundwater recharge to support the water demand in the watershed. Increasing the water availability would also enable the community to convert the 46% of the land area in the watershed that is currently underutilised into productive uses. The surface water potential for the upper Thurinjalar watershed was assessed by applying the USDA–NRCS model with daily time steps. This modelling exercise indicated that the annual runoff from the 323 km² area of the watershed is 61 million m³. Groundwater recharge in the watershed was assessed by carrying out daily water balance method and indicated that about 43 million m³ of water from recharge is available on an annual basis or about 14% of annual rainfall. A simple regression model was developed to compute groundwater recharge from rainfall based on water balance computations and this was statistically verified. The modelling indicated that there is sufficient water available in the watershed to support current land uses and to increase the productivity of underutilised land in the area. The study also demonstrates that simple regression models can be used as an effective tool to compute groundwater recharge for ungauged basins with proper calibration.     

Study of geothermal water intrusion due to groundwater exploitation in the Puebla Valley aquifer system, Mexico

Significant intrusion of geothermal water into fresh groundwater takes place in the Puebla Valley aquifer system, Mexico. The decline in the potentiometric surface due to the overexploitation of the groundwater induces this intrusion. This hydrological system comprises three aquifers located in Plio-Quaternary volcanic sediments and Mesozoic calcareous rocks. The hydraulic balance of the aquifer shows that the annual output exceeds the natural inputs by 12 million m³. Between 1973 and 2002, a drop in the potentiometric surface, with an 80 m cone of depression, was identified in a 5-km-wide area located southwest of the city of Puebla. Chemical analyses performed on water samples since 1990 have shown an increase in total dissolved solids (TDS) of more than 500 mg/L, coinciding with the region showing a cone of depression in the potentiometric surface. A three-dimensional flow and transport model, based on the hydrogeological and geophysical studies, was computed by using the MODFLOW and MT3D software. This model reproduces the evolution of the aquifer system during the last 30 years and predicts for 2010 an additional drawdown in the potentiometric surface of 15 m, and an increase in the geothermal water intrusion.     

Hydraulic response of a tidally forced coastal aquifer, Pontal do Paraná, Brazil

The unconfined aquifer supported by the extensive Holocene coastal plain of southern Brazil is found dissected at several places by tidally forced streams. One such unit of coastal plain, falling within the watershed limits of Perequê Stream in Paraná state, was analyzed using MODFLOW to assess the surface and subsurface water components under different stress conditions. The steady-state simulation suggests that ~23,000 m³/day of water forms the total budget of the designated aquifer zone and this zone discharges ~16,000 m³/day across the interface between the freshwater aquifer zone and the sea as direct submarine discharge. The maximum evaporation loss was assumed to be 2,700 m³/day. Any change in effective aquifer storage was insignificant, as there is no large-scale exploitation of local groundwater reserves. A calibration exercise, based on field measurements lasting 1 year, supports validity of the model considered. Simulations representing two different transient conditions point to the dynamic hydraulic response exhibited by the aquifer. A schematic map, showing directions and velocity of subsurface flows, and a three-dimensional model of the groundwater reserve helped to visualize the hydrogeologic changes and to formulate management plans.     

Fractured-karst spring-flow protections: a case study in Jinan, China

Jinan Springs are an important historical heritage in China and have been well known for hundreds of years. Over-abstraction of groundwater in the Jinan area has seriously endangered the hydrological system of the springs, which have stopped flowing for significant periods of time in recent years. A three-dimensional finite-element model programmed primarily at the Hefei University of Technology has been developed to simulate groundwater level change in the large fractured-karst aquifer system in the Jinan Springs field. Various spring protection plans have been explored and their effects on the water table analyzed and compared. It was found that the present rate of groundwater withdrawal from the fractured-karst system in this area was inappropriate for spring protection. The simulated results suggest that decreasing the rate of groundwater pumping from 6.9×10⁵ to 2.7×10⁵ m³/day is needed to protect spring flows. Additional water resource requirements in that area may be met by use of surface water and recycled waste water.     

Geochemical patterns and origin of alkaline thermal waters in Central Greece (Platystomo and Smokovo areas)

This study investigates the origin and chemical composition of the thermal waters of Platystomo and Smokovo areas in Central Greece as well as any possible relationships of them to the neighboring geothermal fields located in the south-eastern part of Sperchios basin. The correlations between different dissolved salts and the temperature indicate that the chemical composition of thermal waters are controlled by, the mineral dissolution and the temperature, the reactions due to CO₂ that originates possibly by diffusion from the geothermal fields of Sperchios basin and the mixing of thermal waters with fresh groundwater from karst or shallow aquifers. Two major groups of waters are recognized on the basis of their chemistry: thermal waters of Na–HCO₃–Cl type and thermal waters mixed with fresh groundwater of Ca–Mg–Na–HCO₃ type. All thermal waters of the study area are considered as modified by water–rock interaction rainwater, heated in depth and mixed in some cases with fresh groundwater when arriving to the surface. Trace elements present low concentrations. Lithium content suggests discrimination between the above two groups of waters. Boron geochemistry confirms all the above remarks. Boron concentration ranges from 60 μg L^?1 to 10 mg L^?1, while all samples’ constant isotopic composition (δ¹¹B ≈ 10 ‰) indicates leaching from rocks. The positive correlation between the chemical elements and the temperature clearly indicates that much of the dissolved salts are derived from water–rock interactions. The application of geothermometers suggests that the reservoir temperature is around 100–110 °C. Chalcedony temperatures are similar to the emergent temperatures and this is typical of convective waters in fault systems in normal thermal gradient areas.     

Groundwater budget for the upper and middle parts of the River Gash Basin, eastern Sudan

The River Gash Basin is filled by the Quaternary alluvial deposits, unconformably overlying the basement rocks. The alluvial deposits are composed mainly of unconsolidated layers of gravel, sand, silt, and clays. The aquifer is unconfined and is laterally bounded by the impermeable Neogene clays. The methods used in this study include the carry out of pumping tests and the analysis of well inventory data in addition to the river discharge rates and other meteorological data. The average annual discharge of the River Gash is estimated to be 1,056?×?10⁶ m³ at El Gera gage station (upstream) and 587?×?10⁶ m³ at Salam-Alikum gage station (downstream). The annual loss mounts up to 40% of the total discharge. The water loss is attributed to infiltration and evapotranspiration. The present study proofs that the hydraulic conductivity ranges from 36 to 105 m/day, whereas the transmissivity ranges from 328 to 1,677 m²/day. The monitoring of groundwater level measurements indicates that the water table rises during the rainy season by 9 m in the upstream and 6 m in the midstream areas. The storage capacity of the upper and middle parts of the River Gash Basin is calculated as 502?×?10⁶ m³. The groundwater input reach 386.11?×?10⁶ m³/year, while the groundwater output is calculated as 365.98?×?10⁶ m³/year. The estimated difference between the input and output water quantities in the upper and middle parts of the River Gash Basin demonstrates a positive groundwater budget by about 20?×?10⁶ m³/year     

Temperature-driven groundwater convection in cold climates

The aim was to study density-driven groundwater flow and analyse groundwater mixing because of seasonal changes in groundwater temperature. Here, density-driven convection in groundwater was studied by numerical simulations in a subarctic climate, i.e. where the water temperature was <4 °C. The effects of soil permeability and groundwater temperature (i.e. viscosity and density) were determined. The influence of impermeable obstacles in otherwise homogeneous ground was also studied. An initial disturbance in the form of a horizontal groundwater flow was necessary to start the convection. Transient solutions describe the development of convective cells in the groundwater and it took 22 days before fully developed convection patterns were formed. The thermal convection reached a maximum depth of 1.0 m in soil of low permeability (2.71 · 10^?9 m²). At groundwater temperature close to its density maximum (4 °C), the physical size (in m) of the convection cells was reduced. Small stones or frost lenses in the ground slightly affect the convective flow, while larger obstacles change the size and shape of the convection cells. Performed simulations show that “seasonal groundwater turnover” occurs. This knowledge may be useful in the prevention of nutrient leakage to underlying groundwater from soils, especially in agricultural areas where no natural vertical groundwater flow is evident. An application in northern Sweden is discussed.     

Quantification of the water balance and hydrogeological processes of groundwater–lake interactions in the Pampa Plain, Argentina

This paper gives an account of the assessment and quantification of the water balance and the hydrogeological processes related to lake–groundwater interaction in the Pampa Plain by using hydrogeochemical, isotopic and flow numerical modeling techniques. La Salada is a permanent shallow lake, with an area of 5.8 km², located on the SE of Buenos Aires Province. A total of 29 lake water samples and 15 groundwater samples were collected for both hydrochemical analysis and environmental stable isotope determination. Water table depths were measured in wells closed to the lake. Groundwater samples appear grouped on the Local Meteoric Water Line, suggesting a well-mixed system and that rainfall is the main recharge source to the aquifer. Water evaporation process within La Salada is also corroborated by its isotopic composition. The model that best adjusts to La Salada Lake hydrochemical processes includes evaporation from groundwater, calcite precipitation with CO₂ release and cationic exchange. The annual water balance terms for the lake basin indicates for each hydrological component the following values: 1.16 E⁰⁸ m³ rainfall, 8.15 E⁰⁷ m³ evapotranspiration, 1.90 E⁰⁶ m³ runoff, 1.55 E⁰⁷ m³ groundwater recharge, 6.01 E⁰⁶ m³ groundwater discharge to the lake, 9.54 E⁰⁶ m³ groundwater discharge to the river, 5.00 E⁰⁵ m³ urban extraction and 4.90 E⁰⁶ m³ lake evaporation. Integrated analysis of hydrochemical and isotopic information helped to calibrate the groundwater flow model, to validate the conceptual model and to quantitatively assess the basin water balance.     

Estimates of heat flow and heat production and a thermal model of the São Francisco craton

An updated analysis of geothermal data from the highland area of eastern Brazil has been carried out and the characteristics of regional variations in geothermal gradients and heat flow examined. The database employed includes results of geothermal measurements at 45 localities. The results indicate that the Salvador craton and the adjacent metamorphic fold belts northeastern parts of the study area are characterized by geothermal gradients in the range of 6–17°C/km. The estimated heat flow values fall in the range of 28–53 mW/m², with low values in the cratonic area relative to the fold belts. On the other hand, the São Francisco craton and the intracratonic São Francisco sedimentary basin in the southwestern parts are characterized by relatively higher gradient values, in the range of 14–42°C/km, with the corresponding heat flow values falling in the range of 36–89 mW/m². Maps of regional variations indicate that high heat flow anomaly in the São Francisco craton is limited to areas of sedimentary cover, to the west of the Espinhaço mountain belt. Crustal thermal models have been developed to examine the implications of the observed intracratonic variations in heat flow. The thermal models take into consideration variation of thermal conductivity with temperature as well as change of radiogenic heat generation with depth. Vertical distributions of seismic velocities were used in obtaining estimates of radiogenic heat production in crustal layers. Crustal temperatures are calculated based on a procedure that makes simultaneous use of the Kirchoff and Generalized Integral Transforms, providing thereby analytical solutions in 2D and 3D geometry. The results point to temperature variations of up to 300°C at the Moho depth, between the northern Salvador and southern São Francisco cratons. There are indications that differences in rheological properties, related to thermal field, are responsible for the contrasting styles of deformation patterns in the adjacent metamorphic fold belts.     

Groundwater recharge mechanism in an integrated tableland of the Loess Plateau,northern China: insights from environmental tracers

Assessing groundwater recharge characteristics (recharge rate, history, mechanisms (piston and preferential flow)) and groundwater age in arid and semi-arid environments remains a difficult but important research frontier. Such assessments are particularly important when the unsaturated zone (UZ) is thick and the recharge rate is limited. This study combined evaluations of the thick UZ with those of the saturated zone and used multiple tracers, such as Cl, NO₃, Br, ²H, ¹⁸O, ¹³C, ³H and ¹⁴C, to study groundwater recharge characteristics in an integrated loess tableland in the Loess Plateau, China, where precipitation infiltration is the only recharge source for shallow groundwater. The results indicate that diffuse recharge beneath crops, as the main land use of the study area, is 55–71 mm yr^?1 based on the chloride mass balance of soil profiles. The length of time required for annual precipitation to reach the water table is 160–400 yrs. The groundwater is all pre-modern water and paleowater, with corrected ¹⁴C age ranging from 136 to 23,412 yrs. Most of the water that eventually becomes recharge originally infiltrated in July–September. The Cl and NO₃ contents in the upper UZ are considerably higher than those in the deep UZ and shallow groundwater because of recent human activities. The shallow groundwater has not been in hydraulic equilibrium with present near-surface boundary conditions. The homogeneous material of the UZ and relatively old groundwater age imply that piston flow is the dominant recharge mechanism for the shallow groundwater in the tableland.     

Modeling of groundwater flow for Mujib aquifer, Jordan

Jordan is an arid country with very limited water resources. Groundwater is the main source for its water supply. Mujib aquifer is located in the central part of Jordan and is a major source of drinking water for Amman, Madaba and Karak cities. High abstraction rates from Mujib aquifer during the previous years lead to a major decline in water levels and deterioration in groundwater quality. Therefore, proper groundwater management of Mujib aquifer is necessary; and groundwater flow modeling is essential for proper management. For this purpose, Modflow was used to build a groundwater flow model to simulate the behavior of the flow system under different stresses. The model was calibrated for steady state condition by matching observed and simulated initial head counter lines. Drawdown data for the period 1985–1995 were used to calibrate the transient model by matching simulated drawdown with the observed one. Then, the transient model was validated by using drawdown data for the period 1996–2002. The results of the calibrated model showed that the horizontal hydraulic conductivity of the B2/A7 aquifer ranges between 0.001 and 40m/d. Calibrated specific yield ranges from 0.0001 to 0.15. The water balance for the steady state condition of Mujib aquifer indicated that the total annual direct recharge is 20.4 × 10⁶m³, the total annual inflow is 13.0 × 10⁶ m³, springs discharge is 15.3 × 10⁶ m³, and total annual outflow is 18.7 × 10⁶ m³. Different scenarios were considered to predict aquifer system response under different conditions. The results of the sensitivity analysis show that the model is highly sensitive to horizontal hydraulic conductivity and anisotropy and with lower level to the recharge rates. Also the model is sensitive to specific yield     

Geochemical,isotopic and hydrological mass balance approaches to constrain the lake water–groundwater interaction in Dal Lake,Kashmir Valley

It is important to have qualitative as well as quantitative understanding of the hydraulic exchange between lake and groundwater for effective water resource management. Dal, a famous urban fresh water lake, plays a fundamental role in social, cultural and economic dynamics of the Kashmir Valley. In this paper geochemical, isotopic and hydrological mass balance approaches are used to constrain the lake water–groundwater interaction of Dal Lake and to identify the sources of lake water. Water samples of precipitation (n = 27), lake water (n = 18) and groundwater (n = 32) were collected across the lake and its catchment for the analysis of δ¹⁸O and δ²H. A total of 444 lake water samples and 440 groundwater samples (springs, tube wells and dug wells) were collected for the analysis of Ca²⁺, Mg²⁺, HCO₃ ^?, SO₄ ^2?, Cl^?, NO₃ ^?, Na⁺ and K⁺. Water table and lake water level were monitored at 40 observation locations in the catchment. Water table map including pH and EC values corroborate and verify the gaining nature of the Dal Lake. Stable isotopes of lake water in Boddal and Gagribal basins showed more deviation from the global meteoric water line than Hazratbal and Nigeen basins, indicating the evaporation of lake water. The isotopic and geochemical mass balance suggested that groundwater contributes a significant proportion (23–40%) to Dal Lake. The estimated average groundwater contribution to Dal Lake ranged from 31.2 × 10³ to 674 × 10³ m³ day^?1 with an average of 276 × 10³ m³ day^?1. The study will be useful to delineate the possible sources of nutrients and pollutants entering the lake and for the management of lake water resources for sustainable development.     

An assessment of the potential and impacts of winter water banking in the Sokh aquifer, Central Asia

The dynamics of artificial recharge of winter surface flows coupled with increased summer groundwater use for irrigation in the Sokh aquifer (Central Asia) have been investigated. Water release patterns from the giant Toktogul reservoir have changed, as priority is now given to hydropower generation in winter in Kyrgyzstan. Winter flows have increased and summer releases have declined, but the Syr Darya River cannot pass these larger winter flows and the excess is diverted to a natural depression, creating a 40?×?10⁹m³ lake. A water balance study of all 18 aquifers feeding the Fergana Valley indicated the feasibility of winter groundwater recharge in storage created by summer abstraction. This modeling study examines the dynamics of the process in one aquifer over a 5-year period, with four scenarios: the current situation; increased groundwater abstraction of around 625 million (M) m³/year; groundwater abstraction with an artificial recharge of 144 Mm³/year, equivalent to the volume available in low flow years in the Sokh River; and with a larger artificial recharge of 268 Mm³/year, corresponding to high flow availability. Summer surface irrigation diversions can be reduced by up to 350 Mm³ and water table levels can be lowered.     

Groundwater-driven temperature changes at thermal springs in response to recent glaciation: Bormio hydrothermal system,Central Italian Alps

Thermal springs are widespread in the European Alps, with hundreds of geothermal sites known and exploited. The thermal circulation and fluid outflows were examined in the area around Bormio (Central Italian Alps), where ten geothermal springs discharge from dolomite bodies located close to the regional Zebrù thrust. Water is heated in deep circulation systems and upwells vigorously at a temperature of about 40 °C. Heat and fluid transport is explored by steady and transient three-dimensional finite-element simulations taking into account the effect of the last glaciation, which in the study area was recognized to end around 11,000–12,000 years ago. The full regional model (ca. 700 km²) is discretized with a highly refined triangular finite-element planar grid. Numerical simulations suggest a reactivation of the system following the end of the Last Glacial Maximum. Results correctly simulate the observed discharge rate of ca. 2,400 L/min and the spring temperatures after ca. 13,000 years from deglaciation, and show a complete cooling of the aquifer within a period of approximately 50,000 years. Groundwater flow and temperature patterns suggest that thermal water flows through a deep system crossing both sedimentary and metamorphic lithotypes along a fracture network associated with the thrust system. This example gives insights into the influences of deep alpine structures and glaciations on groundwater circulation that control the development of many hydrothermal systems not necessarily associated with convective heat flow.     

Stochastic groundwater modeling of a sedimentary aquifer: evaluation of the impacts of abstraction scenarios under conditions of reduced recharge

A transient finite difference groundwater flow model has been calibrated for the Nasia sub-catchment of the White Volta Basin. This model has been validated through a stochastic parameter randomization process and used to evaluate the impacts of groundwater abstraction scenarios on resource sustainability in the basin. A total of 1500 equally likely model realizations of the same terrain based on 1500 equally likely combinations of the data of the key aquifer input parameters were calibrated and used for the scenario analysis. This was done to evaluate model non-uniqueness arising from uncertainties in the key aquifer parameters especially hydraulic conductivity and recharge by comparing the realizations and statistically determining the degree to which they differ from each other. Parameter standard deviations, computed from the calibrated data of the key parameters of hydraulic conductivity and recharge, were used as a yardstick for evaluating model non-uniqueness. All model realizations suggest horizontal hydraulic conductivity estimates in the range of 0.03–78.4 m/day, although over 70 % of the area has values in the range of 0.03–14 m/day. Low standard deviations of the horizontal hydraulic conductivity estimates from the 1500 solutions suggest that this range adequately reflects the properties of the material in the terrain. Lateral groundwater inflows and outflows appear to constitute significant components of the groundwater budgets in the terrain, although estimated direct vertical recharge from precipitation amounts to about 7 % of annual precipitation. High potential for groundwater development has been suggested in the simulations, corroborating earlier estimates of groundwater recharge. Simulation of groundwater abstraction scenarios suggests that the domain can sustain abstraction rates of up to 200 % of the current estimated abstraction rates of 12,960 m³/day under the current recharge rates. Decreasing groundwater recharge by 10 % over a 20-year period will not significantly alter the results of this abstraction scenario. However, increasing abstraction rates by 300 % over the period with decreasing recharge by 10 % will lead to drastic drawdowns in the hydraulic head over the entire terrain by up to 6 m and could cause reversals of flow in most parts of the terrain.