- 1Department of Geology, Carleton College, Northfield, MN, United States
- 2Department of Earth and Planetary Sciences, Yale University, New Haven, CT, United States
- 3Department of Earth and Planetary Sciences, Massachusetts Institute of Technology, Cambridge, MA, United States
- 4Department of Geology and Geophysics, Woods Hole Oceanographic Institution, Woods Hole, MA, United States
- 5Institute of the Environment, University of California, Davis, Davis, CA, United States
- 6Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, CA, United States
- 7Yale Center for Natural Carbon Capture, Yale University, New Haven, CT, United States
Enhanced weathering is an emergent pathway for permanent atmospheric carbon dioxide removal (CDR). However, despite a dramatic increase in academic and commercial research, there remain relatively few published examples of field evidence demonstrating the effectiveness of enhanced weathering. Here, we present results from a three-year field trial that evaluated steel slag and crushed basalt applied as amendments in a conventional agricultural system in the Midwestern United States. Steel slag applied to initially acidic soil increased porewater pH and alkalinity and increased soil pH and Ca-saturation. Together, changes in porewater chemistry and soil properties provide strong evidence for steel slag weathering and CDR. However, steel slag applied to soils with a neutral initial pH did not generate significant changes in soil or porewater chemistry. In addition, coarse-grained crushed basalt did not generate significant change in any of the soils. Strong acid effects were apparent in all 3 years of monitoring soil porewater chemistry. Overall, our results demonstrate clear evidence of CDR from applying steel slag amendments to acidic cropland soils while also highlighting the difficulty of greenhouse gas reduction accounting from enhanced weathering and the variable outcomes that can occur depending on feedstock and soil type.
1 Introduction
Enhanced weathering is a promising carbon dioxide removal (CDR) strategy with co-benefits for crop yields and soil health that could reach gigaton scale targets within the next decade if widely deployed (Beerling et al., 2025; Goll et al., 2021; Smith et al., 2024). In its simplest form, enhanced rock weathering delivers CDR through the neutralization of carbonic acid, which forms as CO2 from the atmosphere or soil dissolves into soil waters. The reactions below illustrate both carbonic acid formation (Equation 1) and the neutralization of carbonic acid through reaction with an example Ca-silicate (here, larnite – Ca2SiO4; Equation 2).
In this simplified case, silicate minerals buffer and consume carbonic acid, release base cations (here Ca2+ but commonly Mg2+ as well), and convert the dissolved carbon into bicarbonate (HCO3−), increasing the alkalinity of the soil porewaters. The bicarbonate and base cations generated through enhanced rock weathering eventually leach through the soil profile and into river networks to reach the ocean. There, bicarbonate can persist for thousands of years, or eventually precipitate calcium carbonate, leading to effectively permanent CDR (Seifritz, 1990; Hartmann and Kempe, 2008; Köhler et al., 2010; Hartmann et al., 2013; Beerling et al., 2018; Beerling et al., 2024).
In croplands, enhanced rock weathering aims to optimize the practice of soil pH management by using materials and deployment strategies designed to maximize CDR, while providing cost effective solutions to farmers with agronomic co-benefits.
Silicate rocks (e.g., basalt) and minerals (e.g., olivine, wollastonite) often offer a similar set of co-benefits as traditional ag lime (often carbonate rocks), including the release of micronutrients, increased crop yields, and reduced emissions of other greenhouse gases such as nitrous oxide (Beerling et al., 2024; Kantola et al., 2023; Skov et al., 2024; Haque et al., 2019; Edwards et al., 2017; Chiaravalloti et al., 2023). Most commonly, enhanced weathering studies target natural silicate rocks like basalt (Beerling et al., 2018; Kantola et al., 2023; Skov et al., 2024). Industrial products like steel slag and returned concrete have also drawn some attention due to their relatively rapid weathering rates, potential for CDR, and agronomic benefits (Das et al., 2020; Das et al., 2019; McDermott et al., 2024). Steel slag could be particularly attractive given the demonstrated benefits of Si fertilizers on crop resilience (Das et al., 2020). However, there has yet to be a field trial testing steel slag in an enhanced weathering context for CDR.
Despite the potential of enhanced weathering, there are still very limited field studies that have published empirical evidence for rates of rock weathering and CDR. Perhaps the most difficult aspect of enhanced weathering as a CDR approach is the direct measurement of mineral weathering and the accompanying CDR in field deployments (Clarkson et al., 2024; Almaraz et al., 2022). Published rates for CDR from enhanced weathering field trials are highly variable, spanning four orders of magnitude (CarbonPlan, n.d.). Given that commercial carbon markets and widespread adoption hinge on adequate proof that enhanced weathering is delivering CDR at scale, robust demonstration of the quantity of CDR is essential (Smith et al., 2024). In the past several years, work on enhanced weathering has rapidly increased across academic and industry groups, with a growing emphasis on developing field trials and measurement, reporting, and verification (MRV) techniques across different soil types, climate patterns, and crop regimes (Beerling et al., 2024; Kantola et al., 2023; Haque et al., 2019; Larkin et al., 2022; Gunnarsen et al., 2022; Guo et al., 2023; Holzer et al., 2023). Current strategies for measurement in enhanced weathering generally focus on tracking weathering rates through measurement of cation-loss from feedstock material (Reershemius et al., 2023; Suhrhoff et al., 2024), direct sampling and measurement of soil porewater chemistry (Holzer et al., 2023; McBride et al., 2025), and local watershed monitoring (Larkin et al., 2022; Taylor et al., 2021).
The collection and measurement of soil porewaters provides arguably the most direct measure of weathering and initial CDR through analysis of total alkalinity (which approximates the concentration of bicarbonate in most natural waters Larkin et al., 2022; Holzer et al., 2023; Neal, 2001), dissolved inorganic carbon, base cations, and non-carbonate system anions. However, porewaters are notoriously difficult to sample, and cannot easily be measured continuously. Watershed scale studies can provide a large-scale integrated measure of CDR, but signals can be dilute and difficult to measure (Larkin et al., 2022). In addition, weathering products have variable mobility in soil, with products interacting with charged exchange sites on soil particles, or becoming entangled in secondary mineral phases that are difficult to detect (Slessarev et al., 2021). Given some of these complexities, solid-phase cation-loss approaches that leverage soil samples to evaluate the loss of base-cations from applied feedstock material is an emerging proxy (Reershemius et al., 2023). However, mass-balance based approaches are a measure of the rate of chemical weathering and base cation release that establishes a potential rate of CDR (McBride et al., 2025). When converting weathering rates to CDR, complications around the transport time of released base cations and potential impacts of strong acid weathering are areas of active research and debate (Kanzaki et al., 2025; Power et al., 2025).
Limited attention has been paid to any associated positive or negative impacts of crushed rock amendments on soil organic matter (SOM), Earth’s largest terrestrial pool of organic C, and a critical source of nitrogen (N) and other nutrients (Scharlemann et al., 2014; Sokol et al., 2024). Given that positive effects on SOM may enhance CDR and boost soil fertility, while negative effects could entirely negate any benefit of inorganic CDR this is an essential area for continued work. The limited available evidence shows contrasting effects of crushed rock on SOM pools in different ecosystems, as well as variable effects on different SOM pools including the slower-cycling mineral-associated soil organic matter (MAOM) and faster-cycling particulate organic matter (POM) (Sokol et al., 2024; Xu et al., 2024; Sohng et al., 2025; Buss et al., 2024). No studies to date have investigated effects on SOM pools in a Midwestern field trial, nor the effects of steel slag on SOM, MAOM, and POM. Put together, there is a clear need for more research and data to inform our collective understanding of enhanced weathering as emerging CDR pathways.
Here, we report results from a three-year field trial in Southeastern Minnesota designed to evaluate the effectiveness of enhanced weathering as a strategy for CDR. We applied crushed basalt and steel slag to a strip-till corn and soybean rotational field typical of a conservation-oriented Midwestern farm. Porewater chemistry, soil properties, and gas fluxes were monitored for three field seasons. We evaluated whether there was direct evidence for mineral weathering and CDR. In addition, we evaluated the role of strong acid weathering and considered the impacts of strong acids on CDR quantification. Finally, we investigated how weathering of crushed basalt and steel slag affected soil organic carbon cycling, through monitoring soil CO2 gas-flux, total SOM, and different soil organic carbon pools (MAOM and POM).
2 Site details and experimental methods
2.1 Field site
The field trial is situated on an eight-hectare agricultural field on the campus of Carleton College in Northfield, Minnesota (Figure 1). Since the 1930s, the field has been in near-constant cultivation, with corn and soybeans rotating as the primary crops for the past several decades. In the mid 2010s the field transitioned from conventional tillage to strip-tillage. Crop rotation during the field trial included soybeans in 2022 and 2024 and corn in 2023. Aerial images of the site show that there has been land-management variation across the site since at least 2003, with the northern portion of the field remaining out of cultivation in a subset of years (see Supplementary Figure S1). In addition, dredged sediments from an artificial lake on campus were spread across the northern portion of the field during the early 2000s (estimated extent of dredging covers the full extent of Block B and all of Block A except the easternmost slag plot). These variations in land use correspond to a gradient in soil properties across the field area (discussed below in Section 3.1), making it ideal to explore the effects of varying soil types on enhanced weathering.
Figure 1. Field map and overview of study area. Location is Northfield Minnesota (inset map in upper left), with the study field (outlined in yellow at left) on campus at Carleton College. Detailed plot layout shown at right with basalt, slag, and control treatments distributed across the field in four experimental blocks.
Soils in the region are dominated by loam developed atop Quaternary glacial sediments or Ordovician sedimentary units (limestones, dolostones, and sandstones). In this field, soils are loam and clay loam developed atop Quaternary loess and fine loamy till mapped primarily as the Hayden loam and Blooming silt loam in the field area (Soil Survey Staff United States Department of Agriculture, n.d.). The Hayden loam is a mesic, superactive, Glossic Hapludalf and The Blooming silt loam is a fine-loamy, mixed, superactive mesic Mollic Hapludalf (Soil Survey Staff United States Department of Agriculture, n.d.).
Climate in the area is temperate, with a mean annual temperature (MAT) of 7.9 °C and mean annual precipitation (MAP) of 902 mm yr.−1 (Supplementary Table S1; NOAA). Growing season climate in the area (May-Sept) has average temperature of 18.5 °C and rainfall of 574 mm (1990–2020 normals for Rice County, NOAA). Conditions during the 3 years of this experiment were generally warmer than normal in 2023 and 2024 and drier than normal in 2022 and 2023 (see Supplementary Table S1).
2.2 Experimental design, baseline sampling, and materials
We used a randomized block design to distribute treatment and control plots (~0.5 acres each) across the study area to minimize impacts of soil variation on our interpretation of treatment effects (Figure 1). A set of baseline soil samples were collected in June 2021 to constrain soil properties across the field area prior to spreading of any feedstock. Baseline soil was collected at four depth intervals: 0–10 cm, 10–20 cm, 20–40 cm, and 40–60 cm. At each depth interval, a total of 12 individual soil cores were collected at semi-random locations across the plot area (~0.5 acres) and homogenized into one sample per plot (see Supplementary Figure S2). All homogenized baseline samples were air dried and sieved to remove any particles > 2 mm prior to analysis.
Each baseline soil sample was analyzed for soil pH, total soil organic matter (SOM) and SOM pools (MAOM and POM; see methods for this below in section 2.3), carbonate minerals, cation exchange capacity (CEC), and 27 element microwave digestion ICP-OES. Soil pH was measured using a Mettler Toledo Seven-Multi pH meter and InLab Routine Pro electrode on a 1:1 soil to water mixture of 10 g dried soil and 10 mL deionized water with a 15-min equilibration period. Weight percentage organic matter (SOM) and carbonate minerals were determined through loss-on-ignition (LOI) at temperatures of 360 °C and 1,000 °C respectively, with a 2-h heating period at each temperature and an intermediate heating at 105 °C to remove any residual water. For CEC and the 27 element ICP-OES, baseline samples were sent to the University of Minnesota Research Analytical Laboratory (RAL) in St. Paul, Minnesota, USA. Complete analytical methods for those analyses are available on the RAL website1 but are briefly summarized here. CEC was calculated by summation of the exchangeable cations Ca, Mg, Na, K, and Al through an ammonium acetate extraction and subsequent evaluation of leachate using an ICP-OES. The 27-element analysis involved a microwave digestion in concentrated nitric acid following EPA procedure 3,051 and measurement using an ICP-OES. The data for the 27 element ICP-OES analysis is not shown in a figure but is available in the Supplementary data.
Feedstock used for treatments in this study was basalt and steel slag. The basalt is a commercially available quarry sand product derived from an aggregate quarry in the Chengwatana Volcanics of the Keweenawan Supergroup, a sequence of lavas, tuffs, and interflow sedimentary rocks associated with the Neoproterozic (~1.1 Ga) failed rifting of the midcontinent (Wirth et al., 1997). The material is dominantly coarse to medium grained sand (~63% passing 0.5 mm sieve) with 20% fine sand (< 0.25 mm) or finer, placing this on the coarse side of feedstocks used for EW trials. Geochemically, the basalt is 15% Ca + Mg and the mineralogy is dominated by plagioclase, diopside/augite, chlorite, and epidote. The presence of epidote and chlorite is related to low-grade metamorphism associated with hydrothermal fluids (Wirth et al., 1997). The steel slag used in this study is a Ca-rich byproduct of steel production. The material is ~50% CaO + MgO with mineralogy dominated by Ca and Mg rich silicate, (hydr)oxides, and amorphous phases. Grain size is finer in comparison to the basalt, with most of the steel slag gradation being medium sand (~65% < 0.6 mm) and ~25% being very fine sand to silt (~27.3% < 0.15 mm). Both the basalt and steel slag were byproducts of industrial processes. Chemical, mineralogical, and grain size data for the feedstocks is available in the Supplementary file.
Spreading of feedstock materials was completed in November 2021. Four treatment plots were amended with basalt at a target rate of 10 tonnes acre−1 and four plots with slag at a target rate of 2 tonnes acre−1 using commercial ag lime spreading equipment. To validate application rates, six pie pans were placed in each plot at random before spreading (for pan test results and average spreading rates see Supplementary Table S2). The lower application rate for steel slag reflects the much higher liming capacity of steel slag relative to basalt, primarily due to the abundance of hydroxide mineral phases in the steel slag. Care should be taken to not over apply steel slag purely for CDR, as application rates should be tailored to properly adjust soil pH depending on local conditions.
2.3 Soil sampling and analysis years 1–3
Soil samples were collected each year of the three-year field trial. Various sampling strategies were adopted in later years to constrain effects of the study treatment and deal with background soil variability. Below, sampling strategies are described for each sampling year: 2022, 2023, and 2024. While all soil data is included in the Supplementary files, we only discuss soil data in detail from the end of the experiment—sampling year 3 (2024). Below, we briefly describe both sampling strategies and conclusions from evaluating soil in the intermediate years of this study.
In October of 2022 (year 1), soil was collected at 0–10 cm and 10–20 cm depth increments for all treatment plots using the same sampling strategy as the baseline soils. These soils were evaluated for soil pH, organic matter, CEC, and the 27 element ICP-OES, all measured at the RAL facility at the University of Minnesota using methods described above.
In 2023 (year 2), soil was collected during October 2023 at a higher resolution in the southern blocks (C and D) of the field area to determine if higher resolution sampling would reveal a clearer picture of treatment effects. Within an experimental plot, 15 sampling points representing a circular area with a 3 m radius (~28.3 m2 area) were established (see Supplementary Figure S2). At each sampling point, fifteen 0–10 cm soil cores were collected and homogenized into a single sampling point sample. Soil pH and exchangeable base cations (Ca and Mg) for year 2 soils significantly increased in both properties in the slag treatment relative to both the control and basalt treatments (see Supplementary Figure S3). The basalt treatment was not statistically different from control for either property (Supplementary Figure S3). Given the improved resolution using higher density sampling, this approach was adopted for evaluating soil properties in year 3 of this project.
In year 3 (2024), all plots were sampled from 0 to 10 cm at 15 sampling points per plot (following scheme outlined above and in Supplementary Figure S2). Year 3 soils were evaluated for soil organic matter (SOM; weight percentage; based on loss-on-ignition at 360 °C), soil pH (1:1 slurry), a series of extractable elements using the Mehlich 3 extraction including P, K, B, S, Fe, Mn, Cu, Zn, and exchangeable base cations Ca, Mg, K, and Na. Cation exchange capacity (CEC) and base saturation is reported for each sample. Analytical work on all soils from year 3 of this study was conducted by Waypoint Analytical, LLC (Champaign, Illinois, USA).
For baseline soils (2021), year 1 soils (2022), and year 3 soils (2024) we also evaluated effects on different SOM pools: mineral-associated organic matter (MAOM) and particulate organic matter (POM). We separated MAOM and POM pools using a standard physical fractionation method (Sokol et al., 2017; Bradford et al., 2008). First, 5-g of air-dried soil was weighed into a 50-mL falcon tube, and 30-mL of 5% sodium hexametaphosphate (NaHMP) was added to disperse aggregates. Samples were then vortexed and placed on a shaker table for 16 h at 200 rpm. After shaking, the soil solution was passed through a 53-μm sieve. The < 53-μm fraction (clay + fine silt fraction; defined as MAOM) and the > 53-μm fraction (sand, coarse silt, and large organic matter fragments; defined as POM) were oven-dried to constant mass at 105 °C. Dried samples were finely ground to a powder-like consistency with a mortar and pestle and analyzed on an elemental analyzer (Lawrence Livermore National Laboratory, Livermore, CA).
We determined the trace and major element composition of the mineral-bound phases to constrain cation loss from the soil surface (e.g., a soil mass balance or TiCAT methodology Reershemius et al., 2023). This approach uses a mass-balanced based mixing model that constrains mineral weathering rates in enhanced weathering feedstocks. As minerals in applied feedstock are weathered, there is a preferential release of mobile base cations relative to immobile elements like titanium. It can be feasible to disentangle the loss of mobile base cations from a soil-feedstock mixture that has weathered relative to initial baseline measurements (Reershemius et al., 2023). All soil and feedstock samples evaluated using this approach were analyzed for isotope dilution inductively coupled mass spectrometry (ID-ICP-MS) at Yale University following methods in Reershemius et al. (2023).
2.4 Porewater chemistry
In each growing season of the study, lysimeters were used to collect shallow (0–20 cm) soil porewater. In 2022, porewater monitoring occurred in each of the 12 experimental plots. Porewater in the northern experimental plots with high initial soil pH (> 7; all of Block B and Block A basalt and control plots) had considerably elevated porewater alkalinity compared to soils with more acidic baseline soil pH, likely due to soil carbonates dissolving (see Supplementary Figure S4). Due to difficulties in acquiring sufficient water for analysis, all additional porewater monitoring was conducted on experimental plots with acidic initial soil pH, which includes all the southern plots in blocks C and D along with the slag-treatment plot in the easternmost portion of Block A (see Figure 1 for plot and block layouts). In addition, a set of lysimeters was placed in the buffer strip area between block C and D to serve as additional porewater control given the lack of replicates for control in acidic soil areas.
Porewater was collected weekly to bi-weekly during a 10-week monitoring period between mid-June and mid-August 2022, 2023, and 2024. In each plot, a minimum of 3 suction cup lysimeters were installed at depths of ~15 cm using a soil slurry. A plastic apron was installed around the surface of the lysimeter soil connection to avoid any direct flow paths of water adjacent to the lysimeters. Following installation, the first round of water collected was discarded to avoid any contamination caused by the installation soil slurry. Dry periods during each collection season created gaps in data due to lack of available water. Collected water was combined for all lysimeters within a single experimental plot and refrigerated immediately after collection. Prior to analysis, all water was filtered using a 0.45 μm filter. For elemental analysis, water samples were spiked to 1% nitric acid for preservation shortly after collection.
Total alkalinity and pH were measured on all samples with sufficient sample volume using a Hach AT1000 automated titration system with hydrochloric acid. Units for total alkalinity were initially reported as mg CaCO3 L−1. Reported units here are in milliequivalents (meq) L−1. We make the conservative assumption that for samples with pH < 8.3, the total alkalinity is equivalent to the bicarbonate alkalinity (Holzer et al., 2023; Neal, 2001). In 2023, to mitigate limited sample volumes, we conducted alkalinity titrations on diluted samples. Based on regular measurements of alkalinity in reverse osmosis (RO) water used in dilutions, we calculated the samples total alkalinity assuming linear mixing. Because RO waters used in dilutions had relatively low pH (~6), there likely was a minor consumption of alkalinity in samples that were diluted—however this does not appear to be a large effect. For samples that were diluted, we often were able to measure the same sample 2–3 times and use an average value for reporting. When not diluting it was often not possible to obtain duplicate or triplicate measurements of the same sample because of minimum sample volumes (~50 mL) when using pure samples. In comparing alkalinity data, we evaluate measurements throughout the season rather than overanalyzing any one collection interval in part due to these complications.
Major anions, cations, and metals were evaluated throughout the three seasons of monitoring. In all three summers, all water samples were evaluated for a suite of 27 elements using an ICP-OES at the RAL at the University of Minnesota. Here, we report data from this analysis for Ca and Mg, but all elemental data is included within our Supplementary dataset. Notably, across all treatments and monitoring years, there are rarely detectable metals (Cr, Pb, Ni, etc) in soil porewaters (see Supplementary datafile). In 2022, nitrate was measured using colorimetric methods at the RAL at the University of Minnesota. In 2023 and 2024, a set of major anions, including nitrate, sulfate, phosphate, and chloride were measured using an Aquion Ion Chromatography system from ThermoFischer Scientific.
In 2024, water samples were evaluated for total carbon (TC), total organic carbon (TOC), and total inorganic carbon (TIC) at the RAL at the University of Minnesota as an independent check on dissolved carbon separate from alkalinity. Samples were first measured for TC in an Elementar TOC Select Combustion Analyzer. Following an initial TC measurement, an aliquot of the same sample was acidified to remove inorganic carbon and remeasured to quantify TOC. TIC was quantified by subtraction, where TIC = TC-TOC.
2.5 Soil CO2 gas flux
In each of the first two growing seasons (2022 and 2023), soil gas flux was measured using a combination of instruments, chamber types, and soil collar distributions. In 2022, we measured CO2 gas flux from soils across the study area using an EGM-5 from PP Systems (Amesbury, Massachusetts, USA). Measurements were conducted weekly to bi-weekly at multiple locations throughout each plot. For each measurement, a survey chamber was sealed gently to either the bare soil or to a pre-installed soil collar (PVC pipe) and CO2 concentration was monitored in the closed chamber for 90 seconds. The flux of CO2 was then calculated using the change in temperature, air temperature, pressure, and volume of the chamber using standard methods (PP Systems, 2018). In 2023, we measured CO2 gas flux from soil using both the EGM-5 system and a Gasmet GT5000 Terra multigas FTIR gas analyzer (Gasmet Technologies Inc.; Vantaa, Finland). Functionally, both instruments measure change in CO2 concentration in a closed chamber over time. In 2023, all measurements were made on pre-installed soil collars at standard locations across each plot.
3 Results
3.1 Soil properties
Baseline soil properties across the different blocks were variable (see Supplementary Figures S5, S6). Blocks A and B in the northern portion of the field trial had higher soil pH, organic matter, carbonate mineral content, and CEC compared with the southern Blocks C and D (Supplementary Figure S5). One exception to this pattern was the slag treatment plot on the eastern side of Block A (Figure 1) which had soil properties that more closely aligned with Blocks C and D (Supplementary Figure S5).
Comparison between soil from baseline and intermediate sampling years (2022 and 2023) revealed considerable variability (see Supplementary Figures S6, S7 for example with soil pH). Here, we focus on surface soils (0–10 cm) collected at higher sampling density in year 3 (2024, see Supplementary Figure S2). In the more acidic soil in blocks C and D, both soil pH and Ca-saturation were significantly higher for slag treated soils compared against control (Figure 2; p-values < 0.05 for Tukey HSD and ANOVA tests). For soils with more neutral initial pH (blocks A and B), we observed no consistent difference in any soil properties in enhanced weathering plots vs. control plots (Figure 2). For acidic soils in blocks C and D, slag treatments showed shifts in soil pH and Ca-saturation relative to baseline (boxplots compared with grey squares in Figure 2); statistical tests made difficult due to variable sampling strategies.
Figure 2. Soil properties comparing year 3 (colored boxplots) and baseline (2021; grey squares) soil properties (pH, Ca saturation, and organic matter) across treatments and blocks. Each boxplot includes n = 15 data points representing the soil sampling points within each plot. Asterisks above boxplots indicate statistically significant differences between control and slag treatment, where p-values < 0.05 for a Tukey HSD and ANOVA test.
There was no impact of basalt or slag on total SOM (p > 0.05 for Tukey HSD and ANOVA tests; Figure 2). There was also no effect of basalt or slag treatment on the change in MAOM-C (mg C g soil−1; p > 0.05 for Tukey HSD and ANOVA tests) or POM-C (mg C g soil−1; p > 0.05 for Tukey HSD and ANOVA tests) pools relative to baseline for both year 1 (2022) and year 3 (2024) soils (Figure 3).
Figure 3. Change in mineral associated organic matter (ΔMAOM C; upper row) and particulate organic matter (ΔPOM; lower row) for Year 1 (2022; left column) and Year 3 (2024; right column) relative to baseline (2021). Overall, there were no statistically significant differences by treatment for either year (p > 0.05 with a Tukey HSD and ANOVA tests), though ΔMAOM C trended higher in plots A and B in steel slag-amended plots in 2024.
3.2 Porewater chemistry
Variation in soil porewater occurred throughout each growing season, in part due to changes in weather (rainfall and temperature), plant activity, and fertilizer applications. In Supplementary Figures S8–S10, major porewater chemical properties including pH, alkalinity, nitrate, and base cations (Ca and Mg) are shown with time for each growing season.
When taken in aggregate across each growing season, porewater pH and alkalinity were significantly elevated in acidic soils treated with steel slag relative to control (Figure 4; all differences are significant at the p < 0.05 level using a Tukey HSD and ANOVA test). There was an increase in Ca and Mg concentrations in slag treated relative to control porewaters in 2022 and 2023, but the difference was not significant in either year (p = 0.63 and p = 0.14 in 2022 and 2023, respectively; see Figure 4). We observed no difference between porewater chemistry in the acidic soils treated with basalt during the 2022 and 2023 growing season, and lysimeters were repurposed in 2024 to focus only on continued monitoring of steel slag treated soil compared to control soil. There were consistent, measurable concentrations of strong acid anions such as nitrate (NO3−) throughout the trial, but there were no significant differences between non-alkaline anion concentrations observed in any of the monitoring years. However, as with calcium and magnesium concentrations, nitrate concentrations in slag treatments in 2023 trended higher than control plots (p = 0.2; see Figure 4).
Figure 4. Porewater chemistry for acidic soil in southern plots (Blocks C and D). Boxplots represent aggregated data across a summer monitoring season in each indicated sampling year (2022, 2023, and 2024). In 2024 only slag and control plots were monitored due to lack of an observable difference in porewater chemistry between basalt and control during the 2022 and 2023 growing season. Stars indicate statistical significance (p < 0.05; Tukey HSD and ANOVA test) between the indicated treatment relative to control. Note that where basalt is significantly different in 2022 in calcium and magnesium concentrations, they are lower than those observed in control, an effect that is opposite that is expected.
In the cases where porewater was measured for heavy metals including Cr, Pb, and Ni, these metals were often below the analytical detection limits (data for Cr, Ni, and Pb are available in the Supplementary file). Ni is the only metal with enough values for statistics, and there was no significant difference with treatment (p-values > 0.05 a Tukey HSD and ANOVA test).
3.3 CO2 gas flux
There was no significant difference in the CO2 gas flux between basalt-amended or slag-amended soils and control soils when measured using survey chambers with or without soil collars during the 2022 and 2023 monitoring season (p > 0.05 for Tukey HSD and ANOVA test; Supplementary Figure S11).
4 Discussion
Despite an abundance of interest and investment in enhanced rock weathering in recent years, there remains relatively limited published data from monitored field deployments. Overall, we found that after 3 years of monitoring, the steel slag treatment yielded measurable evidence for mineral weathering, decreased acid export, and increased alkalinity export in acidic soils. The effects were less pronounced or not significant in neutral pH soils. The basalt in this study generated small, but not significant, differences relative to untreated controls in all soils. While the outcomes of our study are specific to the feedstocks deployed, the local soil types, and our measurement approach, there are several outcomes that are relevant to enhanced weathering as a CDR strategy more broadly which we discuss below.
4.1 Evidence for carbon dioxide removal
CDR occurs in enhanced weathering through mineral dissolution buffering acidity in soil. In soil porewaters, pH and base cation concentrations should increase as mineral weathering occurs (see Equation 2 above). If carbonic acid is directly buffered, bicarbonate (HCO3−) formation will drive an immediate increase in the alkalinity of the system. In soils amended with steel slag, we observed a significant increase in porewater pH and alkalinity in blocks with acidic soil (Blocks C and D) across 3 years of monitoring (Figure 4; p < 0.05 for Tukey HSD and ANOVA). While there was an apparent increase in base cation concentration that accompanied the changes in pH and alkalinity, this effect was only significant in 2024 (Figure 4). In addition to soil porewater, soil pH and Ca-saturation for steel slag treatments in Blocks C and D increased significantly in Year 3 (2024) relative to control (Figure 2), providing further support for sustained weathering and CDR over the three-year monitoring period. However, we did not observe significant effects of steel slag in soils with neutral pH (Blocks A and B) nor basalt treatments regardless of initial soil pH (Figures 2, 4).
In the acidic soils (Blocks C and D), there was consistent evidence for increased pH and increased alkalinity export in steel slag-amended plots, even at a relatively low application rate (2 tonnes acre −1) compared to the much higher application rate of crushed basalt (10 tonnes per acre−1). This highlights the importance of considering the reactivity of materials in feedstock selection. Steel slag is composed of calcium and magnesium oxides and hydroxides, Ca-silicates, and Ca-rich amorphous phases that are highly weatherable and alkaline, making it attractive from the perspective of pH management and CDR (Das et al., 2020; Das et al., 2019; O’Connor et al., 2021). However, steel slag availability in the United States is limited and a large majority of steel slag produced is already utilized within construction markets. Legacy steel slag piles could represent considerable volumes of material but may have more complicated production histories and subsequent potential for heavy metal contamination.
While many studies have focused on basalt as an attractive feedstock for enhanced weathering, the results presented here fail to demonstrate any impact from the applied basalt. This is likely due to the physical and chemical characteristics of the particular basalt used in this study. The grain size of the quarry sand product was likely too coarse-grained to weather measurably in the field over a 3-year period. In addition, the metamorphic minerals epidote and chlorite present in our basalt suggest low-grade metamorphism that may lessen the materials reactivity. In addition, our application rate for basalt (10 tonnes per acre) was relatively low compared to other studies using basalt, and basalt was only applied a single application in the first year, compared to other enhanced weathering deployments with annual applications at higher rates (Kantola et al., 2023; Larkin et al., 2022). The lower total application volume could have limited our ability to detect a weathering signal. Higher application rates come with agronomic and economic trade-offs, especially considering typical application rates for ag lime are just 0.5–2 tonnes per acre for pH management. Increased trucking costs and soil compaction for spreading large volumes of material, particularly for annual applications, are important considerations. This is an obvious documentation that some fine grained commercially available basalt products will not be suitable for EW deployments. Future work to evaluate new basalt products at finer grain sizes in similar settings should help to determine if weathering signals can be detected at lower application rates with basalt.
There are two plausible explanations for why simple measurements failed to detect impacts from enhanced weathering amendments in soils with neutral pH. First, silicate weathering is pH-sensitive and tends to slow in neutral soils, where lower proton availability reduces mineral dissolution and the release of weathering products. Second, more neutral soils can contain soil carbonates (Supplementary Figure S5) which appear to dominate the soil porewater chemistry (Supplementary Figure S4) and have faster weathering kinetics than silicates. This is consistent with studies monitoring chemical weathering in catchments based on stream water chemistry, where carbonate weathering often predominates the weathering signal. Given the more rapid weatherability of carbonates, they require less rock material to be produced, transported, and spread to achieve a similar buffering effect for soil acidity compared to most silicate feedstocks for enhanced weathering. Prior work that evaluated the carbon removal effects of widespread ag lime usage (Hamilton et al., 2007; Oh and Raymond, 2006), and more recent work on carbonate-rich returned concrete as a feedstock for enhanced weathering (McDermott et al., 2024), both support relatively rapid rates of carbonate weathering that can lead to CDR. Given the requirements of voluntary carbon markets for comparison of project activities against a counterfactual (what would have happened absent the project), quantifying net CDR should consider comparing silicates against carbonates in areas where ag lime use is common (e.g., the Midwest U.S.).
Increased Ca-saturation and soil pH provides direct evidence that steel slag is actively weathering and releasing base cations into the soil environment. Given that Ca-saturation increased (Figure 2), at least a portion of the released calcium is adhering to exchange sites, rather than flushing through the soil column. The timescale and rate at which this calcium will migrate through the soil column is not clear, but will be important for full delivery of CDR (Kanzaki et al., 2025). To further quantify rates of cation release and weathering, we applied the Ti-CAT methodology (Reershemius et al., 2023). However, the results were inconclusive using both the high- and low-resolution soil sampling methodologies described above (see Section 2.3 and Supplementary Figure S2). This could be due in part to our relatively low application rates, low Ti abundance in steel slag, or lack of high-resolution sampling for baseline soil prior to application for a direct comparison to weathered mixtures. This serves as an example of the importance of high-resolution sampling for soil-based mass balance approaches.
4.2 Impact on soil CO2 flux and soil organic matter
Inorganic CDR that results from mineral weathering may be synergistic with processes that encourage the storage of organic carbon in soils, similar to the relationship between rock weathering and soil organic carbon accrual over long time scales (Slessarev et al., 2021). However, it remains unclear how different feedstocks will interact with organic carbon cycling in different soil types and climates, and on different timescales. Here we evaluate the impacts of enhanced weathering on soil organic carbon cycling by evaluating changes in total SOM, MAOM and POM pools, and monitoring soil CO2 gas flux from the surface soil.
Both basalt and steel slag did not have any positive or negative effects on total SOM, or MAOM and POM pools (Figures 2, 3). There was an apparent positive increase in MAOM-C in steel slag-amended treatments in neutral soil (Blocks A and B), but not in acidic soil in southern plots (Blocks C and D; Figure 2). In contrast to most other results where steel slag weathering in more acidic soils in Blocks C and D produce significant results, perhaps steel slag in more neutral soils could have a positive effect on MAOM formation.
Interestingly, while steel slag applied to acidic soils (Blocks C and D) led to significant changes to porewater chemistry that indicates weathering and CDR, the increase in MAOM-C in more neutral soils (Blocks A and B; see Figure 3) treated with steel slag could highlight differing responses of inorganic and organic C pools to enhanced weathering feedstocks based on initial soil pH. In the more neutral soils, slag-derived Ca may have promoted MAOM formation through mechanisms like cation bridging, aggregate stabilization or by enhancing microbial transformation of organic inputs (Kang et al., 2024; Rowley et al., 2018; Shabtai et al., 2023). Even if indicators of weathering in soil porewaters were masked by high background carbonate weathering, Ca released from slag could still interact with native SOM. Additionally, higher initial SOM could support MAOM accrual through organo-organic associations (Kang et al., 2024; King, 2025). Future work should aim to determine the effects on SOM pools with higher resolution sampling and better replicated trials that vary in soil pH.
Soil CO2 flux could be influenced by enhanced weathering directly or indirectly. The direct effect of increasing consumption of soil CO2 through mineral weathering could lead to a decrease in the amount of CO2 released from the soil surface. In addition, enhanced weathering-induced increases to soil pH may alter soil microbial activity and influence rates of gas flux in soils (Malik et al., 2018; Vermeiren et al., 2025; Fernández-Calviño and Bååth, 2010). However, there were no significant differences in CO2 flux for any treatments in this study. The lack of any significant difference between basalt and steel slag treatments and control plots suggests there are not any significant losses of carbon through increase soil CO2 flux resulting from enhanced weathering. Although these results are inconclusive, they do not support the hypothesis that mineral addition leads to soil organic carbon loss.
As for any direct impact, the background variability in soil CO2 flux from biological respiration likely masks any signal from inorganic CDR via enhanced weathering. Diurnal cycles can show variation of ~3–5 μmol CO2 m2 s−1, or even more. Spatially, there is also considerable variation even within a given area or experimental plot. Even a low-end rate of CO2 flux in soils from this study of just ~2 μmol CO2 m2 s−1 implies an annual CO2 flux of ~7.5 tCO2 acre−1 yr.−1 being respired out of soils over a frost-free 9-month period. For enhanced weathering amendments to have a chance at being detectable via gas flux measurements, annual rates of carbon removal would need to be perhaps 10–25% of this rate, implying 0.75–1.9 tCO2 acre−1 yr.−1, or higher depending on resolution of measurements. Automated chambers with appropriate spatial resolution and replication could perhaps help to reduce the signal-to-noise issues we encountered here, although existing commercial chambers are expensive and unlikely to be applicable at scale without development of low-cost sensor alternatives. Notably, recent work used an automated array of CO2 sensors embedded in the soil profile to directly detect decreases in CO2 flux based on application of basalt in an enhanced weathering trial in Viriginia (Milliken et al., 2025), suggesting this approach could be developed into an effective monitoring strategy.
4.3 Influence of strong acids on weathering
Soil pH management in agriculture is dominantly a response to the generation of strong acids from the degradation of fertilizers. As an example, ammonia (NH4+), a common ingredient in nitrogen fertilizers, releases acidity during oxidation and forms nitric acid (HNO3), which dissociates into acidity (H+) and nitrate ions (NO3−) in most soil environments (Equation 3).
Like carbonic acid, strong acids weather silicate minerals and release base cations (Equation 4 below). Instead of achieving charge balance with bicarbonate, strong acid anions (e.g., nitrate, NO3−; phosphate, PO42−; sulfate, SO42−) in porewaters charge balance cations. To illustrate this process, consider the stoichiometric reactions below (Equations 4–6):
Given that carbonic acid neutralization and bicarbonate generation is typically thought of as a key mechanism for delivering CDR in enhanced weathering, strong acid weathering introduces complexity for quantifying CDR in enhanced weathering deployments.
Pure carbonic acid weathering should yield a 1:1 equivalence ratio between base cations and alkalinity (bicarbonate; Equation 4), while strong acid weathering will act to increase this ratio by generating an abundance of base cations compared to alkalinity (McDermott et al., 2024; Hamilton et al., 2007). In each year of this study, base cations in soil porewaters exceeded measured alkalinity (see Figure 5). Nonetheless, in all cases, the base cation to alkalinity ratios exceeded 1:1 and in many cases exceeded 2:1 (expected for a mixed acid weathering regime that includes carbonic acid, Equation 6; Figure 5). The apparent strong acid effects were most pronounced in 2023 when the field was planted with corn, requiring the addition of N-based fertilizers. In 2022 and 2024, the field was planted with soybeans, which do not require additional N-fertilizer. Nitrate in soil porewaters often better aligns with base cation concentration (Figure 5), indicating the importance of N cycling in agricultural soil acid base balances.
Figure 5. Porewater alkalinity (left panel set) and nitrate (right panel set) against base cations calcium and magnesium for each year of field trial monitoring for acidic soils only. Solid lines represent a 1:1 relationship expected in the case of pure weathering with carbonic acid (alkanity plots) and nitric acid (nitrate plots). Dashed lines represent 2:1 relationship expected with a mixed carbonic strong acid weathering (Hamilton et al., 2007). Soybeans were planted in 2022 and 2024, which do not require direct fertilization with N-based fertilizers. Corn was planted in 2023 and N-based fertilizers were applied.
Alkalinity loss (ΔHCO3−) is defined as the reduction in alkalinity due to the charge balance of base cations with strong acid anions rather than bicarbonate (Perrin et al., 2008). We assume titrated alkalinity is equivalent to the concentration of HCO3− (in units of meq L−1; referred to below as HCO3−measured). The theoretical concentration of HCO3− (HCO3−theoretical) is assumed to be equal to the concentration of Ca and Mg (in units of meq L−1) and ΔHCO3− is calculated as in Equation 7:
Comparison of alkalinity loss to nitrate concentrations reveals a strong linear relationship that is close to a 1:1 line (see Figure 6), suggesting that N-fertilizers are dominantly responsible for alkalinity losses. The observation of this correlation in controls, as well as treated soils, highlights that acid addition in many soils likely leads to alkalinity loss as bicarbonate neutralizes acidity.
Figure 6. Alkalinity loss (∆HCO3−) plotted against porewater nitrate across all years of the study. The solid line represents a 1:1 ratio showing strong association of nitrate with alkalinity loss.
In acidic soils with low pH, the use of silicate minerals to buffer acidity and bring soils back towards neutral would be particularly attractive relative to carbonate rocks as there would be no risk of releasing fossil CO2 from carbonates during neutralization.
The weathering of silicate minerals with strong acids does not lead to a direct removal of CO2 and subsequent formation of bicarbonate. However, the buffering effect of a silicate mineral on strong acids like nitric acid should decrease CO2 evasion due to acidification and help to retain bicarbonate alkalinity. For instance, consider the carbonate system present in soils (Equation 8):
The net effect of acidification due to nitrogen fertilizers would be to add H+, pushing the equilibrium towards CO2 (Equation 8). From a purely inorganic carbon perspective then, acidified systems will degas CO2. Buffering acidity (H+) from any source of acid (carbonic acid or, in this case, nitric acid) should limit CO2 degassing and drive the carbonate system towards bicarbonate thereby encouraging the stability and preservation of alkalinity. The results here support the inclusion of robust accounting for strong acid effects and acid neutralization in frameworks for greenhouse gas accounting for enhanced weathering deployments.
5 Conclusion
Our three-year field trial of enhanced weathering in the Midwestern U.S. provides evidence that steel slag can deliver measurable carbon dioxide removal (CDR) in acidic soils (initial pH < 6.8) with no carbonate mineral contents. In these settings, steel slag amendments raised soil pH and Ca saturation and elevated soil porewater pH and alkalinity indicating effective weathering leading to CDR at modest application rates (2 tonnes acre−1). In contrast, no CDR or weathering signal was detected from basalt amendments, even at higher application rates (10 tonnes acre−1), or from steel slag in soils with neutral pH (~7) and pre-existing soil carbonate minerals.
Strong acid driven weathering and/or re-equilibration of the carbonic acid system, primarily from N-fertilizers, was apparent in porewater chemistry during the entire three-year period of monitoring. Strong acid production drives alkalinity loss and blocks the initial removal of CO2 to bicarbonate. But, buffering acid likely has downstream impacts on the carbonate system that eventually leads to alkalinity preservation (reduced CO2 degassing), an effect that will likely be difficult to measure.
Measurement techniques for enhanced weathering require high temporal and spatial resolution to detect a relatively small signal in an otherwise noisy system. For porewater chemistry, the dataset presented here was limited due to the labor-intensive nature of sampling and the limited capacity for water recovery in non-saturated media. Soil sampling densities required to detect significant effects were greater than originally anticipated and could indicate substantial time and cost to apply at larger deployment scales. For the measurement techniques that failed to detect a significant signal (chamber-based gas flux measurement and cation-loss mass balance approaches), density of sampling and measurement arrays likely limited the utility of these approaches in this study. Our results together emphasize the importance of high density measurement for reliable single field CDR estimates.
Enhanced weathering continues to emerge as a promising pathway for CDR, but it is critical to better constrain CDR efficiencies in field trials where variable feedstocks are tested in combination with differing soil types and climates. This study provided clear documentation that not all commercially available feedstocks are suitable for deployments. In addition, continued development of sampling and measurement strategies that robustly detect weathering, quantify CDR to develop scalable strategies for monitoring, reporting, and verification are key as enhanced weathering become more widespread.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material.
Author contributions
DM: Supervision, Methodology, Investigation, Data curation, Conceptualization, Writing – review & editing, Funding acquisition, Resources, Formal analysis, Writing – original draft, Project administration, Visualization. EM: Methodology, Investigation, Writing – review & editing, Conceptualization. JY: Methodology, Conceptualization, Investigation, Writing – review & editing. EW: Investigation, Formal analysis, Methodology, Writing – original draft. RG: Writing – review & editing, Investigation, Methodology. LS: Methodology, Investigation, Writing – review & editing. JS: Investigation, Writing – review & editing, Methodology. NS: Investigation, Methodology, Funding acquisition, Writing – review & editing. NP: Methodology, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. Funding for this project was provided by Carleton College, the Keck Geology Consortium, and the National Science Foundation under Grant Nos. 2050697 and 2208133. Work at LLNL by NWS was funded by the LLNL LDRD Program (22-LW-022, 24-SI-002) and supported by LLNL’s ‘Terraforming Soil’ Energy Earthshot Research Center (EERC), funded by the U.S. Department of Energy Office of Science, Award# SCW1841, and was performed under the auspices of the DOE, Contract DE-AC52-07NA27344.
Acknowledgments
The authors thank several individuals who have contributed to this project, including Adele Fredericks, Demetrius Blackmon-Jimenez, Sebastian Forero Escovar., Jingwen (Crystal) Wu, Anna Greenlee, Ella Cunningham, Fiona Anstey, Sarah Leibovitz, Sophie Naylor, Mike and Shane Peterson, Jonathon Cooper, and Kenzo Esquivel.
Conflict of interest
NP was a co-founder of the CDR company Lithos Carbon but does not have any direct financial tie or equity stake in the company.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The authors declare that no Gen AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fclim.2025.1657058/full#supplementary-material
Footnotes
References
Almaraz, M., Bingham, N. L., Holzer, I. O., Geoghegan, E. K., Goertzen, H., Sohng, J., et al. (2022). Methods for determining the CO2 removal capacity of enhanced weathering in agronomic settings. Front. Clim. 4:970429. doi: 10.3389/fclim.2022.970429
Beerling, D. J., Epihov, D. Z., Kantola, I. B., Masters, M. D., Reershemius, T., Planavsky, N. J., et al. (2024). Enhanced weathering in the US Corn Belt delivers carbon removal with agronomic benefits. Proc. Natl. Acad. Sci. USA 121:e2319436121. doi: 10.1073/pnas.2319436121,
Beerling, D. J., Kantzas, E. P., Lomas, M. R., Taylor, L. L., Zhang, S., Kanzaki, Y., et al. (2025). Transforming US agriculture for carbon removal with enhanced weathering. Nature 638, 425–434. doi: 10.1038/s41586-024-08429-2,
Beerling, D. J., Leake, J. R., Long, S. P., Scholes, J. D., Ton, J., Nelson, P. N., et al. (2018). Farming with crops and rocks to address global climate, food and soil security. Nat Plants 4, 138–147. doi: 10.1038/s41477-018-0108-y
Bradford, M. A., Fierer, N., and Reynolds, J. F. (2008). Soil carbon stocks in experimental mesocosms are dependent on the rate of labile carbon, nitrogen and phosphorus inputs to soils. Funct. Ecol. 22, 964–974. doi: 10.1111/j.1365-2435.2008.01404.x
Buss, W., Hasemer, H., Sokol, N. W., Rohling, E. J., and Borevitz, J. (2024). Applying minerals to soil to draw down atmospheric carbon dioxide through synergistic organic and inorganic pathways. Commun. Earth Environ. 5:602. doi: 10.1038/s43247-024-01771-3
CarbonPlan Does enhanced weathering work? We’re still learning. Available online at: https://carbonplan.org/research/enhanced-weathering-fluxes (Accessed May 14, 2025).
Chiaravalloti, I., Theunissen, N., Zhang, S., Wang, J., Sun, F., Ahmed, A. A., et al. (2023). Mitigation of soil nitrous oxide emissions during maize production with basalt amendments. Front. Clim. 5:1203043. doi: 10.3389/fclim.2023.1203043
Clarkson, M. O., Larkin, C. S., Swoboda, P., Reershemius, T., Suhrhoff, T. J., Maesano, C. N., et al. (2024). A review of measurement for quantification of carbon dioxide removal by enhanced weathering in soil. Front. Clim. 6:6. doi: 10.3389/fclim.2024.1345224,
Das, S., Galgo, S. J., Alam, M. A., Lee, J. G., Hwang, H. Y., Lee, C. H., et al. (2020). Recycling of ferrous slag in agriculture: potentials and challenges. Crit. Rev. Environ. Sci. Technol. 1, 1–35. doi: 10.1080/10643389.2020.1853458
Das, S., Kim, G. W., Hwang, H. Y., Verma, P. P., and Kim, P. J. (2019). Cropping with slag to address soil, environment, and food security. Front. Microbiol. 10:1320. doi: 10.3389/fmicb.2019.01320,
Edwards, D. P., Lim, F., James, R. H., Pearce, C. R., Scholes, J., Freckleton, R. P., et al. (2017). Climate change mitigation: potential benefits and pitfalls of enhanced rock weathering in tropical agriculture. Biol. Lett. 13:20160715. doi: 10.1098/rsbl.2016.0715,
Fernández-Calviño, D., and Bååth, E. (2010). Growth response of the bacterial community to pH in soils differing in pH: growth response of the bacterial community to pH. FEMS Microbiol. Ecol. doi: 10.1111/j.1574-6941.2010.00873.x
Goll, D. S., Ciais, P., Amann, T., Buermann, W., Chang, J., Eker, S., et al. (2021). Potential CO2 removal from enhanced weathering by ecosystem responses to powdered rock. Nat. Geosci. 14, 545–549. doi: 10.1038/s41561-021-00798-x
Gunnarsen, K. C., Schjoerring, J. K., Gómez-Muñoz, B., de Neergaard, A., and Jensen, L. S. (2022). Can silicon in glacial rock flour enhance phosphorus availability in acidic tropical soil? Plant Soil 477, 241–258. doi: 10.1007/s11104-022-05399-0
Guo, F., Wang, Y., Zhu, H., Zhang, C., Sun, H., Fang, Z., et al. (2023). Crop productivity and soil inorganic carbon change mediated by enhanced rock weathering in farmland: a comparative field analysis of multi-agroclimatic regions in Central China. Agric. Syst. 210:103691. doi: 10.1016/j.agsy.2023.103691
Hamilton, S. K., Kurzman, A. L., Arango, C., Jin, L., and Robertson, G. P. (2007). Evidence for carbon sequestration by agricultural liming. Glob. Biogeochem. Cycles 21:21. doi: 10.1029/2006GB002738,
Haque, F., Santos, R. M., Dutta, A., Thimmanagari, M., and Chiang, Y. W. (2019). Co-benefits of wollastonite weathering in agriculture: CO2 sequestration and promoted plant growth. ACS Omega 4, 1425–1433. doi: 10.1021/acsomega.8b02477,
Hartmann, J., and Kempe, S. (2008). What is the maximum potential for CO2 sequestration by “stimulated” weathering on the global scale? Naturwissenschaften 95, 1159–1164. doi: 10.1007/s00114-008-0434-4,
Hartmann, J., West, A. J., Renforth, P., Köhler, P., De La Rocha, C. L., Wolf-Gladrow, D. A., et al. (2013). Enhanced chemical weathering as a geoengineering strategy to reduce atmospheric carbon dioxide, supply nutrients, and mitigate ocean acidification. Rev. Geophys. 51, 113–149. doi: 10.1002/rog.20004
Holzer, I. O., Nocco, M. A., and Houlton, B. Z. (2023). Direct evidence for atmospheric carbon dioxide removal via enhanced weathering in cropland soil. Environ. Res. Commun. 5:101004. doi: 10.1088/2515-7620/acfd89
Kang, J., Qu, C., Chen, W., Cai, P., Chen, C., and Huang, Q. (2024). Organo–organic interactions dominantly drive soil organic carbon accrual. Glob. Change Biol. 30:e17147. doi: 10.1111/gcb.17147,
Kantola, I. B., Blanc-Betes, E., Masters, M. D., Chang, E., Marklein, A., Moore, C. E., et al. (2023). Improved net carbon budgets in the US Midwest through direct measured impacts of enhanced weathering. Glob. Chang. Biol. 29, 7012–7028. doi: 10.1111/gcb.16903,
Kanzaki, Y., Planavsky, N., Zhang, S., Jordan, J., Suhrhoff, T. J., and Reinhard, C T. (2025). Soil cation storage as a key control on the timescales of carbon dioxide removal through enhanced weathering. Environ. Res. Lett. 20:074055. doi: 10.1088/1748-9326/ade0d5
King, A. (2025). Soil carbon formation is promoted by saturation deficit and existing mineral-associated carbon, not by microbial carbon-use efficiency. Sci. Adv. 11:eadv9482. doi: 10.1126/sciadv.adv9482,
Köhler, P., Hartmann, J., and Wolf-Gladrow, D. A. (2010). Geoengineering potential of artificially enhanced silicate weathering of olivine. Proc. Natl. Acad. Sci. USA 107, 20228–20233. doi: 10.1073/pnas.1000545107,
Larkin, C. S., Andrews, M. G., Pearce, C. R., Yeong, K. L., Beerling, D. J., Bellamy, J., et al. (2022). Quantification of CO2 removal in a large-scale enhanced weathering field trial on an oil palm plantation in Sabah, Malaysia. Front. Clim. 4:959229. doi: 10.3389/fclim.2022.959229
Malik, A. A., Puissant, J., Buckeridge, K. M., Goodall, T., Jehmlich, N., Chowdhury, S., et al. (2018). Land use driven change in soil pH affects microbial carbon cycling processes. Nat. Commun. 9:3591. doi: 10.1038/s41467-018-05980-1,
McBride, A. L., Skov, K., Wade, P., Betz, J., Stubbs, A., Albahri, T., et al. Quantifying potential carbon dioxide removal via enhanced weathering using porewater from a field trial in Scotland. EarthArXiv 2025 [Preprint]. doi: 10.31223/X5NX5H.
McDermott, F., Bryson, M., Magee, R., and Van Acken, D. (2024). Enhanced weathering for CO2 removal using carbonate-rich crushed returned concrete; a pilot study from SE Ireland. Appl. Geochem. 169:106056. doi: 10.1016/j.apgeochem.2024.106056
Milliken, E., Woodley, A., and Planavsky, N. (2025). Direct measurement of carbon dioxide removal due to enhanced weathering. CDRXIV [Preprint]. Environmental Science and Technology. 12, 970–976. doi: 10.1021/acs.estlett.5c00441
Neal, C. (2001). Alkalinity measurements within natural waters: towards a standardised approach. Sci. Total Environ. 265:99. doi: 10.1016/S0048-9697(00)00652-5,
O’Connor, J., Nguyen, T. B. T., Honeyands, T., Monaghan, B., O’Dea, D., Rinklebe, J., et al. (2021). Production, characterisation, utilisation, and beneficial soil application of steel slag: a review. J. Hazard. Mater. 419:126478. doi: 10.1016/j.jhazmat.2021.126478,
Oh, N., and Raymond, P. A. (2006). Contribution of agricultural liming to riverine bicarbonate export and CO2 sequestration in the Ohio River basin. Glob. Biogeochem. Cycles 20:2005GB002565. doi: 10.1029/2005GB002565
Perrin, A.-S., Probst, A., and Probst, J.-L. (2008). Impact of nitrogenous fertilizers on carbonate dissolution in small agricultural catchments: implications for weathering CO2 uptake at regional and global scales. Geochim. Cosmochim. Acta 72, 3105–3123. doi: 10.1016/j.gca.2008.04.011
Power, I. M., Hatten, V. N. J., Guo, M., Schaffer, Z. R., Rausis, K., and Klyn-Hesselink, H. (2025). Are enhanced rock weathering rates overestimated? A few geochemical and mineralogical pitfalls. Front. Clim. 6:1510747. doi: 10.3389/fclim.2024.1510747
PP Systems 2018 EGM-5 portable CO2 gas analyzer operations manual version 1.03. Amesbury, MA, USA: PP Systems.
Reershemius, T., Kelland, M. E., Jordan, J. S., Davis, I. R., D’Ascanio, R., Kalderon-Asael, B., et al. (2023). Initial validation of a soil-based mass-balance approach for empirical monitoring of enhanced rock weathering rates. Environ. Sci. Technol. 57, 19497–19507. doi: 10.1021/acs.est.3c03609,
Rowley, M. C., Grand, S., and Verrecchia, É. P. (2018). Calcium-mediated stabilisation of soil organic carbon. Biogeochemistry 137, 27–49. doi: 10.1007/s10533-017-0410-1
Scharlemann, J. P., Tanner, E. V., Hiederer, R., and Kapos, V. (2014). Global soil carbon: understanding and managing the largest terrestrial carbon pool. Carbon Manag. 5, 81–91. doi: 10.4155/cmt.13.77
Shabtai, I. A., Wilhelm, R. C., Schweizer, S. A., Höschen, C., Buckley, D. H., and Lehmann, J. (2023). Calcium promotes persistent soil organic matter by altering microbial transformation of plant litter. Nat. Commun. 14:6609. doi: 10.1038/s41467-023-42291-6,
Skov, K., Wardman, J., Healey, M., McBride, A., Bierowiec, T., Cooper, J., et al. (2024). Initial agronomic benefits of enhanced weathering using basalt: a study of spring oat in a temperate climate. PLoS One 19:e0295031. doi: 10.1371/journal.pone.0295031,
Slessarev, E. W., Chadwick, O. A., Sokol, N. W., Nuccio, E. E., and Pett-Ridge, J. (2021). Rock weathering controls the potential for soil carbon storage at a continental scale. Biogeochemistry 157, 1–13. doi: 10.1007/s10533-021-00859-8
Smith, S., Geden, O., Gidden, M. J., Lamb, W. F., Nemet, G. F., Minx, J. C., et al. (eds.) (2024). The state of carbon dioxide removal 2024 - 2nd edition. doi: 10.17605/OSF.IO/F85QJ
Sohng, J., Sokol, N. W., Whiteaker, S., Schmidt, R., Holzer, I., Goertzen, H., et al. (2025). Combining organic amendments with enhanced rock weathering shifts soil carbon storage in croplands. Sci. Total Environ. 998:180179. doi: 10.1016/j.scitotenv.2025.180179
Soil Survey Staff United States Department of Agriculture. Web Soil Survey - Home. Available online at: https://websoilsurvey.nrcs.usda.gov/app/homepage.htm (Accessed May 20, 2025).
Sokol, N. W., Kuebbing, S. E., and Bradford, M. A. (2017). Impacts of an invasive plant are fundamentally altered by a co-occurring Forest disturbance. Ecology 98, 2133–2144. doi: 10.1002/ecy.1906,
Sokol, N. W., Sohng, J., Moreland, K., Slessarev, E., Goertzen, H., Schmidt, R., et al. (2024). 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. Biogeochemistry 167, 989–1005. doi: 10.1007/s10533-024-01160-0
Suhrhoff, T. J., Reershemius, T., Wang, J., Jordan, J. S., Reinhard, C. T., and Planavsky, N. J. (2024). A tool for assessing the sensitivity of soil-based approaches for quantifying enhanced weathering: a US case study. Front. Clim. 6:6. doi: 10.3389/fclim.2024.1346117,
Taylor, L. L., Driscoll, C. T., Groffman, P. M., Rau, G. H., Blum, J. D., and Beerling, D. J. (2021). Increased carbon capture by a silicate-treated forested watershed affected by acid deposition. Biogeosciences 18, 169–188. doi: 10.5194/bg-18-169-2021
Vermeiren, C., Ceulemans, J., Thiry, Y., and Smolders, E. (2025). Increased soil pH and enhanced microbial activity stimulate the gradual immobilisation of selenate added to soils. Soil Biol. Biochem. 202:109688. doi: 10.1016/j.soilbio.2024.109688
Wirth, K. R., Naiman, Z. J., and Vervoort, J. D. (1997). The Chengwatana volcanics, Wisconsin and Minnesota: petrogenesis of the southernmost volcanic rocks exposed in the midcontinent rift. Can. J. Earth Sci. 34, 536–548. doi: 10.1139/e17-043
Keywords: enhanced weathering, carbon dioxide removal, basalt, steel slag, soil carbon
Citation: Maxbauer DP, Milliken E, Yambing JR, Watson E, Gregg RB, Swanson L, Sohng J, Sokol NW and Planavsky NJ (2026) Evidence for carbon dioxide removal via enhanced rock weathering with steel slag, though not basalt, in a midwestern U.S. field trial. Front. Clim. 7:1657058. doi: 10.3389/fclim.2025.1657058
Edited by:
Yiwen Pan, Zhejiang University, ChinaReviewed by:
Nataliya Yurkevich, Institute of Petroleum Geology and Geophysics (RAS), RussiaLianghan Cong, University of Vaasa, Finland
Copyright © 2026 Maxbauer, Milliken, Yambing, Watson, Gregg, Swanson, Sohng, Sokol and Planavsky. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Daniel P. Maxbauer, ZG1heGJhdWVyQGNhcmxldG9uLmVkdQ==
Jahmaine Renzo Yambing1,3,4