Paddy surface water serves as the primary source to artificial drainage and rainfall runoff leading to phosphorus (P) loss from paddy fields. The quantification of P dynamics in paddy surface water on a large scale is challenging due to the spatiotemporal heterogeneity of influencing factors and the limitations of field measurements. Based on 1,226 data sets from 33 field sites covering the three main rice-growing regions of China (the Southeast Coast, the Yangtze River Basin, and the Northeast Plain), we analyzed the spatiotemporal characteristics of P attenuation in paddy surface water and its influencing factors. A new multi-site and long-term phosphorus estimation model for paddy (MLEpaddy-P) was proposed to evaluate the total phosphorus (TP) dynamics at national scale by improving the initial concentration (C0) and attenuation coefficient (k) of the first-order kinetic model (). Our study showed that: (1) Fertilizer amounts, soil organic matter content, soil Olsen-P content, soil pH, and soil total phosphorus are the primary factors affecting the variation of C0 and k; (2) Yangtze River Basin possessed highest C0 (6.87±12.97 mg/L) and high k≤7 (0.262 in 1-7 days after fertilization), followed by Southeast Coast (4.15±5.33 mg/L; 0.263) and Northeast Plain (1.33±1.50 mg/L; 0.239), respectively; (3) MLEpaddy-P performed well in daily TP dynamics estimation at national scale with R2 of 0.74-0.85; (4) Middle and lower reaches of the Yangtze River Basin were the critical regions with high TP concentration due to high fertilizer amount and soil Olsen-P content. The new universal model realizes the multi-site and long-term inversion of P dynamics while greatly saving multi-site monitoring costs. This study provides a basis for early warning and targeted control of P loss from paddies.