Models are needed to understand how plant-soil nutrient stores and fluxes have responded to the last two centuries of widespread anthropogenic nutrient pollution and predict future change. These models need to integrate across carbon, nitrogen, and phosphorus (C, N, and P) cycles and simulate changes over suitable timescales using available driving data. It is also vital that they are constrainable against observed data to provide confidence in their outputs. To date, no models address all of these requirements. To meet this need, a new model, N14CP, is introduced, which is initially applied to Northern Hemisphere temperate and boreal ecosystems over the Holocene. N14CP is parameterized and tested using 88 northern Europe plot-scale studies, providing the most robust test of such a model to date. The model simulates long-term P weathering, based on the assumption of a starting pool of weatherable P (P weath0 , g m À2 ), which is gradually transformed into organic and sorbed pools. Nitrogen fixation (and consequently primary production) is made dependent on available P. In the absence of knowledge about the spatial variability of P weath0 , N14CP produces good average soil and plant variables but cannot simulate variations among sites. Allowing P weath0 to vary between sites improves soil C, N, and P results greatly, suggesting that contemporary soil C, N, and P are sensitive to long-term P weathering. Most sites were found to be N limited. Anthropogenic N deposition since 1800 was calculated to have increased plant biomass substantially, in agreement with observations and consequently increased soil carbon pools.