Next Article in Journal
Effect of Polyethylene Glycol-Simulated Drought Stress on Stomatal Opening in “Modern” and “Ancient” Wheat Varieties
Next Article in Special Issue
The Long-Term Effects of Barren Land Afforestation on Plant Productivity, Soil Fertility, and Soil Moisture in China: A Meta-Analysis
Previous Article in Journal
Effects of Long-Term Sod Culture Management on Soil Fertility, Enzyme Activities, Soil Microorganisms, and Fruit Yield and Quality in “Jiro” Sweet Persimmon Orchard
Previous Article in Special Issue
Effects of Straw Returning and New Fertilizer Substitution on Rice Growth, Yield, and Soil Properties in the Chaohu Lake Region of China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Integrating NDVI, SPAD, and Canopy Temperature for Strategic Nitrogen and Seeding Rate Management to Enhance Yield, Quality, and Sustainability in Wheat Cultivation

1
Research Unit of Plant Breeding and Genetic Resources Conservation, Regional Center of Agricultural Research of Settat, National Institute of Agricultural Research, P.O. Box 589, Settat 26000, Morocco
2
Laboratory of Agrifood and Health, Faculty of Sciences and Techniques, Hassan First University of Settat, P.O. Box 577, Settat 26000, Morocco
3
Dryland Research Center, National Institute of Agricultural Research, P.O. Box 589, Settat 26000, Morocco
*
Authors to whom correspondence should be addressed.
Plants 2024, 13(11), 1574; https://doi.org/10.3390/plants13111574
Submission received: 10 May 2024 / Revised: 29 May 2024 / Accepted: 2 June 2024 / Published: 6 June 2024
(This article belongs to the Special Issue Soil Fertility, Plant Nutrition and Nutrient Management)

Abstract

:
This study explores the interplay between nitrogen doses and seeding rates on wheat yield, biomass, and protein content. Utilizing tools such as the Normalized Difference Vegetation Index (NDVI), Soil Plant Analysis Development (SPAD) measurements, and canopy temperature (CT), we conducted experiments over five growing seasons. The treatments included three nitrogen levels (0, 60, 120 kg/ha) and three seeding rates (300, 400, 500 seeds/m2) in a split-plot design with 90 plots and two replications. Our results show that an intermediate nitrogen dose (60 kg/ha) combined with a moderate seed rate (400 seeds/m2) enhances wheat yield by 22.95%. Reduced nitrogen levels increased protein content, demonstrating wheat’s adaptive mechanisms under nitrogen constraints. NDVI analysis highlighted significant growth during the tillering phase with high nitrogen, emphasizing early-stage nutrient management. SPAD measurements showed that early nitrogen applications boost chlorophyll content, essential for vigorous early growth, while CT data indicate that optimal nitrogen and seed rates can effectively modulate plant stress responses. As crops mature, the predictive capacity of NDVI declines, indicating the need for adjusted nitrogen strategies. Collectively, these findings advocate for refined management of nitrogen and seeding rates, integrating NDVI, SPAD, and CT assessments to enhance yields and promote sustainable agricultural practices while minimizing environmental impacts.

1. Introduction

In wheat farming, providing a detailed understanding of crop management practices, particularly nitrogen use and seeding strategies, is critical to advancing sustainable agricultural practices. These factors are integral to enhancing crop yields, plant health, and nutritional quality, profoundly affecting both the productivity and ecological sustainability of farming operations. Developments in precision agriculture technologies such as the Normalized Difference Vegetation Index (NDVI), Soil Plant Analysis Development (SPAD) measurements, and canopy temperature assessments have dramatically improved our ability to refine these vital agricultural inputs [1,2,3].
Durum wheat (Triticum durum Desf.) is a staple crop globally, known for its adaptability and nutritional value. Effective cultivation of durum wheat requires careful management of nitrogen (N) and seeding rates to optimize yield and quality. Nitrogen management is particularly crucial because of its essential role in plant growth and the associated environmental risks, including water pollution and greenhouse gas emissions. Carefully managing nitrogen applications helps not only to improve plant health, but also to minimize ecological damage [4,5]. Additionally, selecting appropriate seeding rates is vital as it directly impacts how plants utilize the available space and resources, influencing their growth patterns, yield outcomes, and overall quality. Advanced strategies for optimizing nitrogen management can significantly enhance the yield and environmental performance of wheat crops while minimizing negative impacts, such as nitrate leaching and the emission of nitrous oxide [4,5]. Moreover, the integration of seeding rates with nitrogen management practices is essential for achieving optimal agronomic and environmental outcomes in wheat production [6,7].
Employing NDVI, we gain real-time insights into plant health, enabling precise adjustments to farming practices that cultivate optimal growth conditions. NDVI provides detailed feedback on plant light absorption and reflection, indicating plant health and the effectiveness of our agricultural strategies [8]. This technique has been widely used for various crops, including wheat, to estimate yield and monitor growth stages [9,10]. Studies have demonstrated NDVI’s effectiveness in optimizing nitrogen management and reducing environmental impacts [11]. Additionally, NDVI measurements have shown strong correlations with wheat yield variability, supporting its use in developing new fertilization strategies aligned with crop needs throughout the growing season [12]. SPAD measurements complement these insights by providing an estimate of chlorophyll content, which correlates with nitrogen levels in plants, shedding light on their photosynthetic activity [13]. These technologies, including multispectral remote sensing and machine learning, contribute significantly to precision agriculture by enabling more accurate monitoring of crop status and productivity. This allows for more targeted and efficient management of resources, ultimately leading to improved yield and sustainability of wheat production [12,14].
This research employed a split-plot experimental design over five growing seasons, involving five durum wheat cultivars and 90 plots to examine the effects of different seeding rates (300, 400, 500 seeds/m2) under varying nitrogen levels (0, 60, 120 kg/ha), with each plot replicated twice. This design mimicked real-world agricultural conditions, providing robust data on crop responses. Our study explores the relationship between nitrogen dosages and seeding rates to identify optimal strategies for managing wheat crops effectively. By integrating data from NDVI, SPAD, and canopy temperature measurements, we aim to refine farming practices that not only boost yields, but also do so sustainably. This approach leverages the capabilities of precision agriculture tools such as UAV multispectral imagery and machine learning algorithms to precisely manage nitrogen application, thus enhancing both productivity and environmental stewardship. Investigations into the effects of nitrogen management practices on SPAD values and NDVI readings have demonstrated the potential of these tools in optimizing wheat production under various environmental conditions [1,14,15,16].
The goal of this study was to investigate how different nitrogen levels and seed rates affect wheat yield, biomass, and protein content, utilizing precision agriculture tools. We aimed to identify and refine agricultural practices that enhance crop performance at various growth stages while supporting sustainable farming objectives. Through this research, we intend to offer actionable recommendations that can help improve wheat farming efficiency and environmental sustainability, contributing to better food security and agricultural practices.

2. Materials and Methods

2.1. Study Site and Environmental Conditions

The research covered five growing seasons (S): 2016–2017, 2017–2018, 2018–2019, 2019–2020, and 2020–2021, and took place at the Sidi El Aidi facility of the National Institute of Agricultural Research (INRA) in Morocco. Precipitation data over the five growing seasons are provided in Table 1. Further soil and site-specific parameters are summarized in Table 2. Soil samples were collected and analyzed annually before planting, with the results averaged to provide an overall view of soil conditions. The wheat was grown in a rotation with a fallow period each year to avoid additional nitrogen input from previous legume crops in the experimental station and to maintain soil health.

2.2. Experimental Design, Agronomic Practices, and Genetic Material

Each experimental plot covered an area of 2.7 square meters, measuring 2.5 m in length and 1.08 m in width. The study employed a split-plot design with five durum wheat cultivars and 90 plots to examine the effects of different seed rates (300, 400, 500 seeds/m2) under varying nitrogen levels (0, 60, 120 kg/ha), with each plot replicated twice. Standard agronomic practices, including soil preparation, weed control, and irrigation, were implemented to ensure optimal growth conditions. Sowing occurred in mid-November using a Wintersteiger plot seeder (Wintersteiger, Ried im Innkreis, Austria). Nitrogen was applied in two stages: half at the tillering stage and the remaining half at the stem extension stage. Ammonium nitrate fertilizer with a 33.5% nitrogen content was applied according to the treatment levels specified. Pest management included the application of Vitavax Ultra (Dhanuka, Haryana, India) for seed treatment and Roundup® TURBO (St. Louis, MI, USA) for pre-emergence weed control. Manual weeding was performed as needed to maintain optimal crop health and to eliminate other factors that can affect the main factor of the study.
Five durum wheat cultivars were used in the experiments: Karim, Nassira, Faraj, Luiza, and Itri. These cultivars were selected for their adaptability and performance under Moroccan conditions (Table 3). The data used in this study are the mean of varieties for each combination of nitrogen and seed rate to remove the varietal effect and the mean of all the growing seasons to remove the environmental effect and the year-to-year variation.

2.3. Methodology for Data Acquisition and Trait Analysis

Yield measurements were quantified by harvesting the entire grain from each plot, measured in grams over an area of 2.7 square meters, later standardized to kilograms per hectare to ensure uniformity across the study. Biomass assessments occurred shortly before harvest, with the entire above-ground biomass from each plot recorded in kilograms to reflect productivity. For the evaluation of the critical quality trait of protein content, we utilized the capabilities of Chopin Technologies’ Infraneo, a near-infrared spectroscopy (NIRS) instrument (Chopin Technologies, Villeneuve-la-Garenne, France). This device underwent routine calibration at the National Institute of Agricultural Research facility in Rabat to maintain measurement precision. We strengthened the validity of the NIRS assessments through cross-validation processes, which involved comparing the Infraneo readings with those from a FOSS Infratec NIR analyzer (Infratec, Dresden, Germany), carefully calibrated and operated at INRA Settat. The protein content analysis was conducted using 800 g of grain collected from the harvested yield of each plot.
This dual-analytical approach not only enhanced the robustness of our protein content analysis, but also allowed us to cross-verify our NIRS data against the gold-standard Kjeldahl method. Through this validation process, we ensured that our protein content measurements were both accurate and reliable. This strategy of employing dual verification methods reinforces the integrity and consistency of our findings, providing a solid foundation for our research conclusions.
To augment our understanding of the plant physiological responses and enhance the precision of our agronomic evaluations, we incorporated several measurement techniques.
Normalized Difference Vegetation Index (NDVI) measurements were performed using a FieldSpec HandHeld FSHH 325-3075P Spectris multispectral radiometer (Artisan Technology Group, Champaign, IL, USA). This tool measures reflectance at specific wavelengths to provide detailed insights into plant health and vigor. NDVI was calculated using the formula N D V I = R 900 R 680 R 900 + R 680 , where R900 is the reflectance at the near-infrared wavelength of 900 nm, and R680 is the reflectance at the visible light wavelength of 680 nm. Measurements were conducted during key phenological stages: tillering (Zadoks 20–26), stem extension (Zadoks 30–39), and heading (Zadoks 50–58) [17]. Each measurement was repeated ten times within each plot to account for variability, and the mean values were used for analysis. These measurements were taken under clear sky conditions between 10:00 a.m. and 2:00 p.m. to ensure consistency and minimize the impact of fluctuating light conditions Figure 1a.
Chlorophyll content, an indicator of plant health and nitrogen status, was measured using a SPAD-502 Plus meter (Konika Minolta, Tokyo, Japan). This non-destructive method allowed for rapid in-field assessment of chlorophyll levels providing a direct measure of the photosynthetic potential and nitrogen content of the plants. Measurements were conducted during the same key phenological stages as NDVI (tillering, stem extension, and heading). During each stage, measurements were taken from ten randomly selected plants per plot. From each plant, chlorophyll content was measured on the penultimate fully expanded leaf Figure 1b. Each leaf was measured three times, and the average value was recorded to ensure accuracy and consistency.
Canopy temperature was measured using an infrared thermometer to assess crop water stress and physiological responses. Measurements were taken under clear sky conditions between 10:00 a.m. and 2:00 p.m., providing data critical for understanding the thermal dynamics of the crop canopy and its correlation with water use efficiency and stress tolerance. For each plot, two measurements were taken from the part most exposed to the sun, ensuring the operator’s shadow and neighboring plots’ shadows were avoided. Measurements were consistently recorded from the same end of each plot with the sun behind the operator. The trigger was held for 3–5 s to allow the infrared thermometer to average the temperature readings (Figure 1c). Water stress in plants was assessed by comparing canopy temperatures between plots. Higher canopy temperatures indicate that the plants are experiencing more water stress due to decreased transpiration rates, leading to reduced cooling effects on the leaf surface. This method detects subtle temperature changes indicating stress levels based on plant physiology principles.
These methodologies, integrated into our study, enhance the robustness and precision of our agronomic assessments, enabling detailed analysis of how management practices influence crop performance.

2.4. Statistical Analysis Techniques

Data management began with the use of Microsoft Excel (version 2108, Build 14332.20324) for collecting all data annually for each trait studied and performing preliminary calculations, such as computing means, to prepare the dataset for further analysis.
For inferential statistics, we utilized Minitab 18 to conduct a two-way analysis of variance (ANOVA), which helped assess the significance of our findings and calculate standard deviations and F-values. This facilitated a deeper understanding of data variability and effect significance.
Following ANOVA, Tukey’s Honest Significant Difference (HSD) method was applied for post hoc tests to control the family-wise error rate and verify the distinctiveness of group means, ensuring rigorous comparison standards. The tests were performed at a significance level of 0.05. Minitab 18 was used to conduct the analysis, evaluating the effects of varying nitrogen doses and seeding rates. The results indicated mean values for each combination, accompanied by their respective letter annotations indicating statistical significance.
In addition, we employed R software (version 4.4.0) to perform Pearson’s correlation analyses, creating matrices that helped elucidate linear relationships between variables.
To visually display our results, we used R to generate heatmaps, which visually represented the data’s variability through color gradients. Additionally, R was used to generate trend figures showing the NDVI and SPAD values across key growth stages under nitrogen doses. Origin Pro (version 2024b, 10.1.5.132) was used to produce boxplots that graphically depicted data distribution and central tendencies, offering insights into data spread. For the boxplot figures, significance levels were determined by one-way ANOVA and the subsequent Tukey’s test using Minitab 18.

3. Results

3.1. Analysis of Yield Responses to Nitrogen and Seed Rates

The agronomic experiment conducted provided an evaluation of wheat yield influenced by defined nitrogen doses (N1: 120 kg/ha, N2: 60 kg/ha, N3: 0 kg/ha) and seed rates (S1: 500 seeds/m2, S2: 400 seeds/m2, S3: 300 seeds/m2). Interactions between these variables on yield are graphically represented in Figure 2’s heatmap and quantitatively delineated in Table 4.
ANOVA results, presented in Table 4, highlight that nitrogen doses significantly influenced yield (F = 4.65, p = 0.012), while the seed rate alone did not have a significant impact (F = 0.92, p = 0.402). Moreover, the interaction between nitrogen dose and seed rate was significant (F = 4.38, p = 0.003), illustrating the crucial role of these combined factors in determining yield. As visually illustrated in the heatmap of Figure 2 and supported by the boxplot in Figure 3, this relationship was explored in more depth. The peak mean yield was achieved with the N2 dose (60 kg/ha) at a moderate nitrogen application rate, combined with the S2 seed rate of 400 seeds/m2, resulting in an optimal yield of 2978 kg/ha. This yield significantly surpassed the mean yield for the N2 dose by 22.95%, highlighting the effectiveness of this treatment combination. Conversely, the yield substantially decreased at the S3 seed rate with the same nitrogen dose (N2), yielding a mean of 2008 kg/ha—indicative of a 17.09% decline from the N2 average and a significant 32.58% fall from the highest observed yield.
The delineation of an optimal yield zone within the intermediate nitrogen dose range and a specific seed rate underscores the critical balance required in agronomic management. This balance is essential in directing the optimization of fertilization and seeding strategies for attaining maximum yield. These strategic insights have practical implications for promoting sustainable agricultural practices.

3.2. Analysis of Biomass Responses to Nitrogen and Seed Rates

The agronomic investigation conducted evaluated the dynamics of wheat yield and biomass, offering a global view of plant productivity. Both Figure 4 and Table 4 demonstrate that the agronomic treatments optimized yield and significantly increased biomass. The treatment conditions that produced the highest mean biomass further validated the effectiveness of the agronomic practices utilized. Specifically, the N2 dose paired with the S2 seed rate resulted in the highest biomass measurement at 9742 kg/ha, a 12.08% increase over the average biomass for the N2 treatments. This uniformity across different plant growth metrics implies the potential for a harmonized agronomic management approach. The results presented in Figure 4 and Figure 5 and Table 4 illustrate the impact of these agronomic strategies on productivity.
Figure 6 offers additional clarity on the influence of seed rates on biomass within the context of the N2 nitrogen dose. Here, the S2 seed rate under the N2 dose demonstrates a considerably higher mean biomass than S1 and S3, paralleling the yield trend. The mean biomass at the S2 seed rate was 9742 kg/ha, clearly indicating this seed density’s efficacy in maximizing biomass, akin to its impact on yield. However, the influence of seed rate on biomass is not significant when considered independently, as evidenced by an F-value of (F = 0.62, p = 0.539) It is within the interaction with the N2 nitrogen dose that the seed rate’s effect becomes significant, underscoring the importance of considering these factors in combination rather than isolation (F = 3.53, p = 0.01).
The findings portrayed in the boxplots are robustly supported by statistical evidence from Table 4, particularly the ANOVA results, which affirm the significant interplay between nitrogen doses and seed rates on biomass. While the impact of seed rate alone on biomass did not reach statistical significance, the interaction effects are substantial, reflecting the dynamics at play in crop production systems.

3.3. Analysis of Protein Content Responses to Nitrogen and Seed Rates

The analysis expanded from wheat yield and biomass to protein content to explore the variability of this essential quality parameter under varying nitrogen doses and seeding rates. Figure 7’s heatmap reveals contrasts in protein content across nitrogen doses, with the most pronounced protein content identified in the N3 group at an average of 18.25%. The ANOVA did not detect a significant interaction effect between nitrogen dose and seed rate on protein content, indicating that the observed differences are primarily due to the nitrogen dose alone.
Further investigating the protein distribution, Figure 8 demonstrates the effect of nitrogen doses, with the N3 treatment exhibiting a higher mean protein percentage than the N1 and N2 treatments. This trend suggests an adaptive response of wheat to limited nitrogen availability, which may represent a compensatory mechanism in protein concentration. Conversely, Figure 9 examines seed rate impacts and indicates that no particular seed rate significantly influences protein content.
The protein content analysis is statistically reinforced by Table 4 of the ANOVA results, which confirm a significant effect of nitrogen dose on protein content (F = 4.11, p = 0.02). While the seed rate alone did not show a significant impact (F = 1.79, p = 0.173), and neither did the interaction between nitrogen and seed rate (F = 1.88, p = 0.122), these outcomes hint at a latent complexity in how these factors influence protein content, indicating potential areas for future research.
The dataset indicates that optimizing protein content in wheat may align with reduced nitrogen input. This scenario contrasts with the conventional agronomic strategy of increasing nitrogen to boost yield. The data suggest that lower nitrogen doses result in higher protein content, challenging existing paradigms.

3.4. Evaluating the Impact of Nitrogen and Seed Rates on NDVI through Key Growth Stages of Wheat

The comprehensive investigation into wheat cultivation reported here expands the scope of traditional agronomic research by employing the Normalized Difference Vegetation Index (NDVI) to gauge plant health and photosynthetic capability at critical growth stages, thus offering essential insights for precision agriculture enhancements.
Statistical analyses from Table 5 confirm significant differences in nitrogen doses at both the tillering (F = 9.64, p < 0.001) and stem extension stages (F = 9.42, p < 0.001). However, the heading stage shows non-significant results (F = 3.09, p = 0.051), underlining a detailed relationship between plant development and nitrogen supply that demands additional investigation.
Notably, seed rate alone did not significantly impact NDVI, suggesting that it may have a less dominant effect compared to nitrogen. However, the interaction effects observed at the tillering stage (F = 2.56, p = 0.045) point to a dynamic interplay between nitrogen and seeding rate, with the potential for diminishing influence as the crop matures.
During the tillering phase, represented by NDVI1 in Figure 10, plants receiving the highest nitrogen dose (N1) of 120 kg/ha showed substantial growth across all seeding rates, particularly at the S3 seed rate (300 seeds/m2), underscoring the significant role of nitrogen in early plant development. This is further illustrated in Figure 11, where the distribution of NDVI1 values clarifies how nitrogen dose impacts plant health. Plants under the N1 treatment exhibited the highest NDVI1 values, indicating robust growth, whereas lower nitrogen doses showed decreased NDVI1 values, reflecting lesser vegetative vigor.
As the study transitioned into the stem extension phase, reflected in NDVI2 with Figure 12 and Figure 13, the highest nitrogen dose group (N1) at 120 kg/ha maintained its position with the most consistent plant health, as evidenced by an average NDVI2 value of 0.827. This value suggests a continuation of robust vegetative growth from the tillering phase, indicating the sustaining effects of the highest nitrogen application during this crucial growth stage. Contrary to expectations of a reduction in vigor, this higher nitrogen availability underpins a stable vegetative health as the plants prepare for reproductive development. The data compel a reassessment of nitrogen management strategies to ensure plant health is optimized throughout all stages of growth, acknowledging the crucial role of nitrogen, especially in the early to mid-phases of the crop cycle.
The data from Table 5, corroborated by Figure 14’s heatmap of mean NDVI3 values by nitrogen doses and seed rates, reveals a declining trend in NDVI3 across all nitrogen treatments as the wheat progresses to the heading stage. The highest nitrogen treatment (N1) group shows a mean NDVI3 value of 0.801, indicating a decrease from the tillering and stem extension phases. Similarly, the intermediate (N2) group exhibits a mean NDVI3 value of 0.784, and the non-fertilized (N3) group has the lowest mean value of 0.771, both confirming a general reduction in NDVI as the crops approach reproductive maturity. This trend reflects the anticipated shift in plant physiology from leafy, vegetative growth to grain filling, and suggests that the influence of nitrogen on vegetative vigor becomes less pronounced as wheat enters the reproductive phase. Table 5’s statistics provide a quantitative basis for this observation, underscoring the importance of adjusting nitrogen management strategies to align with the changing nutritional requirements of the wheat throughout its growth cycle.
In this study of wheat physiology, Figure 15 presents an insightful analysis of the interrelations among plant traits, particularly the Normalized Difference Vegetation Index (NDVI) during tillering (NDVI1), stem extension (NDVI2), and heading stages (NDVI3), and their influence on vital agronomic outputs such as yield, biomass, and protein content.
The Pearson’s correlation coefficients showcased in Figure 15 reveal that NDVI1 has a strong positive correlation with biomass (r = 0.65, p < 0.001) and a moderate positive correlation with yield (r = 0.38, p < 0.001). These correlations confirm that a strong start in plant growth is indicative of higher biomass accumulation and a substantial influence on yield. Additionally, the positive relationship between NDVI1 and NDVI2 (r = 0.64, p < 0.001) suggests that early growth vigor is sustained into the subsequent stem extension phase, emphasizing the importance of initial nitrogen management.
In contrast, as plants progress to the heading stage, NDVI3 displays only a weak association with biomass (r = 0.20, p < 0.05) and an even weaker relationship with yield (r = 0.11, p < 0.05), indicating that NDVI’s predictive capacity diminishes as the wheat matures and transitions from vegetative to reproductive growth phases.
Protein percentage demonstrates a lack of significant correlation with NDVI at all stages, suggesting that protein synthesis in wheat may be regulated by a set of factors distinct from those impacting vegetative growth. This finding is in line with the ANOVA results that show no significant effect of seeding rate or its interaction with nitrogen dose on protein content, revealing complex interactions that dictate this quality attribute.
The NDVI data derived from our study reinforce the critical need for an approach for nitrogen management throughout the wheat growth cycle. Our findings support a model of integrated management that optimizes both crop health and yield, promoting sustainable agricultural practices. This model emphasizes strategic nitrogen timing to enhance productivity efficiently while safeguarding environmental health.
The early-stage predictive capabilities of NDVI prove indispensable for agronomists, offering a method for the precise adjustment of agronomic inputs. These early indicators of plant health allow for the proactive management of resources, maximizing growth potential and ensuring resource use efficiency. This application of NDVI data not only improves yield outcomes, but also reduces the environmental footprint of agricultural practices by targeting input use where and when it is most effective.
Our research lays the groundwork for the next generation of farming practices that combine real-time monitoring with established agricultural methods to elevate crop yield. Looking ahead, leveraging NDVI within precision agriculture frameworks is set to transform crop management practices, ensuring they meet the dual objectives of enhanced productivity and environmental sustainability in an era of increasing global food demand.

3.5. Evaluating the Impact of Nitrogen and Seed Rates on SPAD Chlorophyll Measurements across Key Growth Stages of Wheat

Building upon our NDVI findings, we extend our agronomic analysis to the Soil Plant Analysis Development (SPAD) chlorophyll measurements, providing another dimension to our understanding of wheat physiology across developmental stages. The SPAD indices, like NDVI, serve as non-destructive proxies for assessing plant health, particularly chlorophyll content, which is closely related to nitrogen status and photosynthetic capacity.
Table 6’s statistical analysis reveals significant findings related to the influence of nitrogen dose and the interaction between nitrogen dose and seed rate on SPAD measurements, which are indicative of chlorophyll content in wheat leaves, a proxy for plant health and vigor.
For SPAD1, which reflects the chlorophyll content at the tillering stage, the nitrogen dose presents a statistically significant impact (F = 3.22, p = 0.045). This suggests that nitrogen application at this early growth stage is crucial for establishing the chlorophyll levels that support photosynthetic activity and initial plant development. The highest nitrogen treatment (N1) group exhibits the highest mean SPAD1 value, supporting the notion that adequate nitrogen supply is beneficial early in the plant’s life cycle and could guide targeted management practices to enhance early growth.
However, the significance of the nitrogen effect does not persist into the SPAD3 measurement stage (F = 1.65, p = 0.198), indicating that as the wheat plant progresses towards reproductive maturity, the direct influence of nitrogen on chlorophyll content becomes less pronounced. This observation suggests the need for a reassessment of nitrogen application strategies to ensure optimal resource utilization and support sustainable outcomes such as reduced nitrogen runoff and enhanced soil health.
Interestingly, the interaction between nitrogen dose and seed rate becomes statistically significant in the SPAD3 stage (F = 4.43, p = 0.003). This significant interaction effect at the heading stage implies that the combined influence of the amount of nitrogen available and the density of the seeding rate plays a more complex role in determining chlorophyll content during this later stage of development. It suggests that optimizing chlorophyll content, indicative of plant health during the heading stage, requires a balance between these factors, aligning agronomic practices with the plant’s changing nutrient demands.
During the tillering stage, SPAD1 measurements indicate that plants receiving the highest nitrogen dose (N1) exhibit a pronounced chlorophyll content, reflected in an average SPAD1 value of 51.49. The heatmap in Figure 16 emphasizes that within the N1 group, it is the S3 seed rate (300 seeds/m2) that displays the highest mean SPAD1 reading of 52.47. This is consistent with the boxplot shown in Figure 17, which further confirms the superior performance of the N1 nitrogen dose within the medium seed rate category, highlighting the nuanced effects of nitrogen applications on chlorophyll levels during early growth stages. These findings align with NDVI1 observations, demonstrating the critical role of appropriate nitrogen provision in establishing vigorous plant growth, as posited by precision agriculture’s targeted nutrient management approach.
As illustrated by Figure 18, SPAD2 measurements at the stem extension stage indicate a consistent retention of chlorophyll content across nitrogen treatments. Specifically, the average SPAD2 value for the highest nitrogen dose (N1) slightly declines to 50.62, reflecting a continued healthy chlorophyll presence in the leaves. This persistence from the tillering to the stem extension phase underscores the critical role of early nitrogen application in maintaining chlorophyll density and, by implication, plant health throughout the vital early growth stages.
Advancing to the heading stage, Table 6 and Figure 19 indicates a decline in SPAD3 values across all nitrogen doses, mirroring the pattern observed with NDVI3. The highest nitrogen dose group (N1) shows a decrease to an average SPAD3 value of 47.06. This downward trend across the nitrogen treatments highlights a general reduction in chlorophyll content as the crop matures, transitioning resources from vegetative growth to grain filling.
To better visualize the trends of NDVI and SPAD values throughout the key growth stages, Figure 20 presents these indices over time. The NDVI values, as shown in Figure 20a, display a consistent decline from tillering to heading. The highest nitrogen dose (N1) maintains superior values throughout all stages, with the mean NDVI values for N1 starting at 0.830 during tillering and decreasing to 0.801 by heading. This trend emphasizes the importance of early-stage nutrient management, with statistically significant differences observed among nitrogen treatments (F = 9.64, p = 0.000 for NDVI 1; F = 9.42, p = 0.000 for NDVI 2; Table 5).
Similarly, the SPAD values, illustrated in Figure 20b, follow a declining trend but highlight significant differences among the nitrogen treatments, particularly at the tillering stage. The mean SPAD values for N1 start at 51.49 during tillering and decrease to 47.06 by heading, showing significant variation among treatments (F = 3.22, p = 0.045 for SPAD 1; Table 6).
In synthesizing these SPAD results with our earlier findings on NDVI, yield, biomass, and protein content, it becomes clear that nitrogen application must be managed strategically across the wheat’s lifecycle to promote not just yield, but also physiological health and efficiency. This approach advocates for precision in agronomic decision-making, emphasizing sustainable practices that optimize the various facets of crop development for high-quality agricultural production.

3.6. Impact of Nitrogen and Seed Rates on Canopy Temperature Dynamics in Wheat

Building on the insights from our NDVI and SPAD metrics analysis, we further explored the physiological responses of wheat to varying agronomic treatments by examining canopy temperature (CT) dynamics.
The relationship between canopy temperature (CT), seed rate, and nitrogen dose during the tillering stage reveals fascinating insights into wheat’s physiological response to varying agronomic practices. Table 7 highlights a significant interaction effect between nitrogen dose and seed rate on CT1 (F = 8.26, p < 0.001), elucidated by the corresponding heatmap in Figure 21. This interaction points to the balance between nitrogen availability and seed density that modulates canopy temperature, a proxy for water stress and photosynthetic efficiency.
The CT1 results demonstrate the effects of nitrogen and seeding rates on canopy temperature during wheat’s tillering stage. Specifically, the highest seed rate (S1) combined with the highest nitrogen dose (N1) registers a CT1 mean temperature of 20.95 °C, which aligns with expectations of a warmer canopy due to denser planting. On the other hand, a significant observation is that the same high nitrogen dose (N1) paired with the lowest seed rate (S3) results in the coolest mean CT1 temperature of 19.72 °C, as highlighted in Figure 21. This unexpected result may suggest a higher water use efficiency or a different physiological adaptation to thermal stress under these conditions. The distinct thermal footprint evident in the heatmap adds a layer of depth to our agronomic understanding and indicates that canopy temperature, much like NDVI and SPAD indices, can provide valuable insights into the optimal balance of agronomic inputs for healthy plant development.
As the study progressed to CT2 and CT3 measurements during the stem extension and heading stages, the statistical significance of the interaction between nitrogen and seed rate decreased (F = 2.07, p = 0.093 for CT2 and F = 0.33, p = 0.856 for CT3), as shown in Table 7. This trend underscores the necessity of adjusting nitrogen and seed rates in the early growth stages to effectively modulate canopy temperature, which could influence water use efficiency and overall plant health.
The findings from NDVI, SPAD, and CT measurements demonstrate that strategic management of nitrogen application and seed density throughout the wheat’s lifecycle is crucial for optimizing yield and ensuring plant health. This comprehensive approach promotes precision in agronomic decision-making and supports sustainable practices that enhance crop development and quality. By integrating precise agronomic data, we advocate a management model that enhances crop yield and sustainability, establishing a framework for agronomic excellence and environmental stewardship.

4. Discussion

In wheat agriculture, the significant interaction between nitrogen doses and seeding rates, as demonstrated in our study, corroborates the growing consensus on the need for nuanced nutrient management. Such a management strategy not only maximizes yield, but also enhances other aspects of crop performance. Specifically, we found that an intermediate nitrogen dose (N2) paired with a medium seed rate (S2) leads to optimal wheat yield, aligning with sustainable agronomy principles that advocate for precision rather than excess input use. This synergy reflects a complex balance, whereby nitrogen availability and seed density act together to optimize growth conditions, potentially minimizing environmental impacts. Our findings are supported by the work of Arduini et al. [18], who evaluated grain yield, and dry matter and nitrogen accumulation and remobilization, in durum wheat as affected by variety and seeding rate, illustrating how different seeding rates can impact wheat’s nitrogen use and yield under various conditions. Furthermore, the research by Malkarnekar Saharsh et al. [19] on the response of nitrogen and plant growth regulators on the growth and yield of wheat highlights the importance of appropriate nitrogen doses in conjunction with other agronomic practices to enhance yield and crop performance.
Additionally, Iqbal et al. [6] demonstrated that the interaction between seeding rates and nitrogen levels significantly affected grain yield, emphasizing the critical role of balanced nutrient management in achieving high yield. This finding supports our results and highlights the importance of optimizing both nitrogen doses and seeding rates for sustainable and productive wheat agriculture. Moreover, Ecco et al. [20] found that specific nitrogen doses positively influenced the number of tillers and spike length up to a certain threshold, although it did not significantly affect grain yield. These insights underscore the complexity of nitrogen’s impact on different wheat growth parameters and further reinforce the need for precise nutrient management strategies to optimize both nitrogen doses and seeding rates for sustainable wheat production.
Extending beyond yield, our study underscores the importance of considering crop biomass as a key output parameter. The highest biomass was recorded with the same N2 and S2 combination, suggesting that factors optimizing yield also benefit the overall plant growth, a relationship that has been previously observed. This dual benefit is critical, with biomass serving not only as a measure of crop productivity but also as an essential resource for livestock nutrition. Trentin et al. [21] demonstrate the significant impact of nitrogen fertilizer dosing on durum wheat biomass production, underscoring the importance of nitrogen availability in enhancing biomass and its potential benefits for livestock nutrition (Trentin et al. [21]). Additionally, the study by Souissi et al. [22] evaluates the effects of nitrogen fertilization on the agronomic and economic performances of durum wheat in rainfed semi-arid environments similar to those of our study. Their findings indicate that specific nitrogen dosages can significantly enhance N-use efficiency and grain yield, emphasizing the importance of nuanced nutrient management strategies in optimizing wheat biomass for enhanced livestock nutrition (Souissi et al. [22]).
Furthermore, the research by Latiri-Souki et al. [23] highlights that nitrogen fertilizer can increase dry matter, grain production, and radiation and water use efficiencies for durum wheat under semi-arid conditions, reinforcing the critical role of nitrogen in biomass production. Additionally, the study by Mon et al. [24] found that the interaction between nitrogen fertilization and irrigation significantly influences grain yield, canopy temperature, and nitrogen use efficiency, illustrating the complex dynamics between nitrogen management and environmental factors in enhancing biomass and overall crop performance.
The relationship between nitrogen levels and seeding rates sheds light on the nuances of agronomic practices. Although the combined effect of nitrogen dose and seed rate on yield was not statistically significant, it was observed that lower nitrogen rates, particularly when coupled with the highest seed rate, profoundly affected protein content, with the lower nitrogen dose (N3) surprisingly correlating with higher grain protein. This finding aligns with recent research by Melash et al. [25] and Banach et al. [26], suggesting that, under certain conditions, reduced nitrogen may lead to more efficient protein synthesis and accumulation. This phenomenon might be indicative of a nitrogen-sparing effect where the crop adapts to limited nitrogen by diverting it towards grain protein synthesis, presenting a case for re-examining nitrogen recommendations for protein optimization.
Contrary to our findings, Subedi et al. [27] reported that grain protein concentration generally increases with higher nitrogen application, even when grain yield did not benefit from split nitrogen application. This underscores the commonly accepted importance of nitrogen management in achieving high protein content. Similarly, Ghimire et al. [28] found that nitrogen application significantly increased grain yield and protein content, with spring and split applications showing better results than fall application in years with a risk of nitrogen loss. These studies suggest that, under typical conditions, higher nitrogen levels are expected to enhance protein content, whereas our study presents an interesting deviation where lower nitrogen levels coupled with higher seed rates resulted in higher protein content.
Moving beyond the tangible outputs of yield and biomass, our study leverages advanced agricultural tools to offer deeper insights into plant health and management strategies, such as NDVI, which has emerged as an invaluable tool in agronomy for assessing plant health and vigor across growth stages. Our study’s NDVI results offer compelling evidence of the role of early growth stages in determining the overall health and productivity of wheat crops. The pronounced growth observed during the tillering phase, under the highest nitrogen dose (N1), underscores the foundational importance of adequate nitrogen application for setting the stage for healthy plant development. This early vegetative vigor, as indicated by NDVI1, establishes a robust foundation that is strongly correlated with biomass accumulation (r = 0.65, p < 0.001) and has a significant impact on yield (r = 0.38, p < 0.001). Such findings corroborate the hypothesis that efficient nitrogen management at initial growth stages is paramount for optimizing subsequent crop outputs. Insights from relevant studies further reinforce the significance of precise nitrogen application as supported by NDVI data. The work of Vian et al. [29], demonstrates how NDVI, measured by an active optical canopy sensor, can be used to optimize nitrogen topdressing doses, enabling variable-rate nitrogen fertilization and affirming NDVI’s role in facilitating efficient nitrogen management. Additionally, Kizilgeci et al. [1] support the utility of NDVI in identifying nitrogen deficiency and enhancing precision nitrogen management in durum wheat cultivars under semi-arid conditions, showcasing its potential in ensuring food security amid climate change scenarios.
As the crop advances to the stem extension phase, the sustained high NDVI2 values reflect the continued vegetative health of the plants, highlighting the enduring effects of the initial nitrogen application. This phase marks a critical juncture where the maintained vigor from tillering supports the plant through its preparation for reproductive development. It is a testament to the necessity of reevaluating nitrogen application strategies to maintain plant health throughout the crop cycle, ensuring that the growth is not only sustained but also optimized for the reproductive phase. Li et al. [30] provide evidence supporting this notion. Their research indicates that adjusting nitrogen usage can significantly enhance grain yield and nitrogen utilization efficiency in wheat crops. Such findings suggest a crucial role of informed nitrogen management in sustaining plant health and optimizing production as the crop progresses through various developmental stages. Furthermore, Gezahegn et al. [31] demonstrated that split applications of nitrogen at specific growth stages significantly improved the yield and quality of durum wheat, reinforcing the importance of strategic nitrogen management throughout the crop cycle.
However, NDVI’s predictive capacity appears to wane as the plant transitions to the heading stage, evidenced by the declining NDVI3 values across all nitrogen treatments. This trend is indicative of the plant’s physiological shift from vegetative growth to reproductive development and grain filling, as demonstrated by Zhang et al. [32], who explored the impacts of optimizing fertilization on soil nitrogen cycling and wheat nitrogen utilization. Their findings underline the importance of adjusting nitrogen application to enhance wheat’s nitrogen uptake efficiency and reduce loss, shedding light on the relationship between nitrogen management and plant maturity stages [32]. The diminished correlation between NDVI3 with both yield and biomass as plants advance to the heading stage signifies a transition in plant physiology. This is further exemplified by Hussain et al. [33], who investigated how different nitrogen fertilizers impact grain yield and agronomic nitrogen use efficiency in wheat cultivars. Their work highlights how nitrogen source selection can significantly influence yield parameters and efficiency, emphasizing the necessity to tailor nitrogen management strategies to the specific developmental needs of the plant as it matures [33]. Additionally, Ibarra-Villarreal et al. [34] found that specific nitrogen rates can influence NDVI values and canopy temperature, further illustrating the nuanced impacts of nitrogen management on crop development.
Interestingly, the lack of a significant correlation between NDVI and protein content across all growth stages sheds light on the complex dynamics that govern grain quality traits. This disconnect suggests that factors influencing protein synthesis may extend beyond the vegetative growth parameters captured by NDVI, possibly encompassing a broader range of genetic and environmental influences. Studies by Dashkevich et al. [35], on the genetic potential of spring durum wheat, and by De Santis et al. [36], reviewing the influence of drought and abiotic stress on grain quality in Mediterranean environments, suggest such complexities. Dashkevich et al. revealed significant variations in protein content across durum wheat genotypes, underscoring the genetic underpinnings of quality traits, while De Santis et al. highlighted the role of environmental stressors in affecting protein content and composition, further emphasizing the multifactorial nature of grain quality. These insights indicate the profound implications for breeding and management practices, revealing the potential for targeted agronomic interventions to enhance grain quality without adversely affecting yield. Additionally, Mefleh et al. [37] found that genotypic variation in grain nitrogen content correlated with the content of different protein fractions, further emphasizing the genetic factors in protein synthesis. Afzal et al. [38] also highlighted the variability in protein profiles among wheat cultivars and the significant impact of environmental conditions on protein expression.
In synthesizing these insights, our study proposes a holistic approach to nitrogen and seed rate management that transcends traditional yield and biomass maximization objectives. By integrating NDVI-based assessments into agronomic decision-making, we can leverage early indicators of plant health to inform targeted interventions that enhance crop quantity. This strategy not only aligns with sustainable agricultural practices but also promotes environmental stewardship by optimizing input use and minimizing waste. The substantial economic benefits demonstrated by Finco et al. [39], through a decrease in labor and pesticide costs while enhancing nitrogen and seed distribution efficiencies, highlight the tangible outcomes of integrating precision agriculture in durum wheat production. Moreover, the potential for using NDVI as an effective tool for predicting yield and managing nitrogen applications in durum wheat is further supported by Wang et al. [40]. Their findings on crop models and canopy reflectance index, specifically NDVI at flowering time, as beneficial methods for managing nitrogen applications, echo our stance on the value of precision agriculture technologies. Together, these studies underscore the transformative potential of NDVI assessments and precision agriculture in fostering more nuanced crop management decisions, thereby advancing agronomic efficiency and sustainability in wheat production. Additionally, Toscano et al. [8] demonstrated the significant and positive linear relationships between NDVI and yield monitoring data, explaining most of the within-field variability in durum wheat, and highlighting the effectiveness of NDVI in precision agriculture. Vian et al. [29] also found that NDVI can be used to develop models for estimating shoot biomass and nitrogen content, which can help determine optimal nitrogen topdressing doses, thus improving nitrogen use efficiency.
Collectively, the findings underscore the transformative potential of integrating NDVI data into crop management strategies. This approach not only enhances our understanding of the dynamic interactions between nitrogen application, seed density, and plant growth, but also propels us toward more sustainable and efficient agricultural practices. Future research should continue to explore the multifaceted relationships between these agronomic factors and crop performance, expanding the scope of precision agriculture to encompass a broader range of crops and environmental conditions. The study by Santaga et al. [41] further supports this, showing that advanced precision nitrogen fertilization models limited yield losses and reduced intra-field variability under less favorable conditions. Mitra et al. [42] emphasized the effectiveness of NDVI sensor-based nitrogen management strategies in significantly improving wheat yield, nitrogen use efficiencies, and economic returns, validating the integration of advanced technologies in agronomic practices.
As we further explore the physiological aspects of wheat growth, the SPAD and CT results shed light on the complexity of plant health beyond conventional metrics. The SPAD results are indicative of chlorophyll content and, by proxy, plant health. The findings from Fiorentini et al. [43], who evaluated the relationship between SPAD readings, chlorophyll concentration, and Nitrogen Nutrition Index in durum wheat under various agricultural practices, demonstrated a strong association between SPAD readings, chlorophyll concentration, and efficient nitrogen management, supporting the essential role of nitrogen in establishing initial plant vigor, which correlates with overall plant productivity. Additionally, research by Sharma et al. [44] assessed the impact of different nitrogen and phosphorus fertilizer rates on chlorophyll content in winter wheat varieties, reinforcing the importance of nutrient management strategies for optimal plant growth reflected in SPAD values.
Moreover, the interaction between nitrogen dose and seed rate significantly impacts SPAD3 measurements, underscoring the nuanced nutrient utilization strategies as plants approach reproductive maturity. The findings from Skudra and Ruza [45] demonstrate that adjusted levels of nitrogen and seed rate markedly affect wheat growth and productivity, indicating the importance of these factors in enhancing chlorophyll content for sustained photosynthesis and grain development. This adaptation in agronomic practices, informed by precise nitrogen management and seeding strategies, is essential for maintaining plant vigor and optimizing wheat yield. Additionally, Jhanji and Sekhon [46] highlighted that chlorophyll meters could effectively quantify chlorophyll and nitrogen content in wheat leaves, showing significant correlations with the SPAD index, thereby enhancing the accuracy of nitrogen management practices. Furthermore, Wang et al. [47] found that multiple SPAD measurements on the same plants can improve the estimation of crop nitrogen status, further supporting the use of SPAD meters for effective nitrogen management in wheat.
As we examine the physiological intricacies of durum wheat development, the CT results from this study offer a fascinating glimpse into the thermal dynamics influencing plant growth. The observed coolest mean CT1 temperatures under conditions of high nitrogen application and low seed rate suggest an ideal thermal environment that may enhance growth due to reduced intra-species competition and optimized water usage. This observation resonates with the findings of Al-Karaki et al. [48]. Their work emphasizes the influence of thermal conditions on durum wheat, particularly how optimal temperatures can drive growth efficiency. Furthermore, as the wheat matures, the interaction between nitrogen and seeding rates diminishes, indicating a shift in the plant’s physiological needs. This phase of growth correlates with insights from Marti and Slafer et al. [49], who explored how durum wheat yields respond to a broad spectrum of environmental factors, highlighting the adaptive nature of plant responses to abiotic stress. Additionally, Cossani and Sadras [50] found that the interplay of nitrogen and water availability can influence how elevated temperature impacts wheat yield. Hou et al. [51] investigated the impact of experimental warming on nitrogen uptake in winter wheat, highlighting the significant effects of temperature and nitrogen management on crop performance.
In light of the relationships between nitrogen management, crop performance, and the incorporation of precision agriculture technologies, our study propels the discourse towards redefining conventional nitrogen management paradigms. The observed inverse relationship between reduced nitrogen availability and increased protein content necessitates a critical reassessment of nitrogen strategies. Such a phenomenon suggests that wheat possesses an inherent adaptive mechanism, potentially reallocating limited nitrogen towards grain protein synthesis under nitrogen-limited conditions. This insight, echoing findings from contemporary agronomic research, underscores the necessity for a recalibrated approach to fertilization that harmonizes yield optimization with grain quality enhancement. Abedi et al. [52] report that both the rate and timing of nitrogen fertilization play crucial roles in enhancing wheat protein content and its quality, alongside high wheat production, thus supporting a strategic rethink in nitrogen application strategies. Additionally, Li et al. [53] demonstrated that applying remedial nitrogen prior to low-temperature stress was more effective in enhancing wheat morphology and nitrogen uptake efficiency compared to post-low-temperature stress applications, highlighting the importance of strategic nitrogen management in varying environmental conditions.
The consideration of trade-offs in agricultural management—particularly between yield enhancement and protein concentration—highlights the complex balance required in crop management. This balance, reflecting a broader physiological response observed across crop species to varying nitrogen levels, advocates for a strategic approach to fertilization. Such an approach aims to balance nitrogen inputs to optimize both yield and grain quality, contributing to the evolving narrative on sustainable crop management practices. Recent studies, like that conducted by Karatay et al. [54], show the profitability and risk considerations of site-specific nitrogen management strategies in wheat production, which not only aim at enhancing grain quality and yield, but also consider price premiums associated with higher quality grains under different nitrogen management scenarios.
Furthermore, the integration of technologies such as NDVI and SPAD into our agronomic framework illuminates the path towards refined crop management strategies. By harnessing real-time insights into plant health and nitrogen status, these precision agriculture tools offer the potential to tailor agronomic interventions more closely to the plant’s lifecycle needs, thereby increasing both yield and quality outcomes. For instance, research by Kizilgeci et al. [1] highlights the effectiveness of SPAD and NDVI in assessing and optimizing nitrogen management strategies in durum wheat, demonstrating their critical role in improving grain yield and nutrient efficiency under semi-arid conditions. This supports the importance of these technologies in adapting agronomic practices to meet specific crop needs, enhancing sustainability and productivity in durum wheat cultivation.
Equally significant is our exploration of CT, which underscores the critical role of water use efficiency within the context of crop management. The strategic manipulation of agronomic factors to maintain optimal canopy temperatures shows the importance of adaptive management strategies in ensuring crop productivity, especially in water-scarce regions like the location of our study, Sidi El Aidi in Morocco. This facet of our findings reinforces the utility of precision agriculture in navigating the modern farming challenges, highlighting the importance of agronomic decisions in enhancing crop resilience and sustainability. In this context, the research by Devkota et al. [55] provides an analysis of genotype and agronomic management interactions that enhance wheat yield and water use efficiency in the Mediterranean rainfed environment of Morocco. Their study specifically examines how different genotypes perform under varied agronomic practices, shedding light on the optimal strategies to improve both yield and water use efficiency under the challenging conditions typical of Moroccan agriculture.
Our study underscores a crucial shift towards a more nuanced, stage-specific approach to crop management but also emphasizes the significance of integrating precision agriculture technologies to advance sustainable farming practices. By examining the adaptive responses of wheat to nitrogen availability and leveraging new technologies, we pave the way for future research aimed at fostering agricultural systems that address the pressing needs of food security while upholding environmental stewardship. As we continue to unravel the interplay between agronomic factors and crop performance, our collective efforts will be instrumental in shaping the future of sustainable agriculture. In particular, the study by Finco et al. [39] explores the economic outcomes of precision agriculture investments in durum wheat production, demonstrating how sustainable farming practices, combined with precision farming techniques, significantly enhance the sustainability and profitability of durum wheat cultivation, thereby supporting the broader goals of sustainable agricultural development.

5. Conclusions

This study examined the effects of various nitrogen doses and seeding rates on wheat yield, biomass, and protein content, employing precision agriculture technologies such as NDVI, SPAD, and canopy temperature measurements. Strategic nitrogen management, especially during the tillering phase as revealed by NDVI, is crucial for initiating robust early plant growth, which significantly impacts overall biomass and yield. The analysis demonstrated that an intermediate nitrogen dose combined with a moderate seed rate optimally boosts yield and biomass, showcasing the effectiveness of precision agriculture in refining crop management practices. NDVI data also highlighted a decline in vegetative vigor as plants matured, necessitating adjustments in nitrogen strategies to accommodate the plants’ evolving physiological needs. Furthermore, SPAD measurements confirmed that higher nitrogen levels substantially enhance chlorophyll content during early growth stages, supporting plant health. However, this influence wanes in later stages, indicating a physiological shift toward reproductive development. Canopy temperature data underscored that proper management of nitrogen and seed rates early in the growth cycle plays a vital role in modulating plant stress responses and enhancing water use efficiency. Together, these findings advocate for finely tuned nitrogen and seeding strategies to not only enhance crop yield, but also to promote sustainable agricultural practices. This integrated approach underscores the potential of precision agriculture to significantly improve both the efficiency and sustainability of wheat production, aligning advanced agronomic tools with sustainable farming objectives.

Author Contributions

The following contributions were made by the authors to the study and preparation of this manuscript: O.H. was responsible for the Writing—Original Draft Preparation and Visualization of the data, A.B. was responsible for Project Administration, Supervision, and providing Resources, I.K. contributed to Validation of the study, A.L. contributed to Validation and Funding Acquisition, N.N. was responsible for the Methodology and Validation of the study and A.A. was responsible for the Writing—Review and Editing, Supervision, Data Curation, and providing Resources. All authors have read and agreed to the published version of the manuscript.

Funding

The paper is part of the study conducted on Agricultural Innovation, Resilience, and Income Improvement in the Case of Wheat Value Chain (RESIFARM) Project. The project is funded by UM6P under the APRD Program and managed in collaboration with CNRST, MESRSI and technically coordinated by INRA Morocco.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to ongoing analyses and intended use in future publications.

Acknowledgments

The authors express gratitude to their institutional affiliates for their continual support, as well as their colleagues for their insightful contributions to this study. Appreciation is also extended to the reviewers and the MDPI Plants editorial team for their guidance in manuscript improvement.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Kizilgeci, F.; Yildirim, M.; Islam, M.S.; Ratnasekera, D.; Iqbal, M.A.; Sabagh, A.E. Normalized difference vegetation index and chlorophyll content for precision nitrogen management in durum wheat cultivars under semi-arid conditions. Sustainability 2021, 13, 3725. [Google Scholar] [CrossRef]
  2. Yin, Q.; Zhang, Y.; Li, W.; Wang, J.; Wang, W.; Ahmad, I.; Zhou, G.; Huo, Z. Better inversion of wheat canopy SPAD values before heading stage using spectral and texture indices based on UAV multispectral imagery. Remote Sens. 2023, 15, 4935. [Google Scholar] [CrossRef]
  3. Cavalaris, C.; Megoudi, S.; Maxouri, M.; Anatolitis, K.; Sifakis, M.; Levizou, E.; Kyparissis, A. Modeling of durum wheat yield based on sentinel-2 imagery. Agronomy 2021, 11, 1486. [Google Scholar] [CrossRef]
  4. Pampana, S.; Mariotti, M. Durum wheat yield and N uptake as affected by N source, timing, and rate in two mediterranean environments. Agronomy 2021, 11, 1299. [Google Scholar] [CrossRef]
  5. Tedone, L.; Ali, S.A.; Verdini, L.; De Mastro, G. Nitrogen management strategy for optimizing agronomic and environmental performance of rainfed durum wheat under Mediterranean climate. J. Clean. Prod. 2018, 172, 2058–2074. [Google Scholar] [CrossRef]
  6. Iqbal, J.; Hayat, K.; Hussain, S.; Ali, A.; Bakhsh, M.A.A.H.A. Effect of seeding rates and nitrogen levels on yield and yield components of wheat (Triticum aestivum L.). Pak. J. Nutr. 2012, 11, 531. [Google Scholar] [CrossRef]
  7. Tedone, L.; Verdini, L.; Grassano, N.; Tarraf, W.; De Mastro, G. Optimising nitrogen in order to improve the efficiency, eco-physiology, yield and quality on one cultivar of durum wheat. Ital. J. Agron. 2014, 9, 49–54. [Google Scholar] [CrossRef]
  8. Toscano, P.; Castrignanò, A.; Di Gennaro, S.F.; Vonella, A.V.; Ventrella, D.; Matese, A. A precision agriculture approach for durum wheat yield assessment using remote sensing data and yield mapping. Agronomy 2019, 9, 437. [Google Scholar] [CrossRef]
  9. Messina, G.; Badagliacca, G.; Praticò, S.; Preiti, G.; Monti, M.; Modica, G. Multispectral UAV-Based Monitoring of Behavior of Different Wheat and Barley Varieties. In Proceedings of the Conference of the Italian Society of Agricultural Engineering, Palermo, Italy, 19–22 September 2022; Springer International Publishing: Cham, Switzerland, 2022; pp. 1173–1181. [Google Scholar]
  10. Sultana, S.R.; Ali, A.; Ahmad, A.; Mubeen, M.; Zia-Ul-Haq, M.; Ahmad, S.; Ercisli, S.; Jaafar, H.Z. Normalized difference vegetation index as a tool for wheat yield estimation: A case study from Faisalabad, Pakistan. Sci. World J. 2014, 2014, 725326. [Google Scholar] [CrossRef]
  11. Vizzari, M.; Santaga, F.; Benincasa, P. Sentinel 2-based nitrogen VRT fertilization in wheat: Comparison between traditional and simple precision practices. Agronomy 2019, 9, 278. [Google Scholar] [CrossRef]
  12. Aranguren, M.; Castellón, A.; Aizpurua, A. Wheat yield estimation with NDVI values using a proximal sensing tool. Remote Sens. 2020, 12, 2749. [Google Scholar] [CrossRef]
  13. Rodrigues, F.A., Jr.; Blasch, G.; Defourny, P.; Ortiz-Monasterio, J.I.; Schulthess, U.; Zarco-Tejada, P.J.; Taylor, J.A.; Gérard, B. Multi-temporal and spectral analysis of high-resolution hyperspectral airborne imagery for precision agriculture: Assessment of wheat grain yield and grain protein content. Remote Sens. 2018, 10, 930. [Google Scholar] [CrossRef] [PubMed]
  14. Badagliacca, G.; Messina, G.; Praticò, S.; Lo Presti, E.; Preiti, G.; Monti, M.; Modica, G. Multispectral Vegetation Indices and Machine Learning Approaches for Durum Wheat (Triticum durum Desf.) Yield Prediction across Different Varieties. AgriEngineering 2023, 5, 2032–2048. [Google Scholar] [CrossRef]
  15. Didal, V.; Vidyasagar, G.; Mahender Kumar, R.; Surekha, K.; Narender Reddy, S.; Bhooshan, B. Effect of nitrogen management practices on SPAD values and NDVI readings of rice crop. Pharma Innov. J 2022, 11, 367–371. [Google Scholar]
  16. Vannoppen, A.; Gobin, A. Estimating farm wheat yields from NDVI and meteorological data. Agronomy 2021, 11, 946. [Google Scholar] [CrossRef]
  17. 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]
  18. Arduini, I.; Masoni, A.; Ercoli, L.; Mariotti, M. Grain yield, and dry matter and nitrogen accumulation and remobilization in durum wheat as affected by variety and seeding rate. Eur. J. Agron. 2006, 25, 309–318. [Google Scholar] [CrossRef]
  19. Saharsh, M.; Singh, R.; Indu, T. Response of Nitrogen and Plant Growth Regulators on Growth and Yield of Wheat (Triticum aestivum L.). Int. J. Plant Soil Sci. 2023, 35, 66–73. [Google Scholar] [CrossRef]
  20. Ecco, M.; Perin, R.F.; Duarte, R.P.; Barbosa, A.P.; Gabriel, A.M. Biometric responses of wheat plants submitted to different does of nitrogen. Rev. Cient. Rural. 2020, 22, 169–184. [Google Scholar]
  21. Trentin, C.; Bredemeier, C.; Vian, A.L.; Drum, M.A.; Santos, F.L.d. Biomass production and wheat grain yield and its relationship with NDVI as a function of nitrogen availability. Rev. Bras. Ciências Agrárias Recife. 2021, 16, e34. [Google Scholar] [CrossRef]
  22. Souissi, A.; Bahri, H.; Cheikh M’hamed, H.; Chakroun, M.; Benyoussef, S.; Frija, A.; Annabi, M. Effect of tillage, previous crop, and N fertilization on agronomic and economic performances of durum wheat (Triticum durum desf.) under rainfed semi-arid environment. Agronomy 2020, 10, 1161. [Google Scholar] [CrossRef]
  23. Latiri-Souki, K.; Nortcliff, S.; Lawlor, D. Nitrogen fertilizer can increase dry matter, grain production and radiation and water use efficiencies for durum wheat under semi-arid conditions. Eur. J. Agron. 1998, 9, 21–34. [Google Scholar] [CrossRef]
  24. Mon, J.; Bronson, K.; Hunsaker, D.; Thorp, K.; White, J.; French, A. Interactive effects of nitrogen fertilization and irrigation on grain yield, canopy temperature, and nitrogen use efficiency in overhead sprinkler-irrigated durum wheat. Field Crops Res. 2016, 191, 54–65. [Google Scholar] [CrossRef]
  25. Melash, A.A.; Bogale, A.A.; Bytyqi, B.; Nyandi, M.S.; Ábrahám, É.B. Nutrient management: As a panacea to improve the caryopsis quality and yield potential of durum wheat (Triticum turgidum L.) under the changing climatic conditions. Front. Plant Sci. 2023, 14, 1232675. [Google Scholar] [CrossRef]
  26. Banach, J.K.; Majewska, K.; Żuk-Gołaszewska, K. Effect of cultivation system on quality changes in durum wheat grain and flour produced in North-Eastern Europe. PLoS ONE 2021, 16, e0236617. [Google Scholar] [CrossRef]
  27. Subedi, K.; Ma, B.; Xue, A. Planting date and nitrogen effects on grain yield and protein content of spring wheat. Crop Sci. 2007, 47, 36–44. [Google Scholar] [CrossRef]
  28. Ghimire, D.; Das, S.; Mueller, N.D.; Creech, C.F.; Santra, D.; Baenziger, P.S.; Easterly, A.C.; Maust, B.; Maharjan, B. Effects of cultivars and nitrogen management on wheat grain yield and protein. Agron. J. 2021, 113, 4348–4368. [Google Scholar] [CrossRef]
  29. Vian, A.L.; Bredemeier, C.; Turra, M.A.; Giordano, C.P.d.S.; Fochesatto, E.; Silva, J.A.d.; Drum, M.A. Nitrogen management in wheat based on the normalized difference vegetation index (NDVI). Ciência Rural. 2018, 48, e20170743. [Google Scholar] [CrossRef]
  30. Li, S.; Li, Z.; Bi, X.; Feng, B.; Wang, Z.; Wang, F.; Si, J.; Shi, J.; Liu, K. Nitrogen fertilizer management on wheat yield and nitrogen utilization. J. Plant Nutr. 2022, 45, 1953–1960. [Google Scholar] [CrossRef]
  31. Gezahegn, A.M.; Tafes Desta, B.; Takel, A.; Eshetu, S. Split application of nitrogen fertilizer for optimum yield of durum wheat (Triticum turgidum l. var. durum) at central ethiopia. Commun. Soil Sci. Plant Anal. 2022, 53, 1809–1822. [Google Scholar] [CrossRef]
  32. Zhang, Z.; Yu, Z.; Zhang, Y.; Shi, Y. Impacts of fertilization optimization on soil nitrogen cycling and wheat nitrogen utilization under water-saving irrigation. Front. Plant Sci. 2022, 13, 878424. [Google Scholar] [CrossRef]
  33. Hussain, M.; Cheema, S.A.; Abbas, R.Q.; Ashraf, M.F.; Shahzad, M.; Farooq, M.; Jabran, K. Choice of nitrogen fertilizer affects grain yield and agronomic nitrogen use efficiency of wheat cultivars. J. Plant Nutr. 2018, 41, 2330–2343. [Google Scholar] [CrossRef]
  34. Ibarra-Villarreal, A.L.; Villarreal-Delgado, M.F.; Parra-Cota, F.I.; Yepez, E.A.; Guzmán, C.; Gutierrez-Coronado, M.A.; Valdez, L.C.; Saint-Pierre, C.; Santos-Villalobos, S.d.L. Effect of a native bacterial consortium on growth, yield, and grain quality of durum wheat (Triticum turgidum L. subsp. durum) under different nitrogen rates in the Yaqui Valley, Mexico. Plant Signal. Behav. 2023, 18, 2219837. [Google Scholar] [CrossRef]
  35. Dashkevich, S.; Utebayev, M.; Kradetskaya, O.; Chilimova, I.; Zhylkybaev, R.; Babkenov, A. The Geneti c Potential of Spring Durum Wheat Grain Quality in the Nort h of Kazakhstan. OnLine J. Biol. Sci 2022, 22, 347–355. [Google Scholar] [CrossRef]
  36. De Santis, M.A.; Soccio, M.; Laus, M.N.; Flagella, Z. Influence of drought and salt stress on durum wheat grain quality and composition: A review. Plants 2021, 10, 2599. [Google Scholar] [CrossRef]
  37. Mefleh, M.; Motzo, R.; Samson, M.-F.; Morel, M.-H.; Giunta, F. N partitioning between gluten fractions in a set of Italian durum wheat cultivars: Role of the grain n content. Foods 2020, 9, 1684. [Google Scholar] [CrossRef]
  38. Afzal, M.; Sielaff, M.; Curella, V.; Neerukonda, M.; El Hassouni, K.; Schuppan, D.; Tenzer, S.; Longin, C.F.H. Characterization of 150 wheat cultivars by LC-MS-based label-free quantitative proteomics unravels possibilities to design wheat better for baking quality and human health. Plants 2021, 10, 424. [Google Scholar] [CrossRef]
  39. Finco, A.; Bucci, G.; Belletti, M.; Bentivoglio, D. The economic results of investing in precision agriculture in durum wheat production: A case study in central Italy. Agronomy 2021, 11, 1520. [Google Scholar] [CrossRef]
  40. Wang, G.S.; Ottman, M.; Thorp, K. Late N management in durum wheat using crop models and canopy reflectance. Forage Grain A Coll. Agric. Life Sci. Rep. 2012, 65–72. [Google Scholar]
  41. Santaga, F.S.; Benincasa, P.; Toscano, P.; Antognelli, S.; Ranieri, E.; Vizzari, M. Simplified and advanced Sentinel-2-based precision nitrogen management of wheat. Agronomy 2021, 11, 1156. [Google Scholar] [CrossRef]
  42. Mitra, B.; Singha, P.; Roy Chowdhury, A.; Sinha, A.K.; Skalicky, M.; Laing, A.M.; Alamri, S.; Hossain, A. Optical sensor-based nitrogen management: An environmentally friendly and cost-effective approach for sustainable wheat (Triticum aestivum L.) production on Eastern plains of India. Front. Sustain. Food Syst. 2023, 7, 1153575. [Google Scholar] [CrossRef]
  43. 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]
  44. Sharma, P.K.; Gurjar, R.; Saini, R.; Meena, D.; Meena, S.; Sabal, S.K. Response of Different Physiochemical Parameters of Wheat Genotypes (Triticum aestivum L. & Triticum durum L.) in Nutrient Stress Condition. Int. J. Plant Soil Sci 2022, 34, 1074–1083. [Google Scholar]
  45. Skudra, I.; Ruza, A. Effect of nitrogen and sulphur fertilization on chlorophyll content in winter wheat. Rural. Sustain. Res. 2017, 37, 29–37. [Google Scholar] [CrossRef]
  46. Jhanji, S.; Sekhon, N.K. Evaluation of potential of portable chlorophyll meter to quantify chlorophyll and nitrogen contents in leaves of wheat under different field conditions. Indian J. Exp. Biol. 2018, 56, 750–758. [Google Scholar]
  47. Wang, G.; Bronson, K.F.; Thorp, K.R.; Mon, J.; Badaruddin, M. Multiple leaf measurements improve effectiveness of chlorophyll meter for durum wheat nitrogen management. Crop Sci. 2014, 54, 817–826. [Google Scholar] [CrossRef]
  48. Al-Karaki, G.N. Phenological development-yield relationships in durum wheat cultivars under late-season high-temperature stress in a semiarid environment. Int. Sch. Res. Not. 2012, 2012, 456856. [Google Scholar] [CrossRef]
  49. Marti, J.; Slafer, G.A. Bread and durum wheat yields under a wide range of environmental conditions. Field Crops Res. 2014, 156, 258–271. [Google Scholar] [CrossRef]
  50. Cossani, C.M.; Sadras, V.O. Nitrogen and water supply modulate the effect of elevated temperature on wheat yield. Eur. J. Agron. 2021, 124, 126227. [Google Scholar] [CrossRef]
  51. Hou, R.; Xu, X.; Ouyang, Z. Effect of experimental warming on nitrogen uptake by winter wheat under conventional tillage versus no-till systems. Soil Tillage Res. 2018, 180, 116–125. [Google Scholar] [CrossRef]
  52. Abedi, T.; Alemzadeh, A.; Kazemeini, S.A. Wheat yield and grain protein response to nitrogen amount and timing. Aust. J. Crop Sci. 2011, 5, 330–336. [Google Scholar]
  53. Li, C.; Liu, M.; Dai, C.; Zhu, Y.; Zhu, M.; Ding, J.; Zhu, X.; Zhou, G.; Guo, W. Morphology and nitrogen uptake and distribution of wheat plants as influenced by applying remedial urea prior to or post low-temperature stress at seedling stage. Agronomy 2022, 12, 2338. [Google Scholar] [CrossRef]
  54. Meyer-Aurich, A.; Karatay, Y.N.; Nausediene, A.; Kirschke, D. Effectivity and cost efficiency of a tax on nitrogen fertilizer to reduce GHG emissions from agriculture. Atmosphere 2020, 11, 607. [Google Scholar] [CrossRef]
  55. Devkota, M.; Devkota, K.P.; Karrou, M.; Nangia, V. Genotype and agronomic management interaction to enhance wheat yield and water use efficiency in the Mediterranean rainfed environment of Morocco: I. Field data analysis. Eur. J. Agron. 2023, 151, 126972. [Google Scholar] [CrossRef]
Figure 1. Implementation of precision measurement techniques in our field research. (a) Use of a multispectral radiometer for NDVI measurement. (b) SPAD-502 Plus meter used for assessing chlorophyll content. (c) Infrared thermometer for canopy temperature analysis.
Figure 1. Implementation of precision measurement techniques in our field research. (a) Use of a multispectral radiometer for NDVI measurement. (b) SPAD-502 Plus meter used for assessing chlorophyll content. (c) Infrared thermometer for canopy temperature analysis.
Plants 13 01574 g001
Figure 2. Heatmap of wheat mean yield variation by nitrogen dose (N1: 120 kg/ha, N2: 60 kg/ha, N3: 0 kg/ha) and seed rate (S1: 500 seeds/m2, S2: 400 seeds/m2, S3: 300 seeds/m2).
Figure 2. Heatmap of wheat mean yield variation by nitrogen dose (N1: 120 kg/ha, N2: 60 kg/ha, N3: 0 kg/ha) and seed rate (S1: 500 seeds/m2, S2: 400 seeds/m2, S3: 300 seeds/m2).
Plants 13 01574 g002
Figure 3. Boxplot of wheat yield distribution by nitrogen application levels (N1: 120 kg/ha, N2: 60 kg/ha, N3: 0 kg/ha) at the mean seed rate (S2: 400 seeds/m2). The black diamond symbol indicates outliers. Significance levels were determined by one-way ANOVA and subsequent Tukey test. Different letters indicate significant differences among the nitrogen application levels.
Figure 3. Boxplot of wheat yield distribution by nitrogen application levels (N1: 120 kg/ha, N2: 60 kg/ha, N3: 0 kg/ha) at the mean seed rate (S2: 400 seeds/m2). The black diamond symbol indicates outliers. Significance levels were determined by one-way ANOVA and subsequent Tukey test. Different letters indicate significant differences among the nitrogen application levels.
Plants 13 01574 g003
Figure 4. Heatmap of wheat biomass distribution by nitrogen dose (N1: 120 kg/ha, N2: 60 kg/ha, N3: 0 kg/ha) and seed rate (S1: 500 seeds/m2, S2: 400 seeds/m2, S3: 300 seeds/m2).
Figure 4. Heatmap of wheat biomass distribution by nitrogen dose (N1: 120 kg/ha, N2: 60 kg/ha, N3: 0 kg/ha) and seed rate (S1: 500 seeds/m2, S2: 400 seeds/m2, S3: 300 seeds/m2).
Plants 13 01574 g004
Figure 5. Boxplot of wheat biomass distribution by nitrogen application levels (N1: 120 kg/ha, N2: 60 kg/ha, N3: 0 kg/ha) at the mean seed rate (S2: 400 seeds/m2). The black diamond symbol indicates outliers. Significance levels were determined by one-way ANOVA and subsequent Tukey test. Different letters indicate significant differences among the nitrogen application levels.
Figure 5. Boxplot of wheat biomass distribution by nitrogen application levels (N1: 120 kg/ha, N2: 60 kg/ha, N3: 0 kg/ha) at the mean seed rate (S2: 400 seeds/m2). The black diamond symbol indicates outliers. Significance levels were determined by one-way ANOVA and subsequent Tukey test. Different letters indicate significant differences among the nitrogen application levels.
Plants 13 01574 g005
Figure 6. Boxplot of wheat biomass distribution by seed rate levels (S1: 500 seeds/m2, S2: 400 seeds/m2, S3: 300 seeds/m2) at a mean nitrogen application (N2: 60 kg/ha). Significance levels were determined by one-way ANOVA and subsequent Tukey test. Different letters indicate significant differences among the nitrogen application levels.
Figure 6. Boxplot of wheat biomass distribution by seed rate levels (S1: 500 seeds/m2, S2: 400 seeds/m2, S3: 300 seeds/m2) at a mean nitrogen application (N2: 60 kg/ha). Significance levels were determined by one-way ANOVA and subsequent Tukey test. Different letters indicate significant differences among the nitrogen application levels.
Plants 13 01574 g006
Figure 7. Heatmap of mean protein content in wheat by nitrogen dose (N1: 120 kg/ha, N2: 60 kg/ha, N3: 0 kg/ha) and seed rate (S1: 500 seeds/m2, S2: 400 seeds/m2, S3: 300 seeds/m2).
Figure 7. Heatmap of mean protein content in wheat by nitrogen dose (N1: 120 kg/ha, N2: 60 kg/ha, N3: 0 kg/ha) and seed rate (S1: 500 seeds/m2, S2: 400 seeds/m2, S3: 300 seeds/m2).
Plants 13 01574 g007
Figure 8. Variation of protein content in wheat as influenced by different nitrogen dose levels (N1: 120 kg/ha, N2: 60 kg/ha, N3: 0 kg/ha) at the mean seed rate (S2: 400 seeds/m2). The black diamond symbol indicates outliers. Significance levels were determined by one-way ANOVA and subsequent Tukey test. Different letters indicate significant differences among the nitrogen application levels.
Figure 8. Variation of protein content in wheat as influenced by different nitrogen dose levels (N1: 120 kg/ha, N2: 60 kg/ha, N3: 0 kg/ha) at the mean seed rate (S2: 400 seeds/m2). The black diamond symbol indicates outliers. Significance levels were determined by one-way ANOVA and subsequent Tukey test. Different letters indicate significant differences among the nitrogen application levels.
Plants 13 01574 g008
Figure 9. Protein content variability in wheat by seed rate levels (S1: 500 seeds/m2, S2: 400 seeds/m2, S3: 300 seeds/m2) at a mean nitrogen application (N2: 60 kg/ha). Significance levels were determined by one-way ANOVA and subsequent Tukey test. Different letters indicate significant differences among the nitrogen application levels.
Figure 9. Protein content variability in wheat by seed rate levels (S1: 500 seeds/m2, S2: 400 seeds/m2, S3: 300 seeds/m2) at a mean nitrogen application (N2: 60 kg/ha). Significance levels were determined by one-way ANOVA and subsequent Tukey test. Different letters indicate significant differences among the nitrogen application levels.
Plants 13 01574 g009
Figure 10. Heatmap of mean Normalized Difference Vegetation Index (NDVI1) by nitrogen dose (N1: 120 kg/ha, N2: 60 kg/ha, N3: 0 kg/ha) and seed rate (S1: 500 seeds/m2, S2: 400 seeds/m2, S3: 300 seeds/m2).
Figure 10. Heatmap of mean Normalized Difference Vegetation Index (NDVI1) by nitrogen dose (N1: 120 kg/ha, N2: 60 kg/ha, N3: 0 kg/ha) and seed rate (S1: 500 seeds/m2, S2: 400 seeds/m2, S3: 300 seeds/m2).
Plants 13 01574 g010
Figure 11. Distribution of Normalized Difference Vegetation Index (NDVI1) values by nitrogen dose levels (N1: 120 kg/ha, N2: 60 kg/ha, N3: 0 kg/ha) at the mean seed rate (S2: 400 seeds/m2). The black diamond symbol indicates outliers. Significance levels were determined by one-way ANOVA and subsequent Tukey test. Different letters indicate significant differences among the nitrogen application levels.
Figure 11. Distribution of Normalized Difference Vegetation Index (NDVI1) values by nitrogen dose levels (N1: 120 kg/ha, N2: 60 kg/ha, N3: 0 kg/ha) at the mean seed rate (S2: 400 seeds/m2). The black diamond symbol indicates outliers. Significance levels were determined by one-way ANOVA and subsequent Tukey test. Different letters indicate significant differences among the nitrogen application levels.
Plants 13 01574 g011
Figure 12. Heatmap of mean Normalized Difference Vegetation Index (NDVI2) by nitrogen dose (N1: 120 kg/ha, N2: 60 kg/ha, N3: 0 kg/ha) and seed rate (S1: 500 seeds/m2, S2: 400 seeds/m2, S3: 300 seeds/m2).
Figure 12. Heatmap of mean Normalized Difference Vegetation Index (NDVI2) by nitrogen dose (N1: 120 kg/ha, N2: 60 kg/ha, N3: 0 kg/ha) and seed rate (S1: 500 seeds/m2, S2: 400 seeds/m2, S3: 300 seeds/m2).
Plants 13 01574 g012
Figure 13. Distribution of Normalized Difference Vegetation Index (NDVI2) values by nitrogen dose levels (N1: 120 kg/ha, N2: 60 kg/ha, N3: 0 kg/ha) at the mean seed rate (S2: 400 seeds/m2). Significance levels were determined by one-way ANOVA and subsequent Tukey test. Different letters indicate significant differences among the nitrogen application levels.
Figure 13. Distribution of Normalized Difference Vegetation Index (NDVI2) values by nitrogen dose levels (N1: 120 kg/ha, N2: 60 kg/ha, N3: 0 kg/ha) at the mean seed rate (S2: 400 seeds/m2). Significance levels were determined by one-way ANOVA and subsequent Tukey test. Different letters indicate significant differences among the nitrogen application levels.
Plants 13 01574 g013
Figure 14. Heatmap of mean Normalized Difference Vegetation Index (NDVI3) by nitrogen dose (N1: 120 kg/ha, N2: 60 kg/ha, N3: 0 kg/ha) and seed rate (S1: 500 seeds/m2, S2: 400 seeds/m2, S3: 300 seeds/m2).
Figure 14. Heatmap of mean Normalized Difference Vegetation Index (NDVI3) by nitrogen dose (N1: 120 kg/ha, N2: 60 kg/ha, N3: 0 kg/ha) and seed rate (S1: 500 seeds/m2, S2: 400 seeds/m2, S3: 300 seeds/m2).
Plants 13 01574 g014
Figure 15. Pearson’s correlation matrix for wheat yield, biomass, Normalized Difference Vegetation Index (NDVI), and protein content.
Figure 15. Pearson’s correlation matrix for wheat yield, biomass, Normalized Difference Vegetation Index (NDVI), and protein content.
Plants 13 01574 g015
Figure 16. Heatmap of mean SPAD1 values by nitrogen dose (N1: 120 kg/ha, N2: 60 kg/ha, N3: 0 kg/ha) and seed rate (S1: 500 seeds/m2, S2: 400 seeds/m2, S3: 300 seeds/m2).
Figure 16. Heatmap of mean SPAD1 values by nitrogen dose (N1: 120 kg/ha, N2: 60 kg/ha, N3: 0 kg/ha) and seed rate (S1: 500 seeds/m2, S2: 400 seeds/m2, S3: 300 seeds/m2).
Plants 13 01574 g016
Figure 17. Distribution of SPAD1 values by nitrogen dose levels (N1: 120 kg/ha, N2: 60 kg/ha, N3: 0 kg/ha) at the mean seed rate (S2: 400 seeds/m2). Significance levels were determined by one-way ANOVA and subsequent Tukey test. Different letters indicate significant differences among the nitrogen application levels.
Figure 17. Distribution of SPAD1 values by nitrogen dose levels (N1: 120 kg/ha, N2: 60 kg/ha, N3: 0 kg/ha) at the mean seed rate (S2: 400 seeds/m2). Significance levels were determined by one-way ANOVA and subsequent Tukey test. Different letters indicate significant differences among the nitrogen application levels.
Plants 13 01574 g017
Figure 18. Heatmap of mean SPAD2 values by nitrogen dose (N1: 120 kg/ha, N2: 60 kg/ha, N3: 0 kg/ha) and seed rate (S1: 500 seeds/m2, S2: 400 seeds/m2, S3: 300 seeds/m2).
Figure 18. Heatmap of mean SPAD2 values by nitrogen dose (N1: 120 kg/ha, N2: 60 kg/ha, N3: 0 kg/ha) and seed rate (S1: 500 seeds/m2, S2: 400 seeds/m2, S3: 300 seeds/m2).
Plants 13 01574 g018
Figure 19. Heatmap of mean SPAD3 values by nitrogen dose (N1: 120 kg/ha, N2: 60 kg/ha, N3: 0 kg/ha) and seed rate (S1: 500 seeds/m2, S2: 400 seeds/m2, S3: 300 seeds/m2).
Figure 19. Heatmap of mean SPAD3 values by nitrogen dose (N1: 120 kg/ha, N2: 60 kg/ha, N3: 0 kg/ha) and seed rate (S1: 500 seeds/m2, S2: 400 seeds/m2, S3: 300 seeds/m2).
Plants 13 01574 g019
Figure 20. (a) Trends of NDVI values across key growth stages under different nitrogen doses (N1: 120 kg/ha, N2: 60 kg/ha, N3: 0 kg/ha). (b) Trends of SPAD values across key growth stages under different nitrogen doses (N1: 120 kg/ha, N2: 60 kg/ha, N3: 0 kg/ha).
Figure 20. (a) Trends of NDVI values across key growth stages under different nitrogen doses (N1: 120 kg/ha, N2: 60 kg/ha, N3: 0 kg/ha). (b) Trends of SPAD values across key growth stages under different nitrogen doses (N1: 120 kg/ha, N2: 60 kg/ha, N3: 0 kg/ha).
Plants 13 01574 g020
Figure 21. Heatmap of mean canopy temperature (CT1) by nitrogen dose (N1: 120 kg/ha, N2: 60 kg/ha, N3: 0 kg/ha) and seed rate (S1: 500 seeds/m2, S2: 400 seeds/m2, S3: 300 seeds/m2).
Figure 21. Heatmap of mean canopy temperature (CT1) by nitrogen dose (N1: 120 kg/ha, N2: 60 kg/ha, N3: 0 kg/ha) and seed rate (S1: 500 seeds/m2, S2: 400 seeds/m2, S3: 300 seeds/m2).
Plants 13 01574 g021
Table 1. Research station’s soil characteristics, annual rainfall, geographic details, and elevation over five growing seasons.
Table 1. Research station’s soil characteristics, annual rainfall, geographic details, and elevation over five growing seasons.
Research StationSoil TypeAnnual Rainfall (mm) per Growing Season Geographic
Coordinates
Altitude (Meters)
2016–20172017–20182018–20192019–20202020–2021Mean PrecipitationsLatitudeLongitude
Sidi El Aidi (SEA)Vertisol29050521024246734333.12218°N7.63315°W235
Table 2. Soil and Site-Specific Parameters.
Table 2. Soil and Site-Specific Parameters.
ParameterValue Classification
pH8.25Alkaline soil
Conductivity (C.E)0.45 dS/mLow salinity
Nitrate (NO3-N)12.83 ppmLow nitrogen
Phosphorus (P)13.40 ppmHigh phosphorus
Potassium (K)182.68 ppmHigh potassium
Organic Matter (MO)2.13%Moderate organic matter
Table 3. Durum wheat cultivars used in the experiment.
Table 3. Durum wheat cultivars used in the experiment.
CultivarYear of RegistrationPedigree
Karim1985Bittern ‘S’ «JO’S’. AA”:S’//FG’S’»
Nassira2003INRA Selection on CIMMYT EII, 12 TA14/BD3//Isly # CF41530–1548
Faraj2007Hybrid Nassira, Qarmal, Lahn (ICARDA)
Luiza2011RASCON_39/TILO_1
Itri2017RISSA/GAN//POHO_1/3/PLATA_3//CREX/ALLA/x Karim
Table 4. Impact of nitrogen doses and seeding rates on wheat yield, biomass, and protein content across different treatment combinations.
Table 4. Impact of nitrogen doses and seeding rates on wheat yield, biomass, and protein content across different treatment combinations.
Nitrogen DoseSeed RateMean Yield (kg/ha)Standard DeviationMean Biomass (kg/ha)Standard DeviationMean Protein %Standard Deviation
N1S12841 ab3359800 a165116.91 ab2.55
S22622 abc6149017 ab119715.26 b2.27
S32761 abc5139367 ab126817.82 ab1.45
Mean 27414879395137216.662.09
N2S12279 abc4438825 ab107516.53 ab1.69
S22978 a10359742 a274217.16 ab1.74
S32008 c4627508 b112217.79 ab0.91
Mean 24226478692164617.161.45
N3S12349 abc4537542 b121619.12 a3.24
S22203 bc3637450 b69617.86 ab3.16
S32470 abc1908233 ab102717.76 ab1.47
Mean 2341335774298018.252.62
Source F-Value (Yield)p-Value (Yield)F-Value (Biomass)p-Value (Biomass)F-Value (Protein)p-Value (Protein)
Nitrogen Dose 4.650.012 *9.920.000 ***4.110.02 *
Seed Rate 0.920.402 ns0.620.539 ns1.790.173 ns
Nitrogen × Seed Rate 4.380.003 **3.530.01 *1.880.122 ns
Notes: ANOVA: * = p ≤ 0.05, ** = p ≤ 0.01, *** = p ≤ 0.001, and ns = not significant. Means with identical letters are not significantly different at the 95% confidence interval (Tukey method).
Table 5. Effects of nitrogen doses and seeding rates on Normalized Difference Vegetation Index (NDVI) measurements at three key growth stages of wheat.
Table 5. Effects of nitrogen doses and seeding rates on Normalized Difference Vegetation Index (NDVI) measurements at three key growth stages of wheat.
Nitrogen DoseSeed RateNDVI 1 MeanStandard DeviationNDVI 2 MeanStandard DeviationNDVI 3 MeanStandard Deviation
N1S10.834 a0.031 0.818 ab0.0450.813 a0.034
S20.817 a0.038 0.807 ab0.0440.81 a0.044
S30.84 a0.028 0.857 a0.0260.781 a0.074
Mean 0.8300.0320.8270.0380.8010.051
N2S10.815 a0.022 0.796 ab0.0290.767 a0.06
S20.818 a0.06 0.802 ab0.0670.784 a0.059
S30.805 ab0.028 0.793 b0.0370.8 a0.027
Mean 0.8130.0370.7970.0440.7840.049
N3S10.757 b0.04 0.772 b0.0210.761 a0.036
S20.806 ab0.023 0.788 b0.0360.772 a0.037
S30.799 ab0.053 0.777 b0.0660.779 a0.039
Mean 0.7870.0390.7790.0410.7710.037
Source F-Value (NDVI 1)p-Value (NDVI 1)F-Value (NDVI 2)p-Value (NDVI 2)F-Value (NDVI 3)p-Value (NDVI 3)
Nitrogen Dose 9.640.000 ***9.420.000 ***3.090.051 ns
Seed Rate 1.030.361 ns0.760.472 ns0.240.785 ns
Nitrogen × Seed Rate 2.560.045 *1.610.181 ns1.30.278 ns
Notes: ANOVA: * = p ≤ 0.05, *** = p ≤ 0.001, and ns = not significant. Means with identical letters are not significantly different at the 95% confidence interval (Tukey method).
Table 6. Influence of nitrogen dose and seed rate on SPAD chlorophyll content at key developmental stages of wheat.
Table 6. Influence of nitrogen dose and seed rate on SPAD chlorophyll content at key developmental stages of wheat.
Nitrogen DoseSeed RateSPAD 1 MeanStandard DeviationSPAD 2 MeanStandard DeviationSPAD 3 MeanStandard Deviation
N1S150.03 ab1.7849.88 a2.3145.98 ab4.76
S251.97 ab2.2550.62 a1.6348.98 a4.92
S352.47 a2.8150.25 a2.546.23 ab6.57
Mean 51.492.2850.252.1547.065.42
N2S149.73 ab2.9350.69 a2.9446.79 ab3.12
S251.08 ab2.0248.92 a2.4242.11 b4.42
S352.15 ab2.3249.21 a2.0947.75 ab5.19
Mean 50.992.4249.612.4845.554.24
N3S150.9 ab2.1250.79 a2.7948.02 ab3.66
S248.81 b3.0449.05 a2.9645.12 ab3.7
S350.12 ab2.0249.51 a3.0841.7 b4.72
Mean 49.942.3949.782.9444.954.03
Source F-Value (SPAD 1)p-Value (SPAD 1)F-Value (SPAD 2)p-Value (SPAD 2)F-Value (SPAD 3)p-Value (SPAD 3)
Nitrogen Dose 3.220.045 *0.50.608 ns1.650.198 ns
Seed Rate 2.530.086 ns1.140.324 ns1.210.304 ns
Nitrogen × Seed Rate 2.40.056 ns0.840.507 ns4.430.003 **
Notes: ANOVA: * = p ≤ 0.05, ** = p ≤ 0.01 and ns = not significant. Means with identical letters are not significantly different at the 95% confidence interval (Tukey method).
Table 7. Effects of nitrogen dose and seed rate on canopy temperature dynamics during wheat growth stages.
Table 7. Effects of nitrogen dose and seed rate on canopy temperature dynamics during wheat growth stages.
Nitrogen DoseSeed RateCT 1 MeanStandard DeviationCT 2 MeanStandard DeviationCT 3 MeanStandard Deviation
N1S120.95 a0.6024.05 a1.0823.58 a1.64
S220.72 ab0.2223.64 a0.8423.63 a0.80
S319.72 c0.8522.92 a1.4223.60 a1.01
Mean 20.460.5623.531.1123.601.15
N2S120.26 abc0.6123.16 a0.6223.93 a0.85
S220.00 bc0.4523.61 a0.5723.59 a1.05
S320.48 abc0.5523.49 a0.8524.11 a0.93
Mean 20.240.5423.420.6823.880.94
N3S120.24 abc0.4223.07 a0.6723.78 a0.89
S220.03 bc0.3823.04 a0.7024.00 a0.81
S320.54 ab0.6223.06 a0.6623.90 a0.87
Mean 20.270.4723.050.6723.890.86
Source F-Value (CT 1)p-Value (CT 1)F-Value (CT 2)p-Value (CT 2)F-Value (CT 3)p-Value (CT 3)
Nitrogen Dose 1.450.24 ns2.540.085 ns0.780.464 ns
Seed Rate 1.840.166 ns1.020.365 ns0.140.869 ns
Nitrogen × Seed Rate 8.260.000 ***2.070.093 ns0.330.856 ns
Notes: ANOVA: *** = p ≤ 0.001, and ns = not significant. Means with identical letters are not significantly different at the 95% confidence interval (Tukey method).
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Hnizil, O.; Baidani, A.; Khlila, I.; Nsarellah, N.; Laamari, A.; Amamou, A. Integrating NDVI, SPAD, and Canopy Temperature for Strategic Nitrogen and Seeding Rate Management to Enhance Yield, Quality, and Sustainability in Wheat Cultivation. Plants 2024, 13, 1574. https://doi.org/10.3390/plants13111574

AMA Style

Hnizil O, Baidani A, Khlila I, Nsarellah N, Laamari A, Amamou A. Integrating NDVI, SPAD, and Canopy Temperature for Strategic Nitrogen and Seeding Rate Management to Enhance Yield, Quality, and Sustainability in Wheat Cultivation. Plants. 2024; 13(11):1574. https://doi.org/10.3390/plants13111574

Chicago/Turabian Style

Hnizil, Oussama, Aziz Baidani, Ilham Khlila, Nasserelhaq Nsarellah, Abdelali Laamari, and Ali Amamou. 2024. "Integrating NDVI, SPAD, and Canopy Temperature for Strategic Nitrogen and Seeding Rate Management to Enhance Yield, Quality, and Sustainability in Wheat Cultivation" Plants 13, no. 11: 1574. https://doi.org/10.3390/plants13111574

APA Style

Hnizil, O., Baidani, A., Khlila, I., Nsarellah, N., Laamari, A., & Amamou, A. (2024). Integrating NDVI, SPAD, and Canopy Temperature for Strategic Nitrogen and Seeding Rate Management to Enhance Yield, Quality, and Sustainability in Wheat Cultivation. Plants, 13(11), 1574. https://doi.org/10.3390/plants13111574

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