Impacts of Pleistocene glacial loading on abnormal pore‐water pressure in the eastern Michigan Basin

The hydromechanical effects of Pleistocene glacial loading on the Michigan Basin are assessed using numerical analysis based on coupled stress and pore‐water pressure. The two‐dimensional model domain included the Basin cross section and extended 10 km into the Precambrian. In the analysis, we considered the effects of the number of glacial loading cycles, the presence and connectedness of a deep Cambrian aquifer, the direction of glacial advance, the effect of a wet versus dry glacier/soil interface, topographic effects, density‐driven flow effects, and lithosphere flexure on the development of anomalous pressures. The modeling results were compared with data collected from highly instrumented wells completed in the eastern margin of the Basin. The present‐day results define regions of significant underpressure in the upper Ordovician and lower Silurian formations characterized by very low hydraulic conductivity and regions of overpressure where hydraulic conductivity is higher. To achieve the degree of underpressure observed in the instrumented wells using the model, a specific loading cycle applied over 100 000 years was repeated four times. As the number of loading cycles increased, the Paleozoic formations reached a state where the underpressures remain constant. Our results also illustrate the difference in poroelastic modeling between the application of mechanical loads on the land surface and the application of an equivalent hydraulic head, where the latter developed overpressures rather than the observed underpressures. The modeling also shows that the overpressures observed in the Cambrian formations are most likely to be the result of density‐driven flow. Finally, the simulations show that the effects of lithosphere flexure in the hydromechanical models results in the development of lateral stresses that generate overpressures rather than underpressures in the southern half of the domain. As there are no suitable observation points, these results remain unconfirmed, and further study is warranted.     

Changes in fluid pathways in a calcite vein mesh (Natih Formation,Oman Mountains): insights from stable isotopes

We present a structural, microstructural, and stable isotope study of a calcite vein mesh within the Cretaceous Natih Formation in the Oman Mountains to explore changes in fluid pathways during vein formation. Stage 1 veins form a mesh of steeply dipping crack‐seal extension veins confined to a 3.5‐m‐thick stratigraphic interval. Different strike orientations of Stage 1 veins show mutually crosscutting relationships. Stage 2 veins occur in the dilatant parts of a younger normal fault interpreted to penetrate the stratigraphy below. The δ¹⁸O composition of the host rock ranges from 21.8‰ to 23.7‰. The δ¹³C composition ranges from 1.5‰ to 2.3‰. This range is consistent with regionally developed diagenetic alteration at top of the Natih Formation. The δ¹⁸O composition of vein calcite varies from 22.5‰ to 26.2‰, whereas δ¹³C composition ranges from ?0.8‰ to 2.1‰. A first trend observed in Stage 1 veins involves a decrease of δ¹³C to compositions nearly 1.3‰ lower than the host rock, whereas δ¹⁸O remains constant. A second trend observed in Stage 2 calcite has δ¹⁸O values up to 3.3‰ higher than the host rock, whereas the δ¹³C composition is similar. Stable isotope data and microstructures indicate an episodic flow regime for both stages. During Stage 1, formation of a stratabound vein mesh involved bedding‐parallel flow, under near‐lithostatic fluid pressures. The ¹⁸O fluid composition was host rock‐buffered, whereas ¹³C composition was relatively depleted. This may reflect reaction of low ¹³C CO₂ derived by fluid interaction with organic matter in the limestones. Stage 2 vein formation is associated with fault‐controlled fluid flow accessing fluids in equilibrium with limestones about 50 m beneath. We highlight how evolution of effective stress states and the growth of faults influence the hydraulic connectivity in fracture networks and we demonstrate the value of stable isotopes in tracking changes in fluid pathways.     

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.