Search NASA⌕ Search

SEARCH · Search NASA

Results for “Geoscience”

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.

At least 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Contrasting stream nitrate and sulfate response to recovery from experimental watershed acidification

Improvements in air quality have led to ecosystem recovery from acidic deposition, but the mechanisms and trajectories of this recovery are not fully understood. Here, we present long-term stream response and recovery data for paired watersheds at the Bear Brook Watershed in Maine (BBWM) during declining ambient SO 4 and NO 3 in precipitation. East Bear (EB) received ambient deposition from 1989 to 2018; West Bear (WB) received artificially elevated N+S from 1989 to 2016. The WB treatment was discontinued after 2016, the beginning of the recovery from both the experimental N+S and ambient decline. Stream SO 4 in WB gradually declined after the treatment ended, from ~147 µeq L -1 in 2010-16 to ~126 µeq L -1 in 2017-18. The declining S inputs induced desorption of SO 4 from soil phase surfaces, with stream loss far exceeding precipitation input. At the current rate of recovery, it will be many decades before the WB stream returns to pre-treatment SO 4 concentrations. In contrast, NO 3 is only weakly adsorbed in soil, and WB stream NO 3 concentrations rapidly declined from ~39 µeq L -1 in 2010-16 to ~5 µeq L -1 in 2017-18, comparable to the N-limited EB stream. The acid anions are strongly coupled to base cation chemistry in streams, and there was a distinct hysteretic response of Ca and Mg to the chronic acidification, as (Ca+Mg) increased rapidly during the initial years, followed by declining values due to depletion of the soil exchange complex. This 30-year monitoring study (1989-2019) provides insights into recovery mechanisms from acidic deposition and highlights the role of abiotic processes in soil that mediate nutrient cycling and retention. Finally, documenting the rapid response of N alongside the slower recovery for S identifies the temporal resolution necessary for other whole-watershed recovery studies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Evidence linking calcium to increased organo-mineral association in soils

Geochemical indicators are emerging as important predictors of soil organic carbon (SOC) dynamics, but evidence concerning the role of calcium (Ca) is scarce. This study investigates the role of Ca prevalence in SOC accumulation by comparing otherwise similar sites with (CaCO 3 -bearing) or without carbonates (CaCO 3 -free). We measured the SOC content and indicators of organic matter quality (C stable isotope composition, expressed as δ 13 C values, and thermal stability) in bulk soil samples. We then used sequential sonication and density fractionation (DF) to separate two occluded pools from free and mineral-associated SOC. The SOC content, mass, and δ 13 C values were determined in all the fractions. X-ray photoelectron spectroscopy was used to investigate the surface chemistry of selected fractions. Our hypothesis was that occlusion would be more prevalent at the CaCO 3 -bearing site due to the influence of Ca on aggregation, inhibiting oxidative transformation, and preserving lower δ 13 C values. Bulk SOC content was twice as high in the CaCO 3 -bearing profiles, which also had lower bulk δ 13 C values, and more occluded SOC. Yet, contrary to our hypothesis, occlusion only accounted for a small proportion of total SOC (< 10%). Instead, it was the heavy fraction (HF), containing mineral-associated organic C, which accounted for the majority of total SOC and for the lower bulk δ 13 C values. Overall, an increased Ca prevalence was associated with a near-doubling of mineral-associated SOC content. Future investigations should now aim to isolate Ca-mediated complexation processes that increase organo-mineral association and preserve organic matter with lower δ 13 C values.

58 GEOSCIENCES↗
DOE OSTI2021

Differential effects of redox conditions on the decomposition of litter and soil organic matter

Soil redox conditions exert substantial influence on biogeochemical processes in terrestrial ecosystems. Humid tropical forests are often characterized by fluctuating redox dynamics, yet how redox fluctuations affect patterns in soil versus litter decomposition and associated CO2 fluxes are not well understood in these ecosystems. We used a 13C-labeled litter addition to explicitly follow the decomposition of litter-derived vs. native soil-derived organic matter in response to four different soil redox regimes—static oxic or anoxic, and two oscillating treatments—in soil from the Luquillo Experimental Forest, Puerto Rico. We coupled this incubation experiment with high-resolution mass spectrometry analysis to characterize the preferential decomposition of specific classes of organic molecules. CO2 production from litter and soil organic matter (SOM) showed distinctly different responses to redox manipulation. The cumulative production of SOM-derived CO2 was positively correlated with the length of soil exposure to an oxic headspace (r = 0.89, n = 20), whereas cumulative 13C-litter-derived CO2 production was not linked to oxygen availability. The CO2 production rate from litter was highest under static anoxic conditions in the first half of the incubation period, and later dropped to the lowest among all redox treatments. In the consistently anoxic soils, we observed the depletion of more oxidized water-extractable organic matter (especially amino sugars, carbohydrates, and proteins) over time, suggesting that under anaerobic conditions, microbes preferentially used more oxidized litter-derived compounds during the early stages of decomposition. Results from kinetic modeling showed that more frequent anoxic exposure limited the decomposition of a slow-cycling C pool, but not a fast-cycling pool. Overall, our results demonstrate that substrate source—freshly added litter vs. native organic matter—plays an important role in the redox sensitivity of organic matter decomposition. In soil environments that regularly experience redox fluctuations, anaerobic heterotrophs can be surprisingly effective in degrading fresh plant litter.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Rock weathering controls the potential for soil carbon storage at a continental scale

Abstract As rock-derived primary minerals weather to form soil, they create reactive, poorly crystalline minerals that bind and store organic carbon. By implication, the abundance of primary minerals in soil might influence the abundance of poorly crystalline minerals, and hence soil organic carbon storage. However, the link between primary mineral weathering, poorly crystalline minerals, and soil carbon has not been fully tested, particularly at large spatial scales. To close this knowledge gap, we designed a model that links primary mineral weathering rates to the geographic distribution of poorly crystalline minerals across the USA, and then used this model to evaluate the effect of rock weathering on soil organic carbon. We found that poorly crystalline minerals are most abundant and most strongly correlated with organic carbon in geographically limited zones that sustain enhanced weathering rates, where humid climate and abundant primary minerals co-occur. This finding confirms that rock weathering alters soil mineralogy to enhance soil organic carbon storage at continental scales, but also indicates that the influence of active weathering on soil carbon storage is limited by low weathering rates across vast areas.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Influence of contrasting redox conditions on iron (oxyhydr)oxide transformation and associated phosphate sorption

Iron (oxyhydr)oxides strongly adsorb phosphate and limit its bioavailability, but interactions between phosphate and various Fe (oxyhydr)oxides are poorly constrained in natural systems. An in-situ incubation experiment was conducted to explore Fe (oxyhydr)oxide transformation and effects on phosphate sorption in soils with contrasting saturation and redox conditions. Synthetic Fe (oxyhydr)oxides (ferrihydrite, goethite and hematite) were coated onto quartz sand and either pre-sorbed with phosphate or left phosphate-free. The oxide-coated sands were mixed with natural organic matter, enclosed in mesh bags, and buried in and around a vernal pond for up to 12 weeks. Redox conditions were stable and oxic in the upland soils surrounding the vernal pond but largely shifted from Fe reducing to Fe oxidizing in the lowland soils within the vernal pond as it dried during the summer. Iron (oxyhydr)oxides lost more Fe (- 41% ± 10%) and P (- 43 ± 11%) when incubated in the redox-dynamic lowlands compared to the uplands (- 18% ± 5% Fe and - 24 ± 8% P). Averaged across both uplands and lowlands, Fe losses from crystalline goethite and hematite (- 38% ± 6%) were unexpectedly higher than losses from short range ordered ferrihydrite (- 12% ± 10%). We attribute losses of Fe and associated P from goethite and hematite to colloid detachment and dispersion but losses from ferrihydrite to reductive dissolution. Iron losses were partially offset by retention of solubilized Fe as organic-bound Fe(III). Iron (oxyhydr)oxides that persisted during the incubation retained or even gained P, indicating low amounts of phosphate sorption from solution. In conclusion, these results demonstrate that hydrologic variability and Fe (oxyhydr)oxide mineralogy impact Fe mobilization pathways that may regulate phosphate bioavailability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Reduced accrual of mineral-associated organic matter after two years of enhanced rock weathering in cropland soils, though no net losses of soil organic carbon

Abstract Enhanced rock weathering (ERW), the application of crushed silicate rock to soil, can remove atmospheric carbon dioxide by converting it to (bi) carbonate ions or solid carbonate minerals. However, few studies have empirically evaluated ERW in field settings. A critical question remains as to whether additions of crushed rock might positively or negatively affect soil organic matter (SOM)—Earth’s largest terrestrial organic carbon (C) pool and a massive reservoir of organic nitrogen (N). Here, in three irrigated cropland field trials in California, USA, we investigated the effect of crushed meta-basalt rock additions on different pools of soil organic carbon and nitrogen (i.e., mineral-associated organic matter, MAOM, and particulate organic matter, POM), active microbial biomass, and microbial community composition. After 2 years of crushed rock additions, MAOM stocks were lower in the upper surface soil (0–10 cm) of plots with crushed rock compared to unamended control plots. At the 2 sites where baseline pre-treatment data were available, neither total SOC nor SON decreased over the 2 years of study in plots with crushed rock or unamended control plots. However, the accrual rate of MAOM-C and MAOM-N at 0–10 cm was lower in plots with crushed rock vs. unamended controls. Before ERW is deployed at large scales, our results suggest that field trials should assess the effects of crushed rock on SOM pools, especially over multi-year time scales and in different environmental contexts, to accurately assess changes in net C and understand the mechanisms driving interactions between ERW and SOM cycling.

58 GEOSCIENCES↗
DOE OSTI2026

Observations of the Marine Atmospheric Boundary Layer’s Response to a Solar Eclipse

The atmospheric response to the solar eclipse of 8 April 2024 in North America is investigated with a specific focus on the marine atmospheric boundary layer (MABL). We leverage measurements collected during the Third Wind Forecast Improvement Project (WFIP3), including Doppler lidars, sonic anemometers, and thermodynamic profiler data to investigate the atmospheric response across sites that experienced partial eclipse conditions with nearly 90% obscuration. Using these measurements, we examine eclipse-induced changes in key meteorological parameters, such as temperature, wind speed, and turbulent fluxes. Most previous eclipse studies have been conducted over land, whereas this study provides new observations for both coastal and marine environments, offering additional insight into eclipse-driven variability in the MABL. The findings confirm a notable decrease in downwelling shortwave radiation during the eclipse, which results in rapid cooling of surface air. The temperature reduction ranges from $1.2^\circ \text {C}$ to $1.4^\circ \text {C}$ in coastal regions and from $0.3^\circ \text {C}$ to $0.5^\circ \text {C}$ over the ocean. This analysis suggests that the MABL’s higher thermal inertia compared to coastal regions moderates the temperature decrease during the eclipse. Wind speed exhibits a more complex behavior, as it is influenced by both the MABL and preexisting synoptic conditions. Although a reduction in wind speed is observable up to approximately 140 m above ground level (AGL) at more inland sites, at other locations closer to the coast, this reduction is constrained to the lowest 100 m AGL. Turbulence parameters retrieved from sonic anemometers, such as turbulence kinetic energy, turbulent heat flux, and friction velocity, decrease during the eclipse at coastal sites, accompanied by a brief transition of atmospheric stability from unstable to neutral or weakly stable conditions. For the open-ocean sites, the variability in turbulence statistics and atmospheric stability is minimal during the occurrence of the eclipse.

16 TIDAL AND WAVE POWER↗
DOE OSTI2020

Simulation of mineral dissolution at the pore scale with evolving fluid-solid interfaces: review of approaches and benchmark problem set

This manuscript presents a benchmark problem for the simulation of single-phase flow, reactive transport, and solid geometry evolution at the pore scale. The problem is organized in three parts that focus on specific aspects: flow and reactive transport (part I), dissolution-driven geometry evolution in two dimensions (part II), and an experimental validation of three-dimensional dissolution-driven geometry evolution (part III). Five codes are used to obtain the solution to this benchmark problem, including Chombo-Crunch, OpenFOAM-DBS, a lattice Boltzman code, Vortex, and dissolFoam. These codes cover a good portion of the wide range of approaches typically employed for solving pore-scale problems in the literature, including discretization methods, characterization of the fluid-solid interfaces, and methods to move these interfaces as a result of fluid-solid reactions. A short review of these approaches is given in relation to selected published studies. Results from the simulations performed by the five codes show remarkable agreement both quantitatively--based on upscaled parameters such as surface area, solid volume, and effective reaction rate--and qualitatively--based on comparisons of shape evolution. This outcome is especially notable given the disparity of approaches used by the codes. Therefore, these results establish a strong benchmark for the validation and testing of pore-scale codes developed for the simulation of flow and reactive transport with evolving geometries. They also underscore the significant advances seen in the last decade in tools and approaches for simulating this type of problem.

58 GEOSCIENCES↗
DOE OSTI2020

Stable and radioactive carbon isotope partitioning in soils and saturated systems: a reactive transport modeling benchmark study

This benchmark provides the first rigorous test of a three-isotope system [ 12 C, 13 C, and 14 C] subject to the combined effects of radioactive decay and both stable equilibrium and kinetic fractionation. We present a series of problems building in complexity based on the cycling of carbon in both organic and inorganic forms. The key components implement (1) equilibrium fractionation between multiple coexisting carbon species as a function of pH, (2) radioactive decay of radiocarbon with associated mass-dependent speciation demonstrating appropriate correction of the Δ 14 C value in agreement with reporting convention, and (3) kinetic stable isotope fractionation due to the oxidation of organic carbon to inorganic forms as a function of time and space in an open, through-flowing system. Participating RTM codes are CrunchTope, ToughReact, Hytec, and The Geochemist’s Workbench. Across all problem levels, simulation results from all RTMs demonstrate good agreement.

58 GEOSCIENCES↗
DOE OSTI2020

Three-dimensional fracture continuum characterization aided by surface time-domain electromagnetics and hydrogeophysical joint inversion—proof-of-concept

Efficient and safe production of hydraulically fractured reservoirs benefits from the prediction of their geometrical attributes. Geophysical methods have the potential to provide data that are sensitive to fracture geometries, alleviating the typically sparse nature of in situ reservoir observations. Moreover, surface-based methods can be logistically and economically attractive since they avoid operational interference with the injection well infrastructure. Here, we investigate the potential of the surface-based time-domain electromagnetic (EM) method. EM methods can play an important role owing to their sensitivity to injection-induced fluid property changes. Two other advantageous factors are the EM signal-enhancing effect of vertical steel-cased wells and the fact that injected proppants can be enhanced to produce a stronger electrical conductivity contrast with the reservoir’s connate fluid. Nevertheless, an optimal fracture characterization will no doubt require the integration of EM and reservoir injection and production data. We hence carry out our investigations within a hydrogeophysical parameter estimation framework where EM data and injection flow rates are combined in a fully coupled way. Given the interdisciplinary nature of coupled hydrogeophysical inverse modeling, we dedicate one section to laying out key aspects in a didactic manner.

58 GEOSCIENCES↗
DOE OSTI2020

Approximate inverse-based block preconditioners in poroelasticity

We focus on the fully implicit solution of the linear systems arising from a three-field mixed finite element approximation of Biot’s poroleasticity equations. The objective is to develop algebraic block preconditioners for the efficient solution of such systems by Krylov subspace methods. In this work, we investigate the use of approximate inverse-based techniques to decouple the native system of equations and obtain explicit sparse approximations of the Schur complements related to the physics-based partitioning of the unknowns by field type. Here, the proposed methods are tested in various numerical experiments including real-world applications dealing with petroleum and geotechnical engineering.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Microscale mechanical modeling of deformable geomaterials with dynamic contacts based on the numerical manifold method

Abstract Micromechanical modeling of geomaterials is challenging because of the complex geometry of discontinuities and potentially large number of deformable material bodies that contact each other dynamically. In this study, we have developed a numerical approach for micromechanical analysis of deformable geomaterials with dynamic contacts. In our approach, we detect contacts among multiple blocks with arbitrary shapes, enforce different contact constraints for three different contact states of separated, bonded, and sliding, and iterate within each time step to ensure convergence of contact states. With these features, we are able to simulate the dynamic contact evolution at the microscale for realistic geomaterials having arbitrary shapes of grains and interfaces. We demonstrate the capability with several examples, including a rough fracture with different geometric surface asperity characteristics, settling of clay aggregates, compaction of a loosely packed sand, and failure of an intact marble sample. With our model, we are able to accurately analyze (1) large displacements and/or deformation, (2) the process of high stress accumulated at contact areas, (3) the failure of a mineral cemented rock samples under high stress, and (4) post-failure fragmentation. The analysis highlights the importance of accurately capturing (1) the sequential evolution of geomaterials responding to stress as motion, deformation, and high stress; (2) large geometric features outside the norms (such as large asperities and sharp corners) as such features can dominate the micromechanical behavior; and (3) different mechanical behavior between loosely packed and tightly packed granular systems.

58 GEOSCIENCES↗
DOE OSTI2020

Prediction of gas hydrate saturation using machine learning and optimal set of well-logs

We report resistivity and acoustic logs are widely used to estimate gas hydrate saturation in various sedimentary systems using one of the two popular methods ((1) acoustic velocity and (2) electrical resistivity), but the limitations of these two methods are often overlooked, which include (i) well-specific calibration of empirical exponents in the electrical resistivity method, (ii) assumption of known pore morphology for gas hydrates in the acoustic velocity method, and (iii) presence of unknown mineralogy and bulk modulus terms in the acoustic velocity method. NMR-density porosity-derived gas hydrate saturation based on the analysis of the transverse magnetization relaxation time (T2) is considered the most precise method, but acquisition of NMR-based logs is limited at relatively recent drilled sites; additionally, its use in conventional oil and gas reservoirs is not that common due to higher cost and operational deployment limitations associated with acquiring NMR well-logs. This study proposes a new method that predicts gas hydrate saturation (S h ) for any well using porosity, bulk density, and compressional wave (P wave) velocity well-logs with neural network (or stochastic gradient descent regression) without any well-specific calibration and/or other aforementioned shortcomings of the existing methods. The method is developed by examining the underlying dependency between S h and different combinations of well-logs, chosen from 6 routine logs, with 12 different machine learning (ML) algorithms. The accuracy of the proposed method in predicting S h is ~ 84%, which is better than the accuracy of seismic and electrical resistivity methods (≤ 75%) per the results reported by three different studies. The robustness of the method in the specific case of permafrost-associated gas hydrates is demonstrated with well-log data from two wells drilled on the Alaska North Slope.

58 GEOSCIENCES↗
DOE OSTI2020

A model for discrete fracture-clay rock interaction incorporating electrostatic effects on transport

Abstract A model based on the code CrunchClay is presented for a fracture-clay matrix system that takes electrostatic effects on transport into account. The electrostatic effects on transport include those associated with the development of a diffusion potential as captured by the Nernst-Planck equation, and the formation of a diffuse layer bordering negatively charged clay particles within which partial anion exclusion occurs. The model is based on a dual continuum formulation that accounts for diffuse layer and bulk water pore space, providing a more flexible framework than is found in the classical mean electrostatic potential models. The diffuse layer model is obtained by volume averaging ion concentrations in the Poisson-Boltzmann equation, but also includes the treatment of longitudinal transport within this continuum. The calculation of transport within the bulk and diffuse layer porosity is based on a new formulation for the Nernst-Planck equation that considers averaging of diffusion coefficients and accumulation factors at grid cell interfaces. Equations for function residuals and the associated Jacobian matrix are presented such that the system of nonlinear differential-algebraic equations can be solved with Newton’s method. As an example, we consider a 2D system with a single discrete fracture within which flow and advective transport occurs that is coupled to diffusion in the clay-rich matrix. The simulation results demonstrate the lack of retardation for anions (e.g., 36 Cl − ) of the contaminant plume within the fracture flow system because they are largely excluded from the charged clay rock, while the migration of cations (e.g., 90 Sr ++ ) is more strongly attenuated. The diffusive loss of divalent cations in particular from the fracture is accentuated by their accumulation in the diffuse layer within the clay-rich matrix.

58 GEOSCIENCES↗
DOE OSTI2022

Application of machine learning to characterize gas hydrate reservoirs in Mackenzie Delta (Canada) and on the Alaska north slope (USA)

Here, artificial neural network-trained models were used to predict gas hydrate saturation distributions in permafrost-associated deposits in the Eileen Gas Hydrate Trend on the Alaska North Slope (ANS), USA and at the Mallik research site in the Beaufort-Mackenzie Basin, Northwest Territories, Canada. The database of Logging-While-Drilling (LWD) and wireline logs collected at five wells (Mount Elbert, Ignik Sikumi, and Kuparuk 7–11–12 wells at ANS, plus 2L-38 and 5L-38 wells at the Mallik research site) includes more than 10,000 depth points, which were used for training, validation, and testing the machine learning (ML) models. Data used in training the ML models include the well logs of density, porosity, electrical resistivity, gamma radiation, and acoustic wave velocity measurements. Combinations of two or three out of these five well logs were found to reliably predict the gas hydrate saturation with accuracy varying between 80 and 90% when compared to the gas hydrate saturations derived from Nuclear Magnetic Resonance (NMR)-based technique. The ML models trained on data from three ANS wells achieved high fidelity predictions of gas hydrate saturation at the Mallik site. The results obtained in this study indicate that ML models trained on data from one geological basin can successfully predict key reservoir parameters for permafrost-associated gas hydrate accumulations within another basin. A generalized approach for selecting a well log combination that can improve model accuracy is discussed. Overall, the study outcome supports earlier work demonstrating that ML models trained on non-NMR well logs are a viable alternative to physics-driven methods for predicting gas hydrate saturations.

58 GEOSCIENCES↗
DOE OSTI2022

TOUGH3-FLAC3D: a modeling approach for parallel computing of fluid flow and geomechanics

The recent development of the TOUGH3 code allows for a faster and more reliable fluid flow simulator. At the same time, new versions of FLAC3D are released periodically, allowing for new features and faster execution. In this paper, we present the first implementation of the coupling between TOUGH3 and FLAC3Dv6/7, maintaining parallel computing capabilities for the coupled fluid flow and geomechanical codes. We compare the newly developed version with analytical solutions and with the previous approach, and provide some performance analysis on different meshes and varying the number of running processors. Finally, we present two case studies related to fault reactivation during CO 2 sequestration and nuclear waste disposal. The use of parallel computing allows for meshes with a larger number of elements, and hence more detailed understanding of thermo-hydro-mechanical processes occurring at depth.

58 GEOSCIENCES↗
DOE OSTI2023

Impact of artificial topological changes on flow and transport through fractured media due to mesh resolution

Abstract We performed a set of numerical simulations to characterize the interplay of fracture network topology, upscaling, and mesh refinement on flow and transport properties in fractured porous media. We generated a set of generic three-dimensional discrete fracture networks at various densities, where the radii of the fractures were sampled from a truncated power-law distribution, and whose parameters were loosely based on field site characterizations. We also considered five network densities, which were defined using a dimensionless version of density based on percolation theory. Once the networks were generated, we upscaled them into a single continuum model using the upscaled discrete fracture matrix model presented by Sweeney et al. (2019). We considered steady, isothermal pressure-driven flow through each domain and then simulated conservative, decaying, and adsorbing tracers using a pulse injection into the domain. For each simulation, we calculated the effective permeability and solute breakthrough curves as quantities of interest to compare between network realizations. We found that selecting a mesh resolution such that the global topology of the upscaled mesh matches the fracture network is essential. If the upscaled mesh has a connected pathway of fracture (higher permeability) cells but the fracture network does not, then the estimates for effective permeability and solute breakthrough will be incorrect. False connections cannot be eliminated entirely, but they can be managed by choosing appropriate mesh resolution and refinement for a given network. Adopting octree meshing to obtain sufficient levels of refinement leads to fewer computational cells (up to a 90% reduction in overall cell count) when compared to using a uniform resolution grid and can result in a more accurate continuum representation of the true fracture network.

58 GEOSCIENCES↗
DOE OSTI2024

A mineral precipitation model based on the volume of fluid method

A novel volume of fluid method is presented for mineral precipitation coupled with fluid flow and reactive transport. The approach describes the fluid-solid interface as a smooth transitional region, which is designed to provide the same precipitation rate and viscous drag force as a sharp interface. Specifically, the governing equation of mineral precipitation is discretized by an upwind scheme, and a rigorous effective viscosity model is derived around the interface. The model is validated against analytical solutions for mineral precipitation in channel and ring-shaped structures. It also compares well with interface tracking simulations of advection-diffusion-reaction problems. Here, the methodology is finally employed to model mineral precipitation in fracture networks, which is challenging due to the low porosity and complex geometry. Compared to other approaches, the proposed model has a concise algorithm and contains no free parameters. In the modeling, only the pore space requires meshing, which improves the computational effciency especially for low-porosity media.

58 GEOSCIENCES↗