<p><strong>Abstract.</strong> The PROFILE model, now incorporated in the ForSAFE model can accurately reproduce the chemical and mineralogical evolution of the soil unsaturated zone. However, in deeper soil layers and in groundwater systems, it appears to overestimate weathering rates. This overestimation has been corrected by improving the kinetic expression describing mineral dissolution by adding or upgrading <q>breaking functions</q>. The base cation and aluminium brakes have been strengthened, and an additional silicate brake has been developed, improving the ability to describe mineral-water reactions in deeper soils. These brakes are developed from a molecular-level model of the dissolution mechanisms. Equations, parameters and constants describing mineral dissolution kinetics have now been obtained for 102 different minerals from 12 major structural groups, comprising all types of minerals encountered in most soils. The PROFILE and ForSAFE weathering sub-model was extended to cover two-dimensional catchments, both in the vertical and the horizontal direction, including the hydrology. Comparisons between this improved model and field observations is available in Erlandsson Lampa et al. (2019, This special issue). The results showed that the incorporation of a braking effect of silica concentrations was necessary and helps obtain more accurate descriptions of soil evolution rates at greater depths and within the saturated zone.</p>

Abstract. Soil and water acidification was internationally recognised as a severe environmental problem in the late 1960s. The interest in establishing “critical loads” led to a peak in weathering research in the 1980s and 1990s, since base cation weathering is the long-term counterbalance to acidification pressure. Assessments of weathering rates and associated uncertainties have recently become an area of renewed research interest, this time due to demand for forest residues to provide renewable bioenergy. Increased demand for forest fuels increases the risk of depleting the soils of base cations produced in situ by weathering. This is the background to the research programme Quantifying Weathering Rates for Sustainable Forestry (QWARTS), which ran from 2012 to 2019. The programme involved research groups working at different scales, from laboratory experiments to modelling. The aims of this study were to (1) investigate the variation in published weathering rates of base cations from different approaches in Sweden, with consideration of the key uncertainties for each method; (2) assess the robustness of the results in relation to sustainable forestry; and (3) discuss the results in relation to new insights from the QWARTS programme and propose ways to further reduce uncertainties. In the study we found that the variation in estimated weathering rates at single-site level was large, but still most sites could be placed reliably in broader classes of weathering rates. At the regional level, the results from the different approaches were in general agreement. Comparisons with base cation losses after stem-only and whole-tree harvesting showed sites where whole-tree harvesting was clearly not sustainable and other sites where variation in weathering rates from different approaches obscured the overall balance. Clear imbalances appeared mainly after whole-tree harvesting in spruce forests in southern and central Sweden. Based on the research findings in the QWARTS programme, it was concluded that the PROFILE/ForSAFE family of models provides the most important fundamental understanding of the contribution of weathering to long-term availability of base cations to support forest growth. However, these approaches should be continually assessed against other approaches. Uncertainties in the model approaches can be further reduced, mainly by finding ways to reduce uncertainties in input data on soil texture and associated hydrological parameters but also by developing the models, e.g. to better represent biological feedbacks under the influence of climate change.

Soil mineral weathering is one of the major sources of base cations (BC), which play a dual role for a forest ecosystem; they function both as plant nutrients, and for buffering against acidification of catchment runoff. On a long-term basis, the soil weathering rates will determine the highest sustainable forest productivity without causing acidification. It is believed that the hydrologic residence time plays a key role in determining weathering rates on a landscape scale.The weathering model PROFILE has been used for almost 30 years to calculate weathering rates in the rooting zone of forest 25 soils. However, the mineral dissolution equations in PROFILE are not adapted for the unsaturated zone, and employing these equations on a catchment scale results in a significant over-prediction of base cation release rates to surface waters. In this study we use a revised set of PROFILE equations which, among other features, include retardation from silica concentrations.Relationships between the water transit time (WTT) and soil water concentrations were derived for each base cation, by simulating the soil water chemistry along a one-dimensional flowpath, using the mineralogy from a glacial till soil. We show 30 how the revised PROFILE equations are able to reproduce patterns in BC-and Si-concentrations, as well as BC-ratios (Ca 2+ /BC, Mg 2+ /BC and Na + /BC), observed in soil water profiles and catchment runoff. As opposed to the original set of PROFILE equations, the revised set of equations could reproduce how increasing WTT led to decreasing Na + /BC, as well as increasing Ca 2+ /BC and Mg 2+ /BC. Furthermore, the total release of base cations from a hillslope was calculated using a mixing model, where water of different WTT was mixed according to an externally modelled WTT-distribution. The revised set of 35 SOIL Discuss., https://doi.equations gave a 50% lower base cation release (0.23 eq·m-2·yr-1) than the original PROFILE equations, and are in better agreement with mass balance calculations of weathering rates. The results from this study thus demonstrate that the revised mineral dissolution equations for PROFILE are a major step forward in modelling weathering rates on a catchment scale.

Abstract. Soil mineral weathering is one of the major sources of base cations (BC), which play a dual role in forest ecosystems: they function as plant nutrients and buffer against the acidification of catchment runoff. On a long-term basis, soil weathering rates determine the highest sustainable forest productivity that does not cause acidification. It is believed that the hydrologic residence time plays a key role in determining the weathering rates at the landscape scale. The PROFILE weathering model has been used for almost 30 years to calculate weathering rates in the rooting zone of forest soils. However, the mineral dissolution equations in PROFILE are not adapted for the saturated zone, and employing these equations at the catchment scale results in a significant overprediction of base cation release rates to surface waters. In this study, we use a revised set of PROFILE equations which, among other features, include retardation due to silica concentrations. Relationships between the water transit time (WTT) and soil water concentrations were derived for each base cation, by simulating the soil water chemistry along a one-dimensional flow path, using the mineralogy from a glacial till soil. We show how the revised PROFILE equations are able to reproduce patterns in BC and silica concentrations as well as BC ratios (Ca2+/BC, Mg2+/BC and Na+/BC) that are observed in the soil water profiles and catchment runoff. In contrast to the original set of PROFILE equations, the revised set of equations could reproduce the fact that increasing WTT led to a decreasing Na+/BC ratio and increasing Ca2+/BC and Mg2+/BC ratios. Furthermore, the total release of base cations from a hillslope was calculated using a mixing model, where water with different WTTs was mixed according to an externally modeled WTT distribution. The revised set of equations gave a 50 % lower base cation release (0.23 eq m−2 yr−1) than the original PROFILE equations and are in better agreement with mass balance calculations of weathering rates. Thus, the results from this study demonstrate that the revised mineral dissolution equations for PROFILE are a major step forward in modeling weathering rates at the catchment scale.