Cultivar Dependent Impact on Yield and Its Components of Young Almond Trees under Sustained-Deficit Irrigation in Semi-Arid Environments

Almond (Prunus dulcis Mill. (D.A. Webb)) plantations in irrigated semi-arid areas need to successfully face the new scenarios of climate change combining sustainable irrigation strategies and tolerant cultivars to water stress. This work examines the response of young almond (cvs. Guara, Marta, and Lauranne) subjected to different irrigation doses under semi-arid conditions (South-West Spain). The trial was conducted during two seasons (2018–2019) with three irrigation strategies: A full-irrigated treatment (FI), which received 100% of the irrigation requirements (IR), and two sustained-deficit irrigation strategies that received 75% (SDI75) and 65% (SDI65) of IR. Crop water status was assessed by leaf water potential (Ψleaf ) and stomatal conductance (gs) measurements, determining the yield response at the end of each season. Different physiological responses for the studied cultivars were observed, especially considering the Ψleaf measurements. In this way, cv. Marta behaved more tolerant, while cvs. Guara and Lauranne maintained higher gs rates in response to water stress. These differences were also observed in terms of yield. The cv. Lauranne did not reflect yield losses, and the opposite trend was observed for cv. Guara, in which reductions on fruit numbers per tree were detected. On overall, effective irrigation water savings (≈2100 m3·ha−1 in SDI65) could be feasible, although these responses are going to be substantially different, depending on the used cultivar.


Introduction
Today, deficit irrigation (DI) strategies cannot be considered as a novelty, and their effects and consequences have been widely studied in order to improve water resources management under water shortage conditions. With them, irrigation can be reduced until certain levels, with the main aim of maximizing the water savings keeping the yields within an acceptable range, close to those obtained under non-water restrictions [1]. Several authors have reported positive results when these divergent in response to SDI strategies, which ultimately effects on yield and its components. Thus, the aim of this study was to assess the nut yield and physiological response of three commercial almond cultivars namely Guara, Marta, and Lauranne subjected to sustained-deficit irrigation strategies, elucidating the tolerant cultivar under these strategies in a semi-arid Mediterranean environment.

Experimental Site
The trial was conducted during two consecutive years (2018-2019) in a commercial almond orchard (Prunus dulcis Mill. (D.A. Webb) cvs. Guara, Marta, and Lauranne ( Figure 1); grafted onto GN15 rootstock), and located in the Guadalquivir river basin (SW Spain, 37 • 30 27.4" N, 5 • 55 48.7" W). Trees were planted in 2013, spaced 8 × 6 m, and drip irrigated using two pipelines with emitters of 2.3 L·h −1 , at 0.75 m intervals. Canopy volumes were very similar within each cultivar. That is, for cv. Marta, canopy volumes ranged from 64 to 65 m 3 ; for cv. Guara between 65 and 66 m 3 ; and for the case of cv. Lauranne, some higher than the previous cultivars, with canopy volumes between 72 and 74 m 3 .
Agronomy 2020, 10, x FOR PEER REVIEW 3 of 15 components. Thus, the aim of this study was to assess the nut yield and physiological response of three commercial almond cultivars namely Guara, Marta, and Lauranne subjected to sustaineddeficit irrigation strategies, elucidating the tolerant cultivar under these strategies in a semi-arid Mediterranean environment.

Experimental Site
The trial was conducted during two consecutive years (2018-2019) in a commercial almond orchard (Prunus dulcis Mill. (D.A. Webb) cvs. Guara, Marta, and Lauranne ( Figure 1); grafted onto GN15 rootstock), and located in the Guadalquivir river basin (SW Spain, 37°30′27.4", 5°55′48.7" W). Trees were planted in 2013, spaced 8 × 6 m, and drip irrigated using two pipelines with emitters of 2.3 L·h −1 , at 0.75 m intervals. Canopy volumes were very similar within each cultivar. That is, for cv. Marta, canopy volumes ranged from 64 to 65 m 3 ; for cv. Guara between 65 and 66 m 3 ; and for the case of cv. Lauranne, some higher than the previous cultivars, with canopy volumes between 72 and 74 m 3 . The soil is a silty loam typical Fluvisol, more than 2 m deep, with organic matter <1.5%. Roots were located predominately in the first 50 cm soil depth, corresponding to the intended wetting depth. Soil water content values at field capacity (−0.033 MPa) and permanent wilting point (−1.5 MPa) were close to 0.40 and 0.15 m 3 ·m −3 , respectively. The climatic classification of the study area is attenuated meso-Mediterranean, with an annual ET0 rate of 1400 mm and accumulated rainfall of 540 mm (average data corresponding to the last 15 years; obtained from the Andalusian Weather information Network [27]).
Regarding to the experimental conditions registered in the irrigation period (from March to October) during the monitored years, in 2018 the average temperature was 20.5 °C with minimum and maximum average temperatures of 7.3 and 37.7 °C, respectively. Similarly, during the second year (2019), the average temperature was 25 °C with minimum and maximum average values of 7.2 and 36.5 °C, respectively. In relation to the average vapor pressure deficit (VPD) in 2018 it ranged from 1.20 to 3.36 KPa, while during 2019 it was between 1.10 and 2.49 KPa. Reference evapotranspiration (ET0) and rainfall during the irrigation period (March-October) amounted to 1102 and 326 mm for 2018, meanwhile, in the second year ET0 and rainfall were 1221 and 85 mm, respectively. Crop evapotranspiration (ETC) rates for both seasons averaged 860 mm, which was in line with previous reports by López-López et al. [28] or García-Tejero et al. [11] who estimated water requirements of 800 and 900 mm for mature almond trees in Guadalquivir river basin. The soil is a silty loam typical Fluvisol, more than 2 m deep, with organic matter <1.5%. Roots were located predominately in the first 50 cm soil depth, corresponding to the intended wetting depth. Soil water content values at field capacity (−0.033 MPa) and permanent wilting point (−1.5 MPa) were close to 0.40 and 0.15 m 3 ·m −3 , respectively. The climatic classification of the study area is attenuated meso-Mediterranean, with an annual ET 0 rate of 1400 mm and accumulated rainfall of 540 mm (average data corresponding to the last 15 years; obtained from the Andalusian Weather information Network [27]).
Regarding to the experimental conditions registered in the irrigation period (from March to October) during the monitored years, in 2018 the average temperature was 20.5 • C with minimum and maximum average temperatures of 7.3 and 37.7 • C, respectively. Similarly, during the second year (2019), the average temperature was 25 • C with minimum and maximum average values of 7.2 and 36.5 • C, respectively. In relation to the average vapor pressure deficit (VPD) in 2018 it ranged from 1.20 to 3.36 KPa, while during 2019 it was between 1.10 and 2.49 KPa. Reference evapotranspiration (ET 0 ) and rainfall during the irrigation period (March-October) amounted to 1102 and 326 mm for 2018, meanwhile, in the second year ET 0 and rainfall were 1221 and 85 mm, respectively. Crop evapotranspiration (ET C ) rates for both seasons averaged 860 mm, which was in line with previous reports by López-López et al. [28] or García-Tejero et al. [11] who estimated water requirements of 800 and 900 mm for mature almond trees in Guadalquivir river basin.
In both seasons, irrigation was applied from the middle of March to the end of October, the IR being estimated according to the methodology proposed by Allen et al. [29] (Equations (1) and (2)), obtaining the values of reference evapotranspiration (ET 0 ) by using a weather station installed in the same experimental orchard (Davis Advance Pro2, Davis Instruments, Valencia, Spain).
where ET C is the crop evapotranspiration; K C is the single-crop coefficient; K r is the crop reduction coefficient, which depends on the percentage of shaded area cast by the tree canopy; ET 0 is the reference evapotranspiration; and IR is the irrigation requirements. The local crop coefficients used during the experimental period ranged from 0.4 to 1.2 according to García-Tejero et al. [30]. Additionally, the IR was reduced for SDI 75 and SDI 65 by multiplying it by 0.75 and 0.65, respectively.

Plant Measurements
During the maximum evapotranspirative demand period, coinciding with the kernel-filling and pre-harvest stages (from early June to mid-August; 160-225 DOY in 2018; and 160-219 DOY in 2019), crop-water monitoring was done by means of measurements of leaf water potential (Ψ leaf ) and the stomatal conductance to water vapor (g s ); these readings being taken between 12:00 and 13:30 GTM, and with a periodicity of 10-15 days. The Ψ leaf was measured using a pressure chamber (Soil Moisture Equipment Corp., Sta. Barbara, CA, USA), monitoring 8 trees per irrigation treatment and two leaves per tree, located in the north side of the tree and being totally mature, fresh and shaded; at 1.5 m of height, approximately and NW exposed. In these same trees, we measured the g s , using a porometer SC-1 (Decagon Devices, INC, Pullman, WA, USA), these measurements being done in two leaves per monitored tree, completely exposed to the sun, at 1.5 m of height, and preferably with south-eastern facing.
Additionally, with the aim of quantifying total water stress supported by the crop, the stress integral (SI) was calculated from Ψ leaf and g s data, following the methodology proposed by Myers et al. [31] (Equations (3) and (4)): where SI lea f is the stress integral in terms of leaf-water potential values, Ψ av lea f is the average leaf water potential for any interval; Ψ max lea f is the maximum value of leaf-water potential registered during the experimental period; SI gs is the stress integral in terms of stomatal conductance values, g av s is the average stomatal conductance for any interval; g min s is the minimum value of stomatal conductance during the experimental period; and n is the days numbers within each interval.
According to these indexes, higher water stress gathered by the crop would be related to higher values of SI lea f and lesser values of SI gs . In this sense, these indexes provide information to quantify the crop water stress accumulated along a period, allowing quantifying the effect of water restriction beyond its temporal distribution, integrating the global stress supported by the crop in comparison to the punctual measurements.
At the end of each season it was monitored the almond yield in terms of nut and kernel weight. Harvesting was carried out by using a mechanical vibrator with a mechanical peeling to remove the hull. Once cleaned, almonds were left to air dry and weighed once they reached a humidity content of around 6%. Finally, almonds were processed with shelling machine, obtaining the kernel yield for each irrigation treatment and cultivar.
In relation to the harvesting dates these were different for each cultivar and season. In 2018, the almond harvest labors were done at 232 DOY for cv. Guara and 239 DOY for cvs. Marta and Lauranne; meanwhile in 2019, these were done at 219, 221, and 235 DOY for Guara, Marta, and Lauranne, respectively.
Taking into consideration that two components determines the yield (the almond size and number of almonds per tree); the first one was obtained by weighing 100 almonds per monitored tree (n = 8); obtaining the kernel unit weight and ratio between kernel and nut (kernel + shell). After this, the second component (number of almonds per tree) was estimated by dividing "the kernel yield per monitored tree" by "the kernel unit weight", tracing the most affected yield component by water stress imposed in each cultivar.
Finally, and considering the total irrigation applied in each treatment, it was estimated the irrigation water productivity (IWP; kg·m −3 ), defined as the ratio between kernel yield and irrigation, water applied.

Experimental Design and Stadistical Analysis
The experimental design was of randomized blocks, with four replications per irrigation treatment and cultivar. Each replication had 12 trees (3 rows and 4 trees per row), and the two central trees for each replication were monitored. Thus, eight trees per treatment of irrigation strategy were monitored (n = 8).
Statistical analysis was developed by using the Sigma Plot statistical software (version 12.5, Systat Software, Inc., San Jose, CA, USA) and the SPSS software (SPSS Inc., 15.0 Statistical packages; Chicago, IL, USA). Year-to-year, an exploratory descriptive analysis of the whole of physiological measurements (Ψ leaf and g s ) for each treatment and cultivar was done; applying a Levene's test to check the variance homogeneity of the variables studied. After this, for each cultivar, an ANOVA for repeated measures was developed (three treatments and 2 freedom degrees), applying a Bonferroni test to compare pairs of treatments when significant differences in the ANOVA were detected. Moreover, and with the aim of identifying those days in which differences between irrigation treatments were detected, a one-way analysis of variance (ANOVA) and a Tukey's test were done for each measurement day.
Additionally, considering the results provided by SI gs and SI lea f ; for each irrigation treatment and cultivar, a two-way ANOVA was developed, followed by a Tukey's multiple range test. These results allow development of a simpler statistical analysis (considering that all measurements were transformed into an stress index that summarizes the total water stress accumulated by the crop) and to find a single value to evaluate the water stress supported by each cultivar and irrigation treatment and to analyze the hypothetical differences between cultivars and treatments in physiological terms with ease. Moreover, the linear regressions between the average values of SI lea f and SI gs registered for each treatment in both years were obtained (n = 6), applying an analysis of covariance (ANCOVA) to evaluate the differences in the interception points and slopes of these regressions. These analyses would allow the determination of whether these relationships are cultivar dependent and if this dependence is accompanied by similar yield responses.
Finally, and year-to-year, the kernel yield and its components were analyzed (kernel unit weight, the ratio between kernel weight vs. almond weight (kernel + shell), and fruits number per tree); by applying a Levene's test to check the variance homogeneity and ANOVA with a Tukey's test, considering as factors, the irrigation treatment, the cultivar, and their interactions.

Physiological Response to Water Stress
In the course of the first experimental season (2018) the ANOVA for repeated measures did not show significant differences between irrigation treatments within each cultivar. That is, on overall as Ψ leaf as g s evidenced similar results for the studied treatments; although afterwards, the analysis done independently for each monitoring day reflected significant differences between treatments. These differences between treatments were especially noticeable for the case Ψ leaf (Figure 2), in comparison to the higher similarities reported by g s measurements ( Figure 3). When confronting the Ψ leaf values registered in each cultivar, these ranged from −0.9 to −2.3 MPa for cvs. Guara and Lauranne, meanwhile for cv. Marta these ranged between −0.8 and −2.0 MPa.
In relation to g s , the obtained values between treatments and cultivars were very similar, which was very evident by the absence of significant differences. Comparing the g s values for each cultivar, cv. Guara showed values between 80 and 175 mmol·m −2 ·s −1 , cv. Marta between 75 and 180 mmol·m −2 ·s −1 , and cv. Lauranne registered g s rates between 90 and 175 mmol·m −2 ·s −1 .
As happened in 2018, during 2019 the differences among treatments were higher in Ψ leaf than in g s . (Figure 2). Again, the ANOVA for repeated measures did not reflect significant differences between irrigation treatments. For cv. Guara, Ψ leaf values ranged from −1.4 to −1.9 MPa in FI, whereas in SDI treatments these values oscillated from −1.7 to −2.5 MPa. Somewhat lower were the values for cv. Lauranne with Ψ leaf between −1.1 and −1.7 MPa in FI, and for both SDI 75 and SDI 65 between −1.5 and −2.0 MPa. Also, some differences for cv. Marta were observed with Ψ leaf with values between −1.1 and −1.6 MPa in FI; and for both SDI treatments between −1.3 and −1.9 MPa. In relation to g s in 2019, no significant differences among treatments were observed for any studied cultivar (Figure 3), very similar to that reported in the previous season. In particular, for cv. Guara, the g s rates were between 163 and 314 mmol·m −2 ·s −1 whereas cv. Lauranne displayed similar rates between 161 and 275 mmol·m −2 ·s −1 . Finally, cv. Marta registered g s values that ranged from 150 to 310 mmol·m −2 ·s −1 , these being similar in all studied irrigation treatments.
In summary and considering the observed values during the two studied years, cv. Marta showed a higher capability to register higher values for Ψ leaf with respect to the other remaining studied cultivars, especially under SDI conditions.
More perceptible findings were achieved in relation to the water stress integral (Table 1). During the first experimental season (2018) significant effects in relation to cultivar (p < 0.05) were found in SI Ψleaf and SI gs . In this sense, cvs. Guara and Lauranne showed similar values of SI Ψleaf and SI gs , and significant (p < 0.05) higher than those calculated in cv. Marta. Particularly noticeable were the values of SI gs for cv. Marta, which were, globally, 38% and 30% lower than those registered in cvs. Lauranne and Guara, respectively. Respect to the irrigation doses effects, no significant differences were observed for SI Ψleaf and SI gs , this being in agreement with the absence of differences reported by the ANOVA for repeated measured previously discussed.
During 2019, the highest differences were observed again among cultivars, as for SI Ψleaf as for SI gs . Thus, the lowest and highest values of SI Ψleaf were reached by cvs. Marta and Guara, respectively, very similar to the response observed the previous season. Moreover, the highest differences in 2019 were determined in SI gs among cultivars, with values for cvs. Guara and Lauranne higher that those obtained for cv. Marta (p < 0.01). That is, irrigation doses did not promote differences in SI Ψleaf (as the previous season), different to the observed in SI gs with values for SDI 65 20% higher than those detected in FI.
Taking into account the interaction between both considered factors (irrigation × cultivar), these were detected only the second studied year, as for SI Ψleaf as for SI gs , the cultivar being the main factor of this interaction. More interesting were the relationships between the SI leaf and SI gs , which were estimated for each cultivar (Figure 4). Whereas cvs. Guara and Lauranne showed similar relationships, and the ANCOVA did not reveal differences for the intercepts and slopes between both cultivars; significant differences were registered with the interception point obtained in cv. Marta, although the slopes were very similar.  More perceptible findings were achieved in relation to the water stress integral (Table 1). During the first experimental season (2018) significant effects in relation to cultivar (p < 0.05) were found in SIΨleaf and SIgs. In this sense, cvs. Guara and Lauranne showed similar values of SIΨleaf and SIgs, and significant (p < 0.05) higher than those calculated in cv. Marta. Particularly noticeable were the values of SIgs for cv. Marta, which were, globally, 38% and 30% lower than those registered in cvs. Lauranne and Guara, respectively. Respect to the irrigation doses effects, no significant differences were   According to these results, when physiological response was analyzed by using SI le  af and SIgs, these differences were more evident, this being mostly relevant in comparing the type of cultivars. Considering previous results reported by Gutiérrez-Gordillo et al. [26], cv. Marta was more sensitive to water stress in physiological terms, evidencing a stronger stomatal control under RDI strategy. On the contrary, cv. Guara triggered physiological mechanisms that were able to maximize the gasexchange rates by increasing gs and decreasing leaf values. By contrasting these findings with those outlined by García-Tejero et al. [32], declines of leaf encourage to lessen the carbon assimilation rate, According to these results, when physiological response was analyzed by using SI leaf and SI gs , these differences were more evident, this being mostly relevant in comparing the type of cultivars. Considering previous results reported by Gutiérrez-Gordillo et al. [26], cv. Marta was more sensitive to water stress in physiological terms, evidencing a stronger stomatal control under RDI strategy. On the contrary, cv. Guara triggered physiological mechanisms that were able to maximize the gas-exchange rates by increasing g s and decreasing leaf values. By contrasting these findings with those outlined by García-Tejero et al. [32], declines of leaf encourage to lessen the carbon assimilation rate, disclosing the almond capability in maintaining high g s values even when leaf is close to −2.5 MPa [33]. Moreover, before reaching important exhaustion in g s rates, reductions in leaf are not accompanied by relevant depletions in g s [32] . This might happen because almond is competent at holding optimum g s and carbon assimilation rates, improving the water-use efficiency, at least under moderate water-stress situations. Therefore, almond reaction to water stress would demand different regulation mechanisms to counteract the adverse impact of water stress as was corroborated with the results of the present study. In addition, according to Fernández et al. [34] and Fu et al. [35], the almond has a lower capability to regulate the stoma under mild water stress situations. Moreover, these findings corroborate with those revealed by García-Tejero et al. [11] who defined two threshold values for Ψ leaf ; the first about −1.4 MPa without reductions for g s , and a second of −2.0 MPa, when significant depletions in g s are observed. Also, other authors reported that under moderate water stress conditions the almond tree decreased the leaf rates much more than g s , and only when a certain threshold value is reached, significant depletions on photosynthesis rates are detected [36][37][38].

Nut Yield and Irrigation-Water Productivity (IWP)
Yield and its related components showed results that agreed with the physiological responses observed in the three studied cultivars (Table 2).
During the first experimental season, significant differences on kernel yield were observed in cv. Guara, with yield reductions around 13% in SDI 65 in comparison to FI treatment. These differences were exclusively reflected in fruits number per tree, with a fruit number depletion similar to that reported on total yield. Something similar happened with cv. Marta, with yield reductions around 11% in SDI 65 comparing to FI treatment, as consequence of fruits number reductions of 13.5% on average. However, these depletions were partially corrected by increasing the kernel unit weight, especially in SDI 75 . Moreover, no differences on kernel yield were observed for cv. Lauranne, with similar results among treatments in all the yield components. FI, Full irrigated treatment; SDI 75 , Sustained-deficit irrigation at 75% RI; SDI 65 , Sustained-deficit irrigation at 65% RI. Values followed by the same letter within the same row and factor were not significantly different (p < 0.05) by Tukey's test.
During the second year, cv. Guara showed again significant differences between SDI 65 and FI conditions, although in this case SDI 75 offered similar productions to those detected under FI; with yield reductions around 8% and 17% on SDI 75 and SDI 65 , respectively. These depletions were associated with fruit number reductions roughly 9% on SDI 75 and 32% on SDI 65 . It is noticeable that these reductions in the fruits number were partially corrected because of a significant increasing of kernel unit weight on SDI 65 (23% higher than the obtained value of FI). Something different were the obtained values for cv. Marta during the second year, which did not evidence relevant yield losses by effect of SDI strategies. Finally, as it was previously discussed for the first studied season, cv. Lauranne registered the best response to SDI strategies, without relevant differences in the studied yield components.
Analysing these results regarding to the irrigation water productivity (IWP, kg·m −3 ) (Table 3), in the first studied season and within each treatment, cvs. Guara and Marta evidenced similar results, these being 20%, 40%, and 31% higher in cv. Lauranne for FI, SDI 75 , and SDI 65 ; respectively. Regarding to the effects of water stress, all treatments fixed significant improvements in comparison to the results detected under FI. Thus, comparing the IWP obtained in SDI 65 with that registered under FI, the average values for both studied seasons offered improvements on IWP of 31%, 36%, and 43% for Guara, Marta, and Lauranne; respectively. These results were confirmed during the second season. All cultivars offered similar values of IWP for FI and SDI 75 , meanwhile, for the case of SDI 65 relevant improvements were detected in cvs. Marta and Lauranne in comparison to cv. Guara. Regarding to the effects of water stress, all cultivars registered significant increasing trend under SDI strategies, cvs. Lauranne and Marta offering the best response versus cv. Guara when sustained water withholdings of 35% were imposed. These values would be comparable to those reported by Egea et al. [39] for cv. Marta (0.25-0.40 kg·m −3 ) or Phogat et al. [18] who highlighted that water productivity increased substantially respect to full irrigated trees when SDI strategies were applied.  Taking into account the whole of data, cv. Guara was the most sensitive cultivar to water stress under sustained-deficit irrigation strategies. This fact is especially noticeable, taking into account that this cultivar evidenced a very positive response when water withholding was applied under RDI strategies during the kernel-filling period [26,28,32]. This contrary response would evidence that the final response to DI strategies would be determined by the added effect of three factors: The cultivar, the water stress, and the irrigation strategy. However, taking into consideration the findings in the present work, when this water stress was applied during the whole irrigation period, it was detected a significant fruit dropping (during the fruit setting period) ( Table 2), and ultimately, this determined the final yield with significant reductions linked to less fruit numbers per tree. These results agree with those reported by other authors. In this sense, for a sustained-deficit irrigation study (2500 m 3 ·ha −1 ), Alegre et al. [40] in Catalonia (North-Eastern Spain) reported productions for seven-year-old almond plantations of cvs. Guara and Lauranne of 1.65 and 2.02 t·ha −1 , respectively. The absolute difference with respect to our findings could mainly be ascribed to the amount of irrigation water applied in water-stressed treatments (4250 and 4730 m 3 ·ha −1 ) and the climatic conditions of South-Western Spain. Likewise, these results would confirm the better response of cv. Lauranne in comparison to cv. Guara under SDI strategies. This is in line with Miarnau et al. [41] who outlined that kernel yield of new almond plantations could be ranged between 1.50 and 2.0 t·ha −1 with water allocations of 2000-3000 and 6000 m 3 ·ha −1 , respectively. In general, the yield potential of almond is highly related to the irrigation amount provided as it was revealed by Miarnau et al. [42]. These authors pointed out that under a SDI strategy with water applications around 2000 m 3 ·ha −1 , the nut yield for cvs. Guara and Marta amounted to 1.20 and 1.85 t·ha −1 , whereas under FI conditions (≈7500 m 3 ·ha −1 ) these values were 2.80 and 3.55 t·ha −1 , respectively. Consequently, this fact would suggest that cv. Marta would be able to activate a physiological prevention mechanism to mitigate the water stress, yielding more than cv. Guara, as was observed in the present work. Moreover, similar agronomical and physiological responses to water stress of cvs. Guara and Lauranne were highlighted by Girona et al. [16].
According to effects of SDI in the fruit unit weight, Alegre et al. [40] reported similar values than those found in the present experiment with kernel unit weights of 1.5 and 1.2 g for cv. Guara and Lauranne. In addition, according to Miarnau et al. [43] kernel unit weight FI conditions for cvs. Guara, Lauranne, and Marta would be around 1.50, 1.20, and 1.50 g, respectively, implying that water stress provoked weight reduction as it was found in the present study. Likewise, it is worth mentioning the improvements in terms of fruit unit weight observed in cv. Marta under SDI 75 during the first experimental season; and something similar in cv. Guara under SDI 65 , during the second year of this experiment. These results would reinforce the possibility of improving the fruit size when SDI is imposed, this being an added value in relation to fruit marketability and consumer acceptance [44].

Conclusions
Combining the type of almond cultivars with water stress through deficit irrigation will be vital to reach an equilibrium between water allocations and sustainable nut yields under climate change scenarios. In the framework of the present experiment, the almond response to SDI strategies was cultivar-dependent, and hence, this fact should be considered before designing a proper DI strategy.
The findings allow for a conclusion on the importance of the cultivar when a DI strategy is being applied because different physiological behaviors will promote different responses in terms of yield and its components. In this way, the cv. Marta exhibited the most conservative behavior to water stress in physiological terms, which allowed it to obtain very similar productions than those registered under FI conditions. Furthermore, cv. Lauranne, despite showing a physiological behavior similar to cv. Guara, it was able to reach the best yield values when a moderate-to-severe SDI was applied. Also, according to our findings, cv. Guara registered the lesser promising results, with significant yield reductions (≈14%) when water restrictions around 35% of irrigation requirements were applied; these being particularly promoted by depletions in the fruit number per tree. That is, SDI 65 would be suitable strategies to cvs. Lauranne and Marta, whereas for the case of cv. Guara we should select a more moderate SDI strategy (as SDI 75 ) or re-consider the application of other more appropriate treatments such as RDI during the kernel-filling period.
Finally, taking into consideration the absence of differences in cvs. Lauranne and Marta, and the differing results observed in cv. Guara; long-term experiment could be advisable; in order to get a deeper knowledge respect to cumulative effects of more severe water stress strategies imposed during several consecutive seasons.