Modeling of global soil organic carbon (SOC) is accompanied by large uncertainties. The heavy computational requirement limits our flexibility in disentangling uncertainty sources especially in high latitudes. We build a structured sensitivity analyzing framework through reorganizing the Organizing Carbon and Hydrology in Dynamic Ecosystems (ORCHIDEE)-aMeliorated Interactions between Carbon and Temperature (MICT) model with vertically discretized SOC into one matrix equation, which brings flexibility in comprehensive sensitivity assessment. Through Sobol's method enabled by the matrix, we systematically rank 34 relevant parameters according to variance explained by each parameter and find a strong control of carbon input and turnover time on long-term SOC storages. From further analyses for each soil layer and regional assessment, we find that the active layer depth plays a critical role in the vertical distribution of SOC and SOC equilibrium stocks in northern high latitudes (>50°N). However, the impact of active layer depth on SOC is highly interactive and nonlinear, varying across soil layers and grid cells. The stronger impact of active layer depth on SOC comes from regions with shallow active layer depth (e.g., the northernmost part of America, Asia, and some Greenland regions). The model is sensitive to the parameter that controls vertical mixing (cryoturbation rate) but only when the vertical carbon input from vegetation is limited since the effect of vertical mixing is relatively small. And the current model structure may still lack mechanisms that effectively bury nonrecalcitrant SOC. We envision a future with more comprehensive model intercomparisons and assessments with an ensemble of land carbon models adopting the matrix-based sensitivity framework.
Keywords: Earth system model; active layer depth; computationally expensive; cryoturbation; soil organic matter.