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.     

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.     

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.     

Hydrothermal fluid flow models of stratiform ore genesis in the McArthur Basin, Northern Territory, Australia

This paper explores the role of basin‐scale fluid migration in stratiform Pb–Zn ore formation in the southern McArthur Basin, Australia. Mathematical models are presented for coupled brine migration and heat transport in the basin. The models account for: (i) topographically driven flow (forced convection) during periods when parts of the McArthur Basin were subaerial and elevated above the central Batten Fault Zone; (ii) density‐driven flow (free convection) during periods when the basin was mostly submarine; and (iii) transient flows associated with fault rupture during periods of transpression. These hydrologic models help to compare and contrast a variety of hypotheses concerning deep fluid migration and the origin of base metal ores in the McArthur Basin. The numerical results exhibit a strong structural control on fluid flow caused by the north‐trending fault systems that characterize the Batten Fault Zone. As a result, fluids descend to depths of a few kilometers along the western side, migrate laterally to the east through the clastic and volcanic aquifers of the upper Tawallah and lowest McArthur Groups, and then ascend along the eastern side of the fault zone. This recharge–discharge pattern dominates all of the hydrogeologic models. The basin‐wide flow pattern suggests that Na–Ca–Cl brines acquired base metals in the deepest levels of the basin stratigraphy as the fluids migrated eastwards through the aquifer system. Upward flow was relatively rapid along the Emu Fault Zone, so much so that fluid temperatures likely approached 130°C in the muddy sediments near the sea floor due to upward flow and venting at the HYC (‘Here’s Your Chance'). Transient pulses of flow characterized periods of transpressional stress and subsequent faulting may have punctuated the basin history. Large‐scale free convection, however, characterized notably long periods of diagenesis and ore mineralization during the Proterozoic in the McArthur Basin.     

Composition and evolution of fluids during skarn development in the Monte Capanne thermal aureole,Elba Island,central Italy

We describe the chemistry of the fluids circulating during skarn formation by focusing on fluids trapped in calcsilicate minerals of the inner thermal aureole of the Late Miocene Monte Capanne intrusion of western Elba Island (central Italy). Primary, CH₄‐dominant, C‐O‐H‐S‐salt fluid inclusions formed during prograde growth of the main skarn‐forming mineral phases: grossular/andradite and vesuvianite. The variable phase ratios attest to heterogeneous entrapment of fluid, with co‐entrapment of an immiscible hydrocarbon–brine mixture. Chemical elements driving skarn metasomatism such as Na, K, Ca, S and Cl, Fe and Mn were dominantly partitioned into the circulating fluid phase. The high salinity (apparent salinity between 58 and 70 wt% NaCl eq.) and the C‐component of the fluids are interpreted as evidence for a composite origin of the skarn‐forming fluids that involves both fluids derived from the crystallizing intrusion and contributions from metamorphic devolatilization. Oxidation of a Fe‐rich brine in an environment dominated by fluctuation in pressure from lithostatic to hydrostatic conditions (maintained by active crack‐sealing) contributed to skarn development. Fluid infiltration conformed to a geothermal gradient of about 100°C km^?1, embracing the transition from high‐temperature contact metamorphism and fluid‐assisted skarn formation (at ca 600°C) to a barren hydrothermal stage (at ca 200°C).     

Cryogenic vitrification and 3D serial sectioning using high resolution cryo‐FIB SEM technology for brine‐filled grain boundaries in halite: first results

The structure of brine films in grain boundaries of halite has been the subject of much controversy over the past 20 years; although a number of innovative methods have been developed to study these structures, much is still unknown and fundamental information is missing. In this study, we investigated different methods of plunge‐freezing to vitrify the brine fill of grain boundaries for natural salt polycrystal. This was followed by a preliminary study of the 3D morphology of a vitrified grain boundary in a natural rock salt sample with a focused ion beam (FIB) excavation system. We have shown that brine‐filled grain boundaries in rock salt can be efficiently well frozen when dimensions are less than about 1 mm. Coupled with an ion beam tool, cryo‐SEM allows 3D observation of the well‐frozen grain boundaries in large volumes and high resolution. Initial results of brine‐filled natural halite grain boundaries show non‐faceted crystal–brine interfaces and unexpectedly low dihedral angles at room temperature and pressure.     

Topography‐driven flow versus buoyancy‐driven flow in the U.S. midcontinent: implications for the residence time of brines

Topography‐driven flow is normally considered to be the dominant groundwater flow system in uplifted sedimentary basins. In the U.S. midcontinent region east of the Rocky Mountains, the presence of brines derived from dissolution of halite suggests that significant topography‐driven flushing has occurred to remove older brines that presumably formed concurrently with Permian evaporites in the basin. However, the presence of evaporites and brines in the modern basin suggests that buoyancy‐driven flow could limit topography‐driven flushing significantly. Here we used numerical models of variable‐density fluid flow, halite dissolution, solute transport, and heat transport to quantify flow patterns and brine migration. Results indicate the coexistence of large‐scale topography‐ and buoyancy‐driven flow. Buoyancy‐driven flow and low permeability evaporites act to isolate brines, and the residence time of the brines was found to be quite long, at least 50 Myr. The modern distribution of salinity appears to reflect near‐steady‐state conditions. Results suggest that flushing of original evaporatively‐concentrated brines occurred tens of millions of years ago, possibly concurrent with maximum uplift ca. 60 Ma. Simulations also suggest that buoyancy‐driven convection could drive chemical exchange with crystalline basement rocks, which could supply significant Ca²⁺, Sr²⁺, and metals to brines.     

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.     

The hydrogeochemistry of subsurface brines in and west of the Jordan–Dead Sea Transform fault

This study is based on 113 analyses of brines with Cl > 0.57 mol l^?1 (modern seawater), which were collected and analysed mostly during several decades of exploration for gas and oil in Israel. Based on critical evaluation of correlations of elements and ionic ratios and on spider patterns, six different brine events or source brines were identified in the Phanerozoic: the Triassic, Lower Cretaceous and the Mio/Pliocene brine families which were identified in boreholes Sdom‐1, Sdom Deep‐1 and Ha'on, and the Holocene Dead Sea brines. The Triassic brines are nowadays also encountered in under‐ and overlying rock units such as the Paleozoic Negev‐Yam Suf and the Jurassic Arad Groups, respectively. The southern Jordan–Dead Sea Transform (also known as the Rift) hosts the Mio‐Pliocene Sdom Deep and Sdom brine families. Brine bodies not sufficiently isolated by impervious sedimentary layers were flushed out during the Pliocene when the southern Valley drained north‐ and westwards through the Yizre'el Valley to the Mediterranean Sea. In the northern Rift Miocene to Pliocene seawater evaporated and infiltrated into the Rift sediments and into adjacent rocks. Further diluted by freshwater, it emerges as the Ha'on brine. Together with its derivatives, they form the Ha'on family. The derivatives of the Holocene Dead Sea brine family occur along the shoreline of the recent Dead Sea. Apart of all these evaporation brines, brines deriving from dissolution of evaporites locally occur in the area. The time‐bound chemical composition of paleoseawater is considered when discussing the ionic ratios of brines generated during different geological periods. Spider patterns of each brine family are compared and, where necessary, the relationship of brines to distinct families of brines is supported by inverse modelling.     

Hydrothermal,multiphase convection of H2O‐NaCl fluids from ambient to magmatic temperatures: a new numerical scheme and benchmarks for code comparison

Thermohaline convection of subsurface fluids strongly influences heat and mass fluxes within the Earth's crust. The most effective hydrothermal systems develop in the vicinity of magmatic activity and can be important for geothermal energy production and ore formation. As most parts of these systems are inaccessible to direct observations, numerical simulations are necessary to understand and characterize fluid flow. Here, we present a new numerical scheme for thermohaline convection based on the control volume finite element method (CVFEM), allowing for unstructured meshes, the representation of sharp thermal and solute fronts in advection‐dominated systems and phase separation of variably miscible, compressible fluids. The model is an implementation of the Complex Systems Modelling Platform CSMP++ and includes an accurate thermodynamic representation of strongly nonlinear fluid properties of salt water for magmatic‐hydrothermal conditions (up to 1000°C, 500 MPa and 100 wt% NaCl). The method ensures that all fluid properties are taken as calculated on the respective node using a fully upstream‐weighted approach, which greatly increases the stability of the numerical scheme. We compare results from our model with two well‐established codes, HYDROTHERM and TOUGH2, by conducting benchmarks of different complexity and find good to excellent agreement in the temporal and spatial evolution of the hydrothermal systems. In a simulation with high‐temperature, high‐salinity conditions currently outside of the range of both HYDROTHERM and TOUGH2, we show the significance of the formation of a solid halite phase, which introduces heterogeneity. Results suggest that salt added by magmatic degassing is not easily vented or accommodated within the crust and can result in dynamic, complex hydrologies.     

Reactive transport and thermo‐hydro‐mechanical coupling in deep sedimentary basins affected by glaciation cycles: model development,verification, and illustrative example

Deep sedimentary basins are complex systems that over long time scales may be affected by numerous interacting processes including groundwater flow, heat and mass transport, water–rock interactions, and mechanical loads induced by ice sheets. Understanding the interactions among these processes is important for the evaluation of the hydrodynamic and geochemical stability of geological CO₂ disposal sites and is equally relevant to the safety evaluation of deep geologic repositories for nuclear waste. We present a reactive transport formulation coupled to thermo‐hydrodynamic and simplified mechanical processes. The formulation determines solution density and ion activities for ionic strengths ranging from freshwater to dense brines based on solution composition and simultaneously accounts for the hydro‐mechanical effects caused by long‐term surface loading during a glaciation cycle. The formulation was implemented into the existing MIN3P reactive transport code (MIN3P‐THCm) and was used to illustrate the processes occurring in a two‐dimensional cross section of a sedimentary basin subjected to a simplified glaciation scenario consisting of a single cycle of ice‐sheet advance and retreat over a time period of 32 500 years. Although the sedimentary basin simulation is illustrative in nature, it captures the key geological features of deep Paleozoic sedimentary basins in North America, including interbedded sandstones, shales, evaporites, and carbonates in the presence of dense brines. Simulated fluid pressures are shown to increase in low hydraulic conductivity units during ice‐sheet advance due to hydro‐mechanical coupling. During the period of deglaciation, Darcy velocities increase in the shallow aquifers and to a lesser extent in deeper high‐hydraulic conductivity units (e.g., sandstones) as a result of the infiltration of glacial meltwater below the warm‐based ice sheet. Dedolomitization is predicted to be the most widespread geochemical process, focused near the freshwater/brine interface. For the illustrative sedimentary basin, the results suggest a high degree of hydrodynamic and geochemical stability.     

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

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

Fluid sinks within the earth's crust

The Urach 3 research borehole in south‐west (SW) Germany has been drilled through the sedimentary cover, and the gneisses of the Variscian crystalline basement at 1600 m below the surface (Black Forest basement) has been reached. An additional 2800 m has been drilled through the fractured crystalline rocks, and the borehole has been used for a number of hydraulic tests in the context of a ‘hot‐dry rock’ (HDR) project exploring for geothermal energy. The fracture system of the basement is saturated with a NaCl brine with about 70 g L^?1 dissolved solids. Water table measurements in the borehole cover a period of 13 years of observation, during which the water table continuously dropped and did not reach a steady‐state level. This unique set of data shows that the hydraulic potential decreases with depth, causing a continuous flow of fluid to the deeper parts of the upper continental crust. The potential decrease and the associated downward migration of fluid is an evidence for the progress of water (H₂O)‐consuming reactions in the crystalline rocks. Computed stability relations among relevant phases at the pressure temperature (PT) conditions in the fracture system and documented fossil fracture coatings in granites and gneisses suggest that the prime candidate for the H₂O‐consuming reaction is the zeolitization of feldspar. The potential of the gneisses to chemically bind H₂O matches the estimated amount of migrating H₂O.     

Methods for calculating brine evaporation rates during salt production

Salt is recognized by archaeologists as an important commodity due to the biological need for sodium and other cultural uses. Numerous studies have described the various techniques used in converting brine to crystallized salt, but few, if any, have attempted to quantify the physical processes of evaporation in pre-industrial societies. Apart from the few areas where salt mining is possible, nearly all forms of salt production require evaporation of water to concentrate brine and ultimately produce salt crystals. This study quantifies three of the most common evaporation techniques and provides insight into the production rates of salt and fuel requirements. Methods of calculation are provided for determining evaporation through (1) direct solar heating of brine, (2) applied external heat to a vessel, and (3) an immersed heated object (e.g., stone). These results provide physical constraints on the evaporation process and provide investigators with techniques for estimating efficiency and total production of prehistoric and historic saltworks.     

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     

Infiltration of basinal fluids into high-grade basement, South Norway: sources and behaviour of waters and brines

Quartz veins hosted by the high‐grade crystalline rocks of the Modum complex, Southern Norway, formed when basinal fluids from an overlying Palaeozoic foreland basin infiltrated the basement at temperatures of c. 220°C (higher in the southernmost part of the area). This infiltration resulted in the formation of veins containing both two‐phase and halite‐bearing aqueous fluid inclusions, sometimes with bitumen and hydrocarbon inclusions. Microthermometric results demonstrate a very wide range of salinities of aqueous fluids preserved in these veins, ranging from c. 0 to 40 wt% NaCl equivalent. The range in homogenization temperatures is also very large (99–322°C for the entire dataset) and shows little or no correlation with salinity. A combination of aqueous fluid microthermometry, halogen geochemistry and oxygen isotope studies suggest that fluids from a range of separate aquifers were responsible for the quartz growth, but all have chemistries comparable to sedimentary formation waters. The bulk of the quartz grew from relatively low δ¹⁸O fluids derived directly from the basin or equilibrated in the upper part of the basement (T < 200°C). Nevertheless, some fluids acquired higher salinities due to deep wall‐rock hydration reactions leading to salt saturation at high temperatures (>300°C). The range in fluid inclusion homogenization temperatures and densities, combined with estimates of the ambient temperature of the basement rocks suggests that at different times veins acted as conduits for influx of both hotter and colder fluids, as well as experiencing fluctuations in fluid pressure. This is interpreted to reflect episodic flow linked to seismicity, with hotter dry basement rocks acting as a sink for cooler fluids from the overlying basin, while detailed flow paths reflected local effects of opening and closing of individual fractures as well as reaction with wall rocks. Thermal considerations suggest that the duration of some flow events was very short, possibly in the order of days. As a result of the complex pattern of fracturing and flow in the Modum basement, it was possible for shallow fluids to penetrate basement rocks at significantly higher temperatures, and this demonstrates the potential for hydrolytic weakening of continental crust by sedimentary fluids.     

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.     

Fluid effect on hydraulic fracture propagation behavior: a comparison between water and supercritical CO2‐like fluid

The initiation of hydraulic fractures during fluid injection in deep formations can be either engineered or induced unintentionally. Upon injection of CO₂, the pore fluids in deep formations can be changed from oil/saline water to CO₂ or CO₂ dominated. The type of fluid is important not only because the fluid must fracture the rock, but also because rocks saturated with different pore fluids behave differently. We investigated the influence of fluid properties on fracture propagation behavior by using the cohesive zone model in conjunction with a poroelasticity model. Simulation results indicate that the pore pressure fields are very different for different pore fluids even when the initial field conditions and injection schemes (rate and time) are kept the same. Low viscosity fluids with properties of supercritical CO₂ will create relatively thin and much shorter fractures in comparison with fluids exhibiting properties of water under similar injection schemes. Two significant times are recognized during fracture propagation: the time at which a crack ceases opening and the later time point at which a crack ceases propagating. These times are very different for different fluids. Both fluid compressibility and viscosity influence fracture propagation, with viscosity being the more important property. Viscosity can greatly affect hydraulic conductivity and the leak‐off coefficient. This analysis assumes the in‐situ pore fluid and injected fluid are the same and the pore space is 100% saturated by that fluid at the beginning of the simulation.     

A simulation of the hydrothermal response to the Chesapeake Bay bolide impact

Groundwater more saline than seawater has been discovered in the tsunami breccia of the Chesapeake Bay Impact Crater. One hypothesis for the origin of this brine is that it may be a liquid residual following steam separation in a hydrothermal system that evolved following the impact. Initial scoping calculations have demonstrated that it is feasible such a residual brine could have remained in the crater for the 35 million years since impact. Numerical simulations have been conducted using the code HYDROTHERM to test whether or not conditions were suitable in the millennia following the impact for the development of a steam phase in the hydrothermal system. Hydraulic and thermal parameters were estimated for the bedrock underlying the crater and the tsunami breccia that fills the crater. Simulations at three different breccia permeabilities suggest that the type of hydrothermal system that might have developed would have been very sensitive to the permeability. A relatively low breccia permeability (1 × 10^?16 m²) results in a system partitioned into a shallow water phase and a deeper superheated steam phase. A moderate breccia permeability (1 × 10^?15 m²) results in a system with regionally extensive multiphase conditions. A relatively high breccia permeability (1 × 10^?14 m²) results in a system dominated by warm‐water convection cells. The permeability of the crater breccia could have had any of these values at given depths and times during the hydrothermal system evolution as the sediments compacted. The simulations were not able to take into account transient permeability conditions, or equations of state that account for the salt content of seawater. Results suggest, however, that it is likely that steam conditions existed at some time in the system following impact, providing additional evidence that is consistent with a hydrothermal origin for the crater brine.     

What do aqueous geothermometers really tell us?

The application of chemical geothermometry to shallow groundwaters or spring discharge assumes that there is minimal mixing or re-equilibration of water as it travels from depth to the surface. In this study, we examine the potential for mixing and re-equilibration by examining heat and fluid flow along crustal-scale faults in tectonic geothermal systems. Numerical modeling results indicate that maximum in situ temperatures could be under-predicted by up to 30% due to mixing of fluids that enter the fault at different depths. This, coupled with the depression of isotherms by downward groundwater flow in the hanging wall, could cause underestimates of maximum circulation depth of greater than 80% in extreme cases. Kinetics does not favor re-equilibration in the shallower portions of faults due to low temperatures and higher fluid velocities. However, in areas of deeper circulation or higher heat flow such reactions are possible.