Phosphorus limitation on CO2 fertilization effect in tropical forests informed by a coupled biogeochemical model

Tropical forests store more than half of the world’s terrestrial carbon (C) pool and account for one-third of global net primary productivity (NPP). Many terrestrial biosphere models (TBMs) estimate increased productivity in tropical forests throughout the 21st century due to CO2 fertilization. However, phosphorus (P) limitations on vegetation photosynthesis and productivity could significantly reduce the CO2 fertilization effect. Here, we used a carbon-nitrogen-phosphorus coupled model (Dynamic Land Ecosystem Model; DLEM-CNP) with heterogeneous maximum carboxylation rates to examine how P limitation has affected C fluxes in tropical forests during 1860–2018. Our model results showed that the inclusion of the P processes enhanced model performance in simulating ecosystem productivity. We further compared the simulations from DLEM-CNP, DLEM-CN, and DLEM-C and the results showed that the inclusion of P processes reduced the CO2 fertilization effect on gross primary production (GPP) by 25% and 45%, and net ecosystem production (NEP) by 28% and 41%, respectively, relative to CN-only and C-only models. From the 1860s to the 2010s, the DLEM-CNP estimated that in tropical forests GPP increased by 17%, plant respiration (Ra) increased by 18%, ecosystem respiration (Rh) increased by 13%, NEP increased by 121% per unit area, respectively. Additionally, factorial experiments with DLEM-CNP showed that the enhanced NPP benefiting from the CO2 fertilization effect had been offset by 135% due to deforestation from the 1860s to the 2010s. Our study highlights the importance of P limitation on the C cycle and the weakened CO2 fertilization effect resulting from P limitation in tropical forests.