Abstract

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.

Highlights

  • The dissolution of minerals in soils is one of the key processes behind the diversity of surface water chemistry around the world

  • We demonstrate how the base cation release rates from mineral dissolution in a hillslope comprised of a range of flow paths with different water transit time (WTT) can be calculated by mixing the water from the mineral dissolution model according to a given WTT distribution

  • In the soil profile S22 (Fig. 2), soil water base cation concentrations were approximately constant with depth within the unsaturated zone

Read more

Summary

Introduction

The dissolution of minerals in soils is one of the key processes behind the diversity of surface water chemistry around the world. In forest ecosystems, mineral dissolution is a major source of base cations. These cations act as tree nutrients in soil water and provide buffering capacity to both soil and the runoff sustaining aquatic ecosystems. Both forest harvest and acid deposition are anthropogenic factors that remove base cations and can potentially create imbalances that threaten soil and aquatic ecosystems (Akselsson et al, 2007). Quantifying the base cation release rates from mineral dissolution is fundamental when determining sustainable rates of forest harvest as well as critical loads of acidifying deposition (Klaminder et al, 2011). Instead key management decisions for controlling acid deposition and forest management practices have been based on approaches that do not explicitly predict weathering rates

Methods
Results
Discussion
Conclusion
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call