Fast, precise and affordable soil analytical techniques are needed for the determination of soil fertility of each zone of a field in site-specific land management. The objective of this study was to develop a reliable calibration model with the visible (VIS) and near infrared (NIR) spectroscopy method to predict the phosphorus (P) content in fresh soil as a preliminary step towards the development of an automated variable rate fertiliser distributor. The Zeiss Coronas 45 visnir fibre-type spectrophotometer was used to measure soil spectra in reflectance mode in the wavelength range of 305–1710 nm. A total of 185 soil samples were collected from arable fields and grassland distributed over Flanders (Belgium) amended with 29 and 14 samples taken from two fields near Leuven, (Belgium) to obtain a balanced P data set. Only samples with texture classes of silt loam, sandy loam and loamy sand texture were held, to avoid possible averse effects of small amounts of soil samples with different textures. Fresh samples were subject to chemical analyses and optical reflectance measurements. Soil spectra were correlated with available P using the partial least-squares (PLS) regression and full cross-validation. After removing the noise, the VIS-NIR wavelengths in the range 401–1663 nm were used. From the same data set, two suitable calibration models were developed using different data preprocessing methods. In the validation stage both models predicted available P in the soil with a coefficient of determination R2 of 075 and 073. On two extra date sets (126 and 36 samples, collected from two fields near Leuven), both models were evaluated. The prediction of P on these samples, provided acceptable results (R2 ¼ 068 and 063). The two P models, developed in this research, can be used for rapid prediction of P in fresh soil samples in the laboratory and in situ in fields in Flanders having texture classes: loamy sand, sandy loam and silt loam. This detailed information could b