May 07 2026

The role of reactive transport models in enhanced weathering today

by Tyler Kukla +Freya Chay

The Modeling Bytes series examines decisions that matter when using models to estimate enhanced weathering outcomes. Most posts in this series are quite technical, so to make the series more accessible we’re starting with a brief introduction to reactive transport models (RTMs) — the family of models we focus on throughout. This post covers what RTMs are, and the role we think they should play in enhanced weathering (EW) today.

So, what are RTMs?

True to their name, reactive transport models simulate chemical reactions and fluid transport. These coupled processes underpin the many systems where RTMs are applied, including: soil formation, chemical weathering, geothermal reservoirs, geological carbon sequestration, bioremediation, industrial waste management, contaminant transport, and groundwater quality. RTMs have been used in these contexts for decades, but we’ve only recently started applying them to EW.

Beyond the basic requirements that RTMs include reaction and transport processes, the best way to understand these models is not by trying to describe what they are, but by unpacking what they can be set up to do. That’s because RTMs are highly flexible. Modelers can make dozens of decisions to use the same model for vastly different applications — across topics, or within EW. Figure 1 shows some of these decisions, grouped by which stage of model setup they relate to: domain, forcings, spinup, reactions, and transport.

Domain
Spatial dimension
Numerical solver
Physical characteristics
+ more
Forcing
Deployment and other interventions
Meteorology
Ambient chemistry
+ more
Spinup
Equilibration timescale
Equilibrium metric
Initial conditions
Target conditions
+ more
Reactions
Species
Geochemical database
Surface complexation and ion exchange
Biochemistry
+ more
Transport
Solute transport
Hydraulic properties and flow regimes
Flow heterogeneity
Gas diffusivity
+ more
figure 1
/

Some examples of consequential RTM set up decisions. Domain refers to the physical space and time conditions the model simulates; forcing refers to the external inputs that the system responds to; spinup refers to how the model represents the initial domain conditions; and reactions and transport refer to how the model represents chemistry and fluid transport, respectively. Click any box for an expanded description.

To date, most EW RTM work has been done with a handful of different models (Table 1), though the list is growing. Critically, the model you choose constrains the decision space laid out above. For example, SMEW is zero-dimensional but the code is easy to modify and customize; CrunchFlow has extensive options for reaction kinetics and uses a very fast and accurate numerical solver; and MIN3P has sophisticated soil hydrology, well-suited for event-scale infiltration and groundwater coupling.

ModelKey EW references
ARTEMIS (with PHREEQC)Taylor et al. preprint
ATS (with PFLOTRAN)Shaheen et al. preprint
CrunchFlowMaher and Rogers, in prep.
Geochemist’s workbenchKhalidy et al. preprint
PHREEQCKelland et al. (2020); Lewis et al. (2021); Vienne et al. (2022)
SCEPTERKanzaki et al. (2022); Kanzaki et al. (2024)
SMEWBertagni et al. (2025); Anand et al. (2026)
T&C SMEWZhang et al. (2025)
TOUGHREACTDeng et al. (2023)
Beerling modelBeerling et al. (2020); Beerling et al. (2025)
Cipolla model
(precursor to SMEW)
Cipolla et al. (2021a); Cipolla et al. (2021b)

In other words, RTMs are less like a single tool you can just pick up and use, and more like a toolbox with many options for tackling different projects. Figuring out how to use them effectively for EW requires an iterative process — applying the models, learning about EW itself, and identifying where modeling approaches need refinement. We are still in the early stages of that process.

Using RTMs to learn

The flexibility inherent to RTMs makes them ideal tools for basic research and learning, which we see as their primary role in EW today. Practitioners can develop and test theories by adding layers of complexity, piece-by-piece, within a single, internally consistent framework. In this way, RTMs can help us build a theoretical foundation for EW that holds across scales — connecting the dissolution of micron-scale dust in soils, to alkalinity export through a catchment, and ultimately the storage of carbon in the ocean.

We don’t have that coherent theoretical framework in EW today, but RTMs are already helping us learn. They’re providing extra context for field data,11M B Bertagni et al. (2025) Advancing enhanced weathering modeling in soils: Critical comparison with experimental data JAMES ,22H Green et al. (2024) Carbon dioxide removal via weathering of sugarcane mill ash under different soil conditions Applied Geochemistry ,33M E Kelland et al. (2020) Increased yield and CO₂ sequestration potential with the C₄ cereal Sorghum bicolor cultivated in basaltic rock dust-amended agricultural soil Global Change Biology ,-44A Vienne et al. (2022) Enhanced weathering using basalt rock powder: Carbon sequestration, co-benefits and risks in a mesocosm study with Solanum tuberosum Frontiers in Climate unpacking how key geochemical interactions play out, 55G Cipolla et al. (2022) Effects of precipitation seasonality, irrigation, vegetation cycle and soil type on enhanced weathering – modeling of cropland case studies across four sites Biogeosciences ,66H Deng et al. (2023) The environmental controls on efficiency of enhanced rock weathering in soils Scientific Reports ,77Y Kanzaki et al. (2025) Soil cation storage is a key control on the carbon removal dynamics of enhanced weathering Environmental Research Letters and helping us understand what it might take to do EW at scale.88S H Baek et al. (2023) Impact of climate on the global capacity for enhanced rock weathering on croplands Earth's Future ,99D J Beerling et al. (2020) Potential for large-scale CO₂ removal via enhanced rock weathering with croplands Nature ,1010D J Beerling et al. (2025) Transforming US agriculture for carbon removal with enhanced weathering Nature But to continue making these learning efforts effective, we also need to better understand the RTMs themselves.

The path forward

If we eventually want to use RTMs to confidently simulate real-world outcomes, there is a lot of work ahead.

First, we need to develop RTMs to capture more of the complexity of EW. Soil biogeochemistry is messy, and EW simulations today make some drastic simplifications. For example, almost all EW RTM studies to date have used spatial domains defined in zero or one dimension — a point or a column of soil — rather than higher dimensions capable of resolving lateral flow paths. They also typically represent just the top ~20-100 cm of the soil column, capturing where alkalinity is generated but not how it’s exported and stored. Studies often make further simplifications, such as fixed rates of baseline alkalinity production1111G Cipolla et al. (2021) The role of hydrology on enhanced weathering for carbon sequestration I. Modeling rock-dissolution reactions coupled to plant, soil moisture, and carbon dynamics Advances in Water Resources or soil respiration,99D J Beerling et al. (2020) Potential for large-scale CO₂ removal via enhanced rock weathering with croplands Nature and tunable parameters that collapse numerous complex processes to a single value.11M B Bertagni et al. (2025) Advancing enhanced weathering modeling in soils: Critical comparison with experimental data JAMES On the one hand, these simplifications are necessary — we learn by adding complexity one layer at a time. But they also show us how far we have to go. Bridging the gap between models and real-world complexity, or even justifying the simplifications, will require significant investment in model development.

Second, we need to scrutinize model decisions and develop field-wide norms to guide how we set up EW simulations. Today, almost all consequential modeling choices rest on individual experts’ judgment. But in more mature fields, they draw from a shared body of work. For example, Earth system model parameterizations for stomatal conductance or atmospheric convection rest on numerous research papers establishing which choices are acceptable for the problem at hand. Such a body of work stress-tests the options and surfaces the tradeoffs and uncertainties that are otherwise implicit, or overlooked. This foundation benefits the whole field: it makes it easier to justify model choices and to contextualize results across studies.

The Modeling Bytes series is one attempt to contribute to this foundation-building effort. As the iterative process of applying, learning, and refining EW models continues, we expect things to change. Mistakes will be discovered and corrected. Model setups will grow more sophisticated and increasingly grounded in emerging field-wide norms. Since this process is only just beginning, we encourage anyone relying on RTM outputs to treat the carbon removal estimates they produce as provisional. We hope this primer, and the posts that follow, give you better tools for evaluating what EW models can and can’t yet tell us.


Credits

Tyler and Freya wrote the article. Kata Martin designed the figure. Grayson Badgley provided feedback on the text.

Please cite as T Kukla and F Chay (2026) “The role of reactive transport models in enhanced weathering today.” CarbonPlan https://carbonplan.org/research/modeling-bytes-01-primer

Terms

CarbonPlan’s work on this article was supported by a grant from the Chan Zuckerberg Initiative DAF, an advised fund of Silicon Valley Community Foundation.

Article text and figures are made available under a CC BY 4.0 International license.


Questions? Interested in collaborating on these problems?
email us
May 07 2026
EMAIL
hello@carbonplan.org
NEWSLETTER
Subscribe
CarbonPlan is a registered nonprofit public benefit corporation in California with 501(c)(3) status.
(c)
CARBONPLAN
TERMS OF USE
/
PRIVACY POLICY
SCROLL: 0.00
aa40475