We developed TECO-CNP Sv1.0, a coupled carbon-nitrogen-phosphorus model with data assimilation for subtropical forests. The model effectively captures observed carbon, nitrogen, and phosphorus pools and fluxes, and significantly improves carbon flux estimates through data assimilation. TECO-CNP provides enhanced biogeochemical cycle representations, enabling more reliable predictions of ecosystem carbon cycle responses to global change through efficient parameter optimization.