ORIGINAL RESEARCH article
Sec. Negative Emission Technologies
Enhanced weathering potentials—the role of in situ CO2 and grain size distribution
- 1Center for Earth System Sciences and Sustainability, Institute for Geology, Universität Hamburg, Hamburg, Germany
- 2Université Grenoble Alpes, CNRS, ISTerre, Grenoble, France
- 3Alfred-Wegener-Institut Helmholtz-Zentrum für Polar- und Meeresforschung Bremerhaven, Bremerhaven, Germany
- 4Potsdam Institute for Climate Impact Research, Member of the Leibniz Association, Potsdam, Germany
The application of rock powder on agricultural land to ameliorate soils and remove carbon dioxide (CO2) from the air by chemical weathering is still subject to many uncertainties. To elucidate the effects of grain size distribution and soil partial pressure of carbon dioxide (pCO2) levels on CO2 uptake rates, two simple column experiments were designed and filled nearly daily with an amount of water that simulates humid tropical conditions, which prevail in areas known for being hotspots of weathering. Multiple materials (dunite, basanite, agricultural oxisol, a combination of the latter two, and loess) were compared under ambient and 100% CO2 atmosphere. In a second series, single material columns (dunite) were filled with three different grain size distributions. Total alkalinity, pH, major ions, and dissolved silica were determined in the outflow water of the columns for about 300 days. Under ambient atmospheric conditions, the CO2 consumption was the lowest in the oxisol column, with 100 t CO2 km−2 year−1, while dunite and basanite showed similar consumption rates (around 220 t CO2 km−2 year−1). The values are comparable to high literature values for ultramafic lithologies. Interestingly, the mixture of basanite and oxisol has a much higher consumption rate (around 430 t CO2 km−2 year−1) than the basanite alone. The weathering fluxes under saturated CO2 conditions are about four times higher in all columns, except the dunite column, where fluxes are increased by a factor of more than eleven. Grain size distribution differences also play a role, with the highest grain surface area normalized weathering rates observed in the columns with coarser grains, which at first seems counterintuitive. Our findings point to some important issues to be considered in future experiments and a potential rollout of EW as a carbon dioxide removal method. Only in theory do small grain sizes of the spread-material yield higher CO2 drawdown potentials than coarser material. The hydrologic conditions, which determine the residence times in the pore space, i.e., the time available for weathering reactions, can be more important than small grain size. Saturated-CO2 column results provide an upper limit for weathering rates under elevated CO2.
Chemical weathering of silicate rocks exposed to meteoric water is a process that removes carbon dioxide (CO2) from the atmosphere (Ebelmen, 1845). To counteract rising CO2 levels by removal of CO2 from the atmosphere, the use of rock products was briefly proposed by Seifritz (1990) and made popular in the climate science community as Enhanced Weathering (EW) by Schuiling and Krijgsman (2006). Since then, assessing the realistic potential of EW as a tool for carbon dioxide removal (CDR) has become the focus of literature reviews (Hartmann et al., 2013) and modeling studies (e.g., Hartmann and Kempe, 2008; Hartmann et al., 2013; Taylor et al., 2015; Beerling et al., 2018; Strefler et al., 2018a). However, these studies have relied on several assumptions for the parameterization of the weathering process. Moreover, robust experimental data are still missing, aside from a few laboratories (Renforth et al., 2015), mesocosm (ten Berge et al., 2012; Dietzen et al., 2018; Amann et al., 2020; Kelland et al., 2020; Vienne et al., 2022), and field trials (Haque et al., 2020).
In addition to CDR, the application of rock products has also been thought to improve soil conditions (via fertilization and physical amelioration effects), thereby increasing agricultural production (Julius, 1894; Leonardos et al., 1987; van Straaten, 2002). The resulting increase in biomass growth leads to additional CO2 uptake and storage as organic carbon, thereby augmenting the overall CDR potential of EW.
The sensitivity of weathering reactions to CO2 concentrations in a system, i.e., the magnitude of chemical weathering effects from CO2 concentration changes, has implications for estimates of weathering rates of EW applications, and therefore CO2 removal. The parameterization of Earth system models targeting long time scales relies on, for example, realistic assumptions of the CO2 and related climate sensitivity (e.g., Royer et al., 2007). In general, the sensitivity of weathering rates to CO2 levels also has consequences for the evolution of the atmospheres on other planets and is important for, e.g., habitability assessments for early Mars (Teitler et al., 2014; Cockell et al., 2016; Foley and Smye, 2018).
Small grain sizes (i.e., with high surface areas) are favorable for fast reaction rates in the system, but the comminution of rock to fine grain sizes implies energy costs that may render EW less attractive. Thus, an ideal grain size range that balances these two constraint needs must be defined for different types of source material and soil targets. This energy constraint may not be relevant if side products from the mining process, often deposited as fines, are used.
To address these issues, column experiments were set up to better understand the effects of the grain size distribution (also referred to as particle size distribution) and partial pressure of CO2 in soils with respect to CO2 uptake rates, these being important controls on weathering rates in the field, and thus the overall CDR potential. One of the objectives of this study was therefore to test realistic upper limits of weathering rates and CDR potential considering CO2 levels as the primary driver of weathering rates.
To explore the CO2 effect, different materials were treated under humid tropical conditions, typical for weathering hotspots, with two soil CO2 conditions, ambient (0.046%) and saturated (100%) CO2, where the latter is an upper limit for this parameter, which is of course not realized in soils. Four rock types were used: dunite, basanite, laterite, and loess. Olivine-rich dunite was selected as a chemically simple model rock type that is primarily composed of olivine, while basaltoids (like the basanite used here) were proposed as a more suitable rock type for EW (Hartmann et al., 2013; Amann and Hartmann, 2019) that is widely available (Börker, 2019). Laterite or oxisol was included as a representative of typically abundant soil types associated with weathering hotspots in humid tropical regions. As a commonly proposed EW substrate, this laterite was used in pure form and as a laterite-basanite mixture. The last material type is loess sediment from aeolian deposits—these can be considered a natural analog for EW. The effect of grain sizes and their distribution was tested using the dunite with three different grain size distributions under fixed ambient CO2 conditions.
Material origin and characterization
Materials used in the experiments were dunite, basanite, laterite, a basanite-laterite mix, and loess. The bulk chemical composition of all materials (Table 1) was determined by Panalytical Magix Pro wavelength dispersive X-ray fluorescence (XRF) analysis. Specific surface areas (Table 2) were measured by BET analyses using Kr adsorption (Brunauer et al., 1938) with a Quantachrome autosorb iQ. Particle size characteristics (Table 2, Figure 1) were determined using Sympatec Helos KFMagic laser granulometry coupled to a Quixel wet dispersing unit, covering a range of 0.5/18–3,200 μm.
Table 2. Specific surface areas of the source material, derived by BET analyses, as well as grain size distribution characteristics.
The dunite rock used here, also deployed in Amann et al. (2020), contains ~93% forsterite (Mg-endmember olivine, Mg2SiO4). The other 7% are comprised of lizardite (Mg-rich serpentine), Cr-bearing chlorite (including chromite or chrome-spinel inclusions), and traces of chabazite (zeolite group) and Mg-hornblende (amphibole). The mineralogical composition was determined by energy-dispersive X-ray spectroscopy (Zeiss LEO 1455 VP coupled with an Oxford Instruments EDX detector). The sample originates from the Almklovdalen peridotite complex (Åheim mineral deposit mined by the North Cape Minerals Company, Norway). More information on the geochemistry of the material can be found in Beyer (2006); Hövelmann et al. (2012) describes the geological setting. The material was provided in two-grain sizes, a fine and a coarse fraction (Figure 1).
The material used is a commercially available rock powder called Eifelgold, produced by Lava Union GmbH (https://www.rpbl.de) and sourced from the Eifel, a low mountain range in western Germany with notable volcanic rock occurrences. It was used in an “as delivered” condition with a grain size similar to the fine dunite, with 44 μm as the dominant grain size fraction (Figure 1). Based on the geochemical composition the calculated normative mineralogy (based on the CIPW norm by Cross et al. (1902) is diopside (32%), plagioclase (16%), leucite (15%), olivine (14%), nepheline (14%), and others (10%). Based on the modal mineralogical QAPF classification of volcanic rocks (Streckeisen, 1980) or the TAS classification (Le Bas et al., 1986), this nepheline-rich material is of basaltic composition, and more specifically, it is classified as basanite.
This material, received from coffee farmers through private contact, originates from the area of Conceição de Ipanema, Minas Gerais, Brazil. According to the USDA soil taxonomy (Soil Survey Staff, 1999), it can be considered to be oxisol or ultisol, yet a more detailed soil description is not available. We use the more generic term laterite here and chose this material based on its cation-depleted character and high iron and aluminum content. As this soil material comes from an actively used agricultural plot, it is known that some fertilizer had been applied to the soil, but more precise information on the fertilizer type and quantities used is not known. There are two types of laterite (Table 1), which were sampled at different locations of the same plot. One type was used for the pure laterite column experiments, and due to insufficient material, the second type was used for the basanite-laterite mix.
The loess material was freshly collected from a well-described Upper Pleistocene loess occurrence in Remagen, western Germany (50.5611°N, 7.2435°E; Klasen et al., 2015) during a field trip in 2017. It was dried at 40°C to remove moisture. The average CaCO3 content is reported as 17–19%, depending on the formation (Schirmer, 2011).
Acrylic cylindrical tubes (outer-ø/inner-ø/length: 60/56/250 mm) were filled with the untreated materials (Table 3) to a height of 20 cm while tapping the cylinder to gently settle the powder. All cylinders were fixed vertically using clamps and clamp stands (Figure 2). The top/inlet was covered to decrease water loss by evaporation. A 20-μm nylon mesh was fixed at the bottom of each column to allow water to pass through while retaining the fill material. The experiments were carried out in a temperature-regulated room at about 23°C (Table 3 for details) during the duration of each run. In total, 75 ml of deionized water was added (all at once) to each cylinder every morning, 5 days a week. This amounts to about 8,000 mm a−1 (= 8 m3 m−2 a−1), to simulate intensive tropical rainfall (e.g., up to 10,000 mm a−1 on La Réunion (Strasberg et al., 2005) or Hawai'i (Giambelluca et al., 2013). In contrast to these natural settings where a variable, but significant amount of water is surface runoff, water is forced through the column material in the experiments.
Figure 2. Experimental setup showing the semi-automated watering system for the CO2 column experiment.
The effluent solution was collected in polyethylene (PE) bottles, and initially, collected and measured daily. The sampling frequency was then decreased in the later stages of the experiment when the fluid composition showed decreasing variability. The experiments ran for about 10 months (Table 3). Samples were not pooled over multiple watering events. It should be noted that each daily-collected fluid only represents the water that flowed through the column starting from the day before. The samples were stored in a cool and dark environment until analyses.
One series of columns was treated with deionized water in equilibrium with ambient air. The saturated-CO2 series used deionized water that was sparged with a pure CO2 stream for 15 min before application. These columns were capped with a gas-tight cover.
Stock solutions had an average pH of 5.34 (5.58; note that all values in parentheses were calculated with PhreeqC) and an average electrical conductivity (EC) of 0.9 ± 0.6 μS cm−1 (1.0) when in equilibrium with the ambient lab atmosphere. Under CO2-saturated conditions, the average pH was 3.93 (3.91) and the average EC was 43.9 ± 5.6 μS cm−1 (48.0).
Grain size experiment
This experiment was run with dunite (Table 3). Three replicates of three-grain size distributions were prepared: a coarse fraction, a fine fraction, and a 1:1 mix by mass of the two-grain size fractions. The mass of the dunite powder in each cylinder was subsequently measured.
Fluid composition analysis
The fluid effluent samples were measured for pH, alkalinity, dissolved silica, and major ion composition. Alkalinity was measured by automated titration to pH 4.3, using a Metrohm Titrando. The titration method determines total alkalinity (TA), defined by Dickson (1981) as:
Measurements were calibrated with certified references by A. Dickson (recovery ≥ 99%). Dissolved silica was measured using the molybdate blue colorimetric method (after Hansen and Koroleff, 1999). The average and standard deviation of the coefficients of variation from multiple measurements were 0.6 ± 0.8% for all measured samples. The major cations, Na+, K+, Ca2+, and Mg2+, were determined by ion chromatography (Metrohm 881 Compact IC Pro system). Measurements were considered valid when the standard recovery (done every ten samples) was ≥ 95%. The pH was measured using a WTW handheld pH meter calibrated with 3 buffer standards. The effluent solution was routinely sampled before new water was added, implying that the solution was in equilibrium with the atmosphere (time for equilibration was generally ~24 h).
To estimate the total cation loss from the column material, data gaps in the obtained time series were approximated by linear interpolations between measured values. In the next step, daily fluxes of each major cation were summed over the entire experimental period. To derive the percentage loss of cations relative to the amount present in the initial fresh material, the following expression was employed:
here Q [L] is the solution outflux volume from the column per time step t [day] and c [mol L−1] is the concentration of the considered cation. The moles of a given cation are initially present, ncat [mol], can be calculated using this expression:
where MCat and MO are the molar mass of the given cation and oxygen, mCatχO is the column-specific available mass [g] of the considered cation oxide (Catx%), as derived by XRF and found in Table 1:
with mc [g] being the initial mass of the solid material in the column.
Weathering rate calculations
TA release rate
To calculate the “average” instantaneous release rate of TA from a particular material, we used the following equation:
As is commonly the case, weathering or elemental release rates are typically reported in terms of moles per unit area of solid reactant per unit time (mol m−2 s−1). In the above equation, Q [L] is the solution outflux volume from the column per time step t, c [mol L−1] is the concentration of total alkalinity (TA) at timestep t, ATSA [m2] is the initial total surface area of the column material (which changes over time but is approximated here as being constant), and τ [s] is the residence time of the solution occupying the pore space of the material.
The total surface area is given by:
where mc [g] is the material mass and ASSA [m2 g−1] is the specific surface area.
The residence time is calculated by:
with Vp as the calculated porespace within the column and Vo(t) is the volume of outflowing water per timestep t, with t being the length of each time step (sampling interval).
The pore space is derived from the density differences as:
where ρs is the material-specific density [g m−3] and Vc [m3] is the volume of the empty column. Column mass and material density data are given in Supplementary Material 1.
Equation 5 is equivalent to the standard equation used to calculate mineral dissolution rates using mixed-flow reactors (e.g., Equation 2 in Hellmann, 1994). Nonetheless, the release rates based on Equation 5 provide what can be considered to be a ‘column-averaged' elemental release rate, which in all cases would be less than the corresponding rate based on the dissolution of the same amount of material in an MFR. The reason for this is that as the solution flows from top to bottom in a column, non-linear chemical gradients will develop along the column axis, such that the rate of dissolution of material at the top of the reactor will be higher than that at the bottom. For the purposes of this study, the use of an “average rate” should be an acceptable proxy.
Land surface-based CO2 consumption rate (Rland)
Based on the weathering reactions for two common silicate minerals, forsterite:
we assume a 1:1 conversion of CO2 into alkalinity in the second step of the reaction, which represents the initial weathering product, while the second steps in Equations 9 and 10 represent the long-term final stages of the reaction.
Based on this conversion, we use the TA concentrations together with outflow volumes to calculate a first-order estimate of land surface-based weathering rates, which are equal to the CO2 consumption rates.
The annual CO2 consumption is estimated based on average TA fluxes during the final stages of the experiments, based on the assumption that dynamic near-equilibrium conditions are representative of long-term fluxes. The values used were t > 225 days for dunite, t > 268 days for basanite-laterite, and t > 289 days for basanite, laterite, and loess experiments.
Individual, daily TA fluxes were calculated according to
based on the measured data for outfluxes Q(t) [L d−1] and TA concentrations cTA(t) [mol L−1] from the aforementioned experiments and corresponding periods.
With these individual fluxes, an average annual CO2 consumption (in t km−2 a−1) was calculated by
where n is the number of available flux values, Ac is the cross-sectional surface area of the column, MCO2 is the molar mass of CO2, and the factor of 365 serves to convert daily to annual CO2 uptake.
To evaluate the additional CO2 uptake over the first year, which is characterized by non-equilibrium fluxes, a simplified approach was chosen. Since TA concentrations were not determined every day due to economic and time constraints, the available data were fitted to a non-linear function in order to derive an estimate of the TA fluxes over the course of the experiments (curve fits shown in Supplementary Material 2). Based on the fitting parameters, modeled daily TA concentrations were summed over 1 year:
The additional CO2 consumption in the first year is then simply calculated by
All saturation indices were calculated using the Python package PhreeqPython (Vitens, 2021), which is based on PhreeqC (Parkhurst and Appelo, 2013). Thermodynamic input was taken from the WATEQ4F database (Ball and Nordstrom, 1991). Calculations consider the major cation concentrations, dissolved silica, phosphate, alkalinity, pH, and temperature as input variables.
Total alkalinity and pH
All experiments show an immediate effect of watering after the start of each experiment, for both ambient CO2 and saturated CO2 conditions (Figure 3). There is an initial phase in which all materials show an immediate spike in TA concentration, followed by a slow exponential decrease, and then the attainment of a roughly linear relation characterized by slowly decreasing concentrations (Figures 3A,B). The weak linear (downward) trends remain, even during the final stages of each experiment, indicating that true steady-state conditions were not yet fully reached.
Figure 3. Development of TA, the sum of cation equivalent concentrations, and pH in the columns treated under ambient (A,B,G) and saturated (C,D,H) CO2 conditions. The ratios of TA and total cations (CO2 saturated/CO2 ambient) are shown in (E,F). pH is given in (G,H), note the starting pH for water under ambient CO2 was on average 5.34, while the starting pH under 100% CO2 was on average 3.93.
Except for the first week, the highest TA concentrations under ambient CO2 were observed for basanite, with values around 2,000 μmol L−1 (Figure 3A), whereas under CO2-saturated conditions, dunite showed the highest concentrations, with values starting >100,000 μmol L−1, and then leveling out further on to around 10,000 μmol L−1 (Figure 3C).
The strongest effect of CO2 saturation on TA was observed for dunite, as evidenced by a TA ratio (CO2-saturated/CO2-ambient) equalling 15–30. The effect was less pronounced for the other materials, but all columns showed higher TA concentrations under saturated conditions (Figure 3E).
The pH measurements in the effluent solutions are quite scattered, but nonetheless, there is an overall trend of decreasing pH with time for CO2-ambient conditions, which does not seem to be the case in saturated-CO2 conditions. In general, it appears that higher pH values occurred during the first half of all experiments for CO2-ambient conditions, but during the second half, the output pH values were slightly higher at saturated-CO2 conditions (Figures 3G,H). The effluent pH values reflect in large part the significant differences in the starting pH values of the input solutions.
The sum of cation equivalents released follows trends similar to those for TA (Figures 3B,D,F), in particular the sharp initial spike followed by an exponential decrease and then the attainment of a weak but negative linear relation (Figures 3B,D).
The relative loss of cations after the experimental period was generally higher under CO2-saturated conditions in all materials (Table 4). The ratio between the losses at ambient and saturated CO2 varies, depending on the cation and source material. On average, the loss ratio is 3.7 ± 3.8. The loss ratio of Mg is highest in the dunite column, where the presence of high CO2 concentrations leads to a nearly 16× higher flux than under ambient-CO2 conditions (Table 4). While the ratio is highest for Mg, it stands out that major cation losses from the dunite are dominated by monovalent K and Na. Also, basanite and basanite/laterite released relatively higher Na and K contents than Mg under saturated-CO2 conditions (Table 4). Losses from the laterite are markedly high for Mg, Ca, and K. The loess column lost predominantly divalent cations.
Table 4. Fraction of major cations (mole %) relative to the cation content of the fresh/unreacted column material released during the experiment phase. The resulting ratios between saturated and ambient CO2 treatment can be interpreted as the maximum achievable ratios in the field since watering volumes were chosen to simulate tropical humid high precipitation rates with enhanced complete drainage.
The relative release of elements is not proportional to their content in the untreated substrates (Figure 4). The major cation release of dunite is dominated by Na and K, even though their abundance is only 1% in the dunite. The basanite columns predominantly show Mg release under saturated CO2 conditions, but Ca, Na, and K contribute 35% under atmospheric conditions (Figure 4). Compared to the original material, laterite shows a pronounced release of K at ambient CO2 and of Ca at saturated CO2. The cation release from loess is proportionally larger for Mg at both CO2 concentrations, while Na and K releases are very low compared to the original composition.
Figure 4. Relative fractions of major cations [from left to right: Mg (green), Ca (light orange), Na (light red), K (red)] in the original, unreacted material (see XRF analyses, Table 1) and the relative cumulative fractions of each cation released over the 300-day reaction time under ambient and saturated-CO2 conditions. (A) Dunite, (B) Basanite, (C) Laterite, (D) Loess.
This series of experiments was designed to shed light on the effects of grain sizes (their distribution and specific surface area) on weathering rates and ion fluxes using dunite. Columns with smaller grain sizes resulted in higher Mg, TA, and effluent pH (Figure 5). The mixed granulometry results approximated those of the fine grain size rather than the coarse grains. Results for Mg (Figure 5B) and TA concentrations (Figure 5C) show similar trends, as was the case for dunite under ambient CO2 in the CO2 experimental series (Figures 3A,C). The effluent pH trends are very similar to TA and show a sharp initial spike, and then an exponential decrease to values of ~8 (Figure 5D). Si concentrations do not follow the generally observed trends of TA, Mg, and pH (Figure 5A). Dissolved silica concentrations are very low (<90 μmol L−1, Figure 5A) in comparison to the released magnesium (>1,500 μmol L−1). Moreover, Si release showed very similar trends for the three-grain sizes, with the mixed case showing three anomalous concentration excursions of unknown origin (Figure 5B). Si values generally fluctuate between 5 and 20 μmol L−1, with four outliers.
Figure 5. Temporal evolution of averaged dissolved silica [Si, (A)], magnesium [Mg, (B)], total alkalinity [TA, (C)] concentrations, and pH (D) in the outflow of the columns filled with dunite of different grain size distributions under ambient CO2 conditions. (E) Evolution of the Mg/Si molar ratio; here the y-axis minimum was chosen to represent the nominal ratio of Mg/Si in olivine, i.e., congruent dissolution (Mg2SiO4 → 2:1), Note that during the final period, values are far higher than 2, indicating a non-stoichiometric alteration process. All the data are differentiated by fine (green squares), coarse (red circles), and mixed (yellow triangles) grain sizes.
The Mg/Si molar ratios are always well above the theoretical value of 2 (Figure 5E), with values >100 at the beginning and an average long-term average of 47.8. The mechanism responsible for this non-stoichiometry can (in theory) result from either preferential Mg release, or preferential Si retention, and can only be elucidated with nanometer-scale solid-state surface analyses (e.g., Hellmann et al., 2012, 2021; Zandanel et al., 2022).
Annual CO2 sequestration based on TA fluxes
The initial period spanning the first few weeks showed substantially elevated TA concentrations (Figures 3, 5), which in turn corresponded to initially higher land surface-based release rates. Over significantly longer time periods (up to 1 year) the TA concentrations, and consequently the CO2 removal rates, decreased until a baseline characterized by a more constant behavior was reached, which can be considered to be representative of the first year (Figure 6). These lower, but constant weathering rates show that, under ambient-CO2 conditions, CO2 consumption is highest in the loess column and lowest in the laterite column, with 730 and 100 t CO2 km−2 a−1, respectively (Figure 6A). Dunite (all-grain sizes under ambient CO2) and basanite have similar consumption rates (around 220 t CO2 km−2 a−1). These values are comparable to the plotted literature values for ultramafic lithologies (Figure 6A). The mix of basanite and laterite has a much higher consumption rate (around 430 t CO2 km−2 a−1) compared to the basanite alone. If the initial intensive weathering period is considered representative of the first year, CO2 consumption rates are significantly higher in all ambient-CO2 experiments (by 20–250%). The CO2-saturated experiments show higher CO2 consumption rates (Figure 6B), compared to the ambient CO2 treatment by roughly a factor of four for the basanite, laterite, and basanite/laterite, but around 11 for the dunite column, excluding the initial flushing phase (Figure 6C).
Figure 6. The long-term annual CO2 consumption normalized to land surface area in t CO2 km−2 a−1 under ambient-CO2 (A) and saturated-CO2 conditions (B). The light blue bars indicate the additional CO2 consumption during the first year, as rates steeply increase during the first weeks of deployment. (C) The ratio between rates of saturated over ambient CO2 conditions is based on the expected annual consumption without the initial flushing phase. 1from Table 1 in Hartmann et al. (2009): 11.77, 13.97, and 2.05 t CO2 km−2 a−1, respectively; 2from Table 1 in Dessert et al. (2003); 3from Table 6 in Schopka et al. (2011).
Mineral saturation states
As a simplified indicator for potential secondary mineral precipitation, saturation indices (SI) were calculated. It should be kept in mind that these are thermodynamic calculations, and therefore may not correctly predict precipitation, given that kinetics may in fact be rate-controlling. The SI of calcite and amorphous silica throughout the experiments was calculated using the composition of the outflow solution from all columns (Figure 7). Under ambient conditions, the saturation indices of the effluent solutions suggest that the precipitation of carbonate minerals (e.g., calcite, Figures 7A–D) was thermodynamically favorable only for a short period at the beginning of the experiment. Beyond this initial period, the solutions were undersaturated with respect to all considered carbonate species (Supplementary Material 3). Under CO2-saturated conditions, water draining the columns was, after equilibration with the atmosphere, most of the time oversaturated or near saturation with respect to calcite for the basanite and basanite-laterite columns (Figures 7A,B). These SI indices change if it is assumed that the effluent water remained 100% CO2-saturated; in this scenario, the solutions were then consistently undersaturated with respect to calcite. When the dunite columns were considered with the same scenario (100% CO2-saturated effluent solution), calcite was also undersaturated, except over the first few weeks. Similar results were calculated for aragonite and dolomite (Supplementary Material 3-1). The hydrous magnesium mineral nesquehonite was always undersaturated, even in the dunite experiment (Supplementary Material 3-2). Magnesite may have precipitated throughout the experiment if equilibration with the atmosphere occurred.
Figure 7. Calculated saturation indices for calcite (A–D) and amorphous silica (E–H) based on the composition of the outflow solution of the column experiments. Values were calculated using PhreeqC. Ambient and saturated values in (A–C,E,F) are calculated based on the outflow solution, which was equilibrated for 24 h or longer with the laboratory atmosphere. In addition, saturation indices were calculated assuming the effluent solutions remained CO2-saturated.
All solutions were generally below saturation with respect to amorphous SiO2 under ambient-CO2 conditions, while at saturated-CO2 conditions values fluctuated around zero, indicating that precipitation could have occurred (Figures 7E,F).
Land surface-based CO2 consumption
The CO2 consumption rates shown in Figure 6 were calculated with a straightforward approach based on TA fluxes from the columns. Two perspectives on achievable CO2 sequestration rates can be drawn from our results. When scaled to 1 year, total CO2 uptake of nearly 1,000 t CO2 km2 a−1 was calculated for loess and dunite. However, this estimate is only valid for the first year of application, which includes a very elevated weathering signal contribution from the first 3 weeks. Further on in time, the systems approach an approximate dynamic steady-state state that leads to markedly lower overall CO2 uptake rates. The possible explanations for the decline in rates are discussed in Section 4.2. The dynamic steady-state in the later stages of the experiments yields significantly lower CO2 uptake rates, which for dunite and basanite are comparable to the highest reported literature values of catchments found for tropical areas with volcanic rocks. This order of magnitude equivalence in CO2 removal validates the applied environmental parameters in this study (temperature and precipitation simulating tropical conditions). Applying the observed rates to the globally available tropical agricultural areas [5.1 × 106 km2, Strefler et al. (2018a), yields a total annual CO2 uptake of 1.1 Gt CO2 a−1]. This is about 14% of the annual minimum required CDR to attain the 1.5° target estimate, stipulated in the Paris Agreement (United Nations, 2015; Strefler et al., 2018b), and indicates that inorganic carbon uptake by EW based on atmospheric CO2 levels in soils might contribute significantly to a portfolio of CDR methods. Moreover, in situ soil pCO2 is normally considerably higher than ambient atmospheric values, where estimates of up to 70,000 ppm in tropical soils are not unusual (Brook et al., 1983; Davidson and Trumbore, 2017; Hashimoto et al., 2017; Romero-Mujalli et al., 2018). Higher CO2 levels increase CO2 drawdown rates, by a factor of more than four under a scenario of 100% CO2-saturated conditions compared to ambient levels, as demonstrated in this study. Counterintuitively, the land surface-normalized CO2 sequestration rate of dunite is roughly equal to that of basanite under ambient-CO2 conditions (Figure 6A), and yet the ratio of CO2-saturated to ambient-CO2 uptake in the dunite columns is more than two times higher than for all other materials (11.3 vs. ~4.5). We hypothesize that the weathering reaction in the dunite column under ambient-CO2 conditions was limited by the available CO2.
To put into perspective the weathering rates that we present, it needs to be re-iterated that the experiments in this study were designed to simulate plausible maximum weathering rates, under specific conditions of (very) humid tropics, and an (unlikely) 100% drainage of water within the columns (water flowing through a bed of rock powder). Although the investigated effects of saturated-CO2 conditions are not realistic for natural soils, the results likely represent the upper limit for inorganic weathering, without considering the potential contributions of plants, microbes, and bioturbation. The CO2-enhancement factor for elevated realistic CO2 levels in soils compared to ambient conditions remains to be determined and might not be linear (c.f. discussion Section The CO2 effect and feedback constraints in Earth system model parameterizations).
Decrease of weathering rates with time
The surface area normalized weathering rates, here represented by TA (Figure 8), are in general lower than far-from-equilibrium rates using mixed-flow or batch reactors (cf., Palandri and Kharaka, 2004). Dissolution experiments using the latter types of reactors are generally characterized by very low solid to solution ratios, pre-treatment of grains to remove fines, uniform grain sizes, and run times commonly on the order of hundreds of hours to several tens of days. In contrast, our column experiments, which have a high rock-to-water volume ratio, provide rates that have maximum values below (ambient) or in the lower range of one standard deviation of fitted reported literature values [basanite columns, Figures 8C,D, comparison with fitted data from Bandstra and Brantley (2008)]. Moreover, in column experiments rates can decrease by up to two orders of magnitude within the first year (Figures 8E,F).
Figure 8. TA based weathering rates plotted against pH. Fitted data for literature weathering rates (gray area and dashed lines) from Bandstra and Brantley (2008); rate values were recalculated with a stoichiometric factor to get TA to release rates. Four moles TA are produced from one mole forsterite (Equation 9). Rates for basalt are reported as mol Si m-2 s-1; here, the conversion to TA release rates is based on a generalized basalt composition (from Navarre-Sitchler and Brantley, 2007), where the equivalent of one mole Si yields 2.5 moles TA. Gradient colors indicate elapsed time of the experiment. Note that in configurations e.g., there are more rate values available, yet no corresponding pH values were recorded. (A) Dunite, ambient CO2, (B) Dunite, saturated CO2, (C) Basanite, ambient CO2, (D) Basanite, saturated CO2, (E) Dunite, fine, (F) Dunite, mixed, (G) Dunite, coarse.
The decrease in reaction rates with time can be attributed to several processes. Initially, sample preparation conditions can lead to hardly quantifiable decreases in observed dissolution rates (Eggleston et al., 1989). Furthermore, a natural continuous decrease of available material may also influence the observed rate to a minor extent.
More importantly, dissolution rate declines may be attributed to the very rapid dissolution of fine particles, and/or surface sites having high surface energies (Holdren and Berner, 1979; Brantley, 2008). The crushed materials used in this study contain fines (no effort was made to remove them, e.g., see SEM images of dunite in Supplementary Material 4 and Table 2, column “share of smallest class”), and most likely also grain surfaces with sites of high energy density that were created during comminution.
Reaction rates reported in the literature based on mixed-flow and batch reactors, typically with agitators stirring a few grams of the minerals in a fluid medium, mitigate against transport-limited conditions of mineral dissolution and rather favor surface-controlled reactions (Awad et al., 2000; Martinez et al., 2014).
In a column reactor, the somewhat fast and preferential flow of water through channels, canals, or pore interconnections will be similar and also surface reaction controlled. Water with dissolved CO2 is moving through the pore system, and not all pores will exchange water at the same rate. We observed slight compaction of the material during the experimental runs, decreasing the pore space and the surface area available for solid-fluid contact, which also could have contributed to the reduction of dissolution rates.
In addition, the material will be in contact with stagnant fluids trapped within very tiny fluid pore spaces where it is unlikely that the reactions are transport controlled.
Past experiments on forsterite dissolution have observed non-stochiometric release of Mg/Si in the initial phase of dissolution, becoming stoichiometric after a few 100 h of reaction (Pokrovsky and Schott, 2000b). This was attributed to the formation of an amorphous silica-rich layer (ASL) on the mineral surface and has been reported for many other silicate minerals including albite, labradorite, diopside, wollastonite, and garnet (Wogelius and Walther, 1992; Pokrovsky and Schott, 2000b; Hellmann et al., 2012). At steady-state, the layer thickness varies and depends on the mineral, as well as the fluid dynamics, solution chemistry and pH, and temperature. Reported thicknesses range from <1 nm for laboratory altered forsterite [Zandanel et al., 2022; see also evidence for slightly thicker layers in Pokrovsky and Schott (2000a)] to 50 nm on naturally weathered K-feldspars and 200 nm on serpentine (Hellmann et al., 2012). The high Mg/Si ratio at the end of the dunite experiment (Figure 5E) suggests that Si is actively precipitating, even after 250 days. These results indicate that the system was still far from the steady-state, unlike in experiments with low mineral to water ratios, where steady-state and stoichiometric dissolution were often observed within a few 100 h (Wogelius and Walther, 1992; Pokrovsky and Schott, 2000b; Hellmann et al., 2012).
If it is assumed that the formation of the ASL is governed by a coupled interfacial dissolution-reprecipitation (CIDR) mechanism (Hellmann et al., 2003, 2012, 2021), and the ASL is porous, it should not impede dissolution of the silicates (Putnis and Putnis, 2007). However, if the ASL is non-porous and impermeable to fluid exchange between the dissolving silicate surface and the surrounding bulk fluid, the dissolution reaction may become kinetically hindered, or in a worst-case scenario, the reaction rate would evolve to such low values so as to be unmeasurable (Putnis and Putnis, 2007). Even if the ASL is porous and permeable, the distance between the exterior bulk fluid and the reacting silicate surface might result in a decline in rates over the course of the reaction (Pedrosa et al., 2016), at least for low pressure and temperature conditions (Pedrosa et al., 2017). Moreover, recent studies suggest that ASLs are heterogeneous and can structurally evolve with time, e.g., an increase in density. As a consequence, dissolution rates may decline further over time (Daval et al., 2011, 2017; Maher et al., 2016). Considering all of our experiments, after 1 year of reaction there was no evidence that the dissolution of the minerals stopped, even if some of the release curves are noisy and do not display clear trends (Figure 3).
Si can also precipitate with different morphologies. In an experiment by Lazaro et al. (2012), it was found that amorphous silica precipitated during intense olivine dissolution as nano-silica colloids, with a high SSA of 100–300 m2g−1. The observed nano-silica has primary particles of 10–25 nm that were agglomerated in clusters with an average pore size of ~20 nm. An increase in mesopores could (further) increase the water holding capacity of the soil (Schaller et al., 2020). This probably important aspect in the context of drought conditions needs further research to assess application consequences from the evolution of amorphous silica in quality and quantity, depending on the applied amount of silicate minerals, their grain sizes, and shapes.
Comparing pH-controlled dissolution experiments from the literature and the column-derived data from this study is difficult as the experiment conditions are different, yet they can provide a reference for the discussion. In the ambient-CO2 column experiments, dunite rates increase with pH (Figure 8) with values >7, as opposed to the trends shown by literature data fits (from Bandstra and Brantley, 2008). While in the literature, rates are measured at fixed pH in a system far from equilibrium, the column experiment is a system that reacts dynamically with water and CO2 input by proton production or consumption. Under ambient conditions, the capacity to consume protons declines toward the end of the experiment, indicated by lower pH and decreased weathering rates. Under CO2-saturated conditions, rates also decline with time, but the pH range is greater and not directly correlated with alkalinity production (Figures 8B,D). Weathering rates are still higher at the end of the experiment compared to rates under ambient conditions.
Effect of grain size (distribution)
It is often assumed that finer rock material grain sizes purposed for enhanced weathering would yield higher CO2 capture rates (e.g., Strefler et al., 2018a), simply due to the higher surface area available for reaction. This hypothesis is not supported by our experiment results for EW.
The experimental series designed to evaluate the effect of different grain sizes using the same dunite material shows higher TA fluxes from the fine grain size column, especially over the first 100 days. However, when normalized to the total grain surface areas, this result changes: while the total grain surface area of the fine column material is on average 8.4 ± 0.2 times higher than the surface area of the material in the coarse columns, the CO2 consumption of the fine material is only ~7% higher. Although the porosity is much higher in the fine grain size column (Table 1-3 in Supplementary Material), which would indicate a larger potential contact area for rock-water interaction, it is assumed that the characteristics of this larger pore space do not allow for rate-stimulating effects. The pore sizes may be much smaller, reflected in decreased hydraulic conductivity, qualitatively observed by longer residence times of the supernatant water (which percolated more slowly through the column). Moreover, significantly smaller pores likely did not contribute much to the transport of water through the column, which was more likely characterized by preferential fluid channels. In contrast, fluid flow through the coarse material with larger and inter-connected pores likely enabled a greater material-water contact area, leading to enhanced dissolution.
A comparable effect may explain the higher weathering fluxes of the laterite/basanite mix, compared to basanite alone, although the mix is only comprised of 50% basanite. The plain tropical soil, laterite, naturally dominated by very small grains (p80 < 30 μm) showed little hydraulic conductivity and long retention times of the supernatant water. Mixing with the coarse basanite yielded a grain size distribution more suitable for enhanced weathering, so the effect of hydrology, an increase of hydraulic conductivity, and enhanced grain-water contact time may lead to the enhanced observed weathering fluxes.
Both observed effects underline the importance of the hydrologic conditions within the material, which are controlled by the grain size distribution (as also hypothesized in Amann et al., 2020; c.f., modeled hydrological effects from EW in de Oliveira Garcia et al., 2020) and the nature of grain packing, which in the experiment was also influenced by manually filling the columns. Preferential flow along the column walls, or clogging and missing pore interconnections (as could be visually observed in the columns) may decrease apparent weathering rates, while a grain size distribution optimized for water residence times, with a balance between sufficient reaction time and oversaturation, can maximize the weathering effect and thus the CO2 uptake efficiency (de Oliveira Garcia et al., 2020).
For a real-world application, simply testing the rock flour material alone, without considering the soils it is applied to, will not be an effective strategy. Besides testing material mixed with the target soil, it seems imperative that hydrological effects should also be considered. Moreover, the effects of bioturbation, plant roots, or uptake of elements by plants (Akter and Akagi, 2010; Taylor et al., 2015; Krahl, 2020; Verbruggen et al., 2021) affect the transferability of the column experiment results to real-world soil-plant systems. Such effects on EW must be studied in more detail to enable enhanced predictions by models in the future, backed up of course, by experimental data.
Studies of major cation leaching are relevant for the assessment of different rock types to act as potential CO2 sinks via weathering. Our results (Table 4) suggest that the total dissolution of any of the evaluated materials under the conditions in this study would take about a century or longer under atmospheric pCO2 while being accelerated in the CO2-saturated environment. It is relevant for a CDR application scenario to understand the cation release rate and alkalinity production as a function of CO2 levels in the amended soils.
The relative release of major cations from the substrate diverges from its relative content in the parent mineral (Figure 4). Cation release is largely determined by the dissolution kinetics of dominant minerals. Mg dominates the cation flux from the dunite column, which is expected, as forsterite (Mg2SiO4) dominates the source rock. In the dunite experiment, notable amounts of K were released, which might be unexpected, given the general composition of the material. This behavior may be traced back to the occurrence of chabazite (Supplementary Material 5), a mineral from the zeolite group, which can contain K and has rather rapid weathering rates (Schofield et al., 2015; Moravec et al., 2021). K is relevant for achieving additional carbon sequestration via EW through biological processes (Goll et al., 2021), as K is an important nutrient to most plants. The release of Na and K from the basanite columns could be explained by the abundance of nepheline [(Na,K)AlSiO4], a feldspathoid mineral that commonly occurs in the silica-depleted magmas of the Eifel (Mertes and Schmincke, 1985). The dissolution rate of nepheline is significantly higher than other Na-dominated rock-forming minerals, i.e., albite (Lasaga et al., 1994; Franke, 2009). The high dissolution rate indicates that rocks rich in this mineral may be a preferable source for EW, even if a monovalent ion is associated with ~50% less CO2 drawdown compared to a divalent cation.
Target materials for EW application should be tested for available water (rainfall) and CO2 concentration gradients in the soil, as these two external parameters predominantly govern CO2 uptake. However, since K, but also Mg and Ca, are taken up by plants, their removal from the soil-water system is likely to decrease the instantaneous CO2 sequestration rate because plant uptake of cations is counterbalanced by the release of protons (Britto and Kronzucker, 2008), and cation concentrations in the water phase leaving the system should be lower due to their temporal storage in plants. The quantity of these elements that remain in agricultural fields after crop/tree harvest will affect the final long-term field-based CO2 removal rate. At the same time, given sufficient time, the leftover organic matter could re-release cations. These processes call for further investigations to provide a better picture of full inorganic CO2 sequestration through EW. To assess the full EW potential, inorganic and organic CO2 sequestration should be studied in parallel, focusing additionally on the effects of P release, an important micro-nutrient for plants that are commonly used in fertilizers, and which is abundant in the basanite (Table 1) but was not the focus of this study.
Riverine cation elemental ratios in watersheds (e.g., Ca/Na) can be used to distinguish between carbonate and silicate weathering (Gaillardet et al., 1999). Results here suggest that EW application with high amounts of nepheline might change Ca/Na ratios due to the non-stoichiometry of the alteration process, which potentially affects the interpretation of weathering sources, and also can alter the geochemical baselines of rivers (Hartmann et al., 2007). Under a large-scale EW deployment scenario, the influence of altered geochemical fluxes in natural systems on scientific interpretations has to be considered as a side effect not addressed so far.
The CO2 effect and feedback constraints in earth system model parameterizations
The CO2 feedback on weathering rates and CO2 consumption in Earth system or planetary models is often represented by the factor (Walker et al., 1981). In an abiotic world, as in our column experiments, this feedback would represent the direct effect of CO2 dissolution and consumption in water, in addition to hydrological feedback, as the weathering flux is largely controlled by water throughflow, given a fixed lithological composition. In general, it is not possible to distinguish the proportions each of these effects contribute to the feedback. The exponent n is often modulated to fit an assumed feedback-strength range (Uchikawa and Zeebe, 2008; Lehmer et al., 2020). We can compare our results with the resulting feedback strength for n = 0.1 to n = 0.5 using a pCO2 of 400 ppm as a baseline, instead of the pre-industrial value. This comparison would only consider the CO2 effect and no further hydrological feedback due to elevated CO2 in the atmosphere. Solving the above equation for n with α = 4.5, representative of the results in this experiment (Figure 6C) leads to n = 0.2, a value that matches the assumed default conditions in the simulations of Uchikawa and Zeebe (2008). This value might be representative of conditions on Mars or other planets with a high pCO2 atmosphere. With n = 0.3 the feedback is comparable to results from the dunite column, where it must be assumed that observed weathering fluxes are diffusion-limited (see above). The case of n = 0.5 would result in an even higher CO2-enhancement factor of 50. This would imply that in a basalt-dominated world, the additional hydrological effect would be stronger by a factor of 10. As the early, primitive surfaces of Earth or Mars were more likely to be comparable to a basaltic than to an olivine-dominated surface, this comparison invites a discussion about the combined effects of CO2 and hydrology in a planetary system, and how they can be best represented. If the CO2 feedback effect due to EW is modeled at the planetary scale, one should distinguish CO2 effects on mineral weathering and on hydrology separately.
Results suggest that n values in the lower range of values considered are more realistic. As only two CO2 conditions were tested, a simple scaling law would suggest that for humid tropical areas a doubling to tripling in weathering rates can be assumed for elevated (>70,000 ppm) CO2 conditions compared to ambient CO2 conditions (Table 5).
Table 5. Weathering rate factor , normalized to the baseline condition of 400 ppm CO2 as applied in common Earth system models to simulate weathering rate changes with increased pCO2, depending on the choice of exponent n.
Enhancement of weathering by rhizosphere processes
As the column experiment was primarily abiotic, the influence of biological processes on weathering rates was not explicitly considered, although microbial components may have exerted a minor influence on the columns filled with soils, as these were not sterilized. Processes in the rhizosphere can speed up dissolution processes and therefore increase EW potentials (Krahl, 2020; Verbruggen et al., 2021). Beerling et al. (2020) deduced a biotic weathering enhancement factor normalized to net primary productivity in the range of 1–7 (from work of Akter and Akagi, 2005, 2010; Quirk et al., 2012, 2014). This factor incorporates the effects of rhizosphere biologically accelerated weathering. Given the overall scarcity of data, but the relevant influence of plants and the rhizosphere, these processes need to be considered to determine realistic constraints on the boundaries of EW for diverse environmental settings in the future.
Precipitation of carbonates
The secondary precipitation of carbonates as the final long-term step in the weathering process can release half of the captured CO2 counterbalanced by Ca and to some proportion of Mg. Thus, carbonate precipitation should be monitored to evaluate the overall carbon sequestration potential of an agricultural system. However, carbonate mineral saturation indices indicate that carbonate formation only occurred during the initial phases of our experiments, since long-term oversaturation was not observed. Consequently, the precipitation of carbonates may not be expected in agricultural soils under long-term treatment with rock material, given the humid, tropical conditions. Precipitation downstream of the application area might occur after water is re-equilibrating with the atmosphere (Figure 7). See further discussion on efficacy of enhancing carbonate weathering for carbon dioxide sequestration in Knapp and Tipper (2022).
Our counter-intuitive finding that the smallest grain sizes may not lead to the fastest CO2 uptake (Figure 6) calls for a more detailed investigation of the influence of grain size distributions and packing density of soils in individual application cases. Larger target grain sizes will decrease the energy requirements for comminution. This will lower costs and the potential CO2 penalty from fossil energy use.
The pronounced alkalinity production during the initial phases of the experiments suggests that a modified application scheme might be considered: rock flour replenishment at shorter intervals (e.g., annually) may be preferable over a one-time application of a larger amount of material to achieve higher CO2 drawdown rates. This, however, needs to be tested with further experiments, as biological processes were not simulated in this study.
We have also shown that the individual mineral compositions making up applied rock material is important for the evaluation of CO2 uptake speeds and side effects. A site-specific material selection process, going further than a simple search for local “basalt,” can optimize dissolution rates (and thus CO2 uptake potential), as well as plant nutrient potential. Looking beyond basaltic material, loess can sustain rather high dissolution rates over millennia (Börker et al., 2020), which also showed high dissolution rates in the experiments. As such, it may even be a well-suited material to produce alkalinity—and remove CO2–in the context of CDR as it is often regionally available, easily extractable, and fertile. Material-specific aging processes and the formation of secondary precipitates will be issues that also need to be understood in more detail.
Overall, future models and projections for EW-related CDR should consider the dissolution rates obtained from column experiments under 100% CO2, as they will represent an upper boundary in an inorganic system. Further studies should assess the soil pCO2 gradients, simulating conditions in soils of different climates and environments with up to more than 100,000 ppm. Biologic processes should further accelerate the observed rates, which demands tailored experiments, distinguishing biological and inorganic geochemical processes, to improve EW-CDR-potential assessments.
A further observation from the experiments, the promising utilization of electrical conductivity as a predictor for TA production, a proxy for the CO2 sequestration effect, is discussed in a complementary publication (Amann and Hartmann, 2022).
Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.
JH, TA, and AM designed the study. AM ran the first experiment comparing grain sizes and contributed to the discussion. TA wrote the manuscript with input from JH, RH, EP, and AM. EP contributed to the discussion. RH provided input throughout the manuscript and especially to section decrease of weathering rates with time. All authors approved the submitted version.
This work was funded by Germany's Excellence Strategy—EXC 2037 Climate, Climatic Change, and Society—project number 390683824, contribution to the Center for Earth System Research and Sustainability (CEN) of the Universität Hamburg.
We acknowledge the help of Tom Jäppinen, Peggy Bartsch, Rhiannon Breider, Mark Brosell, Marvin Keitzel, Eric Marques, and Walid Karimi for valuable contributions from the wet lab, and Sebastian Lindhorst for providing granulometric analyses (all from the Institute for Geology, Universität Hamburg). We thank Stephan Jung and Joachim Ludwig (from the Institute for Mineralogy and Petrography, Universität Hamburg) for contributing to the XRF and XRD analyses. Discussions were further stimulated by work with the Carbon Drawdown Initiative (Project Carbdown).
Conflict of interest
The 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.
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.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fclim.2022.929268/full#supplementary-material
Amann, T., Hartmann, J., Struyf, E., de Oliveira Garcia, W., Fischer, E. K., Janssens, I., et al. (2020). Enhanced weathering and related element fluxes – a cropland mesocosm approach. Biogeosciences 17, 103–119. doi: 10.5194/bg-17-103-2020
Awad, A., Koster van Groos, A. F., and Guggenheim, S. (2000). Forsteritic olivine: effect of crystallographic direction on dissolution kinetics. Geochim. Cosmochim. Acta 64, 1765–1772. doi: 10.1016/S0016-7037(99)00442-1
Ball, J. W., and Nordstrom, D. K. (1991). “User's manual for WATEQ4F, with revised thermodynamic data base and text cases for calculating speciation of major, trace, and redox elements in natural waters,” in Open-File Report. Version 2.0 (Menlo Park, CA). doi: 10.3133/ofr91183
Bandstra, J. Z., and Brantley, S. L. (2008). “Data fitting techniques with applications to mineral dissolution kinetics,” in Kinetics of Water-Rock Interaction, eds S. L. Brantley, J. D. Kubicki, and A. F. White (New York, NY: Springer), 211–257. doi: 10.1007/978-0-387-73563-4_6
Beerling, D. J., Kantzas, E. P., Lomas, M. R., Wade, P., Eufrasio, R. M., Renforth, P., et al. (2020). Potential for large-scale CO2 removal via enhanced rock weathering with croplands. Nature 583, 242–248. doi: 10.1038/s41586-020-2448-9
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
Beyer, E. E. (2006). Transformation of archaean lithospheric mantle by refertilization: evidence from exposed peridotites in the western Gneiss Region, Norway. J. Petrol. 47, 1611–1636. doi: 10.1093/petrology/egl022
Börker, J., Hartmann, J., Amann, T., Romero-Mujalli, G., Moosdorf, N., and Jenkins, C. (2020). Chemical weathering of loess and its contribution to global alkalinity fluxes to the coastal zone during the last glacial maximum, mid-holocene, and present. Geochem. Geophys. Geosyst. 21, e2020GC008922. doi: 10.1029/2020GC008922
Brantley, S. L. (2008). “Kinetics of mineral dissolution,” in Kinetics of Water-Rock Interaction, eds S. L. Brantley, J. D. Kubicki, and A. F. White (New York, NY: Springer), 151–210. doi: 10.1007/978-0-387-73563-4_5
Cross, W. I., Joseph, P., Pirsson, L. V., and Washington, H. S. (1902). A quantitative chemico-mineralogical classification and nomenclature of igneous rocks. J. Geol. 10, 555–690. doi: 10.1086/621030
Daval, D., Bernard, S., Rémusat, L., Wild, B., Guyot, F., Micha, J. S., et al. (2017). Dynamics of altered surface layer formation on dissolving silicates. Geochim. Cosmochim. Acta 209, 51–69. doi: 10.1016/j.gca.2017.04.010
Daval, D., Sissmann, O., Menguy, N., Saldi, G. D., Guyot, F., Martinez, I., et al. (2011). Influence of amorphous silica layer formation on the dissolution rate of olivine at 90 °C and elevated pCO2. Chem. Geol. 284, 193–209. doi: 10.1016/j.chemgeo.2011.02.021
de Oliveira Garcia, W., Amann, T., Hartmann, J., Karstens, K., Popp, A., Boysen, L. R., et al. (2020). Impacts of enhanced weathering on biomass production for negative emission technologies and soil hydrology. Biogeosciences 17, 2107–2133. doi: 10.5194/bg-17-2107-2020
Dessert, C., Dupré, B., Gaillardet, J., François, L. M., and Allègre, C. J. (2003). Basalt weathering laws and the impact of basalt weathering on the global carbon cycle. Chem. Geol. 202, 257–273. doi: 10.1016/j.chemgeo.2002.10.001
Dickson, A. G. (1981). An exact definition of total alkalinity and a procedure for the estimation of alkalinity and total inorganic carbon from titration data. Deep Sea Res. Part A Oceanogr. Res. Pap. 28, 609–623. doi: 10.1016/0198-0149(81)90121-7
Dietzen, C., Harrison, R., and Michelsen-Correa, S. (2018). Effectiveness of enhanced mineral weathering as a carbon sequestration tool and alternative to agricultural lime: an incubation experiment. Int. J. Greenh. Gas Con. 74, 251–258. doi: 10.1016/j.ijggc.2018.05.007
Ebelmen, J. J. (1845). Sur les Produits de la Décomposition des Espèces Minérales de la Famille des Silicates. Paris: Carilian-G.ury et Vor Dalmont, Libraires des Corps Royaux des Ponts et Chausees et dus Mines.
Eggleston, C. M., Hochella, M. F., and George, P. A. (1989). Sample preparation and aging effects on the dissolution rate and surface composition of diopside. Geochim. Cosmochim. Acta 53, 797–804. doi: 10.1016/0016-7037(89)90026-4
Gaillardet, J., Dupré, B., Louvat, P., and Allègre, C. J. (1999). Global silicate weathering and CO2 consumption rates deduced from the chemistry of large rivers. Chem. Geol. 159, 3–30. doi: 10.1016/S0009-2541(99)00031-5
Giambelluca, T. W., Chen, Q., Frazier, A. G., Price, J. P., Chen, Y.-L., Chu, P.-S., et al. (2013). Online rainfall atlas of Hawai‘i. Bull. Am. Meteorol. Soc. 94, 313–316. doi: 10.1175/BAMS-D-11-00228.1
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
Haque, F., Santos, R. M., and Chiang, Y. W. (2020). CO2 sequestration by wollastonite-amended agricultural soils – an Ontario field study. Int. J. Greenh. Gas Con. 97, 103017. doi: 10.1016/j.ijggc.2020.103017
Hartmann, J., Jansen, N., Dürr, H. H., Kempe, S., and Köhler, P. (2009). Global CO2-consumption by chemical weathering: what is the contribution of highly active weathering regions? Glob. Planet. Change 69, 185–194. doi: 10.1016/j.gloplacha.2009.07.007
Hartmann, J., Jansen, N., Kempe, S., and Dürr, H. (2007). Geochemistry of the river rhine and the upper danube: recent trends and lithological influence on baselines. J. Environ. Sci. Sustain. Soc. 1, 39–46. doi: 10.3107/jesss.1.39
Hartmann, J., and Kempe, S. (2008). What is the maximum potential for CO(2) 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
Hashimoto, S., Tanaka, N., Suzuki, M., Inoue, A., Takizawa, H., Kosaka, I., et al. (2017). Soil respiration and soil CO2 concentration in a tropical forest, Thailand. J. For. Res. 9, 75–79. doi: 10.1007/s10310-003-0046-y
Hellmann, R., Penisson, J. M., Hervig, R. L., Thomassin, J. H., and Abrioux, M. F. (2003). An EFTEM/HRTEM high-resolution study of the near surface of labradorite feldspar altered at acid pH: evidence for interfacial dissolution-reprecipitation. Phys. Chem. Min. 30, 192–197. doi: 10.1007/s00269-003-0308-4
Hellmann, R., Wirth, R., Daval, D., Barnes, J.-P., Penisson, J.-M., Tisserand, D., et al. (2012). Unifying natural and laboratory chemical weathering with interfacial dissolution–reprecipitation: a study based on the nanometer-scale chemistry of fluid–silicate interfaces. Chem. Geol. 294–295, 203–216. doi: 10.1016/j.chemgeo.2011.12.002
Hellmann, R., Zhai, Y., Robin, E., Findling, N., Mayanna, S., Wirth, R., et al. (2021). The hydrothermal alkaline alteration of potassium feldspar: a nanometer-scale investigation of the orthoclase interface. Chem. Geol. 569. doi: 10.1016/j.chemgeo.2021.120133
Hövelmann, J., Austrheim, H., and Jamtveit, B. (2012). Microstructure and porosity evolution during experimental carbonation of a natural peridotite. Chem. Geol. 334, 254–265. doi: 10.1016/j.chemgeo.2012.10.025
Kelland, M. E., Wade, P. W., Lewis, A. L., Taylor, L. L., Sarkar, B., Andrews, M. G., et al. (2020). Increased yield and CO2 sequestration potential with the C4 cereal Sorghum bicolor cultivated in basaltic rock dust-amended agricultural soil. Glob. Chang. Biol. 26, 3658–3676. doi: 10.1111/gcb.15089
Klasen, N., Fischer, P., Lehmkuhl, F., and Hilgers, A. (2015). Luminescence dating of loess deposits from the Remagen-Schwalbenberg site, Western Germany. Geochronometria 42. doi: 10.1515/geochr-2015-0008
Krahl, L. L. (2020). Mineral Formation and Element Release From Aluminosilicate Rocks Promoted by Maize Rhizosphere (Ph.D). Available online at: https://repositorio.unb.br/bitstream/10482/38678/1/2020_LuiseLotticiKrahl.pdf
Lasaga, A. C., Soler, J. M., Ganor, J., Burch, T. E., and Nagy, K. L. (1994). Chemical weathering rate laws and global geochemical cycles. Geochim. Cosmochim. Acta 58, 2361–2386. doi: 10.1016/0016-7037(94)90016-7
Lazaro, A., Brouwers, H. J. H., Quercia, G., and Geus, J. W. (2012). The properties of amorphous nano-silica synthesized by the dissolution of olivine. Chem. Eng. J. 211–212, 112–121. doi: 10.1016/j.cej.2012.09.042
Le Bas, M. J., Maitre, R. W. L., Streckeisen, A., and Zanettin, B. (1986). A chemical classification of volcanic rocks based on the total alkali-silica diagram. J. Petrol. 27, 745–750. doi: 10.1093/petrology/27.3.745
Lechler, P. J., and Desilets, M. O. (1987). A review of the use of loss on ignition as a measurement of total volatiles in whole-rock analysis. Chem. Geol. 63, 341–344. doi: 10.1016/0009-2541(87)90171-9
Lehmer, O. R., Catling, D. C., and Krissansen-Totton, J. (2020). Carbonate-silicate cycle predictions of earth-like planetary climates and testing the habitable zone concept. Nat. Commun. 11, 6153. doi: 10.1038/s41467-020-19896-2
Leonardos, O. H., Fyfe, W. S., and Kronberg, B. I. (1987). The use of ground rocks in laterite systems: an improvement to the use of conventional soluble fertilizers? Chem. Geol. 60, 361–370. doi: 10.1016/0009-2541(87)90143-4
Maher, K., Johnson, N. C., Jackson, A., Lammers, L. N., Torchinsky, A. B., Weaver, K. L., et al. (2016). A spatially resolved surface kinetic model for forsterite dissolution. Geochim. Cosmochim. Acta 174, 313–334. doi: 10.1016/j.gca.2015.11.019
Martinez, R. E., Weber, S., and Bucher, K. (2014). Quantifying the kinetics of olivine dissolution in partially closed and closed batch reactor systems. Chem. Geol. 367, 1–12. doi: 10.1016/j.chemgeo.2013.12.017
Moravec, B. G., Keifer, V., Root, R. A., White, A. M., Wang, Y., Olshansky, Y., et al. (2021). Experimental weathering of a volcaniclastic critical zone profile: Key role of colloidal constituents in aqueous geochemical response. Chem. Geol. 559, 119886. doi: 10.1016/j.chemgeo.2020.119886
Palandri, J. L., and Kharaka, Y. K. (2004). A Compilation of Rate Parameters of Water-Mineral Interaction Kinetics for Application to Geochemical Modeling. Menlo Park, CA: U.S. Geological Survey. doi: 10.3133/ofr20041068
Parkhurst, D. L., and Appelo, C. A. J. (2013). “Description of input and examples for PHREEQC version 3—a computer program for speciation, batch-reaction, one-dimensional transport, and inverse geochemical calculations,” in U.S. Geological Survey Techniques and Methods (Denver, CO). doi: 10.3133/tm6A43
Pedrosa, E. T., Putnis, C. V., Renard, F., Burgos-Cara, A., Laurich, B., and Putnis, A. (2016). Porosity generated during the fluid-mediated replacement of calcite by fluorite. Cryst. Eng. Comm. 18, 6867–6874. doi: 10.1039/C6CE01150K
Pokrovsky, O. S., and Schott, J. (2000a). Forsterite surface composition in aqueous solutions: a combined potentiometric, electrokinetic, and spectroscopic approach. Geochim. Cosmochim. Acta 64, 3299–3312. doi: 10.1016/S0016-7037(00)00435-X
Quirk, J., Andrews, M. Y., Leake, J. R., Banwart, S. A., and Beerling, D. J. (2014). Ectomycorrhizal fungi and past high CO2 atmospheres enhance mineral weathering through increased below-ground carbon-energy fluxes. Biol. Lett. 10, 20140375. doi: 10.1098/rsbl.2014.0375
Quirk, J., Beerling, D. J., Banwart, S. A., Kakonyi, G., Romero-Gonzalez, M. E., and Leake, J. R. (2012). Evolution of trees and mycorrhizal fungi intensifies silicate mineral weathering. Biol. Lett. 8, 1006–1011. doi: 10.1098/rsbl.2012.0503
Renforth, P., Pogge von Strandmann, P. A. E., and Henderson, G. M. (2015). The dissolution of olivine added to soil: implications for enhanced weathering. Appl. Geochem. 61, 109–118. doi: 10.1016/j.apgeochem.2015.05.016
Romero-Mujalli, G., Hartmann, J., Börker, J., Gaillardet, J., and Calmels, D. (2018). Ecosystem controlled soil-rock pCO2 and carbonate weathering – constraints by temperature and soil water content. Chem. Geol. 527, 118634. doi: 10.1016/j.chemgeo.2018.01.030
Schaller, J., Cramer, A., Carminati, A., and Zarebanadkouki, M. (2020). Biogenic amorphous silica as main driver for plant available water in soils. Sci. Rep. 10, 2424. doi: 10.1038/s41598-020-59437-x
Schofield, R. E., Hausrath, E. M., and Gainey, S. R. (2015). “Zeolite weathering in laboratory and natural settings, and implications for mars,” in 46th Lunar and Planetary Science Conference (The Woodlands, TX).
Schopka, H. H., Derry, L. A., and Arcilla, C. A. (2011). Chemical weathering, river geochemistry and atmospheric carbon fluxes from volcanic and ultramafic regions on Luzon Island, the Philippines. Geochim. Cosmochim. Acta 75, 978–1002. doi: 10.1016/j.gca.2010.11.014
Soil Survey Staff (1999). “Soil taxonomy: a basic system of soil classification for making and interpreting soil surveys,” in Natural Resources Conservation Service. U.S. Department of Agriculture Handbook 436, 2nd edn (Soil Survey Staff). Available online at: https://www.nrcs.usda.gov/wps/portal/nrcs/main/soils/survey/class/taxonomy/
Strasberg, D., Rouget, M., Richardson, D. M., Baret, S., Dupont, J., and Cowling, R. M. (2005). An Assessment of habitat diversity and transformation on la réunion island (Mascarene Islands, Indian Ocean) as a basis for identifying broad-scale conservation priorities. Biodivers. Conserv. 14, 3015–3032. doi: 10.1007/s10531-004-0258-2
Streckeisen, A. (1980). Classification and nomenclature of volcanic rocks, lamprophyres, carbonatites and melilitic rocks IUGS Subcommission on the systematics of igneous rocks. Geol. Rundschau 69, 194–207. doi: 10.1007/BF01869032
Strefler, J., Amann, T., Bauer, N., Kriegler, E., and Hartmann, J. (2018a). Potential and costs of carbon dioxide removal by enhanced weathering of rocks. Environ. Res. Lett. 13, 034010. doi: 10.1088/1748-9326/aaa9c4
Strefler, J., Bauer, N., Kriegler, E., Popp, A., Giannousakis, A., and Edenhofer, O. (2018b). Between scylla and charybdis: delayed mitigation narrows the passage between large-scale CDR and high costs. Environ. Res. Lett. 13, 044015. doi: 10.1088/1748-9326/aab2ba
Taylor, L. L., Quirk, J., Thorley, R. M. S., Kharecha, P. A., Hansen, J., Ridgwell, A., et al. (2015). Enhanced weathering strategies for stabilizing climate and averting ocean acidification. Nat. Clim. Chang. 6, 402–406. doi: 10.1038/nclimate2882
Teitler, Y., Le Hir, G., Fluteau, F., Philippot, P., and Donnadieu, Y. (2014). Investigating the Paleoproterozoic glaciations with 3-D climate modeling. Earth Planet. Sci. Lett. 395, 71–80. doi: 10.1016/j.epsl.2014.03.044
ten Berge, H. F. M., van der Meer, H. G., Steenhuizen, J. W., Goedhart, P. W., Knops, P., and Verhagen, J. (2012). Olivine weathering in soil, and its effects on growth and nutrient uptake in ryegrass Lolium perenne L.: a pot experiment. PLoS ONE 7, e42098. doi: 10.1371/journal.pone.0042098
Vienne, A., Ibanez, S. P., Portillo-Estrada, M., Hartmann, J., Ijiehon, S., Wade, P., et al. (2022). Enhanced weathering using basalt rock powder: carbon sequestration, co-benefits and risks in a mesocosm study with Solanum tuberosum. Front. Clim. 4, 869456. doi: 10.3389/fclim.2022.869456
Vitens (2021). PhreeqPython. Vitens. Available online at: https://github.com/Vitens/phreeqpython).
Walker, J. C. G., Hays, P. B., and Kasting, J. F. (1981). A negative feedback mechanism for the long-term stabilization of earth's surface temperature. J. Geophys. Res. 86, 9776–9782. doi: 10.1029/JC086iC10p09776
Keywords: climate change, negative emissions, carbon dioxide removal, enhanced weathering, column experiment
Citation: Amann T, Hartmann J, Hellmann R, Pedrosa ET and Malik A (2022) Enhanced weathering potentials—the role of in situ CO2 and grain size distribution. Front. Clim. 4:929268. doi: 10.3389/fclim.2022.929268
Received: 26 April 2022; Accepted: 27 June 2022;
Published: 22 July 2022.
Edited by:Fatima Haque, National Taiwan University, Taiwan
Reviewed by:Davide Ciceri, Agroplantae, United States
Debanjan Chandra, Indian Institute of Technology Bombay, India
Copyright © 2022 Amann, Hartmann, Hellmann, Pedrosa and Malik. 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: Thorben Amann, email@example.com
†ORCID: Thorben Amann https://orcid.org/0000-0001-9347-0615
Jens Hartmann https://orcid.org/0000-0003-1878-9321
Roland Hellmann https://orcid.org/0000-0002-9540-6817
Elisabete Trindade Pedrosa https://orcid.org/0000-0003-3148-5349
Aman Malik https://orcid.org/0000-0002-7310-8448