Rate measurements and detection of gas microseepage to the atmosphere from an enhanced oil recovery/sequestration project,Rangely, Colorado,USA

One of the proposals for large-scale sequestration of fossil fuel-derived CO₂ is deep geologic disposal in depleted oil/gas reservoirs or deep aquifers. Previously published scenarios for this inadequately proven technology have either ignored or dismissed the possibility of vertical migration of gases caused by overpressure. Overpressuring of a reservoir or aquifer will be necessary in order to have acceptable rates for dispersal of injected CO₂. This research describes methodology and the results of measurement of microseepage of CO₂ and CH₄ at a large-scale CO₂-enhanced oil recovery (EOR) operation at Rangely, Colorado, USA. Shallow and deep soil gas concentrations, and direct transport of CO₂ and CH₄ into the atmosphere were measured. The interpretation of the measurements was complemented by both stable and radiogenic isotopic measurements of C. The results have demonstrated an estimated microseepage to the atmosphere of approximately 400 metric tonnes of CH₄/a from the 78 km² area of the Rangely field. Preliminary estimates of deep-sourced CO₂ losses are <3800 tonnes/a, based on stable isotope measurements of soil gases. Several holes up to 10 m deep were drilled on, and off the field for nested gas sampling of composition and stable C isotopic ratios for CO₂ and CH₄. Carbon-14 measurements on CO₂ from these holes indicate that deep-sourced CO₂ microseepage losses were approximately 170 tonnes/a.     

Sensitivity analysis of hydrological parameters in modeling flow and transport in the unsaturated zone of Yucca Mountain, Nevada, USA

The unsaturated fractured volcanic deposits at Yucca Mountain in Nevada, USA, have been intensively investigated as a possible repository site for storing high-level radioactive waste. Field studies at the site have revealed that there exist large variabilities in hydrological parameters over the spatial domain of the mountain. Systematic analyses of hydrological parameters using a site-scale three-dimensional unsaturated zone (UZ) flow model have been undertaken. The main objective of the sensitivity analyses was to evaluate the effects of uncertainties in hydrologic parameters on modeled UZ flow and contaminant transport results. Sensitivity analyses were carried out relative to fracture and matrix permeability and capillary strength (van Genuchten α) through variation of these parameter values by one standard deviation from the base-case values. The parameter variation resulted in eight parameter sets. Modeling results for the eight UZ flow sensitivity cases have been compared with field observed data and simulation results from the base-case model. The effects of parameter uncertainties on the flow fields were evaluated through comparison of results for flow and transport. In general, this study shows that uncertainties in matrix parameters cause larger uncertainty in simulated moisture flux than corresponding uncertainties in fracture properties for unsaturated flow through heterogeneous fractured rock.     

Thermal history of the unsaturated zone at Yucca Mountain,Nevada, USA

Secondary calcite, silica and minor amounts of fluorite deposited in fractures and cavities record the chemistry, temperatures, and timing of past fluid movement in the unsaturated zone at Yucca Mountain, Nevada, the proposed site of a high-level radioactive waste repository. The distribution and geochemistry of these deposits are consistent with low-temperature precipitation from meteoric waters that infiltrated at the surface and percolated down through the unsaturated zone. However, the discovery of fluid inclusions in calcite with homogenization temperatures (T_h) up to ∼80 °C was construed by some scientists as strong evidence for hydrothermal deposition. This paper reports the results of investigations to test the hypothesis of hydrothermal deposition and to determine the temperature and timing of secondary mineral deposition. Mineral precipitation temperatures in the unsaturated zone are estimated from calcite- and fluorite-hosted fluid inclusions and calcite δ¹⁸O values, and depositional timing is constrained by the ²⁰⁷Pb/²³⁵U ages of chalcedony or opal in the deposits. Fluid inclusion T_h from 50 samples of calcite and four samples of fluorite range from ∼35 to ∼90 °C. Calcite δ¹⁸O values range from ∼0 to ∼22‰ (SMOW) but most fall between 12 and 20‰. The highest T_h and the lowest δ¹⁸O values are found in the older calcite. Calcite T_h and δ¹⁸O values indicate that most calcite precipitated from water with δ¹⁸O values between −13 and −7‰, similar to modern meteoric waters.     

The use of tracers to assess leakage from the sequestration of CO2 in a depleted oil reservoir,New Mexico,USA

Geological sequestration of CO₂ in depleted oil reservoirs is a potentially useful strategy for greenhouse gas management and can be combined with enhanced oil recovery. Development of methods to estimate CO₂ leakage rates is essential to assure that storage objectives are being met at sequestration facilities. Perfluorocarbon tracers (PFTs) were added as three 12 h slugs at about one week intervals during the injection of 2090 tons of CO₂ into the West Pearl Queen (WPQ) depleted oil formation, sequestration pilot study site located in SE New Mexico. The CO₂ was injected into the Permian Queen Formation. Leakage was monitored in soil–gas using a matrix of 40 capillary adsorbent tubes (CATs) left in the soil for periods ranging from days to months. The tracers, perfluoro-1,2-dimethylcyclohexane (PDCH), perfluorotrimethylcyclohexane (PTCH) and perfluorodimethylcyclobutane (PDCB), were analyzed using thermal desorption, and gas chromatography with electron capture detection. Monitoring was designed to look for immediate leakage, such as at the injection well bore and at nearby wells, and to develop the technology to estimate overall CO₂ leak rates based on the use of PFTs. Tracers were detected in soil–gas at the monitoring sites 50 m from the injection well within days of injection. Tracers continued to escape over the following years. Leakage appears to have emanated from the vicinity of the injection well in a radial pattern to about 100 m and in directional patterns to 300 m. Leakage rates were estimated for the 3 tracers from each of the 4 sets of CATs in place following the start of CO₂ injection. Leakage was fairly uniform during this period. As a first approximation, the CO₂ leak rate was estimated at about 0.0085% of the total CO₂ sequestered per annum.     

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.     

Isocon analysis of migmatization in the Front Range, Colorado, USA

Isocon analysis has been applied to five sets of leucosome, mafic selvages and immediately adjacent mesosome in the migmatites from a 15-m outcrop in the Colorado Front Range. The results show: (i) mafic selvages formed from the adjacent mesosome by loss of felsic components and therefore the mesosomes are indeed palaeosomes or protoliths; (ii) the leucosomes did not form in a closed system from the palaeosome (in which case the material lost from the palaeosome during selvage formation would become the leucosome). The observed volumes and compositions of leucosomes require that the present leucosome must contain some material in addition to the felsic components lost from the selvages. The materials that must be added are leucotonalitic to granitic in composition, varying greatly in K/(Na + Ca) ratio. The trend in leucosome composition can be reproduced by assuming that a metasomatic exchange, KNa + Ca, modified originally leucotonalitic leucosomes to more K-rich compositions. These leucosomes most likely formed by injection of silicate melts accompanied, or followed, by metasomatism. The trend of leucosome compositions in this study reflects the general trend in the leucosome compositions which have been published from other areas, indicating that the proposed mechanism can be applicable to other regional migmatites.     

Transport of elemental mercury in the unsaturated zone from a waste disposal site in an arid region

Mercury contained in buried landfill waste may be released via upward emission to the atmosphere or downward leaching to groundwater. Data from the US Geological Survey’s Amargosa Desert Research Site (ADRS) in arid southwestern Nevada reveal another potential pathway of Hg release: long-distance (10² m) lateral migration of elemental Hg (Hg⁰) through the unsaturated zone. Gas collected from multiple depths from two instrumented boreholes that sample the entire 110-m unsaturated zone thickness and are located 100 and 160 m away from the closest waste burial trench exhibit gaseous Hg concentrations of up to 33 and 11 ng m⁻³, respectively. The vertical distribution of gaseous Hg in the borehole closest to the disposal site shows distinct subsurface peaks in concentration at depths of 1.5 and 24 m that cannot be explained by radial diffusive transport through a heterogeneous layered unsaturated zone. The inability of current models to explain gaseous Hg distribution at the ADRS highlights the need to advance the understanding of gas-phase contaminant transport in unsaturated zones to attain a comprehensive model of landfill Hg release.     

Movement of water through the thick unsaturated zone underlying Oro Grande and Sheep Creek Washes in the western Mojave Desert,USA

Previous studies indicate that a small quantity of recharge occurs from infiltration of streamflow in intermittent streams in the upper Mojave River basin, in the western Mojave Desert, near Victorville, California. Chloride, tritium, and stable isotope data collected in the unsaturated zone between 1994 and 1998 from boreholes drilled in Oro Grande and Sheep Creek Washes indicate that infiltration of streamflow occurs to depths below the root zone, and presumably to the water table, along much of Oro Grande Wash and near the mountain front along Sheep Creek Wash. Differences in infiltration at sites along each wash are the result of hydrologic variables such as proximity to the mountain front, quantity of streamflow, and texture of the subsurface deposits. Differences in infiltration between the washes are the result of large-scale geomorphic processes. For example, Oro Grande wash is incised into the Victorville fan and infiltration has occurred at approximately the same location over recent geologic time. In contrast, Sheep Creek Wash overlies an active alluvial fan and the stream channel can move across the fan surface through time. Infiltration does not occur to depths below the root zone at control sites outside of the washes. Electronic Publication     

Peridotite xenoliths in silica-rich,potassic latite from the transition zone of the Colorado Plateau in north-central Arizona

Potassic latite in the transition zone of the Colorado Plateau near Chino Valley, Arizona, contains abundant eclogite and amphibolite xenoliths and minor websterite and pyroxenite xenoliths. One unit contains peridotite xenoliths; analyzed samples have mg-ratios of 68 and 71, 58 and 63 wt% SiO₂, and are enriched in potassium and other large ion lithophile (LIL) elements. Rare earth element (REE) patterns are light REE enriched with La greater than 100 times chondritic abundance. The peridotite xenoliths are partly to totally altered, but contain remnant olivine, orthopyroxene, and clinopyroxene; one harzburgite nodule also contains spinel. Mineral compositions from the xenoliths are relatively refractory and similar to those in other spinel peridotite xenoliths from the Colorado Plateau. Geothermometry on olivine-spinel and two-pyroxene pairs indicates equilibration temperatures of less than 800° C for the peridotite nodules. The relatively low temperatures calculated from mineral equilibria are consistent with temperature estimates for other mantle nodules from under the Colorado Plateau.Peridotite xenoliths, mg-ratios, and Ni contents are evidence that the latite magma was derived from mantle peridotite. The potassic nature of the magma probably accounts for its silica-rich composition. The potassic, silica-rich nature of the latite and its enrichment in LREE and other LIL elements are consistent with a source which was metasomatically enriched in these elements either before or during partial melting. The source could have been either spinel or garnet peridotite.     

Geologic assessment and injection design for a pilot CO<Subscript>2</Subscript>-enhanced oil recovery and sequestration demonstration in a heterogeneous oil reservoir: Citronelle Field,Alabama, USA

CO₂ pilot injection studies, with site-specific geologic assessment and engineering reservoir design, can be instrumental for demonstrating both incremental enhanced oil recovery and permanent geologic storage of greenhouse gases. The purpose of this paper is to present the geologic and reservoir analyses in support of a field pilot test that will evaluate the technical and economic feasibility of commercial-scale CO₂-enhanced oil recovery to increase oil recovery and extend the productive life of the Citronelle Oil Field, the largest conventional oil field in Alabama (SE USA). Screening of reservoir depth, oil gravity, reservoir pressure, reservoir temperature, and oil composition indicates that the Cretaceous-age Donovan sand, which has produced more than 169 × 10⁶ bbl in Citronelle Oil Field, is amenable to miscible CO₂ flooding. The project team has selected an 81 ha (200 ac) 5-spot test site with one central gas injector, two producers, and two initially temporarily abandoned production wells that are now in production. Injection is planned in two separate phases, each consisting of 6,804 t (7,500 short tons) of food-grade CO₂. The Citronelle Unit B-19-10 #2 well (Permit No. 3232) is the CO₂ injector for the first injection test. The 14-1 and 16-2 sands of the upper Donovan are the target zones. These sandstone units consist of fine to medium-grained sandstone that is enveloped by variegated mudstone. Both of these sandstone units were selected based on the distribution of perforated zones in the test pattern, production history, and the ability to correlate individual sandstone units in geophysical well logs. The pilot injections will evaluate the applicability of tertiary oil recovery to Citronelle Field and will provide a large volume of information on the pressure response of the reservoirs, the mobility of fluids, time to breakthrough, and CO₂ sweep efficiency. The results of the pilot injections will aid in the formulation of commercial-scale reservoir management strategies that can be applied to Citronelle Field and other geologically heterogeneous oil fields and the design of similar pilot injection projects.     

Mixed-layer illite/smectite as a paleotemperature indicator in the Amethyst vein system,Creede district,Colorado, USA

Ordered illite/smectite is the most abundant and widespread product of argillic alteration associated with the Oligocene, epithermal, Ag-Pb-Zn-Cu-Au Amethyst vein system, southwest Colorado. Hydrothermal illite/smectites exhibit all Reichweite from 1 to 3 and span the composition range from about 25% to about 3% illite.— The composition and Reichweite of illite/smectite vary smoothly with distance from the Amethyst vein. With increasing distance from the vein, percent illite and the Reichweite decrease. When composition-temperature relationships of illite/smectite from modern geothermal systems are applied to the fossil Amethyst system, isotherms describing the thermal regime at the time of argillization can be estimated. Temperatures near 240° C appear to have existed near the Amethyst vein; these temperatures agree with homogenization temperatures obtained from fluid inclusions in vein minerals. The most distal illite/smectites, sampled about 260 m from the vein, indicate temperatures near 110° C. Estimated thermal gradients are on the order of 0.4 to 1° C per meter. — Although illite/smectite composition and structure vary systematically with distance from the Amethyst vein, there are no systematic trends associated with the numerous, smaller veins and veinlets in the hanging wall of the system. This indicates that temperatures of both wallrock and the fluids in all but the major Amethyst vein were nearly the same during clay formation. Apparently, the hydrothermal system had reached a fairly steady, mature, thermal state at the time of argillization.     

Geoarchaeological investigation of lithic resources in the central Mesa Verde region,Colorado, USA

An important question in Southwestern archaeological research is how social and economic ties between prehistoric groups changed through time. The chemical and mineralogical matching of stone tools and ceramics to source areas is one approach to addressing this question. Here we challenge the assumption that igneous rocks used for manufacturing stone tools and ceramics in the Mesa Verde region of Colorado were originally derived from the Chuska Mountains of New Mexico. Recent geological surveys suggest that small igneous features present in and around Mesa Verde National Park and Ute Mountain Ute Tribal Park contain similar lithologies and may represent a local source. We systematically visited seven of these igneous exposures and nine downstream gravel deposits, looking for geological sources and quarry sites matching the lithology of analyzed debitage from this region. Our results suggest that the ancestral Puebloans could have procured such igneous rocks as aphanitic minette and trachybasalt pottery temper from local areas. Because the trachybasalt is usually thought to have been imported from the Chuska Mountains, this discovery may require a revision of interpreted trade networks based on pottery temper and lithic raw material distinctions. © 2009 Wiley Periodicals, Inc.     

Groundwater flow patterns in the San Luis Valley,Colorado, USA revisited: an evaluation of solute and isotopic data

Groundwater systems in the San Luis Valley, Colorado, USA have been re-evaluated by an analysis of solute and isotopic data. Existing stream, spring, and groundwater samples have been augmented with 154 solute and isotopic samples. Based on geochemical stratification, three groundwater regimes have been identified within 1,200 m of the surface: unconfined, upper active confined, and lower active confined with maximum TDS concentrations of 35,000, 3,500 and 600 mg/L, respectively. The elevated TDS of northern valley unconfined and upper active confined systems result from mineral dissolution, ion exchange and methanogenesis of organic and evaporate lake sediments deposited in an ancient lake, herein designated as Lake Sipapu. Chemical evolutions along flow paths were modeled with NETPATH. Groundwater ages, and δ¹³C, δ²H and δ¹⁸O compositions and distributions, suggest that mountain front recharge is the principle recharge mechanism for the upper and lower confined aquifers with travel times in the northern valley of more than 20,000 and 30,000 ¹⁴C years, respectively. Southern valley confined aquifer travel times are 5,000 ¹⁴C years or less. The unconfined aquifer contains appreciable modern recharge water and the contribution of confined aquifer water to the unconfined aquifer does not exceed 20%.     

Estimation of actual evapotranspiration by numerical modeling of water flow in the unsaturated zone: a case study in Buenos Aires, Argentina

A method is presented to estimate actual evapotranspiration (ET_A) from potential evapotranspiration (ET_P) by numerical modeling of water flow in the unsaturated zone. Water flow is described by the Richards equation with a sink term representing the root water uptake. Evaporation is included in the model as a Neumann boundary condition at the soil surface. The Richards equation is solved in a one-dimensional domain using a mixed finite element method. The values of ET_A are obtained by applying a water stress factor to ET_P to account for soil moisture changes during the simulation period. The proposed numerical model is used to estimate ET_A in an experimental plot located in a flatland area in Buenos Aires (Argentina). Numerical results show that the proposed model is a useful tool for evaluating evapotranspiration under different scenarios.     

美国Paradox盆地盐构造新生代的拉伸活化

美国科罗拉多高原Paradox盆地中的裂口盐墙和盐背斜,主要是在古生界和中生界沉积物的差异负载作用下发育起来的。盐构造顶部的平缓褶皱,通常被认为是晚白垩纪至古新世拉腊米挤压运动的产物,而其脊顶地堑与山谷,则被认为是挤压后的松弛拉伸及盐溶作用所造成的。基于野外调查、物理模拟以及现代盐构造学理论,研究认为这些平缓褶皱和脊顶地堑,主要是新生代始新世至渐新世北北东向区域拉伸作用的产物。因为岩盐较其围岩软弱,拉伸形变主要集中在早先形成的盐墙、盐背斜及其较薄的顶板中,使盐构造发生活化并形成新的拉伸构造。顶板中的拉伸构造主要包括正断层、地堑、断层滑移(断滑)褶皱以及滚动褶皱。由于区域拉伸方向斜交盐构造走向,断层分布多呈羽列状。在母盐层较厚地区,拉伸使盐刺穿在正断层之下次动上升;而在母盐层较薄地区,拉伸使盐刺穿下降。在初始盐刺穿较高地区,盐构造的侧翼下降;而在初始盐刺穿较低地区,盐构造整体下降。拉伸使盐刺穿变宽,脊顶地堑下沉,形成拉伸型下降盐刺穿特征性的尖角盐体构造。局部地区,盐体从破碎的顶板间溢出,覆盖了下沉的脊顶地堑。     

Late Pleistocene and Holocene paleoseismology of an intraplate seismic zone in a large alluvial valley, the New Madrid seismic zone, Central USA

The New Madrid seismic zone (NMSZ) is an intraplate right-lateral strike-slip and thrust fault system contained mostly within the Mississippi Alluvial Valley. The most recent earthquake sequence in the zone occurred in 1811–1812 and had estimated moment magnitudes of 7–8 (e.g., [Johnston, A.C., 1996. Seismic moment assessment of stable continental earthquakes, Part 3: 1811–1812 New Madrid, 1886 Charleston, and 1755 Lisbon. Geophysical Journal International 126, 314–344; Johnston, A.C., Schweig III, E.S, 1996. The enigma of the New Madrid earthquakes of 1811–1812. Annual Reviews of Earth and Planetary Sciences 24, 339–384; Hough, S.E., Armbruster, J.G., Seeber, L., Hough, J.F., 2000. On the modified Mercalli intensities and magnitudes of the New Madrid earthquakes. Journal of Geophysical Research 105 (B10), 23,839–23,864; Tuttle, M.P., 2001. The use of liquefaction features in paleoseismology: Lessons learned in the New Madrid seismic zone, central United States. Journal of Seismology 5, 361–380]). Four earlier prehistoric earthquakes or earthquake sequences have been dated A.D. 1450 ± 150, 900 ± 100, 300 ± 200, and 2350 B.C. ± 200 years using paleoliquefaction features, particularly those associated with native American artifacts, and in some cases surface deformation ([Craven, J. A. 1995. Paleoseismology study in the New Madrid seismic zone using geological and archeological features to constrain ages of liquefaction deposits. M.S thesis, University of Memphis, Memphis, TN, U.S.A.; Tuttle, M.P., Lafferty III, R.H., Guccione, M.J., Schweig III, E.S., Lopinot, N., Cande, R., Dyer-Williams, K., Haynes, M., 1996. Use of archaeology to date liquefaction features and seismic events in the New Madrid seismic zone, central United States. Geoarchaeology 11, 451–480; Guccione, M.J., Mueller, K., Champion, J., Shepherd, S., Odhiambo, B., 2002b. Stream response to repeated co-seismic folding, Tiptonville dome, western Tennessee. Geomorphology 43(2002), 313–349; Tuttle, M.P., Schweig, E.S., Sims, J.D., Lafferty, R.H., Wolf, L.W., Haynes, M.L., 2002. The earthquake potential of the New Madrid seismic zone, Bulletin of the Seismological Society of America, v 92, n. 6, p. 2080–2089; Tuttle, M.P., Schweig III, E.S., Campbell, J., Thomas, P.M., Sims, J.D., Lafferty III, R.H., 2005. Evidence for New Madrid earthquakes in A.D. 300 and 2350 B.C. Seismological Research Letters 76, 489–501]). The two most recent prehistoric and the 2350 B.C. events were probably also earthquake sequences with approximately the same magnitude as the historic sequence.Surface deformation (faulting and folding) in an alluvial setting provides many examples of stream response to gradient changes that can also be used to date past earthquake events. Stream responses include changes in channel morphology, deviations in the channel path from the regional gradient, changes in the direction of flow, anomalous longitudinal profiles, and aggradation or incision of the channel ([Merritts, D., Hesterberg, T, 1994. Stream networks and long-term surface uplift in the New Madrid seismic zone. Science 265, 1081–1084.; Guccione, M.J., Mueller, K., Champion, J., Shepherd, S., Odhiambo, B., 2002b. Stream response to repeated co-seismic folding, Tiptonville dome, western Tennessee. Geomorphology 43 (2002), 313–349]). Uplift or depression of the floodplain affects the frequency of flooding and thus the thickness and style of vertical accretion or drowning of a meander scar to form a lake. Vegetation may experience trauma, mortality, and in some cases growth enhancement due to ground failure during the earthquake and hydrologic changes after the earthquake ([VanArdale, R.B., Stahle, D.W., Cleaveland, M.K., Guccione, M.J., 1998. Earthquake signals in tree-ring data from the New Madrid seismic zone and implications for paleoseismicity. Geology 26, 515–518]). Identification and dating these physical and biologic responses allows source areas to be identified and seismic events to be dated.Seven fault segments are recognized by microseismicity and geomorphology. Surface faulting has been recognized at three of these segments, Reelfoot fault, New Madrid North fault, and Bootheel fault. The Reelfoot fault is a compressive stepover along the strike-slip fault and has up to 11 m of surface relief ([Carlson, S.D., 2000. Formation and geomorphic history of Reelfoot Lake: insight into the New Madrid seismic zone. M.S. Thesis, University of Arkansas, Fayetteville, Arkansas, U.S.A]) deforming abandoned and active Mississippi River channels ([Guccione, M.J., Mueller, K., Champion, J., Shepherd, S., Odhiambo, B., 2002b. Stream response to repeated co-seismic folding, Tiptonville dome, western Tennessee. Geomorphology 43 (2002), 313–349]). The New Madrid North fault apparently has only strike-slip motion and is recognized by modern microseismicity, geomorphic anomalies, and sand cataclasis ([Baldwin, J.N., Barron A.D., Kelson, K.I., Harris, J.B., Cashman, S., 2002. Preliminary paleoseismic and geophysical investigation of the North Farrenburg lineament: primary tectonic deformation associated with the New Madrid North Fault?. Seismological Research Letters 73, 393–413]). The Bootheel fault, which is not identified by the modern microseismicity, is associated with extensive liquefaction and offset channels ([Guccione, M.J., Marple, R., Autin, W.J., 2005, Evidence for Holocene displacements on the Bootheel fault (lineament) in southeastern Missouri: Seismotectonic implications for the New Madrid region. Geological Society of America Bulletin 117, 319–333]). The fault has dominantly strike-slip motion but also has a vertical component of slip. Other recognized surface deformation includes relatively low-relief folding at Big Lake/Manila high ([Guccione, M.J., VanArdale, R.B., Hehr, L.H., 2000. Origin and age of the Manila high and associated Big Lake “Sunklands”, New Madrid seismic zone, northeastern Arkansas. Geological Society of America Bulletin 112, 579–590]) and Lake St. Francis/Marked Tree high ([Guccione, M.J., VanArsdale, R.B., 1995. Origin and age of the St. Francis Sunklands using drainage patterns and sedimentology. Final report submitted to the U. S. Geological Survey, Award Number 1434-93-G-2354, Washington D.C.]), both along the subsurface Blytheville arch. Deformation at each of the fault segments does not occur during each earthquake event, indicating that earthquake sources have varied throughout the Holocene.     

The noble gas geochemistry of natural CO2 gas reservoirs from the Colorado Plateau and Rocky Mountain provinces, USA

Identification of the source of CO₂ in natural reservoirs and development of physical models to account for the migration and interaction of this CO₂ with the groundwater is essential for developing a quantitative understanding of the long term storage potential of CO₂ in the subsurface. We present the results of 57 noble gas determinations in CO₂ rich fields (>82%) from three natural reservoirs to the east of the Colorado Plateau uplift province, USA (Bravo Dome, NM., Sheep Mountain, CO. and McCallum Dome, CO.), and from two reservoirs from within the uplift area (St. John’s Dome, AZ., and McElmo Dome, CO.). We demonstrate that all fields have CO₂/³He ratios consistent with a dominantly magmatic source. The most recent volcanics in the province date from 8 to 10 ka and are associated with the Bravo Dome field. The oldest magmatic activity dates from 42 to 70 Ma and is associated with the McElmo Dome field, located in the tectonically stable centre of the Colorado Plateau: CO₂ can be stored within the subsurface on a millennia timescale.The manner and extent of contact of the CO₂ phase with the groundwater system is a critical parameter in using these systems as natural analogues for geological storage of anthropogenic CO₂. We show that coherent fractionation of groundwater ²⁰Ne/³⁶Ar with crustal radiogenic noble gases (⁴He, ²¹Ne, ⁴⁰Ar) is explained by a two stage re-dissolution model: Stage 1: Magmatic CO₂ injection into the groundwater system strips dissolved air-derived noble gases (ASW) and accumulated crustal/radiogenic noble gas by CO₂/water phase partitioning. The CO₂ containing the groundwater stripped gases provides the first reservoir fluid charge. Subsequent charges of CO₂ provide no more ASW or crustal noble gases, and serve only to dilute the original ASW and crustal noble gas rich CO₂. Reservoir scale preservation of concentration gradients in ASW-derived noble gases thus provide CO₂ filling direction. This is seen in the Bravo Dome and St. John’s Dome fields. Stage 2: The noble gases re-dissolve into any available gas stripped groundwater. This is modeled as a Rayleigh distillation process and enables us to quantify for each sample: (1) the volume of groundwater originally ‘stripped’ on reservoir filling; and (2) the volume of groundwater involved in subsequent interaction. The original water volume that is gas stripped varies from as low as 0.0005 cm³ groundwater/cm³ gas (STP) in one Bravo Dome sample, to 2.56 cm³ groundwater/cm³ gas (STP) in a St. John’s Dome sample. Subsequent gas/groundwater equilibration varies within all fields, each showing a similar range, from zero to ∼100 cm³ water/cm³ gas (at reservoir pressure and temperature).