ATTENUATION RELATIONSHIP FOR LOW MAGNITUDE EARTHQUAKES USING STANDARD SEISMOMETRIC RECORDS

Northwestern Italian weak-motion data were used to study attenuation characteristics of horizontal peak ground acceleration (PGA) and horizontal peak ground velocity (PGV) from earthquakes of local magnitudes (M _l ) up to 5.1. Data have been provided by the RSNI (Regional seismic network of Northwestern Italy) and RSLG (Regional seismic network of Lunigiana-Garfagnana) waveform database. The database consists of more than 14000 horizontal components recorded in the period 1999-2002 by both broadband and enlarged band seismometers. The accuracy of the procedure used to extract PGA values from the velocity recordings was verified comparing observed and derived PGA values at station STV2, which was equipped with both a temporary K2 Kinemctrtcs accelerometer and Guralp CMG40 broadband sensor. The attenuation of both peak ground acceleration and peak ground velocity was found to be logarithmically distributed with a strong attenuation for low distances (less than 50 km) and low M _l values (<3.0). The resulting equations are:

Log(PGA)=?3.19+0.87M?0.042M ²?1.92 Log(R)+0.249S,

Log(PGA)=?4.23+0.76M?0.018M²?1.56 Log(R)+0.230S,

where PGA is expressed in g, PGV is expressed in m/s, M is local magnitude, R is the hypocentral distance in kilometers and S is a dummy variable assuming values of 0 and 1 for rock and soil respectively. For increasing distance and magnitude, both PGA and PGV values show a linear distribution. The validity range of the obtained attenuation relationships is 0–200 km for distances and M _l up to 4.5. Sensitivity studies performed by analysis of residuals, showed that predicted PGA and PGV values are stable with respect to reasonable variations of the model and distances providing the data. Comparisons with attenuation relationships proposed for Italian region, derived from strong motion records, are also presented.     

Tectonic control over large-scale diffuse degassing in eastern Sicily (Italy)

Eastern Sicily (southern Italy) is characterised by the presence of many natural gas emissions (mofettes, mud volcanoes). These gases are mostly carbon dioxide and methane, with minor amounts of helium, hydrogen, carbon monoxide and hydrocarbons. In this study, the extent and orientation of soil gas anomalies (He and CO₂) were investigated on a wide area (approximately 110 km²) located just SW of Mt. Etna. From a structural point of view, this area lays on a typical foredeep–foreland system that marks the boundary between the southern part of the Eurasian plate and the northern part of the African plate in the central Mediterranean. No tectonic structure was revealed in this area by surface geological surveys. Very high soil emissions were found, and their spatial pattern reveals the existence of some active faults all directed about N50°E. This direction coincides with that of two major fault systems that cut eastern Sicily and are evident, respectively, NE and SW of the study area. Soil gas data suggest that these fault systems are the expression of a single continuous structural line which is probably responsible for the past and present magma uprise in eastern Sicily. Isotopic values of carbon of CO₂ suggest a minor contribution of organic carbon. Moreover, in the highest degassing sites the isotopic values of He found in association with CO₂ (He abundance = 11–70 p.p.m.; R/R_a between 6.0 and 6.2) suggest that both gases are mantle derived. The extent of the areas affected by high gas emissions and the amounts of deep CO₂ emitted in the investigated area (several hundred tonnes per day) may provide additional supporting evidence of a mantle upwelling taking place beneath this region.     

Hydrogeology of the Krafla geothermal system,northeast Iceland

The Krafla geothermal system is located in Iceland's northeastern neovolcanic zone, within the Krafla central volcanic complex. Geothermal fluids are superheated steam closest to the magma heat source, two‐phase at higher depths, and sub‐boiling at the shallowest depths. Hydrogen isotope ratios of geothermal fluids range from ?87‰, equivalent to local meteoric water, to ?94‰. These fluids are enriched in ¹⁸O relative to the global meteoric line by +0.5–3.2‰. Calculated vapor fractions of the fluids are 0.0–0.5 wt% (~0–16% by volume) in the northwestern portion of the geothermal system and increase towards the southeast, up to 5.4 wt% (~57% by volume). Hydrothermal epidote sampled from 900 to 2500 m depth has δD values from ?127 to ?108‰, and δ¹⁸O from ?13.0 to ?9.6‰. Fluids in equilibrium with epidote have isotope compositions similar to those calculated for the vapor phase of two‐phase aquifer fluids. We interpret the large range in δD_EPIDOTE and δ¹⁸O_EPIDOTE across the system and within individual wells (up to 7‰ and 3.3‰, respectively) to result from variable mixing of shallow sub‐boiling groundwater with condensates of vapor rising from a deeper two‐phase reservoir. The data suggest that meteoric waters derived from a single source in the northwest are separated into the shallow sub‐boiling reservoir, and deeper two‐phase reservoir. Interaction between these reservoirs occurs by channelized vertical flow of vapor along fractures, and input of magmatic volatiles further alters fluid chemistry in some wells. Isotopic compositions of hydrothermal epidote reflect local equilibrium with fluids formed by mixtures of shallow water, deep vapor condensates, and magmatic volatiles, whose ionic strength is subsequently derived from dissolution of basalt host rock. This study illustrates the benefits of combining phase segregation effects in two‐phase systems during analysis of wellhead fluid data with stable isotope values of hydrous alteration minerals when evaluating the complex hydrogeology of volcano‐hosted geothermal systems.     

Geothermometry and geobarometry of overpressured environments in Qiongdongnan Basin,South China Sea*

We demonstrate the use of PVT fluid inclusion modelling in the calculation of palaeofluid formation pressures, using samples from the YC21‐1‐1 and YC21‐1‐4 wells in the YC21‐1 structural closure, Qiongdongnan Basin, South China Sea. Homogenisation temperatures and gas/liquid ratios were measured in aqueous fluid inclusions, and associated light hydrocarbon/CO₂‐bearing inclusions, and their compositions were determined using a crushing technique. The vtflinc software was used to construct P–T phase diagrams that enabled derivation of the minimum trapping pressure for each order of fluid inclusion. Through the projection of average homogenisation temperatures (155, 185.5 and 204.5°C) for three orders of fluid inclusion on the thermal‐burial history diagram of the Oligocene Yacheng and Lingshui formations, their trapping times were constrained at 4.3, 2.1 and 1.8 Ma, respectively. The formation pressure coefficient, the ratio of fluid pressure/hydrostatic pressure established by PVT modelling coupled with DST data, demonstrates that one and a half cycles of pressure increase–discharge developed in the Yacheng and Lingshui formations for about 4.3 Ma. In comparison, the residual formation pressure determined by 2D numerical modelling in the centre of LeDong depression shows two and a half pressure increase–discharge cycles for about 28 Ma. The two different methods suggest that a high fluid potential in the Oligocene reservoir of the YC21‐1 structure developed at two critical stages for regional oil and natural gas migration and accumulation (5.8 and 2.0 Ma, respectively). Natural gas exploration in this area is therefore not advisable.     

Fluidized sandstone intrusions as an indicator of Paleostress orientation, Santa Cruz, California

The late Miocene sandstone intrusions of northern Santa Cruz County, California, are the largest subaerial exposures of clastic intrusions on earth. The intrusions are sourced from a sandstone, underlying mudstone, accumulated in an outer shelf to upper slope environment. Dikes are the most frequent intrusion type, reach the greatest thickness and tend to strike north‐east and dip steeply. One giant dike is more than 150 m wide. Sills are least frequent, locally > 8 m thick and have no clear preferred geographical distribution. Clustered intrusions are commonly < 10 cm thick and mostly composed of dikes of various attitudes. The majority of the intrusions probably were injected shallowly as some extrude onto the seafloor. The local seafloor extrusion also indicates injection during the deposition of the Santa Cruz Mudstone (7–9 Ma). The intrusions are concentrated at the basin margin. Fluid pressure at the centre of the basin and perhaps hydrocarbons were communicated to the basin margin through the then sand, causing fluid overpressures that contributed to the fluidization and intrusion into the overlying mudstone. Primarily north‐east‐striking, steeply dipping dikes and secondarily, shallowly dipping sills are most significant in terms of regional connectivity of intrusions and physical dilation of the formation. The orientation of the dikes and sills indicates a regional stress field with a horizontal NE–SW maximum and a NW–SE minimum compressive direction. The simultaneous development of dikes and sills suggests similar magnitudes of the minimum and intermediate principal stresses. Preferential weakness along bedding contributed to the development of sills. Palaeomagnetic data indicate no significant block rotation around a vertical axis. The maximum principal stress direction indicated by the intrusions is about 55° to the San Gregorio Fault and about 70° to the San Andreas Fault during the late Miocene. This stress field is similar to the modern stress field and suggests moderate fault weakness.     

An Investigation of Characteristic Periods of Seismic Ground Motions

The design seismic base shear was obtained from the spectral elastic acceleration S_a divided by a system behavior factor R, accounting for ductility and overstrength. The behavior factor is currently taken as a constant for a given type of structures in various codes regardless of structural periods. In fact, the behavior factor is also a spectrum varying with the natural periods of structures. In order to understand the relationship between the spectral values and the corresponding characteristic periods in these two spectra, S_a and R_μ, this article carries out an investigation into the characteristic periods of 370 seismic ground motions from 4 site types. It is found that the periods T_ga at which the peak values appear in the S_a spectra are much less than the periods T _gR at which the R_μ spectra take a maximum value. Two characteristic periods are necessary to determine the seismic action if a more elaborate procedure is required in practice. Statistical study on these two periods is carried out for the 370 records, and results are presented. For site types A–D, the ratio of these two periods has a statistically averaged value of 5.5–6.7.

The maximum input energy S _EI , relative velocity S _v , power density P _SD , and the Fourier amplitude F _S spectra were constructed to determine their characteristic periods, respectively. These four spectra predict similar characteristic periods to T _gR . T _gR is very close to the characteristic period T _gd of the elastic displacement spectra.

Analysis of SDOF systems under combined harmonic excitations shows that the S_a spectrum is more sensitive to high-frequency excitations, while the displacement spectrum is more sensitive to long period excitations. For the elastic-plastic S_a spectra, peak values tend to appear at shorter periods even the amplitudes of the longer periods are greater than that of the shorter period. This provides an explanation on different characteristic periods in the S_a and R_μ spectra.     

The role of the stress regime on microseismicity induced by overpressure and cooling in geologic carbon storage

Fluid injection in deep geological formations usually induces microseismicity. In particular, industrial‐scale injection of CO₂ may induce a large number of microseismic events. Since CO₂ is likely to reach the storage formation at a lower temperature than that corresponding to the geothermal gradient, both overpressure and cooling decrease the effective stresses and may induce microseismicity. Here, we investigate the effect of the stress regime on the effective stress evolution and fracture stability when injecting cold CO₂ through a horizontal well in a deep saline formation. Simulation results show that when only overpressure occurs, the vertical total stress remains practically constant, but the horizontal total stresses increase proportionally to overpressure. These hydro‐mechanical stress changes result in a slight improvement in fracture stability in normal faulting stress regimes because the decrease in deviatoric stress offsets the decrease in effective stresses produced by overpressure. However, fracture stability significantly decreases in reverse faulting stress regimes because the size of the Mohr circle increases in addition to being displaced towards failure conditions. Fracture stability also decreases in strike slip stress regimes because the Mohr circle maintains its size and is shifted towards the yield surface a magnitude equal to overpressure minus the increase in the horizontal total stresses. Additionally, cooling induces a thermal stress reduction in all directions, but larger in the out‐of‐plane direction. This stress anisotropy causes, apart from a displacement of the Mohr circle towards the yield surface, an increase in the size of the Mohr circle. These two effects decrease fracture stability, resulting in the strike slip being the least stable stress regime when cooling occurs, followed by the reverse faulting and the normal faulting stress regimes. Thus, characterizing the stress state is crucial to determine the maximum sustainable injection pressure and maximum temperature drop to safely inject CO₂.     

Metamorphic and basin fluids in quartz–carbonate–sulphide veins in the SW Scottish Highlands: a stable isotope and fluid inclusion study

Metalliferous (Fe–Cu–Pb–Zn) quartz–carbonate–sulphide veins cut greenschist to epidote–amphibolite facies metamorphic rocks of the Dalradian, SW Scottish Highlands, with NE–SW to NW–SE trends, approximately parallel or perpendicular to regional structures. Early quartz was followed by pyrite, chalcopyrite, sphalerite, galena, barite, late dolomite–ankerite and clays. Both quartz–sulphide and carbonate vein mineralisation is associated with brecciation, indicating rapid release of fluid overpressure and hydraulic fracturing. Two distinct mineralising fluids were identified from fluid inclusion and stable isotope studies. High temperature (>350°C) quartz‐precipitating fluids were moderately saline (4.0–12.7 wt.% NaCl equivalent) with low (approximately 0.05). Quartz δ¹⁸O (+11.7 to +16.5‰) and sulphide δ³⁴S (?13.6 to ?1.1‰) indicate isotopic equilibrium with host metasediments (rock buffering) and a local metasedimentary source of sulphur. Later, low‐temperature (T_H = 120–200°C) fluids, probably associated with secondary carbonate, barite and clay formation, were also moderately saline (3.8–9.1 wt.% NaCl equivalent), but were strongly enriched in ¹⁸O relative to host Dalradian lithologies, as indicated by secondary dolomite–ankerite (δ¹⁸O = +17.0 to +29.0‰, δ¹³C = ?1.0 to ?3.0‰). Compositions of carbonate–forming fluids were externally buffered. The veins record the fluid–rock interaction history of metamorphic host rocks during cooling, uplift and later extension. Early vein quartz precipitated under retrograde greenschist facies conditions from fluids probably derived by syn‐metamorphic dehydration of deeper, higher‐grade rocks during uplift and cooling of the Caledonian metamorphic complex. Veins are similar to those of mesothermal veins in younger Phanerozoic metamorphic belts, but are rare in the Scottish Dalradian. Early quartz veins were reactivated by deep penetration of low‐temperature basin fluids that precipitated carbonate and clays in veins and adjacent Dalradian metasediments throughout the SW Highlands, probably in the Permo‐Carboniferous. This event is consistent with paragenetically ambiguous barite with δ³⁴S characteristic of late Palaeozoic basinal brines.     

3D numerical modeling of hydrothermal processes during the lifetime of a deep geothermal reservoir

Understanding hydrothermal processes during production is critical to optimal geothermal reservoir management and sustainable utilization. This study addresses the hydrothermal (HT) processes in a geothermal research doublet consisting of the injection well E GrSk3/90 and production well Gt GrSk4/05 at the deep geothermal reservoir of Groß Schönebeck (north of Berlin, Germany) during geothermal power production. The reservoir is located between ?4050 to ?4250 m depth in the Lower Permian of the Northeast German Basin. Operational activities such as hydraulic stimulation, production (T = 150°C; Q = ?75 m³ h^?1; C = 265 g l^?1) and injection (T = 70°C; Q = 75 m³ h^?1; C = 265 g l^?1) change the HT conditions of the geothermal reservoir. The most significant changes affect temperature, mass concentration and pore pressure. These changes influence fluid density and viscosity as well as rock properties such as porosity, permeability, thermal conductivity and heat capacity. In addition, the geometry and hydraulic properties of hydraulically induced fractures vary during the lifetime of the reservoir. A three‐dimensional reservoir model was developed based on a structural geological model to simulate and understand the complex interaction of such processes. This model includes a full HT coupling of various petrophysical parameters. Specifically, temperature‐dependent thermal conductivity and heat capacity as well as the pressure‐, temperature‐ and mass concentration‐dependent fluid density and viscosity are considered. These parameters were determined by laboratory and field experiments. The effective pressure dependence of matrix permeability is less than 2.3% at our reservoir conditions and therefore can be neglected. The results of a three‐dimensional thermohaline finite‐element simulation of the life cycle performance of this geothermal well doublet indicate the beginning of thermal breakthrough after 3.6 years of utilization. This result is crucial for optimizing reservoir management. Geofluids (2010) 10 , 406–421     

The effect of deformation on permeability development in anhydrite and implications for caprock integrity during geological storage of CO2

Geological storage of CO₂ in depleted oil and gas reservoirs is one of the most promising options to reduce atmospheric CO₂ concentrations. Of great importance to CO₂ mitigation strategies is maintaining caprock integrity. Worldwide many current injection sites and potential storage sites are overlain by anhydrite‐bearing seal formations. However, little is known about the magnitude of the permeability change accompanying dilatation and failure of anhydrite under reservoir conditions. To this extent, we have performed triaxial compression experiments together with argon gas permeability measurements on Zechstein anhydrite, which caps many potential CO₂ storage sites in the Netherlands. Our experiments were performed at room temperature at confining pressures of 3.5–25 MPa. We observed a transition from brittle to semi‐brittle behaviour over the experimental range, and peak strength could be described by a Mogi‐type failure envelope. Dynamic permeability measurements showed a change from ‘impermeable’ (<10^?21 m²) to permeable (10^?16 to 10^?19 m²) as a result of mechanical damage. The onset of measurable permeability was associated with an increase in the rate of dilatation at low pressures (3.5–5 MPa), and with the turning point from compaction to dilatation in the volumetric versus axial strain curve at higher pressures (10–25 MPa). Sample permeability was largely controlled by the permeability of the shear faults developed. Static, postfailure permeability decreased with increasing effective mean stress. Our results demonstrated that caprock integrity will not be compromised by mechanical damage and permeability development. Geofluids (2010) 10 , 369–387     

Mineralogy and fluid inclusion gas chemistry of production well mineral scale deposits at the Dixie Valley geothermal field,USA

At the Dixie Valley geothermal field, Nevada, USA, fluid boiling triggered the precipitation of carbonate scale minerals in concentric bands around tubing inserted into production well 28–33. When the tubing was removed, this mineral scale was sampled at 44 depth intervals between the wellhead and 1227 m depth. These samples provide a unique opportunity to evaluate the effects of fluid boiling on the scale mineralogy and geochemistry of the vapor and liquid phase. In this study, the mineralogy of the scale deposits and the composition of the fluid inclusion gases trapped in the mineral scales were analyzed. The scale consists mainly of calcite from 670–1112 m depth and aragonite from 1125 to 1227 m depth, with traces of quartz and Mg‐smectite. Mineral textures, including hopper growth, twinning, and fibrous growth in the aragonite and banded deposits of fine grained calcite crystals, are the result of progressive boiling. The fluid inclusion noncondensable gas was dominated by CO₂. However, significant variations in He relative to N₂ and Ar provide evidence that the geothermal reservoir consists of mixed source deeply circulating reservoir water and shallow, air saturated meteoric water. Gas analyses for many inclusions also showed higher CH₄ and H₂ relative to CO₂ than measured in gas sampled from this well, other production wells, and fumaroles. These inclusions are interpreted to have trapped CH₄‐ and H₂‐enriched gas resulting from early stages of boiling.     

Geothermal fluids circulation at Caldas do Moledo area, Northern Portugal: geochemical and isotopic signatures

A hydrogeological conceptual model of the Caldas do Moledo geothermal site is proposed that shows mixing between geothermal waters and local shallow groundwaters. Stable isotope values of Caldas do Moledo geothermal waters indicate recharge areas located at relatively high altitudes (850–1250 m a.s.l.). The NW–SE Vigo–Régua shear zone plays an important role in fluid recharge and circulation towards the NNE–SSW Régua–Verin fault system, forming a path for ascent of geothermal fluids. The apparent ¹⁴C age of geothermal fluids (15.66 ± 2.86 ka BP) was estimated in the total dissolved inorganic carbon (TDIC). Geothermometer calculations indicate that, assuming a conductive temperature gradient of 32°C per kilometer for northern Portugal, the maximum depth of circulation is roughly 1.8 ± 0.4 km. The K, Ca and SO₄ concentrations found in some Caldas do Moledo geothermal spring waters show mixing between deep geothermal and shallow groundwater systems. Local shallow groundwaters showing the highest SO₄ concentrations were found at low elevation areas, originating from fertilisers and pesticides applied to the Port wine vineyards in the Douro River valley. Geothermal waters from boreholes AC1 and AC2 do not show evidences of direct pollution from the spreading of such agrochemicals.     

Temperature‐dependent Li isotope ratios in Appalachian Plateau and Gulf Coast Sedimentary Basin saline water

Lithium (Li) concentrations of produced water from unconventional (horizontally drilled and hydraulically fractured shale) and conventional gas wells in Devonian reservoirs in the Appalachian Plateau region of western Pennsylvania range from 0.6 to 17 mmol kg^?1, and Li isotope ratios, expressed as in δ⁷Li, range from +8.2 to +15‰. Li concentrations are as high as 40 mmol kg^?1 in produced waters from Plio‐Pleistocene through Jurassic‐aged reservoirs in the Gulf Coast Sedimentary Basin analyzed for this study, and δ⁷Li values range from about +4.2 to +16.6‰. Because of charge‐balance constraints and rock buffering, Li concentrations in saline waters from sedimentary basins throughout the world (including this study) are generally positively correlated with chloride (Cl), the dominant anion in these fluids. Li concentrations also vary with depth, although the extent of depth dependence differs among sedimentary basins. In general, Li concentrations are higher than expected from seawater or evaporation of seawater and therefore require water–mineral reactions that remove lithium from the minerals. Li isotope ratios in these produced waters vary inversely with temperature. However, calculations of temperature‐dependent fractionation of δ⁷Li between average shale δ⁷Li (?0.7‰) and water result in δ⁷Li_water that is more positive than that of most produced waters. This suggests that aqueous δ⁷Li may reflect transport of water from depth and/or reaction with rocks having δ⁷Li lighter than average shale.     

Permeability of the Mercia Mudstone: suitability as caprock to carbon capture and storage sites

The Upper Triassic Mercia Mudstone is the caprock to potential carbon capture and storage (CCS) sites in porous and permeable Lower Triassic Sherwood Sandstone reservoirs and aquifers in the UK (primarily offshore). This study presents direct measurements of vertical (k_v) and horizontal (k_h) permeability of core samples from the Mercia Mudstone across a range of effective stress conditions to test their caprock quality and to assess how they will respond to changing effective stress conditions that may occur during CO₂ injection and storage. The Mercia samples analysed were either clay‐rich (muddy) siltstones or relatively clean siltstones cemented by carbonate and gypsum. Porosity is fairly uniform (between 7.4 and 10.7%). Porosity is low either due to abundant depositional illite or abundant diagenetic carbonate and gypsum cements. Permeability values are as low as 10^?20 m² (10nD), and therefore, the Mercia has high sealing capacity. These rocks have similar horizontal and vertical permeabilities with the highest k_h/k_v ratio of 2.03 but an upscaled k_h/k_v ratio is 39, using the arithmetic mean of k_h and the harmonic mean of k_v. Permeability is inversely related to the illite clay content; the most clay‐rich (illite‐rich) samples represent very good caprock quality; the cleaner Mercia Mudstone samples, with pore‐filling carbonate and gypsum cements, represent fair to good caprock quality. Pressure sensitivity of permeability increases with increasing clay mineral content. As pore pressure increases during CO₂ injection, the permeability of the most clay‐rich rocks will increase more than carbonate‐ and gypsum‐rich rocks, thus decreasing permeability heterogeneity. The best quality Mercia Mudstone caprock is probably not geochemically sensitive to CO₂ injection as illite, the cause of the lowest permeability, is relatively stable in the presence of CO₂–water mixtures.     

Interpretation of observed fluid potential patterns in a deep sedimentary basin under tectonic compression: Hungarian Great Plain, Pannonian Basin

The ≈ 40 000 km² Hungarian Great Plain portion of the Pannonian Basin consists of a basin fill of 100 m to more than 7000 m thick semi‐ to unconsolidated marine, deltaic, lacustrine and fluviatile clastic sediments of Neogene age, resting on a strongly tectonized Pre‐Neogene basement of horst‐and‐graben topography of a relief in excess of 5000 m. The basement is built of a great variety of brittle rocks, including flysch, carbonates and metamorphics. The relatively continuous Endr?d Aquitard, with a permeability of less than 1 md (10^?15 m²) and a depth varying between 500 and 5000 m, divides the basin's rock framework into upper and lower sequences of highly permeable rock units, whose permeabilities range from a few tens to several thousands of millidarcy. Subsurface fluid potential and flow fields were inferred from 16 192 water level and pore pressure measurements using three methods of representation: pressure–elevation profiles; hydraulic head maps; and hydraulic cross‐sections. Pressure–elevation profiles were constructed for eight areas. Typically, they start from the surface with a straight‐line segment of a hydrostatic gradient (γ_st = 9.8067 MPa km^?1) and extend to depths of 1400–2500 m. At high surface elevations, the gradient is slightly smaller than hydrostatic, while at low elevations it is slightly greater. At greater depths, both the pressures and their vertical gradients are uniformly superhydrostatic. The transition to the overpressured depths may be gradual, with a gradient of γ_dyn = 10–15 MPa km^?1 over a vertical distance of 400–1000 m, or abrupt, with a pressure jump of up to 10 MPa km^?1 over less than 100 m and a gradient of γ_dyn > 20 MPa km^?1. According to the hydraulic head maps for 13 100–500 m thick horizontal slices of the rock framework, the fluid potential in the near‐surface domains declines with depth beneath positive topographic features, but it increases beneath depressions. The approximate boundary between these hydraulically contrasting regions is the 100 m elevation contour line in the Duna–Tisza interfluve, and the 100–110 m contours in the Nyírség uplands. Below depths of ≈ 600 m, islets of superhydrostatic heads develop which grow in number, areal extent and height as the depth increases; hydraulic heads may exceed 3000 m locally. A hydraulic head ‘escarpment’ appears gradually in the elevation range of ? 1000 to ? 2800 m along an arcuate line which tracks a major regional fault zone striking NE–SW: heads drop stepwise by several hundred metres, at places 2000 m, from its north and west sides to the south and east. The escarpment forms a ‘fluid potential bank’ between a ‘fluid potential highland’ (500–2500 m) to the north and west, and a ‘fluid potential basin’ (100–500 m) to the south and east. A ‘potential island’ rises 1000 m high above this basin further south. According to four vertical hydraulic sections, groundwater flow is controlled by the topography in the upper 200–1700 m of the basin; the driving force is orientated downwards beneath the highlands and upwards beneath the lowlands. However, it is directed uniformly upwards at greater depths. The transition between the two regimes may be gradual or abrupt, as indicated by wide or dense spacing of the hydraulic head contours, respectively. Pressure ‘plumes’ or ‘ridges’ may protrude to shallow depths along faults originating in the basement. The basement horsts appear to be overpressured relative to the intervening grabens. The principal thesis of this paper is that the two main driving forces of fluid flow in the basin are gravitation, due to elevation differences of the topographic relief, and tectonic compression. The flow field is unconfined in the gravitational regime, whereas it is confined in the compressional regime. The nature and geometry of the fluid potential field between the two regimes are controlled by the sedimentary and structural features of the rock units in that domain, characterized by highly permeable and localized sedimentary windows, conductive faults and fracture zones. The transition between the two potential fields can be gradual or abrupt in the vertical, and island‐like or ridge‐like in plan view. The depth of the boundary zone can vary between 400 and 2000 m. Recharge to the gravitational regime is inferred to occur from infiltrating precipitation water, whereas that to the confined regime is from pore volume reduction due to the basement's tectonic compression.     

Fluids,basin analysis,and mineral deposits

Sedimentary basins are the largest structures on the surface of our planet and the most significant sources of energy‐related commodities. With time, sedimentary successions in basins normally are subjected to increasingly intense diagenesis that results in differential evolution of basin hydrology. This hydrologic structure is in turn vitally important in determining how and where deposition of metals may occur. Fluids in all basins originate and flow as a result of sedimentological and tectonic events, so that fluid histories should reflect the control of both lithology and tectonism on ore deposition. Sandstone lithologies, in particular, reflect fluid‐flow events because they are normally the major aquifers in basins. However, early cementation results in occlusion of primary permeability in some facies (diagenetic aquitards) whereas in others, permeability develops due to the dissolution of unstable grains (diagenetic aquifers). Particularly for ore deposits in Precambrian basins, identification of paleohydrologic systems during basin evolution requires the integration of data derived from tectonics, sedimentology, stratigraphy, diagenesis, geochemistry and geology. Assessment of all these data is a prerequisite for the ‘holistic basin analysis’ needed to guide the search for basin‐hosted ores. Recent results from the Paleoproterozoic Mt Isa and McArthur basins in northern Australia serve as a template for exploring for mineral deposits in basins. Basinal fluids were saline, 200–300°C and evolved primarily from meteoric water in the Mt Isa Basin and from seawater in the McArthur Basin during burial to depths of 4–12 km. The δD_fluid and δ¹⁸O_fluid values in these brines were isotopically identical to those in the Zn‐Pb, Cu and U deposits. Geochemical changes of various lithologies during alteration support detrital minerals as the major source of the U, and volcanic units proximal to diagenetic aquifers as a source for the transition metals. Ages of diagenetic phases extracted from aquifer lithologies reveal that fluid migration from the diagenetic aquifers effectively covers the period of formation for U, Zn‐Pb and Cu mineralization, and that the deposits formed in response to tectonic events reflected in the apparent polar wandering path for the area. Sequence stratigraphic analysis and models of fluid flow also indicate that basinal reservoirs were likely sources for the mineralizing fluids. Thus, diagenetic aquifer lithologies were being drained of fluids at the same time as the deposits were forming from fluids that were chemically and isotopically similar, linking diagenesis and fluid events within the basin to the formation of the deposits.     

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.     

Fluid inclusion and isotopic constraints on the origin of ore‐forming fluid of the Jinman–Liancheng vein Cu deposit in the Lanping Basin,western Yunnan,China

Extensive quartz–carbonate–Cu sulfide veins occur in clastic rocks and are spatially related to Paleocene granites in the western border of the Lanping Basin, western Yunnan, China. Abundant aqueous‐carbonic fluid inclusions occur in these veins but their origin is debated. In the Jinman–Liancheng deposit, individual primary inclusion groups contain either exclusively liquid‐rich inclusions (Gl), or coexisting liquid‐rich and vapor‐rich inclusions (Glv). Microthermometry and estimate of CO₂ content indicate that type Gl inclusions either have homogenization temperatures (Th) 238–263°C and contain c. 3.9–5.5 mole % CO₂, or have Th 178–222°C and contain c. 1.6–3.2 mole % CO₂. Type Glv inclusions are thought to represent samples of fluid unmixing that occurred at 183–218°C. At that time, the liquid phase in the unmixing fluid may contain c. 2.0–3.3 mole % CO₂. As such, the correlation of CO₂ content with Th for type Gl inclusions is thought to be caused by fluid unmixing with decreasing temperature and subsequent CO₂ escape. δ¹⁸O and δD values of the parent water mainly fall in the field below that of primary magmatic water, indicative of fluid derivation from degassed (in open system) magmatic water, with no contributions from basinal or meteoric water. Initial Sr isotopic compositions of hydrothermal carbonates suggest that the fluid was magmatic, probably derived from the Paleogene granites. δ¹³C_PDB values (?4‰ to ?7‰) of the hydrothermal carbonates and δ³⁴S_V‐CDT values of sulfides (mainly ?11‰ to +5‰) indicate that the carbon and sulfur can be derived from (degassed) magma and/or nonmagmatic sources. The CO₂‐rich and magmatic‐water‐derived fluid at Jinman–Liancheng differs from the CO₂‐poor and basinally derived fluid in sediment‐hosted stratiform Cu (SSC) deposits, which suggests that there are no genetic linkages between the vein Cu and SSC deposits in the Lanping Basin.     

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.     

Hydrothermal palaeofluid circulation in the fracture network of the Baksa Gneiss Complex of SW Pannonian Basin,Hungary

A well‐developed fracture‐filling network is filled by dominantly Ca‐Al‐silicate minerals that can be found in the polymetamorphic rock body of the Baksa Gneiss Complex (SW Hungary). Detailed investigation of this vein network revealed a characteristic diopside→epidote→sphalerite→albite ± kfeldspar→chlorite1 ± prehnite ± adularia→chlorite2→chlorite3→pyrite→calcite1→calcite2→calcite3 fracture‐filling mineral succession. Thermobarometric calculations (two feldspar: 230–336°C; chlorites: approximately 130–300°C) indicate low‐temperature vein formation conditions. The relative succession of chlorites in the mineral sequence combined with the calculated formation temperatures reveals a cooling trend during precipitation of the different chlorite phases (T_chlorite1: 260 ± 32°C →T_chlorite2: 222 ± 20°C →T_chlorite3: 154 ± 13°C). This cooling trend can be supported by the microthermometry data of primary fluid inclusions in diopside (T_h: 276–362°C) and epidote (T_h: 181–359°C) phases. The identical chemical character (0.2–1.5 eq. wt% NaCl) of these inclusions mean that vein mineralization occurred in a same fluid environment. The high trace element content (e.g. As, Cu, Zn, Mn) and Co/Ni ratio approximately 1–5 of pyrite grains support the postmagmatic hydrothermal origin of the veins. The vein microstructure and identical fluid composition indicate that vein mineralization occurred in an interconnected fracture system where crystals grew in fluid filled cracks. Vein system formed at approximately <200 MPa pressure conditions during cooling from approximately 480°C to around 150°C. The rather different fluid characteristics (T_h: 75–124°C; 17.5–22.6 eq. wt% CaCl₂) of primary inclusions of calcite1 combining with the special δ¹⁸O signature of fluid from which this mineral phase precipitated refer to hydrological connection between the crystalline basement and the sedimentary cover.