Hydraulic and hydrochemical properties of deep sedimentary reservoirs of the Upper Rhine Graben,Europe

Hydraulic and hydrochemical data from several hundred wells mostly drilled by the oil and gas industry within the four deep carbonate and siliciclastic reservoirs of the Upper Rhine Graben area in France and Germany have been compiled, examined, validated and analysed with the aim to characterize fluids and reservoir properties. Due to enhanced temperatures in the subsurface of the Upper Rhine Graben, this study on hydraulic and hydrochemical properties has been motivated by an increasing interest in deep hydrogeothermal energy projects in the Rhine rift valley. The four examined geothermal reservoir formations are characterized by high hydraulic conductivity reflecting the active tectonic setting of the rift valley and its fractured and karstified reservoirs. The hydraulic conductivity decreases only marginally with depth in each of the reservoirs, because the Upper Rhine Graben is a young tectonically active structure. The generally high hydraulic conductivity of the reservoir rocks permits cross‐formation advective flow of thermal water. Water composition data reflect the origin and hydrochemical evolution of deep water. Shallow water to 500 m depth is, in general, weakly mineralized. The chemical signature of the water is controlled by fluid–rock geochemical interactions. With increasing depth, the total of dissolved solids (TDS) increases. In all reservoirs, the fluids evolve to a NaCl‐dominated brine. The high salinity of the reservoirs is partly derived from dissolution of halite in evaporitic Triassic and Cenozoic formations, and partly from the fluids residing in the crystalline basement. Water of all four reservoirs is saturated with respect to calcite and other minerals including quartz and barite.     

Permeability of the continental crust: dynamic variations inferred from seismicity and metamorphism

The variation of permeability with depth can be probed indirectly by various means, including hydrologic models that use geothermal data as constraints and the progress of metamorphic reactions driven by fluid flow. Geothermal and metamorphic data combine to indicate that mean permeability ( k ) of tectonically active continental crust decreases with depth ( z ) according to log  k  ≈ −14–3.2 log  z , where k is in m² and z in km. Other independently derived, crustal-scale k – z relations are generally similar to this power-law curve. Yet there is also substantial evidence for local-to-regional-scale, transient, permeability-generation events that entail permeabilities much higher than these mean k – z relations would suggest. Compilation of such data yields a fit to these elevated, transient values of log  k  ≈ −11.5–3.2 log  z , suggesting a functional form similar to that of tectonically active crust, but shifted to higher permeability at a given depth. In addition, it seems possible that, in the absence of active prograde metamorphism, permeability in the deeper crust will decay toward values below the mean k – z curves. Several lines of evidence suggest geologically rapid (years to 10³ years) decay of high-permeability transients toward background values. Crustal-scale k – z curves may reflect a dynamic competition between permeability creation by processes such as fluid sourcing and rock failure, and permeability destruction by processes such as compaction, hydrothermal alteration, and retrograde metamorphism.     

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.     

Deep geothermal groundwater flow in the Seferihisar–Balçova area,Turkey: results from transient numerical simulations of coupled fluid flow and heat transport processes

The Seferihisar–Balçova Geothermal system (SBG) is characterized by complex temperature and hydrochemical anomalies. Previous geophysical and hydrochemical investigations suggest that hydrothermal convection in the faulted areas of the SBG and recharge flow from the Horst may be responsible for the observed patterns. A numerical model of coupled fluid flow and heat transport processes has been built in order to study the possible fluid dynamics of deep geothermal groundwater flow in the SBG. The results support the hypothesis derived from interpreted data. The simulated scenarios provide a better understanding of the geophysical conditions under which the different fluid dynamics develop. When recharge processes are weak, the convective patterns in the faults can expand to surrounding reservoir units or below the seafloor. These fault‐induced drag forces can cause natural seawater intrusion. In the Melange of the Seferihisar Horst, the regional flow is modified by buoyant‐driven flow focused in the series of vertical faults. As a result, the main groundwater divide can shift. Sealing caprocks prevent fault‐induced cells from being overwhelmed by vigorous regional flow. In this case, over‐pressured, blind geothermal reservoirs form below the caprocks. Transient results showed that the front of rising hot waters in faults is unstable: the tip of the hydrothermal plumes can split and lead to periodical temperature oscillations. This phenomenon known as Taylor–Saffman fingering has been described in mid‐ocean ridge hydrothermal systems. Our findings suggest that this type of thermal pulsing can also develop in active, faulted geothermal systems. To some extent, the role of an impervious fault core on the flow patterns has been investigated. Although it is not possible to reproduce basin‐scale transport processes, this first attempt to model deep groundwater geothermal flow in the SBG qualitatively supported the interpreted data and described the different fluid dynamics of the basin. Geofluids (2010) 10 , 388–405     

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.     

Evolution of porosity, permeability and fluid pressure in dilatant faults post-failure: implications for fluid flow and mineralization

Mineralization of brittle fault zones is associated with sudden dilation, and the corresponding changes in porosity, permeability and fluid pressure, that occur during fault slip events. The resulting fluid pressure gradients cause fluid to flow into and along the fault until it is sealed. The volume of fluid that can pass through the deforming region depends on the degree of dilation, the porosity and permeability of the fault and wall rocks, and the rate of fault sealing. A numerical model representing a steep fault cutting through a horizontal seal is used to investigate patterns of fluid flow following a dilatant fault slip event. The model is initialized with porosity, permeability and fluid pressure representing the static mechanical state of the system immediately after such an event. Fault sealing is represented by a specified evolution of porosity, coupled to changes in permeability and fluid pressure, with the rate of porosity reduction being constrained by independent estimates of the rate of fault sealing by pressure solution. The general pattern of fluid flow predicted by the model is of initial flow into the fault from all directions, followed by upward flow driven by overpressure beneath the seal. The integrated fluid flux through the fault after a single failure event is insufficient to account for observed mineralization in faults; mineralization would require multiple fault slip events. Downward flow is predicted if the wall rocks below the seal are less permeable than those above. This phenomenon could at least partially explain the occurrence of uranium deposits in reactivated basement faults that cross an unconformity between relatively impermeable basement and overlying sedimentary rocks.     

Salt tectonics and shallow subseafloor fluid convection: models of coupled fluid‐heat‐salt transport

Thermohaline convection associated with salt domes has the potential to drive significant fluid flow and mass and heat transport in continental margins, but previous studies of fluid flow associated with salt structures have focused on continental settings or deep flow systems of importance to petroleum exploration. Motivated by recent geophysical and geochemical observations that suggest a convective pattern to near‐seafloor pore fluid flow in the northern Gulf of Mexico (GoMex), we devise numerical models that fully couple thermal and chemical processes to quantify the effects of salt geometry and seafloor relief on fluid flow beneath the seafloor. Steady‐state models that ignore halite dissolution demonstrate that seafloor relief plays an important role in the evolution of shallow geothermal convection cells and that salt at depth can contribute a thermal component to this convection. The inclusion of faults causes significant, but highly localized, increases in flow rates at seafloor discharge zones. Transient models that include halite dissolution show the evolution of flow during brine formation from early salt‐driven convection to later geothermal convection, characteristics of which are controlled by the interplay of seafloor relief and salt geometry. Predicted flow rates are on the order of a few millimeters per year or less for homogeneous sediments with a permeability of 10^?15 m², comparable to compaction‐driven flow rates. Sediment permeabilities likely fall below 10^?15 m² at depth in the GoMex basin, but such thermohaline convection can drive pervasive mass transport across the seafloor, affecting sediment diagenesis in shallow sediments. In more permeable settings, such flow could affect methane hydrate stability, seafloor chemosynthetic communities, and the longevity of fluid seeps.     

Modeling fluid flow in a heterogeneous, fault-controlled hydrothermal system

Previous studies have shown that most hydrothermal systems discharging at the land surface are associated with faulting, and that the location, temperature and rate of discharge of these systems are controlled by the geometry and style of the controlling fault(s). Unfortunately, the transport of heat and fluid in fault-controlled hydrothermal systems is difficult to model realistically; although heterogeneity and anisotropy are assumed to place important controls on flow in faults, few data or observations are available to constrain the distribution of hydraulic properties within active faults. Here, analytical and numerical models are combined with geostatistical models of spatially varying hydraulic properties to model the flow of heat and fluid in the Borax Lake fault of south-east Oregon, USA. A geometric mean permeability within the fault of 7 × 10⁻¹⁴ m² with 2× vertical/horizontal anisotropy in correlation length scale is shown to give the closest match to field observations. Furthermore, the simulations demonstrate that continuity of flow paths is an important factor in reproducing the observed behavior. In addition to providing some insight into possible spatial distributions of hydraulic properties at the Borax Lake site, the study highlights one potential avenue for integrating field observations with simulation results in order to gain greater understanding of fluid flow in faults and fault-controlled hydrothermal and petroleum reservoirs.     

Linked thermal convection of the basement and basinal fluids in formation of the unconformity‐related uranium deposits in the Athabasca Basin,Saskatchewan, Canada

World‐class unconformity‐related U deposits in the Athabasca Basin (Saskatchewan, Canada) are generally located within or near fault zones that intersect the unconformity between the Athabasca Group sedimentary basin rocks and underlying metamorphic basement rocks. Two distinct subtypes of unconformity‐related uranium deposits have been identified: those hosted primarily in the Athabasca Group sandstones (sediment‐hosted) and those hosted primarily in the underlying basement rocks (basement‐hosted). Although significant research on these deposits has been carried out, certain aspects of their formation are still under discussion, one of the main issues being the fluid flow mechanisms responsible for uranium mineralization. The intriguing feature of this problem is that sediment‐hosted and basement‐hosted deposits are characterized by oppositely directed vectors of fluid flow via associated fault zones. Sediment‐hosted deposits formed via upward flow of basement fluids, basement‐hosted deposits via downward flow of basinal fluids. We have hypothesized that such flow patterns are indicative of the fluid flow self‐organization in fault‐bounded thermal convection (Transport in Porous Media, 110, 2015, 25). To explore this hypothesis, we constructed a simplified hydrogeologic model with fault‐bounded thermal convection of fluids in the faulted basement linked with fluid circulation in the overlying fault‐free sandstone horizon. Based on this model, a series of numerical experiments was carried out to simulate the hypothesized fluid flow patterns. The results obtained are in reasonable agreement with the concept of fault‐bounded convection cells as an explanation of focused upflow and downflow across the basement/sandstone unconformity. We then discuss application of the model to another debated problem, the uranium source for the ore‐forming basinal brines.     

Geophysical and geochemical evidence of large scale fluid flow within shallow sediments in the eastern Gulf of Mexico,offshore Louisiana

We analyse the fluid flow regime within sediments on the Eastern levee of the modern Mississippi Canyon using 3D seismic data and downhole logging data acquired at Sites U1322 and U1324 during the 2005 Integrated Ocean Drilling Program (IODP) Expedition 308 in the Gulf of Mexico. Sulphate and methane concentrations in pore water show that sulphate–methane transition zone, at 74 and 94 m below seafloor, are amongst the deepest ever found in a sedimentary basin. This is in part due to a basinward fluid flow in a buried turbiditic channel (Blue Unit, 1000 mbsf), which separates sedimentary compartments located below and above this unit, preventing normal upward methane flux to the seafloor. Overpressure in the lower compartment leads to episodic and focused fluid migration through deep conduits that bypass the upper compartment, forming mud volcanoes at the seabed. This may also favour seawater circulation and we interpret the deep sulphate–methane transition zones as a result of high downward sulphate fluxes coming from seawater that are about 5–10 times above those measured in other basins. The results show that geochemical reactions within shallow sediments are dominated by seawater downwelling in the Mars‐Ursa basin, compared to other basins in which the upward fluid flux is controlling methane‐related reactions. This has implications for the occurrence of gas hydrates in the subsurface and is evidence of the active connection between buried sediments and the water column.     

Hydrothermal convection in and around mineralized fault zones: insights from two- and three-dimensional numerical modeling applied to the Ashanti belt, Ghana

Among hydrogeological processes, free convection in faults has been cited as a possible cause of gold mineralization along major fault zones. Here, we investigate the effects of free convection to determine whether it can cause giant orogenic gold deposits and their regular spatial distribution along major fault/shear zones. The approach comprises: (i) coupled two- and three-dimensional numerical heat- and fluid-flow simulations of simplified geological models; and (ii) calculation of the rock alteration index (RAI) to delineate regions where precipitation/dissolution can occur. Then, comparing the deduced alteration patterns with temperature distribution, potential areas of gold mineralization, defined by T  > 200°C and RAI < 0, are predicted. The models are based on the orogenic Paleoproterozoic ore deposits of the Ashanti belt in western Africa. These deposits occur in the most permeable parts of the fault zone, where the lateral permeability contrast is the highest. For a simple geometry, with a fault zone adjacent to a sedimentary basin half as permeable, we note a transition from three-dimensional circulation within the fault to a two-dimensional convective pattern in the basin far from the fault. Moreover, whereas two-dimensional undulated isotherms dominate in the basin, three-dimensional corrugated isotherms result from the preferred convective pattern within the fault, thus enhancing a periodic distribution of thermal highs and lows. In our most elaborate three-dimensional model with an imposed lateral permeability gradient, the RAI distribution indicates that fluid circulation in fault zones gives rise to a spatial periodicity of alteration patterns consistent with field data.     

Episodic fluid flow within continental rift basins: some insights from field data and mathematical models of the Rhinegraben

Vitrinite reflectance data from a petroleum exploration well in the northern Upper Rhinegraben show an unusual vertical maturity trend. Above and below a 500 m thick marl layer the vitrinite reflectance levels are consistent with modern, conductive, geothermal gradients. Between about 1000 and 1500 m depth, however, vitrinite reflectance levels are significantly elevated (about 0.6%Ro). This anomaly cannot be explained with one‐dimensional conductive or conductive–convective heat transfer models, and thermal effects of sedimentation or igneous intrusion seem implausible for this geological setting. The thermal anomaly that formed this maturation anomaly must have been hydrothermal in origin, two‐dimensional in nature, and persisted long enough to elevate the vitrinite reflectance values within this marl unit, yet it must have dissipated before the thermal perturbation would have altered the organic matter below and above the unit. In this study, we propose that the vitrinite reflectance anomalies were caused by a transient thermal inversion induced by episodic, lateral flow of hot (130–160°C) groundwater along conductive fractures and bedding planes. Heat flow constraints suggest that fluids must have moved rapidly up a vertical feeder fault from a depth of at least 3.6 km before migrating laterally. To test this hypothesis, we present a suite of simple, idealized mathematical models of groundwater flow, heat transfer, thermal degradation of kerogen and vitrinite systematics to explore the episodic flow that could have produced the observed thermal anomaly. In these simulations, a single, horizontal aquifer is sandwiched between two less permeable units: the total dimensions of the vertical section model are 4 km thick by 10 km long. The top of the aquifer coincides with the position of the observed thermal maturity anomaly in the Rhinegraben. Boundary conditions along the left edge of this aquifer were varied through time to allow for the migration of hot fluids out into the basin. Inflow temperature, horizontal velocity, duration and frequency of flow and thickness of the aquifer were varied. We found that a thermal maturity anomaly could only be produced by a rather restrictive set of hydrothermal conditions. It was possible to produce the observed vitrinite reflectance anomaly by a single hydrothermal flow event of 130°C fluid migrating laterally into the aquifer at a rate of 1 m a^?1 for about 10 000 years. The anomaly is spatially confined to near the left edge of the basin, near the feeder fault. If the flow event lasted longer than 100 000 years, then the maturation anomaly disappeared as the lower confining unit approached steady‐state thermal conditions. It is possible that such an event occurred about 5 million years ago in response to increases in fault permeability associated with far field Alpine tectonism.     

Computing geochemical mass transfer and water/rock ratios in submarine hydrothermal systems: implications for estimating the vigour of convection

A method of calculating chemical water/rock ratios is presented that enables the estimation of fluid velocities in open, flow‐through hydrologic systems. The approach is based on relating the gain/loss of a chemical species per kilogram of solid phase to the loss/gain of that species in the fluid phase, integrated across a specified length of the flowpath. After examining the underlying approximations of the approach using a one‐dimensional model of seawater moving through a basalt under nonisothermal conditions, the method is applied to representative zones within a two‐dimensional hydrothermal convective system. The method requires that regions within the flow system can be identified in which the direction of flow is steady for an extended period of time. Estimates of fluid velocity are spatial and temporal averages for the length of the flowpath used in the calculation. The location within the flow system and the nature of the alteration reactions determine which species can provide reliable values of the chemical water/rock ratio and useful estimates of fluid velocities. Over the length of the flowpath considered, the calculation of water/rock ratios works best when a species is controlled by a single reaction. Accurate estimates are favoured if the concentration profile of a species along the flowpath increases or decreases monotonically. If the length of the flowpath extends over more than one reaction zone, then erroneous estimates of the water/rock ratio and fluid velocity are more likely. Model calculations suggest that the quartz/silica system should provide reliable estimates for fluid velocity under a wide range of temperature and flow conditions, in particular in those regions of a system at or near quartz equilibrium, so that the aqueous silica concentration is buffered by quartz and correlated with the temperature distribution.     

Landscape and Soil Development on Monaro Basalt West of Nimmitabel,New South Wales

Prolonged stability of the southern tableland of NSW, developed on early Tertiary Monaro basalt, is evident from weak landscape dissection and numerous small lakes. However, poorly developed soils and weathered zones are inconsistent with prolonged Cainozoic weathering or formation under climates substantially different to the present. Thin Holocene sediments are also at variance with the sizes of the lake basins. Low depositional landforms, termed ‘lake shadows’, occur on the eastern perimeters of most lakes, and comprise clay pellets blown from basin floors during dry phases. Deflation, dominantly during Quaternary glacial periods, and reactivated today, is responsible for the thin soils and sediment cover. The lake basins are partly solutional and partly deflational features. The modern soil cover reflects the modern climate.     

Maintaining high fluid pressures in older sedimentary basins

High fluid pressures in old geologic basins, where the mechanisms that generate high fluid pressure have ceased to operate, pose the problem of how high fluid pressures are maintained through geologic time. Recent oil and gas exploration reveals that low permeability shales, the source beds for oil and gas, contain large quantities of gas that are now being exploited in many sedimentary basins in North America. No earlier analyses of how to maintain high fluid pressure in older sedimentary basins included a shale bed as a source of adsorbed gas; this is a new conceptual element that will fundamentally change the analysis. Such a large fluid source can compensate for a low rate of bleed off in a dynamic system. If the fluid source is large enough, as the gas within these shale source beds appears to be, there will no appreciable drop in pressure accompanying a low rate of leakage from the basin for long periods. For the dynamic school of basin analysts this may provide the missing piece in the puzzle, explaining how high fluid pressures are maintained for long periods of geologic time in a crust with finite, non-zero permeability. This is a hypothesis which needs to be tested by new basin analyses.     

Fluid flow in shear zones: insights from the geometry and evolution of ore bodies at Renco gold mine,Zimbabwe

The geometry of mineral deposits can give insights into fluid flow in shear zones. Lode gold ore bodies at Renco Mine, in the Limpopo Belt, Zimbabwe, occur as siliceous breccias and mylonites within amphibolite facies shear zones that dip either gently or steeply. The two sets of ore bodies formed synchronously from hydrothermal fluids. The ore bodies are oblate, but have well‐defined long axes. Larger ore bodies are more oblate. High‐grade gold ore shoots have long axes that plunge down dip; this direction is perpendicular to the long axes of the low‐grade ore bodies. The centres of the high‐grade ore bodies align within the low‐grade ore bodies along strike in both gently and steeply dipping groups. The range of sizes and shapes of the ore bodies are interpreted as a growth sequence. Geometrical models are proposed for the gently and steeply dipping ore bodies, in which individual ore bodies grow with long axes plunging down dip, and merge to form larger, more oblate ore bodies. The models show that when three or more ore bodies coalesce, the long axis of the merged ore body is perpendicular to the component ore bodies, and that ore bodies in the deposit may have a range of shapes due to both growth of individual ore bodies, and their coalescence. The long axes of the high‐grade ore bodies are parallel to the shear directions of both the gently and steeply dipping dip slip shear zones, which were the directions of greatest permeability and fluid flow. The larger, lower grade bodies, which may have formed by coalescence, are elongate perpendicular to these directions.     

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.     

Hydrologic and geochemical controls on soluble benzene migration in sedimentary basins

The effects of groundwater flow and biodegradation on the long‐distance migration of petroleum‐derived benzene in oil‐bearing sedimentary basins are evaluated. Using an idealized basin representation, a coupled groundwater flow and heat transfer model computes the hydraulic head, stream function, and temperature in the basin. A coupled mass transport model simulates water washing of benzene from an oil reservoir and its miscible, advective/dispersive transport by groundwater. Benzene mass transfer at the oil–water contact is computed assuming equilibrium partitioning. A first‐order rate constant is used to represent aqueous benzene biodegradation. A sensitivity study is used to evaluate the effect of the variation in aquifer/geochemical parameters and oil reservoir location on benzene transport. Our results indicate that in a basin with active hydrodynamics, miscible benzene transport is dominated by advection. Diffusion may dominate within the cap rock when its permeability is less than 10^?19 m². Miscible benzene transport can form surface anomalies, sometimes adjacent to oil fields. Biodegradation controls the distance of transport down‐gradient from a reservoir. We conclude that benzene detected in exploration wells may indicate an oil reservoir that lies hydraulically up‐gradient. Geochemical sampling of hydrocarbons from springs and exploration wells can be useful only when the oil reservoir is located within about 20 km. Benzene soil gas anomalies may form due to regional hydrodynamics rather than separate phase migration. Diffusion alone cannot explain the elevated benzene concentration observed in carrier beds several km away from oil fields.     

The influence of poorly interconnected fault zone flow paths on spring geochemistry

Thermal springs commonly occur along faults because of the enhanced vertical permeability afforded by fracture zones. Field and laboratory studies of fault zone materials document substantial heterogeneities in fracture permeabilities. Modeling and field studies of springs suggest that spatial variations in permeability strongly influence spring locations, discharge rates and temperatures. The impact of heterogeneous permeability on spring geochemistry, however, is poorly documented. We present stable isotope and water chemistry data from a series of closely spaced thermal springs associated with the Hayward Fault, California. We suggest that substantial spatial variations observed in δ¹⁸O and chloride values reflect subsurface fluid transport through a poorly connected fracture network in which mixing of subsurface waters remains limited. Our measurements provide insight into the effect of fracture zone heterogeneities on spring geochemistry, offer an additional tool to intuit the nature of tectonically induced changes in fault zone plumbing, and highlight the need to consider local variations when characterizing fracture zone fluid geochemistry from spring systems with multiple discharge sites.     

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.