Submitted as: model description paper
24 Aug 2023
Submitted as: model description paper |  | 24 Aug 2023
Status: this preprint is currently under review for the journal GMD.

In-silico calculation of soil pH by SCEPTER v1.0

Yoshiki Kanzaki, Isabella Chiaravalloti, Shuang Zhang, Noah J. Planavsky, and Christopher T. Reinhard

Abstract. One of the soil properties most commonly measured to describe agronomic and biogeochemical conditions of soils is “soil pH”. Soil pH measures the concentration of exchangeable H+ that resides in bulk soil samples taken from the field, through aqueous H+ measurements of extractants (e.g., deionized water or electrolyte solutions) added to dried bulk soil samples in the laboratory. Therefore, “soil pH” differs from “porewater pH”, the latter of which we define here as an in-situ measure of porewater H+ concentration in soil/weathering profiles. The difference between the two pH measurements is often not fully known for a given system but could lead to a misunderstanding of soil conditions if the two measurements are directly compared. Agricultural soils are one of the targeted loci for application of the “Enhanced Rock Weathering” (ERW), a technique aimed at counteracting increasing anthropogenic carbon dioxide from burning fossil fuels, and an increase in pH is thought to be one of key advantages of ERW as this can mitigating soil acidification and secure crop yields. As a result, fully evaluating the biogeochemical and agronomic consequences of ERW approaches requires accurate simulation of both soil pH (pHs) and porewater pH (pHpw). This paper presents an updated version of the reactive transport code SCEPTER (Soil Cycles of Elements simulator for Predicting TERrestrial regulation of greenhouse gases), which enables simulation of bulk soil pH measurement in the laboratory in addition to porewater pH as measured in the field along with a more comprehensive representation of cation exchange with solid-phase constituents of bulk soil. We first describe the implementation of cation exchange in the SCEPTER model, then introduce conceptual modelling frameworks enabling the calculation of bulk pHs. The validity of the model is examined through comparison of model results with soil pH measurements from mesocosm experiments of maize production with crushed basalt amendments. Finally, illustrative example simulations are shown demonstrating that a difference between pHs and pHpw can lead to significantly different estimates of carbon capture by ERW for a given targeted pH in cropland systems.

Yoshiki Kanzaki et al.

Status: open (until 28 Oct 2023)

Comment types: AC – author | RC – referee | CC – community | EC – editor | CEC – chief editor | : Report abuse

Yoshiki Kanzaki et al.

Yoshiki Kanzaki et al.


Total article views: 223 (including HTML, PDF, and XML)
HTML PDF XML Total BibTeX EndNote
163 55 5 223 3 2
  • HTML: 163
  • PDF: 55
  • XML: 5
  • Total: 223
  • BibTeX: 3
  • EndNote: 2
Views and downloads (calculated since 24 Aug 2023)
Cumulative views and downloads (calculated since 24 Aug 2023)

Viewed (geographical distribution)

Total article views: 209 (including HTML, PDF, and XML) Thereof 209 with geography defined and 0 with unknown origin.
Country # Views %
  • 1
Latest update: 23 Sep 2023
Short summary
Soil pH is one of the most commonly measured agronomical and biogeochemical indices, mostly reflecting exchangeable acidity. Explicit simulation of both porewater pH and bulk soil pH is thus crucial to accurate evaluation of alkalinity required to counteract soil acidification and resulting capture of anthropogenic carbon dioxide through the Enhanced Rock Weathering technique. This has been enabled by the updated reactive-transport SCEPTER code and newly developed framework to simulate soil pH.