Recent observations indicate that long-term N additions can slow decomposition, leading to C accumulation in soils, but this process has received limited consideration by models. To address this, we developed a model of soil organic matter (SOM) dynamics to be used with the PnET model and applied it to simulate N addition effects on soil organic carbon (SOC) stocks. We developed the model’s SOC turnover times and responses to experimental N additions using measurements from the Harvard Forest, Massachusetts. We compared model outcomes to SOC stocks measured during the 20th year of the Harvard Forest Chronic Nitrogen Amendment Study, which includes control, low (5 g N m⁻² yr⁻¹) and high (15 g N m⁻² yr⁻¹) N addition to hardwood and red pine stands. For unfertilized stands, simulated SOC stocks were within 10 % of measurements. Simulations that used measured changes in decomposition rates in response to N accurately captured SOC stocks in the hardwood low N and pine high N treatment, but greatly underestimated SOC stocks in the hardwood high N and the pine low N treatments. Simulated total SOC response to experimental N addition resulted in accumulation of 5.3–7.9 kg C per kg N following N addition at 5 g N m⁻² yr⁻¹and 4.1–5.3 kg C per kg N following N addition at 15 g N m⁻² yr⁻¹. Model simulations suggested that ambient atmospheric N deposition at the Harvard Forest (currently 0.8 g N m⁻² yr⁻¹) has led to an increase in cumulative O, A, and B horizons C stocks of 211 g C m⁻²(3.9 kg C per kg N) and 114 g C m⁻²(2.1 kg C per kg N) for hardwood and pine stands, respectively. Simulated SOC accumulation is primarily driven by the modeled decrease in SOM decomposition in the Oa horizon.
