Intra-field heterogeneity of soil properties, such as soil organic carbon (SOC), nitrogen (N), phosphorous (P), exchangeable cations, pH, or soil texture, is a function of complex interactions between biological factors, physical factors, and historic agricultural management. Mapping the crop growth and final yield heterogeneity and quantifying their link with soil properties can contribute to an optimization of amendment/fertilizer application and crop yield in a management variable zones (MVZ) approach. To this end, we studied a field of 17 ha consisting of four former fields that were merged in early 2017 and cropped with winter wheat in 2018. Historical management practices data were collected. The topsoil characteristics were analyzed by grid-based sampling and kriged to create maps. We tested the capacity of a multispectral MicaSense^® RedEdge-M^TM camera sensor embedded on an unmanned aerial vehicle (UAV) to map in-season growth of winter wheat. Relating several vegetation indices (VIs) to the plant area index (PAI) measured in the field highlighted the red-edge NDVI (RENDVI) as the most suitable to follow the crop growth throughout the growing season. The georeferenced final grain yield of the winter wheat was measured by a combine harvester. The spatial patterns in RENDVI at three phenological stages were mapped and analyzed together with the yield map. For each of these images a conditional inference forest (CI-forest) algorithm was used to identify the soil properties significantly influencing these spatial patterns. Historical management practices of the four former fields have induced significant heterogeneity in soil properties and crop growth. The spatial patterns of RENDVI are rather constant over time and their Spearman rank correlation with yield is similar along the growing season (r ≃ 0.7). Soil properties explain between 87% (mid-March) to 78% (mid-May) of the variance in RENDVI throughout the growing season, as well as 66% of the variance in yield. The pH and exchangeable K are the most significant factors explaining from 15 to 26% of the variance in crop growth. The methodology proposed in this paper to quantify the importance of soil parameters based on the CI-forest algorithm can contribute to a better management of amendment/fertilizer inputs by stressing the most important parameters to take into consideration for site-specific management. We also showed that heterogeneity induced by the soil properties can be described by a crop map early in the season and that this crop map can be used to optimize soil sampling and thus amendment/fertilizer management.

Intra-field heterogeneity of soil properties is function of complex interaction between biological, physical factors and historic agricultural management. Quantifying the spatial patterns of soil properties such as soil organic carbon (SOC), nitrogen (N), phosphorous, exchangeable cations, pH, soil texture will contribute to an optimization of fertilizer application and crop yields. We tested the capacity of the multispectral Micasense rededge camera mounted on a UAV in order to map the development of winter wheat and related the Red-Edge NDVI (RENDVI) from the sensor to the Plant Area Index (PAI) measured in the field. The geo-referenced grain yield of the winter wheat was measured by a combine harvester and the top soil characteristics analysed by a grid based sampling. The spatial patterns in RENDVI at three phenological stages were mapped together with the yields. For each of these images conditional inference trees were used to derive the soil properties that significantly influenced these spatial patterns. Within-field variation in PAI (cv 41 % in March, 27% in April and 9 % in May) and yield (cv 4%) can be observed. The spatial patterns of RENDVI are rather constant and their correlation with yields is highest in March and April (r=0.64). Soil properties explain between 67 to 79 % of the variance in vegetation index throughout the growing season as well as 67 % of the variance in yield.  Legacy effects of land consolidation two years earlier reflected in the field lay-out, pH and exchangeable K are significant factors explaining around 12-18 % of the variance in crop yield each. The SOC contents were overall low (8-15 g kg-1). Hence, the N supply resulting from SOC mineralization throughout the growing season covered less than 10% of the crop needs and the yield patterns did not reflect the variation in SOC contents.