Comparative Study between Urea and Biogas Digestate Application towards Enhancing Sustainable Fertilization Management in Olive ( Olea europaea L., cv. ‘Koroneiki’) Plants

: Organic fertilization is a promising strategy to decrease N mineralization rates and high N losses via leaching and denitriﬁcation, thus synchronizing N application with N uptake for crops. A 230-day experiment with olive plants was realized under greenhouse conditions to compare urea and biogas digestate (BD) application on the growth, nutrient uptake, and physiological performance of Olea europaea L. plants. The following treatments were applied: (i) UREA, (ii) UREA + DCD (nitriﬁcation inhibitor), (iii) BD, (iv) BD + DCD, (v) acidiﬁed BD, (vi) acidiﬁed BD + DCD, (vii) CONTROL (typical native Marl soil type, without fertilizer application). Under BD application, signiﬁcantly higher levels of organic matter, Olsen P, and exchangeable K compared to UREA were found. Signiﬁcantly lower main shoot length and total plant biomass were recorded in the CONTROL soil compared to BD. Foliar N was lower in the ACID.BD and ACID.BD + DCD treatments, while leaf P was lower in UREA. The highest PSII activity was recorded in UREA + DCD, while the highest photosynthetic rate and intercellular CO 2 concentration were determined in UREA. It is expected that these data will constitute a ﬁrst comparative approach between urea and BD application in olive plants, which should be carefully considered, towards boosting sustainable fertilization in the frame of circular economy strategy.


Introduction
Olea europaea L. is the most important fruit tree crop of the Mediterranean basin. Despite its economic and ecological importance for the Mediterranean countries, many aspects on tree crop physiology and the field factors influencing olive productivity (e.g., temperature, as well as fertilization and irrigation) still remain unknown. For example, only very few published studies [1] exist with regard to the influence of fertilization on olive trees' nutrition, physiology, and yields. In addition, according to our knowledge, only Chatzistathis et al. (2020) [2] performed a comparative study between organic (manure application) and inorganic fertilization (controlled release fertilizer) on the growth, nutrient uptake, and physiology of Olea europaea L. plants.
Although N is the most important nutrient in fertilization programs [3,4] and especially for the nutrition of olive trees (it contributes to creation of new vegetation, flowering) [5], less than 20% of the applied N to orchards seems to be recovered by fruit trees. Thus, used in the hyper-intensive olive management systems of hedgerows (receiving high N fertilization inputs).

Plant Material, Biogas Digestate Characteristics, and Treatments
Olive plants (one-year olive rooted cuttings, coming from a commercial nursery, approximately 1 m in height), after randomization (based on their initial height), they were transplanted into 50 L pots (each one containing 20 kg soil) and were divided into 7 similar groups. The plants were grown for 230 days (from the 7th of November until the 25th of June) inside an experimental greenhouse (under natural light conditions, with minimum, maximum and average temperature of 10, 32 and 25 • C, respectively, while the relative humidity was ranged from 60 to 78%) of the Department of Forestry and Natural Environment, in the Aristotle University of Thessaloniki. The following 7 treatments (coinciding into one of the 7 groups mentioned above) were applied, with the rate of 261 mL/pot UREA or BD (corresponding to an application of 1402 mg N/pot, either in BD or UREA form): (a) plants fertilized with urea, (b) plants fertilized with urea with supplementary application of the nitrification inhibitor DCD (2.42 g), (c) plants fertilized with the organic material (bio-fertilizer) biogas digestate (BD), (d) plants fertilized with BD with supplementary application of DCD (2.42 g), (e) plants fertilized with acidified BD (pH adjusted to 6.3 from 7.8, with addition of 0.8 mL H 2 SO 4 /100 mL BD), (f) plants fertilized with acidified BD with supplementary application of DCD (2.42 g), and (g) plants grown on a Marl native soil type that receives no fertilization (CONTROL). The dose of 261 mL/pot corresponds to an addition of 150 kg N/ha, which is the typical N dose in olive groves. All the fertilizer additions were applied once at the beginning of the experiment. During the whole experimental period, all the plants were irrigated with an automated spray irrigation system with 200 mm of high quality water to achieve soil moisture in approximately 70% of the water holding capacity. In each of the 7 treatments, 5 plant-replicates were included; thus, the total number of experimental plants was 35. Biogas digestate (BD), used as bio-fertilizer, was a nonfiltered liquid substrate that contained both liquid and solid fractions (generated during anaerobic processes) of domestic animal (cow) origin and all the concentrations of toxic metals (lead, arsenic, nickel, chromium) were below the detectable limits. More specifically, BD had the following nutrient content: N: 0.55% D.W., P: 0.16% D.W. and K: 0.30% D.W., Ca: 2.24% D.W., Mg: 0.58% D.W., Fe: 3.7% D.W., Mn: 512 mg kg −1 D.W., Zn: 500 mg kg −1 D.W., and Cu approximately 200 mg kg −1 D.W. These nutrient concentrations were determined according to standard international protocols existing for vegetal tissue nutrients and metals' concentrations (such as Kjeldahl for N and ICP for the other nutrients and metals, which are described in detail in the Section 2.4).

Soil Sampling and Lab Analyses
A typical native Marl soil type was chosen as a substrate for olive plants' growth from the region of Epanomi, Thessaloniki, Central Macedonia (northern Greece). This soil type was chosen because is one of the 3 main soil types for olive culture in the region of Thessaloniki and Chalkidiki, Central Macedonia [20]. Soil sampling was performed from the upper 60 cm, since most of the active olive trees' root system exists [5]. Ten soil sub-samples were randomly selected, collected and afterwards mixed (before filling pots and performing soil analyses) to achieve homogenization [2]. After mixing, the unified (homogenized) soil sample was transferred to the experimental greenhouse to fill the 50 L pots.
The unified (homogenized) soil sample was dried inside the greenhouse and the stones were removed; afterwards, it was sieved to pass a 10-mesh screen and was transferred to lab for chemical analyses; its main properties are presented in the Table 1. Lab analyses included pH, particle size analysis, organic matter, % CaCO 3 , NO 3 -N, available P, exchangeable cations (Ca, Mg and K) and micronutrients (Fe, Mn, Zn and Cu). The above-mentioned parameters were determined as follows: pH in a soil-distilled water paste (1:1) [21], particle size analysis according to the Bouyoucos method [22], organic matter with potassium dichromate (K 2 Cr 2 O 7 ) [23], NO 3 -N concentration according to the VCl 3 /Griess method [24], while available P was determined with the Olsen method [25]. Finally, the exchangeable cations were determined with the ammonium acetate (CH 3 COONH 4 ) method [26] and micronutrients' concentrations after the extraction of 10 g soil with DTPA, pH 7.3 [27]. The concentrations of Ca, Mg, K, Fe, Mn, Zn, and Cu were determined by ICP (OPTIMA 2100 DV, Optical Emission Spectrometer, Perkin Elmer, Waltham, MA, USA) [28].

Plant Growth Data
During experimentation (from day 0, i.e., in the beginning of the experiment until the end of it-day 230), the main shoot elongation (cm/day) of all the olive plants was measured. In addition, at the end of the experimental period, apart from the main shoot length, the following plant growth parameters were determined: weight of basal and apical leaves, weight of basal and apical stems, weight of root system and the ratio shoot/root. More specifically, the main shoot of each experimental plant was divided into two equal parts (i.e., the basal and apical shoot), and afterwards root, leaves, and stems of each part were separated from each other; before being washed, the fresh weights (F.W.) of the root, stem, and leaves were measured. By adding the F.W. of all the plant parts, the total plant F.W. was calculated. After all the plant tissues were carefully washed (once with tap and twice with distilled water), they were dried at 75 • C, for 48 h. Then, the dry weights (D.W.) of all the parts (root, basal, and apical leaves; basal and apical stems) were determined. By adding the D.W. of all the plant parts, the total plant D.W. was calculated. Finally, the ratios between the shoot (i.e., leaves + stem) and root were determined, both on the F.W. and D.W. basis.

Tissue Nutrient Analyses, Total Plant Nutrient Content, and N Use Efficiency (NUE)
At the end of the experimental period, after the plants were harvested and the tissues were separated, washed, and dried, they were ground to a fine powder to pass a 30-mesh screen. A portion of 0.5 g. of the fine powder of each sample was dry-ashed in a muffle furnace at 515 • C for 5 h. Then, the ash was dissolved with 3 mL of 6 N HCl, diluted with double distilled water up to 50 mL, and the concentrations of P, K, Ca, Mg, Fe, Mn, Zn, and Cu were determined by ICP (Perkin Elmer-Optical Emission Spectrometer, OPTIMA 2100 DV) [28].
A 0.25 g portion of sample powder was used for N determination. Nitrogen was determined by the Kjeldahl method [29]. Macronutrient (N, P, K, Ca and Mg) concentrations were expressed in % D.W., while those of micronutrients (Fe, Mn, Zn and Cu) were expressed in mg kg −1 D.W. Multiplying the concentration of each nutrient (µg or mg g −1 D.W.) found in each plant part by the D.W. of the corresponding part, the content (absolute quantity) of each nutrient contained in each plant part at the end of this experiment was calculated. Through the addition of the nutrient contents of different plant parts, the total nutrient content (mg or g) per plant, and thus the total nutrient uptake per plant, was computed. Finally, N use efficiency (NUE) was defined as the amount of total biomass (d.w.) produced per unit of N [30].

Chlorophyll Fluorescence and Gas Exchange Measurements
In the middle (February), as well as at the end of the experimental period (June), the chlorophyll fluorescence parameters F v /F m and Performance Index (PI) were deter- mined. Particularly, the ratio F v /F m shows the maximum quantum yield of photosystem II (PSII) (F 0 : minimum fluorescence, F m : maximum fluorescence, F v = F m − F 0 : variable fluorescence) and PI reflects the functionality of both photosystems I and II and provides quantitative information on the current state of plant performance under stress conditions [31]. These parameters were determined in the youngest, fully expanded, mature (taken their final size) leaves from the upper part of the main shoot of plants using the PAM-2000 fluorometer (Heinz Walz GmbH, Effeltrich, Germany), after preconditioning the leaves in the dark for 20 min [32].
For the gas exchange measurements (i.e., photosynthetic rate, stomatal conductance, and intercellular CO 2 concentration), the LCPRO portable gas exchange measuring system (ADC Bioscientific Ltd., Hoddesdon, UK) was used. The measurements were performed in the youngest, mature, fully expanded leaves from the upper part of the main shoot. Measurements were performed both in the middle (February) and at the end of the experiment (June) from 10:00 to 12:00 a.m. at natural, full light intensity.

Statistical Analysis
The experimental design consisted of a 7 × 1 completely randomized factorial, with 7 fertilization treatments and 1 olive cultivar. In each of the seven treatments, five plantsreplicates (one per pot) were included (the total number of the experimental plants was 7 × 5 = 35). The data were statistically analyzed by the SPSS statistical program (ONE-WAY ANOVA); particularly, for the comparison of the mean values among the treatments, the Duncan's multiple range test for p ≤ 0.05 was used. The data for plant physiology measurements (chlorophyll fluorescence parameters, photosynthetic rate, stomatal conductance, and intercellular CO 2 concentration) were statistically analyzed by TWO-WAY ANOVA, and for the comparison of the mean values among the treatments and the time of measurements (February or June), the Tukey's test for p ≤ 0.05, was used. Table 1 shows the results of initial soil fertility (i.e., the parameters describing the fertility of soil substrate before the beginning of the experiment, i.e., before application of urea or biogas digestate-BD). As it is clear from Table 1, texture of the original soil was sandyclay-loam (SCL), pH was slightly above 7 (7.26), organic matter content was relatively poor (1.32%), NO 3 -N concentration was approximately 15 mg kg −1 , while extractable-according to Olsen method-P was poor (6.5 mg kg −1 ) ( Table 1). Exchangeable cations' concentrations (those of K, Ca and Mg) were within the normal levels of sufficiency, while from the micronutrients, only Fe concentration was insufficient (approximately 3.4 mg kg −1 ); in contrast, those of Mn, Zn, and Cu were sufficient (approximately 9.3, 1.0, and 0.7 mg kg −1 , respectively) ( Table 1).

Soil Fertility
With regard to final soil fertility (i.e., the parameters describing the fertility of soil mixtures at the end of the experiment), Table 2 concludes that pH did not significantly differ among the treatments. In contrast, organic matter content was significantly higher in the BD and BD + DCD treatments (>1.7%) compared to UREA and CONTROL (approximately 1.3%) ( Table 2). The highest concentrations of NO 3 -N were determined in UREA and BD (>13 mg kg −1 ), which were significantly higher to those found in the ACID BD + DCD and CONTROL. Olsen P was 3-5 times lower in the UREA treatments (UREA, UREA + DCD) and control (about 4.5-5.5 mg kg −1 ) compared to the organic treatments (BD, BD + DCD, ACID. BD, ACID. BD + DCD) (approximately 15.5-28 mg kg −1 ) ( Table 2). From the exchangeable cations (K, Ca and Mg), only K concentrations were significantly different among the treatments; more specifically, exchangeable K was significantly higher in the ACID. BD and ACID. BD + DCD, compared to the CONTROL, UREA, UREA + DCD and BD treatments ( Table 2). From the micronutrients (Fe, Mn, Zn, Cu and B), only Zn, Cu, and B concentrations were found to be significantly higher in the organic treatments  Table 2).

Plant Growth (Main Shoot Length, Plant Tissue Weights, and Total Plant Biomass)
For many growth parameters, significantly lower values were recorded in the control plants compared to other treatments; the ratio shoot/root was significantly lower in the ACID. BD and CONTROL plants ( Table 3). The elongation rate showed significant fluctuations among the treatments and time period (Tables 4 and 5). However, at the end of the experiment, significant differences were found in the main shoot length only between (i) the BD and (ii) the ACID. BD + DCD and CONTROL; all other differences were insignificant ( Figure 1). Finally, the d.w. of the total plant biomass was significantly lower in the CONTROL plants, compared to those treated with BD, ACID. BD, BD + DCD, and ACID. BD + DCD ( Figure 2).

Foliar Nutrition, Total Plant Nutrient Content and N Use Efficiency (NUE)
The highest foliar N (>2% D.W., both in the basal and apical leaves) concentrations were found in the UREA + DCD treatment; in contrast, the lowest leaf N levels were determined in the ACID.BD, ACID.BD + DCD, and CONTROL plants; with regard to the other macronutrients (P, K, Ca and Mg), only in the case of P, in the basal leaves, was found significantly lower concentration (0.08% D.W.) in UREA, compared to the BD + DCD treatment (0.12% D.W.) (Table 6). Regarding micronutrients, significantly lower foliar Mn (for apical leaves) were found in the ACID.BD and ACID.BD + DCD treatments. For Zn, significantly lower foliar levels (both for basal and apical leaves) were found in BD + DCD. Significantly higher leaf Cu concentrations (both in basal and apical leaves)

Foliar Nutrition, Total Plant Nutrient Content and N Use Efficiency (NUE)
The highest foliar N (>2% D.W., both in the basal and apical leaves) concentrations were found in the UREA + DCD treatment; in contrast, the lowest leaf N levels were determined in the ACID.BD, ACID.BD + DCD, and CONTROL plants; with regard to the other macronutrients (P, K, Ca and Mg), only in the case of P, in the basal leaves, was found significantly lower concentration (0.08% D.W.) in UREA, compared to the BD + DCD treatment (0.12% D.W.) (Table 6). Regarding micronutrients, significantly lower foliar Mn (for apical leaves) were found in the ACID.BD and ACID.BD + DCD treatments. For Zn, significantly lower foliar levels (both for basal and apical leaves) were found in BD + DCD. Significantly higher leaf Cu concentrations (both in basal and apical leaves) were determined in the UREA + DCD, compared to the organic treatments and CONTROL (Table 7). Numbers are means of five replicates (N = 5); the different letters in the same column symbolize statistically significant differences among the treatments according to the Duncan's multiple range test for p ≤ 0.05.
The highest total plant N content (approximately 2.5 g) was found in the UREA treatment; significant differences in N content were determined between: (i) UREA and (ii) BD + DCD, ACID.BD, ACID.BD + DCD, and CONTROL plants ( Figure 3A). In contrast, the highest P content was found in the BD + DCD, ACID.BD, and ACID.BD + DCD treatments, compared to the UREA, UREA + DCD, and CONTROL plants ( Figure 3B). With regard to K content, significant differences were determined only between: (i) the treatments ACID.BD and ACID.BD + DCD and (ii) the UREA + DCD and CONTROL plants ( Figure 3C). Finally, with regard to total plant Ca and Mg contents, significant differences were recorded only between UREA and CONTROL, and between the (i) ACID.BD and (ii) UREA + DCD, BD + DCD and CONTROL plants, respectively ( Figure 3D,E).   Similarly, significantly higher Fe content was recorded in the BD, ACID.BD, and ACID.BD + DCD treatments compared to UREA + DCD and CONTROL ( Figure 4A). Significantly higher Zn content was found in the ACID.BD + DCD compared to UREA, UREA + DCD, BD, BD + DCD, and CONTROL ( Figure 4C), while the lowest Cu contents were determined in the BD, BD + DCD, and ACID.BD treated plants (which were significantly lower to those found in the UREA + DCD and ACID.BD + DCD) ( Figure 4D).   Figure 5 describes the photosystem II activity (as expressed by the parameters Fv/Fm and Performance Index-PI). The highest (and significantly higher) value of the ratio Fv/Fm was found at the end of the experimental period (June) in the UREA + DCD, compared to the other treatments ( Figure 5A). Similar to Fv/Fm, the Performance Index (PI)  3.4. Photosystem II Activity (F v /F m and Performance Index-PI) Figure 5 describes the photosystem II activity (as expressed by the parameters F v /F m and Performance Index-PI). The highest (and significantly higher) value of the ratio F v /F m was found at the end of the experimental period (June) in the UREA + DCD, compared to the other treatments ( Figure 5A). Similar to F v /F m , the Performance Index (PI) was significantly higher in the UREA + DCD in June compared to all the other treatments, both in the middle (February) and at the end of the experiment (June). The lowest PI value was recorded in the CONTROL plants in June, compared to the other treatments ( Figure 5B).

Photosynthetic Rate, Stomatal Conductance, Transpiration Rate, and Intercellular CO2 Concentration
The highest photosynthetic rate in February was recorded in the UREA treatment; this value (approximately 11 μmol CO2 m −2 s −1 ) was significantly higher to those determined in the UREA + DCD and BD treatments for the same period ( Figure 6A). In June, significant differences were found only between the UREA and UREA + DCD treatments. Generally, significantly higher photosynthetic rates were recorded in February, than in June ( Figure 6A). In contrast, stomatal conductance, in most treatments, was significantly lower in February, than in June ( Figure 6B). Intercellular CO2 concentration in February was significantly higher in UREA compared to most of the other treatments; in contrast, for the same treatment in June, the lowest intercellular CO2 concentration was recorded compared to all the other treatments ( Figure 6C).

Photosynthetic Rate, Stomatal Conductance, Transpiration Rate, and Intercellular CO 2 Concentration
The highest photosynthetic rate in February was recorded in the UREA treatment; this value (approximately 11 µmol CO 2 m −2 s −1 ) was significantly higher to those determined in the UREA + DCD and BD treatments for the same period ( Figure 6A). In June, significant differences were found only between the UREA and UREA + DCD treatments. Generally, significantly higher photosynthetic rates were recorded in February, than in June ( Figure 6A). In contrast, stomatal conductance, in most treatments, was significantly lower in February, than in June ( Figure 6B). Intercellular CO 2 concentration in February was significantly higher in UREA compared to most of the other treatments; in contrast, for the same treatment in June, the lowest intercellular CO 2 concentration was recorded compared to all the other treatments ( Figure 6C). ments. Generally, significantly higher photosynthetic rates were recorded in February, than in June ( Figure 6A). In contrast, stomatal conductance, in most treatments, was significantly lower in February, than in June ( Figure 6B). Intercellular CO2 concentration in February was significantly higher in UREA compared to most of the other treatments; in contrast, for the same treatment in June, the lowest intercellular CO2 concentration was recorded compared to all the other treatments ( Figure 6C).

Discussion
Final pH (after the end of the experiment) was 0.8-1 unit higher, compared to the initial pH (before the beginning of the experiment) after BD and UREA application. Manure application increases soil pH and salt accumulation; thus, the risk of soil salinization is increased by the continuous application of animal manure [33]. The fact that organic matter was higher in the BD applications compared to the UREA treatment and CONTROL (Table 2), is ascribed to the organic C input after BD addition. These data agree with those of Chatzistathis et al. (2020) [2], who found that, after cow and goat manure application, soil organic matter was significantly higher in the organic treatments (3.49-3.91% d.w.) compared to the control soil (1.55% d.w.). The lowest NO3-N concentration determined in the CONTROL soil (4.85 mg kg −1 d.w.) (Tables 1 and 2), which may be ascribed to the fact that in the CONTROL soil no N input was received, while in the UREA and BD treatments N inputs, were applied. In the BD and acidified BD treatments, the fact that higher P, K, Zn, and Cu concentration were recorded compared to UREA and CONTROL (Table 2), should be ascribed to the fact that BD contains approximately 0.16% d.w. P, 0.30% d.w. K, 500 mg/kg d.w. Zn, and approximately 200 mg/kg d.w. Cu (described in the Section 2.1.), while UREA contains only N (its type is 46-0-0).
The fact that insignificant differences in the main shoot length were recorded: (i) between UREA and BD (or even higher total plant biomass was determined in ACID. BD compared to UREA) (Figure 1)), together with: (ii) the significantly lower total plant biomass (in a d.w. basis) recorded in the CONTROL compared to all the BD treatments (Figure 2), means that BD application was likely the most beneficial strategy to support plant growth; thus, it needs more attention as organic fertilizer on the part of olive growers for further use under field conditions to sustain high crop yields and enhance the sustainable fertilization management of olive groves. In addition, within its advan-

Discussion
Final pH (after the end of the experiment) was 0.8-1 unit higher, compared to the initial pH (before the beginning of the experiment) after BD and UREA application. Manure application increases soil pH and salt accumulation; thus, the risk of soil salinization is increased by the continuous application of animal manure [33]. The fact that organic matter was higher in the BD applications compared to the UREA treatment and CONTROL (Table 2), is ascribed to the organic C input after BD addition. These data agree with those of Chatzistathis et al. (2020) [2], who found that, after cow and goat manure application, soil organic matter was significantly higher in the organic treatments (3.49-3.91% d.w.) compared to the control soil (1.55% d.w.). The lowest NO 3 -N concentration determined in the CONTROL soil (4.85 mg kg −1 d.w.) (Tables 1 and 2), which may be ascribed to the fact that in the CONTROL soil no N input was received, while in the UREA and BD treatments N inputs, were applied. In the BD and acidified BD treatments, the fact that higher P, K, Zn, and Cu concentration were recorded compared to UREA and CONTROL (Table 2), should be ascribed to the fact that BD contains approximately 0.16% d.w. P, 0.30% d.w. K, 500 mg/kg d.w. Zn, and approximately 200 mg/kg d.w. Cu (described in the Section 2.1), while UREA contains only N (its type is 46-0-0).
The fact that insignificant differences in the main shoot length were recorded: (i) between UREA and BD (or even higher total plant biomass was determined in ACID. BD compared to UREA) (Figure 1)), together with: (ii) the significantly lower total plant biomass (in a d.w. basis) recorded in the CONTROL compared to all the BD treatments (Figure 2), means that BD application was likely the most beneficial strategy to support plant growth; thus, it needs more attention as organic fertilizer on the part of olive growers for further use under field conditions to sustain high crop yields and enhance the sustainable fertilization management of olive groves. In addition, within its advantageous aspects of application should be included the significantly higher NUE values in the ACID. BD treatments compared to UREA and UREA + DCD (data not-shown). However, in order to include a more thorough approach of the advantages of BD application (in the frame of circular economy) on the enhancement of sustainable fertilization management, long-term field studies are needed in the near future.
The fact that the highest leaf N (for both basal and apical leaves) concentrations were found in the UREA + DCD treatment (approximately 2.20 and 2.40% dw, respectively) ( Table 6) probably means that the addition of DCD in urea restricted N losses via denitrification, volatilization, and leaching and improved N uptake. One of the most problematic issues in modern agriculture is the high quantity of ammonia being released, resulting in a partial loss of the N, which is contained in urea fertilizers [34]; according to the same authors, the best results were obtained with the use of the urease inhibitor NBPT, which can decrease ammonia emissions from urea fertilizers by 30-70% [34]. Corrochano-Monsalve et al. (2021) [14] found that application of urea with double inhibitor diminished N losses in the form of NH 3 volatilization by more than 50% compared to single urea application, while the co-application of biochar and ammonium sulphate, with the nitrification inhibitor DMPP, increased N uptake up to 34% and N use efficiency up to 13.7% [35]. Centeno et al. (2017) [13], who studied the effects of N fertilization and nitrification inhibitors on vegetative growth, production, and olive oil quality of 'Arbequina' hedgerow and 'Picual' vase-trained orchards, found that the effect of nitrification inhibitor (DMPP) on olive growth or production was not consistent, and N fertilization was not always necessary, since, when leaf N was higher than 2% d.w., fruit d.w., oil content, and oil phenolic content were reduced. Although no direct comparisons of our data to those of Centeno et al. (2017) [13] can be done (due to the different conditions-field vs. large pot experimentationof the two studies), it is clear that in the UREA treatments (especially in the UREA + DCD treatment) foliar N levels exceeded 2%. The optimum N levels from mature non-fruiting olive leaves from the middle of annual vegetation should vary from 1.6% to 1.8% d.w. [5].
The lowest P concentration in basal leaves (0.08% d.w.) and the lowest total P content were found in UREA (Table 6, Figure 3) and this is probably ascribed to the lack of P in urea (46-0-0). The differences among the treatments in content for the other nutrients (K, Ca, Mg, and micronutrients) may be ascribed either to the differences in nutrient uptake or to the differences in biomass production (or to both of them). Similar differences in total plant nutrient content among the inorganic/organic fertilization treatments (ascribed either to differences in nutrient uptake or to differences in total plant dry weight) were also found in the study of Chatzistathis et al. (2020) [2]. With regard to micronutrients, it is interesting to note that the highest foliar Cu concentrations were recorded in UREA + DCD and UREA (Table 6), despite the fact that the lowest soil Cu availability occurred in these two treatments ( Table 2). This may be ascribed to a possible synergistic effect between Cu and Zn in uptake by olive plants (since similar behavior was also observed for soil and foliar Zn-Tables 2 and 6). The kind of interaction effect (synergistic or antagonistic) between divalent cations (e.g., Zn 2+ and Cu 2+ ) is highly dependent on plant species and/or genotypes [36]; in contrast to our results, other researchers found inhibitory or no effects from influence between Zn 2+ and Cu 2+ for other plant species [37].
Photosystem II activity (as expressed by the parameters F v /F m and Performance Index-PI) was significantly higher in the UREA + DCD treatment in June ( Figure 5). These findings may be ascribed to the highest foliar Mn concentrations (63 and 52 mg/kg for basal and apical leaves, respectively) determined in UREA + DCD (Table 6). Manganese is mainly involved in photosystem II activity [38] and foliar Mn levels varying from 50 to 150 mg/kg in the youngest mature olives leaves are the optimum ones for PSII activity [5]. In February, the maximum photosynthetic rate was recorded in UREA ( Figure 6A); this was ascribed to the significantly higher stomata opening recorded in this treatment (Fig-ure 6B). In June, significantly higher photosynthetic rates occurred in UREA + DCD and BD + DCD, compared to UREA ( Figure 6A); however, high stomatal conductance values in June were recorded not only in UREA + DCD and BD + DCD, but also in UREA and ACID.BD treatments ( Figure 6B). This could likely be ascribed to the differences in CO 2 efficiency (photosynthetic rate/intercellular CO 2 ) among the treatments, i.e., in some cases, although high stomata openings were recorded, the lower CO 2 efficiency might lead to lower photosynthetic rates.

Conclusions
From all the data presented in this study, it is clear that BD application improved soil properties (organic matter, Olsen P, exchangeable K, and DTPA extractable Cu) compared to UREA, which proves its beneficial role as soil amendment. Regarding the influence of N treatment and nitrification inhibitor on N uptake, it was found that the highest foliar N was determined in the UREA + DCD and the lowest leaf N levels were found in the ACID.BD and ACID.BD + DCD treatments compared to the other ones, something which probably shows the lower N mineralization rate and uptake under these two treatments. In the case of UREA + DCD, the higher leaf N concentration may be ascribed to the possible effect of the nitrification inhibitor on the decrease of N losses via leaching and denitrification, thus causing the increase of N uptake by the olive plants. However, further research, including N losses via denitrification and leaching, should be realized in the near future to test this hypothesis. For the other nutrients, the lowest total plant contents were recorded in the CONTROL soil, which is mainly ascribed to the lowest total plant biomass, rather than to lower nutrient uptake. Significantly higher NUE values were found in ACID. BD and ACID. BD + DCD compared to the UREA treatments (UREA, UREA + DCD), which shows that the optimum treatments for NUE were those of ACID. BD. Finally, the highest PSII efficiency at the end of the experiment in the UREA + DCD treatment is likely owed to the optimum foliar Mn levels determined in this treatment. It is expected that the above data will constitute a first step towards a comparative approach between UREA and BD application in super-intensive olive groves, with the aim to decrease high (often unnecessary) N fertilization rates and boost sustainable nutrient management in the frame of circular economy strategy. However, long-term field studies are needed in the near future to obtain stable conclusions.