Groundwater age,brine migration,and large‐scale solute transport in the Alberta Basin,Canada

There is a great contrast in geochemical and hydrogeologic estimates of the residence times of pore fluids in sedimentary basins. This contrast is particularly evident in the Alberta Basin, Canada, which has served as the study area for important studies of long‐term fluid flow and transport. To address these differences, we developed two‐dimensional simulations of groundwater age, constrained by both hydrogeologic and geochemical observations, to estimate the residence time of fluids and the amount and timing of flushing by meteoric waters in the Alberta Basin. Results suggest that old, residual brines have been retained in the deepest parts of the basin since their formation ca. 400 Ma, but significant dilution by younger waters has reduced the age of these pore waters to no more than approximately 200 My. Shallower formations have been flushed extensively by fresh, young waters. Loss of brines and dilution of older pore waters occurred primarily after the uplift of the Rockies with the introduction of the gravity‐driven flow regime. Despite these large changes in flow regime, solute exchange between deep saline aquifers and the overlying vigorous freshwater flow system was found to be consistent with long‐term dispersive mixing across subhorizontal concentration gradients rather than by direct flushing. Sensitivity studies using an analytic solution supported the use of 100 m for transverse dispersivity in large‐scale numerical models. These simulations confirm that the age and origin of brines are in many cases poor indicators of long‐term solute transport rates in sedimentary basins, but the geochemical indicators that are used to determine the origin of brines can provide useful constraints for calculating groundwater age and are far more commonly available than isotopic groundwater age tracers.     

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.     

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.     

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.     

Palaeo‐formation water evolution in the Latrobe aquifer,Gippsland Basin,south‐eastern Australia continental shelf

Offshore fresh or brackish groundwater has been observed around the globe and represents an interesting but unusual freshwater reserve. Formation waters in sedimentary basins evolve at geological time through fluid–rock interactions and water movements in aquifers. However, the mechanism and timing of freshwater displacing and mixing with pre‐existing formation water offshore under the seafloor has not been investigated in many cases. The growing need for developing freshwater resources in deeper parts of sedimentary basins that have not been economic or technically feasible in the past, may potentially lead to an increasing conflict with petroleum production or injection of carbon dioxide. For being able to assess and mitigate possible impacts of fluid production or injection on groundwater flow and quality, a better understanding of the natural history of the interaction between fresh meteoric water and deep basin formation water is necessary. A low‐salinity wedge of meteoric origin with less than 5000 ppm currently extends to about 20 km offshore in the confined Latrobe aquifer in the Gippsland Basin (Australia). The Latrobe aquifer is a freshwater resource in the onshore, hosts major petroleum reservoirs and has been considered for carbon dioxide storage in the offshore parts of the basin. The objective of this study is to constrain the evolution of formation water in the Latrobe aquifer by investigating the water naturally trapped in fluid inclusions during burial. The measured palaeo‐salinities from onshore and offshore rock samples have a minimum of about 12 500 ppm (NaCl equivalent) and a maximum of about 50 000 ppm. Most of the salinities are in the 32 000–35 000 ppm range. There is no evidence for freshwater in fluid inclusions and the variation in palaeo‐salinity across the basin is consistent with the palaeogeography of deposition of the sedimentary rocks. The current low‐salinity water wedge must have started to form recently after most of the diagenetic processes that led to the trapping of water in fluid inclusions happened. The minimum homogenisation temperatures (T_h) recorded are consistent with current formation temperature. However, they are generally higher than present day suggesting that hotter temperatures were attained in the past. The T_h and salinity data together suggest that the fluid inclusions record the diagenetic modification of connate water to higher salinities over a time period that was accompanied by an increase in temperature, consistent with a westward palaeo‐fluid flow from the deeper part of the basin through the aquifer. Subsequent pore‐water evolution from palaeo‐ to current day conditions is consistent with an influx of fresher and cooler meteoric water into the Latrobe Group. The meteoric recharge originates from the area of the Baragwanath anticline in the onshore part of the basin where the Latrobe Group subcrops at high elevations.     

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).     

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.     

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.     

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.     

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.     

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

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

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

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

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.     

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.     

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.     

Numerical models of petroleum migration via buoyancy-driven porosity waves in viscously deformable sediments

The apparent ability of petroleum to migrate rapidly through low permeability sediments in sedimentary basins has led to many questions about the manner of its transport. One possible explanation is suggested from observations of the compaction of viscously deformable porous media. These systems have been found in some cases to give rise to regions of locally elevated liquid fraction, in the form of fluid‐filled porosity waves that can ascend at rates much greater than that of the background flow. Previous research on the phenomenon has focused on its implications for magma transport, but recognition of the fact that the compaction of viscous porous media can take place in sedimentary basins has suggested the possibility that porosity waves could similarly be important for hydrocarbon transport. The purpose of the present study was to test this hypothesis by quantifying the transport that would occur as the result of porosity waves initiated during the conversion of kerogen to petroleum. A one‐dimensional numerical model was constructed solving equations for the mechanics of viscous compaction and for the kinetics of reactions describing the formation of petroleum from kerogen. The results showed that porosity waves would develop readily in viscously deformable regions of sedimentary basins, but would not necessarily provide enhanced transport over that of the background flow regime. In order for the waves to achieve this enhanced transport, they must develop high amplitudes, i.e. high porosities relative to the background porous medium. To achieve the high wave amplitudes, the background porosity must be very low in absolute terms. In addition, high kerogen contents are needed in the source layer, and the source layer needs to be buried rapidly to a high temperature region of the oil window. Considerable uncertainty exists as to the value of the matrix shear viscosity of sediments in basins. However, the wave volumetric transport capacity was not found to be significantly altered as a result of variations in the value of this parameter. The physical form of the waves was strongly altered by the matrix shear viscosity, with higher values leading to lower amplitudes and generation frequency, but higher wavelengths. Thus the waves become less recognizable physically at higher values of the matrix shear viscosity. As the waves ascend to higher stratigraphic levels, where the porosity is higher, they gradually lose their physical definition and become absorbed into the background.     

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

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

Boron isotope geochemistry of Na‐bicarbonate,Na‐chloride,and Ca‐chloride waters from the Northern Apennine Foredeep basin: other pieces of the sedimentary basin puzzle

The boron stable isotope ratio δ¹¹B of 12 water samples representative of three chemical facies (fresh Na‐bicarbonate, brackish Na‐chloride, saline, and brine Ca‐chloride) has been analyzed. Interpretation of the δ¹¹B data, along with the chemical compositions, reveals that Na‐carbonate waters from the Northern Apennine are of meteoric origin, with boron contributions from clay desorption and mixing with seawater‐derived fluids of Na‐chloride or Ca‐chloride compositions. The comparison of our new results with the literature data on other sedimentary basins of Mediterranean, and worldwide, confirms the contribution of Na‐bicarbonate waters to the genesis of mud volcano fluids. The Na‐chloride sample of Salvarola (SAL), which may represent the end‐member of the mud volcanoes, and the Ca‐chloride brine water from Salsomaggiore (SM) indicate boron release from clays compatible with the diagenetic process. The empirical equation: relating boron concentration and the stable isotope composition of the fluids observed in this study and the literature is proposed to trace the effect of diagenesis in sedimentary basins. A geothermometer associated to the diagenetic equation is also proposed: The application of this equation to obtain reservoir temperatures from δ¹¹B compositions of waters should be carefully evaluated against the results obtained from other chemical and isotopic geothermometers from other basins around the world.     

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.     

A hydrogeologic model of stratiform copper mineralization in the Midcontinent Rift System, Northern Michigan, USA

This paper presents a suite of two‐dimensional mathematical models of basin‐scale groundwater flow and heat transfer for the middle Proterozoic Midcontinent Rift System. The models were used to assess the hydrodynamic driving mechanisms responsible for main‐stage stratiform copper mineralization of the basal Nonesuch Formation during the post‐volcanic/pre‐compressional phase of basin evolution. Results suggest that compaction of the basal aquifer (Copper Harbor Formation), in response to mechanical loading during deposition of the overlying Freda Sandstone, generated a pulse of marginward‐directed, compaction‐driven discharge of cupriferous brines from within the basal aquifer. The timing of this pulse is consistent with the radiometric dates for the timing of mineralization. Thinning of the basal aquifer near White Pine, Michigan, enhanced stratiform copper mineralization. Focused upward leakage of copper‐laden brines into the lowermost facies of the pyrite‐rich Nonesuch Formation resulted in copper sulfide mineralization in response to a change in oxidation state. Economic‐grade mineralization within the White Pine ore district is a consequence of intense focusing of compaction‐driven discharge, and corresponding amplification of leakage into the basal Nonesuch Formation, where the basal aquifer thins dramatically atop the Porcupine Mountains volcanic structure. Equilibrium geochemical modeling and mass‐balance calculations support this conclusion. We also assessed whether topography and density‐driven flow systems could have caused ore genesis at White Pine. Topography‐driven flow associated with the Ottawan orogeny was discounted because it post‐dates main‐stage ore genesis and because recent seismic interpretations of basin inversion indicates that basin geometry would not be conductive to ore genesis. Density‐driven flow systems did not produce focused discharge in the vicinity of the White Pine ore district.