Hydrology [H]

H23C  MS:Exh Hall B   Tuesday
Advances in Hydrogeomechanics II Posters
Presiding: T A Ghezzehei, Lawrence Berkeley National Laboratory; M Berli, Desert Research Institute

H23C-1509 

Using Earthquake-triggered Groundwater Fluctuations to Estimate the Anomalous Stress- source Area

Wang, S (n4891110@ccmail.ncku.edu.tw), Department of Resources Engineering, NCKU, No.1, University Road, Tainan, 701, Taiwan * Hsu, K (kchsu@mail.ncku.edu.tw), Department of Resources Engineering, NCKU, No.1, University Road, Tainan, 701, Taiwan Wang, C (clwang@mail.ncku.edu.tw), Department of Resources Engineering, NCKU, No.1, University Road, Tainan, 701, Taiwan Lai, W (laiwenji@hotmail.com), Department of Resources Engineering, NCKU, No.1, University Road, Tainan, 701, Taiwan Lai, W (laiwenji@hotmail.com), Disaster Prevention Research Center, NCKU, 3F, No. 500, Sec.3, An-Ming Road, Tainan, 709, Taiwan Chiu, F (n4894122@mail.ncku.edu.tw), Department of Resources Engineering, NCKU, No.1, University Road, Tainan, 701, Taiwan

Earthquakes are usually triggered in the stress-concentrated area. To explain the stress-pore pressure relationship, poroelastic theory is one of the commonly used models. In this model, stress and pore pressure are coupled. We applied an analytical solution of the poroelastic model with point force source to estimate the stress-concentrated area based on earthquake-triggered groundwater fluctuations. A classification system of groundwater level anomalies is developed and step groundwater fluctuations are used for data analysis. The collected anomalous data shows that the detectable distance of groundwater level anomaly increases as the earthquake magnitude increases. The maximum detectable distance is about 250 km and the minimum detectable earthquake magnitude is about 5.0 based on the groundwater fluctuations. Semi-analytical and empirical models are also developed to construct the relations among earthquake magnitude, epicentral distance, and groundwater level anomaly. Using the principle of superposition, the stress-concentrated areas can be estimated from two earthquake events in Taiwan. The results show that earthquake-triggered groundwater fluctuations can be a potential tool to estimate the stress-concentrated area and may be used to indicate the possible area of an epicenter.

H23C-1510 

Detection and Localization of Hydromechanical Disturbances in a Sandbox using the Self-potential Method

* Crespy, A (crespy@cerege.fr), CNRS-CEREGE, Universite Aix Marseille III, Europole de l'Arbois BP 80, Aix en Provence, 13545, France Revil, A (arevil@mines.edu), Colorado School of Mines, Green Center, dept of Geophysics 1500 Illinois street, Golden, CO 80401, United States Linde, N (linde@aug.ig.erdw.ethz.ch), Swiss Federal Institute of Technology, Institute of Geophysics HPP O 2 Schafmattstr. 30, Zurich, 8093, Switzerland Byrdina, S (lana@ipgp.jussieu.fr), IPGP, 4 Place Jussieu, Paris, 75252, France Jardani, A (Abderrahim.Jardani@univ-rouen.fr), Département de Géologie et Laboratoire M2C (Rouen), Universite de Rouen Bât.IRESE A Place Emile Blondel, Mont-Saint-Aignan, 76821, France Boleve, A (boleve@cerege.fr), SOBESOL, Savoie Technolac BP 230, Le Bourget du Lac, 73370, France Henry, P (henry@cdf.u-3mrs.fr), Collège de France-Chaire geodynamique, Europôle de l'Arbois Bat Le Trocadéro - Aile Sud BP 80, Aix en Provence, 13545, France

The self-potential response during hydromechanical disturbance of water-infiltrated porous material was investigated in the laboratory. Foursandbox experiments were performed to understand the electrokinetic response associated with the injection of a pulse of water and the abstraction of a small volume of pore water. The resulting self-potential signals are measured using 32 Ag/AgCl very sensitive and non-polarising medical electrodes. The injected/abstracted volumes of water are responsible for hydro-mechanical disturbances. In turn, these disturbances generate dipolar electrical anomalies of electrokinetic nature with a distinct electrical signature and an amplitude of few V. The source function is a product of a dipolar Green's function and a source intensity function that depends solely on the product between the streaming potential coupling coefficient of the sand and the pore fluid overpressure with respect to the hydrostatic pressure. Numerical modeling is performed to solve the coupled hydro-mechanical problem and to determine the distribution of the resulting streaming potential in the course of the experiments. We use 2D and 3D algorithms based on the intercorrelation method and wavelet analysis of potential fields to show that the source corresponds to a vertical dipole. These methods are also used to localize the position of the source of the hydromechanical disturbance from the self-potential collected at the top surface of the tank. Applications concern the detection and localization of fracturing in active volcanoes.

H23C-1511 

The Effects of Gravity-Induced Stresses on Hydrogeologic Properties Near Ridges: Examples From two Rift Basins

Savage, W Z (savage@usgs.gov), U. S. Geological Survey, Denver Federal Center, Denver, CO 80225, United States * Morin, R H (rhmorin@usgs.gov), U. S. Geological Survey, Denver Federal Center, Denver, CO 80225, United States Matter, J M (jmatter@ldeo.columbia.edu), Lamont-Doherty Earth Observatory of Columbia University, 61 Route 9W, Palisades, NY 10964, United States Rivard, C (crivard@nrcan.gc.ca), Geological Survey of Canada, 490, rue de la Couronne, Quebec, QC G1K 9A9, Canada Goldberg, D S (goldberg@ldeo.columbia.edu), Lamont-Doherty Earth Observatory of Columbia University, 61 Route 9W, Palisades, NY 10964, United States

Twenty Mesozoic rift basins have been identified along the northeast coast of North America. At two of these, the Minas Basin in Nova Scotia and the Newark Basin in New York, steep topographic ridges formed initially by extensional tectonic processes have been evaluated for groundwater potential by obtaining geophysical logs in nearby water wells. The two locations have common lithologic and topographic characteristics, whereby a basalt cap overlying sedimentary rocks is exposed along a ridge flank. Despite these similarities, the basalts in Nova Scotia serve as a local water supply whereas the wells in New York are non-producing. Analyses of acoustic televiewer logs and measurements from flowmeter/pumping tests locate and distinguish the permeable fractures from the general fracture population. In addition, a complementary finite-element stress model is used to simulate ridge geometry from topographic maps and incorporates rock elastic properties as determined from processing full-waveform sonic logs. Results derived from field data and modeling indicate that gravity-induced stresses near ridges result from a lack of buttressing along the free faces and are amplified by spreading effects due to a contrast in Poisson's ratio between the basalts and the sedimentary rocks. These conditions act to enhance the permeability of fractures striking subparallel to the ridge axis in Nova Scotia. Conversely, the predominant fracture population in New York strikes perpendicular to the ridge axis and to the maximum horizontal principal stress. Consequently, there are few fractures preferentially aligned with the stress field and no hydraulically conductive features are recognized in these rocks. The hydrogeologic properties of the basalts near the ridge flanks at these two sites are markedly different because of the orientation of pre-existing fractures and their response to imposed stresses as derived from topography and lithostratigraphy.

H23C-1512 

Evaluation of the Effects of Cementation on Specific Storage of Granular Porous Media Using Discrete Element Models

* Plourde, K E (kplourde@geo.umass.edu), University of Massachusetts - Amherst, Geoscience Department, Morrill Science Center, 611 North Pleasant Street, Amherst, MA 01003, United States Boutt, D F (dboutt@geo.umass.edu), University of Massachusetts - Amherst, Geoscience Department, Morrill Science Center, 611 North Pleasant Street, Amherst, MA 01003, United States Goodwin, L B (laurel@geology.wisc.edu), University of Wisconsin - Madison, Department of Geology and Geophysics University of Wisconsin - Madison 1215 W. Dayton St., Madison, WI 53706, United States Buchheit, T E (tebuchh@sandia.gov), Sandia National Laboratories, MS 0751, Albuquerque, NM 87123-0751, United States

In recent years, population growth, climate change and surface water contamination have resulted in increased development of deep groundwater resources. As this trend continues, quantitative data of aquifer properties (i.e. specific storage) will be crucial for ensuring the sustainability of groundwater resources. Specific storage is the volume of water expelled due to the compressibility of the surrounding granular skeleton. The compressibility is believed to be controlled by the microscale properties of the surrounding granular skeleton (e.g. grain size, shape and mineralogy) and cementation. In general, cementation is assumed to reduce specific storage by reducing compressibility. However, data which adequately quantifies the relationship between cementation and storage does not exist due to unknown variations in mineralogy, amount and distribution of the cement. The goal of this research is to address this problem by quantifying the effects of cementation of granular porous media on specific storage using discrete-based numerical models. We modeled the modification of specific storage by varying microscale parameters (cement stiffness and percent of cementation) in a series of discrete element models (DEM). Modifications to microscale properties produce bulk-scale responses that are represented in a series of stress/strain curves for each model. The stress/strain curves are used to determine the stiffness for each model and to calculate percent of change in specific storage as a function of the microscale properties. Additionally, we are fabricating analog granular assemblages which will be used to better understand fluid storage in weakly cemented clastic materials. The assemblages will be used to evaluate the differences between natural and synthetic systems and to improve and interpret the DEM models. Preliminary results suggest that cementation has a strong control on the elastic storage properties of cemented granular media. We present results which quantify the proportion of cement needed to effectively reduce the specific storage. We anticipate being able to use the bulk-scale results from the DEM models as an input parameter for continuum-based coupled fluid-solid deformation models. The coupled models will allow us to further investigate elastic and inelastic poroelastic properties, (specific storage and aquifer compaction, respectively), with implications for carbon sequestration during which reactive transport (i.e. precipitation) is assumed to be important for fluid storage.

H23C-1513 

Conditions of Fissuring in a Pumped-Faulted Aquifer System

* HERNANDEZ-MARIN, M (mhmarin@vt.edu), DEPT. OF GEOLOGICAL SCIENCES VIRGINIA TECH, 4044 DERRING HALL, BLACKSBURG, VA 24061, United States BURBEY, T J (tjburbey@vt.edu), DEPT. OF GEOLOGICAL SCIENCES VIRGINIA TECH, 4044 DERRING HALL, BLACKSBURG, VA 24061, United States

Earth fissuring associated with subsidence from groundwater pumping is problematic in many arid-zone heavily pumped basins such as Las Vegas Valley. Long-term pumping at rates considerably greater than the natural recharge rate has stressed the heterogeneous aquifer system resulting in a complex stress-strain regime. A rigorous artificial recharge program coupled with increased surface-water importation has allowed water levels to appreciably recover, which has led to surface rebound in some localities. Nonetheless, new fissures continue to appear, particularly near basin-fill faults that behave as barriers to subsidence bowls. The purpose of this research is to develop a series of computational models to better understand the influence that structure (faults), pumping, and hydrostratigraphy has in the generation and propagation of fissures. The hydrostratigraphy of Las Vegas Valley consists of aquifers, aquitards and a relatively dry vadoze zone that may be as thick as 100m in much of the valley. Quaternary faults are typically depicted as scarps resulting from pre- pumping extensional tectonic events and are probably not responsible for the observed strain. The models developed to simulate the stress-strain and deformation processes in a faulted pumped aquifer-aquitard system of Las Vegas use the ABAQUS CAE (Complete ABAQUS Environment) software system. ABAQUS is a sophisticated engineering industry finite-element modeling package capable of simulating the complex fault- fissure system described here. A brittle failure criteria based on the tensile strength of the materials and the acting stresses (from previous models) are being used to understand how and where fissures are likely to form. , Hypothetical simulations include the role that faults and the vadose zone may play in fissure formation