Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Silicates”

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 1,729 records · Page 96
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Ab Initio Structures and Energetics of Hydrated Flat and Terrace-Step Surfaces of Forsterite (Mg 2 SiO 4 )

Forsterite (Mg 2 SiO 4 ), a model divalent metal silicate mineral, has been extensively studied in the context of mineral carbonation. Although dissolution is a key step in this process, the mechanisms by which forsterite dissolves under high CO 2 conditions remain poorly understood. Atomistic simulations could aid in exploring these mechanisms, but it is essential first to understand the structures and energetics of the relevant forsterite surfaces. We present an ab initio study of the structure and surface energy at 0 K of the flat $(010), (110), (001), (111), (021), (101)$ and $(120)$ faces of forsterite using the density functional PBE Hamiltonian and a plane-wave basis set. Dry surfaces became stabilized upon hydration through the formation of bonds between surface Mg and O from water, as well as by the formation of hydrogen bonds. According to surface energy values, the stability order of the hydrated forsterite faces was found to be $(120) < (101) < (021) < (111) < (001) < (110) < (010)$. We also investigated the energetics of the terrace-step $(0\bar{41})$ surface as a model site for forsterite dissolution. Among all the facets, the $(0\bar{41})$ surface is the least stable termination in water. Hydration of Mg atoms on the $(0\bar{41})$ surface increases their susceptibility to dissolution. The presence of a step and its hydration destabilizes the terraces, making step retreat more likely than a dissolution front advancing along the [010] direction. This research will support future simulations to investigate forsterite dissolution in water under CO 2 -rich conditions.

PBE Hamiltonian↗
DOE OSTI2026

Modeling an S N 2 Reaction Mechanism for Hydrolysis of Siloxane Linkages with Density Functional Theory under Basic Conditions and Implications for Dissolution of Quartz

An extended silicate molecular cluster was used to perform density functional theory calculations to determine if an S N 2 mechanism, where OH – attacks a Q 1 Si coupled with Na + charge balancing an adjacent siloxane bridging O atom, could explain the observed energy of activation of quartz dissolution under basic conditions.

chemical calculations↗
DOE OSTI2026

Decadal-Scale Hydrothermal Alteration of Basalt: Implications for Long-Term Reactivity in Subsurface Technologies

Global demands to secure critical mineral supply chains and managing water resources require innovative resource management strategies. Subsurface basalt reservoirs offer a dual-purpose application, where CO2 enhanced critical metal recovery can be coupled with water remediation strategies for the extraction of economically important and environmentally relevant elements. This study utilizes static batch experiments conducted on Columbia River Basalt Group samples to evaluate long-term fluid–rock interactions relevant to coupled subsurface applications. Samples were reacted under elevated temperature and pressure conditions, for 296 and 4,061 days, providing a unique dataset to assess geochemical evolution, element partitioning, and aqueous chemistry over decadal timescales. The alteration was dominated by silicate, aluminosilicate, and oxide phases. Progressive hydrothermal alteration increased the cation exchange capacity through the formation of secondary minerals (e.g. clays and zeolites), influencing solute retention, ion exchange behavior, and aqueous chemistry. Grain-scale and bulk analyses revealed selective cation release and retention patterns, establishing a mechanistic baseline for predicting element mobility. These geochemical constraints inform the prospective use of CO2-enhanced metal recovery for critical resources (e.g. Co, Ni, Mn, Mg, Ti, Fe) and associated water-quality impacts, providing a foundation for the design and management of engineered fluid–rock systems in basalt reservoirs.

Gibson, Marie N.↗
DOE OSTI2025

Electrocoagulation Combined with Ultrafiltration Membranes as Pretreatment for RO Desalination of Synthetic Cooling Tower Blowdown Water

Electrocoagulation (EC), an electrochemical water treatment process, is commonly used to remove particulate and colloidal matter from water. Here, in this study, we demonstrate that EC, when coupled with membrane filtration, is also capable of removing dissolved species, such as Ca + , Mg 2+ , and SiO x . The removal of such species is important for downstream membrane-based desalination treatment of the water that can suffer from reduced performance due to membrane scaling. Here, we describe how EC can be combined with a low-pressure membrane (LPM) system to offer efficient (and potentially universal) pretreatment for downstream membrane desalination. Synthetic water, simulating cooling tower blowdown (CTBD) with elevated concentrations of hardness and silicates (Ca: 418 ppm, Mg: 63 ppm, SiO 2 : 50 ppm) is treated using EC coupled to ultrafiltration (UF) to remove Mg up to 30 ± 1%, Ca up to 29 ± 1%, and silica up to 99 ± 1%. We evaluated the effectiveness of the EC-UF pretreatment system in reducing downstream RO scaling using thermodynamic modeling to predict the saturation index (SI) at the RO membrane/water interface. An SI value below zero (SI < 0) indicates under-saturated conditions (with respect to a particular mineral) where mineral scaling does not take place, which correlates with improved water recovery. Our findings suggest that the EC-UF pretreatment system was able to increase water recovery by up to 30%, compared to 0% recovery without pretreatment, under optimal conditions (feed solution pH of 7 and an EC charge loading of 1800 C/L). Finally, we conducted an economic analysis showing that implementing an EC-UF system for CTBD water could yield a cost benefit of up to $\$$14.13 per m 3 compared to direct brine disposal.

cooling tower blowdown↗
DOE OSTI2023

In Situ Imaging of Faujasite Surface Growth Reveals Unique Pathways of Zeolite Crystallization

Zeolite crystallization occurs by complex processes involving a variety of possible mechanisms. The sol gel media used to prepare zeolites leads to heterogeneous mixtures of solution and solid states with diverse solute species. At later stages of zeolite synthesis when growth occurs predominantly from solution, classical two-dimensional nucleation and spreading of layers on crystal surfaces via the addition of soluble species is the dominant pathway. At earlier stages, these processes occur in parallel with nonclassical pathways involving crystallization by particle attachment (CPA). The relative roles of solution- and solid-state species in zeolite crystallization have been a subject of debate. Here, in this work, we investigate the growth mechanism of a commercially relevant zeolite, faujasite (FAU). In situ atomic force microscopy (AFM) measurements reveal that supernatant solutions extracted from a conventional FAU synthesis at various times do not result in growth, indicating that FAU growth predominantly occurs from the solid state through a disorder-to-order transition of amorphous precursors. Elemental analysis shows that supernatant solutions are significantly more siliceous than both the original growth mixture and the FAU zeolite product; however, in situ AFM studies using a dilute clear solution with a lower Si/Al ratio revealed three-dimensional growth of surfaces that is distinct from layer-by-layer and CPA pathways. This unique mechanism of growth differs from those observed in studies of other zeolites. Given that relatively few zeolite frameworks have been the subject of mechanistic investigation by in situ techniques, these observations of FAU crystallization raise the question whether its growth pathway is characteristic of other zeolite structures.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Effects of Cofeeding Hydrogen on Propane Dehydrogenation Catalyzed by Isolated Iron Sites Incorporated into Dealuminated BEA

Iron sites dispersed on nonacidic siliceous supports have been reported to be catalytically active for propane dehydrogenation (PDH), yet the precise relationship between site structure and catalytic activity remains elusive. Here, this study provides a comprehensive understanding of the catalytic performance of iron supported on dealuminated BEA (DeAlBEA) zeolites for PDH. Using XAS, UV-vis, and IR spectroscopy of adsorbed pyridine and deuterated acetonitrile, it was found that, at an Fe/Al 0 of 0.04, isolated Fe sites form. These isolated sites exhibit a forward rate of PDH of 213 mol propene/mol Fe·h at 823 K and a feed containing 15 kPa propane. When 15 kPa of H 2 is added to the feed, the forward rate of PDH rises to 391 mol of propene/mol of Fe·h. In both cases, the propene selectivity is over 99%. IR spectroscopy of d 3 -acetonitrile suggests that the open Lewis acid site ((-Si-O-) 2 Fe 3+ -OH) serves as the active site responsible for PDH, while Brønsted acid sites (≡Fe 3+ -O(H)-Si≡) contribute to propane cracking with increasing Fe/Al 0 ratios. Kinetic analysis of the effects of H 2 addition to the propane feed on PDH kinetics shows that H 2 enhances the activity of 0.04FeDeAlBEA primarily by enhancing the strength of the propane adsorption.

Alghannam, Afnan [University of California, Berkel↗
DOE OSTI2025

Interactions of Polar and Nonpolar Groups of Alcohols in Zeolite Pores

Understanding the quantitative interactions among zeolite pore walls, Bro̷nsted acid sites, and molecules with both polar and nonpolar regions is essential for scoping out the potential of zeolites as sorbents and catalysts. Purely siliceous zeolites (MFI and Beta in the present study) are hydrophobic, whereas those containing aluminum are considered hydrophilic, preferentially adsorbing organic molecules even in aqueous environments. To characterize these interactions, we use primary alcohols of increasing molecular weight, quantifying their specific interactions in the confined pore space of the alkyl (CH x ) and OH groups. Three types of interactions were identified: (i) alkyl CH x groups interacting with the zeolite pore walls (approximately 10 kJ mol −1 per carbon), (ii) alcohol OH groups interacting with the pore walls (30−35 kJ mol −1 ), and (iii) alcohol OH groups interacting with Bro̷nsted acid sites (37 kJ mol −1 ). All three interactions were well mirrored by computational simulations. The contribution of the alkyl CH x groups was inferred from the incremental increase in sorption enthalpy with increasing molecular weight; the interaction strength of the OH groups was determined by extrapolating the global adsorption enthalpy of the alcohols to a hypothetical OH group without an alkyl group. This value was identical to the adsorption enthalpy of water. The experiments demonstrated that only water has an adsorption enthalpy on zeolite pore walls lower than its condensation enthalpy (30−35 kJ mol −1 vs 45 kJ mol −1 ), limiting the concentration of water that can be adsorbed.

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

Quantifying Fluid-Wettable Effective Pore Space in the Utica and Bakken Oil Shale Formations

Combined (ultra-) small angle neutron scattering measurements [(U)SANS] and a contrast matching technique were employed to quantify the porosity and pore size distribution from 1 nm to 10 μm and to differentiate accessible (open) pores and inaccessible (closed) pores with respect to organophilic and hydrophilic fluids for two Utica and two Bakken shale samples. Overall, the results indicate that around 40–70% of the pores in the Utica oil shales (mixed carbonate mudstone) are accessible to oil, and 34–37% of the pore surfaces are water-wet. In contrast, the Bakken oil shales (mixed siliceous mudstone and carbonate/siliceous mudstone), which have high total organic carbon contents, have a higher proportion of isolated pore space that is not preferentially wet by oil or water, with less than 36% of the pores accessible to both fluids. In addition, for both formations, pores less than 3 nm in diameter are not oil accessible (organic matter related) but water accessible (clay tactoids related).

58 GEOSCIENCES↗
DOE OSTI2020

Direct Observation of Shock‐Induced Disordering of Enstatite Below the Melting Temperature

Abstract We report in situ structural measurements of shock‐compressed single crystal orthoenstatite up to 337 ± 55 GPa on the Hugoniot, obtained by coupling ultrafast X‐ray diffraction to laser‐driven shock compression. Shock compression induces a disordering of the crystalline structure evidenced by the appearance of a diffuse X‐ray diffraction signal at nanosecond timescales at 80 ± 13 GPa on the Hugoniot, well below the equilibrium melting pressure (>170 GPa). The formation of bridgmanite and post‐perovskite have been indirectly reported in microsecond‐scale plate‐impact experiments. Therefore, we interpret the high‐pressure disordered state we observed at nanosecond scale as an intermediate structure from which bridgmanite and post‐perovskite crystallize at longer timescales. This evidence of a disordered structure of MgSiO 3 on the Hugoniot indicates that the degree of polymerization of silicates is a key parameter to constrain the actual thermodynamics of shocks in natural environments.

Hernandez, J. ‐A.↗
DOE OSTI2021

A Long-Lived Planetesimal Dynamo Powered by Core Crystallization

The existence of numerous iron meteorite groups indicates that some planetesimals underwent melting that led to metal-silicate segregation, sometimes producing metallic cores. Meteorite paleomagnetic records suggest that crystallization of these cores generated dynamo magnetic fields. In this work, we describe the magnetic history of the partially differentiated IIE iron meteorite parent body. This is the first planetesimal for which we have a time-resolved paleomagnetic record constrained by 40 Ar/ 39 Ar chronometry spanning several tens of million years (Ma). We find that the core of the IIE parent body generated a dynamo, likely powered by core crystallization, starting before 78 ± 13 Ma after solar system formation and lasting at least 80 Ma. Such extended core crystallization suggests that the core composed a substantial fraction of the body (≳ 13%–19% core-to-body radius ratio depending on the body’s radius), indicating efficient core formation within some partially differentiated planetesimals.

58 GEOSCIENCES↗
DOE OSTI2020

Femtosecond X-ray Diffraction of Laser-shocked Forsterite (Mg 2 SiO 4 ) to 122 GPa

The response of forsterite, Mg 2 SiO 4 , under dynamic compression is of fundamental importance for understanding its phase transformations and high-pressure behavior. Here, we have carried out an in situ X-ray diffraction study of laser-shocked polycrystalline and single-crystal forsterite from 19 to 122 GPa using the Matter in Extreme Conditions end-station of the Linac Coherent Light Source. Under laser-based shock loading, forsterite does not transform to the high-pressure equilibrium assemblage of MgSiO 3 bridgmanite and MgO periclase, as has been suggested previously. Instead, we observe forsterite and forsterite III, a metastable polymorph of Mg 2 SiO 4 , coexisting in a mixed-phase region from 33 to 75 GPa for both polycrystalline and single-crystal samples. Densities inferred from X-ray diffraction are consistent with earlier gas-gun shock data. At higher stress, the response is sample-dependent. Polycrystalline samples undergo amorphization above 79 GPa. For [010]- and [001]-oriented crystals, a mixture of crystalline and amorphous material is observed to 108 GPa, whereas the [100]-oriented forsterite adopts an unknown phase at 122 GPa. The first two sharp diffraction peaks of amorphous Mg 2 SiO 4 show a similar trend with compression as those observed for MgSiO 3 in both recent static- and laser-driven shock experiments. This study provides new insight into the transformation of forsterite under nanosecond-duration shock loading. This work emphasizes the importance of formation of metastable phases along the Hugoniot and adds to evidence that the 300-K single-crystal diamond anvil cell experiments have relevance for understanding structures formed under shock compression. In particular, the metastable phase forsterite III has now been shown to form under dynamic compression from 10s to 100s of nanoseconds as well as under 300-K static compression. Upon compression to higher pressures, Mg 2 SiO 4 transforms to an amorphous phase. These results have broad relevance for understanding the behavior of silicates under dynamic compression.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Enhanced Visibility of Subduction Slabs by the Formation of Dense Hydrous Phase A

Abstract Phase A (Mg 7 Si 2 O 8 (OH) 6 ) is one of the important dense hydrous magnesium silicates in subducting slab, because it forms after the breakdown of antigorite serpentine and could be the dominant hydrous phase in the upper‐mantle deep slab for the water transportation into deep Earth. In this study, the compressional ( P ) and shear ( S ) wave velocities of phase A were measured at simultaneous pressure and temperature conditions up to 11 GPa and 1073 K. Combined with elastic properties of olivine and pyroxenes, we calculate the hydration effect on the velocities of harzburgite lithology in cold subduction zones throughout the depth ranges where phase A is thermodynamically stable. Our calculations suggest that the hydration increases both P ‐ and S ‐wave velocities of harzburgite; at ∼5 wt% hydration, its seismic detectability is enhanced by 1%–1.5% in velocity contrasts relative to its anhydrous counterpart.

Cai, Nao↗
DOE OSTI2023

Corrosive Influence of Carbon Dioxide on Crack Initiation in Quartz: Comparison With Liquid Water and Vacuum Environments

The stimulation of crack growth in quartz and siliceous materials by injecting carbon dioxide (CO 2 ) represents a key technology in long-term carbon storage and in the development of natural gas wells. While this technology is widely used, the molecular impact of CO 2 interactions on the solid matrix is only incompletely understood. In this work, we employ reactive molecular dynamics simulations to study how the CO 2 fluid environment affects the mechanical properties of pre-cracked single-crystal quartz. The thermodynamic conditions of interest are those relevant to subsurface reservoirs. We report how structural properties of quartz—bond length distribution and crack tip shape—evolve upon introduction of a fluid. These properties are directly related to macroscopic quantities of the global stress-strain curves, thus reaffirming the inherent coupling across multiple scales for fluid-solid interactions in the subsurface. Here, we find that CO 2 reduces the fracture toughness of quartz by 12.1% compared to that of quartz in vacuum, thereby promoting crack growth and enhancing fluid transport in the subsurface.

58 GEOSCIENCES↗
DOE OSTI2024

Dynamic Rupture Simulations of Caldera Collapse Earthquakes: Effects of Wave Radiation, Magma Viscosity, and Evidence of Complex Nucleation at Kilauea 2018

All instrumented basaltic caldera collapses have generated M w > 5 very long period earthquakes. However, previous studies of source dynamics have been limited to lumped models treating the caldera block as rigid, leaving open questions related to how ruptures initiate and propagate around the ring fault, and the seismic expressions of those dynamics. Here, we present the first 3D numerical model capturing the nucleation and propagation of ring fault rupture, the mechanical coupling to the underlying viscoelastic magma, and the associated seismic wavefield. We demonstrate that seismic radiation, neglected in previous models, acts as a damping mechanism reducing coseismic slip by up to half, with effects most pronounced for large magma chamber volume/ring fault radius or highly compliant crust/compressible magma. Viscosity of basaltic magma has negligible effect on collapse dynamics. In contrast, viscosity of silicic magma significantly reduces ring fault slip. We use the model to simulate the 2018 Kilauea caldera collapse. Three stages of collapse, characterized by ring fault rupture initiation and propagation, deceleration of the downward-moving caldera block and magma column, and post-collapse resonant oscillations, in addition to chamber pressurization, are identified in simulated and observed (unfiltered) near-field seismograms. A detailed comparison of simulated and observed displacement waveforms corresponding to collapse earthquakes with hypocenters at various azimuths of the ring fault reveals a complex nucleation phase for earthquakes initiated on the northwest. Our numerical simulation framework will enhance future efforts to reconcile seismic and geodetic observations of caldera collapse with conceptual models of ring fault and magma chamber dynamics.

58 GEOSCIENCES↗
DOE OSTI2024

A New Ecosystem Model for Arctic Phytoplankton Phenology From Ice‐Covered to Open‐Water Periods: Implications for Future Sea Ice Retreat Scenarios

Abstract To predict ecosystem change in the Arctic Ocean, understanding Arctic phytoplankton phenology is essential. We develop a marine ecosystem model focusing on phytoplankton dynamics and the competition for multiple resources based on knowledge from in situ data obtained in the Chukchi Sea. The model is designed to include the unique ecological characteristics of the Arctic Ocean. To the best of our knowledge, this is the first model successfully simulating the current paradigm for Arctic phytoplankton phenology, including both under‐ice blooms and succession of phytoplankton groups. Sensitivity experiments show that the dominance of diatoms can be sustained by lateral transport of a high‐silicate water mass (Pacific Winter Water). Experiments for future scenarios show that sea ice retreat drives a decrease in large diatoms and an increase in the other relatively small non‐diatom groups. These results suggest that dominant phytoplankton groups can be shifted from diatoms to non‐diatoms under ongoing Arctic Ocean sea ice retreat.

Choi, Jang‐Geun↗
DOE OSTI2025

The Effect of Hydrothermal Alteration and Microcracks on Hydraulic Properties and Poroelastic Deformation: A Case Study of the Blue Mountain Geothermal Field

Abstract Geothermal energy plays a vital role in decarbonizing electricity and heat supply. Effective utilization of geothermal resources hinges on identifying or generating permeable reservoir zones and understanding how effective pressure variations affect fluid circulation and reservoir properties by poroelastic deformation. Hydrothermal alteration can modify the petrophysical properties of geothermal reservoir rocks, which may increase or decrease its productivity. Understanding these alteration effects is essential to predict and optimize long‐term sustainable geothermal operations. Here, we investigate the impact of hydrothermal alteration on poroelastic and hydraulic properties of diverse lithologies in a series of deformation tests performed at several confining (0–80 MPa) and pore pressure (10–30 MPa) levels. Experimental results of hydrothermally altered dikes and phyllites obtained from the Blue Mountain geothermal field (Nevada, USA) are compared to thermally cracked La Peyratte granite (France) and correlated with petrophysical properties, mineral composition, and microstructures. Argillic alteration of dikes increases porosity and storage capacity but lowers thermal conductivity and increases pore compressibility. Conversely, silicate precipitation in phyllites increases stiffness and thermal conductivity but also reduces porosity and permeability. Experimentally determined effective pressure coefficients range from 0.1 to 0.9, differ for permeability and volumetric strain and decrease with increasing effective pressure. The presence of compliant microcracks and crack‐like pores significantly increases the stress sensitivity of La Peyratte granite and silicified phyllites. This study demonstrates how thermal and chemical alteration impacts poromechanical and petrophysical characteristics of geothermal targets, which ultimately govern reservoir stability and subsidence, induced seismicity as well as fluid and heat extraction efficiency during geothermal operations.

Schuster, Valerian [Helmholtz Centre Potsdam GFZ G↗
DOE OSTI2025

Composition and Plume Gas Interaction Control Iron Fractional Solubility More Than Particle Size in Volcanic Ash: Implications for Fertilization of the North Atlantic

Deposition of volcanic ash is thought to impact marine biogeochemical cycling by adding soluble iron (Fe) to the surface ocean. The magnitude of this input is a function of the amount of ash deposited, the total Fe content in the ash, and ash‐derived Fe's fractional solubility. However, the relative importance of chemical composition, acidic processing by the volcanic plume, and ash particle size in determining solubility is unclear. We paired an aerosol leach meant to provide an upper limit for fractional Fe solubility with chemical analyses of ash from the Cumbre Vieja (CV) and La Soufrière eruptions, which both impacted the North Atlantic in 2021. Fe in the ash samples is <6% soluble, but Fe fractional solubility in CV ash is approximately triple that of La Soufrière ash. Compared to La Soufrière, a larger proportion of the Fe in CV ash is in silicate rather than oxide minerals, which release more soluble Fe. Elevated levels of surficial fluorine (F) also suggest that CV ash was subjected to a more fluorine‐rich eruption plume and underwent more acidic processing. Particle size does not appear to be a primary control on Fe release. We estimate that the CV eruption had a much larger impact on dissolved Fe (DFe) concentration in the surface ocean than the La Soufrière eruption because of differences in soluble Fe content and particle deposition velocity. These differences may help explain why some eruptions elicit a biological response in the ocean while others do not.

Elliott, Hope E. [Univ. of Miami, FL (United State↗
DOE OSTI2026

Model‐Based Interpretation of Solute Exports and Carbon Partitioning During Shale Weathering in a Mountainous Hillslope

The weathering of sedimentary rocks in high-elevation catchments influences freshwater quality and the global carbon cycle. While individual biogeochemical mechanisms involved in this process are relatively well understood, quantifying their contributions to solute export and carbon fluxes under natural, transient conditions remains challenging. Here, we implement a numerical multidimensional and multiphase model to simulate coupled hydrological and biogeochemical processes in a shale-underlain, snow-dominated hillslope in the Rocky Mountains, Colorado. The model captures the dynamic interplay between soil respiration, mineral weathering, and climate-driven hydrological forcing, reproducing observed soil CO 2 dynamics, groundwater chemistry, and subsurface flow. Our results reveal that seasonal snowmelt enhances carbonate weathering by promoting the infiltration of CO 2 -rich water to depth, while pyrite oxidation is primarily sensitive to low water saturation that facilitates O 2 diffusion through the regolith. Topography modulates the spatial distribution of shale weathering, as steeper slopes enhance lateral drainage, favoring the delivery of reactants to greater depths. While shale weathering at our site acts as a transient carbon sink, with silicates and carbonates buffering acidity and promoting atmospheric CO 2 consumption (1% of soil-derived CO 2 ), the exported dissolved inorganic carbon is predominantly geogenic (∼73%). Consequently, when accounting for long-term marine carbonate precipitation. The current weathering regime represents a net source of carbon to the atmosphere. The oxidation of pyrite and petrogenic organic carbon together release approximately 0.9 mol·m −2 ·yr −1 of CO 2 . Our findings highlight the role of topography, hydroclimate, and the coupling between acid-base reactions in shaping the carbon balance and the solute exports in mountainous critical zones.

carbon cycling↗