In this study, a phosphorus (P) module containing the biogeochemical P cycle has been developed and integrated into the forest ecosystem model ForSAFE. The model was able to adequately reproduce the measured soil water chemistry, tree biomass (wood and foliage), and the biomass nutrient concentrations at a spruce site in southern Sweden. Both model and measurements indicated that the site showed signs of P limitation at the time of the study, but the model predicted that it may return to an N-limited state in the future if N deposition declines strongly. It is implied by the model that at present time, the plant takes up 0.50 g P m(-2) y(-1) of which 80% comes from mineralization and the remainder comes from net inputs, i.e. deposition and weathering. The sorption/desorption equilibrium of P contributed marginally to the supply of bioavailable P, but acted as a buffer, particularly during disturbances.