Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geochemical Modeling”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

Hydrothermal systems as environments for the emergence of life

Analysis of the chemical disequilibrium provided by the mixing of hydrothermal fluids and seawater in present-day systems indicates that organic synthesis from CO2 or carbonic acid is thermodynamically favoured in the conditions in which hyperthermophilic microorganisms are known to live. These organisms lower the Gibbs free energy of the chemical mixture by synthesizing many of the components of their cells. Primary productivity is enormous in hydrothermal systems because it depends only on catalysis of thermodynamically favourable, exergonic reactions. It follows that hydrothermal systems may be the most favourable environments for life on Earth. This fact makes hydrothermal systems logical candidates for the location of the emergence of life, a speculation that is supported by genetic evidence that modern hyperthermophilic organisms are closer to a common ancestor than any other forms of life. The presence of hydrothermal systems on the early Earth would correspond to the presence of liquid water. Evidence that hydrothermal systems existed early in the history of Mars raises the possibility that life may have emerged on Mars as well. Redox reactions between water and rock establish the potential for organic synthesis in and around hydrothermal systems. Therefore, the single most important parameter for modelling the geochemical emergence of life on the early Earth or Mars is the composition of the rock which hosts the hydrothermal system.

Non-NASA Center↗
DOE OSTI2025

Development of a kinetic-thermodynamic model for lime-stabilization of Na-bentonite

This study presents the first kinetic model to predict the solid and pore solution composition of Na-bentonite clay reacting with slaked lime over a period of 720 days. The model successfully accounts for most experimental data using a single kinetic rate constant. The following sequence of reactions was predicted by the model: initial rapid dissolution of portlandite within the first 7 days, leading to a decrease in pH and dissolved calcium, and concurrent formation of calcium silicate hydrates (C-S-H: jennite), calcium aluminate hydrate (C-A-H: C₄AH₁₃), calcium aluminosilicate hydrates (stratlingite) and hydrotalcite. After 7 days, jennite and stratlingite are predicted to transform into tobermorite-II, contributing to strength development up to 28 days. From 28 to 90 days, continued montmorillonite dissolution is predicted, along with minor formation of ettringite, partial tobermorite-II dissolution, and precipitation of secondary phases such as albite and talc. Experimentally, portlandite dissolution was confirmed by TGA and XRD and found to be complete within 7 days, in agreement with model predictions. However, other predicted solid-phase transformations (e.g., tobermorite-II formation and dissolution, ettringite, albite, and talc formation) could not be conclusively verified through experimental techniques. Aqueous phase measurements confirmed that the pH and Ca trends in solution, and that equilibrium was reached by 90 days.

Chemical kinetics↗
DOE OSTI2024

Effect of Anoxic Iron Corrosion on WIPP Brine Geochemistry FY23 Final Report (U)

A 280-day study was completed to evaluate the effect of zero-valent iron (Fe 0 ) on the Waste Isolation Pilot Plant (WIPP) brine geochemistry under anticipated reducing conditions. Hydrogen (H 2 ) gas is expected to be present in the repository after closure due to the anoxic corrosion of a vast quantity of iron contained in the waste forms disposed at WIPP; therefore, a background argon atmosphere containing H 2 was chosen for this study. WIPP groundwater brine pH and E h will impact the mobility and fate of plutonium within the repository. Modeling and laboratory results for Castile WIPP brine indicate that equilibrium fa values relative to the standard hydrogen electrode (SHE) are 40 mV more reducing (i.e., more negative) than those for Salado WIPP brine (-480 mV vs. -440 mV, respectively) because of the higher pH of the Castile brine (pH 9 .3 for Castile vs. pH 8.8 for Salado). The E h and pH data were corrected for the effects of high ionic strength. The experimental results for both brines are consistent with thermodynamic predictions using OLI Systems' Mixed Solvent Electrolyte chemical equilibrium model. The measured and corrected pH and E h data from this study are provided in Table ES-I and Table ES-2, respectively. The experimental study, with four test conditions in triplicate, was performed in a dual glovebox with a nominally 3 vol.% H 2 in argon atmosphere (target H 2 range: 3 ± I vol.%). Simulants containing MgO only ( experimental control) and MgO+Fe 0 (WIPP base case) were prepared for both the Salado and Castile brines. MgO was included in all simulants to account for the use of bulk magnesium oxide in the WIPP repository. Fe 0 was included in some simulants to incorporate the effects of the anoxic corrosion of iron and in-situ hydrogen generation in the study. The brine compositions were developed by Sandia National Laboratory (SNL; Xiong, 2008) and have been used in previous WIPP evaluations. The test method (agitation, etc.) is partially based on ASTM D3987-12. Twelve rounds of periodic measurements of pH and E h were performed over the course of the study. Chemical analysis results for liquids and solids (ICP-MS, ICP-ES, IC Anion, TIC, SEM-EDX) are consistent with the pH, E h , and thermodynamic modeling results. This study included the following conditions that deviate from anticipated post-closure conditions following brine intrusion, but were selected to facilitate bench-scale testing to validate modeling of pH and E h for the post-closure WIP P repository: an anoxic glove box atmosphere containing ≤ 4 vol. % H 2 vs. substantially higher H 2 gas concentrations assumed in the WIPP Performance Assessment (PA); a significantly higher liquid-to-solid test ratio compared to the much lower phase ratio anticipated in the WIP P repository; agitation of the simulant bottles to maximize mass transfer; and finally the use of Fe 0 reagents having a much greater surface area than expected in the WIP P repository. Non-representative conditions were chosen for various reasons such as: to provide bounding conservative results, to provide a margin of safety for testing, or to facilitate simulant sub-sampling and analysis. In a parallel effort, aqueous electrolyte thermodynamic models were developed for the synthetic Salado and Castile brines to inform the experimental design, facilitate laboratory data interpretation, and allow extension of evaluations beyond the parameters tested. Thermodynamic modeling simulations including the MgO and Fe 0 additives that are directly relevant to the experimental measurements (e.g., pH calibration curve, ORP corrections) are included in this report. The measured fa of the simulants was close to the OLI model predictions for both brines and was largely controlled by the background H 2 partial pressure in the vapor phase as well as H 2 generated in situ in the aqueous phase by the Fe 0 corrosion. The H 2 gas-phase concentration tested and thermodynamically evaluated was much lower than is assumed in the WIPP PA; however, H 2 (g) concentrations significantly below this level are still predicted to result in very reducing conditions. In conclusion: • The experimental results are consistent with thermodynamic model predictions for fa, pH, and the effects of high ionic strength. • Evidence to date suggests that the H2 concentration in the glovebox atmosphere ultimately determined the final E h values of the simulants and resulted in highly reducing conditions. As a result, little difference was observed between the control simulants containing only MgO and the WIPP base-case simulants that contained MgO and Fe 0 . • This test methodology is recommended for future studies evaluating WIPP repository conditions. The methodology includes: (1) background H 2 in argon with agitation ( or could alternatively include in-situ-generated H 2 in sealed bottles); (2) carefully measured and corrected ORP data ( with much effort focused on allowing the probes to fully stabilize); and (3) ionic-strength-corrected pH data. Other best practices, such as simulant sparging/handling, ORP probe replacement, etc., should also be considered. • The coupling of experimental studies and thermodynamic modeling is also highly recommended because these methods inform and direct one another leading to greater confidence in and understanding of the results.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Integrated Disposal Facility FY2011 Glass Testing Summary Report [Erratum]

Pacific Northwest National Laboratory was contracted by Washington River Protection Solutions, LLC to provide the technical basis for estimating radionuclide release from the engineered portion of the disposal facility (e.g., source term). Vitrifying the low-activity waste at Hanford is expected to generate over 1.6 x 10 5 m 3 of glass (Certa and Wells 2010). The volume of immobilized low-activity waste (ILAW) at Hanford is the largest in the DOE complex and is one of the largest inventories (approximately 8.9 x 10 14 Bq total activity) of long-lived radionuclides, principally 99 Tc (t 1/2 = 2.1 x 10 5 ), planned for disposal in a low-level waste (LLW) facility. Before the ILAW can be disposed, DOE must conduct a performance assessment (PA) for the Integrated Disposal Facility (IDF) that describes the long-term impacts of the disposal facility on public health and environmental resources. As part of the ILAW glass testing program PNNL is implementing a strategy, consisting of experimentation and modeling, in order to provide the technical basis for estimating radionuclide release from the glass waste form in support of future IDF PAs. The purpose of this report is to summarize the progress made in fiscal year (FY) 2011 toward implementing the strategy with the goal of developing an understanding of the long-term corrosion behavior of low-activity waste glasses.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Solubility controls on plutonium and americium release in subsurface environments exposed to acidic processing wastes

To identify the role of waste composition and sediment interactions in controlling Pu and Am mobility in contaminated sediments at Hanford, a legacy nuclear site, Pu and Am concentrations in solutions equilibrated with contaminated sediments from beneath the 216-Z-9 (Z-9) Trench were compared to the solubilities of PuO 2 materials, synthesized by methods representative of the disposed wastes, in the absence of sediments. Furthermore, this work shows that the solubilities of PuO 2 materials synthesized by different methods, and with varying particle sizes, agree with PuO 2 (am,hyd), although dissolution kinetics differed between materials. According to saturation index (SI) calculations, PuO 2 (am,hyd) is likely also controlling Pu release from sediments under conditions where phosphate concentrations are low. However, both Pu-phosphate and Am-phosphate phases, identified in SI calculations and by high resolution transmission electron microscopy, play roles in controlling release in low pH, high phosphate, shallow sediments just below the Z-9 Trench. The elevated phosphate is likely due to decomposition of tributyl phosphate from waste solutions over time. Sediments from deeper in the subsurface beneath the Z-9 Trench are less acidic and contain less phosphate, with Pu solubility likely controlled by PuO 2 (am,hyd) that precipitated following neutralization of the acidic waste stream. Controls on Am concentrations in deeper sediments are more complex and potentially involve sediment adsorption and/or release from Pu (1-x) Am x O 2 following Am in-growth. The concentrations of both Pu and Am were elevated in the colloidal fraction associated with shallow sediments, but not in PuO 2 experiments, suggesting the presence of Pu/Am pseudocolloids (e.g., Pu/Am associated with mineral colloids). However, Pu and Am association with the colloidal size fraction was not observed in deeper sediments, suggesting transport of Pu and Am to these depths beneath the Z-9 Trench was not due to colloidal transport.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2025

Reactive Transport Modeling of Hydrogen Production from Serpentinization of Olivine

Hydrogen production from serpentinization of ultramafic rocks represents a promising natural pathway for generating carbon-free energy, yet its kinetics and controlling factors remain incompletely understood. A key challenge in advancing serpentinization research lies in the heterogeneity of porosity and permeability in rocks, which leads to nonuniform fluid velocity fields, as well as uncertainties in estimating reactive surface area and identifying appropriate mineral reaction equilibria. Additional complexities arise from the role of dissolved SiO 2 , Fe 2+ /Fe 3+ partitioning, and the limited effect of pH variations within the strongly alkaline regime on hydrogen yields. These challenges hinder straightforward extrapolation from laboratory tests to practical applications of hydrogen production from natural rocks. Here, in this work, we address these questions using a simulation-based reactive transport modeling framework calibrated against controlled laboratory experiments reported elsewhere. The model couples geochemical kinetics, multiphase flow, and mineralogical feedbacks, enabling systematic evaluation of how surface area, dissolved silica concentration, Fe redox state, temperature, and pressure govern serpentinization and H2 generation. We find that surface area exerts the strongest control on reaction rates and hydrogen yields, while Fe 2+ /Fe 3+ ratios act as secondary modulators. Elevated dissolved silica concentrations suppress hydrogen production but accelerate serpentine precipitation, whereas increasing pH beyond 12 within the strongly alkaline regime produces only marginal gains. Finally, we demonstrate that integrating targeted experiments with calibrated simulations offers a powerful and efficient approach for predicting hydrogen yields and assessing parameter trade-offs in industrial-scale applications. This integration can substantially reduce the experimental burden while improving predictive capability, thereby enhancing both the mechanistic understanding and the practical feasibility of hydrogen production from serpentinization.

08 HYDROGEN↗
DOE OSTI2021

Thin Water Films Enable Low-Temperature Magnesite Growth Under Conditions Relevant to Geologic Carbon Sequestration

Injecting supercritical CO2 (scCO2) into basalt formations for long-term storage is a promising strategy for mitigating CO2 emissions. Mineral carbonation can result in permanent entrapment of CO2; however, carbonation kinetics in thin H2O films in humidified scCO2 is not well understood. We investigated forsterite (Mg2SiO4) carbonation to magnesite (MgCO3) via amorphous magnesium carbonate (AMC; MgCO3·xH2O, 0.5 < x < 1), with the goal to establish the fundamental controls on magnesite growth rates at low H2O activity and temperature. Experiments were conducted at 25, 40, and 50 °C in 90 bar CO2 with a H2O film thickness on forsterite that averaged 1.78 ± 0.05 monolayers. In situ infrared spectroscopy was used to monitor forsterite dissolution and the growth of AMC, magnesite, and amorphous SiO2 as a function of time. Geochemical kinetic modeling showed that magnesite was supersaturated by two to three orders of magnitude and grew according to a zero-order rate law. The results indicate that the main drivers for magnesite growth are sustained high supersaturation coupled with low H2O activity, a combination of thermodynamic conditions not attainable in bulk aqueous solution. This improved understanding of reaction kinetics can inform subsurface reactive transport models for better predictions of CO2 fate and transport.

Kerisit, Sebastien N.↗
DOE OSTI2025

Geo-Responsive Chemomechanics in Aluminum Oxyhydroxide via Alkali-Driven Dehydroxylation for Supercritical Geothermal Systems

Widespread use of enhanced geothermal systems can revolutionize global renewable electrical power access, yet its advancement is hindered by the inherent instability of Portland cement-based chemistries for geothermal well construction under high-temperature corrosive conditions. Here, in this work, we demonstrate the tunable mechanical performance of aluminum oxyhydroxides as cementitious materials through an alkali-controlled dehydroxylation reaction pathway for long-term applications under supercritical geothermal environments. Notably, the synthesized aluminum oxyhydroxides demonstrate remarkable stability, maintaining superior mechanical performance under supercritical conditions for over 30 days. Synchrotron X-ray diffraction, spectroscopy measurements and geochemical thermodynamic modeling uncover that the gibbsite dehydroxylation pathway functions as a key dial for tuning the chemomechanics, rendering the aluminum oxyhydroxide a strong cementitious material. By uncovering the mechanistic role of alkali-driven dehydroxylation, this work proposes a cementitious chemistry distinct from conventional Portland cement and geopolymer-dominated alkali-activated systems, laying the groundwork for developing next-generation cementitious materials for supercritical geothermal energy exploitation.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

A Geo‐Structurally Based Correction Factor for Apparent Dissolution Rates in Fractured Media

Abstract Field measurements of apparent geochemical weathering reaction rates in subsurface fractured porous media are known to deviate from laboratory measurements by multiple orders of magnitude. To date, there is no geologically based explanation for this discrepancy that can be used to predict reaction rates in field systems. Proposed correction factors are typically based on ad hoc characterizations related to geochemical kinetic models. Through a series of high‐fidelity reactive transport simulations of mineral dissolution within explicit 3D discrete fracture networks, we are able to link the geo‐structural attributes with reactive transport observations. We develop a correction factor to linear transition state theory for the prediction of the apparent dissolution rate based on measurable geological properties. The modified rate law shows excellent agreement with numerical simulations, indicating that geological structure could be a primary reason for the discrepancy between laboratory and field observations of apparent dissolution rates in fractured media.

58 GEOSCIENCES↗
DOE OSTI2024

Weathering Incongruence in Mountainous Mediterranean Climates Recorded by Stream Lithium Isotope Ratios

Lithium isotope ratios (δ 7 Li) of rivers are increasingly serving as a diagnostic of the balance between chemical and physical weathering contributions to overall landscape denudation rates. Here, we show that intermediate weathering intensities and highly enriched stream δ 7 Li values typically associated with lowland floodplains can also describe small upland watersheds subject to cool, wet climates. This behavior is revealed by stream δ 7 Li between +22.4 and +23.5‰ within a Critical Zone observatory located in the Cévennes region of southern France, where dilute stream solute concentrations and significant atmospheric deposition otherwise mask evidence of incongruence. The water-rock reaction pathways underlying this behavior are quantified through a multicomponent, isotope-enabled reactive transport model. Using geochemical characterization of soil profiles, bedrock, and long-term stream samples as constraints, we evolve the simulation from an initially unweathered granite to a steady state weathering profile which reflects the balance between (a) fluid infiltration and drainage and (b) bedrock uplift and soil erosion. Enriched stream δ 7 Li occurs because Li is strongly incorporated into actively precipitating secondary clay phases beyond what prior laboratory experiments have suggested. Chemical weathering incongruence is maintained despite relatively slow reaction rates and moderate clay accumulation due to a combination of two factors. First, reactive primary mineral phases persist across the weathering profile and effectively “shield” the secondary clays from resolubilization due to their greater solubility. Second, the clays accumulating in the near-surface profile are relatively mature weathering byproducts. These factors promote characteristically low total dissolved solute export from the catchment despite significant input of exogenous dust.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1984

An integrated dynamical and geochemical approach to lunar origin modelling

The three major categories of models of lunar origin which explain the Moon's properties are complete and a more general scenario is presented. The model presented is as the Earth grew by planetesimal bombardment, a circumterrestrial cloud of particles was created from a combination of impact ejected mantle material and planetismals captured directly into orbit around the Earth. The compositional properties are explained two ways: (1) a few big late planetismals of diverse composition are captured in orbit and/or hit the Earth; and (2) the circumterrestrial swarm acts as a filter, preferentially capturing small weak silicate bodies, while passing large iron planetismals.

Greenberg, R.↗
NASA NTRS1985

Seafloor hydrothermal activity and spreading rates - The Eocene carbon dioxide greenhouse revisited

A suggestion has been made that enhanced rates of hydrothermal activity during the Eocene could have caused a global warming by adding calcium to the ocean and pumping CO2 into the atmosphere (Owen and Rea, 1984). This phenomenon was purported to be consistent with the predictions of the CO2 geochemical cycle model of Berner, Lasaga and Garrels (1983) (henceforth BLAG). In fact, however, the BLAG model predicts only a weak connection between hydrothermal activity and atmospheric CO2 levels. By contrast, it predicts a strong correlation between seafloor spreading rates and pCO2, since the release rate of CO2 from carbonate metamorphism is assumed to be proportional to the mean spreading rate. The Eocene warming can be conveniently explained if the BLAG model is extended by assuming that the rate of carbonate metamorphism is also proportional to the total length of the midocean ridges from which the spreading originates.

Kasting, J. F.↗
NASA NTRS1989

Iridium in the Vredefort Bronzite Granophyre - Impact melting and limits on a possible extraterrestrial component

The Bronzite Granophyre was analyzed for Ir and other elements to examine whether its origin was due to impact melting or conventional igneous assimilation and to search for a possible extraterrestrial component in the melt. No such component is identified. Ir values range from 50-130 pg/g, equivalent to less than 0.05 percent C1 material. Further evidence against an extraterrestrial source of siderophiles are the low Ir/Au ratios and the absence of systematic correlations between Ir and other elements. However, the Granophyre is significantly enriched in Ir over other Vredefort granitic rocks, implying a fundamental difference between them. Geochemical mixing models indicate that the Ir observed in the Granophyre can be supplied from high-Ir local shales and quartzites during an impact melting event. The results are consistent with the view that the Granophyre is an impact melt and that Vredefort is an impact structure.

French, B. M.↗
NASA NTRS1993

Characterization of Minnesota lunar simulant for plant growth

Processing of lunar regolith into a plant growth medium is crucial in the development of a regenerative life support system for a lunar base. Plants, which are the core of such a system, produce food and oxygen for humans and, at the same time, consume carbon dioxide. Because of the scarcity of lunar regolith, simulants must be used to infer its properties and to develop procedures for weathering and chemical analyses. The Minnesota Lunar Simulant (MLS) has been identified to date as the best available simulant for lunar regolith. Results of the dissolution studies reveal that appropriately fertilized MLS can be a suitable medium for plant growth. The techniques used in conducting these studies can be extended to investigate the suitability of actual lunar regolith as a plant growth medium. Dissolution experiments were conducted using the MLS to determine its nutritional and toxicity characteristics for plant growth and to develop weathering and chemical analysis techniques. Two weathering regimes, one with water and one with dilute organic acids simulating the root rhizosphere microenvironment, were investigated. Elemental concentrations were measured using inductively-coupled-plasma (ICP) emission spectrometry and ion chromatography (IC). The geochemical speciation model, MINTEQA2, was used to determine the major solution species and the minerals controlling them. Acidification was found to be a useful method for increasing cation concentrations to meaningful levels. Initial results indicate that MLS weathers to give neutral to slightly basic solutions which contain acceptable amounts of the essential elements required for plant nutrition (i.e., potassium, calcium, magnesium, sulfur, zinc, sodium, silicon, manganese, copper, chlorine, boron, molybdenum, and cobalt). Elements that need to be supplemented include carbon, nitrogen, and perhaps phosphorus and iron. Trace metals in solution were present at nontoxic levels.

Oglesby, James P.↗
NASA NTRS1994

Lunar and martian meteorite delivery services

Launch mechanisms for lunar and martian meteorites have been investigated, by integrating physical modeling constraints, geochemical cosmic-ray exposure (CRE) constraints, and petrologic constraints. The potential source region for lunar meteorites is remarkably small compared to the final crater volume. CRE constraints indicate that most launches start at depths of less than or equal to 3.2 m, and cratering theory implies derivation of suitably accelerated objects from a subvolume with diameter only about 0.3 x the final crater diameter. The shallow depth provenance is probably related to shock-wave interference, enhanced by the lunar regolith's extremely low compressional wave velocity. CRE constraints alone imply that four to five separate launch events are represented among the eight well-studied lunar meteorites. Most of the lunar meteorites are regolith breccias, which tend to show only limited compositional diversity within any kilometer-scale region of the Moon. Several others are polymict breccias, which also show relatively subdued compositional diversity, compared to igneous rocks. The observed diversity among these samples in terms of abundances of mare basalt and KREEP, and in Mg/(Mg + Fe) ratio, implies that among eight well-studied lunar meteorites only two potential source craters pairings are plausible: between Asuka-881757 + Y-793169 (most probable) and between Y-793274 + EET875721. Altogether, these eight lunar meteorites apparently represent at least six separate source craters, including three in the past 10(exp 5) years and five in the past 10(exp 6) years. CRE constraints imply that SNC meteorites are launched from systematically greater than lunar meteorites. SNCs are also systematically bigger, and all nine well-studied SNCs are uncommonly young (by martian standards) mafic igneous rocks. Comparison between Viking and Apollo results reveals that rocks the size of common meteorites are remarkably scarce in the martian regolith, probably due to pervasive weathering. A plausible explanation for these trends is that most old, small, and shallow rocks on Mars have been weakened by a two-stage process of brecciation followed by pervasive weathering, to the point where they seldom survive the stresses of spallation off the planet. The scarcity of source-crater pairing among the lunar meteorites implies that these objects can be launched from craters much smaller than previously estimated and tends to support the suggestion of Rabinowitz (1993) that present-day cratering rates for the Earth-Moon region may be higher than previously estimated.

Warren, Paul H.↗
NASA NTRS1985

Seafloor hydrothermal activity and spreading rates: the Eocene carbon dioxide greenhouse revisted

A suggestion has been made that enhanced rates of hydrothermal activity during the Eocene could have caused a global warming by adding calcium to the ocean and pumping CO2 into the atmosphere (Owen and Rea, 1984). This phenomenon was purported to be consistent with the predictions of the CO2 geochemical cycle model of Berner, Lasaga and Garrels (1983) (henceforth BLAG). In fact, however, the BLAG model predicts only a weak connection between hydrothermal activity and atmospheric CO2 levels. By contrast, it predicts a strong correlation between seafloor spreading rates and pCO2, since the release rate of CO2 from carbonate metamorphism is assumed to be proportional to the mean spreading rate. The Ecocene warming can be conveniently explained if the BLAG model is extended by assuming that the rate of carbonate metamorphism is also proportional to the total length of the midocean ridges from which the spreading originates.

Seawater/chemistry↗