Accurate estimation and spatial allocation of economic optimum nitrogen (N) rates (EONR) can support sustainable crop production systems by reducing chemical compounds to be applied to the ground while preserving the optimum yield and profitability Smart Farming (SF) techniques such as historical precision agriculture (PA) machinery data, satellite multispectral imagery, and on-machine nitrogen adjustment sensors can bring together state-of-the-art precision in determining EONR. The novelty of this study is in introducing an efficient optimization framework using SF technology to enable real-time and prescription based EONR application execution. An optimization strategy called response surface modelling (RSM) was implemented to support decision making by fusing multiple sources of information while keeping the underlying computation simple and interpretable. Here, a field of winter wheat with an area of 7 ha was used to prove the proposed concept of determining EONR for each location in the field using auxiliary variables called multispectral indices (MSIs) derived from Sentinel 2. Three different image acquisition dates before the actual N application were considered to find the best time combination of MSIs along with the best MSIs to model yield. The best MSIs were filtered out through three phases of feature selection using analysis of variance (ANOVA), Lasso regression, and model reduction of RSM. For the date 2020.03.25, 14 out of 21 MSIs exhibited a significant interaction with the N applied as determined through an on-machine N sensor. For dates 2020.03.30 and 2020.04.04, the numbers of significant indices were identified as 6 and 10, respectively. Some of the MSIs were no longer significant after five days of the growth period (5-day interval between Sentinel 2 revisits). The best model demonstrated an average prediction error of 14.5%. Utilizing the model’s coefficients, the EONR was computed to be between 43 kg/ha and 75 kg/ha for the target field. By incorporating MSIs into the fitted model for a given N range, it was demonstrated that the shape of the yield-N relation (RSM) varied due to field heterogeneity. The proposed analytical approach integrates farmer engagement by participatory annual post-mortem analysis. Using the determined RSM approach, retrospective assessment compares economically optimal N input, based on observed MSIs values to each location, with the actual applied rates.
{"title":"In season estimation of economic optimum nitrogen rate with remote sensing multispectral indices and historical telematics field-operation data","authors":"Morteza Abdipourchenarestansofla, Hans-Peter Piepho","doi":"10.1007/s11119-025-10224-6","DOIUrl":"https://doi.org/10.1007/s11119-025-10224-6","url":null,"abstract":"<p>Accurate estimation and spatial allocation of economic optimum nitrogen (N) rates (EONR) can support sustainable crop production systems by reducing chemical compounds to be applied to the ground while preserving the optimum yield and profitability Smart Farming (SF) techniques such as historical precision agriculture (PA) machinery data, satellite multispectral imagery, and on-machine nitrogen adjustment sensors can bring together state-of-the-art precision in determining EONR. The novelty of this study is in introducing an efficient optimization framework using SF technology to enable real-time and prescription based EONR application execution. An optimization strategy called response surface modelling (RSM) was implemented to support decision making by fusing multiple sources of information while keeping the underlying computation simple and interpretable. Here, a field of winter wheat with an area of 7 ha was used to prove the proposed concept of determining EONR for each location in the field using auxiliary variables called multispectral indices (MSIs) derived from Sentinel 2. Three different image acquisition dates before the actual N application were considered to find the best time combination of MSIs along with the best MSIs to model yield. The best MSIs were filtered out through three phases of feature selection using analysis of variance (ANOVA), Lasso regression, and model reduction of RSM. For the date 2020.03.25, 14 out of 21 MSIs exhibited a significant interaction with the N applied as determined through an on-machine N sensor. For dates 2020.03.30 and 2020.04.04, the numbers of significant indices were identified as 6 and 10, respectively. Some of the MSIs were no longer significant after five days of the growth period (5-day interval between Sentinel 2 revisits). The best model demonstrated an average prediction error of 14.5%. Utilizing the model’s coefficients, the EONR was computed to be between 43 kg/ha and 75 kg/ha for the target field. By incorporating MSIs into the fitted model for a given N range, it was demonstrated that the shape of the yield-N relation (RSM) varied due to field heterogeneity. The proposed analytical approach integrates farmer engagement by participatory annual post-mortem analysis. Using the determined RSM approach, retrospective assessment compares economically optimal N input, based on observed MSIs values to each location, with the actual applied rates.</p>","PeriodicalId":20423,"journal":{"name":"Precision Agriculture","volume":"8 1","pages":""},"PeriodicalIF":6.2,"publicationDate":"2025-02-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143435032","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"农林科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-02-17DOI: 10.1007/s11119-025-10225-5
Manushi B. Trivedi, Terence R. Bates, James M. Meyers, Nataliya Shcherbatyuk, Pierre Davadant, Robert Chancia, Rowena B. Lohman, Justine Vanden Heuvel
The ability to reduce sampling distance or time is crucial for growers to monitor vineyard nutrients more frequently. Extension specialists often recommend collecting large random samples, but this is frequently overlooked, leading to inaccurate fertilizer recommendations. A novel, one-location square grid area-based sampling method called “box” sampling was developed to capture the overall nutrient distribution within a block, providing guidance for growers on sample collection in vineyards for nutrient monitoring. Box sampling was compared with random and stratified sampling methods at both bloom and veraison for grapevine foliar nitrogen (N%), phosphorus (P%), potassium (K%), magnesium (Mg%), and calcium (Ca%). Box and stratified sampling locations were determined based on Synthetic Aperture Radar (SAR) from Sentinel-1 and Sentinel-2 Normalized Difference Vegetation Index (NDVI) images. SAR and NDVI images were stratified into three variability zones using the k-means + + algorithm. Representative pixels from each zone were sampled using the stratified method, while the junction of these variability zones (30mx30m sampling window) was sampled using the new box method. In 2021 and 2022, these methods were compared against nutrient population parameters in two vineyard blocks. Both methods showed marginal differences in mean, median, and standard deviation, with box sampling consistently capturing a broader range of variations. This was evidenced by the Bhattacharya coefficient, which indicates the overlap between two probability distributions (with values closer to 1 for greater overlap). The coefficient was > 0.80 for N%, P%, and Mg%, and > 0.60 for K% and Ca% at both bloom and veraison. For 14 different commercial vineyards in 2022 and 2023, box sampling accurately captured random nutrient variability for N%, P% and Mg% at both bloom and veraison. However, for K% (at veraison) and Ca% box sampling performed poorly due to high spatial variability. Box sampling reduced the sampling distance and time by 75% compared to random sampling.
{"title":"Box sampling: a new spatial sampling method for grapevine macronutrients using Sentinel-1 and Sentinel-2 satellite images","authors":"Manushi B. Trivedi, Terence R. Bates, James M. Meyers, Nataliya Shcherbatyuk, Pierre Davadant, Robert Chancia, Rowena B. Lohman, Justine Vanden Heuvel","doi":"10.1007/s11119-025-10225-5","DOIUrl":"https://doi.org/10.1007/s11119-025-10225-5","url":null,"abstract":"<p>The ability to reduce sampling distance or time is crucial for growers to monitor vineyard nutrients more frequently. Extension specialists often recommend collecting large random samples, but this is frequently overlooked, leading to inaccurate fertilizer recommendations. A novel, one-location square grid area-based sampling method called “box” sampling was developed to capture the overall nutrient distribution within a block, providing guidance for growers on sample collection in vineyards for nutrient monitoring. Box sampling was compared with random and stratified sampling methods at both bloom and veraison for grapevine foliar nitrogen (N%), phosphorus (P%), potassium (K%), magnesium (Mg%), and calcium (Ca%). Box and stratified sampling locations were determined based on Synthetic Aperture Radar (SAR) from Sentinel-1 and Sentinel-2 Normalized Difference Vegetation Index (NDVI) images. SAR and NDVI images were stratified into three variability zones using the <i>k</i>-means + + algorithm. Representative pixels from each zone were sampled using the stratified method, while the junction of these variability zones (30mx30m sampling window) was sampled using the new box method. In 2021 and 2022, these methods were compared against nutrient population parameters in two vineyard blocks. Both methods showed marginal differences in mean, median, and standard deviation, with box sampling consistently capturing a broader range of variations. This was evidenced by the Bhattacharya coefficient, which indicates the overlap between two probability distributions (with values closer to 1 for greater overlap). The coefficient was > 0.80 for N%, P%, and Mg%, and > 0.60 for K% and Ca% at both bloom and veraison. For 14 different commercial vineyards in 2022 and 2023, box sampling accurately captured random nutrient variability for N%, P% and Mg% at both bloom and veraison. However, for K% (at veraison) and Ca% box sampling performed poorly due to high spatial variability. Box sampling reduced the sampling distance and time by 75% compared to random sampling.</p>","PeriodicalId":20423,"journal":{"name":"Precision Agriculture","volume":"49 1","pages":""},"PeriodicalIF":6.2,"publicationDate":"2025-02-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143435072","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"农林科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}