Microbially-explicit soil biogeochemistry (BGC) models are thought to be more biologically realistic than conventional models, and improved model structures can potentially better represent field- and lab-scale observations. The microbe-enabled BGC module will be implemented in the PFLOTRAN-BGC framework, which has been developed under the NGEE-Arctic project. PFLOTRAN can solve a system of nonlinear partial differential equations describing multi-phase, multi-component and multi-scale 3-D flow and reactive-transport in porous media. We have developed a generic interface in the Energy Exascale Earth System Model (E3SM) through the CMDV project. The interface facilitates the coupling of PLFOTRAN into the E3SM Land Model (ELM). The development of this interface is to enable flexible and fast development and evaluation of soil BGC modules and their coupling to various thermal-hydrology (TH) and aboveground vegetation modules. The interface includes a generic data-structure to pass data between submodels (i.e., vegetation, TH and BGC) and allows users to select a specific submodel from multiple options (e.g., ELM-BGC or PFLOTRAN-BGC for BGC). We evaluate ELM-PFLOTRAN and compare it to the original ELM and observations at a field site near Barrow, AK. The ELM-PFLOTRAN simulations (e.g., GPP, LAI, and total soil organic carbon) were improved by modifying the nitrogen (N) uptake profile with N fixation profile (representing root distribution). ELM-PFLOTRAN using the nutrient-limited kinetics can reproduce the relative-demand behavior in ELM. A complete evaluation of ELM-PFLOTRAN at the Barrow site is underway and the global run of ELM-PFLOTRAN BGC is being tested.