Molybdenum speciation tracking hydrocarbon migration in fine-grained sedimentary rocks
Explore the source record for details and available documents.
Ardakania, Omid H.↗
SEARCH · Search NASA
Results for “SEDIMENTARY ROCK”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Mineralogical, magnetic and geochemical data constrain the pathways and extent of weathering of mineralized sedimentary rocks
The oxidative weathering of sulfidic rock can profoundly impact watersheds through the resulting export of acidity and metals. Weathering leaves a record of mineral transformation, particularly involving minor redox-sensitive phases, that can inform the development of conceptual and quantitative models. In sulfidic sedimentary rocks, however, variations in depositional history, diagenesis and mineralization can change or overprint the distributions of these trace minerals, complicating the interpretation of weathering signatures. Here we show that a combination of bulk mineralogical and geochemical techniques, micrometer-resolution X-ray fluorescence microprobe analysis and rock magnetic measurements, applied to drill core samples and single weathered fractures, can provide data that enable the development of a geochemically consistent weathering model. This work focused on one watershed in the Upper Colorado River Basin sitting within the Mesaverde Formation, a sedimentary sandstone bedrock with disseminated sulfide minerals, including pyrite and sphalerite, that were introduced during diagenesis and subsequent magmatic-hydrothermal mineralization. Combined analytical methods revealed the pathways of iron (Fe), carbonate and silicate mineral weathering and showed how pH controls element retention or release from the actively weathering fractured sandstone. Drill core logging, whole rock X-ray diffraction, and geochemical measurements document the progression from unweathered rock at depth to weathered rock at the surface. X-ray microprobe analyses of a 1-cm size weathering profile along a fracture surface are consistent with the mobilization of Fe(II) and Fe(III) into acidic pore water from the dissolution of primary pyrite, Fe-sphalerite, chlorite, and minor siderite and pyrrhotite. These reactions are followed by the precipitation of secondary minerals such as of goethite and jarosite, a Fe-(oxyhydr)oxide and hydrous Fe(III) sulfate, respectively. Microscale analyses also helped explain the weathering reactions responsible for the mineralogical transformations observed in the top and most weathered section of the drill core. For example, dissolution of feldspar and chlorite neutralizes the acidity generated by Fe and sulfide mineral oxidation, oversaturating the solution in both Fe-oxides. The combination of X-ray spectromicroscopy and magnetic measurements show that the Fe(III) product is goethite, mainly present either as a coating on fracture surfaces in the actively weathering region of the core or more homogeneously contained within the unconsolidated regolith at the top of the core. Low-temperature magnetic data reveal the presence of ferromagnetic Fe-sulfide pyrrhotite that, although it occurs at trace concentrations, could provide a qualitative proxy for unweathered sulfide minerals because the loss of pyrrhotite is associated with the onset of oxidative weathering. Pyrrhotite loss and goethite formation are detectable through room-temperature magnetic coercivity changes, suggesting that rock magnetic measurements can determine weathering intensity in rock samples at many scales. In conclusion, this work contributes evidence that the weathering of sulfidic sedimentary rocks follows a geochemical pattern in which the abundance of sulfide minerals controls the generation of acidity and dissolved elements, and the pH-dependent mobility of these elements controls their export to the ground- and surface-water.
58 GEOSCIENCES↗
Core to Local Scale Assessment of REE and Critical Minerals in Sedimentary Rocks from the Powder River Basin
For presentation at the Geological Society of America 2023 Annual Meeting, Pittsburgh, PA, October 15-18, 2023
Montross, Scott↗
A visco-plastic framework for interface processes in sedimentary reservoir rocks at HPHT conditions
As energy operations face the challenge of reservoirs at ever-increasing depths, modelling the response of reservoir rocks at high pressure and high temperature (HPHT) conditions is a crucial step for successfully unlocking new resources. At these conditions, rocks can experience thermal and pressure sensitivity, as well as admit internal phase changes at the pore-solid interfaces that determine their macroscopic response to external loading. It is therefore expected that constitutive laws of Geomechanics should be enriched to accommodate such effects. In this work, a multi-physics constitutive theory for sedimentary rocks is proposed within the general framework of viscoplasticity. The viscosity of the material is assumed to be a function of the temperature, pore-pressure and energy required to alter the inter-granular interfaces. This energy is expressed through the chemical potentials of the phases involved during processes like chemical dissolution/precipitation of the interfaces or mechanical debonding. The resulting flow law and corresponding stress equilibrium are coupled to the energy and mass conservation laws, constituting a closed system of equations, which is solved using the Finite Element simulator REDBACK. A series of numerical tests for performance and calibration is then successfully performed against different types of reservoir rocks, confining pressures and temperatures (from room temperature to over 800K). Here, we show that the mechanical response of sedimentary porous rocks at strains usually achieved in laboratory testing can be explained by accounting for the interface processes taking place at the cementitious material bonding the grains, a process that can also determine the brittle-to-ductile transition of sedimentary rocks.
58 GEOSCIENCES↗
Global patterns of radiocarbon depletion in subsoil linked to rock-derived organic carbon
Organic matter stored in sedimentary rocks is one of the largest stocks of carbon at Earth’s surface. The fate of this rock organic carbon (OC petro ) during weathering in soils influences the geological carbon cycle, and impacts soil radiocarbon content that is used to quantify soil carbon turnover. Here, we assess the potential contribution of OC petro to soils, using a mixing model generated by a global dataset of soil radiocarbon measurements ( 14 C). Soils developed on sedimentary rocks (rather than on igneous substrate) have a paired OC content and 14 C values consistent with OC petro input, giving rise to apparent increase in soil residence time. We call for renewed assessment of OC petro input to soils, in terms of its impact on soil radiocarbon inventories, and its potential to release carbon dioxide.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Pore-Scale Study on the Positive Feedback Between Stress and Porosity Caused by Pressure Solution in Porous Media
Pressure solution is an important process in the evolution of sedimentary rocks, which provide storage space for most of our petroleum resources. It directly influences the generation, migration, and storage of petroleum fluids in subsurface sedimentary rocks. Here, in this paper, we develop a pore-scale, mechanochemical model to demonstrate a possible positive feedback between the local porosity and pore surface stress, in which a higher local porosity causes a higher local pore surface stress, thus enhancing pressure solution and consequently further increasing the local porosity. Pore surface stress represents stress on a solid grain adjacent to a pore. Specifically, the pore-scale, mechanochemical model directly simulates the stress distribution over solid and pore surfaces using a finite element model. The dissolution of solids at the solid-pore interfaces under a far-from-equilibrium condition is simulated using a first-order kinetics model that accounts for the local stress distribution. The updated pore geometry, caused by pore surface dissolution, is then used in the stress simulation in the next numerical iteration. Two types of porous media, the Oriskany sandstone and an artificial porous medium with spherical pores, were tested in the mechanochemical simulation. The positive stress-porosity feedback during pressure solution was observed in both samples. In addition, the model quantitatively illustrated the distribution of local mineral dissolution rates on all pore surfaces, as well as its relation to the effective mineral dissolution rate of the entire sample. Based on the comparison between the two porous media, the local mineral dissolution was regulated by pore space distribution, geometry, and coalescence during pressure solution. This work is the first that uses direct, pore-scale numerical simulation to demonstrate the positive stress-porosity feedback during pressure solution, which has the potential to advance the understanding of the mechanical-chemical (MC) coupling in many geological processes that are relevant to subsurface energy systems, such as the recovery of petroleum hydrocarbons and geothermal energy.
CT scanning↗
High-frequency and Vertical Variations of Soil Moisture and their Impact on Ecosystem Dynamics
The vertical distribution of subsurface moisture is a major control on ecosystem dynamics and climate. The Mediterranean climate of California, with wet winters and dry summers, provides a testbed for understanding how transpiration is sustained through extended dry periods. The project developed and tested a new subsurface hydrology model that includes a parameterization of stochastic hydraulic conductivity to capture fracture flow, a diagnostic model of rooting profile that depends on the subsurface hydrology, and a representation of hydraulic redistribution by plant roots. A sensitivity study shows that it is plants with rooting depths >8 m that could access the subsurface moisture. The insight is applied to understand differential tree mortality during the 2012-2016 drought in California. The analysis shows that greater mortality trees on plutonic rocks than on sedimentary rocks. The study demonstrates, for the first time, the importance of lithology for the distribution and health of ecosystems.
54 ENVIRONMENTAL SCIENCES↗
Concentration‐Discharge Relationships of Dissolved Rhenium in Alpine Catchments Reveal Its Use as a Tracer of Oxidative Weathering
Abstract Oxidative weathering of sedimentary rocks plays an important role in the global carbon cycle. Rhenium (Re) has been proposed as a tracer of rock organic carbon (OC petro ) oxidation. However, the sources of Re and its mobilization by hydrological processes remain poorly constrained. Here, we examine dissolved Re as a function of water discharge, using samples collected from three alpine catchments that drain sedimentary rocks in Switzerland (Erlenbach and Vogelbach) and Colorado, USA (East River). The Swiss catchments reveal a higher dissolved Re flux in the catchment with higher erosion rates, but have similar [Re]/[Na + ] and [Re]/[SO 4 2− ] ratios, which indicate a dominance of Re from OC petro . Despite differences in rock type and hydro‐climatic setting, the three catchments have a positive correlation between river water [Re]/[Na + ] and [Re]/[SO 4 2− ] and water discharge. We propose that this reflects preferential routing of Re from a near‐surface, oxidative weathering zone. The observations support the use of Re as a proxy to trace rock‐organic carbon oxidation, and suggest it may be a hydrological tracer of vadose zone processes. We apply the Re proxy and estimate CO 2 release by OC petro oxidation of 5.7 +6.6 / −2.0 tC km −2 yr −1 for the Erlenbach. The overall weathering intensity was ∼40%, meaning that the corresponding export of unweathered OC petro in river sediments is large, and the findings call for more measurements of OC petro oxidation in mountains and rivers as they cross floodplains.
54 ENVIRONMENTAL SCIENCES↗
Validating the rhenium proxy for rock organic carbon oxidation using weathering profiles
Chemical weathering over geological timescales acts as a source or sink of atmospheric carbon dioxide (CO 2 ), while influencing long-term redox cycling and atmospheric oxygen (O 2 ) at Earth's surface. There is a growing recognition that the oxidative weathering of rock organic carbon (OC petro ) can release more CO 2 than is locally drawn down by silicate weathering, and may vary due to changes in erosion and climate. The element rhenium (Re) has emerged as a proxy to track the oxidative weathering of OC petro , yet uncertainties in its application remain namely that we lack a systematic assessment of the comparative mobility of Re and OC petro during sedimentary rock weathering. Here we measure Re and OC petro loss across gradients in rock weathering at 9 global sites, spanning a range of initial OC petro values from ~0.2 % to >10 %. We use titanium to account for volume changes during weathering and assess Re and OC petro loss alongside major elements that reflect silicate (Na, Mg), carbonate (Ca, Mg) and sulfide (S) weathering. Across the dataset, Re loss is correlated with OC petro loss but not with loss of any other major element. Further, across the weathering profiles, the average molar ratio of OC petro to Re loss was 0.84 ± 0.15, with 8 out of 9 sites having a ratio >0.74. At one site (Marcellus Shale), the average ratio was lower at 0.58 ± 0.11. The excess loss of Re matches expectations that, typically, between ~0 and 20 % of the Re liberated by sedimentary rock weathering derives from silicate or sulfide phases, while some OC petro may be physically or chemically protected from weathering. Overall, our measurements provide validation for the Re proxy of OC petro oxidation and allow future work to further improve our knowledge of regional and global-scale rates of this important source of CO 2 in the geochemical carbon cycle.
58 GEOSCIENCES↗
Permeability and Induced Polarization of Mudstones
Electrical measurements can be used to estimate hydraulic properties such as permeability ( k ) in sedimentary rocks. Previous work has focused on sandstones, siltstones, and carbonates, while investigations on mudstones have rarely been reported. In this study, we report on electrical geophysical measurements for 23 mudstone samples using an experimental approach designed to reliably saturate these low permeability mudstones. The modified Hagen-Poiseuille model linking permeability to the formation factor ( F ) and an effective pore radius ( r ) provides an excellent fit to the data set with a near-constant pore radius, indicating that the effective porosity (1/ F ) is the controlling factor on k . In these samples, the surface area normalized to pore volume ( S por ), frequently used in permeability estimation models, varies by 1–2 orders of magnitude and is thus not a reliable proxy of the inverse effective hydraulic radius. The formation factor also exerts the primary control on induced polarization (IP) parameters, whereas Spor shows no relation to the IP parameters. A strong linear relationship is found between IP parameters (imaginary conductivity and normalized chargeability) and surface conductivity, although the proportionality factor is significantly lower than those observed in more permeable rocks and sediments. Apparent relationships between the polarization strength-derived and time constant-derived geophysical length scales and the effective hydraulic radius appear to be driven by variations in the electrochemical parameters (i.e., specific polarizability and diffusion coefficient). Overall, these findings emphasize that predicting hydraulic properties from electrical measurements in fine-grained rocks remains challenging and requires further investigation into the electrochemical properties involved.
58 GEOSCIENCES↗
Experimental study of compaction localization in carbonate rock and constitutive modeling of mechanical anisotropy
Abstract Sedimentary rocks are inherently anisotropic and prone to strain localization. While the influence of rock anisotropy on the brittle/dilative regime has been studied extensively, its influence on the ductile/compactive regime is much less explored. This paper discusses the anisotropic behavior of a high‐porosity carbonate rock from central Europe (the Maastricht Tuffeau). A set of triaxial tests with concurrent x‐ray tomography has been performed at different confining pressures. The anisotropic characteristics of this rock have been investigated by testing samples cored at different inclinations of the bedding, thus revealing non‐negligible effects of the coring direction on yielding and compaction behavior. Specifically, samples cored perpendicular to bedding display higher strength and longer stages of post‐yielding deformation before manifesting re‐hardening. Despite such alterations of the inelastic response, Digital Image Correlation has revealed that the strain localization mode is independent of the coring direction, thus being primarily affected by the confinement level. To capture the observed interaction between material anisotropy and compaction behavior at the continuum‐scale, an elastoplastic constitutive law has been proposed. For this purpose, a set of tensorial bases has been introduced to replicate how the oriented rock fabric modulates the yielding and plastic flow characteristics of the material. The analyses show that the impact of the coring direction on yield function and plastic flow rule is fundamentally different, thus requiring the use of distinct projection strategies (a strategy here defined heterotopic mapping). The performance of the model, studied through parametric analyses and by calibrating the experimental results, illustrates the improved capability of the proposed constitutive approach when applied to strongly anisotropic porous rocks.
Shahin, Ghassan↗
Characterizing Stress Roughness at Utah FORGE Through Simulation of Hydraulic Fracture Growth During 16A Stage 3 Stimulation
In most geologic formations, the vertical gradient of the minimum horizontal stress (𝜎 h ) exceeds the hydrostatic gradient, often driving hydraulic fractures to propagate upward. However, observations indicate that the upward growth of real hydraulic fractures is less pronounced than theoretical predictions based on smoothly varying 𝜎 h fields. In fact, the layered structure of sedimentary rocks hinders fracture propagation across layers, partially explaining the limited height growth observed in practice (Zoback et al., 2022). While crystalline rocks lack the pervasive layered fabric of sedimentary formations, they still possess structural fabric formed over their geologic history, resulting in inherently "rough" in situ stress fields. Recent studies have identified another key factor: the "roughness" of in situ stress, characterized by temporal, relatively short-wavelength fluctuations superimposed on the overall stress gradient. This stress roughness leads to apparent toughness anisotropy, where the vertical toughness appears significantly larger than the horizontal toughness (P. Fu et al., 2019; Dontsov & Suarez-Rivera, 2021). Neglecting the effects of rock fabric and stress roughness can yield inaccurate predictions of hydraulic fracture geometry and growth rates at larger length scales.
58 GEOSCIENCES↗
Identifying Direct SP -Converted Waves Constrains Local Induced Earthquake Depths
Abstract Seismicity in southern Kansas and northern Oklahoma in the past decade has been associated with fluid injections. In southcentral Kansas, the Wellington earthquake catalog is primarily composed of local, low-magnitude events. Approximately 22% of recorded earthquakes over a 2.5 yr period exhibit a seismic phase arriving between the direct P phase and direct S phase with particle motion similar to the P wave. This intermediate phase was identified as an S to P conversion (SP phase) occurring in the sedimentary rocks instead of the hypothesized basement to sedimentary section transition. We exploit the SP-converted phases to improve the depth accuracy of shallow earthquakes and to constrain VP/VS. The revised depth calculations further confirm that these local induced earthquakes are occurring in the shallow crystalline basement, below the sedimentary section in which fluids are injected.
Geochemistry & Geophysics↗
Large-Scale Molecular Dynamics Simulation of the Dehydration of a Suspension of Smectite Clay Nanoparticles
Fine-grained sediments and sedimentary rocks play important roles in a variety of modern energy technologies from petroleum geology to geological carbon sequestration and radioactive waste management. However, despite their utility and ubiquity, many of their properties remain poorly understood. In particular, the ability to predict the permeability and mechanics of these media remains a persistent fundamental challenge in the geosciences. In the present work, we show how large-scale classical molecular dynamics (MD) simulations can help interpret the properties of fine-grained sedimentary material. All-atom simulations containing 30 discrete clay particles are utilized to understand the evolution of a clay nanoparticle suspension during its progressive dehydration. Microstructural (pore size distribution, tortuosity, anisotropy), thermodynamic (enthalpy and free energy of hydration, anion exclusion), mechanical (total suction), and transport properties (diffusion coefficient tensors of water and sodium) are calculated and compared to the experiment. Overall, our results provide new insight into the coupled chemistry, mechanics, and transport properties of disordered nanoparticle assemblages and shed light upon the important role of water films in controlling these properties.
02 PETROLEUM↗
(U)SAXS characterization of porous microstructure of chert: insights into organic matter preservation
This study characterizes the microstructure and mineralogy of 132 (ODP sample), 1000 and 1880 million-year-old chert samples. By using ultra-small-angle X-ray scattering (USAXS), wide-angle X-ray scattering and other techniques, the preservation of organic matter (OM) in these samples is studied. The scarce microstructural data reported on chert contrast with many studies addressing porosity evolution in other sedimentary rocks. The aim of this work is to solve the distribution of OM and silica in chert by characterizing samples before and after combustion to pinpoint the OM distribution inside the porous silica matrix. The samples are predominantly composed of alpha quartz and show increasing crystallite sizes up to 33 ± 5 nm (1σ standard deviation or SD). In older samples, low water abundances (~0.03%) suggest progressive dehydration. (U)SAXS data reveal a porous matrix that evolves over geological time, including, from younger to older samples, (1) a decreasing pore volume down to 1%, (2) greater pore sizes hosting OM, (3) decreasing specific surface area values from younger (9.3 ± 0.1 m 2 g -1 ) to older samples (0.63 ± 0.07 m 2 g -1 , 1σ SD) and (4) a lower background intensity correlated to decreasing hydrogen abundances. The pore-volume distributions (PVDs) show that pores ranging from 4 to 100 nm accumulate the greater volume fraction of OM. Raman data show aromatic organic clusters up to 20 nm in older samples. Raman and PVD data suggest that OM is located mostly in mesopores. Observed structural changes, silica–OM interactions and the hydrophobicity of the OM could explain the OM preservation in chert.
(U)SAXS↗
Rock-Sourced Nitrogen in Semi-Arid, Shale-Derived California Soils
Models suggest that rock-derived nitrogen (N) inputs are of global importance to ecosystem N budgets; however, field studies demonstrating the significance of rock N inputs are rare. We examined rock-derived N fluxes in soils derived from sedimentary rocks along a catena formed under a semi-arid climate. Our measurements demonstrate that there are distinct and traceable pools of N in the soil and bedrock and that the fraction of rock-derived N declines downslope along the catena. We used geochemical mass balance weathering flux measurements to estimate a rock-derived N flux of 0.145 to 0.896 kg ha -1 yr -1 at the ridgecrest. We also developed independent N flux estimates using a 15 N-based isotope mixing model. While geochemical mass-balance-based estimates fell within the 95% confidence range derived from the isotope mixing model (-1.1 to 44.3 kg ha -1 yr -1 ), this range was large due to uncertainty in values for atmospheric 15 N deposition. Along the catena, N isotopes suggest a diminishing effect of rock-derived N downslope. Overall, we found that despite relatively large N pools within the saprolite and bedrock, slow chemical weathering and landscape denudation limit the influence of rock-derived N, letting atmospheric N deposition (7.1 kg ha -1 yr -1 ) and N fixation (0.9-3.1 kg ha -1 yr -1 ) dominate N inputs to this grassland ecosystem.
54 ENVIRONMENTAL SCIENCES↗
Helium Release During Fracture and Granular Fragmentation of Rocks
Geogenic Helium-4 ( 4 He) in-situ increases locally in regions of large deformation generated naturally or anthropogenically. This gas release by deformation is a potential geochemical precursor signal for subsurface deformation. To evaluate the applicability of 4 He degassing for correlating deformation in different lithologies, we conducted high force crush tests, up to 97,800 N axial load, to assess the total 4 He released during fragmentation of the rocks. We observed that the highest 4 He released occurred in the sedimentary rocks and that release correlated strongly with lithologic age and U/Th content. Microstructural changes of the pre- and post-test rocks indicate that the degree of grain size reduction relates directly to the total 4 He released during crushing. The range of in-place 4 He was calculated based XRF measurements of uranium and thorium in each lithology, with the results indicating that the majority of the trapped 4 He was not released. However, the 4 He released by deformation depended upon how the each rock deformed during deformation and the degree of grain size reduction. We postulate that 4 He precursor signals can be used to understand subsurface deformation only if geomechanical and geochemical conditions for 4 He enrichment in a lithology are met.
Deformation signals↗
Sulfur isotopes from the Paleoproterozoic Francevillian Basin record multigenerational pyrite formation, not depositional conditions
Bulk-rock sulfur isotope data from pyrite in the ~2.1 billion-year sedimentary rocks of the Francevillian Basin, Gabon, have underpinned ideas about initial oxygenation of Earth’s surface environments and eukaryote evolution. In this paper, we show, using micro-scale analytical methods, that the bulk sulfur isotope record represents progressive diagenetic modification. Our findings indicate no significant change in microbial sulfur cycling processes and seawater sulfate composition throughout that initial phase of atmosphere-ocean oxygenation of Paleoproterozoic time. This offers an alternative view of Earth system evolution during the transition from an anoxic to an oxic state and highlights the need for a judicious reappraisal of conceptual models using sulfur isotope data as primary depositional signals linked to global-scale biogeochemical processes.
58 GEOSCIENCES↗