Temporal change in groundwater level following the 1999 (Mw = 7.5) Chi-Chi earthquake, Taiwan

We examine the post‐seismic change in the groundwater level following the 1999 (M_w = 7.5) Chi‐Chi earthquake in central Taiwan, as recorded by a network of 70 evenly distributed hydrological stations over a large alluvial fan near the epicenter. Four types of post‐seismic responses may be distinguished. In type 1, the groundwater level declined exponentially with time following a coseismic rise. This was the most common response in the study area and occurred in unconsolidated sediments on the Choshui River fan. In type 2, the groundwater level rose exponentially with time following a coseismic fall. This occurred in the deformed and fractured sedimentary rocks in the foothills near the Chelungpu fault that ruptured in the Chi‐Chi earthquake. In type 3, the groundwater level continued to decline with time following a coseismic fall. This also occurred in the deformed and fractured sedimentary rocks near the ruptured fault. Finally, in type 4, the groundwater level, following a coseismic rise, stayed at the same level or even rose with time before it eventually declined. This occurred mostly in unconsolidated sediments along the coast of central Taiwan and along the Peikang Stream. We analyze these post‐seismic responses by using a one‐dimensional model. Together with the results from well test, the analysis show that the type 1 response may be explained by an aquifer model with coseismic recharge and post‐seismic subhorizontal discharge across a length of 500–5000 m; the type 2 response may be explained by a model of coseismic discharge and post‐seismic recharge from surface water; the type 3 response may be explained by a model of coseismic discharge and post‐seismic subhorizontal discharge across a length of 500–5000 m; and the type 4 response may be explained by a model of coseismic recharge and sustained post‐seismic recharge from surface water. The characteristic time for the post‐seismic changes is similar to that for the groundwater‐level decline during dry seasons before the earthquake, suggesting that there was no earthquake‐induced changes in the aquifer properties (i.e. hydraulic conductivity), confirming the earlier results from recession analyses of the post‐seismic streamflow elsewhere after several earthquakes.     

The timescale and mechanisms of fault sealing and water‐rock interaction after an earthquake

Hydrogeochemical monitoring of a basalt‐hosted aquifer, which contains Ice Age meteoric water and is situated at 1220 m below sea level in the Tjörnes Fracture Zone, northern Iceland, has been ongoing since July 2002. Based on hydrogeochemical changes following an earthquake of magnitude (M_w) 5.8 on 16 September 2002, we constrained the timescales of post‐seismic fault sealing and water–rock interaction. We interpret that the earthquake ruptured a hydrological barrier, permitting a rapid influx of chemically and isotopically distinct Ice Age meteoric water from a second aquifer. During the two subsequent years, we monitored a chemical and isotopic recovery towards pre‐earthquake aquifer compositions, which we interpret to have been mainly facilitated by fault‐sealing processes. This recovery was interrupted in November 2004 by a second rupturing event, which was probably induced by two minor earthquakes and which reopened the pathway to the second aquifer. We conclude that the timescale of fault sealing was approximately 2 years and that the approach to isotopic equilibrium (from global meteoric water line) was approximately 18% after >10⁴ years.     

Equations of state for basin geofluids: algorithm review and intercomparison for brines

Physical properties of formation waters in sedimentary basins can vary by more than 25% for density and by one order of magnitude for viscosity. Density differences may enhance or retard flow driven by other mechanisms and can initiate buoyancy‐driven flow. For a given driving force, the flow rate and injectivity depend on viscosity and permeability. Thus, variations in the density and viscosity of formation waters may have or had a significant effect on the flow pattern in a sedimentary basin, with consequences for various basin processes. Therefore, it is critical to correctly estimate water properties at formation conditions for proper representation and interpretation of present flow systems, and for numerical simulations of basin evolution, hydrocarbon migration, ore genesis, and fate of injected fluids in sedimentary basins. Algorithms published over the years to calculate water density and viscosity as a function of temperature, pressure and salinity are based on empirical fitting of laboratory‐measured properties of predominantly NaCl solutions, but also field brines. A review and comparison of various algorithms are presented here, both in terms of applicability range and estimates of density and viscosity. The paucity of measured formation‐water properties at in situ conditions hinders a definitive conclusion regarding the validity of any of these algorithms. However, the comparison indicates the versatility of the various algorithms in various ranges of conditions found in sedimentary basins. The applicability of these algorithms to the density of formation waters in the Alberta Basin is also examined using a high‐quality database of 4854 water analyses. Consideration is also given to the percentage of cations that are heavier than Na in the waters.     

Response of Alum Rock springs to the October 30, 2007 Alum Rock earthquake and implications for the origin of increased discharge after earthquakes

The origins of increased stream flow and spring discharge following earthquakes have been the subject of controversy, in large part because there are many models to explain observations and few measurements suitable for distinguishing between hypotheses. On October 30, 2007 a magnitude 5.5 earthquake occurred near the Alum Rock springs, California, USA. Within a day we documented a several‐fold increase in discharge. Over the following year, we have monitored a gradual return towards pre‐earthquake properties, but for the largest springs there appears to be a permanent increase in discharge. The Alum Rock springs discharge waters that are a mixture between modern (shallow) meteoric water and old (deep) connate waters expelled by regional transpression. After the earthquake, there was a small and temporary decrease in the fraction of connate water in the largest springs. Accompanying this geochemical change was a small (1–2°C) temperature decrease. Combined with the rapid response, this implies that the increased discharge has a shallow origin. Increased discharge at these springs occurs both for earthquakes that cause static volumetric expansion and for those that cause contraction, supporting models in which dynamic strains are responsible for the subsurface changes that cause flow to increase. We make a quantitative comparison between the observed changes and model predictions for three types of models: (i) a permanent increase in permeability; (ii) an increase in permeability followed by a gradual decrease to its pre‐earthquake value; and (iii) an increase of hydraulic head in the groundwater system discharging at the springs. We show that models in which the permeability of the fracture system feeding the springs increases after the earthquake are in general consistent with the changes in discharge. The postseismic decrease in discharge could either reflect the groundwater system adjusting to the new, higher permeability or a gradual return of permeability to pre‐earthquake values; the available data do not allow us to distinguish between these two scenarios. However, the response of these springs to another earthquake will provide critical constraints on the changes that occur in the subsurface and should permit a test of all three types of models.     

Permeability evolution in sorbing media: analogies between organic‐rich shale and coal

Shale gas reservoirs like coalbed methane (CBM) reservoirs are promising targets for geological sequestration of carbon dioxide (CO₂). However, the evolution of permeability in shale reservoirs on injection of CO₂ is poorly understood unlike CBM reservoirs. In this study, we report measurements of permeability evolution in shales infiltrated separately by nonsorbing (He) and sorbing (CO₂) gases under varying gas pressures and confining stresses. Experiments are completed on Pennsylvanian shales containing both natural and artificial fractures under nonpropped and propped conditions. We use the models for permeability evolution in coal (Journal of Petroleum Science and Engineering, Under Revision) to codify the permeability evolution observed in the shale samples. It is observed that for a naturally fractured shale, the He permeability increases by approximately 15% as effective stress is reduced by increasing the gas pressure from 1 MPa to 6 MPa at constant confining stress of 10 MPa. Conversely, the CO₂ permeability reduces by a factor of two under similar conditions. A second core is split with a fine saw to create a smooth artificial fracture and the permeabilities are measured for both nonpropped and propped fractures. The He permeability of a propped artificial fracture is approximately 2‐ to 3fold that of the nonpropped fracture. The He permeability increases with gas pressure under constant confining stress for both nonpropped and propped cases. However, the CO₂ permeability of the propped fracture decreases by between one‐half to one‐third as the gas pressure increases from 1 to 4 MPa at constant confining stress. Interestingly, the CO₂ permeability of nonpropped fracture increases with gas pressure at constant confining stress. The permeability evolution of nonpropped and propped artificial fractures in shale is found to be similar to those observed in coals but the extent of permeability reduction by swelling is much lower in shale due to its lower organic content. Optical profilometry is used to quantify the surface roughness. The changes in surface roughness indicate significant influence of proppant indentation on fracture surface in the shale sample. The trends of permeability evolution on injection of CO₂ in coals and shales are found analogous; therefore, the permeability evolution models previously developed for coals are adopted to explain the permeability evolution in shales.     

Indications of formation water flow and compartmentalization on the flank of a salt structure derived from salinity and seismic data

Vertical and lateral variations in lithology, salinity, temperature, and pressure determined from wireline LAS logs, produced water samples, and seismic data on the south flank of a salt structure on the continental shelf, offshore Louisiana indicate three hydrogeologic zones in the study area: a shallow region from 0 to 1.1 km depth with hydrostatically pressured, shale‐dominated Pleistocene age sediments containing pore waters with sea water (35 g l^?1) or slightly above sea water salinity; a middle region from 1.1 to 3.2 km depth with near hydrostatically pressured, sand‐dominated Pliocene age sediments that contain pore waters that range from seawater salinity to up to 5 times sea water salinity (180 g l^?1); and a deep section below 3.2 km depth with geopressured, shale‐dominated Miocene age sediments containing pore waters that range from sea water salinity to 125 g l^?1. Salt dissolution has generated dense, saline waters that appear to be migrating down dip preferentially through the thick Pliocene sandy section. Sand layers that come in contact with salt contain pore waters with high salinity. Isolated sands have near sea water salinity. Salinity information in conjunction with seismic data is used to infer fluid compartmentalization. Both vertical and lateral lithologic barriers to fluid flow at tens to hundreds of meters scale are observed. Fluid compartmentalization is also evident across a supradomal normal fault. Offset of salinity contours are consistent with the throw of the fault, which suggests that saline fluids migrated before fault formation.     

Carbon dioxide controlled earthquake distribution pattern in the NW Bohemian swarm earthquake region,western Eger Rift,Czech Republic – gas migration in the crystalline basement

The ascent of magmatic carbon dioxide in the western Eger (Oh?e) Rift is interlinked with the fault systems of the Variscian basement. In the Cheb Basin, the minimum CO₂ flux is about 160 m³ h^?1, with a diminishing trend towards the north and ceasing in the main epicentral area of the Northwest Bohemian swarm earthquakes. The ascending CO₂ forms Ca‐Mg‐HCO₃ type waters by leaching of cations from the fault planes and creates clay minerals, such as kaolinite, as alteration products on affected fault planes. These mineral reactions result in fault weakness and in hydraulically interconnected fault network. This leads to a decrease in the friction coefficient of the Coulomb failure stress (CFS) and to fault creep as stress build‐up cannot occur in the weak segments. At the transition zone in the north of the Cheb Basin, between areas of weak, fluid conductive faults and areas of locked faults with frictional strength, fluid pressure can increase resulting in stress build‐up. This can trigger strike‐slip swarm earthquakes. Fault creep or movements in weak segments may support a stress build‐up in the transition area by transmitting fluid pressure pulses. Additionally to fluid‐driven triggering models, it is important to consider that fluids ascending along faults are CO₂‐supersaturated thus intensifying the effect of fluid flow. The enforced flow of CO₂‐supersaturated fluids in the transitional zone from high to low permeability segments through narrowings triggers gas exsolution and may generate pressure fluctuations. Phase separation starts according to the phase behaviour of CO₂‐H₂O systems in the seismically active depths of NW Bohemia and may explain the vertical distribution of the seismicity. Changes in the size of the fluid transport channels in the fault systems caused, or superimposed, by fault movements, can produce fluid pressure increases or pulses, which are the precondition for triggering fluid‐induced swarm earthquakes.     

Geodynamically induced variations in the emission of CO2 gas at San Faustino (Central Apennines,Italy)

The Central Apennines are affected by frequent earthquakes of moderate magnitude that occur mainly within the upper part of the crust at depths of <15 km. A large number of cold gas emissions that are rich in CO₂ are also found in the region. One particular vent with a high rate of degassing was equipped with a sensor to measure flow rates, which were recorded for a number of different periods between 2005 and 2010. Factors that could affect potentially CO₂ flow rates include barometric pressure, atmospheric temperature, precipitation and local seismicity. Our analysis indicates that the periods of anomalous flow rate were related not to the environmental factors but probably to the deformative processes of the crust associated with the local seismicity. Local seismic events as expression of geodynamic processes occurred always before and during these anomalous gas flow periods. This correlation exists only for events that occurred eastwards of the gas emission site close to the Martana fault zone. We herein consider this correlation as indication for a continuous interaction between the field of static strain and the deep fluid pressure. An approximation of the fluid pressure transmission towards the gas emission site gives reasonable values of 1–10 m² sec^?1. To make comparisons with the long‐term effects of the static strain, we also recorded the short‐term effects of the dynamic release of strain induced by the series of strong earthquakes that took place in L’Aquila in 2009. We detected a significant anomalous flow rate that occurred at the same time as this seismic sequence, during which widespread degassing was induced around the focal zone.     

Transient fluid flow in and around a fault

We present the results of simple numerical experiments in which we study the evolution with time of fluid flow around and within a permeable fault embedded in a less permeable porous medium. Fluid movement is driven by an imposed vertical pressure gradient. The results show that fluid flow is controlled by two timescales: τ_f = Sl²/κ_F and τ_F = Sl²/κ_M, where S is the specific storage of the porous material, l the length of the fault, and κ_M and κ_F are the hydraulic conductivities of the porous material and the fault, respectively. Fluid flow and the associated fluid pressure field evolve through three temporal stages: an early phase [t < τ_f] during which the initial fluid pressure gradient within the fault is relaxed; a second transient stage [τ_f < t < τ_F] when fluid is rapidly expelled at one end of the fault and extracted from the surrounding rocks at the other end leading to a reduction in the pressure gradient in the intact rock; a third phase [t < τ_F] characterized by a steady‐state flow. From the numerical experiments we derive an expression for the steady‐state maximum fluid velocity in the fault and the values of the two timescales, τ_f and τ_F. A comparison indicates excellent agreement of our results with existing asymptotic solutions. For km‐scale faults, the model results suggest that steady‐state is unlikely to be reached over geological timescales. Thus, the current use of parameters such as the focusing ratio defined under the assumption of steady‐state conditions should be reconsidered.     

Analysis of strain‐induced ground‐water fluctuations at Devils Hole,Nevada

The pool in Devils Hole is a sensitive indicator of crustal strain and fluctuates in response to changes in atmospheric pressure, earth tides, earthquakes, large‐scale tectonic activity and ground‐water development. Short‐term and cyclic water‐level fluctuations caused by atmospheric pressure and earth tides were found to be on the order of millimeters to centimeters. The 1992 Landers/Little Skull Mountain earthquake sequence and the 1999 Hector Mine earthquake induced water‐level offsets of greater than ?12 and ?3.6 cm, respectively. The results of a dislocation model used to compute volumetric strain for each earthquake indicates that the coseismic water‐level offsets are consistent in magnitude and sense with poroelastic responses to earthquake‐induced strain. Theoretical postseismic fluid‐flow modeling indicates that the diffusivity of the system is on the order of 0.03 m² sec^?1, and identified areas of anomalous water‐level fluctuations. Interpretation of model results suggests that while the persistent post‐Landers rise in water‐level can be attributed to deformation‐induced channeling of fluid to the Devils Hole fault zone, the cause of the pre‐Hector Mine water‐level rise may be related to postseismic excess fluid pressures or preseismic strain accumulation.     

Stress‐induced temperature variations in groundwater of the Monferrato area (north‐western Italy)

Although characterized by low seismicity, the Monferrato area of north‐western Italy was affected by earthquakes, of magnitude M5.1 and M4.8, in 2000 and 2001. At the same time, marked changes were recorded in water temperature and chemistry in several wells within the epicentral area. In May 2004, an automatic network for the continuous monitoring of groundwater was installed in selected wells to study the phenomenon. Here, we report on data collected during a 3‐year period of groundwater monitoring. During the first year, episodes of water heating (by up to 20°C) were observed in one monitored well. The temporal analysis of the seismic activity recorded in the area revealed as almost all seismic events occurred during the period of elevated water temperatures. The similar timing of earthquakes and groundwater‐temperature anomalies suggests that both may be triggered by the same processes acting in the crust.     

Fracture‐focused fluid flow in an acid and redox‐influenced system: diagenetic controls on cement mineralogy and geomorphology in the Navajo Sandstone

Differential cement mineralogy is influenced by depositional textures, structural deformation, pore fluid chemistry, and ultimately influences landscape evolution by introducing heterogeneities in erodibility. In Southern Utah, the region West of the Kaibab uplift known as Mollies Nipple (Mollies) in Grand Staircase‐Escalante National Monument exhibits a complex history of fluid–sediment interactions, which has resulted in a localized zone of anomalous diagenetic iron sulfate (jarosite) mineralogy in a well‐exposed dune–interdune deposit within the Navajo Sandstone. Mineralogy and geochemistry of cements within this region are examined using reflectance and imaging spectroscopy, field investigations, microscopy, and whole‐rock geochemical analyses. These data show that the in‐situ jarosite cement is localized to a plane along the highest ridge of the butte, providing an armor along with other secondary cements, which controls the butte's geomorphic evolution. The jarosite cement is associated with other mineralogies suggesting that the sulfate was one of the latest fluid‐related precipitates in the paragenetic sequence. It was preceded by a regional bleaching event, precipitation of clay cements, some localized concretionary iron oxide precipitation, and formation of deformation bands. At least one generation of dense iron oxide mineralization is associated with cataclastic brittle deformation predating the sulfate precipitation. Trace element geochemistry of cements shows certain metal oxide populations associated with extremely high (>2000 ppm) arsenic values. We interpret the combination of spatial mineral distribution, observed paragenetic sequence, and trace element geochemistry to suggest this region experienced acid sulfate diagenesis along fracture‐controlled fluid conduits related to weathering of proximal, unidentified, sulfides, or H₂S associated with deep source beds. Jarosite is highly soluble, and its presence suggests that abundant fluid flow has not occurred in this region since its formation. This terminal cement‐forming event must have occurred prior to sandstone exhumation and erosion to form the current extreme landscape at Mollies. This site exhibits the influence that fluid geochemistry, sedimentary mineralogy, and structural fabric have on geomorphic evolution.     

Geophysical signatures associated with fluid flow and gas hydrate occurrence in a tectonically quiescent sequence,Qiongdongnan Basin,South China Sea

Interpretation of high‐resolution two‐dimensional (2D) and three‐dimensional (3D) seismic data collected in the Qiongdongnan Basin, South China Sea reveals the presence of polygonal faults, pockmarks, gas chimneys and slope failure in strata of Pliocene and younger age. The gas chimneys are characterized by low‐amplitude reflections, acoustic turbidity and low P‐wave velocity indicating fluid expulsion pathways. Coherence time slices show that the polygonal faults are restricted to sediments with moderate‐amplitude, continuous reflections. Gas hydrates are identified in seismic data by the presence of bottom simulating reflectors (BSRs), which have high amplitude, reverse polarity and are subparallel to seafloor. Mud diapirism and mounded structures have variable geometry and a great diversity regarding the origin of the fluid and the parent beds. The gas chimneys, mud diapirism, polygonal faults and a seismic facies‐change facilitate the upward migration of thermogenic fluids from underlying sediments. Fluids can be temporarily trapped below the gas hydrate stability zone, but fluid advection may cause gas hydrate dissociation and affect the thickness of gas hydrate zone. The fluid accumulation leads to the generation of excess pore fluids that release along faults, forming pockmarks and mud volcanoes on the seafloor. These features are indicators of fluid flow in a tectonically‐quiescent sequence, Qiongdongnan Basin. Geofluids (2010) 10 , 351–368     

Simulation of the impact of faults on CO2 injection into sandstone reservoirs

Numerical simulations of multiphase CO₂ behavior within faulted sandstone reservoirs examine the impact of fractures and faults on CO₂ migration in potential subsurface injection systems. In southeastern Utah, some natural CO₂ reservoirs are breached and CO₂‐charged water flows to the surface along permeable damage zones adjacent to faults; in other sites, faulted sandstones form barriers to flow and large CO₂‐filled reservoirs result. These end‐members serve as the guides for our modeling, both at sites where nature offers ‘successful’ storage and at sites where leakage has occurred. We consider two end‐member fault types: low‐permeability faults dominated by deformation‐band networks and high‐permeability faults dominated by fracture networks in damage zones adjacent to clay‐rich gouge. Equivalent permeability (k) values for the fault zones can range from <10^?14 m² for deformation‐band‐dominated faults to >10^?12 m² for fracture‐dominated faults regardless of the permeability of unfaulted sandstone. Water–CO₂ fluid‐flow simulations model the injection of CO₂ into high‐k sandstone (5 × 10^?13 m²) with low‐k (5 × 10^?17 m²) or high‐k (5 × 10^?12 m²) fault zones that correspond to deformation‐band‐ or fracture‐dominated faults, respectively. After 500 days, CO₂ rises to produce an inverted cone of free and dissolved CO₂ that spreads laterally away from the injection well. Free CO₂ fills no more than 41% of the pore space behind the advancing CO₂ front, where dissolved CO₂ is at or near geochemical saturation. The low‐k fault zone exerts the greatest impact on the shape of the advancing CO₂ front and restricts the bulk of the dissolved and free CO₂ to the region upstream of the fault barrier. In the high‐k aquifer, the high‐k fault zone exerts a small influence on the shape of the advancing CO₂ front. We also model stacked reservoir seal pairs, and the fracture‐dominated fault acts as a vertical bypass, allowing upward movement of CO₂ into overlying strata. High‐permeability fault zones are important pathways for CO₂ to bypass unfaulted sandstone, which leads to reduce sequestration efficiency. Aquifer compartmentalization by low‐permeability fault barriers leads to improved storativity because the barriers restrict lateral CO₂ migration and maximize the volume and pressure of CO₂ that might be emplaced in each fault‐bound compartment. As much as a 3.5‐MPa pressure increase may develop in the injected reservoir in this model domain, which under certain conditions may lead to pressures close to the fracture pressure of the top seal.     

Gas breakthrough experiments on pelitic rocks: comparative study with N2, CO2 and CH4

The capillary‐sealing efficiency of intermediate‐ to low‐permeable sedimentary rocks has been investigated by N₂, CO₂ and CH₄ breakthrough experiments on initially fully water‐saturated rocks of different lithological compositions. Differential gas pressures up to 20 MPa were imposed across samples of 10–20 mm thickness, and the decline of the differential pressures was monitored over time. Absolute (single‐phase) permeability coefficients (k_abs), determined by steady‐state fluid flow tests, ranged between 10^?22 and 10^?15 m². Maximum effective permeabilities to the gas phase k_eff(max), measured after gas breakthrough at maximum gas saturation, extended from 10^?26 to 10^?18 m². Because of re‐imbibition of water into the interconnected gas‐conducting pore system, the effective permeability to the gas phase decreases with decreasing differential (capillary) pressure. At the end of the breakthrough experiments, a residual pressure difference persists, indicating the shut‐off of the gas‐conducting pore system. These pressures, referred to as the ‘minimum capillary displacement pressures’ (P_d), ranged from 0.1 up to 6.7 MPa. Correlations were established between (i) absolute and effective permeability coefficients and (ii) effective or absolute permeability and capillary displacement pressure. Results indicate systematic differences in gas breakthrough behaviour of N₂, CO₂ and CH₄, reflecting differences in wettability and interfacial tension. Additionally, a simple dynamic model for gas leakage through a capillary seal is presented, taking into account the variation of effective permeability as a function of buoyancy pressure exerted by a gas column underneath the seal.     

Seismic constraints on the effects of gas hydrate on sediment physical properties and fluid flow: a review

The formation of gas hydrates in marine sediments changes their physical properties and hence influences fluid flow. Here, we review seismic indicators of gas hydrates and relate these indicators to gas hydrate formation and fluid migration. Analyses of seismic data from sediments containing gas and gas hydrates in a variety of locations have shown that the characteristic bottom‐simulating reflector (BSR), which commonly marks the hydrate phase boundary is caused mainly by the presence of gas beneath the gas hydrate stability zone (GHSZ). The amplitude of the BSR is also dependent on the hydrate concentration and on the porosity of the sediment. The presence of gas hydrate alters the elastic properties of sediments, particularly if it cements sediment grains. However, multifrequency studies in various geological provinces show that any loss of reflectivity or blanking observed within the GHSZ is dependent on both the nature of the sediments and concentration of hydrate present. Gas beneath the BSR may cause amplitude anomalies and may result in bright spots and enhanced reflections. The presence of gas beneath the BSR is the primary cause of observed amplitude versus offset (AVO) anomalies, but the amplitude of these anomalies is also dependent on the amount of cementation brought by the gas hydrates within the GHSZ. Fluid migration appears to play an important role in the formation and dissociation of gas hydrates in both active and passive margin settings. Fluid migration in accretionary prisms influences hydrate accumulation and may therefore control the spatial distribution of BSRs. Fluid migration may influence also the type of hydrate formed by bringing thermogenic gas containing higher order hydrocarbons to the GHSZ from below. Fluid advection may cause local dissociation of gas hydrates by bringing heat from below, thus shifting the gas hydrate phase boundary. Fluid flow within the GHSZ is limited by the formation of hydrate in the pore space, which reduces the permeability of the sediment. Features such as pockmarks, acoustic masking and acoustic turbidity are indirect indicators of fluid flow and identification of these features in seismic sections within and beneath the GHSZ may also suggest the formation of gas hydrate.     

Multiphase flow simulation through porous media with explicitly resolved fractures

Accurate simulation of multiphase flow in fractured porous media remains a challenge. An important problem is the representation of the discontinuous or near discontinuous behaviour of saturation in real geological formations. In the classical continuum approach, a refined mesh is required at the interface between fracture and porous media to capture the steep gradients in saturation and saturation‐dependent transport properties. This dramatically increases the computational load when large numbers of fractures are present in the numerical model. A discontinuous finite element method is reported here to model flow in fractured porous media. The governing multiphase porous media flow equations are solved in the adaptive mesh computational fluid dynamics code IC‐FERST on unstructured meshes. The method is based on a mixed control volume – discontinuous finite element formulation. This is combined with the P_N+1DG‐P_NDG element pair, which has discontinuous (order N+1) representation for velocity and discontinuous (order N) representation for pressure. A number of test cases are used to evaluate the method's ability to model fracture flow. The first is used to verify the performance of the element pair on structured and unstructured meshes of different resolution. Multiphase flow is then modelled in a range of idealised and simple fracture patterns. Solutions with sharp saturation fronts and computational economy in terms of mesh size are illustrated.     

Geodynamic processes in the NW Bohemian swarm earthquake region,Czech Republic,identified by continuous gas monitoring

The time series of two continuously operating gas monitoring stations at Oldřišská and Nový Kostel located along seismoactive faults in the epicentral area of the NW Bohemian swarm earthquakes (Czech Republic) are compared with water level fluctuations in two boreholes positioned along these faults and with gas flux variations of a mofette at the Soos mofette field at 9 km distance. The seasonal trend of the monitored CO₂ concentration with a maximum in November and a minimum in March/April is governed by groundwater temperatures, superimposed in spring by soil temperatures. CO₂ concentration variations identified at Oldřišská are also reflected in gas flux variations in the Soos mofette and/or water level fluctuations of two boreholes. Variations in the gas monitoring recordings of station at Nový Kostel are also linked with variations at Oldřišská. In all data sets, diurnal variations generated by earth tides occur, reflecting a daily stress – fault permeability cycle. Additional stress interferes with this cycle. Significant, abrupt changes are attributed to geodynamic processes linked with seismic events, as revealed by local seismicity or by the transient of waves of a strong remote earthquake. Simultaneous variations of the gas concentrations in the Nový Kostel area and in the gas flux in the Soos point to an interconnected hydraulic conductive fault systems present in the northern part of the Cheb Basin. Sharp falls in gas concentration, during or subsequent to, earthquake swarms may reflect fault compression associated with impeded gas migration. However, gas variations also occur in periods without seismic activity, indicating changes in fault permeability were caused by local aseismic fault movements, as revealed by events with opposite trends in the gas recordings at Oldřišská, Nový Kostel and the Soos. Therefore, a mathematical approach to establish a correlation between seismicity and gas geochemical variations is not possible.     

Syntectonic infiltration by meteoric waters along the Sevier thrust front,southwest Montana

Structural, petrographic, and isotopic data for calcite veins and carbonate host‐rocks from the Sevier thrust front of SW Montana record syntectonic infiltration by H₂O‐rich fluids with meteoric oxygen isotope compositions. Multiple generations of calcite veins record protracted fluid flow associated with regional Cretaceous contraction and subsequent Eocene extension. Vein mineralization occurred during single and multiple mineralization events, at times under elevated fluid pressures. Low salinity (T_m = ?0.6°C to +3.6°C, as NaCl equivalent salinities) and low temperature (estimated 50–80°C for Cretaceous veins, 60–80°C for Eocene veins) fluids interacted with wall‐rock carbonates at shallow depths (3–4 km in the Cretaceous, 2–3 km in the Eocene) during deformation. Shear and extensional veins of all ages show significant intra‐ and inter‐vein variation in δ¹⁸O and δ¹³C. Carbonate host‐rocks have a mean δ¹⁸O_V‐SMOW value of +22.2 ± 3‰ (1σ), and both the Cretaceous veins and Eocene veins have δ¹⁸O ranging from values similar to those of the host‐rocks to as low as +5 to +6‰. The variation in vein δ¹³C_V‐PDB of ?1 to approximately +6‰ is attributed to original stratigraphic variation and C isotope exchange with hydrocarbons. Using the estimated temperature ranges for vein formation, fluid (as H₂O) δ¹⁸O calculated from Cretaceous vein compositions for the Tendoy and Four Eyes Canyon thrust sheets are ?18.5 to ?12.5‰. For the Eocene veins within the Four Eyes Canyon thrust sheet, calculated H₂O δ¹⁸O values are ?16.3 to ?13.5‰. Fluid–rock exchange was localized along fractures and was likely coincident with hydrocarbon migration. Paleotemperature determinations and stable isotope data for veins are consistent with the infiltration of the foreland thrust sheets by meteoric waters, throughout both Sevier orogenesis and subsequent orogenic collapse. The cessation of the Sevier orogeny was coincident with an evolving paleogeographic landscape associated with the retreat of the Western Interior Seaway and the emergence of the thrust front and foreland basin. Meteoric waters penetrated the foreland carbonate thrust sheets of the Sevier orogeny utilizing an evolving mesoscopic fracture network, which was kinematically related to regional thrust structures. The uncertainty in the temperature estimates for the Cretaceous and Eocene vein formation prevents a more detailed assessment of the temporal evolution in meteoric water δ¹⁸O related to changing paleogeography. Meteoric water‐influenced δ¹⁸O values calculated here for Cretaceous to Eocene vein‐forming fluids are similar to those previously proposed for surface waters in the Eocene, and those observed for modern‐day precipitation, in this part of the Idaho‐Montana thrust belt.