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.     

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.     

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.     

Evaluation of recoverable energy potential from enhanced geothermal systems: a sensitivity analysis in a poro‐thermo‐elastic framework

This study presents application of an efficient approach to simulate fluid flow and heat transfer in naturally fractured geothermal reservoirs. Fluid flow is simulated by combining single continuum and discrete fracture approaches. The local thermal nonequilibrium approach is used to simulate heat transfer by conduction in the rock matrix and convection (including conduction) in the fluid. Fluid flow and heat transfer models are integrated within a coupled poro‐thermo‐elastic framework. The developed model is used to evaluate the long‐term response of a geothermal reservoir with specific boundary conditions and injection/production schedule. A comparative study and a sensitivity analysis are carried out to evaluate the capability of the integrated approach and understand the degree by which different reservoir parameters affect thermal depletion of Soultz geothermal reservoir, respectively. Also observed, there exists an optimum fracture permeability after which the reservoir stimulation does not change the recovery factor significantly. Estimation of fluid temperature by the assumption of local thermal nonequilibrium heat transfer between the fracture fluid and the rock matrix gives fluid temperature of about 3°C less than that of estimated by thermal equilibrium heat transfer at early stage of hot water production.     

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     

High‐temperature magmatic fluid exsolved from magma at the Duobuza porphyry copper–gold deposit,Northern Tibet

The Duobuza porphyry copper–gold deposit (proven Cu resources of 2.7 Mt, 0.94% Cu and 13 t gold, 0.21 g t^?1 Au) is located at the northern margin of the Bangong‐Nujiang suture zone separating the Qiangtang and Lhasa Terranes. The major ore‐bearing porphyry consists of granodiorite. The alteration zone extends from silicification and potassic alteration close to the porphyry stock to moderate argillic alteration and propylitization further out. Phyllic alteration is not well developed. Sericite‐quartz veins only occur locally. High‐temperature, high‐salinity fluid inclusions were observed in quartz phenocrysts and various quartz veins. These fluid inclusions are characterized by sylvite dissolution between 180 and 360°C and halite dissolution between 240 and 540°C, followed by homogenization through vapor disappearance between 620 and 960°C. Daughter minerals were identified by SEM as chalcopyrite, halite, sylvite, rutile, K–feldspar, and Fe–Mn‐chloride. They indicate that the fluid is rich in ore‐forming elements and of high oxidation state. The fluid belongs to a complex hydrothermal system containing H₂O – NaCl – KCl ± FeCl₂ ± CaCl₂ ± MnCl₂. With decreasing homogenization temperature, the fluid salinity tends to increase from 34 to 82 wt% NaCl equiv., possibly suggesting a pressure or Cl/H₂O increase in the original magma. No coexisting vapor‐rich fluid inclusions with similar homogenization temperatures were found, so the brines are interpreted to have formed by direct exsolution from magma rather than trough boiling off of a low‐salinity vapor. Estimated minimum pressure of 160 MPa imply approximately 7‐km depth. This indicates that the deposit represents an orthomagmatic end member of the porphyry copper deposit continuum. Two key factors are proposed for the fluid evolution responsible for the large size of the gold‐rich porphyry copper deposit of Duobuza: (i) ore‐forming fluids separated early from the magma, and (ii) the hydrothermal fluid system was of magmatic origin and highly oxidized.     

Hydrogeology of the Krafla geothermal system,northeast Iceland

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

Relationship of brines in the Kinnarot Basin,Jordan‐Dead Sea Rift Valley

Element ratios and water stable isotopes reveal the presence of only two independent deep brines in the Kinnarot Basin, Israel: the evaporite dissolution brine of Zemah‐1 and the inferred Ha’on mother brine (HMB) with low and high Br/Cl ratios, respectively. HMB is considered to be a representative of the Late Pliocene evaporated Sedom Sea. The freshwater‐diluted evaporation brine emerges as Ha’on brine on the eastern shore of Lake Tiberias and is also identified in the pore water of lake sediments. HMB is converted into Tiberias mother brine (TMB) by dolomitization of limestones and alteration of abundant volcanic rocks occurring along the western side of the lake. The Ha’on and Tiberias brines, both characterized by high δD and δ¹⁸O values, are similar in Na/Cl and Br/Cl ratios but are dissimilar in Br/K ratios because these brines were subjected to different degrees of interactions with rocks and sediments. Excepting the brine from KIN 8, all brines from the Tabigha area including the nearby off‐shore Barbutim brine are related to the TMB. The brine KIN 8 and all brines from the Fuliya and Hammat Gader areas are related to the HMB. The brine encountered in wildcat borehole Zemah‐1 is generated by halite‐anhydrite/gypsum dissolution and is independent from the HMB system.     

Pore pressure evolution and fluid flow during visco‐elastic single‐layer buckle folding

Pore pressure and fluid flow during the deformational history of geologic structures are directly influenced by tectonic deformation events. In this contribution, 2D plane strain finite element analysis is used to study the influence of different permeability distributions on the pore pressure field and associated flow regimes during the evolution of visco‐elastic single‐layer buckle folds. The buckling‐induced fluid flow regimes indicate that flow directions and, to a lesser degree, their magnitudes vary significantly throughout the deformation and as a function of the stratigraphic permeability distribution. The modelling results suggest that the volumetric strain and the permeability distribution significantly affect the resulting flow regime at different stages of fold development. For homogeneous permeability models (k > 10^?21 m²), low strain results in a mostly pervasive fluid flow regime and is in agreement with previous studies. For larger strain conditions, fluid focusing occurs in the buckling layer towards the top of the fold hinge. For low permeabilities (<10^?21 m²), local focused flow regimes inside the buckling layer emerge throughout the deformation history. For models featuring a low‐permeability layer embedded in a high‐permeability matrix or sandwiched between high‐permeability layers, focused flow regimes inside the folded layer result throughout the deformation history, but with significant differences in the flow vectors of the surrounding layers. Fluid flow vectors induced by the fold can result in different, even reversed, directions depending on the amount of strain. In summary, fluid flow regimes during single‐layer buckling can change from pervasive to focused and fluid flow vectors can be opposite at different strain levels, that is the flow vectors change significantly through time. Thus, a complete understanding of fluid flow regimes associated with single‐layer buckle folds requires consideration of the complete deformation history of the fold.     

Diverse origins and timing of formation of basinal brines in the Gulf of Mexico sedimentary basin

A review of five different field areas in the Gulf of Mexico sedimentary basin (GOM) illustrates some of the potentially diverse chemical and physical processes which have produced basinal brines. The elevated salinities of most of the formation waters in the GOM are ultimately related to the presence of the Middle Jurassic Louann Salt. Some of these brines likely inherited their salinity from evaporated Mesozoic seawater, while other saline fluids have been produced by subsequent dissolution of salt, some of which is occurring today. The timing of the generation of brines has thus not been restricted to the Middle Jurassic. The mechanisms of solute transport that have introduced brines throughout much of the sedimentary section of the GOM are not entirely understood. Free convection driven by spatial variations in formation water temperature and salinity is undoubtedly occurring around some salt structures. However, the driving mechanisms for the broad, diffusive upward solute transport in the northern Gulf rim of Arkansas and northern Louisiana are not known. In the Lower Cretaceous of Texas, fluid flow was much more highly focused, and perhaps episodic. It is clear that many areas of the Gulf basin are hydrologically connected and that large‐scale fluid flow, solute transport, and dispersion have occurred. The Na‐Mg‐Ca‐Cl compositions of brines in the areas of the Gulf Coast sedimentary basin reviewed in this article are products of diagenesis and do not reflect the composition of the evaporated marine waters present at the time of sediment deposition. Large differences in Na, Ca, and Mg trends for waters hosted by Mesozoic versus Cenozoic sediments may reflect differences in: (i) the sources of salinity (evaporated seawater for some of the Mesozoic sediments, dissolution of salt for some of the Cenozoic sediments); (ii) sediment lithology (dominantly carbonates for much of the Mesozoic sediments, and dominantly siliciclastics for the Cenozoic sediments); or (iii) residence times of brines associated with these sediments (tens of millions of years versus perhaps days).     

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.     

Long‐term chemical evolution and modification of continental basement brines – a field study from the Schwarzwald,SW Germany

Highly saline, deep‐seated basement brines are of major importance for ore‐forming processes, but their genesis is controversial. Based on studies of fluid inclusions from hydrothermal veins of various ages, we reconstruct the temporal evolution of continental basement fluids from the Variscan Schwarzwald (Germany). During the Carboniferous (vein type i), quartz–tourmaline veins precipitated from low‐salinity (<4.5wt% NaCl + CaCl₂), high‐temperature (≤390°C) H₂O‐NaCl‐(CO₂‐CH₄) fluids with Cl/Br mass ratios = 50–146. In the Permian (vein type ii), cooling of H₂O‐NaCl‐(KCl‐CaCl₂) metamorphic fluids (T ≤ 310°C, 2–4.5wt% NaCl + CaCl₂, Cl/Br mass ratios = 90) leads to the precipitation of quartz‐Sb‐Au veins. Around the Triassic–Jurassic boundary (vein type iii), quartz–haematite veins formed from two distinct fluids: a low‐salinity fluid (similar to (ii)) and a high‐salinity fluid (T = 100–320°C, >20wt% NaCl + CaCl₂, Cl/Br mass ratios = 60–110). Both fluids types were present during vein formation but did not mix with each other (because of hydrogeological reasons). Jurassic–Cretaceous veins (vein type iv) record fluid mixing between an older bittern brine (Cl/Br mass ratios ~80) and a younger halite dissolution brine (Cl/Br mass ratios >1000) of similar salinity, resulting in a mixed H₂O‐NaCl‐CaCl₂ brine (50–140°C, 23–26wt% NaCl + CaCl₂, Cl/Br mass ratios = 80–520). During post‐Cretaceous times (vein type v), the opening of the Upper Rhine Graben and the concomitant juxtaposition of various aquifers, which enabled mixing of high‐ and low‐salinity fluids and resulted in vein formation (multicomponent fluid H₂O‐NaCl‐CaCl₂‐(SO₄‐HCO₃), 70–190°C, 5–25wt% NaCl‐CaCl₂ and Cl/Br mass ratios = 2–140). The first occurrence of highly saline brines is recorded in veins that formed shortly after deposition of halite in the Muschelkalk Ocean above the basement, suggesting an external source of the brine's salinity. Hence, today's brines in the European basement probably developed from inherited evaporitic bittern brines. These were afterwards extensively modified by fluid–rock interaction on their migration paths through the crystalline basement and later by mixing with younger meteoric fluids and halite dissolution brines.     

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.     

Factors controlling free thermal convection in faults in sedimentary basins: implications for the formation of zinc–lead mineral deposits

Stratiform sediment‐hosted Zn–Pb–Ag mineral deposits constitute about 40% of the Earth's zinc resources ( Allen 2001 ), and in most cases their genesis involves the discharge of basinal brines near or on the seafloor through syndepositional faults ( Sangster 2002 ). From the point of view of base metal exploration, it is therefore essential to identify all possible faults that formerly carried the upwelling ore‐forming solutions during mineralising events. This paper presents a numerical investigation of the relative importance of various physical parameters in controlling fluid discharge, recharge and heat transport in faults. A two‐dimensional, free convection of pure water, hydrogeological model is developed for the McArthur basin in northern Australia based on the surface geology, known stratigraphic and structural relationships and regional geophysical interpretations. Numerical experiments and sensitivity analyses reveal that faults with strong initial heat input, due to depth of penetration or magmatic activity, are the most likely candidates to carry discharge fluids to the sites of metal precipitation. Deeper, wider and more permeable faults are more likely to behave as the fluid discharge pathways, whereas shallow, narrow or less permeable faults act as marine water recharge pathways. Compared with these fault‐related factors, aquifer physical properties are less important in determining fluid flow patterns and the geothermal regime. These results are an important step in understanding hydrothermal fluid flow in sedimentary basins in order to develop effective exploration criteria for the location of stratiform Zn–Pb–Ag deposits.     

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.     

Numerical modeling of the origin of calcite mineralization in the Refugio‐Carneros fault,Santa Barbara Basin,California

Many faults in active and exhumed hydrocarbon‐generating basins are characterized by thick deposits of carbonate fault cement of limited vertical and horizontal extent. Based on fluid inclusion and stable isotope characteristics, these deposits have been attributed to upward flow of formation water and hydrocarbons. The present study sought to test this hypothesis by using numerical reactive transport modeling to investigate the origin of calcite cements in the Refugio‐Carneros fault located on the northern flank of the Santa Barbara Basin of southern California. Previous research has shown this calcite to have low δ¹³C values of about ?40 to ?30‰PDB, suggesting that methane‐rich fluids ascended the fault and contributed carbon for the mineralization. Fluid inclusion homogenization temperatures of 80–125°C in the calcite indicate that the fluids also transported significant quantities of heat. Fluid inclusion salinities ranging from fresh water to seawater values and the proximity of the Refugio‐Carneros fault to a zone of groundwater recharge in the Santa Ynez Mountains suggest that calcite precipitation in the fault may have been induced by the oxidation of methane‐rich basinal fluids by infiltrating meteoric fluids descending steeply dipping sedimentary layers on the northern basin flank. This oxidation could have occurred via at least two different mixing scenarios. In the first, overpressures in the central part of the basin may have driven methane‐rich formation waters derived from the Monterey Formation northward toward the basin flanks where they mixed with meteoric water descending from the Santa Ynez Mountains and diverted upward through the Refugio‐Carneros fault. In the second scenario, methane‐rich fluids sourced from deeper Paleogene sediments would have been driven upward by overpressures generated in the fault zones because of deformation, pressure solution, and flow, and released during fault rupture, ultimately mixing with meteoric water at shallow depth. The models in the present study were designed to test this second scenario, and show that in order for the observed fluid inclusion temperatures to be reached within 200 m of the surface, moderate overpressures and high permeabilities were required in the fault zone. Sudden release of overpressure may have been triggered by earthquakes and led to transient pulses of accelerated fluid flow and heat transport along faults, most likely on the order of tens to hundreds of years in duration. While the models also showed that methane‐rich fluids ascending the Refugio‐Carneros fault could be oxidized by meteoric water traversing the Vaqueros Sandstone to form calcite, they raised doubts about whether the length of time and the number of fault pulses needed for mineralization by the fault overpressuring mechanism were too high given existing geologic constraints.     

Morphological evolution of three‐dimensional chemical dissolution front in fluid‐saturated porous media: a numerical simulation approach

This paper is concerned with the morphological evolution of three‐dimensional chemical dissolution fronts that occur in fluid‐saturated porous media. A fully coupled system between porosity, pore‐fluid flow and reactive chemical species transport is considered to describe this phenomenon. Using the newly presented concept of the generalized dimensionless pore fluid pressure‐gradient, which can be used to represent the interaction between solute advection, solute diffusion, chemical kinetics and the shape factor of the soluble mineral, a theoretical criterion has been established to assess the likelihood of instability at a chemical dissolution front in the reactive transport system. To simulate the chemical dissolution front evolution in a three‐dimensional fluid‐saturated porous medium, a numerical procedure combining both the finite difference method and the finite element method has been proposed. As the problem belongs to a complex system science problem, a small randomly generated perturbation of porosity is added to the initial porosity of a three‐dimensional homogeneous domain to trigger instability of a planar chemical dissolution front during its propagation within the fluid‐saturated porous medium. To test the correctness and accuracy of the proposed numerical procedure, a three‐dimensional benchmark problem has been constructed and the related analytical solution has been derived. This enables using the proposed numerical procedure for simulating the morphological evolution of a three‐dimensional chemical dissolution front from a stable, planar state into an unstable, fingering state. The related numerical results demonstrate that the proposed numerical procedure is useful for, and capable of, simulating the morphological instability of a three‐dimensional chemical dissolution front within a fluid‐saturated porous medium.     

Modeling thermal convection in supradetachment basins: example from western Norway

The juxtaposition of fault‐bounded sedimentary basins, above crustal‐scale detachments, with warmer exhumed footwalls can lead to thermal convection of the fluids in the sediments. The Devonian basins of western Norway are examples of supradetachment basins that formed in the hanging wall of the Nordfjord‐Sogn Detachment Zone. In the central part of the Hornelen and Kvamshesten basins, the basin‐fill is chiefly represented by fluvial sandstones and minor lacustrine siltstones, whereas the fault margins are dominated by fanglomerates along the detachment contact. Prominent alteration and low‐greenschist facies metamorphic conditions are associated with the peak temperature estimates of the sediments close to the detachment shear zone. Fluid circulation may have been active during the burial of the sediments, and we quantify the potential role played by thermal convection in redistributing heat within the basins. Different models are tested with homogeneous and layered basin‐fill and with material transport properties corresponding to sandstones and siltstones. We found that thermally driven fluid flow is expected in supradetachment basins as a transient process during the exhumation of warmer footwalls. We demonstrate that the fluid flow may have significantly affected the temperature distribution in the upper five kilometers of the Devonian basins of western Norway. The temperature anomaly induced by the flow may locally reach about 80°C. The sedimentary layering formed by sand‐ and siltstones strata does not inhibit fluid circulation at the scale of the basin. The presence of fluid pathways along the detachment has an important impact on the flow and allows an efficient drainage of the basin by channelizing fluids upward along the detachment.     

Hydrothermal dolomitization,mixing corrosion and deep burial porosity formation: numerical results from 1‐D reactive transport models

Major corrosion has been found at depth in carbonate hydrocarbon reservoirs from different geologic provinces. Fluid inclusion microthermometry and stable isotopic compositions of carbonate cements, predating major corrosion, constrain the interpretation of the evolution of parental fluids during progressive burial and prior to the major corrosion event. Post‐major corrosion mineral paragenesis includes pyrite (‐marcasite), anhydrite, kaolinite (dickite) and fluorite. Although the post‐corrosion mineral paragenesis represents minor volumes of rock, it may provide valuable insights into the post‐corrosion brine chemistry. Using reactive transport numerical models, the roles of cooling and/or mixing of brines on corrosion have been evaluated as controls for dolomitization, deep burial corrosion and precipitation of the post‐corrosion mineral paragenesis. Modelling results show that cooling of deep‐seated fluids moving upward along a fracture may cause minor calcite dissolution and porosity generation. Significant dolomitization along a fracture zone and nearby host‐rock only occurs when deep‐seated fluids have high salinities (4 mol Cl kg^?1 of solution) and Darcian flow rates are relatively high (1 m³ m^?2 year^?1). Only minor volumes of quartz and fluorite precipitate in the newly formed porosity. Moreover, modelling results cannot reproduce the authigenic precipitation of kaolinite (dickite at high temperatures) by cooling. As an alternative to cooling as a cause of corrosion, mixing between two brines of different compositions and salinities is represented by two main cases. One case consists of the flow up along a fracture of deep‐seated fluids with higher salinities than the fluid in the wall rock. Dolomite does not precipitate at a fracture zone. Nevertheless, minor volumes of dolomite are formed away from the fracture. The post‐corrosion mineral paragenesis can be partly reproduced, and the results are comparable to those obtained from cooling calculations. Minor volumes of quartz and fluorite are formed, and kaolinite‐dickite does not precipitate. The major outputs of this scenario are calcite dissolution and slight net increase in porosity. A second case corresponds to the mixing of low salinity deep‐seated fluids, flowing up along fractures, with high salinity brines within the wall rock. Calculations predict major dissolution of calcite and precipitation of dolomite. The post‐corrosion mineral paragenesis can be reproduced. High volumes of quartz, fluorite and kaolinite‐dickite precipitate and may even completely occlude newly formed porosity.     

Structural controls on hydrothermal flow in a segmented rift system, Taupo Volcanic Zone, New Zealand

Vigorous hydrothermal convection transfers 10 times the average continental heat flow through the central Taupo Volcanic Zone (TVZ), a region of active extension (approximately 8 mm year^?1) and productive rhyolitic volcanism. Over 20 high‐temperature (>250°C) geothermal fields occur within Quaternary pyroclastic basins, with convective circulation to depths of 7–8 km presumably extending through basement rocks. Parallel‐striking normal faults, fractures and dikes dissect the convective regime, interacting with fluids to either enhance or restrict flow according to the relative permeability of structure and host rock. In the basement, high bulk permeability is maintained by focussed flow through faults and associated fractures well oriented for reactivation in the prevailing stress field. In contrast, distributed flow through fault‐bounded compartments prevails within Quaternary basins, masking any signal of deeper structural control. Exceptions occur where more competent rocks are exposed at the surface. As in narrow magmatic rifts elsewhere, the extensional fabric is partitioned into discrete rift segments linked along strike by accommodation zones. Eighty per cent of TVZ geothermal fields correlate spatially with rift architecture, with 60% located in accommodation zones. We suggest that segmented rift fabrics generate bulk permeability anisotropy that is to some extent predictable, with rift segments characterized by enhanced axial flow, and accommodation zones characterized by locally enhanced vertical permeability that is tectonically maintained. This provides a plausible explanation for the common occurrence of geothermal fields within accommodation zones and their notable absence within densely faulted rift segments. Maintenance of structural permeability in zones of active hydrothermal precipitation necessarily requires repeated brittle failure. Geothermal plumes therefore exploit tectonically maintained permeability within accommodation zones, with rift segments functioning mostly as drawdown regions. The influence of rift architecture on flow paths has important implications for geothermal extraction and epithermal mineral exploration within the TVZ and other structurally segmented hydrothermal systems, both active and extinct.