Abstract. Reliable soil biogeochemical modeling is a prerequisite for credible projections of climate change and associated ecosystem feedbacks. This recognition has called for frameworks that can support flexible and efficient development and application of new or alternative soil biogeochemical modules in earth system models (ESMs). The BeTR-v1 code (i.e., CLM4-BeTR) is one such framework designed to accelerate the development and integration of new soil biogeochemistry formulations into ESMs, and to analyze structural uncertainty in ESM simulations. With a generic reactive transport capability, BeTR-v1 can represent multi-phase (e.g., gaseous, aqueous, and solid), multi-tracer (e.g., nitrate and organic carbon), and multi-organism (e.g., plants, bacteria and fungi) dynamics. Here, we describe the new version BeTR-v2, which adopts more robust numerical algorithms and improves on structural design over BeTR-v1. BeTR-v2 better supports different mathematical formulations in a hierarchical manner by allowing the resultant model be run either for a single topsoil layer, a vertically resolved soil column, or fully coupled with the land component of the Energy Exascale Earth System Model (E3SM). We demonstrate the BeTR-v2 capability with benchmark cases and example soil BGC implementations. By taking advantage of BeTR-v2’s generic structure integrated in E3SM, we then found that calibration could not resolve biases introduced by different numerical coupling strategies of plant-soil biogeochemistry. These results highlight the importance of numerically robust implementation of soil biogeochemistry and coupling with hydrology, thermal dynamics, and plants— capabilities that the open-source BeTR-v2 provides.