Salinity structure of the central North Slope foreland basin, Alaska, USA: implications for pathways of past and present topographically driven regional fluid flow

Previous studies of the areal variation in heat flow in the National Petroleum Reserve Alaska (NPRA) support the existence of an active topographically driven regional fluid flow regime in this central part of the North Slope foreland basin. Drilling records and wireline logs for over 30 wells drilled in the NPRA provide additional field information, which can be used to further constrain interpretation of the pattern of regional flow of basinal waters within the NPRA. Hydraulic heads estimated from drilling mud weights show that ground water flow occurs generally from south to north, but with divergence to the north‐east and north‐west away from the central part of the NPRA towards coastal areas of elevated shallow heat flow. Salinities calculated from SP logs range from less than 1 g L^?1, to marine values of 35 g L^?1, to hypersaline values of over 150 g L^?1. The entire upper sedimentary section to a depth of 2 km or more in the eastern part of the NPRA has been preferentially flushed with meteoric water through an area corresponding to the sandiest portion of the Nanushuk group. Deeper areas of low salinity occur within the Sadlerochit and Lisburne sections. The pattern of regional flow in the east is complicated, however, by the presence of a large mass of hypersaline water at depth. It is not known whether these brines are being displaced laterally and upward towards the discharge end of the basin or whether fresher waters are simply riding up over the top. Deep, hypersaline waters also occur in fault slices in the Brooks Range and have survived meteoric flushing. The brines were probably formed at the time of deposition of the Lisburne carbonates. The fluid flow regime to the west is different. Low‐salinity waters may be flowing northward underneath this section through the Ellesmerian section and discharging upward nearer the coast. However, sparse well log control severely limits what can be deduced about the details of flow paths in the central and western parts of the NPRA.     

Impacts of Pleistocene glaciation on large‐scale groundwater flow and salinity in the Michigan Basin

Pleistocene melting of kilometer‐thick continental ice sheets significantly impacted regional‐scale groundwater flow in the low‐lying stable interiors of the North American and Eurasian cratons. Glacial meltwaters penetrated hundreds of meters into the underlying sedimentary basins and fractured crystalline bedrock, disrupting relatively stagnant saline fluids and creating a strong disequilibrium pattern in fluid salinity. To constrain the impact of continental glaciation on variable density fluid flow, heat and solute transport in the Michigan Basin, we constructed a transient two‐dimensional finite‐element model of the northern half of the basin and imposed modern versus Pleistocene recharge conditions. The sag‐type basin contains up to approximately 5 km of Paleozoic strata (carbonates, siliciclastics, and bedded evaporites) overlain by a thick veneer (up to 300 m) of glacial deposits. Formation water salinity increases exponentially from <0.5 g l^?1 total dissolved solids (TDS) near the surface to >350 g l^?1 TDS at over 800 m depth. Model simulations show that modern groundwater flow is primarily restricted to shallow glacial drift aquifers with discharge to the Great Lakes. During the Pleistocene, however, high hydraulic heads from melting of the Laurentide Ice Sheet reversed regional flow patterns and focused recharge into Paleozoic carbonate and siliciclastic aquifers. Dilute waters (<20 g l^?1 TDS) migrated approximately 110 km laterally into the Devonian carbonate aquifers, significantly depressing the freshwater‐saline water mixing zones. These results are consistent with ¹⁴C ages and oxygen isotope values of confined groundwaters in Devonian carbonates along the basin margin, which reflect past recharge beneath the Laurentide Ice Sheet (14–50 ka). Constraining the paleohydrology of glaciated sedimentary basins, such as the Michigan Basin, is important for determining the source and residence times of groundwater resources, in addition to resolving geologic forces responsible for basinal‐scale fluid and solute migration.     

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.     

Large‐scale chemical stratification of fluids in the crust: hydraulic and chemical data from the geothermal research site Urach,Germany

The Urach 3 research borehole in SW Germany has been drilled through a sedimentary cover sequence and reached gneisses of the Variscan crystalline basement at 1604 m below surface. An additional 2840 m has been drilled through fractured basement rocks. The borehole has been used for hydraulic tests in the context of a ‘hot dry rock’ (HDR) project. The sedimentary cover ranges from the Carboniferous to the Middle Jurassic (Dogger) in age and comprises mostly clastic sediments in the Paleozoic and limestone and shale in the Mesozoic. Water composition data from 10 different depths include samples from all major lithological units. The total dissolved solids (TDS) increases from the surface to about 650 m where it reaches 4.1 g l^?1 in Triassic limestone. In lower Triassic sandstones, TDS increases very sharply to 28.5 g l^?1 and the water is saturated with pure CO₂ gas. With increasing depth, TDS does not change much in the clastic sediments of the Permian and Carboniferous. The crystalline basement is marked by a very sharp increase in TDS to 55.5 g l^?1 at about 1770 m depth. TDS increases within the basement to more than 78.5 g l^?1 at about 3500 m depth. The data suggest that there is limited vertical chemical communication over long periods of time. The CO₂ gas cap in the lower Triassic sandstones requires a gastight cover. The chemical stratification of the fluids relates to the permeability structure of the crust at the Urach site and fits well with hydraulic and thermal data from the site.     

Upper crustal fluid flow in the outboard region of the Southern Alps,New Zealand

The currently active fluid regime within the outboard region of the Southern Alps, New Zealand was investigated using a combination of field observations, carbon‐ and oxygen‐stable isotopes from fault‐hosted calcites and interpretation of magnetotelluric (MT) data. Active faulting in the region is dominated by NE striking and N striking, oppositely dipping thrust fault pairs. Stable isotopic analyses of calcites hosted within these fault zones range from 10 to 25‰δ¹⁸O and from ?33 to 0‰δ¹³C. These values reflect mixing of three parent fluids: meteoric water, carbon‐exchanged groundwater and minor deeper rock‐exchanged fluids, at temperatures of 10–90°C in the upper 3–4 km of the crust. A broad, ‘U‐shaped’ zone of high electrical conductivity (maximum depth c. 28 km) underlies the central Southern Alps. In the ductile region of the crust, the high‐conductivity zone is subhorizontal. Near‐vertical zones of high conductivity extend upward to the surface on both sides of the conductive zone. On the outboard side of the orogen, the conductive zone reaches the surface coincident with the trace of the active Forest Creek Faults. Near‐surface flow is shown to dominate the outboard region. Topographically driven meteoric water interacts, on a kilometre scale, with either carbon‐exchanged groundwater or directly with organic material within Pliocene gravels, resulting in a distinctive low ¹³C signal within fault‐hosted calcites of the outboard region. The high‐strain zone in the lower crust focuses the migration of deeply sourced fluids upward to the base of the brittle–ductile transition. Interconnected fluid is imaged as a narrow vertical zone of high conductivity in the upper crust, implying continuous permeability and possibly buoyancy‐driven flow of deeply sourced fluids to higher levels of the crust where they are detected by the isotopic analysis of the fault‐hosted calcites.     

Abnormally pressured beds as barriers to diffusive solute transport in sedimentary basins

Diffusion can drive significant solute transport over millions of years, but ancient brines and large salinity gradients are still observed in deep sedimentary basins. Fluid flow within abnormally pressured beds may prevent diffusive transfer over geologically significant periods, if the abnormally pressured bed is surrounded by normally pressured beds. Analytic solutions based on sediment loading and unloading demonstrate that this effect should be considered in beds with a compressibility exceeding 10^?8 Pa^?1, with a thickness of 100 m or more, or a sedimentation rate exceeding 10^?5 m year^?1. Conditions favourable for our model of abnormally pressured beds appear common in sedimentary basins. Large salinity gradients associated with clay beds have previously been attributed to membrane effects, but flow patterns associated with abnormally pressured beds appear more robust in the presence of heterogeneity and discontinuities than membrane effects. Calculations suggest that thick underpressured shales in the Alberta basin may have allowed ancient evaporatively concentrated brines to be preserved beneath a vigorous topography‐driven flow system over the last 60 My. In the Illinois basin, drained overpressured beds may have limited solute transport across the New Albany shale until approximately 250 Ma. It is unlikely, however, that overpressures could have persisted long enough to explain concentration gradients observed in the modern basin. These gradients may instead reflect relatively recent halite dissolution above the New Albany shale.     

The upper continental crust, an aquifer and its fluid: hydaulic and chemical data from 4 km depth in fractured crystalline basement rocks at the KTB test site

Detailed information on the hydrogeologic and hydraulic properties of the deeper parts of the upper continental crust is scarce. The pilot hole of the deep research drillhole (KTB) in crystalline basement of central Germany provided access to the crust for an exceptional pumping experiment of 1‐year duration. The hydraulic properties of fractured crystalline rocks at 4 km depth were derived from the well test and a total of 23100 m³ of saline fluid was pumped from the crustal reservoir. The experiment shows that the water‐saturated fracture pore space of the brittle upper crust is highly connected, hence, the continental upper crust is an aquifer. The pressure–time data from the well tests showed three distinct flow periods: the first period relates to wellbore storage and skin effects, the second flow period shows the typical characteristics of the homogeneous isotropic basement rock aquifer and the third flow period relates to the influence of a distant hydraulic border, probably an effect of the Franconian lineament, a steep dipping major thrust fault known from surface geology. The data analysis provided a transmissivity of the pumped aquifer T = 6.1 × 10^?6 m² sec^?1, the corresponding hydraulic conductivity (permeability) is K = 4.07 × 10^?8 m sec^?1 and the computed storage coefficient (storativity) of the aquifer of about S = 5 × 10^?6. This unexpected high permeability of the continental upper crust is well within the conditions of possible advective flow. The average flow porosity of the fractured basement aquifer is 0.6–0.7% and this range can be taken as a representative and characteristic values for the continental upper crust in general. The chemical composition of the pumped fluid was nearly constant during the 1‐year test. The total of dissolved solids amounts to 62 g l^?1 and comprise mainly a mixture of CaCl₂ and NaCl; all other dissolved components amount to about 2 g l^?1. The cation proportions of the fluid (X_Ca approximately 0.6) reflects the mineralogical composition of the reservoir rock and the high salinity results from desiccation (H₂O‐loss) due to the formation of abundant hydrate minerals during water–rock interaction. The constant fluid composition suggests that the fluid has been pumped from a rather homogeneous reservoir lithology dominated by metagabbros and amphibolites containing abundant Ca‐rich plagioclase.     

Spatial variations in the salinity of pore waters in northern deep water Gulf of Mexico sediments: implications for pathways and mechanisms of solute transport

Spatial variations in the salinity of pore waters in sedimentary basins can provide important insight into basin-scale hydrogeologic processes. Although there have been numerous studies of brine seeps in the deep water Gulf of Mexico, much less is known about porewater salinities in the vast areas between seeps. A study has been made of spatial variation in pore water salinities in sediments in an approximately 500 km by 200 km area of the northern deep water (water depth >500 m) Gulf of Mexico sedimentary basin (GOM) to provide insight into pathways and mechanisms of solute transport in this portion of the basin. A second objective was to document salinities in the upper 500 m of the sedimentary section, the approximate depth to which methane hydrates, a potential future energy resource, may be stable. Elevated salinities would reduce the P – T stability range of hydrates. Salinities were calculated from borehole logs using a dual electrical conductivity model. Even though much of the northern GOM is underlain by allochthonous salt most of the undisturbed shallow sedimentary section has not been permeated by hypersaline waters. These waters are limited to areas near brine seeps. Hypersaline waters having salinities in excess of 100 g l⁻¹ become more common at subseafloor depths of 2 km and greater. A field study at Green Canyon 65 and published numerical simulations of fluid flow above tabular salt bodies suggest that brines produced by salt dissolution migrate laterally and pond above salt and/or within minibasins and that the dominant mechanism of vertical solute transport is a combination of compaction-driven advection and diffusion, not large-scale thermohaline overturn. Superimposed on this diffuse upward flux of dissolved salt is the more focused and localized expulsion of saline fluids up faults.     

Thermal convection in faulted extensional sedimentary basins: theoretical results from finite-element modeling

Thermal convection has the potential to be a significant and widespread mechanism of fluid flow, mass transport, and heat transport in rift and other extensional basins. Based on numerical simulation results, large‐scale convection can occur on the scale of the basin thickness, depending on the Rayleigh number for the basin. Our analysis indicates that for syn‐rift and early post‐rift settings with a basin thickness of 5 km, thermal convection can occur for basal heat flows ranging from 80 to 150 mW m^?2, when the vertical hydraulic conductivity is on the order of 1.5 m year^?1 and lower. The convection cells have characteristic wavelengths and flow patterns depending on the thermal and hydraulic boundary conditions. Steeply dipping extensional faults can provide pathways for vertical fluid flow across large thicknesses of basin sediments and can modify the dynamics of thermal convection. The presence of faults perturbs the thermal convective flow pattern and can constrain the size and locations of convection cells. Depending on the spacing of the faults and the hydraulic properties of the faults and basin sediments, the convection cells can be spatially organized to align with adjacent faults. A fault‐bounded cell occurs when one convection cell is constrained to occupy a fault block so that the up‐flow zone converges into one fault zone and the down‐flow zone is centred on the adjacent fault. A fault‐bounded cell pair occurs when two convection cells occupy a fault block with the up‐flow zone located between the faults and the down‐flow zones centred on the adjacent faults or with the reverse pattern of flow. Fault‐bounded cells and cell pairs can be referred to collectively as fault‐bounded convective flow. The flow paths in fault‐bounded convective flow can be lengthened significantly with respect to those of convection cells unperturbed by the presence of faults. The cell pattern and sense of circulation depend on the fault spacing, sediment and fault permeabilities, lithologic heterogeneity, and the basal heat flow. The presence of fault zones also extends the range of conditions for which thermal convection can occur to basin settings with Rayleigh numbers below the critical value for large‐scale convection to occur in a basin without faults. The widespread potential for the occurrence of thermal convection suggests that it may play a role in controlling geological processes in rift basins including the acquisition and deposition of metals by basin fluids, the distribution of diagenetic processes, the temperature field and heat flow, petroleum generation and migration, and the geochemical evolution of basin fluids. Fault‐bounded cells and cell pairs can focus mass and heat transport from longer flow paths into fault zones, and their discharge zones are a particularly favourable setting for the formation of sediment‐hosted ore deposits near the sea floor.     

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.     

Effects of episodic fluid flow on hydrocarbon migration in the Newport‐Inglewood Fault Zone,Southern California

Fault permeability may vary through time due to tectonic deformations, transients in pore pressure and effective stress, and mineralization associated with water‐rock reactions. Time‐varying permeability will affect subsurface fluid migration rates and patterns of petroleum accumulation in densely faulted sedimentary basins such as those associated with the borderland basins of Southern California. This study explores the petroleum fluid dynamics of this migration. As a multiphase flow and petroleum migration case study on the role of faults, computational models for both episodic and continuous hydrocarbon migration are constructed to investigate large‐scale fluid flow and petroleum accumulation along a northern section of the Newport‐Inglewood fault zone in the Los Angeles basin, Southern California. The numerical code solves the governing equations for oil, water, and heat transport in heterogeneous and anisotropic geologic cross sections but neglects flow in the third dimension for practical applications. Our numerical results suggest that fault permeability and fluid pressure fluctuations are crucial factors for distributing hydrocarbon accumulations associated with fault zones, and they also play important roles in controlling the geologic timing for reservoir filling. Episodic flow appears to enhance hydrocarbon accumulation more strongly by enabling stepwise build‐up in oil saturation in adjacent sedimentary formations due to temporally high pore pressure and high permeability caused by periodic fault rupture. Under assumptions that fault permeability fluctuate within the range of 1–1000 millidarcys (10^?15–10^?12 m²) and fault pressures fluctuate within 10–80% of overpressure ratio, the estimated oil volume in the Inglewood oil field (approximately 450 million barrels oil equivalent) can be accumulated in about 24 000 years, assuming a seismically induced fluid flow event occurs every 2000 years. This episodic petroleum migration model could be more geologically important than a continuous‐flow model, when considering the observed patterns of hydrocarbons and seismically active tectonic setting of the Los Angeles basin.     

Chlorine stable isotope ratios as tracer for pore‐water advection rates in a submarine gas‐hydrate field: implication for hydrate concentration

At Ocean Drilling Program Site 997 in Blake Ridge gas‐hydrate field in West Atlantic, pore‐water studies revealed a pronounced downward depletion of the heavy chlorine isotope to nearly ?4‰δ³⁷Cl at approximately 750 m below sea floor (mbsf) associated with a 10% downward chlorinity decrease. This is one of the stronger ³⁷Cl depletions hitherto reported for marine pore waters. Chlorinity reductions in hydrate‐bearing sediments commonly result from fresh‐water release by hydrate melting. However, in situ measurements at Site 997 suggest that >50% of the chlorinity reduction occurred prior to hydrate dissociation. Modeling the chlorinity profile shows that advection of a strongly ³⁷Cl‐depleted, low‐chlorinity water (506 mm ) from below the drilled sequence can explain the reduction prior to sampling. Fitting the model to the δ³⁷Cl curve yielded an advection rate of 0.18 mm year^?1. Diffusive mixing with near‐0‰‐δ³⁷Cl paleo‐seawater with maximum chlorinity at shallow subsurface depths (561 mm at approximately 20 mbsf) produced the smooth, steady trend. Separating that part of the freshening caused by advection and diffusion from that due to hydrate dissociation allowed estimation of average hydrate concentrations of 3.8% of the pore space (up to 24.5% in hydrate‐rich layers; near‐100% in rare massive hydrate layers). The deep‐seated source of the ³⁷Cl‐depleted, low‐chlorinity water remains unknown and might be located below the sedimentary section in the oceanic basement.     

Origin of palaeo‐waters in the Ordovician carbonates in Tahe oilfield,Tarim Basin: constraints from fluid inclusions and Sr,C and O isotopes

Petrographic features, isotopes, and trace elements were determined, and fluid inclusions were analyzed on fracture‐filling, karst‐filling and interparticle calcite cement from the Ordovician carbonates in Tahe oilfield, Tarim basin, NW China. The aim was to assess the origin and evolution of palaeo‐waters in the carbonates. The initial water was seawater diluted by meteoric water, as indicated by bright cathodoluminescence (CL) in low‐temperature calcite. The palaeoseawater was further buried to temperatures from 57 to 110°C, nonluminescent calcite precipitated during the Silurian to middle Devonian. Infiltration of meteoric water of late Devonian age into the carbonate rocks was recorded in the first generation of fracture‐ and karst‐filling dull red CL calcite with temperatures from <50°C to 83°C, low salinities (<9.0 wt%), high Mn contents and high ⁸⁶Sr/⁸⁷Sr ratios from 0.7090 to 0.7099. During the early Permian, ⁸⁷Sr‐rich hydrothermal water may have entered the carbonate rocks, from which precipitated a second generation of fracture‐filling and interparticle calcite and barite cements with salinities greater than 22.4 wt%, and temperatures from 120°C to 180°C. The hydrothermal water may have collected isotopically light CO₂ (possibly of TSR‐origin) during upward migration, resulting in hydrothermal calcite and the present‐day oilfield water having δ¹³C values from ?4.3 to ?13.8‰ and showing negative relationships of ⁸⁷Sr/⁸⁶Sr ratios to δ¹³C and δ¹⁸O values. However, higher temperatures (up to 187°C) and much lower salinities (down to 0.5 wt%) measured from some karst‐filling, giant, nonluminescent calcite crystals may suggest that hydrothermal water was deeply recycled, reduced (Fe‐bearing) meteoric water heated in deeper strata, or water generated from TSR during hydrothermal water activity. Mixing of hydrothermal and local basinal water (or diagenetically altered connate water) with meteoric waters of late Permian age and/or later may have resulted in large variations in salinity of the present oilfield waters with the lowest salinity formation waters in the palaeohighs.     

Analytical and numerical models of hydrothermal fluid flow at fault intersections

Fault intersections are the locus of hot spring activity and Carlin‐type gold mineralization within the Basin and Range, USA. Analytical and numerical solutions to Stokes equation suggest that peak fluid velocities at fault intersections increase between 20% and 47% when fracture apertures have identical widths but increase by only about 1% and 8% when aperture widths vary by a factor of 2. This suggests that fault zone intersections must have enlarged apertures. Three‐dimensional finite element models that consider intersecting 10‐ to 20‐m wide fault planes resulted in hot spring activity being preferentially located at fault zone intersections when fault zones were assigned identical permeabilities. We found that the onset of convection at the intersections of the fault zones occurred in our hydrothermal model over a narrow permeability range between 5 × 10^?13 and 7 × 10^?13 m². Relatively high vertical fluid velocities (0.3–3 m year^?1) extended away from the fault intersections for about 0.5–1.5 km. For the boundary conditions and fault plane dimensions used, peak discharge temperatures of 112°C at the water table occurred with an intermediate fault zone permeability of 5 × 10^?13 m². When fault plane permeability differed by a factor of 2 or more, the locus of hot spring activity shifted away from the intersections. However, increasing the permeability at the core of the fault plane intersection by 40% shifted the discharge back to the intersections. When aquifer units were assigned a permeability value equal to those of the fault planes, convective rolls developed that extend about 3 km laterally along the fault plane and into the adjacent aquifer.     

Origin of deep saline groundwaters in the Vienne granitic rocks (France): constraints inferred from boron and strontium isotopes

As part of a preliminary geological characterization programme to assess the feasibility of an underground laboratory in granitic rock, a series of 17 deep boreholes (maximum depth, 900 m) was drilled by ANDRA in the Vienne district, France. A salinity gradient was demonstrated in the granitic waters with concentrations varying from approximately 1 g L^?1 at 150 m depth at the top of the basement (beneath the sedimentary cover) to 10 g L^?1 in the deeper part (from 400 to 600 m depth). Sr and B isotope ratios were measured in order to better understand the origin of the salinity and to evaluate the degree of water–rock interaction in the system. The results obtained were compared to those of mineral spring waters emerging from the granitic basement in the Massif Central. Evidence in support of a significant marine contribution include: (i) the Cl–Br investigations agree with a marine origin for the saline groundwaters without evolution from seawater; (ii) the ⁸⁷Sr/⁸⁶Sr ratio of the Vienne deep groundwaters (0.7078–0.7084) is in agreement with a palaeo‐seawater isotopic signature; (iii) measured δ¹¹B values for the deepest brine samples are enriched in ¹¹B (up to 36.1‰) relative to the granitic springs. The combined use of δ¹¹B, Cl, B, Br, Sr contents and ⁸⁷Sr/⁸⁶Sr ratios makes it possible to define and quantify a mixing model between marine and crustal end‐members in order to explain the origin of the deep saline groundwaters in the Vienne granitic rocks.     

Geochemical fingerprinting of hydraulic fracturing fluids from Qusaiba Hot Shale and formation water from Paleozoic petroleum systems,Saudi Arabia

Provenance studies of produced water are essential to trace flow dynamics and reservoir compartmentalization in petroleum systems and to quantify fluid recovery rates from unconventional fracturing. Produced water from a hydraulically fractured well in the Qusaiba Hot Shale in the Northern Exploration Area, Saudi Arabia, was daily monitored and analyzed for water chemistry, and environmental (δ²H, δ¹³C, δ¹⁸O_H2O, δ¹⁸O_SO4, δ³⁴S_SO4, δ³⁷Cl, ⁸⁷Sr/⁸⁶Sr) and cosmogenic isotopes (³H, ¹⁴C, ³⁶Cl), to differentiate from reference fluids of supply water, fracturing fluids, and formation water from adjacent Paleozoic units. Initially, recovered water is composed of fracturing fluids and subsequently replaced by a homogeneous cut of pristine formation water. Formation water is composed of dominant meteoric water (approximately 84 vol%) and minor fossil evaporated seawater. The young ¹⁴C‐apparent age between 6000 and 6700 years BP and depleted δ¹⁸O/δ²H values for the meteoric component confirm the infiltration of surface water into the Qusaiba Hot Shale interval or adjacent units during the Early Holocene Pluvial Period under cooler and wetter climatic conditions than present, which suggest the presence of a very recent, dynamic hydraulic flow system. ³⁶Cl/Cl ratios between 102 × 10^?15 and 31 × 10^?15 are ambiguous and can be attributed to atmospheric recharge close to the coast, mixing of ³⁶Cl‐enriched Quaternary meteoric recharge with ³⁶Cl‐depleted fossil seawater, and/or hypogene production by U‐Th‐enriched host rock. Produced waters from Qusaiba Hot Shale are within the compositional range of Na‐Cl‐type formation water from Paleozoic reservoir units in northern Saudi Arabia with salinities from 30 000 to 130 000 mg l^?1. As a novel technological approach for exploration wells in Northern Saudi Arabia, multi‐isotopic methods were successfully implemented to quantify flowback volumes from hydraulic fracturing, and to fingerprint pristine formation water or pore water in Paleozoic systems on their provenance, residence time, migration pathways, and secondary alteration processes.     

Mobilization of ore fluids during Alpine metamorphism: evidence from hydrothermal veins in the Variscan basement of Western Carpathians,Slovakia

Hydrothermal polymetallic veins of the Gemeric unit of the Western Carpathians are oriented coherently with the foliation of their low‐grade Variscan basement host. Early siderite precipitated from homogeneous NaCl‐KCl‐CaCl₂‐H₂O brines with minor CO₂, while immiscible gas–brine mixtures are indicative of the superimposed barite, quartz–tourmaline and quartz–sulphide stages. The high‐salinity aqueous fluid (18–35 wt%) found in all mineralization stages corresponds to formation water modified by interaction with crystalline basement rocks at temperatures between 140 and 300°C. High brominity (around 1000 ppm in average) resulted from evaporation and anhydrite precipitation in a Permo‐Triassic marine basin, and from secondary enrichment by dissolution of organic matter in the marine sediments at diagenetic temperatures. Sulphate depletion reflects thermogenic reduction during infiltration of the formation waters into the Variscan crystalline basement. Crystallization temperatures of the siderite fill (140–300°C) and oxygen isotope ratios of the parental fluids (4–10‰) increase towards the centre of the Gemeric cleavage fan, probably as a consequence of decreasing water/rock ratios in rock‐buffered hydrothermal systems operating during the initial stages of vein evolution. In contrast, buoyant gas–water mixtures, variable salinities and strongly fluctuating P–T parameters in the successive mineralization stages reflect transition from a closed to an open hydrothermal system and mixing of fluids from various sources. Depths of burial were 6–14 km (1.7–4.4 kbar, in a predominantly lithostatic fluid regime) during the siderite and barite sub‐stages of the north‐Gemeric veins, and up to 16 km (1.6–4.5 kbar, in a hydrostatic to lithostatic fluid regime) in the quartz–tourmaline stage of the south‐Gemeric veins. The fluid pressure decreased down to approximately 0.6 kbar during crystallization of sulphides. U‐Pb‐Th, ⁴⁰Ar/³⁹Ar and K/Ar geochronology applied to hydrothermal muscovite–phengite and monazite, as well as cleavage phyllosilicates in the adjacent basement rocks and deformed Permian conglomerates corroborated the opening of hydrothermal veins during Lower Cretaceous thrusting and their rejuvenation during Late Cretaceous sinistral transpressive shearing and extension.     

Numerical model of pore‐pressure diffusion associated with the initiation of the 2010–2011 Guy–Greenbrier,Arkansas earthquakes

We model pore‐pressure diffusion caused by pressurized waste‐fluid injection at two nearby wells and then compare the buildup of pressure with the observed initiation and migration of earthquakes during the early part of the 2010–2011 Guy–Greenbrier earthquake swarm. Pore‐pressure diffusion is calculated using MODFLOW 2005 that allows the actual injection histories (volume/day) at the two wells to diffuse through a fractured and faulted 3D aquifer system representing the eastern Arkoma basin. The aquifer system is calibrated using the observed water‐level recovery following well shut‐in at three wells. We estimate that the hydraulic conductivities of the Boone Formation and Arbuckle Group are 2.2 × 10^?2 and 2.03 × 10^?3 m day^?1, respectively, with a hydraulic conductivity of 1.92 × 10^?2 m day^?1 in the Hunton Group when considering 1.72 × 10^?3 m day^?1 in the Chattanooga Shale. Based on the simulated pressure field, injection near the relatively conductive Enders and Guy–Greenbrier faults (that hydraulically connect the Arbuckle Group with the underlying basement) permits pressure diffusion into the crystalline basement, but the effective radius of influence is limited in depth by the vertical anisotropy of the hydraulic diffusivity. Comparing spatial/temporal changes in the simulated pore‐pressure field to the observed seismicity suggests that minimum pore‐pressure changes of approximately 0.009 and 0.035 MPa are sufficient to initiate seismic activity within the basement and sedimentary sections of the Guy–Greenbrier fault, respectively. Further, the migration of a second front of seismicity appears to follow the approximately 0.012 MPa and 0.055 MPa pore‐pressure fronts within the basement and sedimentary sections, respectively.     

The geochemical evolution of very dilute CO2‐rich water in Chungcheong Province,Korea: processes and pathways

A geochemical study was carried out on the CO₂‐rich water occurring in granite areas of Chungcheong Province, Korea. In this area, very dilute and acidic CO₂‐rich waters [62–242 mg l^?1 in total dissolved solid (TDS), 4.0–5.3 in pH; group I) occur together with normal CO₂‐rich waters (317–988 mg l^?1 in TDS, 5.5–6.0 in pH; group II). The concentration levels and ages of group I water are similar to those of recently recharged and low‐mineralized groundwater (group III). Calculation of reaction pathways suggests that group I waters are produced by direct influx of CO₂ gas into group III type waters. When the groundwater is injected with CO₂, it develops the capacity to accept dissolved solids and it can evolve into water with very high solute concentrations. Whether the water is open or closed to the CO₂ gases becomes less important in controlling the reaction pathway of the CO₂‐rich groundwater when the initial pco ₂ is high. Our data show that most of the solutes are dissolved in the CO₂‐rich groundwater at pH > 5 where the weathering rates of silicates are very slow or independent of pH. Thus, groundwater age is likely more important in developing high solute concentrations in the CO₂‐rich groundwaters than accelerated weathering kinetics because of acidic pH caused by high pco ₂.     

Exploring the potential linkages between oil‐field brine reinjection,crystalline basement permeability,and triggered seismicity for the Dagger Draw Oil field,southeastern New Mexico,USA, using hydrologic modeling

We used hydrologic models to explore the potential linkages between oil‐field brine reinjection and increases in earthquake frequency (up to M_d 3.26) in southeastern New Mexico and to assess different injection management scenarios aimed at reducing the risk of triggered seismicity. Our analysis focuses on saline water reinjection into the basal Ellenburger Group beneath the Dagger Draw Oil field, Permian Basin. Increased seismic frequency (>M_d 2) began in 2001, 5 years after peak injection, at an average depth of 11 km within the basement 15 km to the west of the reinjection wells. We considered several scenarios including assigning an effective or bulk permeability value to the crystalline basement, including a conductive fault zone surrounded by tighter crystalline basement rocks, and allowing permeability to decay with depth. We initially adopted a 7 m (0.07 MPa) head increase as the threshold for triggered seismicity. Only two scenarios produced excess heads of 7m five years after peak injection. In the first, a hydraulic diffusivity of 0.1 m² s^?1 was assigned to the crystalline basement. In the second, a hydraulic diffusivity of 0.3 m² s^?1 was assigned to a conductive fault zone. If we had considered a wider range of threshold excess heads to be between 1 and 60 m, then the range of acceptable hydraulic diffusivities would have increased (between 0.1–0.01 m² s^?1 and 1–0.1 m² s^?1 for the bulk and fault zone scenarios, respectively). A permeability–depth decay model would have also satisfied the 5‐year time lag criterion. We also tested several injection management scenarios including redistributing injection volumes between various wells and lowering the total volume of injected fluids. Scenarios that reduced computed excess heads by over 50% within the crystalline basement resulted from reducing the total volume of reinjected fluids by a factor of 2 or more.