Spectral Response of Camelina ( Camelina sativa (L.) Cran�) to Diﬀerent Nitrogen Fertilization Regimes under Mediterranean Conditions

: Knowledge about the spectral response of camelina under diﬀerent regimes of nitrogen (N) fertilization is very scarce. Therefore, 2-year open-ﬁeld trials were carried out in the 2021 and 2022 growing seasons with the aim of evaluating the spectral response of spring camelina to four diﬀerent N fertilization regimes by using remote (UAV) and proximal (leaf-clip Dualex) sensing techniques. The tested treatments were: (i) control: no N application (T0); (ii) top dressing: 60 kg N ha −1 before stem elongation (T1); basal dressing: 60 kg N ha −1 at sowing (T2); basal + top dressing combination: 60 kg N ha −1 at sowing + 60 kg N ha −1 before stem elongation (T3). Camelina seed yield and N use eﬃciency were strongly aﬀected by fertilization regimes, with the best results obtained at T2. A reduction in plant development and seed yield was detected in 2022, probably due to the rise in air temperatures. A signiﬁcant eﬀect of both growing season and N fertilization was observed on the photosynthetic pigments content with the T1 highest values in 2022. The highest seed oil content was achieved at T1, while the protein content increased with increasing N, with the best values at T3. Positive and signiﬁcant correlations were observed among several vegetation indices obtained through UAV ﬂights (NDVI, MRS705, FGCC) and seed yield, as well as between FGCC and leaf N concentration. Overall, these ﬁndings demonstrate the feasibility of utilizing remote sensing techniques from UAVs for predicting seed yield in camelina.


Introduction
Nowadays, the great potential arising from technological progress seems to be able to meet the challenges of modern agriculture aimed at increasing the sustainability and quality of cropping systems [1]. The application of new monitoring techniques to manage intensive agriculture is crucial. Indeed, modern technologies allow farmers to keep plants healthy and achieve high and stable yields, and, at the same time, reduce environmental impacts. Particular a ention is generally paid to nitrogen (N) fertilization, as N is an essential mineral nutrient known to enhance photosynthetic capacity and improve crop yield [2,3]. In agriculture, nitrogen use efficiency (NUE) is generally quite low, and nitrogen losses, mainly through leaching (nitrates) and gas emissions (nitrous oxide and ammonia), cause serious environmental concerns such as water contamination and evaluate the suitability of remote and proximal sensing in monitoring nutritional status. In particular, the effects of different N fertilizer regimes (rate and timing) on crop growth (seed yield and yield components) and quality (protein and oil seed content) were assessed, along with traditional (destructive diagnosis) and non-destructive methods (Dualex portable sensor and UAVs equipped with RGB and high-definition multispectral cameras).

Experimental Setup and Site Description
Two-year experimental field trials were conducted during 2021 and 2022 growing seasons at the Experimental Center of the Department of Agriculture, Food, and Environment (DAFE) of the University of Pisa, San Piero a Grado, Pisa, Italy (43°40′ N; 10°19′ E; 5 m above sea level). In both years of cultivation, winter wheat (Triticum turgidum L. subsp. durum (Desf.) Husn.) preceded camelina. An integrated management system was adopted with conventional tillage practices and mineral fertilization. Pre-planting phosphorus and potassium fertilizations were performed at a rate of 80 kg P2O5 ha −1 as triple superphosphate and 50 kg K2O ha −1 as potassium sulfate. Both years, spring sowing of camelina was performed by a Wintersteiger plot drill on 0.15 m spaced rows and a depth of 0.01 m on 8 March. The seeding rate, in both years, was around 6.5 kg ha −1 , considering percent seed germination as well as 1000-seed weight (TSW), in order to reach a target of 500 plants m −2 . For these experiments, the commercial variety Calena was used, and different nitrogen levels (applied as ammonium nitrate, NH4NO3) and timing (basal and top dressing) were compared. The tested treatments were: T0 = control: no N application; T1 = top dressing: 60 kg N ha −1 before stem elongation; T2 = basal dressing: 60 kg N ha −1 at sowing; T3 = basal + top dressing combination: 60 kg N ha −1 at sowing + 60 kg N ha −1 before stem elongation. A completely randomized block design was adopted with three replications for each treatment (plot size: 8 m × 6 m). No chemical treatments for pests and diseases were necessary during the whole experiment, and weeds were controlled by hand weeding.
Daily meteorological data (temperature and rainfall) were measured through an automatic meteorological station located near the experimental site. The climate in the area is Mediterranean (Csa according to Köppen classification), characterized by hot and dry summers, with rainfall concentrated during winter. Long-term (20 years) rainfall in the March-June period was 212 mm, while in the same period, temperature varied from a minimum of 5.8 °C (1st decade of March) to a maximum of 26.6 °C (3rd decade of June).
Soils on which trials were carried out are classified as Typic Xerofluvens according to USDA soil taxonomy. In both growing seasons, soil samples were collected at 0-30 cm depth at the beginning of the trial (spring 2021 and spring 2022) for physical and chemical analysis. According to USDA soil texture classes, in both years (2021 and 2022), field trials were characterized by loamy soils with a sub-alkaline pH, good soil organic ma er (SOM) content, medium-low levels of total nitrogen content and available phosphorus, and medium-high levels of exchangeable potassium (Table 1). Soil pH determination was performed on a 1:2.5 soil:water suspension following McLean procedure [29]. Electrical conductivity (EC) was measured at 20 °C by using a GLP-31 Crison conductimeter (52.93 electrode) (Montepaone s.r.l., San Mauro Torinese, Torino, Italy). Soil total nitrogen was evaluated using the macro-Kjeldahl digestion procedure [30], available phosphorus by colorimetric analysis using the Olsen method [31], and the cation exchange capacity (CEC) following Mehlich method [32]. The exchangeable K was determined using the Thomas method [33]. Soil organic ma er was calculated by multiplying by 1.724 the soil organic carbon concentration, measured using the modified Walkley-Black wet combustion method [34].

Crop Growth Cycle Monitoring
For each growing season, the crop growth cycle was monitored, and the main phenological phases of camelina, including basal rose e (BBCH 19), beginning of stem elongation (BBCH 31), full flowering (BBCH 65), and seed maturation (BBCH 89), were determined using the extended BBCH scale described by Martinelli and Galasso [35]. In addition, number of days and accumulated growing degree days (GDD) for both years of cultivation (2021 and 2022) were calculated using daily maximum air temperature (Tmax), daily minimum air temperature (Tmin), and base temperature (Tbase) of 4 °C [36], as follows: (1)

Non-Destructive Analysis and Destructive Sampling
Proximal and remote sensing methods, in combination with destructive methods, were carried out when camelina plants were at full flowering stage (BBCH 65) (21 and 12 May, in 2021 and 2022, respectively). In particular, aerial surveys were carried out on two sites of camelina cultivation using the DJI MAVIC 2 PRO drone equipped with a Hasselblad L1D-20C camera and a PARROT SEQUOIA multispectral camera (Parrot Drones S.A.S. , Paris, France) with four spectral bands (Green, Red, Red-edge, NIR). The flights were performed at an altitude of 15 m, with an overlap of 80%, and a spatial resolution in terms of pixels on the ground equal to 3.5 mm for the RGB camera and 14 mm for the multispectral camera. The multispectral sensor was calibrated using a calibration panel (provided by the sensor producer), which must be scanned after takeoff and before landing to correct the acquired data based on the incident light conditions during the flight. The survey activities were preceded by the identification on the ground, with the use of the total station, of 13 ground control points (GCPs), with known coordinates, useful for the georeferencing of all the surveys in a common reference system (ETRF2000 UTM 32N). These points were then materialized with clearly visible "targets" from drones, useful for the georeferencing of the models and orthophotos obtained from the processing of aerial images. Data obtained from UAV flights were used for the calculation of some vegetation indices (VIs) reported in Table 2.
After aerial surveys, three leaves randomly chosen from three plants for each plot were clipped with a field-portable leaf-clip sensor Dualex ® Force-A (Orsay, Cedex, France) for the evaluation of chlorophyll content (µg cm −2 ), flavonols (Abs unit), anthocyanins (Abs unit), and nitrogen balance index (NBI). This instrument performs simple, fast, and non-destructive measurements acquiring information from the UV absorbance measurement of the leaf epidermis by double excitation of chlorophyll fluorescence [37]. Two measurements for each leaf per plant and per plot were performed on both the adaxial and abaxial sides, with 72 total measurements detected. Chlorophyll and anthocyanin values were obtained as average between adaxial and abaxial sides, while the flavonol value was the sum, and NBI was the ratio between chlorophyll and flavonol values. Table 2. Vegetation indices (VIs) evaluated in this study.

Vegetation Indices
Formula Normalized difference vegetation index (NDVI) [38] NDVI = (NIR − RED)/(NIR + RED) Ratio vegetation index (RVI) [39] RVI = RED/NIR Transformed normalized difference vegetation index [40] TNDVI = (NIR − RED)/(NIR + RED) + 0.5] 1/2 Renormalize difference vegetation index (RDVI) [41] RDVI = (NIR − RED)/(NIR + RED) 1/2 Improved modified chlorophyll absorption ratio index (MCARI2) [42] MCARI 2 = 1. In order to validate the collected data non-destructively, three fresh leaf discs (Ø 0.8 cm) from each leaf per three plants per plot were sampled immediately after Dualex ® Force-A measurements and processed for chlorophyll and carotenoid content according to Lichtentahler and Buschmann [45]. The extraction was carried out in 4 mL of 80% acetone by placing the samples in a cold chamber to reduce the volatility of acetone and in the dark so as not to interfere with chlorophyll activity. Finally, destructive plant samplings (~0.075 m 2 ) were performed for N determination in the different plant organs (stems, leaves, and inflorescences). Fresh weight was measured, and plants subsequently allowed to dry in a ventilated oven (70 °C) for dry weight determination and evaluated for their moisture content. Soil samples were collected in the same area of destructive plant samplings to analyze soil moisture and nitrate concentration (NO3 − ). The soil moisture was calculated with the gravimetric method (samples were weighed twice before and after oven drying at 60 °C until constant weight). NO3 − content was determined on the same samples taken to monitor soil moisture by pooling the first three (0-30 cm) and the last two (30-50 cm) layers. The ion chromatography method was used (Dx-500 ion chromatograph; Dionex, Sunnyvale, CA, USA).

Image Processing
The RGB and multispectral images of experimental field plots acquired from the UAVs were orthorectified and mosaicked using Agisoft Photoscan Professional Edition 1.1.6 (Agisoft LLC, St. Petersburg, Russia) in order to correct the possible perspective deformations of the single frames. Finally, subsets of the orthomosaics were created by isolating the various test plots using Erdas Imagine software by Hexagon AB. The RGB images of the textures were then processed using Canopeo, which analyzes the images using color values in the red (R), green (G), and blue (B) system and classifies all the pixels according to three selection parameters: R/G, B/G, and the green excess index (GEI = 2G-R-B) [46][47][48][49]. The primary image was transformed into a binary image depending on whether or not the pixels met the selection criteria (R/G < 0.95, B/G < 0.95, and GEI > 20, default values). The final result of the image processing is, therefore, a percentage index of FGCC (fraction green canopy cover), which expresses the number of pixels complying with the aforementioned parameters and which are therefore classified as green, with respect to the number of total pixels. The multispectral images were processed using the Erdas Imagine 2020 software, a tool for geospatial digital image processing developed by Hexagon AB [50].

Crop Production and Yield Component
Camelina plants were harvested by manually sampling 7 rows from the central portion of each plot for a length of 1.5 m (sampling area = 1.58 m 2 ) at full seed ripening, when more than 90% of siliques were dried and turned brown, and most seeds were reddish-brown in color (seed moisture ≤ 12%, BBCH 89). Spring-sown camelina were harvested on 24 and 20 June, in 2021 and 2022, respectively. Plant density, plant height, branching, aboveground biomass dry weight (stems, empty siliques, and seeds), seed yield, and yield components (number of siliques, number of seeds per silique, 1000-seeds weight) were evaluated. To assess seed yield, the plants were threshed by a fixed machine, using sieves suitable for small seeds. Fresh weight was measured, and plants subsequently allowed to dry in a ventilated oven (70 °C) for dry weight determination and evaluated for their moisture content. The harvest index (HI), an important trait associated with crop yields and the successful partitioning of photosynthates to harvestable products, was calculated by dividing dry seed weight by the dry weight of total aboveground biomass.
Thousand seed weight (TSW, g DW) was also assessed at the Seed Research and Testing Laboratory of DAFE on representative seed samples deriving from each plot according to the International Rules for Seed Testing [51]. The content of N (determined by macro-Kjeldahl digestion procedure) in the different plant organs (stems, straw, and seeds) was also evaluated. Fresh weight was measured, and plant organs subsequently allowed to dry in a ventilated oven (70 °C) for dry weight determination and evaluated for their moisture content. Moreover, Agronomic Efficiency of applied Nitrogen (AE, kg seed yield per kg N applied) was calculated as the difference between the seed yield of treatments with N and the achene yield of the non-fertilized control, divided by the N rate applied. Apparent recovery efficiency (ARE) is defined as the difference in N accumulation between plots receiving fertilizer and unfertilized plots and is in proportion to the amount of N fertilizer applied. AE and ARE indices were then calculated with the following equations [52,53]: (2) where N uptake was computed by multiplying biomass yield per N tissue concentration.

Seed Quality
Seed oil and protein contents were determined on representative samples from each plot. Oil content was determined using a Soxhlet extractor apparatus (mod. R 306) from Behr Labor-Technik (Düsseldorf, Germany). Briefly, about 30 g of camelina seeds were finely ground in a coffee grinder for 40 s. An aliquot of 1.5 g of ground material was exactly weighed in a cellulose extraction thimble (22 mm × 80 mm) from Axiva Sichem Biotech (Delhi, India). The thimble was successively inserted in a glass extractor, and oil extraction was carried out for two hours in a Soxhlet extractor using 60 mL of n-hexane as an organic solvent. The extract containing the oily fraction was then filtered with anhydrous sodium sulfate in a 100 mL flat bo om flask and removed under reduced pressure at 30 °C in a rotary evaporator. The residual oil was dried under a gentle nitrogen flow for 5 min with the flask in a water bath (50-55 °C), exactly weighed, transferred by means of 5 mL of nhexane/i-propanol 4/1 (v/v) in a 10 mL Teflon screw-cap glass tube, and stored at −18 °C. Solvents used were of analytical grade. Oil yield (kg DM ha −1 ) was obtained by multiplying seed yield by seed oil content of each individual replicate. Seed crude protein content was calculated by multiplying the total nitrogen percentage by 6.25, and total N content was determined by means of the mini-Kjeldahl method on different plant organs (stems, empty siliques, and seeds).

Statistical Analysis
Data were subjected to statistical analysis using the software COSTAT cohort V6.201 (2002). Before ANOVA, Levene's test was used to verify homoscedasticity, and Shapiro-Wilk test was used to check the normal distribution of residuals. A two-way ANOVA for evaluating the effects of growing season (2021 vs. 2022) and the different N fertilization regimes (control, T0, T1, T2, and T3) on agronomic and qualitative parameters were performed. Means were compared using Tukey's HSD test when ANOVA F-test was significant at p ≤ 0.05. For each year, Pearson's correlation coefficients were computed between VIs and Dualex chlorophyll and seed yield, leaf chlorophyll, and leaf N concentration.

Meteorological Data and Crop Phenology
Weather conditions recorded over the study timespan are reported in Figure 1. The average minimum and maximum temperatures referring to the camelina growing period (from March to June) were very similar for both years, with a minimum of 8.7 °C and 9.6 °C, and a maximum of 21.4 °C and 21.8 °C in spring 2021 and 2022, respectively. In the first growing season, the increase in temperature occurred in the third decade week of April 2021, while in the second growing season, temperatures started to rise significantly from the second decade of April 2022. A shortening of about 10 days in the vegetative growth phase of the second year of experimentation was induced by the differences in the trend of average temperatures (Table 3). In the first growing season (spring 2021), cumulative rainfall and the GDD from sowing to seed maturity were equal to 117.9 mm and 1195.2, respectively, while in the second growing season (2022), they were 140.8 mm and 1241.9 GDD, respectively. A few days of cycle length discern the two spring sowings; in fact, the plants completed their growth cycle in 109 and 105 days in 2021 and 2022, respectively (Table 3).

Vegetation Indices (VIs)
Results of the statistical two-way analysis (F-Test, ANOVA), related to the effects of year (Y) and N fertilization regimes (F) and their interaction (Y × F), are reported in Table  4. The effect of the year (Y) was significant for the indices NDVI, RVI, TNDVI, WvVI, MSR705, SAVI, and FGGC calculated on the images acquired from the drone. While the effect of N fertilization (F) was statistically appreciable only for NDVI, MSR705, and FGGC, the interaction between year and fertilization was significant only for NDVI.

Dualex Detection and Destructive Samplings
The effects of year (Y) and N fertilization (F) and their reciprocal interactions were evaluated on Dualex parameters and on destructive samplings. In particular, flavonols and NBI values detected by the Dualex leaf-cleap sensor were significantly affected by the year, and only flavonols were also affected by the Y x F interaction (Table 5). Specifically, the flavonols showed an average value of the different N fertilizers higher than T0 in the first year (+17.4%). T0, in the second year, instead, showed the highest value (+21.34%) compared to the average of all other N fertilizers, except for T1. Furthermore, the NBI detected in 2021 showed an average value higher (+12.4%) than in 2022.
Regarding destructive samplings, a significant effect of both Y and F and their reciprocal interaction was also observed in the content of chlorophylls and carotenoids deriving from spectrophotometric assays with some exceptions, as reported in Table 5. Specifically, chlorophyll a content was higher at the top dressing in the second year: T2 (+36.2%) and T3 (+32.3%). T0 and T2 in 2022, despite having lower values, were, on average, higher (+39%) than all fertilizers in 2021. The concentration of chlorophyll b, on the other hand, was, on average, higher in the second year (+41%), with T1 and T2 exhibiting higher average values (+12.9%) than T0 and T3. The concentration of carotenoids was, on average, higher in T0 and T2 in 2022 (+13.9%) than in T1 and T3 and compared to the average in 2021 (+69%).
Soil moisture and nitrate concentration were only affected by the year, while total aboveground biomass produced by the crop was significantly affected by both Y and F, as well as by their reciprocal interaction (Table 5). In particular, soil moisture percentage was higher in 2022 (+17%) compared to 2021, and nitrate concentration was also higher in the second year of the experimentation (+42.5%) than in the first one.

Dry Ma er Accumulation and Nitrogen Uptake
At full flowering, the total aboveground biomass was significantly higher in 2021 (+43%) than in 2022 (Figure 2A). Considering the dry ma er partitioned in the different organs, camelina plants produced more stems and inflorescences in 2021 than in 2022, except for leaf biomass production. In fact, T2 (2021) and T3 (2022) produced more leaves than the other treatments. Among these, all N fertilizer treatments were more productive than the control in two years (Figure 2A). At seed maturity, camelina plants produced more total aboveground biomass at harvest when they were treated with top-dressing N fertilizer (T1) and with basal + top dressing combination (T3) in 2021 (+41%) and 2022 (+23.5%) compared to the other respective treatments. Nevertheless, 2021 showed clearly higher average values of biomass produced than those of 2022 (+45.5%) ( Figure 2B).
Total N uptake by camelina plants when they were at full flowering, was higher in 2021 (+33%), with a be er distribution of N uptake among organs (leaves, stems, and inflorescences) than in 2022 ( Figure 2C). At seed maturity, in the first year of cultivation, camelina plants treated with N fertilization showed an average N straw uptake higher (+30%) than in the second one, although the control showed an opposite trend. Additionally, N seed uptake was also higher in 2021 (+37.4%) compared to 2022 if considering the different N fertilizer levels ( Figure 2D).
Agronomic efficiency (AE), described as the economic production obtained per unit of nitrogen applied, was significantly highest in 2021 and when the plants were treated with N-fertilizer at sowing (T2) followed by top-dressing N-fertilizer (T1) ( Figure 2E). Furthermore, the apparent recovery efficiency (ARE), commonly defined as the difference in nutrient uptake in aboveground parts of the plant between the fertilized and unfertilized crop relative to the quantity of nutrients applied, was significantly higher in 2021 when plants were treated with N-fertilizer both at sowing and top-dressing ( Figure  2F). , and on agronomic efficiency (AE) (E) and apparent recovery efficiency (ARE) (F). Means ± standard error followed by identical le ers are not significantly different for p ≤ 0.05, according to Tukey's HSD post-hoc test. Significance is indicated as follows: *, p ≤ 0.05; **, significant at p ≤ 0.01; ***, significant at p ≤ 0.001. Red line in sub-image F: values >1 indicate that N up-taken by crop was higher than that applied as fertilizer.

Yield, Yield Components, and Seed Quality
The effects of year (Y) and N fertilization (F) and their reciprocal interaction were evaluated on different agronomic and seed qualitative characteristics, as reported in Table  6. The results obtained showed that plant height was significantly influenced by both Y and F, with the highest values recorded in camelina plants in the first year of experimentation, with an average of 13 cm taller compared to the second year. No differences in plant height were observed between the T0 and T1 treatments, although both were smaller (−9.2 cm) than the T2 and T3. Furthermore, a significant effect of the year (Y) was highlighted, while N fertilization was significant only for plant density, branching, and number of siliques per plant. For plant density, a statistically significant interaction between the year and N fertilization was found. In addition, except for plant density, a significant decrease in branching (−0.8 stems per plant), siliques per plant (−11.1 siliques), seeds per silique (−0.8 seeds), and TSW (−3.7%) were observed in the second growing season (2022). In N fertilized plots, the application of 60 + 60 kg N ha −1 (T3) increased branching on average of +16.6%, +33.6%, and +48.5% compared to T2, T0, and T1, respectively, while no statistical differences with the application of 60 kg N ha −1 (T2) and control (T0) was observed. In the first year, only T2 had a higher plant density compared to the untreated control, while in the second year, no significant differences emerged between N fertilizer treatments and control ( Figure 3A). Considering the seed yield, a significant effect of both year of cultivation (Y) and N fertilizer treatment (F) and their reciprocal interaction was observed (Table 6). Camelina crop produced a higher seed yield in the first year (+38%) with respect to the second one. All the fertilizer treatments showed higher seed yield (1.9 Mg ha −1 as the average value) compared to the control (1.4 Mg ha −1 ) ( Table 6). By the Y x F interaction, the N fertilizer showed, only in 2021, the highest seed yield (2.4 Mg ha −1 ) compared to the control and all the treatments in 2022 ( Figure 3B).
No statistical difference between the two years of cultivation was observed in the oil content, although it was affected by N fertilization, with the highest value achieved when camelina plants were top-dressing fertilized (T1), as reported in Table 6 and Figure 3C. Regarding oil yield, a significant effect of both Y and F, as well as of their interaction (Y × F), was observed (Table 6). In particular, camelina produced a higher oil yield in 2021 (+38%) compared to 2022, with an average value of the different N fertilizer treatments higher (+29%) than the control ( Figure 3D). Finally, a significant effect of both Y and F was observed on the crude protein content (Table 6). Specifically, a slight increase (+6%) was observed in 2022 compared to 2021, with the highest value obtained when camelina plants were treated with basal + top dressing N fertilizer (T3) (25.3%), followed by T1 (24.8%) and T2 (24.4%). Table 6. Main effects of the year (Y), N fertilization (F), and their interaction on camelina biometric characteristics, seed yield, yield components, and qualitative traits. Values ± standard error followed by identical le ers are not significantly different for p < 0.05, according to Tukey's HSD post-hoc test. The significance of variability factors according to the F-test: ns, not significant; *, significant at p ≤ 0.05; **, significant at p ≤ 0.01; ***, significant at p ≤ 0.001 level.

Correlation Analysis
The vegetation indices (VIs) calculated on the spectral data acquired from the drone and chlorophyll obtained with Dualex leaf-clip sensor were correlated with the seed yield, leaf laboratory chlorophyll content, and leaf N concentration (Table 7). During the first year of experimentation (2021), a significantly high level of positive correlation was obtained by NDVI, TNDVI, MSR705, and FGCC indices (0.81, 0.73, 0.72, and 0.74, respectively), followed by WvVi (r = 0.67) and SAVI (0.64), while only RVI showed negative correlation (−0.69). The MCARI2 index, on the other hand, was not significantly correlated to the seed yield. In the second year of experimentation (2022), however, only MSR705 (r = 0.80) and FGCC (r = 0.82) indices had a good level of correlation, followed by NDVI (r = 0.53). The other indices were not significantly correlated with camelina seed yield. The only indices significantly correlated with seed yield in both years of experimental trials (Figure 4) were NDVI (Figure 4a), MSR705 (Figure 4b), and FGCC (Figure 4c), with the la er presenting higher significance values on average. The correlation analysis between the same vegetation indices previously mentioned and the chlorophyll leaf content obtained by the spectrophotometric method did not reveal significantly good correlations for any of the parameters analyzed (Table 7). Low correlation levels in the first year of camelina cultivation (2021) were recorded on NDVI and MCARI2, which, however, were not statistically acceptable (r < 0.52). In the first growing season, leaf N concentration was significantly and positively correlated with NDVI (r = 0.68) and FGCC (r = 0.72), while in the second growing season, a strong positive correlation was found with MSR705 (r = 0.90) and FGCC (r = 0.83). However, only FGCC was found to have a significant correlation with foliar N at full flowering in both years of the trials (Figure 4d). A significant correlation was highlighted by Dualex with seed yield in the first year of cultivation (2021) (r = 0.50) and with leaf chlorophyll content in 2021 (r = 0.51). A non-significant correlation was confirmed, instead, in 2022. During the first growing season, a significant positive correlation was found between leaf N concentration and Dualex chlorophyll; however, this correlation was not found during the second growing season. Table 7. Pearson's correlation coefficient (r) between VIs, obtained by UAV images, and chlorophyll obtained by Dualex leaf-clip sensor, in relation to camelina seed yield, chlorophyll leaf content, and foliar N % in the two growing seasons (2021 and 2022).

Discussion
The use of remote sensing has proved to be a valid tool for estimating the seed production of camelina. In fact, the NDVI, MSR 705 indices, and FGCC% showed good levels of correlation with the recorded seed yields. In line with what has been reported in other works [54][55][56], the NDVI was significantly correlated with the grain yield in a linear positive way, with a Pearson coefficient equal to 0.81 in the first year of the trial and 0.53 in the second. While the MSR705 index, although not significantly correlated with the leaf chlorophyll content, as reported in other works [43,57,58], is instead significantly correlated with the grain yield of the crop in 2021 (r = 0.72), more than in 2022 (0.80). The

Parameters
Growing Season  [59,60], but there are still few works that use this data to predict grain production [61]. In this trial, a good correlation was found between FGCC and seed yield, in agreement with what was reported by other authors [61,62]. Considering the low cost necessary to calculate this data, from an application point of view, it would seem to be the best among those proven to be reliable. In our study, FGCC was the only VIs capable of providing an indication of camelina nutritional status since it was correlated with leaf N concentration. In the literature, it is reported that the use of FGCC for yield prediction can be obtained with greater precision during the reproductive phase [63], in line with our results. However, for agronomic purposes, the most useful phase for the determination of nutritional stress is during vegetative development in order to adequately calibrate topdressing fertilization. In this regard, further studies are needed to verify whether the identified index (FGCC) can give reliable indications of camelina nutritional status even during early vegetative phases. In our trials, meteorological differences among growing seasons were not remarkable; however, in the second year, the increase in air temperatures during spring occurred almost 10 days earlier compared to 2021. It is likely that the rise in temperatures caused the reduction in vegetative development length, with negative effects on the accumulation of vegetative biomass. Consequently, the decrease in photosynthesizing biomass probably contributed to the reduction in seed yield observed in the second growing season. Anyhow, the overall duration of plant development and GDD accumulation are comparable to those already reported in the literature for spring-sown camelina in Central Italy [26]. In the first-year trial, all yield components of camelina gave be er results, and this was probably partially due to the rainfall that occurred during the reproductive phase. Therefore, water availability in this phase could have contributed to reducing potential water deficit and improving yield component performances. It has been observed that the occurrence of post-anthesis water stress conditions can lead to negative effects on yield components in various crops, such as rapeseed [64], chickpea [65], and wheat [66]. Secchi et al. [67] reported that in Brassica napus L., water stress during pod formation has a higher negative effect on the number and weight of seeds. In camelina, it has been observed that the absence of water availability or irrigation can induce a decrease in the number of siliques and branching per plant [68]. Similarly, in our trials, the number of siliques per branch showed a higher drop in the second year in comparison with the other yield components. Although there may have been an effect of early temperature rise and post-anthesis water stress during the second growing season, it is not possible to establish to what extent each of these factors influenced seed yield or individual yield components. These two stress events occurred in different periods of plant development, and, to the best of our knowledge, no research describes the sequence and timing of yield component formation in camelina. In this regard, further investigations on the production physiology of camelina are needed. In both years of our study, the increase in N dose from 60 to 120 kg N ha −1 did not lead to a significant increase in the seed yield, although it had a positive effect on the branching level. However, it is reported that camelina can be responsive to N inputs up to 300 kg N ha −1 [69]. Probably, in our study, when the environmental conditions were less favorable, as in the second growing season, the effect of nitrogen fertilization on total aboveground biomass, seed yield, and yield components was deleted. Coherently to Jankowsky et al. [70], in our study, the agronomic efficiency of fertilization in camelina varies between the two years, dropping dramatically during the second growing season. Bronson et al. [71] reported that in camelina, fertilizer agronomic efficiency can be influenced by soil water level, so the AE differences observed in our study between the two growing seasons could be a ributed to different rainfall pa erns. In addition, AE was not influenced by the N application rate, differently from that observed by Jankowsky et al. [70], who reported that fertilization agronomic efficiency started to decrease from the distribution of 80 kg N ha −1 . Similar to our findings, Allen et al. [72], in a 5-year trial, reported that in the Brassicaceae taxon, ARE (named nitrogen recovery index) varied among growing seasons, and for camelina, an average value of 0.3 was detected. Similarly, Afchar et al. [73] observed an increase in the apparent efficiency of nitrogen use in the growing seasons characterized by more favorable conditions for rainfed camelina growth. Contrarywise, in Mahli et al. [74], environmental factors (site-year) did not affect ARE (named percent recovery of applied N). During the first year, high values of ARE were reached, even exceeding the threshold value of 1.0, indicating a condition in which N removals exceeded the N distribution. On the contrary, in the second year, environmental conditions reduced plant development, at the same time decreasing N uptake deriving from fertilization in all fertilized plots, and probably leaving some unused N in the soil, likely to be lost in the environment via leaching. A reduction in root uptake during the second year is confirmed by soil analysis carried out at flowering, which showed an increase in the presence of N and moisture observed compared to the first growing season. This further underlines the necessity of identifying new methodologies to assess the crop nutritional status, also in terms of early seed yield prediction, to properly calibrate top-dressing fertilization. Finally, the oil content measured in our trials resulted in slightly lower values compared to those measured for spring-sown camelina in Central Italy [26], but nonetheless, in line with values (range 28 to 49% DM) obtained from camelina spring cultivation in Central Europe [75]. As expected, while the highest seed oil content was achieved at T1, the protein content increased with increasing N with the highest values at T3. With regard to the protein content (%), even if in previous studies values up to 32% [21] have been found, our results showed an average value of around 25%, which, however, was good and comparable.

Conclusions
Our results confirmed the low productive stability of camelina in the Mediterranean areas and highlighted how the effect of nitrogenous fertilization is dependent on meteorological conditions, especially for spring sowings.
Aside from these evaluations, the further value of this study is for demonstrating the possibility of applying remote and proximal techniques to camelina cultivation, with particular a ention to N management. This aspect is of fundamental importance for optimizing agronomic management and for increasing cropping system sustainability and quality. Among the different N fertilization rates and timing tested in the present study, satisfactory seed yields and high NUE, as well as good seed oil content, were reached with the application of 60 kg N ha −1 before sowing in spring camelina.
We find that UAV-derived vegetation indices during flowering are a good seed yield predictor in camelina. In particular, FGCC was the only VIs capable of providing an indication of camelina nutritional status. Anyway, further studies are needed to verify whether the UAV-derived VIs can give reliable indications on camelina nutritional status even during early vegetative phases. Although preliminary, these results can help in developing proper strategies for N application rates for camelina and for site-specific recommendations for its cultivation in Central Italy.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.