Economic and Environmental Feasibility of Variable Rate Nitrogen Application in Potato by Fusion of Online Visible and Near Infrared (Vis-NIR) and Remote Sensing Data

1. Introduction

The essential role of nitrogen (N) is in ensuring the sufficient food supply for 8 billion people is undeniable. However, the overuse of N fertilizers lead to significant release of N into the environment, contaminating air and water. Globally, 40% of total nitrogen applied to the agricultural lands is taken up by crops, yielding a nitrogen use efficiency (NUE) of 40% [1]. On the other hand, about half of N is released to the environment through leaching, and runoff of N compounds into waterbodies and greenhouse gas emissions [1,2]. This not only harms the environment but also poses risk to human health. In addition, ongoing escalation in the N fertilizer use in recent years has not been matched by proportional improvements in crop yield for essential crops. This suggests that the continuous increase in nitrogen pollution is linked to overuse of nitrogen fertilizers. Consequently, agronomists are actively seeking for a strategy that achieves a 'win–win' outcome by improving crop production while simultaneously mitigating nitrogen pollution.

A positive plant response to nitrogen fertilizer application depends on N threshold level, beyond which additional nitrogen input does not increase the crop yield. On the contrary, excessive N input can increase the fertilizer input cost, and decrease yield due to crop lodging problem [3] in case of growing cereal crops. Despite the threshold N limit, farmers frequently apply N fertilization more than plant biological need with the aim to increase crop yield potential [4]. This over-application not only increase unnecessary cost and the environmental pollution but also deteriorate the soil quality by increasing soil acidification [5,6,7]. Having said that it is suggested that advanced and sustainable nutrient management practices can help to enhance crop production and mitigate N pollution through minimizing the N loses [8].

In context of financial constraints faced by farmers and growing environmental concerns, a range of strategies are being implemented to address the challenges associated with overuse of N fertilizers. These strategies include monitoring of manure production, mandating adherence to nitrogen fertilizer recommendations and performing post-harvest analysis of mineral N in soil [9,10]. However, convincing farmers to adopt most efficient N fertilization practices proves challenging, because of crucial significance of N for affecting crop yield. Therefore, it is an urgent need to develop an innovative and effective approach of N fertilizer use for sustainable crop production.

In traditional agricultural practices, spatial variability within a field is ignored and fields are considered as a single unit, whose soil characteristics, topography and environmental factors are spatially homogeneous [11]. Consequently, applying a uniform rate of N fertilization on entire field area may results in over and under-application of N, potentially leading to significant environmental pollution and poor crop yield [12]. Tackling the spatial variability within the field for variable rate fertilization is a key component of precision agriculture and emerges as a sustainable nutrient management approach. In variable rate application, N fertilizers are applied based on spatial variability of soil and crop characteristics within a field [12]. Site-specific application of fertilizer has been demonstrated to increase nutrient use efficiency and net returns, while also addressing the issue of excessive use of fertilizers [13]. Successful implementation of variable rate N fertilization requires high-resolution data on soil and crop characteristics, to enable mapping the spatial variability at sufficient resolution. This can achieved through the use of online proximal soil sensing and remote crop sensing technologies, coupled with chemometrics and machine learning [14]. These measurement techniques should be rapid, cost-effective and convenient and require minimum labore effort. The most reported proximal technologies that fulfill these requirements are online soil sensors, which play a crucial role in the effectiveness of variable rate fertilizers application. A good example is the visible and near infrared (vis-NIR) spectroscopy-based soil sensor [15,16]. However, the potential of these online soil sensors for variable N fertilization was rarely explored, and only for cereal crops, e.g., wheat, barley, oil seed rape [3,17]. However, to the best of our knowledge no works on adopting online soil sensing for variable N fertilization in potato was reported in the literature.

Potato (Solanum tuberosum) is an important food and high cash crop, whose yield quantity and quality is significantly affected by N fertilisation. Therefore, improvement of N use efficiency is of high priority to enhance potato tuber yield and quality. It can be hypothesized in this work that online vis-NIR soil sensing-based variable rate N fertilization in potato increases potato yield at reduced N application rate compared to the uniform rate N fertilization. Therefore, this work evaluates the economic, environmental and agronomic benefits of variable rate nitrogen application in potato, based on management zones, derived from the fusion of multiple soil properties measured with an online vis-NIR sensor [16] with remote sensing measured crop normalized difference vegetation index (NDVI).

2. Materials and Methods

2.1. Experimental Field

The experiment was conducted in one field (Blondel-3) of an area of 6.62 ha in a commercial farm in the northwest of Belgium near the French boarder (51°01'47.3"N 2°33'31.6"E) (Figure 1). This region has a mean annual temperature of 6 to10 ℃ and a mean annual precipitation, ranging from 700 to 850 mm. Soil texture type of experimental field is a sandy loam determined by texture-by-feel method by Soil Service of Belgium. The topography is flat and has almost no undulation. The experiment was carried out in 2023 cropping season, before which the field was cultivated with barley.

2.2. Online Soil Scanning and Sampling

Online measurement of soil spectra was carried out in August 2022 by using the online vis-NIR multi-sensor platform [16] (Figure 2). The platform consists of a subsoiler fitted to a metal frame attachable to tractor’s three-point linkage. The subsoiler creates 15-20 cm deep trenches, whose bottom is smoothened by downward force of acting on the subsoiler chisel. An optical probe hosted by a mild steel lens holder is attached to backside of the subsoiler heel. Soil reflectance spectra are measurement while the tractor is driven at 3.5 km/h. Online spectra with a spectral range of 305-1700 nm and resolution of 1 nm is recorded by a fiber-type spectrophotometer (CompactSpec Tec5 Technology, Germany). A 100% ceramic disk with 50 mm diameter was used as a white reference. Recalibration of this reference was performed at regular intervals of 30 minutes. Position of soil spectra was recorded by a Differential Global Positioning System (DGPS) (Trimble AG., Trimble Navigation Ltd., Sunnyvale, Canada). Each soil spectrum with georeferenced data was recorded at a rate of 1 Hz in a laptop computer using MultiSpec Pro-II software from Tec5 Technology, Germany. The sensing platform was driven along parallel transects spaced at 12 m apart. During online soil scanning, 13 soil samples were collected randomly from the bottom of the 12 trenches. On average, the data logger captured more than 800 online spoil spectra per hectare. These samples were merged with selected samples from our spectral library to produce calibration models for several soil properties, as detailed below.

2.3. Laboratory Scanning and Analysis

Each soil sample was mixed well and cleaned manually to remove visible plant roots and stones then ground to pass through a 2 mm sieve. Employing the quartering method proposed by Mukhopadhyay and Maiti [18], each soil sample reduced to 150g which was used for the chemical analysis in laboratory. Soil moisture content (MC) of each soil sample was determined by using oven-drying method at 105 C for 24 hours. Subsequently, soil chemical analyses were performed at soil laboratory of Soil Service Belgium by using standard protocol (Table 1). Soil pH was measured in a 1:2.5 soil: water suspension. Soil total organic carbon (TOC) was determined by dry combustion according to Dumas principle, followed by removal of inorganic carbon with HCl treatment. The extractable soil K, P, Mg and Ca were determined in an ammonium lactate extract and analyzed by using inductively coupled plasma atomic emission spectroscopy [19].

2.4. Online vis-NIR Spectral Prediction Models and Maps

Vis-NIR calibration models for soil MC, pH, TOC, P, K, Mg, Ca and CEC were developed using partial least squares regression (PLSR) after spectra preprocessing using R studio (RStudio Inc, USA) with free available libraries [20]. The 13 soil samples collected from the experimental field were insufficient to develop the PLSR calibration models. Therefore, 97 soil spectra from our spectral library, collected from different fields in the same farm were used from the development of calibration models. Adding from samples from the experiment yielding a total of 110 soil samples that were used in modeling.

Spectral modeling was started with preprocessing of spectra by using several algorithms and the best performing combination of algorithms was used. Before preprocessing, edge trimming was executed to eliminate the noisy boundaries within the range of 305–429 nm and 1666–1700 nm, resulting in a final spectral range of 430–1665 nm that was used for subsequent analysis. Different spectra range was selected for different soil properties to ensure the best possible prediction accuracy. The correction of spectral jump occurring at junction between the visible and NIR detectors, at 1045 nm was done following the approach proposed by Mouazen et al. [21]. The parameters applied for the different preprocessing steps are shown in Table 2 for moving average (MV), normalization, Savtizky and Golay (SG) derivative, gap segment (GS) and SG smoothing. The processed spectra was randomly divided into a 70% calibration set and a 30% prediction set. The PLSR analysis with cross-validation was conducted on the calibration dataset and output models were validated using the prediction dataset. The prediction accuracy was evaluated by means of coefficient of determination (R²), root mean square error (RMSE), ratio of prediction to deviation (RPD) and ratio of performance to inter-quartile distance (RPIQ). Spectral analysis and PLSR modelling was done by using "pls" package, as provided by Mevik and Wehrens [22], and accessible on R CRAN. Models providing the best accuracy using the prediction set were used for the prediction of all vis-NIR spectra collected during online measurement.

After successful prediction of soil properties using online soil spectra, the prediction dataset containing DGPS coordinates was used to develop maps of individual soil properties. High-resolution maps necessitates fitting a semi-variogram model for each soil property, followed by interpolation using ordinary kriging (OK) available in ArcGIS (ESRI ArcGIs v10.7, Redlands, CA, USA).

2.5. NDVI Data Collection

To account for crop growth characteristics, NDVI was calculated for barley grown in the previous cropping season. It was not possible to use the NDVI of potato using data from the current cropping season as N fertilizer has to be applied during potato seeding. We used the Band 4 and Band 8 of multi-band images collected on 2022-03-26 from Sentinel-2 open data hub to calculate NDVI index [24], using the following equation:

$$NDVI=\frac{(B8-B4)}{B8+B4)}$$

(1)

2.6. Management Zones Delineation and Application Map

Crop NDVI data and online predicted soil properties were used to delineate management zones (MZ) for nitrogen application map. Before MZs delineation, the different layers of online soil parameters and crop NDVI having different resolution have to be resampled into the same resolution. This was done by a raster analysis to develop a common grids of 5 by 5 m resolution. This raster analysis was followed by K-mean cluster analysis, which was performed on the normalization data to divide the field into limited number of clusters by using R software (RStudio Inc, USA). Each zone has similar soil and crop characteristics. Since no historical yield data was available for this field, it was necessary to discuss with farmer the ranking of MZ into different soil fertility levels for each management zone based on soil and crop NDVI maps. The output of this discussion allowed ranking the four MZs into high soil fertility (VR-H), medium high fertility (VR-MH), medium low fertility (VR-ML) and low fertility zone (VR-L).

In the conventional UR farming, farmers estimate nitrogen fertilization rate based on a composite soil sample, analysed by laboratory chemical methods to including nitrate, and texture. For this study N fertilization rate for potato was 100 L/ha standard liquid fertilizer containing 39% of N .., which was used in this study, for the UR N fertilization treatment. To compare between Strip map for nitrogen application was developed based on fertility management zones by using fishnet function in an open source ArcGIS software. Treatments were designed as variable rate (VR) and uniform rate (UR) parallel to strips. Each strip was 20 m wide parallel to tramline. Nitrogen fertilizer rate was assigned according to fertility MZs for variable rate application. The UR treatment received the recommended N rate by the farmer, which was applied uniformly over the entire treatment. For the VR-N treatment, we followed the Robin Hood principle [3,25], which recommend more N to poor zones (feeding the poor principle) and vice versa for the fertile zones. Therefore, the VR-H zone (plot) received 50% less N fertilizer than UR, VR-MH received 25% less N, VR-ML received 25% higher N fertilizer and VR-L received 50% higher N fertilizer rate than the UR treatment (Figure 3).

2.7. Crop Management and Yield Measurement

This experiment followed standard crop management practices routinely adopted by the farmer for potato production. Prior to seeding, compost and chicken manure were spread across the entire field at the rate of 10 t/ha and 15 t/ha, respectively. Standard liquid N fertilizer (containing 39% N) was applied at variable rate according to MZs recommendation after crop germination on May 15, 2023. The applied potato seed rate was 2200 kg/ha, while different pesticides were applied during the cropping season. Overall, the field received 4 irrigations by means of a gun-type sprinkler irrigation system, applying about 30 mm water each time.

Crop yield was assessed by harvesting potato tuber manually, following the protocol suggested by Munnaf et el. [26]. Sampling points for harvesting were allocated in each polygon (treatment plot), from each of which 11 plants were harvested manually. Yields were then calculated per plot based on area covered by 11 plants, then yield data were extrapolated to strip and treatment (3 strips). Potatoes were cleaned manually, and any rotten potatoes were removed prior to weighing.

2.8. Cost-Benefit and Environmental Analysis

Nitrogen fertilizer cost was calculated in euros (€) from the fertilizer current price and percentage of N in the fertilizer. Liquid N fertilizer with 39% N was applied. The cost of nitrogen fertilizer was 400 €/t. Area for each treatment was calculated from the field map. The current market price for potato yield, obtained from the farmer is 200 €/t. The gross margin was calculated from fertilizer price and current market price of potato tuber for consumption potato. The cost of field scanning was not included at this stage because soil scanning cost evaluated at 25 € per year, should be divided among several variable rate application of farming input resources (e.g., manure, seeding, irrigation, tillage). Simulated profit was calculated for the whole field assuming one treatment is applied in the entire field are. The marginal value for N fertilization (MV_N) in agriculture refers to the additional benefit with applying an extra unit of nitrogen fertilizer. This concept involves assessing the impact of a small change in nitrogen application on the output of interest, typically crop yield. MV_N was calculated by following equation:

Marginal value N (MV_N) was calculated as:

$${MV}_N=\frac{\mathsf{\Delta }N}{\mathsf{\Delta }Y}$$

(2)

Whereas, ΔN is change in N input and ΔY is change in yield.

Yield response index (YRI) measures the relative yield response to variable rate N application (VR-N) compared to a uniform rate N (UR-N). YRI was calculated by following equation.

$$YRI=\frac{\mathrm{Y}\mathrm{i}\mathrm{e}\mathrm{l}\mathrm{d} \mathrm{w}\mathrm{i}\mathrm{t}\mathrm{h} \mathrm{V}\mathrm{R}-\mathrm{N} \mathrm{i}\mathrm{n}\mathrm{p}\mathrm{u}\mathrm{t}}{\mathrm{Y}\mathrm{i}\mathrm{e}\mathrm{l}\mathrm{d} \mathrm{w}\mathrm{i}\mathrm{t}\mathrm{h} \mathrm{U}\mathrm{R}-\mathrm{N} \mathrm{i}\mathrm{n}\mathrm{p}\mathrm{u}\mathrm{t}}$$

(3)

Environmental analysis for N use was performed by calculating nitrogen use efficiency (NUE) that quantifies how effectively N was utilized by the crop in response to crop yield. It was calculated by dividing crop yield by N fertilizer applied in respective treatment [27]. Net nitrogen impact was determined by the change in N input in VR-N treatment with respect to the UR treatment.

2.9. Statistical Analysis

Statistical analyses were performed on laboratory measured soil properties to identify the outlier and assess the overall quality and consistency of soil data by using R studio (RStudio Inc, USA). Significant differences among the various treatments on crop yield were assessed using analysis of variance (ANOVA), followed by post hoc Tukey's HSD test at a significance level of P = 0.05. Relative influence of predictors (soil properties) on crop yield was calculated by using randomforest package and importance function in R studio. This analysis provides a measure of variable importance based on the mean decrease in accuracy when each predictor variable is permuted. Random forest model extracts variable importance scores using the importance function and creates a variable importance plot using varImpPlot function.

3. Results

3.1. Accuracy of Prediction Models of Online vis-NIR Sensor

The accuracy PLSR models demonstrated fair to very good performance across both the calibration set and online prediction of soil properties (Table 3). Significant variations were observed in model accuracy for different soil properties. Accuracy in calibration was relatively higher as compared to online prediction accuracy for most of the soil properties. For instance, in the calibration set, the highest accuracy was observed for soil Mg with R² value of 0.94, RPD value of 5.36 and RPIQ value of 5.51. In online predicted set, Mg accuracy exhibited R² value of 0.73, RPD value of 2.00 and RPIQ value of 2.58. Soil MC, pH, TOC, P, K and CEC accuracies in calibration sets exhibited R² values of 0.94, 0.84, 0.70, 0.87, 0.70, and 0.91 and RPD values of 4.18, 2.56, 1.86, 2.85, 1.85 and 3.53, respectively. While, in online prediction sets, MC, pH, TOC, P, K, Mg and CEC accuracies exhibited R² values of 0.87, 0.72, 0.69, 0.80, 0.71 and 0.64 and RPD values of 2.86, 1.93, 1.85, 2.31, 1.92 and 1.73, respectively. These findings highlight the models efficacy in predicting all soil properties, with particularly notable performance in the calibration set. Overall, the online prediction results evaluated as RPD [17], indicated excellent performance (RPD > 2.5) for MC, very good (RPD = 2.5–2.0) for P and Mg, and fair (RPD = 1.9–1.4) for K and TOC, pH, and CEC.

3.2. Spatial Variation in Soil Fertility

Summary statistics of measured soil properties are shown in Table 1. Online soil predicted maps and NDVI maps showed notable spatial variability, necessitating variable rate management of farming input resources using precision agricultural technologies (Figure 4). For instance, soil pH, P and CEC showed similar pattern of fertility variation, indicating strong correlation among these properties which is evidenced in NDVI and yield map. While K and Mg showed a similar pattern, soil TOC demonstrated a unique spatial distributed pattern. These spatial variations emphasize the importance of employing advanced proximal and remote sensing technology to enhance soil management practices.

Four MZ classes were determined by K mean clustering, which were categorized into different fertility levels namely, H, MH, ML and L (Figure 5). The quality of the MZ ma seems convincing, as it copies the spatial variability pattern of few soil properties and NDVI (Figure 5). D the classification was based on soil fertility status and experience of the farmer of this field, it was not a straight forward decision to be made. According to fertility classification, the middle part of field has MH and H fertility levels, the southern part of the field was of L fertility (low P, NDVI in particular), and northern part was classified as ML fertile zone. Nitrogen application rates were estimated based on these fertility management zones.

3.3. Crop Yield Response to Variable Nitrogen Application

Crop yield response was different among the different variable rate nitrogen application rates (Figure 6) within each MZ. In contrast to the UR treatment, the VR-H treatment, which received 50% less nitrogen than the UR treatment, exhibited an 8.1% increase in potato yield. Additionally, the VR-ML and L treatments, which received 25% and 50% higher nitrogen than the UR treatment, demonstrated increases in potato yield by 14.3% and 2%, respectively. The only zone that has a smaller yield than the control UR was the VR-MH. The impact of predictors on crop yield was assessed using a random forest model. The analysis revealed TOC, P, and CEC as the predominant factors influencing crop yield, with relative influences of 18.02%, 17.33%, and 16.81%, respectively. While, Mg, K, and pH exhibited relative influences of 16.53%, 16.05%, and 15.22%, respectively (Figure 7).

3.4. Cost-Benefit and Environmental Assessment of Variable Rate Nitrogen

Results of the cost benefit analysis demonstrated that VR-N application significantly affected yield throughout the different MZs in the field (Table 4). The yield of VR-H, and VR-ML MZs exceeded that of the UR treatment, whereas the yield for the VR-MH and VR-L MZs was smaller. Overall, the VR-N treatment increased the average yield (60.55 t ha^-1) compared to UR treatment (58.66 t ha^-1), which resulted in corresponding increase in overall gross margin. Among different variable rate treatments, the highest yield (67.06 t ha^-1) obtained by the VR-ML MZ treatment followed by VR-H (63.44 t ha^-1) has compensated to the reduced yield in the VR-MH, and VR-L MZs, resulting in an overall higher relative gross margin of the VR-N treatment of 374.83 (€/ha), compared to that produced by the UR treatment. The simulated relative gross margin of the VR-N treatment for the field was 2481.37 €. The highest yield response index (YRI - ratio of yield of VR-N / yield of UR treatments) was in VR-ML MZ (1.14), followed by VR-H MZ (1.08), VR-L (0.98) and VR-MH (0.95), respectively. Overall, YRI for the VR-N treatment was 1.03 indicating the positive yield response to VR treatment application, compared to UR application. Marginal more N of 0.68 t/kg was consumed in the VR-N treatment, which has resulted in overall in potato yield, indicating improved N use efficiency per unit of N fertilizer applied.

Environmental assessment of N use in VR application was performed by determining N use efficiency (NUE) and net N impact (ΔN) (Table 5). Results showed that NUE was observed highest in VR-H (126.87%), followed by VR-MH (73.92%), UR (58.70%), VR-ML (53.64%) and VR-L 38.2%, respectively. On average, the overall NUE for the VR-N treatment was significantly higher (150.31 %), than that of the UR treatment (58.70 %), confirming the more efficient use of N by the VR-N treatment, compared to the UR treatment. Net N impact (ΔN), indicating the difference in N input in the VR-N treatment with respect to the UR N input, demonstrated a negative value of -1.28 (Table 5). This results explains that the VR-N treatment consumed a larger N input, compared to the UR treatment. The values of ΔN were negative in poor fertility MZs (VR-ML and VR-L), indicating a higher N input compared to UR treatment, whereas positive values can be observed in the VR-H and VR-MH zones.

4. Discussion

4.1. Models Accuracy

Despite encountering low moisture levels (4.05-6.09%) that resulted in a higher noise generated by vibration during online soil scanning, the vis-NIR models demonstrated a satisfactorily performance ranging from fair to excellent prediction of various soil properties (Table 3). Models performance in the cross-validation was slightly higher than that in prediction for most of the soil properties, findings which are comparable with previous studies results [28,29]. The model's performance closely resembled that documented in previous studies [29,30]. In this study, the accuracy obtained for the majority of vis-NIR models met the criteria established by Maleki et al. [31], requiring an R2 > 0.70 to be considered suitable for implementation of site-specific applications. Furthermore, current models showed higher prediction accuracies than those models utilized in other variable rate applications employing the same online vis-NIR sensing system. For instance, Zhang et al. [32] and Guerrero et al. [3], successfully implemented VR manure and nitrogen application, respectively using the same online vis-NIR models, reporting R² values ranging from 0.51 to 0.81. Similarly, Mouazen and Kuang [33] reported a successful site-specific P management using the same online vis-NIR sensor, using in advance prediction model for P with R² = 0.60, compared to the P model of this work having an R² = 0.8 (Table 3). Thus, it can be inferred that current vis-NIR models exhibit sufficient accuracy to determine management zones for VR-N application.

4.2. Within Field Spatial Variability

Considerable spatial variability was identified in the online-predicted soil properties maps (Figure 6). Zones in the field characterized by low pH exhibited similarities with reduced P content, affirming the impact of pH on P. This correlation may be attributed to lower Ca levels in this specific field area. Additionally, the areas with low pH demonstrated elevated Mg and K content, indicating a negative correlation between pH and the levels of Mg and K. The influence of soil pH is crucial for P availability, as P accessibility can be limited in both acidic and alkaline soils [34]. In the current field, pH is in the alkaline range, and this should be corrected for a better P uptake by the crop. The distinct spatial pattern observed in NDVI closely resembled the distribution of P, pH and CEC. Consequently, special consideration is directed in this work and perhaps in a future work towards P, CEC, and NDVI maps in making informed crop management decisions. Furthermore, the NDVI maps confirmed the significance of phosphorus, CEC, and pH as pivotal factors limiting yield, as they exhibited similar spatial patterns.

In previous studies, a lack of similarity between crop yield patterns and NDVI was observed [17,35], highlighting the potential unreliability of relying solely on NDVI for ranking fertility classes in soil management decisions. Since other studies reporting a strong correlation between crop yield and NDVI [36], we opted to use NDVI as an alternative to yield (as historical yield data was not available), along with other soil properties in classifying MZ into different fertility classes (Figure 4). The determination of N application rates relied on the fertility level of each zone and the farmer's experience with current field. Similar practices were employed in previous studies to optimize manure, and P management, resulted in successful variable rate applications [32,33].

4.3. Variable Rate Nitrogen Impact on Crop Yield

Applying right amount of N fertilizer at right time on right place is critical for effective management of N to meet the N needs of plants. Variable rate N application significantly affected crop yield (Figure 4). Nitrogen is the most limiting nutrient in crop production and in most of the plants is higher in concentration than all other mineral nutrients [37]. Nitrogen is applied in greater amount compared to other nutrients which indicates its critical role in potato production [38,39]. According to Moller et al. [40], N availability poses a significant limitation to yield in potato cultivation. To optimize the potential of crop in a given location, it is crucial to apply the appropriate and optimal amount of N, ensuring vitality and maximizing crop performance. In the VR-ML fertility zone, using a 25% higher N rate led to a higher crop yield compared to the standard rate (UR). Similarly, in the VR-H zone, a 50% lower N application than the standard rate still resulted in a relatively higher crop yield. These results align with Munnaf et al. [17], who reported that a high-fertility zone exhibited increased crop yield despite lower nitrogen input. This could be attributed to the possibility that there was a high natural availability of N through mineralization. Even with a lower external N application, the soil's inherent fertility was sufficient to support a relatively higher crop yield. But, to confirm if this was the case for this output it would have been necessary to measure N mineralization rate and available N throughout the cropping sensor. However, this was not performed in the current work, as no sensor technology is available to collect the necessary data on nitrate and mineralization rate. Conversely, in the VR-ML zone, N was a limiting factor, and increasing nitrogen application by 25% improved crop yield compared to the standard rate. However, in the VR-L zone, even applying a 50% higher N input than the UR treatment did not improvement crop yield (Table 5). Indeed, the extra N rate used in the VR-L zone was not utilized by the crop. This can be resulted in a risk of increasing residual N in the soil without a corresponding increase in crop yield [41].`Therefore, a future work should indeed consider monitoring nitrate and measurement of mineralization rates for the calculation of the recommendation of VR-N application, rather than relying on % increase or decrease in N by the Robin Hood method adopted in this work. Increasing N levels in the plant tissue can positively affect crop yield by improving the absorption of solar radiation through enhanced tillering, leaf enlargement, and enhanced photosynthesis [42]. However, it's essential to note that excessive N application can have negative consequences, causing problems like lodging in cereal crop [43] and foliar infections [42], which can ultimately reduce overall yield. For potato the latter issue might be the cause of reducing tuber yields when applying N fertilizer beyond the crop response. Moreover, soil properties such as soil TOC, pH, CEC significantly influence the N availability and crop yield. It is well known that TOC serves as a source of N through microbial decomposition. Breakdown of TOC releases N in various form, including ammonium (NH₄⁺) and organic nitrogen compounds [44]. CEC crucially affects N availability by enhancing the retention of ammonium in high CEC soils [45]. In the low-fertility zone (VR-L) shown in Figure 5, having the lowest pH, CEC, and P, the limited crop yield response to increased N applied (+50 % of that of UR) can be attributed to nutrient imbalances and interactions [46]. Phosphorus deficiency constrains the plant's ability to utilize nitrogen effectively, due to its impact on plant root development and nitrogen metabolism regulation. Moreover, the random forest model confirms that TOC, P and CEC were major influencing factor on crop yield (Figure 7). Thus, the study underscores the importance of VR-N application in precision agriculture, emphasizing the need to balance inputs across varied soil fertility zones and adopt a holistic soil management approach for sustainable crop production.

4.4. Cost-Benefit and Environmental Analysis

Implementation of variable rate N application based on management zones delineation through the fusion of high resolution soil and remote sensing data can increase the profitability by increasing crop yield and reducing excessive nitrogen use [12,17]. Variable rate N input significantly enhanced crop yield compared to the uniform rate (UR) treatment, leading to increased gross margins (Table 4). These findings align with the results reported by Basso et al. [47]. In the VR-H zone, despite receiving 50% less N than the UR treatment, there was a higher crop yield response. Conversely, in the VR-ML zone, receiving 25% higher N input than the UR treatment, the highest yield (67.76 t/ha) was achieved. However, in the VR-L zone, where the crop did not respond to a 50% increase in N compared to the UR rate, the excessive N use in this class exceeded the maximum crop response threshold to N fertilizer [48]. Therefore, it is not surprising to observe almost equal or smaller potato yield, compared to the UR treatment, which is most likely attributed to increase crop sensitivity to foliar infection due biotic stresses. This suggests that 50% increase in N than the recommended N fertilizer should be avoided. Consequently, a recommendation with a 25% change in N input rate proved effective and would not exceed N over the crop N response limit. Overall, VR N application increased the gross margin by 373.84 € per hectare, with a marginal N increase of 0.68 (t/kg), signifying an increase in crop yield with each additional kilogram of N input. This study demonstrates that VR-N fertilization can offer additional profits to farmers compared to UR rate application.

Compared to the UR treatment, N use per hectare increased by 3.25% with the VR-N treatment, however, N remains below the environmental threshold level while contributing to increased crop yield and overall profitability (Table 4). Environmental analysis of N use revealed a significant increase in NUE, compared to the UR treatment, attributed to a better crop yield response to site-specific N application. Precision N application based on MZs delineation is shown to be both environmentally and economically profitable [3,49]. However, recommending a 25% change in N input rate, instead of a 50% change, could potentially yield greater profits and environmental benefits compared to the UR treatment.

5. Conclusions

Implementation of VR-N fertilization in potato based on MZ delineation using high resolution data on soil fertility parameters and crop NDVI measured by means of an online vis-NIR spectroscopy sensor, and Sentinel 2 satellite imagery, respectively has increased crop yield and profitability indicated as the gross margin. The highest fertile zone receiving 50% lower N input and the medium-low fertile zone receiving 25% more N rate than the UR treatment noticeably increased the crop yield, hence, increased the gross margin of the VR-N treatment over the UR treatment. Marginal N value of 0.68 (t/kg) applied in the VR-N treatment has resulted in increased yield and gross margin. This suggests that the crop has used this minor extra N to produce more yield, minimizing the chance for excess N leached to the environment. Moreover, a simulated gross margin of 2481 € was observed from the VR-N fertilisation for the whole field. These results need to be validated further for potato by running similar work in several fields to confirm the potential profitability, with minor increase in applied N that has no negative effect on the environment, as long as the extra N applied is consumed by the crop to produce more yield compared to the traditional UR fertilisation. Furthermore, historical yield data should be incorporated into decision-making process to improve VR-N recommendation. Similarly data on nitrate and nitrogen mineralisation rate per MZ should also be monitored and used in the decision making to overcome the Robin Hood ± % change in N rate adopted in this work.