Effects of horizontal strain in estimating specific storage and compaction in confined and leaky aquifer systems

 Mathematical and numerical distinctions are made between developments of governing equations involving groundwater flow and granular displacement where (1) only vertical strain and no change in total normal load are assumed, and (2) such limitations are eliminated. The former is referred to here as the Jacob-Terzaghi method; the latter is referred to as the Biot method. Numerical comparisons of drawdown, compaction, and strain are made for hypothetical confined and leaky aquifer systems for one- and three-dimensional strain problems. Simulation results indicate that incorporation of horizontal strain results in reduced drawdowns. Vertical compaction for a specified stress distribution is greatly reduced when horizontal strain is invoked. The distribution of strain components indicates that approximately two-thirds of the total volume strain originates from horizontal compaction in pumped aquifers. Incorporation of horizontal strain also affects the distribution of flux through the overlying confining layer as a function of distance from the pumping well and may ultimately affect the interpretation of aquifer-test data. Received, January 1999 / Revised, May 1999, August 1999 / Accepted, August 1999     

Scale aspects of groundwater flow and transport systems

Flow-system analysis is based on the concept of hierarchical groundwater flow systems. The topography of the water table, which is strongly related to the topography of the land surface, is a major factor in the hierarchical nesting of gravity-driven groundwater flow, resulting in flow systems of different orders of magnitude in lateral extent and depth of penetration. The concept of flow systems is extremely useful in the analysis of spatial and temporal scales and their mutual relationships. Basic equations on the laboratory scale are extended to larger, regional scales. Making use of Fourier analysis further develops Tóth's original idea of topography-driven flow systems. In this way, the different spatial scales of the water table are separated in a natural way, leading to a simple expression for the penetration depth of a flow system. This decomposition leads also to the relationship between spatial and temporal scales. Analogous to flow systems, water bodies with different water quality may be called 'transport systems.' Field studies, numerical micro-scale modeling over macro-scale domains, and stochastic dispersion theory indicate that between systems with steady transport, the interfaces are relatively thin. The interfaces are much thinner than the relatively large mixing zones predicted by the conventional engineering approach to macrodispersion, in which relatively large, time-independent macrodispersion lengths are applied. A relatively simple alternative engineering approach is presented. For macrodispersion of propagating solute plumes, the alternative dispersion term gives the same results as the conventional engineering approach and gives correct results for steady-state transport.     

Validity of a sharp-interface model in a confined coastal aquifer

 The problem of seawater intrusion is considered for the case of a confined coastal aquifer in which there is steady seaward flow of fresh water. Using the GWCH2O model, the problem is solved first for the case of no dispersion where a distinct interface exists separating the fresh water from the salt water. The problem is solved next by taking into account dispersion and diffusion of the salt-water component, along with the density effect. In this respect, a two-dimensional finite-element model, 2D-VDTRAN, is developed to simulate density-dependent solute transport. To investigate the limitation of the sharp-interface approach in coastal aquifers for conditions of both steady state and unsteady state, the problem is solved twice using the two models with different parameter values. These parameters are combined in dimensionless form, resulting in four named parameters: seepage factor (A); dispersion-to-advection ratio (B); geometry ratio (C); and time-scale factor (T). Using the density-dependent model, the dimensionless width of the transition zone (W/L) is determined for different values of A, B, C, and T. Steady-state simulations show that the sharp-interface approach is valid only when the system is dominated by advection, i.e., when 0<B≤5% for all values of A and C, or when A≥65% for all values of B and C. However, the unsteady-state analysis shows that the applicability of the sharp-interface approach is sufficiently accurate at early times. Received, October 1997 Revised, June 1998, October 1998 Accepted, November 1998     

Reconstruction of natural groundwater flow paths in the multiple aquifer system of the northern Negev (Israel), based on lithological and structural evidence

 The Judea Group, a limestone and dolomite karstic aquifer of late Albian–Turonian age, is one of the most important sources of water in Israel. In the western part of the country, the Judea Group aquifer is also known as the Yarkon–Taninim basin. In the northern Negev, the Judea Group is a recipient for fresh water flowing southward from the Hebron Mountains and of brackish paleowater flowing northward from Sinai. Very little is known of the hydraulic properties of this aquifer. In order to outline assumed natural flow paths that existed in this karstic environment prior to groundwater exploitation, use was made of lithological, structural, and paleomorphological features. A detailed hydrogeological conceptual model of the Judea Group aquifer in northern Negev was established by the geological interpretation of high-resolution seismic reflection and by analysis of lithological evidence from boreholes. Isopach, isolith-contour, and isolith-ratio maps were compiled for the main lithological components. Increase in transmissivity values is inversely proportional with the cumulative thickness of argillaceous components. The lithological and hydraulic evidence provides the basis for subdividing the subsurface into distinctive permeability zones for the upper and lower sections of the aquifer; for outlining possible preferential groundwater flow paths for both subaquifers; and for improving understanding of groundwater-salinty variations that result from lithological variability, direction of groundwater flow paths, groundwater flow rates, and the duration of rock/water interactions. In an earlier conceptual model of the basin, the Judea Group aquifer was regarded as a continuous and undisturbed entity. The present study reveals an intricate groundwater flow pattern that is controlled by lithological and structural factors that create zones of preferential flow. This interpretation bears on the overall evaluation of groundwater resources and their management and exploitation. Received, December 1996 · Revised, October 1997, June 1998 · Accepted, July 1998     

Variations in hydraulic conductivity with scale of measurement during aquifer tests in heterogeneous, porous carbonate rocks

 Previous studies have shown that hydraulic conductivity of an aquifer seems to increase as the portion of the aquifer tested increases. To date, such studies have all relied on different methods to determine hydraulic conductivity at each scale of interest, which raises the possibility that the observed increase in hydraulic conductivity is due to the measurement method, not to the scale. This study analyzes hydraulic conductivity with respect to scale during individual aquifer tests in porous, heterogeneous carbonate rocks in southeastern Wisconsin, USA. Results from this study indicate that hydraulic conductivity generally increases during an individual test as the volume of aquifer impacted increases, and the rate of this increase is the same as the rate of increase determined by using different measurement methods. Thus, scale dependence of hydraulic conductivity during single tests does not depend on the method of measurement. This conclusion is supported by 22 of 26 aquifer tests conducted in porous-flow-dominated carbonate units within the aquifer. Instead, scale dependency is probably caused by heterogeneities within the aquifer, a conclusion supported by digital simulation. All of the observed types of hydraulic-conductivity variations with scale during individual aquifer tests can be explained by a conceptual model of a simple heterogeneous aquifer composed of high-conductivity zones within a low-conductivity matrix. Received, January 1997 Revised, August 1997, November 1997 Accepted, November 1997     

On the solutions of modeling a slug test performed in a two-zone confined aquifer

The slug test is a popular field method for estimating aquifer parameters. As a finite-thickness skin is present at a wellbore, the result of test-data analysis may represent average hydraulic properties of the skin and formation zones. Moench and Hsieh (Int. Assoc. of Hydrogeology, The IAH 17th International Congress on Hydrology of Rocks of Low Permeability, Int. Assoc. of Hydrogeology. Tucson, Arizona, 7–12 January) provided the Laplace-domain solution for a confined aquifer with a finite-thickness skin that satisfies the governing equations, boundary conditions, and continuity requirements at the interface of the skin and formation zones. Two new formulae presented here relate the modified zero-order Bessel functions of the first and second kinds to their first and second derivatives. Their solution is evaluated using the modified Crump approach for the Laplace inversion. The numerical values obtained by the modified Crump approach are compared with others based on finite-element modeling. The results indicate that the values based on finite-element modeling exhibit poor accuracy, attributed to errors in the mass balance equation and interface condition in their mathematical model.

---

Resumen Los ensayos slug o tipo pulso son métodos populares de campo para estimar parámetros de acuíferos. Debido a que en los sondeos existe una membrana de espesor finito, el resultado del análisis de los datos del ensayo pueden representar las propiedades hidráulicas promedio de las zonas de la formación y dicha membrana. Moench AF y Hsieh PA (Asociación Internacional de Hidrogeólogos, 7–12 Enero 1985) Análisis de datos de prueba slug en un pozo con membrana de espesor finito. 17 Congreso Internacional de la Asociación Internacional de Hidrogeólogos sobre Hidrología de Rocas de Baja Permeabilidad. Asociación Internacional de Hidrogeólogos. Tucson, Arizona. 7–12 Enero. han aportado la solución en el ámbito de Laplace para un acuífero confinado con una membrana de espesor finito que satisface los requerimientos de las ecuaciones dominantes, las condiciones restrictivas, y de continuidad en la interfase de las zonas de formación y membrana. Las dos nuevas fórmulas que se presentan aquí relacionan las funciones modificadas Bessel de orden cero de la primera y segunda clase con su primera y segunda derivadas. Se evalúa su solución usando el enfoque modificado de Crump para la inversión Laplace. Los valores numéricos obtenidos mediante el enfoque modificado de Crump se comparan con otros valores numéricos basados en modelos de elementos finitos. Los resultados indican que los valores basados en el modelo de elementos finitos presentan poca precisión, que se atribuye a errores en la ecuación de balance de masa y a las condiciones de interfase en su modelo matemático.

---

Résumé Le slug-test représente une méthode populaire de terrain pour estimer les paramètres d’un aquifère. Du fait de la présence d’une couche à l’épaisseur limitée sur un site de forages, le résultat de l’analyse des données du test pourrait représenté le terrain mais aussi cette couche. Moench AF et Hsieh PA (Association Internationale d’Hydrogéologie, 7–12 Jan 1985) Analyses de données de slug test dans un puits présentant une couche d’épaisseur finie. 17ème congrès international de l’AIH sur l’Hydrologie des roches de faible perméabilité, Association Internationale d’Hydrogéologie. Tucson, Arizona, 7–12 Jan. apportent une solution du domaine de Laplace pour un aquifère confiné avec une épaisseur limitée qui satisfait l’équation en question, les conditions aux limites, et les besoins de continuité à l’interface de la couche et de l’encaissant. Deux nouvelles formules présentées dans cet article met en relation les fonctions de Bessel modifiées d’ordre zéro du premier et second type, avec leur dérivées première et deuxième. Leurs solutions sont évaluées en utilisant une approche de Crump modifiée pour l’inversion de Laplace. Les valeurs numériques obtenues avec l’approche de Crump sont comparées avec d’autres approches utilisant les modèles aux éléments finis. Les résultats indiquent que les valeurs basées sur les modèles aux éléments finis montrent une faible précision, du fait d’erreurs dans le bilan hydrologique et les conditions aux limites dans le modèle.

     

Application of a solute transport model under variable velocity conditions in a conduit flow aquifer: Olalde karst system, Basque Country, Spain

A model based on numerical solutions, which allows for solving the dispersion equation under variable recharge and velocity conditions, is developed to simulate solute transport in conduit flow aquifers during flow recession periods. As an example, the evolution of a tracer in the little known karst conduit that links the sinking stream of Oma valley to the Olalde spring is investigated in the karstic region of Santa Eufemia-Ereñozar (Basque Country, Spain). The model, with different hypothetical structures, allows for obtaining series of tracer breakthrough curves, which are fitted to experimental data using an optimization algorithm. These results, although they can be used to simulate the tracer evolution between the two points considered, do not allow for determining the internal structure and spatial disposition of contributions in the aquifer.     

Geomorphic aspects of groundwater flow

 The many roles that groundwater plays in landscape evolution are becoming more widely appreciated. In this overview, three major categories of groundwater processes and resulting landforms are considered: (1) Dissolution creates various karst geometries, mainly in carbonate rocks, in response to conditions of recharge, geologic setting, lithology, and groundwater circulation. Denudation and cave formation rates can be estimated from kinetic and hydraulic parameters. (2) Groundwater weathering generates regoliths of residual alteration products at weathering fronts, and subsequent exhumation exposes corestones, flared slopes, balanced rocks, domed inselbergs, and etchplains of regional importance. Groundwater relocation of dissolved salts creates duricrusts of various compositions, which become landforms. (3) Soil and rock erosion by groundwater processes include piping, seepage erosion, and sapping, important agents in slope retreat and headward gully migration. Thresholds and limits are important in many chemical and mechanical groundwater actions. A quantitative, morphometric approach to groundwater landforms and processes is exemplified by selected studies in carbonate and clastic terrains of ancient and recent origins. Received, May 1998 · Revised, September 1998 · Accepted, October 1998     

Relation between "terra rossa" from the Apulia aquifer of Italy and the radon content of groundwater: Experimental results and their applicability to radon occurrence in the aquifer

 The radon-222 (²²²Rn) activity in groundwater of the Apulian karstic aquifer in southern Italy is as great as 500 Becquerel per liter (Bq/L) locally. Normal radium-226 (²²⁶Ra) activity in the limestone and calcareous dolomites of the aquifer is not enough to explain such a high level. Laboratory investigations identified high ²²⁶Ra activity in the "terra rossa," the residuum occupying fissures and cavities in the bedrock, and also the relation between (1) ²²⁶Ra-bearing bedrock and "terra rossa" and (2) ²²²Rn in water. The "terra rossa" is the primary source of the radon in the groundwater. The experimental results show the need to characterize the "terra rossa" of Apulia on the basis of ²²⁶Ra activity and also to study the distribution and variations in ²²²Rn activity over time in the aquifer. Received, January 1996 · Revised, March 1997, September 1997, June 1998 · Accepted, July 1998     

Numerical analysis of the hydrogeologic controls in a layered coastal aquifer system, Oahu, Hawaii, USA

 The coastal aquifer system of southern Oahu, Hawaii, USA, consists of highly permeable volcanic aquifers overlain by weathered volcanic rocks and interbedded marine and terrestrial sediments of both high and low permeability. The weathered volcanic rocks and sediments are collectively known as caprock, because they impede the free discharge of groundwater from the underlying volcanic aquifers. A cross-sectional groundwater flow and transport model was used to evaluate the hydrogeologic controls on the regional flow system in southwestern Oahu. Controls considered were: (a) overall caprock hydraulic conductivity; and (b) stratigraphic variations of hydraulic conductivity in the caprock. Within the caprock, variations in hydraulic conductivity, caused by stratigraphy or discontinuities of the stratigraphic units, are a major control on the direction of groundwater flow and the distribution of water levels and salinity. Results of cross-sectional modeling confirm the general groundwater flow pattern that would be expected in a layered coastal system. Groundwater flow is: (a) predominantly upward in the low-permeability sedimentary units; and (b) predominantly horizontal in the high-permeability sedimentary units. Received, October 1996 Revised, August 1997 Accepted, September 1997     

承压含水层非完整反滤回灌井的稳定流计算

反滤回灌井是地下水库中一种有效的回灌设施,它由回灌井和井口反滤层组成。非完整反滤回灌井井流运动由砂反滤层的竖向流和非完整回灌井井流组合而成。文中提出了承压含水层非完整反滤回灌井稳定流的计算方法,分析了反滤层、井损和回灌堵塞等因素对回灌量的影响,并结合实例证明了反滤回灌井回灌量大大低于普通回灌井的回灌量。     

Influence of groundwater on the degradation of irrigated soils in a semi-arid region, the inner delta of the Niger River, Mali

The problem of soil degradation through alkalinization/salinization in an irrigated area with a semi-arid climate was examined in the inner delta of the Niger River, Mali, by the study of groundwater hydraulics and hydrochemistry in an area recharged by irrigation water. On the basis of data analysis on various scales, it is concluded that the current extent of the surface saline soils is due to a combination of three factors: (1) the existence of ancient saline soils (solonchaks) resulting from the creation of a broad sabkha west of the former course of the Niger River, now called the Fala of Molodo. These saline crusts were gradually deposited during the eastward tilting of the tectonic block that supports the Niger River; (2) the irrigation processes during the recent reflooding of the Fala of Molodo (river diversion in 1950). These used very poorly mineralized surface water but reintroduced into the alluvial groundwater system – generally of a low permeability (K=10–6?m?s–1) – salts derived from the ancient solonchaks; and (3) the redeposition of the dissolved salts on the surface due to the intense evapotranspiration linked to the present Sahelian climate. In this context, only efficient artificial draining of subsurface alluvial groundwater can eliminate most of the highly mineralized flow and thus reduce the current saline deposits.     

Groundwater flow rate and contaminant migration in fissure-karstic aquifer of Opole Triassic system southwest Poland

 This paper presents hydrogeological problems occurring during municipal water exploitation and mine dewatering. These activities result in groundwater quantity and quality changes in the fissure-karstic aquifer. Increase of nitrate concentration up to 12 mg NNO₃/l due to intensive fertilizer use, and high tritium concentration, show water system impact up to 100 m depth. Intensive water exploitation produces large cones of depression with over 40-m water-level depletion in the Opole region. Flow rates of major components and isotopes have been verified by chemical migration history. Some aspects of the protection policy of this type of aquifer are also discussed. Received: 7 March 1997 · Accepted: 17 November 1998     

侧向有限越流承压含水层中非完整井非稳定流模型及解析解

工程建设中当距离抽水井r=rb处水位基本没有变化或不受抽水影响时,或当此处存在止水帷幕时,含水层系统视为侧向有限延伸,rb为有限半径。为此,构建更加符合工程实际的侧向有限延伸的典型弱透水层-承压水层系统中非完整井非稳定流计算模型,同时考虑井径和井储效应的影响,应用Laplace变换和分离变量法得到了水位降深在拉氏空间下的解析解,并应用拉氏数值逆变换Stehfest法得到真实空间下的水位降深。新建立的解析解可以进一步退化为诸多已有解,并进一步将其与已知解和有限元数值解进行对比,验证了所得解的正确性和可靠性。基于新建解重点分析了侧向边界和井的完整性对承压水层水位降深的影响。结果表明:含水层系统的侧向有限边界仅对抽水后期的水位降深影响明显,含水层系统侧向无限延伸情况下的水位降深要大于情形1(在r=rb处为定水头边界)且明显小于情形2(在r=rb处为不透水边界)下的水位降深,rb越小,两者之间的误差越大;抽水井的完整性对整个抽水期间不同情形下的水位降深均有明显的影响,承压含水层顶板处的水位降深随着抽水井滤管的长度和埋深的增加而减小。     

Plan for a groundwater monitoring network in Taiwan

 In Taiwan, rapid economic growth, rising standards of living, and an altered societal structure have in recent years put severe demands on water supplies. Because of its stable quantity and quality, groundwater has long been a reliable source of water for domestic, agricultural, and industrial users, but the establishment of a management program that integrates groundwater and surface-water use has been hampered by the lack of groundwater data. In 1992, the Department of Water Resources (DWR) initiated a program entitled "Groundwater Monitoring Network Plan in Taiwan." Under this program, basic groundwater data, including water-level and water-quality data, are being collected, and a reliable database is being established for the purpose of managing total water resources. This paper introduces the goals, implementation stages, and scope of that plan. The plan calls for constructing 517 hydrogeologic survey stations and 990 groundwater monitoring wells within 17 years. Under this program, water-level fluctuations are continuously monitored, whereas water-quality samples are taken for analysis only at the initial drilling stage and, subsequently, at the time when a monitoring well is being serviced. In 1996, the DWR and the Water Resources Planning Commission were merged to form today's Water Resources Bureau. Received, July 1997 · Revised, February 1998 · Accepted, January 1998     

Assessing aquifer susceptibility to and severity of atrazine contamination at a field site in south-central Wisconsin, USA

 A field study from October 1989 through July 1992, conducted on a 4.1-km² area in south-central Wisconsin, USA, examined the distributions of atrazine and its chlorinated metabolites in groundwater and related those distributions to the groundwater flow system. MODFLOW and PATH3D were used to assess bedrock-aquifer susceptibility to contamination. Estimated travel time from water table to bedrock surface ranges from <0.25 to >512 yr. Spatial distribution of the estimates demonstrates that increased travel time to bedrock can result from the presence of shallow surface-water bodies, greater depths to bedrock, and smaller hydraulic conductivities. Estimated travel times to local domestic wells are inversely related to atrazine and desethylated atrazine concentrations observed in water from those wells. The potential impact of long-term atrazine use on aquifer water quality was investigated using MT3D in two best-case scenarios. Uncertainties associated with predicted atrazine concentrations at various depths and times were estimated. For shallow groundwater, widespread violations of Wisconsin's current preventive action limit were predicted, but with large uncertainty stemming from uncertain estimates of input parameter values. The simulations indicate, however, that moderate inputs at the water table are very unlikely to produce violations of Wisconsin's standards deeper in the aquifer. Received, October 1997 Revised, July 1998 Accepted, July 1998     

Analysis of radial movement of a confined aquifer due to well pumping and injection

Three aspects of radial motion of an aquifer solid matrix, driven by hydraulic forces due to well recharge and discharge, have been investigated mathematically. Firstly, the analytic solutions in the velocity and displacement fields were found using the definition of bulk mass flow, conservation of water and solid mass, the Darcy-Gersevanov law, and the drawdown of hydraulic head for unsteady flow. Secondly, the radial velocity and displacement of aquifer solid matrix were investigated with respect to time and location in response to three pumping-injection schemes (i.e., linear, nonlinear, and sinusoidal well flow rates). Aquifer movement was found to be more responsive to linear and nonlinear pumping rates than to constant pumpage. Aquifer movement in response to periodic well discharge and recharge was governed by an integral sinusoidal function. Finally, the potential risk of earth fissuring at a given location caused by tensile stress, in turn induced by well discharge-recharge activity, has been discussed using the analytic solution. The deformation compression and extension zones derived from analysis of the results indicated the possible occurrence of tensile stress that may cause earth fissures. Sample computations, using assumed parameters, provided a first-hand evaluation of the location of potential fissuring in relation to the well.     

Groundwater as a geologic agent: An overview of the causes, processes, and manifestations

 The objective of the present paper is to show that groundwater is a general geologic agent. This perception could not, and did not, evolve until the system nature of basinal groundwater flow and its properties, geometries, and controlling factors became recognized and understood through the 1960s and 1970s. The two fundamental causes for groundwater's active role in nature are its ability to interact with the ambient environment and the systematized spatial distribution of its flow. Interaction and flow occur simultaneously at all scales of space and time, although at correspondingly varying rates and intensities. Thus, effects of groundwater flow are created from the land surface to the greatest depths of the porous parts of the Earth's crust, and from a day's length through geologic times. Three main types of interaction between groundwater and environment are identified in this paper, with several special processes for each one, namely: (1) Chemical interaction, with processes of dissolution, hydration, hydrolysis, oxidation-reduction, attack by acids, chemical precipitation, base exchange, sulfate reduction, concentration, and ultrafiltration or osmosis; (2) Physical interaction, with processes of lubrication and pore-pressure modification; and (3) Kinetic interaction, with the transport processes of water, aqueous and nonaqueous matter, and heat. Owing to the transporting ability and spatial patterns of basinal flow, the effects of interaction are cumulative and distributed according to the geometries of the flow systems. The number and diversity of natural phenomena that are generated by groundwater flow are almost unlimited, due to the fact that the relatively few basic types are modified by some or all of the three components of the hydrogeologic environment: topography, geology, and climate. The six basic groups into which manifestations of groundwater flow have been divided are: (1) Hydrology and hydraulics; (2) Chemistry and mineralogy; (3) Vegetation; (4) Soil and rock mechanics; (5) Geomorphology; and (6) Transport and accumulation. Based on such a diversity of effects and manifestations, it is concluded that groundwater is a general geologic agent. Received, December 1998 · Revised, January 1999 · Accepted, January 1999     

Effects of groundwater flow on mineral diagenesis, with emphasis on carbonate aquifers

 This article provides a critical synopsis of the effects of groundwater flow on mineral diagenesis. Emphasis is placed on those aspects and processes that change porosity and permeability in carbonate aquifers, because they are of particular importance to human societies as sources of supplies of water for human consumption (drinking, irrigation) and of crude oil and natural gas. Diagenetic settings in carbonates as well as clastics are generally ill defined. This paper proposes a new comprehensive classification of diagenetic settings into near-surface, shallow-, intermediate-, and deep-burial diagenetic settings; hydrocarbon-contaminated plumes; and fractures. These settings are defined on the basis of mineralogy, petroleum, hydrogeochemistry, and hydrogeology. This classification is applicable to all sedimentary basins. Diagenesis is governed by various intrinsic and extrinsic factors that include thermodynamic and kinetic constraints, as well as microstructural factors that may override the others. These factors govern diagenetic processes, such as dissolution, compaction, recrystallization, replacement, and sulfate–hydrocarbon redox-reactions. Processes such as cementation, dissolution, and dolomitization require significant flow of groundwater driven by an externally imposed hydraulic gradient. Other processes, such as stylolitization and thermochemical sulfate reduction, commonly take place without significant groundwater flow in hydrologically nearly or completely stagnant systems that are geochemically "closed." Two major effects of groundwater flow on mineral diagenesis are enhancement and reduction of porosity and permeability, although groundwater flow can also leave these rock properties essentially unchanged. In extreme cases, an aquifer or hydrocarbon reservoir rock can have highly enhanced porosity and permeability due to extensive mineral dissolution, or it can be plugged up due to extensive mineral precipitation. Received, April 1998 · Revised, July 1998 · Accepted, September 1998