Hydraulic conductivity fields: Gaussian or not?

Hydraulic conductivity (K) fields are used to parameterize groundwater flow and transport models. Numerical simulations require a detailed representation of the K field, synthesized to interpolate between available data. Several recent studies introduced high-resolution K data (HRK) at the Macro Dispersion Experiment (MADE) site, and used ground-penetrating radar (GPR) to delineate the main structural features of the aquifer. This paper describes a statistical analysis of these data, and the implications for K field modeling in alluvial aquifers. Two striking observations have emerged from this analysis. The first is that a simple fractional difference filter can have a profound effect on data histograms, organizing non-Gaussian ln K data into a coherent distribution. The second is that using GPR facies allows us to reproduce the significantly non-Gaussian shape seen in real HRK data profiles, using a simulated Gaussian ln K field in each facies. This illuminates a current controversy in the literature, between those who favor Gaussian ln K models, and those who observe non-Gaussian ln K fields. Both camps are correct, but at different scales.

The origin of a large (>3km) maar volcano by coalescence of multiple shallow craters: Lake Purrumbete maar, southeastern Australia

Lake Purrumbete maar is located in the intraplate, monogenetic Newer Volcanics Province in southeastern Australia. The extremely large crater of 3000. m in diameter formed on an intersection of two fault lines and comprises at least three coalesced vents. The evolution of these vents is controlled by the interaction of the tectonic setting and the properties of both hard and soft rock aquifers. Lithics in the maar deposits originate from country rock formations less than 300. m deep, indicating that the large size of the crater cannot only be the result of the downwards migration of the explosion foci in a single vent. Vertical crater walls and primary inward dipping beds evidence that the original size of the crater has been largely preserved. Detailed mapping of the facies distributions, the direction of transport of base surges and pyroclastic flows, and the distribution of ballistic block fields, form the basis for the reconstruction of the complex eruption history,which is characterised by alternations of the eruption style between relatively dry and wet phreatomagmatic conditions, and migration of the vent location along tectonic structures. Three temporally separated eruption phases are recognised, each starting at the same crater located directly at the intersection of two local fault lines. Activity then moved quickly to different locations. A significant volcanic hiatus between two of the three phases shows that the magmatic system was reactivated. The enlargement of especially the main crater by both lateral and vertical growth led to the interception of the individual craters and the formation of the large circular crater. Lake Purrumbete maar is an excellent example of how complicated the evolution of large, seemingly simple, circular maar volcanoes can be, and raises the question if these systems are actually monogenetic.

A Web-based system enabling the integration, analysis, and 3D sub-surface visualization of groundwater monitoring data and geological models

The 3D Water Chemistry Atlas is an intuitive, open source, Web-based system that enables the three-dimensional (3D) sub-surface visualization of ground water monitoring data, overlaid on the local geological model (formation and aquifer strata). This paper firstly describes the results of evaluating existing virtual globe technologies, which led to the decision to use the Cesium open source WebGL Virtual Globe and Map Engine as the underlying platform. Next it describes the backend database and search, filtering, browse and analysis tools that were developed to enable users to interactively explore the groundwater monitoring data and interpret it spatially and temporally relative to the local geological formations and aquifers via the Cesium interface. The result is an integrated 3D visualization system that enables environmental managers and regulators to assess groundwater conditions, identify inconsistencies in the data, manage impacts and risks and make more informed decisions about coal seam gas extraction, waste water extraction, and water reuse.

Relationships between major ions in coal seam gas groundwaters: Examples from the Surat and Clarence-Moreton basins

Using a combination of multivariate statistical techniques and the graphical assessment of major ion ratios, the influences on hydrochemical variability of coal seam gas (or coal bed methane) groundwaters from several sites in the Surat and Clarence-Moreton basins in Queensland, Australia, were investigated. Several characteristic relationships between major ions were observed: 1) strong positive linear correlation between the Na/Cl and alkalinity/Cl ratios; 2) an exponentially decaying trend between the Na/Cl and Na/alkalinity ratios; 3) inverse linear relationships between increasing chloride concentrations and decreasing pH for high salinity groundwaters, and; 4) high residual alkalinity for lower salinity waters, and an inverse relationship between decreasing residual alkalinity and increasing chloride concentrations for more saline waters. The interpretation of the hydrochemical data provides invaluable insights into the hydrochemical evolution of coal seam gas (CSG) groundwaters that considers both the source of major ions in coals and the influence of microbial activity. Elevated chloride and sodium concentrations in more saline groundwaters appear to be influenced by organic-bound chlorine held in the coal matrix; a sodium and chloride ion source that has largely been neglected in previous CSG groundwater studies. However, contrastingly high concentrations of bicarbonate in low salinity waters could not be explained, and are possibly associated with a number of different factors such as coal degradation, methanogenic processes, the evolution of high-bicarbonate NaHCO3 water types earlier on in the evolutionary pathway, and variability in gas reservoir characteristics. Using recently published data for CSG groundwaters in different basins, the characteristic major ion relationships identified for new data presented in this study were also observed in other CSG groundwaters from Australia, as well as for those in the Illinois Basin in the USA. This observation suggests that where coal maceral content and the dominant methanogenic pathway are similar, and where organic-bound chlorine is relatively abundant, distinct hydrochemical responses may be observed. Comparisons with published data of other NaHCO3 water types in non-CSG environments suggest that these characteristic major ion relationships described here can: i) serve as an indicator of potential CSG groundwaters in certain coal-bearing aquifers that contain methane; and ii) help in the development of strategic sampling programmes for CSG exploration and to monitor potential impacts of CSG activities on groundwater resources.

Characterization of alluvial formation by stochastic modelling of paleo-fluvial processes: The concept and method

Modelling fluvial processes is an effective way to reproduce basin evolution and to recreate riverbed morphology. However, due to the complexity of alluvial environments, deterministic modelling of fluvial processes is often impossible. To address the related uncertainties, we derive a stochastic fluvial process model on the basis of the convective Exner equation that uses the statistics (mean and variance) of river velocity as input parameters. These statistics allow for quantifying the uncertainty in riverbed topography, river discharge and position of the river channel. In order to couple the velocity statistics and the fluvial process model, the perturbation method is employed with a non-stationary spectral approach to develop the Exner equation as two separate equations: the first one is the mean equation, which yields the mean sediment thickness, and the second one is the perturbation equation, which yields the variance of sediment thickness. The resulting solutions offer an effective tool to characterize alluvial aquifers resulting from fluvial processes, which allows incorporating the stochasticity of the paleoflow velocity.

Groundwater balance research to improve the accuracy of drawdown predictions in the Surat Basin, Queensland

Groundwater modelling studies rely on an accurate determination of inputs and outputs that make up the water balance. Often there is large uncertainty associated with estimates of recharge and unmetered groundwater use. This can translate to equivalent uncertainty in the forecasting of sustainable yields, impacts of extraction, and susceptibility of groundwater dependent ecosystems. In the case of Coal Seam Gas, it is important to characterise the temporal and special distribution of depressurisation in the reservoir and how this may or may not extend to the adjacent aquifers. A regional groundwater flow model has been developed by the Queensland Government to predict drawdown impacts due to Coal Seam Gas activities in the Surat basin. This groundwater model is undergoing continued refinement and there is currently scope to address some of the key areas of uncertainty including better quantification of groundwater recharge and unmetered groundwater extractions. Research is currently underway to improve the accuracy of estimates of both of these components of the groundwater balance in order to reduce uncertainty in predicted groundwater drawdowns due to CSG activities.

Memory of past random wave conditions in submarine groundwater discharge

Submarine groundwater discharge (SGD) is an integral part of the hydrological cycle and represents an important aspect of land-ocean interactions. We used a numerical model to simulate flow and salt transport in a nearshore groundwater aquifer under varying wave conditions based on yearlong random wave data sets, including storm surge events. The results showed significant flow asymmetry with rapid response of influxes and retarded response of effluxes across the seabed to the irregular wave conditions. While a storm surge immediately intensified seawater influx to the aquifer, the subsequent return of intruded seawater to the sea, as part of an increased SGD, was gradual. Using functional data analysis, we revealed and quantified retarded, cumulative effects of past wave conditions on SGD including the fresh groundwater and recirculating seawater discharge components. The retardation was characterized well by a gamma distribution function regardless of wave conditions. The relationships between discharge rates and wave parameters were quantifiable by a regression model in a functional form independent of the actual irregular wave conditions. This statistical model provides a useful method for analyzing and predicting SGD from nearshore unconfined aquifers affected by random waves

Isotopic fingerprints in surficial waters : Stable isotope methods applied in hydrogeological studies

The driving force behind this study has been the need to develop and apply methods for investigating the hydrogeochemical processes of significance to water management and artificial groundwater recharge. Isotope partitioning of elements in the course of physicochemical processes produces isotopic variations to their natural reservoirs. Tracer property of the stable isotope abundances of oxygen, hydrogen and carbon has been applied to investigate hydrogeological processes in Finland. The work described here has initiated the use of stable isotope methods to achieve a better understanding of these processes in the shallow glacigenic formations of Finland. In addition, the regional precipitation and groundwater records will supplement the data of global precipitation, but as importantly, provide primary background data for hydrological studies. The isotopic composition of oxygen and hydrogen in Finnish groundwaters and atmospheric precipitation was determined in water samples collected during 1995 2005. Prior to this study, no detailed records existed on the spatial or annual variability of the isotopic composition of precipitation or groundwaters in Finland. Groundwaters and precipitation in Finland display a distinct spatial distribution of the isotopic ratios of oxygen and hydrogen. The depletion of the heavier isotopes as a function of increasing latitude is closely related to the local mean surface temperature. No significant differences were observed between the mean annual isotope ratios of oxygen and hydrogen in precipitation and those in local groundwaters. These results suggest that the link between the spatial variability in the isotopic composition of precipitation and local temperature is preserved in groundwaters. Artificial groundwater recharge to glaciogenic sedimentary formations offers many possibilities to apply the isotopic ratios of oxygen, hydrogen and carbon as natural isotopic tracers. In this study the systematics of dissolved carbon have been investigated in two geochemically different glacigenic groundwater formations: a typical esker aquifer at Tuusula, in southern Finland and a carbonate-bearing aquifer with a complex internal structure at Virttaankangas, in southwest Finland. Reducing the concentration of dissolved organic carbon (DOC) in water is a primary challenge in the process of artificial groundwater recharge. The carbon isotope method was used to as a tool to trace the role of redox processes in the decomposition of DOC. At the Tuusula site, artificial recharge leads to a significant decrease in the organic matter content of the infiltrated water. In total, 81% of the initial DOC present in the infiltrated water was removed in three successive stages of subsurface processes. Three distinct processes in the reduction of the DOC content were traced: The decomposition of dissolved organic carbon in the first stage of subsurface flow appeared to be the most significant part in DOC removal, whereas further decrease in DOC has been attributed to adsorption and finally to dilution with local groundwater. Here, isotope methods were used for the first time to quantify the processes of DOC removal in an artificial groundwater recharge. Groundwaters in the Virttaankangas aquifer are characterized by high pH values exceeding 9, which are exceptional for shallow aquifers on glaciated crystalline bedrock. The Virttaankangas sediments were discovered to contain trace amounts of fine grained, dispersed calcite, which has a high tendency to increase the pH of local groundwaters. Understanding the origin of the unusual geochemistry of the Virttaankangas groundwaters is an important issue for constraining the operation of the future artificial groundwater plant. The isotope ratios of oxygen and carbon in sedimentary carbonate minerals have been successfully applied to constrain the origin of the dispersed calcite in the Virttaankangas sediments. The isotopic and chemical characteristics of the groundwater in the distinct units of aquifer were observed to vary depending on the aquifer mineralogy, groundwater residence time and the openness of the system to soil CO2. The high pH values of > 9 have been related to dissolution of calcite into groundwater under closed or nearly closed system conditions relative to soil CO2, at a low partial pressure of CO2.

Estimating specific yield and transmissivity with magnetic resonance sounding in an unconfined sandstone aquifer (Niger)

The unconfined aquifer of the Continental Terminal in Niger was investigated by magnetic resonance sounding (MRS) and by 14 pumping tests in order to improve calibration of MRS outputs at field scale. The reliability of the standard relationship used for estimating aquifer transmissivity by MRS was checked; it was found that the parametric factor can be estimated with an uncertainty a parts per thousand currency sign150% by a single point of calibration. The MRS water content (theta (MRS)) was shown to be positively correlated with the specific yield (Sy), and theta (MRS) always displayed higher values than Sy. A conceptual model was subsequently developed, based on estimated changes of the total porosity, Sy, and the specific retention Sr as a function of the median grain size. The resulting relationship between theta (MRS) and Sy showed a reasonably good fit with the experimental dataset, considering the inherent heterogeneity of the aquifer matrix (residual error is similar to 60%). Interpreted in terms of aquifer parameters, MRS data suggest a log-normal distribution of the permeability and a one-sided Gaussian distribution of Sy. These results demonstrate the efficiency of the MRS method for fast and low-cost prospection of hydraulic parameters for large unconfined aquifers.

Understanding Connectivity between Groundwater Chemistry Data and Geological Stratigraphy via 3D Sub-surface Visualization and Analysis

This paper describes the 3D Water Chemistry Atlas - an open source, Web-based system that enables the three-dimensional (3D) sub-surface visualization of ground water monitoring data, overlaid on the local geological model. Following a review of existing technologies, the system adopts Cesium (an open source Web-based 3D mapping and visualization interface) together with a PostGreSQL/PostGIS database, for the technical architecture. In addition a range of the search, filtering, browse and analysis tools were developed that enable users to interactively explore the groundwater monitoring data and interpret it spatially and temporally relative to the local geological formations and aquifers via the Cesium interface. The result is an integrated 3D visualization system that enables environmental managers and regulators to assess groundwater conditions, identify inconsistencies in the data, manage impacts and risks and make more informed decisions about activities such as coal seam gas extraction, waste water extraction and re-use.

Groundwater Management in India:Physical, Institutional and Policy Alternatives

Groundwater constitutes a vital natural resource for sustaining India’s agricultural economy and meeting the country’s social, ecological and environmental goals. It is a unique resource, widely available, providing security against droughts and yet it is closely linked to surface-water resources and the hydrological cycle. Its availability depends on geo-hydrological conditions and characteristics of aquifers, from deep to alluvium, sediment crystalline rocks to basalt formations; and agro-climate from humid to subhumid and semi-arid to arid. Its reliable supply, uniform quality and temperature, relative turbidity, pollution-safe, minimal evaporation losses, and low cost of development are attributes making groundwater more attractive compared to other resources. It plays a key role in the provision of safe drinking water to rural populations. For example, already almost 80% of domestic water use in rural areas in India is groundwater-supplied, and much of it is being supplied to farms, villages and small towns. Inadequate control of the use of groundwater, indiscriminate application of agrochemicals and unrestrained pollution of the rural environment by other human activities make groundwater usage unsustainable, necessitating proper management in the face of the twin demand for water of good quality for domestic supply and adequate supply for irrigation, ensuring equity, efficiency and sustainability of the resource. Groundwater irrigation has overtaken surface irrigation in the early 1980s, supported by well energization. It is estimated that there are about 24 million energised wells and tube wells now and it is driven by demand rather than availability, evident through the greater occurrence of wells in districts with high population densities. Apart from aquifer characteristics, land fragmentation and landholding size are the factors that decide the density of wells. The ‘rise and fall’ of local economies dependent on groundwater can be summarized as: the green revolution of 1980s, groundwaterbased agrarian boom, early symptoms of groundwater overdraft, and decline of the groundwater socio-ecology. The social characteristics and policy interventions typical of each stage provide a fascinating insight into the human-resource dynamics. This book is a compilation of nine research papers discussing various aspects of groundwater management. It attempts to integrate knowledge about the physical system, the socio-economic system, the institutional set-up and the policy environment to come out with a more realistic analysis of the situation with regard to the nature, characteristics and intensity of resource use, the size of the economy the use generates, and the negative socioeconomic consequences. Complex variables addressed in this regard focusing on northern Gujarat are the stock of groundwater available in the region, its hydrodynamics, its net outflows against inflows, the economics of its intensive use (particularly irrigation in semi-arid and arid regions), its criticality in the regional hydroecological regime, ethical aspects and social aspects of its use. The first chapter by Dinesh Kumar and Singh, dwells on complex groundwater socio-ecology of India, while emphasizing the need for policy measures to address indiscriminate over-exploitation of dwindling resources. The chapter also explores the nature of groundwater economy and the role of electricity prices on it. The next chapter on groundwater issue in north Gujarat provides a description of groundwater resource characteristics followed by a detailed analysis of the groundwater depletion and quality deterioration problems in the region and their undesirable consequences on the economy, ecosystem health and the society. Considering water-buyers and wellowning farmers individually, a methodology for economic valuation of groundwater in regions where its primary usage is in agriculture, and as assessment of the groundwater economy based on case studies from north Gujarat is presented in the fourth chapter. The next chapter focuses on the extent of dependency of milk production on groundwater, which includes the water embedded in green and dry fodder and animal feed. The study made a realistic estimate of irrigation water productivity in terms of the physics and economics of milk production. The sixth chapter analyses the extent of reduction in water usage, increase in yield and overall increase in physical productivity of alfalfa with the use of the drip irrigation system. The chapter also provides a detailed synthesis of the costs and benefits associated with the use of drip irrigation systems. A linear programmingbased optimization model with the objective to minimize groundwater use taking into account the interaction between two distinct components – farming and dairying under the constraints of food security and income stability for different scenarios, including shift in cropping pattern, introduction of water-efficient crops, water- saving technologies in addition to the ‘business as usual’ scenario is presented in the seventh chapter. The results show that sustaining dairy production in the region with reduced groundwater draft requires crop shifts and adoption of water-saving technologies. The eighth chapter provides evidences to prove that the presence of adequate economic incentive would encourage farmers to adopt water-saving irrigation devices, based on the findings of market research with reference to the level of awareness among farmers of technologies and the factors that decide the adoption of water-saving technologies. However, now the marginal cost of using electricity for agricultural pumping is almost zero. The economic incentives are strong and visible only when the farmers are either water-buyers or have to manage irrigation with limited water from tube-well partnerships. The ninth chapter explores the socio-economic viability of increasing the power tariff and inducing groundwater rationing as a tool for managing energy and groundwater demand, considering the current estimate of the country’s annual economic loss of Rs 320 billion towards electricity subsidy in the farm sector. The tenth chapter suggests private tradable property rights and development of water markets as the institutional tool for achieving equity, efficiency and sustainability of groundwater use. It identifies the externalities for local groundwater management and emphasizes the need for managing groundwater by local user groups, supported by a thorough analysis of groundwater socio-ecology in India. An institutional framework for managing the resource based on participatory approach that is capable of internalizing the externalities, comprising implementation of institutional and technical alternatives for resource management is also presented. Major findings of the analyses and key arguments in each chapter are summarized in the concluding chapter. Case studies of the social and economic benefits of groundwater use, where that use could be described as unsustainable, are interesting. The benefits of groundwater use are outlined and described with examples of social and economic impacts of groundwater and the negative aspects of groundwater development with the compilation of environmental problems based on up-to-date research results. This publication with a well-edited compilation of case studies is informative and constitutes a useful publication for students and professionals.

Tracing the source of bottled water using stable isotope techniques

The Indian Summer Monsoon (ISM) precipitation recharges ground water aquifers in a large portion of the Indian subcontinent. Monsoonal precipitation over the Indian region brings moisture from the Arabian Sea and the Bay of Bengal (BoB). A large difference in the salinity of these two reservoirs, owing to the large amount of freshwater discharge from the continental rivers in the case of the BoB and dominating evaporation processes over the Arabian Sea region, allows us to distinguish the isotopic signatures in water originating in these two water bodies. Most bottled water manufacturers exploit the natural resources of groundwater, replenished by the monsoonal precipitation, for bottling purposes. The work presented here relates the isotopic ratios of bottled water to latitude, moisture source and seasonality in precipitation isotope ratios. We investigated the impact of the above factors on the isotopic composition of bottled water. The result shows a strong relationship between isotope ratios in precipitation (obtained from the GNIP data base)/bottled water with latitude. The approach can be used to predict the latitude at which the bottled water was manufactured. The paper provides two alternative approaches to address the site prediction. The limitations of this approach in identifying source locations and the uncertainty in latitude estimations are discussed. Furthermore, the method provided here can also be used as an important forensic tool for exploring the source location of bottled water from other regions. Copyright (C) 2011 John Wiley & Sons, Ltd.

Causality analysis of groundwater dynamics based on a Vector Autoregressive model in the semi-arid basin of Gundal (South India)

Causal relationships existing between observed levels of groundwater in a semi-arid sub-basin of the Kabini River basin (Karnataka state, India) are investigated in this study. A Vector Auto Regressive model is used for this purpose. Its structure is built on an upstream/downstream interaction network based on observed hydro-physical properties. Exogenous climatic forcing is used as an input based on cumulated rainfall departure. Optimal models are obtained thanks to a trial approach and are used as a proxy of the dynamics to derive causal networks. It appears to be an interesting tool for analysing the causal relationships existing inside the basin. The causal network reveals 3 main regions: the Northeastern part of the Gundal basin is closely coupled to the outlet dynamics. The Northwestern part is mainly controlled by the climatic forcing and only marginally linked to the outlet dynamic. Finally, the upper part of the basin plays as a forcing rather than a coupling with the lower part of the basin allowing for a separate analysis of this local behaviour. The analysis also reveals differential time scales at work inside the basin when comparing upstream oriented with downstream oriented causalities. In the upper part of the basin, time delays are close to 2 months in the upward direction and lower than 1 month in the downward direction. These time scales are likely to be good indicators of the hydraulic response time of the basin which is a parameter usually difficult to estimate practically. This suggests that, at the sub-basin scale, intra-annual time scales would be more relevant scales for analysing or modelling tropical basin dynamics in hard rock (granitic and gneissic) aquifers ubiquitous in south India. (c) 2012 Elsevier B.V. All rights reserved.

Constraining complex aquifer geometry with geophysics (2-D ERT and MRS measurements) for stochastic modelling of groundwater flow

Stochastic modelling is a useful way of simulating complex hard-rock aquifers as hydrological properties (permeability, porosity etc.) can be described using random variables with known statistics. However, very few studies have assessed the influence of topological uncertainty (i.e. the variability of thickness of conductive zones in the aquifer), probably because it is not easy to retrieve accurate statistics of the aquifer geometry, especially in hard rock context. In this paper, we assessed the potential of using geophysical surveys to describe the geometry of a hard rock-aquifer in a stochastic modelling framework. The study site was a small experimental watershed in South India, where the aquifer consisted of a clayey to loamy-sandy zone (regolith) underlain by a conductive fissured rock layer (protolith) and the unweathered gneiss (bedrock) at the bottom. The spatial variability of the thickness of the regolith and fissured layers was estimated by electrical resistivity tomography (ERT) profiles, which were performed along a few cross sections in the watershed. For stochastic analysis using Monte Carlo simulation, the generated random layer thickness was made conditional to the available data from the geophysics. In order to simulate steady state flow in the irregular domain with variable geometry, we used an isoparametric finite element method to discretize the flow equation over an unstructured grid with irregular hexahedral elements. The results indicated that the spatial variability of the layer thickness had a significant effect on reducing the simulated effective steady seepage flux and that using the conditional simulations reduced the uncertainty of the simulated seepage flux. As a conclusion, combining information on the aquifer geometry obtained from geophysical surveys with stochastic modelling is a promising methodology to improve the simulation of groundwater flow in complex hard-rock aquifers. (C) 2013 Elsevier B.V. All rights reserved.

What mainly controls recession flows in river basins?

The ubiquity of the power law relationship between dQ/dt and Q for recession periods (-dQ/dt kQ(alpha); Q being discharge at the basin outlet at time t) clearly hints at the existence of a dominant recession flow process that is common to all real basins. It is commonly assumed that a basin, during recession events, functions as a single phreatic aquifer resting on a impermeable horizontal bed or the Dupuit-Boussinesq (DB) aquifer, and with time different aquifer geometric conditions arise that give different values of alpha and k. The recently proposed alternative model, geomorphological recession flow model, however, suggests that recession flows are controlled primarily by the dynamics of the active drainage network (ADN). In this study we use data for several basins and compare the above two contrasting recession flow models in order to understand which of the above two factors dominates during recession periods in steep basins. Particularly, we do the comparison by selecting three key recession flow properties: (1) power law exponent alpha, (2) dynamic dQ/dt-Q relationship (characterized by k) and (3) recession timescale (time period for which a recession event lasts). Our observations suggest that neither drainage from phreatic aquifers nor evapotranspiration significantly controls recession flows. Results show that the value of a and recession timescale are not modeled well by DB aquifer model. However, the above mentioned three recession curve properties can be captured satisfactorily by considering the dynamics of the ADN as described by geomorphological recession flow model, possibly indicating that the ADN represents not just phreatic aquifers but the organization of various sub-surface storage systems within the basin. (C) 2014 Elsevier Ltd. All rights reserved.