The mARM3D spatially distributed soil evolution model: Three-dimensional model framework and analysis of hillslope and landform responses

We present a three-dimensional landscape-pedogenesis model, mARM3D (matrices ARMOUR 3D), which simulates soil evolution as a function of erosion and pedogenic processes. The model simulates the discretized soil profile for points on a spatial grid. The approach, using transition matrices, is computationally efficient and allows the simulation of large-scale spatial coupling and long-term soil evolution. We study the effect of the depth-dependent soil-weathering rate (i.e., the rate of soil particle breakdown) and bedrock-lowering rate (i.e., the rate of conversion of bedrock to soil). The difference in depth-dependent weathering functions has a significant effect on the in-profile soil properties through depth, specifically particle size grading. However, the depth dependency has a relatively minor effect on the surface properties of the soil profile, with all weathering functions generating very similar surface properties. The surface properties were a function of the cumulative amount of weathering (i.e., the integral of the weathering function over exhumation) with finer surface grading for higher weathering rates. Soil thickness could be estimated without explicitly modeling soil thickness. Thickness was negatively correlated with surface median grain size. As thickness decreases, the surface grading coarsens. This was driven by surface erosion, where as surface grading coarsens, erosion decreases and the soil deepens. Weathering and erosion interact to spatially organize the surface soil grading with a log-log relationship between surface grading, contributing area, and local slope. This relationship was independent of the weathering function. This relationship might be useful for the spatial description of soil properties in digital soil mapping.