Next Article in Journal
Can Miscanthus Fulfill Its Expectations as an Energy Biomass Source in the Current Conditions of the Czech Republic?—Potentials and Barriers
Previous Article in Journal
Inline Milk Lactose Concentration as Biomarker of the Health Status and Reproductive Success in Dairy Cows
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Remote and Proximal Sensing Applications for Durum Wheat Nutritional Status Detection in Mediterranean Area

Department of Agricultural, Food and Environmental Sciences (D3A), Agronomy and Crop Science Section, Marche Polytechnic University, 60131 Ancona, Italy
*
Author to whom correspondence should be addressed.
Agriculture 2021, 11(1), 39; https://doi.org/10.3390/agriculture11010039
Submission received: 26 November 2020 / Revised: 27 December 2020 / Accepted: 4 January 2021 / Published: 8 January 2021
(This article belongs to the Section Crop Production)

Abstract

:
Combining remote and proximal sensing in agriculture is essential to monitor crop spatial-temporal variability and to provide high-quality prescription maps for the precision agriculture applications. The study showed how different combinations of soil management (no tillage—NT vs. conventional tillage—CT) and nitrogen (N) fertilization levels (0.90 and 180 kg N ha−1) can affect the durum wheat nutritional status and development through vegetation indices computation and proximal sensing tool application. Chlorophyll and N crop content were measured, in addition a proximal sensing tool and multispectral imagery equipped on unmanned aerial vehicle were used. The N input is the key driver for durum wheat development (4.5 ± 0.92 t ha−1 on average), but when it was not provided the NT performed better than CT (2.51 ± 0.22 vs. 1.46 ± 0.28 t ha−1 respectively) in terms of grain yield. This is due to the greater content of organic matter and N availability which characterizes the NT system. The near infrared (NIR) band-based vegetation indices can well detect the durum wheat nutritional status (R2 = 0.70 on average). The showed results can provide an important contribution in the implementation of ago-environmental policies aimed at environmental impact of cereal-based-cropping systems reduction.

1. Introduction

In the European Union, the production of wheat is estimated at 152 million kg 103, about 10% above the harvest of 2018 [1]; in detail in Mediterranean semi-arid climate regions durum wheat (Triticum turgidum L. subsp. durum) is among the most common winter cereal crops [2,3]. In this geographic area drought stress as a combination of erratic rainfall distribution and high temperature during grain filling stages, is the main constraint limiting crop production [4] and determining their instability.
In the near future, the Mediterranean area is expected to face more severe drought and an increase in average temperature, due to climate change [5]; climate change and the simultaneous global population increase, estimated at 2.6 billion by 2050 [6], require a range of environmental measures such as optimizing resources [7], but production must remain high and environmental degradation must be reduced at the same time. If modern agriculture with conventional tillage techniques have allowed high yields, it is also equally true that such soil management induced greater soil erosion [3].
Both the reduced soil fertility and the effects of climatic change strongly affect durum wheat yields and quality traits [8]; to counteract such criticality, conservative agriculture (CA) practices-based on minimum or no soil disturbance, permanent residue cover and planned crop rotations or associations [9] should be applied; these aim to promote water conservation, biological processes, and organic matter accumulation due to reduced soil erosion and CO2 emissions [10,11,12,13].
However, not all studies on CA practices have led to the same results regarding the crop yield. In a dry context, the CA can have a positive effect on crop productivity due to increased soil moisture and nutrient availability [14]. In the literature, there are often contradictory results when the long-term impact of CA practices for durum wheat yield is analysed.
This is due to differences among experimental sites, such as soil and weather conditions, different experimental designs, and different agronomic management of the site [15,16] for example [17] reported 20 years of experiment results on durum wheat grain yield, and concluded that for durum wheat, NT did not allow substantial yield benefits leading to comparable yields with respect to CT in 10 out of 20 years.
However, these studies have always focused on the final grain yield and not on the crop cycle or nutritional status. This can be easily monitored using a precision agriculture (PA) instrument like remote sensing.
In fact, PA allows the use of farm inputs such as fertilizers and herbicides to be optimized [18,19], and thus maximize profit and minimize negative environmental impacts [16] by addressing spatial variability.
PA methods hold great potential as a tool for: high-throughput phenotyping for plant breeding [20,21], decision making for precision agriculture [22,23], predicting yields [24] and predicting spatial field variability in experimental sites [25].
In order to apply PA methods, spatial-temporal crop variability must be measured, and to do so remote and proximal sensing tools can be used and combined [26]. Their usefulness relies on the fact that they are non-destructive, low time-consuming methods, moreover they are also well-correlated with agronomical and important physiological crop traits [27].
An example of those sensing tools is the unmanned aerial vehicle (UAV) equipped with multispectral camera. This allows to monitor the crop with high spatial-temporal resolution during the production season [28].
Spectral remote sensing is a widely used tool for agricultural production and crop management with respect to site-specific fertilizer applications [29,30]; its potential as a vehicle-based high-throughput phenotyping technology under field conditions has been demonstrated [31,32,33,34].
Studies regarding the effect of CA practices and several nitrogen fertilization levels on chlorophyll content and nitrogen status of durum wheat in relation to vegetation indices are lacking. For example, some studies showed that the yield estimated from the normalized difference vegetation index (NDVI) had a strong relationship with grain yield in wheat [35,36] in addition, NDVI has also been used to estimate crop growth status based on the different patterns of reflection of green organs and soil in wheat and other cereals [37,38,39,40].
Following the approach of [41], it is necessary implement the knowledge of the relationships between field data (such as chlorophyll and nitrogen content in the crop, yield components as the number of kernels per spike) and remote-sensing data (such as several vegetation indices) in relation to the type of management operated on durum wheat.
Based on the above considerations, the aims of this study are to evaluate how different soil management combinations and nitrogen fertilization levels can affect:
  • the durum wheat nutritional status, chlorophyll content and yield components;
  • the effectiveness of proximal and remote sensing tools to well correlate the durum wheat chlorophyll content and nutritional status; in addition, the scientific findings emerging from this experiment will be used to provide a contribution in the implementation of the agri-environmental policies adopted in the examined context.

2. Materials and Methods

2.1. Experimental Site

The experimental site is located in the “Pasquale Rosati” experimental farm of the Polytechnic University of Marche in Agugliano, Italy (43°32′ N, 13°22′ E, at an altitude of 100 m above sea level and a slope gradient of 10%) (Figure 1) on a Silty-clay soil classified as Calcaric Gleyic Cambisol [42].
The weather data were provided by ASSAM (Agenzia Servizi al Settore Agroalimentare delle Marche, Osimo stazione, Ancona, Italy). The thermo-pluviometry trend data (Table 1) showed a similar trend in both temperatures and rainfall with a slight difference of 0.4 °C in the mean air temperature and 30 mm of rainfall between the years 2018 and 2019.
Soil sampling was performed with a hand auger (model: Suelo HA-3, Zhejiang Lujian Instrument Equipment Co., Ltd., Zhejiang, China) immediately before sowing. From each subplot, two samples were taken at a depth of 0–0.20 m and geo-referenced using the Leica Zeno 20 (Leica Geosystem, Heerbrugg, Canton St. Gallen, Switzerland). For each soil sample have been measured the content of sand (g kg−1), silt (g kg−1) and clay (g kg−1) by using the hydrometer method, pH with the pH meter method, Organic matter (g kg−1) with the Walkley–Black chromic acid wet oxidation method, the total nitrogen (g kg−1) with the Kjeldahl method and the C/N as the ratio of previous measures (Table 2).

2.2. Experimental Design and Crop Management

The experimental site is a part of a long-term experiment (LTE) established in 1994 but still on-going [41] and consists of a rainfed 2-year rotation with durum wheat (Triticum turgidum L. subsp. durum cv. Grazia, ISEA) and maize (Zea mays L., DK440 hybrid Dekalb Monsanto, FAO Class 300) [17].
In this study, research activities focused on durum wheat plots where two soil management practices (main plot, 1500 m2) and three nitrogen (N) levels (sub plot 500 m2) were repeated in the same plots every year and arranged according to a split plot experimental design with two replications.
Plots with conventional tillage (CT), which is representative of the business as usual tillage practice in the study area, were ploughed along the maximum slope every year by a mouldboard plough (with 2 ploughs) at a depth of 0.40 m in autumn.
The seedbed was prepared with harrowing before the sowing date. The no tillage (NT) soil was left undisturbed and sprayed with herbicides before sowing prior to direct seed drilling.
The three N fertilizer treatments (N0, N1 and N2) corresponded to 0.90 and 180 kg N ha−1 were distributed in two rates. The unfertilized N0 treatment was chosen as control. The N1 treatment was compliant with the agri-environmental measures adopted within the Rural Development Plans at local scale. The N2 treatment was the business as usual N rates in the study area at the start of the experiment. Dates (dd/mm /yyyy) of all agronomic practices are reported in Table 3.

2.3. Measurements

In both growing seasons under analysis (2018 and 2019) have been acquired Soil Plant Analysis Development (SPAD) readings (Konica Minolta Sensing 2003, Osaka, Japan), leaves chlorophyll concentration (mg g−1), N content (%), nitrogen nutrition index (NNI) at the crop tillering, stem elongation and anthesis phenological stages (ZS22, ZS35 and ZS60 respectively were ZS = Zadoks scale [43]).
In addition, multispectral images have been acquired (MAIA S-2, SAL engineering, Russi, Italy) by using an unmanned aerial vehicle (UAV) platform (DJI Matrice 600 pro, SZ DJI Technology Co., Ltd., Shenzen, China) to compute several vegetation index (VI) algorithms.
At crop maturity (ZS92), the yield components have been measured by the number of kernels per spike (KS), the spikes m−2 and the grain yield (t ha−1) for both years.

2.3.1. Soil Plant Analysis Development (SPAD) Readings

For each plot and for each field survey day, three test areas were randomly selected and georeferenced the biomass sample positions by using the Leica Zeno 20 (Figure 1) for a total of 36 ground control points (GCPs). At each test area, 10 fully expanded and intact leaves were chosen to make the readings with the SPAD (mod. Minolta 502). The readings were taken on the central and distal portion of the leaves, at the same time slot, around midday (11:00–13:00 a.m.).
The 10 leaves for each test area, after the SPAD readings, were cut and sealed in a plastic bag, placed in a portable refrigerator and transferred to the laboratory where the leaves’ chlorophyll concentration (mg g−1) was measured.

2.3.2. Leaves Chlorophyll Concentration

The leaf chlorophyll concentration (mg g−1) was determined according to the Arnon’s method [44]. The 10 fully expanded and intact acquired from the previous sampling activities and were weighed until we obtained about 0.1 g with a precision balance of α = 0.0001 g.
The leaf pigments were extracted in 80% acetone and leaf was pulverized completely by using the homogenizer Bio-Gen PRO200 (PRO Scientific Inc., Oxford, MS, USA) to obtain a completely white foliar tissue. Before the analysis, the leaf pigments were centrifuged for 5 min at 7 rpm to separate solid and liquid substance, and then more solution of 80% acetone was added until the achievement of a final volume of 25 mL. The determination of leaves chlorophyll concentration (mg g−1) is carried out by using a spectrophotometer Varian Cary 50 Scan ultraviolet (UV)-visible spectrophotometric (Agilent Technologies, Inc., Santa Clara, CA, USA), where the absorbance’s (A) were measured at 663 and 645 nm. Total chlorophyll contents of each sample were computed from the follow equation Equation (1):
mg   total   chlorophyll / g   foliar   tissue = 20.2   ·   A 645   +   8.02   ·   A 663   ·   V 1000   ·   W
where:
  • A = absorbance at specific wavelengths
  • V = final volume of chlorophyll extract
  • W = fresh weigh of tissue extracted.

2.3.3. Nitrogen Crop Status Determination

At each test area, fresh plant biomass was manually cut and collected in a 0.5 m long-row. The fresh plant biomass was oven-dried at 80 °C for 48 h and weighed (g). Before analyzing for the N content (%), the dry biomass was grounded to pass a 0.5 mm sieve.
The N content (%) was determined by the automated combustion analysis Dumas method [45,46] in an oxygen-enriched atmosphere at a high temperature by using the EA 1110 LECO CHNS-0 analyser (Leico Corporation, St. Joseph, MI, USA).
Starting from the N content (%) results, the NNI was calculated by dividing actual nitrogen concentration by the critical nitrogen concentration using the wheat dilution curve [47,48]. The following equations were used to calculate the NNI Equations (2) and (3):
Nitrogen   Nutrition   Index   NNI =   % N N c
Critical   Nitrogen   Concentration   ( N c ) = a   ·   DM b
where:
  • %N = actual N concentration
  • Nc = critical N concentration
  • a = 5.35
  • DM = dry matter (g)
  • b = 0.442.

2.3.4. Unmanned Aerial Vehicle (UAV) Flights, Image Acquisition and Processing

Multispectral images were acquired in correspondence of the three phenological stages (ZS22, ZS35 and ZS60).
Before each flight, a radiometric calibration with the calibration reflectance panel was performed. In addition, in order to perform better imagery acquisition, the UAV platform was equipped with the incident light sensor (ILS) which can record the variation of the light reflectance during the flight for each 9 spectral bands (wavelengths between 443 and 865 nm) of the MAIA Sentinel 2.
To obtain a high spatial accuracy of the remote-sensing images, 6 GCPs were geo-referenced before each flight operation were used. All the relevant information regarding the flights is given in the following Table 4.
Each image acquired by UAV flight operations required an image process to analyze all the relevant information. The image processing is composed by four main steps: raw to tiff conversion; orthomosaic reflectance map generation; selection and calculation of VIs; data extraction.
After exporting the raw images from the camera, the dedicated software MultiCam Stitcher Pro (SAL engineering and Eoptis, Russi, Italy) was used to convert the raw format image to the multi-channel tiff format.
To generate the orthomosaic reflectance maps the commercial software Agisoft Metashape (Agisoft LLC, St. Petersburg, Russia) was used. The software is based on structure from motion (SfM) algorithm [49]. The newly generated orthomosaic reflectance map was imported in R statistical software [50] and used to complete the remaining two main steps of the image processing.
To select the most relevant VI based on multispectral imagery for precision agriculture application, ten VIs were compared [51]. The VIs analysed in this study are reported in Table 5.
The csv file format contained the georeferenced positions of the biomass sampling was imported in R to extract the vegetation index values. To perform the values extraction, the “extract” function of “raster” R package [58] was used.

2.3.5. Yield Parameters

To characterize the yield for each of compared treatments, the number of kernels per spike (KS), the number of spikes m−2 and the grain yield (t ha−1) were measured at crop maturity (ZS92).
The KS was estimated on 30 spikes randomly collected per plot while the spikes m−2 was determined by counting the number of spikes along 1-m long row. The grain yield (t ha−1) expressed in dry matter, was measured by using a laboratory thresher for the three test areas (1 m long-row) per subplot.

2.4. Statistical Analysis

All statistical analyses were performed with R statistical software [50].
Every good analysis is based on a good model that is biologically relevant and based on realistic assumptions [59].
To select the model that can fit better the crop parameters and vegetation index dataset, three different models by using the “stats” [60] and “nlme” [61] R packages were built:
(1)
A one factor linear model was built by using the “lm” function on which the soil management was considered the main factor.
(2)
A full factorial linear model was also created by using the “lm” function on which the soil management, nitrogen input and year were considered as factors.
(3)
A mixed model was built by using “lme” function on which the soil management, nitrogen input and year were considered as fixed factors, and the Zadoks scale, the block and subplot were considered as random factor, according with [59].
In order to evaluate the best model, an ‘a posteriori’ selection was made, based on those statistics which put a penalty on ‘complexity’, such as the Akaike information criterion (AIC) [62] the Bayesian information criterion (BIC) [63], and likelihood ratio tests (LRTs) [64].
Based on the previous model analyses, the statistical model that fit better the crop parameters, vegetation index and the yield datasets was the mixed model, as reported by several studies [59,65].
Before performing the analysis of variance (ANOVA) it has been verified that the model met the three assumptions of the ANOVA. The Normality distribution of the model residual was checked both graphically (QQ-plot) and by performing the Shapiro-Wilk normality test. Moreover, the homoscedasticity was checked using the Levene test. The last ANOVA assumption was satisfied by the experimental design and the random sampling.
When all the three ANOVA assumptions were met, The ANOVA was applied to the model. Only when the ANOVA showed a significant difference between the factors (p-value < 0.05) the estimated marginal means post-hoc analysis was performed by using the “emmeans” function with the Bonferroni adjustment of the emmeans R package [66].
The linear relationships between the previous variables were analyzed, by using the Pearson correlation parameter (r) of the “cor” function of the R stats package [50], the coefficient of determination (R2) and the root mean square error (RMSE).
To analyze the yield dataset the same statistical approach it was used as previously proposed for crop parameters and vegetation index dataset, but unlike before we set the year as a unique random factor while the soil management and nitrogen input were set as fixed factors.
The soil data were analyzed by using a full factorial linear model on which soil management and nitrogen input were considered as factors. The assumptions of the ANOVA application were checked and applied to the model. The investigation of the significant differences between groups factors were performed as before.

3. Results

3.1. Crop Parameters

The ANOVA applied to the mixed model showed that the soil management (SM) was statistically significant for the NNI and SPAD readings while the nitrogen input (N) was statistically significant for each crop parameters (Table 6).
The interaction between SM and N was statistically significant for all the measured variables (Table 6); this led to perform the estimated marginal means post hoc analysis on this factorial interaction.
The all pair-wise comparisons test showed that the factorial combination CT-N2 and NT-N2 performed statistically better than the other factorial combinations for the leaf chlorophyll concentration (mg g−1), %N, NNI and SPAD readings (Table 7) with a mean value of 2.73 (mg g−1), 2.18 (%), 1.19 and 47.94 respectively.
The factorial combination CT-N1 and NT-N1 showed a statistical superiority over both unfertilized factorial combinations (CT-N0 and NT-N0) for each variable. This confirmed that confirmed that the nitrogen input is the key driver for the durum wheat growth development (Table 7). In detail, it was observed that there was an average percentage increase of +25.10%, +18.04%, +45.05%, +30.52% respectively for the leaf chlorophyll concentration (mg g−1), %N, NNI and SPAD readings.
NT-N0 reported a statistically higher value for the leaf chlorophyll concentration (mg g−1), %N, NNI and SPAD readings than the factorial combination CT-N0 (Table 7), with a difference of +0.44, +0.29%, +0.16 and +8.97 mg g−1, respectively.

3.2. Vegetation Index

The ANOVA applied to the mixed model showed that the SM was statistically significant only for 4 out of the 10 VIs analysed. These indices were the atmospherically resistant vegetation index (ARVI), green atmospherically resistant vegetation index (GARI), normalized difference red edge index (NDRE) and NDVI (Table 8). The N factor had significant impact on all 10 VIs analysed with a p-value < 0.001.
The interaction between SM and N was statistically significant only for the NIR-based VIs.
The all pair-wise comparisons test confirmed that the factorial combination CT-N2 and NT-N2 showed statistically higher values for the ARVI, modified soil-adjusted vegetation index 2 (MSAVI2), GARI, NDRE and NDVI than the other factorial combinations with a mean value of 0.63, 0.56, 0.89, 0.40 and 0.72, respectively (Table 9).
The factorial combination CT-N1 and NT-N1 performed statistically better than both unfertilized factorial combinations for each VI with an average difference of +0.21, +0.22, +0.14, +0.13, +0.17 for ARVI, MSAVI2, GARI, NDRE and NDVI, respectively.
NT-N0 reported a higher value than CT-N0 for those VIs that has the NIR band for the formula calculation, with a percentage increase of +33.33%, +46.88%, +14.86%, +27.27% and +23.08% for ARVI, MSAVI2, GARI, NDRE and NDVI, respectively.
Those VIs analysed, which do not have the NIR band on the formula calculation, reported a statistical difference only between the N levels (Table 10). Both fertilized treatments showed a higher value than the unfertilized treatment for GLI and VDVI with a mean difference of +0.09 and +0.09 respectively. The VARI, VARI Green and Vig showed a superiority of N2, passing through the N1, to the N0 which showed the lowest values.

3.3. Regression Analysis

Both remote (through vegetation indices computation) and proximal sensing (through SPAD readings) data acquisition showed a significant linear relationship with leaves chlorophyll concentration (mg g−1) and NNI. The mean correlation parameter between all the remote and proximal sensing and leaves chlorophyll concentration (mg g−1) was R2 = 0.47 and RMSE = 0.60. Furthermore, the average of the correlation parameters between the remote and proximal sensing with NNI was R2 = 0.68 and RMSE = 0.31(Table 11).
The MSAVI2 vegetation index showed a better correlation to be related to leaves chlorophyll concentration (mg g−1) with a R2 value of 0.68 and RMSE equal to 0.47; for the NNI the NDRE vegetation index (Figure 2) showed a better performance with a R2 of 0.85 and RMSE of 0.21.
The NIR-based VIs showed a better average correlation parameter for both crop parameters compared to the not NIR-based VIs (Table 12).
The linear relationship between leaves chlorophyll concentration (mg g−1) and the VIs with the NIR band obtained a higher correlation parameter than the VIs without NIR band, with a difference of +0.32 for R2 and −0.18 of RMSE (Table 12). A lower difference of +0.04 for R2 and −0.03 RMSE for NNI was reported between VIs with NIR band and the VIs without the NIR band (Table 12).

3.4. Yield Parameters

The ANOVA showed a statistically significant impact of soil management (SM) on the number of kernels per spike (KS) and grain yield (t ha−1). The nitrogen input (N) affected significantly all the yield components and grain yield (t ha−1). The factorial combination SM × N showed a significant difference only on the KS and grain yield (t ha−1) (Table 13).
With regard to the KS and grain yield (t ha−1) the estimated marginal means showed the same statistical difference as the VIs with the NIR band on the formula calculation, on which the NT-N2 and CT-N2 were statistically higher than the other factorial combinations (Table 14) with a mean value of 35 KS and 5.28 grain yield (t ha−1).
The factorial combination CT-N1 and NT-N1 performed statistically better than both unfertilized factorial combinations, with a mean difference of +17 KS and +1.87 grain yield (t ha−1) (Table 14).
The KS and grain yield (t ha−1) observed on the NT-N0 combination were statistically higher than the CT-N0 (Table 14), with a percentage increase of +41.18% and +41.83% respectively.
The number of spike m−2 reported a difference only between the three different nitrogen levels where the higher mean value was obtained by N2 with 567 spike m−2, followed by N1 with 530 spike m−2 and N0 with 400 spike m−2 (Table 14).

4. Discussion

4.1. Crop Parameters

Nitrogen input is the key driver for the durum wheat growth development [67]. As the nitrogen input increase, the overall crop nutritional status improves. This allowed a higher production of leaves’ chlorophyll concentration (mg g−1), biomass nitrogen content (%), NNI and the SPAD readings as previously observed for durum wheat by [68] and [69].
On the unfertilized treatment, the NT performed better for each crop parameters than the CT. This is probably due to: the higher content of organic matter and consequently nitrogen availability present in the NT system 0–20 cm soil layer [70] than the compared CT plots (Table 2).
The positive effect of CA compared with the conventional tillage practices was reported also by [71] who, by using the Decision Support System for Agrotechnology Transfer model (DSSAT) in the same experimental site, showed that NT leads to an increase in soil organic carbon with a rate of 0.43, 0.31 and 0.03 t ha−1 respectively for 180, 90 and 0 kg N ha−1 year−1.
This characteristic of NT plots was observed also by [72] and [73] who recommended to adopt the CA to increase the soil organic matter, especially to mitigate the variations due to climate change [74].
Recently, [75] demonstrated that the use of the ARMOSA process-based crop model that CA can increase the annual soil organic carbon sequestration by 0.4% year−1.

4.2. Vegetation Index

All VIs were susceptible to the different nitrogen input provided, therefore they are good indicators to identify crop variations due to different nitrogen input as observed by [76] where a significant difference was reported on the SPAD readings and NDVI values between 0, 80 and 160 kg N ha−1. Moreover, [77] observed a significant difference on the NDVI values between 0, 55, 110 and 220 kg N ha−1.
The ability of individual vegetation indices to detect different levels of nitrogen fertilization depends on the changes in reflectance recorded by each individual multispectral band. It is widely accepted that when a crop has a better overall nutritional status led to an higher leaves chlorophyll concentration (mg g−1) production [69]. A higher leaf chlorophyll concentration led to higher absorbance in the red and blue visible spectrum (VIS) area and reflect the green VIS area. This spectral characteristic of chlorophyll allowed the development of the Arnon method 1949 [44] with which we measured the leaves’ chlorophyll concentration (mg g−1) and also the development of the optics that we used in this study such as the SPAD readings, multispectral camera and others proximal and remote-sensing tools applied in precision agriculture.
The VIs calculated with the NIR band detected the variation due to the SM within the N treatments, while the others VIs could not report a statistically alteration due the SM × N factorial interaction. The ARVI, MSAVI2, GARI, NDRE and NDVI statistically detected the superiority of NT-N0 over the CT-N0 as observed by [7] which in cotton observed a statistically increase of the NDVI on NT compared to CT.

4.3. Relationship Between SPAD Readings, Vegetation Indices and Crop Parameters

All significant linear relationships between the non-destructive methods and crop parameters were observed. It was reported by [78] that there was a significant linear relationship between 12 VIs and the nitrogen status in an experimental design where 8 different levels of N fertilization were applied on cotton (0, 130, 177, 194, 210, 307, 324 and 340 kg N ha−1).
A significant and high linear relationship between 30 VIs and leaf chlorophyll content (LCC) was reported by [79], where the MDATT index reached a correlation accuracy of R2 = 0.86 and a RMSE of 6.85.
The VI that reported the best correlation with the leaf chlorophyll concentration (mg g−1) was the MSAVI2, which was the same VI to obtain a significant relationship in each factorial combination reported by [68].
The NDRE was the VI to show the best relationship with NNI (Table 11), in accordance with [78], where for three different cotton growth stages NDRE was the best VI to be related with the nitrogen content.
The VIs which use the NIR band showed a higher mean correlation parameter with measured crop parameters than the VIs which do not use the NIR band on the formula calculation. This is probably due to the higher reflectance of the NIR in the plant cell, and in fact it is widely accepted that a plant with a good nutritional status has a higher reflectance in NIR spectral zone than a plant which does not have a good nutritional status [80].

4.4. Yield Components

The nitrogen input affected all the yield components, and this confirms what reported by [17] where the nitrogen was the key driver for the durum wheat growth and production after analyzing 20 years of data.
For the soil management the results found in the bibliography are discordant. [70] where only the unfertilized combination between CT and NT have been considered, reported the superiority of NT over the CT for KS and grain yield (t ha−1). While [17] reported that the NT did not provide any substantial advantage or disadvantage to durum wheat grain yield in comparison to CT. These results encourage in any case the adoption of the NT management because it leads to higher soil organic matter and nitrogen availability as described by [75] and allows farmers to have lower budget cost per ha as reported by [81].
The spike m−2 was affected by the different nitrogen levels as observed by [68]. The higher grain yield (t ha−1) obtained from NT-N0 over CT-N0 was +41.83% which is in line with the results obtained for KS (+41.18%).
The superiority of NT over CT is due to the higher organic matter (g kg−1). total nitrogen availability (g kg−1) and carbon nitrogen that NT plots had respect to CT (Table 1), because of over 25 years of no soil-tillage practices [73].

5. Conclusions

Over 25 years of repeated no tillage practices allowed to NT plots to accumulate a higher soil organic matter content and total nitrogen availability than the CT plots.
In order to obtain valid results from a split plot LTE and interpret them correctly, the use of the mixed model is recommended, because it is the statistical model capable of better fitting the experimental data than the classics linear models.
Nitrogen input is a key driver for durum wheat growth development and production, as the nitrogen input increases. There is an increase also in the crop parameters analyzed. When no mineral nitrogen was supplied to the crop, NT performed better than CT for any parameters measured in this study apart from the spike m−2.
All VIs were capable of detecting variation due to the nitrogen input, but few of them detected the difference due to the interaction between soil management and nitrogen input, which were the VIs with the NIR band on the formula calculation. Moreover, we recommend readers use the MSAVI2 to monitor the leaf chlorophyll concentration while we suggest using the NDRE to monitor the nitrogen durum wheat status.
The presented results can also provide an important contribution in the development and implementation of agro-environmental policies aimed at containing production inputs to reduce the environmental impact of cereal-based cropping systems.

Author Contributions

Conceptualization, R.O.; Data Curation, R.O., M.F. and S.Z.; Formal Analysis, M.F. and S.Z.; Project Administration, R.O.; Supervision, R.O. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing not applicable.

Acknowledgments

Marche Polytechnic University, trough University Strategic Project (PSA 2017) “PFRLab: Setting of a Precision Farming Robotic Laboratory for cropping system sustainability and food safety and security” and Marche Region trough measure 16.1 Action 2 of the Rural Development Plan (PSR 2014/2020) “Precision agriculture: reduction of the environmental impact of production systems” (project ID 29000).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

AbbreviationsDefinitions
NT No tillage
CT Conventional tillage
N Nitrogen
CA Conservative agriculture
PA Precision agriculture
UAV Unmanned aerial vehicle
LTE Long-term experiment
N0 0 kg N ha−1
N1 90 kg N ha−1
N2 180 kg N ha−1
NNI Nitrogen Nutrition Index
ZS Zadoks scale
KS Kernels spike
GCP Ground control points
A Absorbance at specific wavelenghts
Nc Critical nitrogen concentration
DM Dry matter
ILS Incident light sensor
VI Vegetation index
VIs Vegetation indexes
ARVI Atmospherically resistant vegetation index
GARI Green atmospherically resistant vegetation index
MSAVI2 Modified Soil-adjusted Vegetation Index 2
NDRE Normalized difference red edge index
NDVI Normalized difference vegetation index
GLI Green leaf index
VARI Visible atmospherically resistant index
VARI Green Visible atmospherically resistant index green
VDVI Visible-band difference vegetation index
Vig Vegetation index green
AIC Akaike information criterion
BIC Bayesian information criterion
LRT Likelihood ratio tests
ANOVA Analysis of variance
R2 Coefficient of determination
RMSE Root mean square error
SM Soil management
Y Year

References

  1. FAO (Food and Agriculture Organization of the United Nations). Crop Prospects and Food Situation. Quarterly Global Report. September, 2019; FAO (Food and Agriculture Organization of the United Nations): Rome, Italy, 2019. [Google Scholar]
  2. Tsialtas, J.T.; Theologidou, G.S.; Karaoglanidis, G.S. Effects of pyraclostrobin on leaf, leaf physiology, yield and quality of durum wheat under Mediterranean conditions. J. Crop Prot. 2018, 113, 48–55. [Google Scholar] [CrossRef]
  3. Pagnani, G.; Galieni, A.; D’Egidio, S.; Visioli, G.; Stagnari, F.; Pisante, M. Effect of Soil Tillage and Crop Sequence on Grain Yield and Quality of Durum Wheat in Mediterranean Areas. Agronomy 2019, 9, 488. [Google Scholar] [CrossRef] [Green Version]
  4. Araus, J.L.; Slafer, G.A.; Reynolds, M.P.; Royo, C. Plant breeding and drought in C3 cereals: What should we breed for? Ann. Bot. 2002, 89, 925–940. [Google Scholar] [CrossRef]
  5. Giorgi, F.; Lionello, P. Climate change projections for the Mediterranean region. Glob. Planet. Chang. 2008, 63, 90–104. [Google Scholar] [CrossRef]
  6. UN (United Nations) Reports. World Population Prospects: The 2017 Revision; Working Paper No. ESA/P/WP.250; Dep. of Economic and Social Affairs, Population Division: New York, NY, USA, 2017. [Google Scholar]
  7. Ashapure, A.; Jung, J.; Yeom, J.; Chang, A.; Maeda, M.; Maeda, A.; Landivar, J. A novel framework to detect conventional tillage and no-tillage cropping system effect on cotton growth and development using multi-temporal UAS data. Remote Sens. Environ. 2019, 152, 49–64. [Google Scholar] [CrossRef]
  8. Anderson, W.K.; Impiglia, A. Management of Dryland Wheat. In Bread Wheat Improvement and Production; Curtis, B.C., Rajaram, S., Macpherson, H.G., Eds.; Food and Agriculture Organization of the United Nations: Rome, Italy, 2002; pp. 407–432. [Google Scholar]
  9. Pisante, M.; Stagnari, F.; Acutis, M.; Bindi, M.; Brilli, L.; Di Stefano, V.; Carozzi, M. Conservation agriculture and climate change. In Conservation Agriculture; Farooq, M., Siddique, K., Eds.; Springer International Publishing: Cham, Switzerland, 2015; pp. 579–620. [Google Scholar]
  10. Hobbs, P.R. Conservation agriculture: What is it and why is it important for future sustainable food production? J. Agric. Sci. 2007, 145, 127. [Google Scholar] [CrossRef] [Green Version]
  11. Farooq, M.; Flower, K.; Jabran, K.; Wahid, A.; Siddique, K.H. Crop yield and weed management in rainfed conservation agriculture. Soil Tillage Res. 2011, 117, 172–183. [Google Scholar] [CrossRef]
  12. Pisante, M.; Stagnari, F.; Grant, C.A. Agricultural innovations for sustainable crop production intensification. Ital. J. Agron. 2012, 7, e40. [Google Scholar] [CrossRef] [Green Version]
  13. Stagnari, F.; Jan, S.; Angelica, G.; Pisante, M. Sustainable agricultural practices for water quality protection. In Water Stress and Crop Plants: A Sustainable Approach; Ahmad, P., Ed.; John Wiley & Sons: Chichester, UK, 2016; pp. 75–85. [Google Scholar]
  14. López-Garrido, R.; Madejón, E.; Murillo, J.M.; Moreno, F. Short and long-term distribution with depth of soil organic carbon and nutrients under traditional and conservation tillage in a Mediterranean environment (southwest Spain). Soil Use Manag. 2011, 27, 177–185. [Google Scholar] [CrossRef]
  15. Amato, G.; Ruisi, P.; Frenda, A.S.; Di Miceli, G.; Saia, S.; Plaia, A.; Giambalvo, D. Long-term tillage and crop sequence effects on wheat grain yield and quality. Agron. J. 2013, 105, 1317–1327. [Google Scholar] [CrossRef]
  16. Mazzoncini, M.; Di Bene, C.; Coli, A.; Antichi, D.; Petri, M.; Bonari, E. Rainfed wheat and soybean productivity in a long-term tillage experiment in central. Italy. Agron. J. 2008, 100, 1418–1429. [Google Scholar] [CrossRef]
  17. Seddaiu, G.; Iocola, I.; Farina, R.; Orsini, R.; Iezzi, G.; Roggero, P.P. Long-term effects of tillage practices and N fertilization in rainfed Mediterranean cropping systems: Durum wheat, sunflower and maize grain yield. Eur. J. Agron. 2016, 77, 166–178. [Google Scholar] [CrossRef]
  18. Mulla, D.J. Twenty-five years of remote sensing in precision agriculture: Key advances and remaining knowledge gaps. Biosyst. Eng. 2013, 114, 358–371. [Google Scholar] [CrossRef]
  19. Georgi, C.; Spengler, D.; Itzerott, S.; Kleinschmit, B. Automatic delineation algorithm for site-specific management zones based on satellite remote sensing data. Precis. Agric. 2018, 19, 684–707. [Google Scholar] [CrossRef] [Green Version]
  20. Thilakarathna, M.; Raizada, M. Challenges in Using Precision Agriculture to Optimize Symbiotic Nitrogen Fixation in Legumes: Progress, Limitations, and Future Improvements Needed in Diagnostic Testing. Agronomy 2018, 8, 78. [Google Scholar] [CrossRef] [Green Version]
  21. Deery, D.; Jimenez-Berni, J.; Jones, H.; Sirault, X.; Furbank, R. Proximal remote sensing buggies and potential applications for field-based phenotyping. Agronomy 2014, 5, 349–379. [Google Scholar] [CrossRef] [Green Version]
  22. Sankaran, S.; Khot, L.R.; Espinoza, C.Z.; Jarolmasjed, S.; Sathuvalli, V.R.; Vandemark, G.J.; Miklas, P.N.; Carter, A.H.; Pumphrey, M.O.; Knowles, N.R.; et al. Low-altitude, high-resolution aerial imaging systems for row and field crop phenotyping: A review. Eur. J. Agron. 2015, 70, 112–123. [Google Scholar] [CrossRef]
  23. Zhang, C.; Kovacs, J.M. The application of small unmanned aerial systems for precision agriculture: A review. Precis. Agric. 2012, 13, 693–712. [Google Scholar] [CrossRef]
  24. Gago, J.; Douthe, C.; Coopman, R.E.; Gallego, P.P.; Ribas-Carbo, M.; Flexas, J.; Escalona, J.; Medrano, H. UAVs challenge to assess water stress for sustainable agriculture. Agric. Water. Manag. 2015, 153, 9–19. [Google Scholar] [CrossRef]
  25. Son, N.T.; Chen, C.F.; Chen, C.R.; Minh, V.Q.; Trung, N.H. A comparative analysis of multitemporal MODIS EVI and NDVI data for large-scale rice yield estimation. Agric. For. Meteorol. 2014, 197, 52–64. [Google Scholar] [CrossRef]
  26. Zaman-Allah, M.; Vergara, O.; Araus, J.L.; Terekegne, A.; Magorokosho, C.; Zarco-Tejada, P.J.; Hornero, A.; Albà, H.; Das, B.; Craufurd, P.; et al. Unmanned aerial platform-based multi-spectral imaging for field phenotyping of maize. Plant Methods. 2015, 11, 1–10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Maimaitijiang, M.; Sagan, V.; Sidike, P.; Daloye, A.M.; Erkbol, H.; Fritschi, F.B. Crop Monitoring Using Satellite/UAV Data Fusion and Machine Learning. Remote Sens. 2020, 12, 1357. [Google Scholar] [CrossRef]
  28. Reynolds, M.; Tuberosa, R. Translational research impacting on crop productivity in drought-prone environments. Curr. Opin. Plant Biol. 2008, 11, 171–179. [Google Scholar] [CrossRef] [PubMed]
  29. Khot, L.R.; Sankaran, S.; Carter, A.H.; Johnson, D.A.; Cummings, T.F. UAS imaging-based decision tools for arid winter wheat and irrigated potato production management. Int. J. Remote Sens. 2016, 37, 125–137. [Google Scholar] [CrossRef]
  30. Dellinger, A.; Schmidt, J.; Beegle, D. Developing Nitrogen Fertilizer Recommendations for Corn Using an Active Sensor. Agron. J. 2008, 100, 1546–1552. [Google Scholar] [CrossRef]
  31. Barker, D.; Sawyer, J. Using Active Canopy Sensors to Quantify Corn Nitrogen Stress and Nitrogen Application Rate. Agron. J. 2010, 102, 964–971. [Google Scholar] [CrossRef]
  32. Schmidhalter, U.; Glas, J.; Heigl, R.; Manhart, R.; Wiesent, S.; Gutser, R.; Neudecker, E. Application and testing of a crop scanning instrument-Field experiments with reduced crop width, tall maize plants and monitoring of cereal yield. In Proceedings of the 3rd European Conference on Precision Agriculture, Montpellier, France, 18–20 June 2001; pp. 953–958. [Google Scholar]
  33. Thoren, D.; Schmidhalter, U. Nitrogen status and biomass determination of oilseed rape by laser-induced chlorophyll fluorescence. Eur. J. Agron. 2009, 30, 238–242. [Google Scholar] [CrossRef]
  34. Winterhalter, L.; Mistele, B.; Jampatong, S.; Schmidhalter, U. High-Throughput Sensing of Aerial Biomass and Above-Ground Nitrogen Uptake in the Vegetative Stage of Well-Watered and Drought Stressed Tropical Maize Hybrids. Crop Sci. 2011, 51, 479–489. [Google Scholar] [CrossRef]
  35. Erdle, K.; Mistele, B.; Schmidhalter, U. Spectral high-throughput assessments of phenotypic differences in biomass and nitrogen partitioning during grain filling of wheat under high yielding Western European conditions. Field Crops Res. 2013, 141, 16–26. [Google Scholar] [CrossRef]
  36. Raun, W.R.; Solie, J.B.; Johnson, G.V.; Stone, M.L.; Lukina, E.V.; Thomason, W.E.; Schepers, J.S. In-season prediction of potential grain yield in winter wheat using canopy reflectance. Agron. J. 2001, 93, 131–138. [Google Scholar] [CrossRef] [Green Version]
  37. Magney, T.S.; Eitel, J.U.H.; Huggins, D.R.; Vierling, L.A. Proximal NDVI derived phenology improves in-season predictions of wheat quantity and quality. Agric. For. Meteorol. 2016, 217, 46–60. [Google Scholar] [CrossRef]
  38. Morgounov, A.; Gummadov, N.; Belen, S.; Kaya, Y.; Keser, M.; Mursalova, J. Association of digital photo parameters and NDVI with winter wheat grain yield in variable environments. Turk. J. Agric. For. 2014, 38, 624–632. [Google Scholar] [CrossRef]
  39. Lopresti, M.F.; Di Bella, C.M.; Degioanni, A.J. Relationship between MODIS-NDVI data and wheat yield: A case study in northern Buenos Aires province, Argentina. Inf. Process. Agric. 2015, 2, 73–84. [Google Scholar] [CrossRef] [Green Version]
  40. Mekliche, A.; Hanifi-Mekliche, L.; Aïdaoui, A.; Gate, P.H.; Bouthier, A.; Monneveux, P.H. Grain yield and its components study and their association with normalized difference vegetation index (NDVI) under terminal water deficit and well-irrigated conditions in wheat (Triticum durum Desf. and Triticum aestivum L.). Afr. J. Biotechnol. 2015, 14, 2142–2148. [Google Scholar] [CrossRef]
  41. Orsini, R.; Basili, D.; Belletti, M.; Bentivoglio, D.; Bozzi, C.A.; Chiappini, S.; Conti, C.; Fiorentini, M.; Galli, A.; Giorgini, E.; et al. Setting of a precision farming robotic laboratory for cropping system sustainability and food safety and security: Preliminary results. IOP Conf. Ser. Earth Environ. Sci. 2019, 275, 1–8. [Google Scholar] [CrossRef]
  42. FAO (Food and Agriculture Organization of the United Nations). World Soil Resources Report 103; IUSS, ISRIC; FAO (Food and Agriculture Organization of the United Nations): Rome, Italy, 2006. [Google Scholar]
  43. Zadoks, J.C.; Chang, T.T.; Konzak, C.F. A decimal code for the growth stages of cereals. Weed Res. 1974, 14, 415–421. [Google Scholar] [CrossRef]
  44. Arnon, D.J. Copper enzymes in isolated chloroplasts. Polyphenoloxidase in Beta vulgaris. Plant Physiol. 1949, 24, 1. [Google Scholar] [CrossRef] [Green Version]
  45. Dumas, A. No Title. Ann. Chim. 1826, 33, 342. [Google Scholar]
  46. Buckee, G.K. Determination of total nitrogen in Barley, Malt and Beer by Kjeldahl procedures and the Dumas combustion method. J. Inorg. Biochem. 1994, 100, 57–64. [Google Scholar]
  47. Justes, E.; Mary, B.; Meynard, J.M.; Machet, J.M.; Thelier-Huche, L. Determination of critical nitrogen dilution curve for winter wheat crops. Ann. Bot. 2006, 74, 397–407. [Google Scholar] [CrossRef]
  48. Debaeke, P.; Justes, E. Relationship between the Normalized SPAD Index and the Nitrogen Nutrition Index: Application to durum wheat. J. Plant Nutr. 2006, 29, 75–92. [Google Scholar] [CrossRef]
  49. Verhoeven, G. Taking computer vision aloft archaeological three-dimensional reconstructions from aerial photographs with photoscan. Archaeol. Prospect. 2011, 18, 67–73. [Google Scholar] [CrossRef]
  50. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2019. [Google Scholar]
  51. Xue, J.; Su, B. Significant Remote Sensing Vegetation Indeces: A Review of Developments and Applications. J. Sens. 2017, 2017, 1353691. [Google Scholar] [CrossRef] [Green Version]
  52. Korhonen, L.; Ali-Sisto, D.; Tokola, T. Tropical forest canopy cover estimation using satellite imagery and airborne lidar reference data. Silva Fenn. 2015, 49, 1–18. [Google Scholar] [CrossRef]
  53. Gitelson, A.; Stark, R.; Gritis, U.; Rundquist, D.; Kaufman, Y.; Derry, D. Vegetation and Soil Lines in Visible Spectral Space: A Concept and Technique for Remote Estimation of Vegetation Fraction. Int. J. Remote Sens. 2002, 23, 2537–2562. [Google Scholar] [CrossRef]
  54. Leprieur, C.; Kerr, Y.H.; Mastorchio, S.; Meunier, J.C. Monitoring vegetation cover across semi-arid regions: Comparison of remote observations from various scales. Int. J. Remote Sens. 2000, 21, 281–300. [Google Scholar] [CrossRef]
  55. Barnes, E.M.; Clarke, T.R.; Richards, S.E.; Colaizzi, P.D.; Haberland, J.; Kostrzewski, M.; Waller, P. Coincident detection of crop water stress, nitrogen status and canopy density using ground based multispectral data. In Proceedings of the Fifth International Conference on Precision Agriculture, Bloomington, MN, USA, 16–19 July 2000. [Google Scholar]
  56. Rouse, J.; Haas, R.J.S.; Deering, D. Monitoring Vegetation Systems in the Great Plains with ERTS. In Proceedings of the Third ERTS Symposium, NASA, Washington, DC, USA, 10–14 December 1973; Volume 351, pp. 309–317. [Google Scholar]
  57. Wang, X.; Wang, M.; Wang, S.; Wu, Y. Extraction of vegetation information from visible unmanned aerial vehicle images. Nongye Gongcheng Xuebao/Trans. Chin. Soc. Agric. Eng. 2015, 31, 152–159. [Google Scholar]
  58. Raster: Geographic Data Analysis and Modeling, CRAN. Available online: https://cran.r-project.org/web/packages/raster/index.html (accessed on 14 November 2020).
  59. Onofri, A.; Seddaiu, G.; Piepho, H.P. Long-Term Experiments with cropping systems: Case studies on data analysis. Eur. J. Agron. 2016, 77, 223–235. [Google Scholar] [CrossRef]
  60. Chambers, J.M. Linear models. In Statistical Models in S; Chambers, J.M., Hastie, T.J., Eds.; Wadsworth & Brooks/Cole: Belmont, CA, USA, 1992; Chapter 4. [Google Scholar]
  61. Venables, W.N.; Ripley, B.D. Modern Applied Statistics with S, 4th ed.; Springer: Berlin/Heidelberg, Germany, 2002. [Google Scholar]
  62. Akaike, H. A new look at the statistical model identification. IEEE Trans. Autom. Control 1974, 19, 716–723. [Google Scholar] [CrossRef]
  63. Schwarz, G. Estimating the dimension of a model. Ann. Stat. 1978, 6, 461–464. [Google Scholar] [CrossRef]
  64. Wilks, S.S. The large-sample distribution of the likelihood ratio for testing composite hypotheses. Ann. Math. Stat. 1938, 9, 60–62. [Google Scholar] [CrossRef]
  65. Piepho, H.P.; Büchse, A.; Richter, C. A mixed modelling approach for randomized experiments with repeated measures. J. Agron. Crop Sci. 2004, 190, 230–247. [Google Scholar] [CrossRef]
  66. Emmeans: Estimated Marginal Means, Aka Least-Squares Means, CRAN. Available online: https://cran.r-project.org/web/packages/emmeans/index.html (accessed on 9 December 2020).
  67. Acevedo, P.; Silva, H. Wheat Growth and Physiology; F.A.O. Corporate Repository: Rome, Italy, 2009; pp. 1–24. [Google Scholar]
  68. Orsini, R.; Fiorentini, M.; Zenobi, S. Testing vegetation index categories as influenced by soil management and nitrogen fertilization in cereal based cropping systems. In Proceedings of the Workshop International IEEE 2019 SU, Metrology for Agriculture and Foresty, Portici, Italy, 24–26 October 2019. [Google Scholar]
  69. Fiorentini, M.; Zenobi, S.; Giorgini, E.; Basili, D.; Conti, C.; Pro, C.; Monaci, E.; Orsini, R. Nitrogen and chlorophyll status determination in durum wheat as influenced by fertilization and soil management: Preliminary results. PLoS ONE 2019, 14, e0225126. [Google Scholar] [CrossRef] [PubMed]
  70. Orsini, R.; Fiorentini, M.; Zenobi, S. Evaluation of Soil Management Effect on Crop Productivity and Vegetation Indices Accuracy in Mediterranean Cereal-Based Cropping Systems. Sensors 2020, 20, 3383. [Google Scholar] [CrossRef] [PubMed]
  71. De Sanctis, G.; Roggero, P.P.; Seddaiu, G.; Orsini, R.; Porter, C.H.; Jones, J.W. Long-term no tillage increased soil organic carbon content of rain-fed cereal systems in a Mediterranean area. Eur. J. Agron. 2012, 40, 18–27. [Google Scholar] [CrossRef]
  72. Pastorelli, R.; Vignozzi, N.; Landi, S.; Piccolo, R.; Orsini, R.; Seddaiu, G.; Roggero, P.P.; Pagliai, M. Consequences on macroporosity and bacterial diversity of adopting a no-tillage farming system in a clayish soil of Central Italy. Soil Biol. Biochem. 2013, 66, 78–93. [Google Scholar] [CrossRef]
  73. Iocola, I.; Bassu, S.; Farina, R.; Antichi, D.; Basso, B.; Bindi, M.; Dalla Marta, A.; Danuso, F.; Doro, L.; Ferrise, R.; et al. Can conservation tillage mitigate climate change impacts in Mediterranean cereal systems? A soil organic carbon assessment using long term experiments. Eur. J. Agron. 2017, 90, 96–107. [Google Scholar] [CrossRef]
  74. Liebig, M.A.; Franzluebbers, A.J.; Alvarez, C.; Chiesa, T.D.; Lewczuk, N.; Piñeiro, G.; Posse, G.; Yahdjian, L.; Grace, P.; Machado Rodrigues Cabral, O.; et al. MAGGnet: An international network to foster mitigation of agricultural greenhouse gases. Carbon Manag. 2016, 7, 243–248. [Google Scholar] [CrossRef]
  75. Valkama, E.; Kunypiyaeva, G.; Zhapayev, R.; Karabayev, M.; Zhusupbekov, E.; Perego, A.; Schillaci, C.; Sacco, D.; Moretti, B.; Grignani, C.; et al. Can conservation agriculture increase soil carbon sequestration? A modelling approach. Geoderma 2020, 369, 114298. [Google Scholar] [CrossRef]
  76. Weber, C.; Navarrete, F.; Perona, L.; Acciaresi, H. Remote sensing of nitrogen status in wheat by a radiometric response of its canopy. J. Plant Nutr. 2017, 40, 1877–1886. [Google Scholar] [CrossRef]
  77. Sultana, S.R.; Ali, A.; Ahmad, A.; Mubeen, M.; Zia-Ul-Haq, M.; Ahmad, S.; Ercisli, A.; Jaafar, H.Z.E. Normalized Difference Vegetation Index as a Tool for Wheat Yield Estimation: A Case Study from Faisalabad, Pakistan. Sci. World J. 2014, 2014, 8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  78. Ballester, C.; Hornbuckle, J.; Brinkhoff, J.; Smith, J.; Quayle, W. Assessment of In-Season Cotton Nitrogen Status and Lint Yield Prediction from Unmanned Aerial System Imagery. Remote Sens. 2017, 9, 1149. [Google Scholar] [CrossRef] [Green Version]
  79. Lu, S.; Lu, F.; You, W.; Wang, Z.; Liu, Y.; Omasa, K. A robust vegetation index for remotely assessing chlorophyll content of dorsiventral leaves across several species in different seasons. Plant Methods 2018, 14, 15. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  80. Asner, G.P. Biophysical and Biochemical Sources of Variability in Canopy Reflectance. Remote Sens. Environ. 1998, 64, 234–253. [Google Scholar] [CrossRef]
  81. Ribera, L.A.; Hons, F.M.; Richardson, J.W. An Economic Comparison between Conventional and No-Tillage Farming Systems in Burleson County, Texas. Agron. J. 2004, 96, 415–424. [Google Scholar] [CrossRef]
Figure 1. Experimental location (on the left), planimetry and relative georeferenced positions of sampling biomass points during the two-year experimental survey.
Figure 1. Experimental location (on the left), planimetry and relative georeferenced positions of sampling biomass points during the two-year experimental survey.
Agriculture 11 00039 g001
Figure 2. NDRE Vegetation index maps during the 2019 year at stem elongation—Zadoks scale: 35 (on the (left)) and at anthesis—Zadoks scale: 60 (on the (right)).
Figure 2. NDRE Vegetation index maps during the 2019 year at stem elongation—Zadoks scale: 35 (on the (left)) and at anthesis—Zadoks scale: 60 (on the (right)).
Agriculture 11 00039 g002
Table 1. Mean air temperature (°C) and monthly precipitations (mm) observed in the experimental period (years 2018–2019).
Table 1. Mean air temperature (°C) and monthly precipitations (mm) observed in the experimental period (years 2018–2019).
MonthsJanuaryFebruaryMarchAprilMayJuneJulyAugustSeptemberOctoberNovemberDecember
Mean air T (°C) Average
20188.64.88.916.119.122.725.825.721.416.511.96.815.7
20195.38.512.113.515.425.026.026.120.917.513.49.016.1
Precipitation (mm) Total
201829173143379548572519674261796
201970223659165110515112626653765
Table 2. Soil properties of the 0–20 cm layer in the conventional tillage (CT) and no tillage (NT) in 2020.
Table 2. Soil properties of the 0–20 cm layer in the conventional tillage (CT) and no tillage (NT) in 2020.
SM 1SandSiltClaypHC/N 2Organic MatterTotal Nitrogen
(g kg−1)(g kg−1)(g kg−1) (g kg−1)(g kg−1)
CT 3423 (±44) a126 (±20) a451 (±58) a7.8 (±0.2) a7.9 (±0.5) b13.20 (±1.07) b0.98 (±0.04) b
NT 4436 (±42) a121 (±20) a443 (±39) a7.8 (±0.2) a8.6 (±0.5) a21.52 (±6.40) a1.44 (±0.36) a
1 SM: Soil Management; 2 C/N: Ratio Carbon Nitrogen; 3 CT: Conventional Tillage; 4 NT: No-Tillage. Means within column that are followed by the same letter are not significantly different at p < 0.05%.
Table 3. Agronomic management practices adopted during the two-year experimental period (dd/mm/yyyy).
Table 3. Agronomic management practices adopted during the two-year experimental period (dd/mm/yyyy).
Agro-TechniqueSM 12017–20182018–2019
Ploughing (40 cm)CT 202/10/201726/09/2018
Weed control: Glyphosate 3NT 430/10/201728/10/2018
Harrowing and seed bed preparationCT20/11/201701/11/2018
Sowing 5NT and CT21/11/201730/11/2018
Weed control: Pinoxaden 6NT and CT28/03/201808/03/2019
N fertilization 7NT and CT29/03/201818/03/2019
NT and CT30/04/201802/05/2019
Pest control: Azoxystrobin, Cyproconazole 8NT and CT24/04/201822/04/2019
HarvestNT and CT06/07/201707/07/2019
1 SM: Soil Management; 2 CT: conventional tillage; 3 dose: 2.25 kg ha−1 of active ingredient; 4 NT: no-tillage; 5 Seed rate: 220 kg ha−1; row spacing: 0.17 m; 6 dose: 30 g ha−1 of active ingredient; 7 N source: urea (46%); 50% of N distribution for each date; 8 dose: 0.16 l ha−1 of active ingredient.
Table 4. Crop phenological stage (Zadoks scale) and remote sensing images acquisition details.
Table 4. Crop phenological stage (Zadoks scale) and remote sensing images acquisition details.
Phenological StageDate (mm/dd/yyyy)SensorForward-Side OverlapAbove Ground Level (m)Flight Speed (ms−1)Mode of FlightGround Resolution
(mm pixel−1)
ZS22 103/20/2018MAIA S-275%–65%505Automatic23
ZS3504/28/2018MAIA S-275%–65%505Automatic23
ZS6005/27/2018MAIA S-275%–65%505Automatic23
ZS2203/17/2019MAIA S-275%–65%505Automatic23
ZS3504/18/2019MAIA S-275%–65%505Automatic23
ZS6005/24/2019MAIA S-275%–65%505Automatic23
1 ZS: Zadoks scale.
Table 5. Vegetation index categories, calculated vegetation indices and related formula.
Table 5. Vegetation index categories, calculated vegetation indices and related formula.
Vegetation IndexFormulaReferences
ARVI 1 ARVI = NIR RB NIR + RB
Where: RB   =   Red y   ·   Blu Red
[52]
GARI 2 GARI = NIR Green y   ·   Blue Red NIR + Green y   ·   Blue Red
Where: y = 1 , 7
[53]
MSAVI2 3 MSAVI 2 = 2 · NIR + 1 2   ·   NIR + 1 2 8 · NIR Red 2 [54]
NDRE 4 NDRE   =   NIR     Red   Edge NIR   +   Red   Edge [55]
NDVI 5 NDVI = NIR Red NIR + Red [56]
GLI 6 GLI = Green Red + Green Blue 2 · Green + Red + Blue [57]
VARI 7 VARI = Green Red Green + Red Blue [53]
VARI Green 8 VARI   Green = Green Red Green + Red Violet [53]
VDVI 9 VDVI = 2 · Green Red Blue 2 · Green + Red + Blue [57]
Vig 10 Vig = Green Red Green + Red [53]
1 ARVI: Atmospherically resistant vegetation index; 2 GARI: Green atmospherically resistant vegetation index; 3 MSAVI2: Modified soil-adjusted vegetation index 2; 4 NDRE: Normalized difference red edge index; 5 NDVI: Normalized difference vegetation index; 6 GLI: Green leaf index; 7 VARI: Visible atmospherically resistant index; 8 VARI Green: Visible atmospherically resistant index green; 9 VDVI: Visible-band difference vegetation index; 10 Vig: Vegetation index green.
Table 6. Results of the analysis of variance (ANOVA) applied to the mixed model for measured crop parameters.
Table 6. Results of the analysis of variance (ANOVA) applied to the mixed model for measured crop parameters.
Factorsdf 1Chlorophyll (mg g−1)%NNNI 2SPAD
p-Valuep-Valuep-Valuep-Value
SM 310.160.16****
N 42************
Y 510.840.840.150.20
SM x N2*****
SM x Y10.660.660.720.96
N x Y20.790.790.230.37
SM x N x Y20.620.610.500.32
1 df: Degree of freedom; 2 NNI: Nitrogen Nutrition Index; 3 SM: Soil management; 4 N: Nitrogen input; 5 Y: Year; *: Significant at p < 0.05%; **: Significant at p < 0.01%; ***: Significant at p < 0.001%.
Table 7. Estimated marginal means for crop parameters.
Table 7. Estimated marginal means for crop parameters.
SM 1N 2Chlorophyll (mg g−1)%N 3NNI 4SPAD
CT 5N0 61.51 (±0.30) d1.44 (±0.58) d0.42 (±0.19) d25.35 (±6.56) d
NT 7N01.95 (±0.53) c1.73 (±0.49) c0.58 (±0.26) c34.32 (±4.21) c
CTN1 82.34 (±0.71) b1.91 (±0.57) b0.87 (±0.42) b43.10 (±8.14) b
NTN12.28 (±0.76) b1.97 (±0.53) b0.94 (±0.52) b42.79 (±8.57) b
CTN2 92.73 (±0.93) a2.16 (±0.83) a1.18 (±0.63) a47.43 (±9.66) a
NTN22.72 (±0.91) a2.21 (±0.66) a1.20 (±0.65) a48.44 (±9.67) a
1 SM: Soil management; 2 N: Nitrogen input; 3 %N: Nitrogen % content, 4 NNI: Nitrogen Nutrition Index; 5 CT: Conventional Tillage; 6 N0 = 0 kg N ha−1; 7 NT: No-Tillage; 8 N1 = 90 kg N ha−1; 9 N2 = 180 kg N ha−1. Means within column that are followed by the same letter are not significantly different at p < 0.05%.
Table 8. Results of the ANOVA applied to the mixed model for evaluated vegetation indices.
Table 8. Results of the ANOVA applied to the mixed model for evaluated vegetation indices.
Factorsdf 2NIR 1-Based Vegetation IndexNot NIR-Based Vegetation Index
ARVIMSAVI2GARINDRENDVIGLIVDVIVARIVARI GreenVig
SM 31*0.15*****0.530.530.220.870.88
N 42******************************
Y 510.370.820.620.440.410.240.240.280.110.22
SM x N2******0.520.520.570.530.48
SM x Y10.490.680.520.930.870.860.860.440.980.9
N x Y20.370.420.440.080.50.590.590.690.070.82
SM x N x Y20.870.740.540.760.610.760.760.670.690.72
1 NIR: Near infrared; 2 df: Degree of freedom; 3 SM: Soil management; 4 N: Nitrogen input; 5 Y: Year; *: Significant at p < 0.05%; **: Significant at p < 0.01%; ***: Significant at p < 0.001%.
Table 9. Differences between near-infrared (NIR)-based vegetation index on the compared fertilization levels.
Table 9. Differences between near-infrared (NIR)-based vegetation index on the compared fertilization levels.
SM 1N 2ARVIGARIMSAVI2NDRENDVI
CT 3N0 40.28 (±0.20) d0.17 (±0.16) d0.63 (±0.22) d0.16 (±0.07) d0.40 (±0.17) d
NT 5N00.42 (±0.26) c0.32 (±0.21) c0.74 (±0.31) c0.22 (±0.10) c0.52 (±0.21) c
CTN1 60.55 (±0.28) b0.44 (±0.26) b0.83 (±0.28) b0.30 (±0.16) b0.62 (±0.22) b
NTN10.56 (±0.31) b0.48 (±0.29) b0.81 (±0.30) b0.34 (±0.17) b0.64 (±0.25) b
CTN2 70.63 (±0.28) a0.55 (±0.26) a0.89 (±0.27) a0.39 (±0.18) a0.72 (±0.20) a
NTN20.63 (±0.32) a0.56 (±0.28) a0.88 (±0.31) a0.40 (±0.19) a0.71 (±0.23) a
1 SM: Soil management; 2 N: Nitrogen input; 3 CT: Conventional Tillage; 4 N0 = 0 kg N ha−1; 5 NT: No-Tillage; 6 N1 = 90 kg N ha−1; 7 N2 = 180 kg N ha−1. Means within column that are followed by the same letter are not significantly different at p < 0.05%.
Table 10. Differences between the not NIR-based vegetation index on the compared fertilization levels.
Table 10. Differences between the not NIR-based vegetation index on the compared fertilization levels.
Nitrogen InputGLIVARIVDVIVARI GreenVig
N0 10.19 (±0.10) b0.14 (±0.18) c 0.19 (±0.10) b0.12 (±0.18) c0.10 (±0.13) c
N1 20.27 (±0.11) a0.27 (±0.23) b0.27 (±0.11) a0.24 (±0.23) b0.20 (±0.16) b
N2 30.28 (±0.11) a0.32 (±0.25) a0.28 (±0.11) a0.28 (±0.24) a0.23 (±0.16) a
1 N0 = 0 kg N ha−1; 2 N1 = 90 kg N ha−1; 3 N2 = 180 kg N ha−1. Means within column that are followed by the same letter are not significantly different at p < 0.05%.
Table 11. Coefficient of determination and root mean square error (RMSE) of linear relationship between the leaves chlorophyll concentration (mg g−1), NNI and SPAD or vegetation indices.
Table 11. Coefficient of determination and root mean square error (RMSE) of linear relationship between the leaves chlorophyll concentration (mg g−1), NNI and SPAD or vegetation indices.
Chlorophyll (mg g−1)NNI
R2RMSER2RMSE
SPAD0.630.510.720.29
ARVI0.580.540.760.27
GARI0.610.520.740.28
GLI0.280.710.600.35
MSAVI20.680.470.400.43
NDRE0.610.520.850.21
NDVI0.570.550.740.28
VARI Green0.190.750.690.31
VARI0.360.670.710.30
VDVI0.280.710.600.35
Vig0.340.680.700.30
Table 12. The average of coefficient of determination R2, RMSE of linear relationship between the crop parameters and the no-destructive index, and difference between the averages of NIR band and no NIR band correlation parameters.
Table 12. The average of coefficient of determination R2, RMSE of linear relationship between the crop parameters and the no-destructive index, and difference between the averages of NIR band and no NIR band correlation parameters.
Chlorophyll (mg g−1)NNI
R2RMSER2RMSE
NIR band0.610.520.70.29
No NIR band0.290.70.660.32
Difference0.32−0.180.04−0.03
Table 13. Results of the ANOVA applied to the mixed model for measured yield components.
Table 13. Results of the ANOVA applied to the mixed model for measured yield components.
Factorsdf 1KS 2Spikes m−2Grain Yield (t ha−1)
SM 31***0.41*
N 42*********
SM x N2*0.42*
1 df: Degree of freedom; 2 KS: Kernel per spike; 3 SM: Soil management; 4 N: Nitrogen input; *: Significant at p < 0.05%; ***: Significant at p < 0.001%.
Table 14. Yield component differences between compared treatments.
Table 14. Yield component differences between compared treatments.
Soil ManagementNitrogen InputKS 1Spikes m−2Grain Yield (t ha−1)
CT 2N0 310 (±3) d380 (±85) c1.46 (±0.28) d
NT 4N017 (±4) c420 (±57) c2.51 (±0.22) c
CTN1 529 (±10) b522 (±53) b3.62 (±0.42) b
NTN131 (±9) b539 (±53) b4.08 (±0.52) b
CTN2 634 (±9) a574 (±97) a5.29 (±1.49) a
NTN236 (±5) a559 (±76) a5.26 (±1.25) a
1 KS: Kernel per spike; 2 CT: Conventional Tillage; 3 N0 = 0 kg N ha−1; 4 NT: No-Tillage; 5 N1 = 90 kg N ha−1; 6 N2 = 180 kg N ha−1. Means within column that are followed by the same letter are not significantly different at p < 0.05%.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Fiorentini, M.; Zenobi, S.; Orsini, R. Remote and Proximal Sensing Applications for Durum Wheat Nutritional Status Detection in Mediterranean Area. Agriculture 2021, 11, 39. https://doi.org/10.3390/agriculture11010039

AMA Style

Fiorentini M, Zenobi S, Orsini R. Remote and Proximal Sensing Applications for Durum Wheat Nutritional Status Detection in Mediterranean Area. Agriculture. 2021; 11(1):39. https://doi.org/10.3390/agriculture11010039

Chicago/Turabian Style

Fiorentini, Marco, Stefano Zenobi, and Roberto Orsini. 2021. "Remote and Proximal Sensing Applications for Durum Wheat Nutritional Status Detection in Mediterranean Area" Agriculture 11, no. 1: 39. https://doi.org/10.3390/agriculture11010039

APA Style

Fiorentini, M., Zenobi, S., & Orsini, R. (2021). Remote and Proximal Sensing Applications for Durum Wheat Nutritional Status Detection in Mediterranean Area. Agriculture, 11(1), 39. https://doi.org/10.3390/agriculture11010039

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop