Deﬁcit Irrigation Strategies on Tree Physiological and Chemical Properties: Treatment Effects, Prediction Based Model Analyses and Inter-Correlations

: Irrigation is a key factor for different physiological aspects of fruit trees. Therefore, such irrigation protocols that can save water consumption during irrigation and maintain fruit trees productivity are an essential goal especially under semiarid climate conditions. The aim of this 3-year apricot study was to investigate the effect of four deﬁcit irrigation (DI) treatments (control, moderate regulated deﬁcit irrigation: RDIm, severe RDI: RDIs and continuous DI: CDI) on 15 tree physiological properties (chilling requirement—CR, heat requirement—HR, days from end—dormancy until fruit harvest—DEDFH, sum of growing degree days—sGDD, total number of buds—TNB, number of ﬂower buds—NFB, number of vegetative buds—NVB, starting date of ﬂowering—SDF, number of opened ﬂower buds—NOFB, ﬂower bud abscission—FBA, fruit set—FS, seasonal vegetative growth— SVG, fruit number per tree—FNT, fruit weight—FW, fruit yield—FY), and on two tree chemical properties (total soluble carbohydrates—TSC and total proline content—TPC) on apricot cultivars ‘Ninfa’ and ‘Canino’ in Egypt. Results showed that both DI treatments and cultivars signiﬁcantly inﬂuenced the values of CR, HR, TNB, SDF, NOFB, FS, SVG, FNT, FY, TSC, and TPC. Values of FBA were signiﬁcantly affected by years and DI treatments, while sGDD by years and cultivars. Values of DEDFH, NFB, and FW were signiﬁcantly inﬂuenced only by cultivars, while NVB only by DI treatments. The RDIm treatment gave the most acceptable values for most measured properties compared to the fully irrigated control treatment. Prediction based model analysis demonstrated that generalized linear models (GLMs) can be predictors for the measured tree properties in the DI treatments. The best goodness-of-ﬁt of the predicted GLMs was reached for HR, NOFB, FS, SVG, FNT, TSC, and TPC. In all the four DI treatments, 22 pair-variables (TNB versus (vs.) NFB, TNB vs. NOFB, TNB vs. NOFB, NFB vs. NOFB, NFB vs. FNT, NFB vs. FY, NFB vs. FW, NOFB vs. SVG, NOFB vs. FNT, NOFB vs. FY, FS vs. FNT, FS vs. FY, SVG vs. FNT, SVG vs. FY, SVG vs. TSC, FNT vs. FY, FY vs. FW, CR vs. TSC, HR vs. TNB, HR vs. NFB, HR vs. FNT, HR vs. FY, and NOFB vs. FBA) correlated signiﬁcantly in Pearson correlation and regression analyses. Principal component analyses explained 82% of the total variance and PC1, PC2, and PC3 explained 23, 21, and 15% of the total variance and correlated with the HR, TNB, FS, FNT and FY; FBA, SVG, TSC, and TPC; and NFB, NVB and NOFB, respectively, indicating strong connections among tree physiological and chemical


Introduction
Apricot (Prunus armeniaca L.) is one of the fundamental species of deciduous trees grown all over the world. Mediterranean countries produce more than 60% of the global apricot production [1]. Water is a natural resource and it is increasingly scarce and expensive especially in the Mediterranean regions. Therefore, water management is a key element for environmental sustainability especially in those countries that depends on agriculture activities and/or have limited water resources [2]. Limited water resources in Egypt are current and future challenges for Egyptian agriculture as about 13.5 BCM/year is the gap between the needs and availability of water [3]. Therefore, the optimization of irrigation strategies and improving the water use efficiency are hot research topics in the arid and semiarid regions of the world [4,5].
For improving the water use efficiency, a suggested irrigation strategy for fruit trees is the deficit irrigation (DI) practice especially in the arid and semiarid regions [4][5][6]. The DI practice is based on the knowledge of tree phenological stages and their responses to DI [2,7]. In order to study the benefits of DI, the crop water productivity (CWP) was proposed, which is the ratio of the yield and the consumed water by the crop [4]. Previous studies evaluated the effects of DI on tree vegetative growth, flowering, and fruit yield of various fruit tree species such as apples [7], peaches and pears [6,[8][9][10], plums [11,12], apricot [13][14][15][16][17][18][19], mango [20], and orange [21,22]. In these previous studies, the use of DI practices improved water use efficacy, maintained the fruit yield and decreased the cost of tree pruning, vegetative growth, as well as irrigation times and quantities. The sensitivity of apricot to DI depended upon the phenological stages of the trees [2]. The most sensitive phonological stages to DI were the rapid fruit growth (S3) and the early postharvest (S4) stages [2]. Water restrictions negatively affected fruit sizes during the S3 stage while the lack of water reduced the floral bud initiation and differentiation during the S4 stage [23,24]. On the other hand, Intrigliolo and Castel [11] found no effect of DI during pit-hardening (S2) stage and postharvest stage on flowering or fruit set. Therefore, a regulated deficit irrigation (RDI) was proposed during S2 stage of fruit trees e.g., [5,8,23]. Applications of an RDI during the intensive shoot growth stages reduced vegetative growth by decreasing shoot length and trunk diameter [7,11,17]. Therefore, a separation between shoot and fruit growth stages was proposed for applying the RDI strategies without or with minor negative effects on fruit yield [5,8,23]. As a consequence, a DI strategy has to reduce water use during the insensitive shoot growth stages while irrigating the trees well during the critical fruit growth periods. Although several studies evaluated DI treatments on fruit trees including apricot, an overall evaluation of various DI treatments on the measures of tree physiological and chemical properties coupled with prediction based model analyses, and the inter-correlations among these properties are limited under semiarid Mediterranean climate conditions. Providing the best inter-connections among tree properties and suitable prediction based models can help to clarify the background physiological processes of the apricot trees under DI conditions. Therefore, the aim of this 3-year study was to investigate the effect of four DI treatments (control, moderate RDI: RDIm, severe RDI: RDIs and continuous DI: CDI) on 15 measures of tree physiological properties (chilling requirement-CR, heat requirement-HR, days from end-dormancy until fruit harvest-DEDFH, sum of growing degree days-sGDD, total number of buds-TNB, number of flower buds-NFB, number of vegetative buds-NVB, starting date of flowering-SDF, number of opened flower buds-NOFB, flower bud abscission-FBA, fruit set-FS, seasonal vegetative growth-SVG, fruit number per tree-FNT, fruit weight-FW, fruit yield-FY), and on two tree chemical properties (total soluble carbohydrates-TSC and total proline content-TPC) for apricot cultivars Agronomy 2021, 11, 1361 3 of 33 (cvs) 'Ninfa' and 'Canino' under Egyptian climate conditions. Moreover, prediction based model analyses (Generalized Linear Models-GLMs) and inter-correlations (Pearson correlation, regression analyses and principal component analyses) were performed among the 17 measures of tree physiological and chemical properties in order to propose the best relationships of the correlated measures and provide suitable prediction based models under various DI practice.

Plant Material and Experimental Site
The study was conducted during the tree growing periods of three consecutive years from 2016 to 2018 on two apricot cvs 'Ninfa' and 'Canino' in a commercial apricot orchard at Wadi El Natrun (30.583 • N 30.333 • E), El Beheira Governorate, Northern Egypt ( Figure 1). The orchard is located under semiarid Mediterranean climate conditions and it is 23 and 38 m below sea and Nile River levels, respectively. The soil of the orchard is sandy loam and the soil physical and chemical characteristics are listed in Table 1. NOFB, flower bud abscission-FBA, fruit set-FS, seasonal vegetative growth-SVG, fruit number per tree-FNT, fruit weight-FW, fruit yield-FY), and on two tree chemical properties (total soluble carbohydrates-TSC and total proline content-TPC) for apricot cultivars (cvs) 'Ninfa' and 'Canino' under Egyptian climate conditions. Moreover, prediction based model analyses (Generalized Linear Models-GLMs) and inter-correlations (Pearson correlation, regression analyses and principal component analyses) were performed among the 17 measures of tree physiological and chemical properties in order to propose the best relationships of the correlated measures and provide suitable prediction based models under various DI practice.

Plant Material and Experimental Site
The study was conducted during the tree growing periods of three consecutive years from 2016 to 2018 on two apricot cvs 'Ninfa' and 'Canino' in a commercial apricot orchard at Wadi El Natrun (30.583° N 30.333° E), El Beheira Governorate, Northern Egypt ( Figure  1). The orchard is located under semiarid Mediterranean climate conditions and it is 23 and 38 m below sea and Nile River levels, respectively. The soil of the orchard is sandy loam and the soil physical and chemical characteristics are listed in Table 1. The orchard was established in 2009 and trees were grafted on apricot rootstocks and planted with a distance between and within rows of 4 and 4 m, respectively. The two examined apricot cultivars are the most commonly grown cultivars in the Egyptian region and they are reference cultivars for apricot growing. Mean tree height, mean trunk diameter, ground cover at the beginning and ground cover at the end of the study were 2.8 m, 17 cm, 47 and 55%, respectively. Eight trees replicated four times were selected for each cultivar and each irrigation treatment with nearly similar tree growth vigor and tree fruiting. The same trees were used for sampling during the three consecutive seasons. All treatments and both cultivars in all years received the same orchard management practices according to standard fruit farming practices. The trees had received 158 kg N, 769 kg P2O5 and 110 kg K2O ha −1 year −1 in all irrigation treatments and for both cultivars.  The orchard was established in 2009 and trees were grafted on apricot rootstocks and planted with a distance between and within rows of 4 and 4 m, respectively. The two examined apricot cultivars are the most commonly grown cultivars in the Egyptian region and they are reference cultivars for apricot growing. Mean tree height, mean trunk diameter, ground cover at the beginning and ground cover at the end of the study were 2.8 m, 17 cm, 47 and 55%, respectively. Eight trees replicated four times were selected for each cultivar and each irrigation treatment with nearly similar tree growth vigor and tree fruiting. The same trees were used for sampling during the three consecutive seasons. All treatments and both cultivars in all years received the same orchard management practices according to standard fruit farming practices. The trees had received 158 kg N, 769 kg P 2 O 5 and 110 kg K 2 O ha −1 year −1 in all irrigation treatments and for both cultivars.

Estimation of Irrigation Requirements for Apricot Tree
Irrigation water requirements for apricot trees were calculated according to data obtained from the local agrometeorological station at Al-Beharia Governorate, Central Laboratory for Agricultural Climate (C.L.A.C.), Ministry of Agriculture and Land Reclamation by using the following equation: where IR = irrigation water requirements (m 3 ha −1 day −1 ), ETc = the water requirement of a given crop (mm day −1 ), A = area irrigated (m 2 ), Ea = application efficiency (%) under 90% drip irrigation, CF = covering factor (for peach trees 45%), and LR = leaching requirements. The crop evapotranspiration (ETc) was calculated by using the following equation: where kc = crop coefficient factor. The kc varies among growth stages of the given crop. Normally four growth stages are considered: (i) the initial stage, when the crop uses little water; (ii) the crop development stage, when the water consumption increases; (iii) the mid-season stage, when water consumption reaches a peak; and (iv) the late-season stage, when the maturing crop once again requires less water. Crop coefficient factors of apricot were: 0.5 February, 0.75 March, 0.8 April, 0.9 May, 0.6 June, and 0.5 July-November [25].
ETo is the reference crop evapotranspiration (mm day −1 ). ETo was calculated as: where Epan = pan evaporation and Kpan = pan coefficient (values of Kpan vary between 0.35 and 0.85, with an average of 0.70).

Deficit Irrigation Treatments
For the DI treatments, four periods of phenological growth stages of apricot trees were determined from early February until late September according to the study of Ruiz-Sánchez et al. [23]. The four periods were: (i) stage 1 from early February to early March (swelling and opening of flower buds and flowering and fruit set), (ii) stage 2 from mid March to early May (fruit growth stages I and II), stage 3 from mid May to mid June (fruit growth stage III and fruit harvest), and stage 4 from end June to end September (fruit harvest and postharvest) ( Figure 2). All irrigation treatments were calculated up to the seasonal ETc as follows: (i) control, in which trees were irrigated at 100% of seasonal ETc for all phenological growth stages of the trees, (ii) moderate regular deficit irrigation (RDIm), in which trees were irrigated at 100% of ETc from mid-April to the end of July (from fruit growth stages II to fruit harvest), which defined as critical periods (CPs) according to the study of Ruiz-Sánchez et al. [23]; and trees were irrigated at 50% of ETc during the rest of the season (before mid-April and after the end of July), (iii) a severe regular deficit irrigation (RDIs), in which trees were irrigated at 100% of ETc during the CPs and trees Experimental design for the four deficit irrigation treatments (Control, RDIm, RDIs, and CDI) of apricot trees according to the months of the year, tree phenological phases (dormancy, flowering, fruit set, fruit growth, fruit harvest, postharvest), fruit growth stages (I., II., III.), critical periods and the percent of crop evapotranspiration (ETc) of irrigation treatments at Wadi El Natrun, El Beheira Governorate, Egypt (2016-2018). Control = fully irrigated scheme, RDIm = moderate regular deficit irrigation, RDIs = severe regular deficit irrigation and CDI = continuous deficit irrigation. Table 2. The monthly irrigated water (L tree −1 ) and overall annual irrigated water (L tree −1 , m 3 ha −1 ) for the four deficit irrigation treatments (control, RDIm, RDIs, and CDI) of apricot trees at Wadi El Natrun, El Beheira Governorate, Egypt (2016-2018). Explanation for Control, RDIm, RDIs, and CDI are given in Figure 2.

Control
RDIm RDIs CDI  January  150  75  37.5  75  February  950  475  237.5  475  March  1900  950  475  950  April  1900  1425  1187  950  May  1900  1900  1900  950  June  1900  1900  1900  950  July  1900  1900  1900  950  August  1900  950  475  950  September  1900  950  475  950  October  690  345  172.5  345  November  640  320  160  320  December  230  115   Experimental design for the four deficit irrigation treatments (Control, RDIm, RDIs, and CDI) of apricot trees according to the months of the year, tree phenological phases (dormancy, flowering, fruit set, fruit growth, fruit harvest, postharvest), fruit growth stages (I., II., III.), critical periods and the percent of crop evapotranspiration (ETc) of irrigation treatments at Wadi El Natrun, El Beheira Governorate, Egypt (2016-2018). Control = fully irrigated scheme, RDIm = moderate regular deficit irrigation, RDIs = severe regular deficit irrigation and CDI = continuous deficit irrigation. Table 2. The monthly irrigated water (L tree −1 ) and overall annual irrigated water (L tree −1 , m 3 ha −1 ) for the four deficit irrigation treatments (control, RDIm, RDIs, and CDI) of apricot trees at Wadi El Natrun, El Beheira Governorate, Egypt (2016-2018). Explanation for Control, RDIm, RDIs, and CDI are given in Figure 2. January  150  75  37.5  75  February  950  475  237.5  475  March  1900  950  475  950  April  1900  1425  1187  950  May  1900  1900  1900  950  June  1900  1900  1900  950  July  1900  1900  1900  950  August  1900  950  475  950  September  1900  950  475  950  October  690  345  172.5  345  November  640  320  160  320  December  230  115  Calculation of chilling requirement (CR) was based on chilling unit accumulation. According to Razavi et al. [26], a fixed starting date for chilling accumulation was set up. This date was a day in autumn (mid-October) when chilling accumulation become consistent and a negative effect of ambient temperature was scarce. From this date, 30 apricot branches (in three replicates, 10 branches per replicate) were taken randomly from similar positions of six randomly selected trees for each year, treatment and cultivar. The branches Agronomy 2021, 11, 1361 6 of 33 had a length of 40 cm and a diameter of 4.8-5.2 cm. Then each branch was immersed in a fungicide solution for 5 min and then wrapped with wetted papers then placed into plastic bags under cold conditions at temperature of 3 ± 1 • C. The artificial chilling treatments included 100, 200, 300, 400, 500, 600, and 700 h. After chilling treatments, the branches were placed from the proximal end in 0.5 L pots that were filled with 5% of sucrose solution and forced to grow in the laboratory at constant temperature 25 ± 1 • C and 65% RH for 16 h in light and 8 h in dark photoperiods. The water of pots was renewed every 3 days and the base of the branches was re-cut every 3 days. When 50% of the cuttings was sprouted, it was considered as the bud dormancy was released. Chilling requirements were defined by accumulating the chill hours below 7.2 • C according to the study of Weinberger [27] and Razavi et al. [26].

RDIm RDIs CDI
For evaluating heat requirement (HR), the growing degree hours (GDH) were estimated for each treatment and cultivar. The accumulation of GDH, as heat requirements, was considered under room temperature when 50% of the buds of the branches were sprouted after moving cooled branches from the refrigerator [28]. In addition, days from enddormancy until fruit harvest (DEDFH) was calculated from the date of the end-dormancy (when 5% of the flower buds opened) until the date of fruit harvest for each year, treatment and cultivar. The sum of growing degree days (sGDD) was also calculated during the period from the date of end-dormancy until the date of fruit harvest for each year, treatment and cultivar. The sGDD calculation was based on the equation for growing degree day (GDD) according to the study of Mohamed et al. [29]: where, T max = maximum temperature, T min = minimum temperature and T base is base temperature, 4.4 • C. Then sGDD was calculated by summing up each GDD from the date of the end-dormancy until the date of fruit harvest for each year, treatment and cultivar.

Bud Related Measures
In five trees, four × four branches were selected randomly in each tree and then labeled for each year, treatment and cultivar; and the total number of buds (TNB), the number of flower buds (NFB) and the number of vegetative buds (NVB) were counted. When 5% of the flower buds opened, the starting date of flowering (SDF) was determined for each year, treatment, and cultivar. When all flower buds opened, the same labeled branches were used to count the number of opened flower buds (NOFB) for each year, treatment and cultivar. After this, the percentage of flower bud abscission (FBA) was calculated as (TNB − NOFB) × 100 for each year, treatment and cultivar.

Measures for Yield Related Attributes
Fruit set (FS) was calculated as the percentage of the number of fruits to the number of flowers per branch. For SVG (cm), the vegetative growth of shoots, was measured with a veriner caliper at the date of end-dormancy (SVGed) and measured again at fruit harvest (SVDh) for each year, treatment, and cultivar. Then SVG was calculated as SVDh − SVGed. At harvest, the fruit number per tree (FNT) was counted for each selected tree of experiment for each year, treatment and cultivar. The individual fruit weight (FW, g) was measured for 4 × 50 fruit samples. Yield per tree (FY) was measured with a balance for each selected tree of experiment for each year, treatment and cultivar. For measuring the total soluble carbohydrate (TSC) content (i.e., sum of the glucose, fructose, and sucrose) of the tree was determined according to the method of Ranganna [30] with a slight modification. Apricot tree branches in four replicates were collected at 25, 27, and 24 November in 2016, 2017, and 2018, respectively. The collected samples were transported to laboratory and dried at 75 • C until constant weight. The dried samples were grinded and then powdered. A mixture of 0.1 g powdered sample and 13 mL ethanol (90%) were homogenized and then centrifuged at 4000 rpm for 15 min. The supernatant was used to measure the soluble sugar content. A mixture of 1 mL supernatant, 5 mL phenol (5%) and 5 mL sulfuric acid was vortexed for 30 s. Then the solution was adjusted to room temperature and the absorption of the solution was measured spectrophotometrically (BioTek Instruments Inc., VT 05404-0998, Winooski, VT, USA) at the wavelength of 490 nm. The standard curve was determined for a glucose standard solution. The total sugar content was measured as mg soluble carbohydrates g −1 dry weight.

Total Proline Content in Leaf
Total proline content (TPC) in apricot leaves was measured with the method described by the study of Bates and Teare [31]. Apricot leaves in four replicates were collected on 30, 28, and 29 September in 2016, 2017 and 2018, respectively, and frozen immediately and stored at −20 • C. A mixture of frozen leaves (1 g) and sulfosalicylic acid (3%) was homogenized then the tubes were centrifuged at 12,000× g for 10 min. 1 mL of supernatant was mixed with 1 mL of 1:1 (v/v) ratio of glacial acetic acid and ninhydrin reagent. The mixture solution was placed in water bath at 100 • C for 1 h. Then the samples were extracted with toluene. The absorption of the solution was measured by a spectrophotometer (BioTek Instruments Inc., VT 05404-0998, Winooski, VT, USA) at 520nm wavelength. The proline concentration was determined from a standard curve and the proline content was expressed as (µg g −1 ) of fresh matter.
2.6. Statistical Analysis 2.6.1. ANOVA A randomized complete block design (RCBD) was used to prepare the treatments. ANOVA was performed to analyze data set using an SPSS program (SPSS Inc., Chicago, IL, USA). The effects of DI treatments (control, RDIm, RDIs, and CDI), cultivars ('Ninfa' and 'Canino'), growing season i.e., year (2016, 2017, and 2018) and their two-and three-way interactions were evaluated on all the 17 measures of tree physiological and chemical properties. Duncan's multiple range tests were used to separate means at P = 0.05 levels.

Prediction Based Model Analysis
GLMs (Generalized Linear Models) were conducted, as a 3-way factorial ANOVA [32]. We used the models to predict all the 17 measures of tree physiological and chemical properties as the dependent variables (biological parameters), the independent variables were the factorial variables of irrigation treatment, cultivar and growing season i.e., year. A model for each measure (i.e., for all the 17 measures) was determined. The dataset were separated into a training and a testing part in 50-50%. Models were trained with the k-fold cross-validation method: we applied five folds with five repetitions using the training data. k-fold cross-validation starts with splitting randomly the input data into k-folds; in our case it was a split into five folds where four folds were used to train the models and one was used for testing. In the next step, other four folds were considered as training data and another one for testing. The procedure stops when all folds had been used as testing data. The whole procedure was repeated five times with new randomly selected folds [33,34]. Thus, we were able to exclude the possibility to gain a favorable or an inappropriate model by chance, we had 25 repetitions, with 25 R 2 , Mean Absolute Error (MAE), and Root Mean Square Error (RMSE). The great advantage of the k-fold cross-validation is that all models are applied on independent data (i.e., on a fold which was not included in the training), and the several model solutions provide a feedback on data reliability: if ranges of model fit parameters are in a narrow range, all models tend to justify the usability of a model.
Then, we applied the model on the independent part of the dataset, i.e., testing data. Predictions (modelled data) were compared to the observed data, and differences were plotted on the combination of a violin plot vs. boxplot using the factors of cultivars, irrigation deficit, and the years of measurements. This approach ensured to report effectiveness of the factors as predictors. ANOVA usually reveals the efficiency of a factor with justifying the H0: group means do not differ significantly. Our approach used the model and the result showed, beside the group differences, more importantly, how accurate predictions can be gained with the factors. If a GLM model was accurate, differences between observed and predicted values were around 0, while large deviations indicated issues in the predictors. The satisfactory model efficiency was evaluated using the relationship between RMSE of the models and the standard deviation (SD) of the observed data. Then the criterion for goodness-of-fit of the models was determined by the calculation of Nash-Sutcliffe efficiency (NSE) coefficient [35]. An NSE coefficient value above 0.65 indicated acceptable model fit to the observed data.

Correlation and Linear Regression Analyses among Parameters
In order to quantify relationships among the 17 measures of the tree physiological and chemical properties, Pearson's correlation coefficients were determined for the relationships of the measures in all combinations (136 variable pairs). Correlation analyses were done for the total data of DI treatments and also separately for each DI treatment. Pearson's correlation analyses were performed by using Genstat 5 Release 4.1 (Lawes Agricultural Trust, IACR, Rothamsted, UK). Then, the best-correlated variable pairs were further analyzed by linear regression analyses using the equation of f(x) = ax + b. Then, a t-test was used for comparing regression slopes in order to quantify the differences among the four DI treatments (control, RDIm, RDIs, and CDI) at P = 0.05.

Principal Component Analysis
A standardized Principal Component Analysis (PCA) was prepared based on the correlation matrix with the CR, HR, DEDFH, sGDD, TNB, NFB, NVB, SDF, NOFB, FBA, FS, SVG, FNT, FW, FY, TSC, and TPC variables. A standardization of all variables was prepared transforming the values to z-scores. Model fit had been tested with the Root Mean Square Residual (RMSR) [39]. Principal Components (PCs) were plotted in biplot diagrams. PCA was performed in R 4.04 [36] with the psych [40], FactoMiner [41] and factoextra [42] packages.

General Measures for Tree Physiological Properties
Analysis of variance (ANOVA) on CR showed that there was a significant (P < 0.05) effect for irrigation treatments and cultivars, while the year variable was non-significant. There were no significant two-or three-way interactions for the treatments' effects (Table 3). According to ANOVA, the data were presented as the average of 3 years ( Figure 3A). Table 3. The monthly irrigated water (L tree −1 ) and overall annual irrigated water (L tree −1 , m 3 ha −1 ) for the four deficit irrigation treatments (control, RDIm, RDIs, and CDI) of apricot trees at Wadi El Natrun, El Beheira Governorate, Egypt (2016-2018). Bold figures indicate significant differences at P < 0.05. Explanations for Control, RDIm, RDIs and CDI are given in Figure 2.    Control means that trees were irrigated at 100% of seasonal ETc for all phenological growth stages of the trees. Differences among control, RDIm, RDIs and CDI treatments were represented by LSD 0.05 values at P = 0.05, separately for the two cultivars ('Ninfa' and 'Canino'). After each LSD 0.05 value, first, second, third, and fourth letters belongs to the control, RDIm, RDIs and CDI treatments, respectively. Values within the given cultivar coupled with different letters are significantly different among the control, RDIm, RDIs and CDI treatments at P = 0.05 according to LSD t-tests.

Source of
The lowest CR value was 307.9 h in the control irrigation treatment for cv. 'Canino', while the highest one was 391.8 h in the CDI treatments for cv. 'Ninfa' ( Figure 3A). Cultivar 'Canino' showed significantly lower mean CR values compared to cv. 'Ninfa' (means were 326.06 and 364.4 h, respectively). The two cultivars showed different effect of the DI treatments on CR. Treatments of RDIm and RDIs resulted in significant reduction in CR compared to control or CDI treatments for cv. 'Ninfa'. On the other hand, treatments of RDIs and CDI resulted in significant increase in CR compared to control or RDIm treatments for cv. 'Canino'. For this cultivar, treatments of RDIs and CDI as well as treatments of control and RDIm did not differ from each other at P < 0.05.
ANOVA on HR showed that there was a significant (P < 0.05) effect for irrigation treatment and cultivar, while the year variable was non-significant. There were no significant two-or three-way interactions for the treatments' effects (Table 3). According to ANOVA, the data were presented as the average of 3 years ( Figure 3B). The lowest HR value was 4542.2 h in the CDI treatment for cv. 'Ninfa', while the highest one was 5518.9 h in the RDIm treatment for cv. 'Canino' ( Figure 3B). Cultivar 'Ninfa' showed significantly lower mean HR values compared to cv. 'Canino' (means were 4623.4 and 5418.2 h, respectively). The two cultivars showed different effect of the DI treatments on HR. Treatments of RDIm, RDIs and CDI resulted in significant reduction in HR compared to control for cv. 'Ninfa' but the three DI treatments did not differ significantly from each other. Treatments of RDIm and CDI resulted in significant increase in HR compared to control or RDIs treatments for cv. 'Canino'. For this cultivar, treatments of RDIm and CDI as well as treatments of control and RDIs did not differ from each other at P < 0.05.
ANOVA on DEDFH showed that there was a significant (P < 0.05) effect for cultivar, while the irrigation treatment and the year variable were non-significant. There were no significant two-or three-way interactions for the treatments' effects (Table 3). Therefore, data were discussed only for the two cultivars as the average of 3 years and four irrigation treatments. The DEDEH value of cv. 'Ninfa' (63.6 days) was significantly lower than that of cv. 'Canino' (71.1 days).
ANOVA on sGDD showed that there were significant (P < 0.05) effects for year and cultivar, while the irrigation treatment was non-significant. There were no significant twoor three-way interactions for the treatments' effects (Table 3). According to ANOVA, the data were presented as the average of the four irrigation treatments (

Bud Related Measures for Tree Physiological Properties
ANOVA on TNB showed that there was a significant (P < 0.05) effect for irrigation treatment and cultivar, while the year variable was non-significant (Table 4). There were no significant two-or three-way interactions for the treatments' effects. According to ANOVA, the data were presented as the average of 3 years for TNB ( Figure 5A). The lowest TNB value was 104.6 buds in CDI treatment for cv. 'Canino', while the highest one was 128.9 buds in the control treatments for cv. 'Ninfa' ( Figure 5A). Cultivar 'Ninfa' showed significantly higher mean TNB values compared to cv. 'Canino' (means were 120.9 and 109.8 buds, respectively). The two cultivars showed different effect of the DI treatments on TNB. Treatments of RDIm, RDIs and CDI resulted in significant reduction in TNB compared to control treatments for cv. 'Ninfa' but the three DI treatments did not differ significantly from each other. Treatments of RDIs resulted in a significant increase in TNB compared to control or CDI treatments for cv. 'Canino'. or three-way interactions for the treatments' effects (Table 3). According to ANOVA, the data were presented as the average of the four irrigation treatments (Figure 4). The lowest sGDD value was 1101.6 days in 2016 for cv. 'Ninfa', while the highest one was 1368.6 days in 2018 for cv. 'Canino' (Figure 4). Cultivar 'Ninfa' showed significantly lower mean sGDD values compared to cv. 'Canino' (means were 1170.5 and 1315.8 days, respectively). Values of sGDD was significantly higher in 2017 and 2018 than in 2016 for cv. 'Ninfa' while the years of 2018 and 2016 differed significantly for cv. 'Canino'.    Table 3.   Number of opened flower buds (NOFB) ))))))))  ANOVA on NFB showed a significant effect for cultivar, while the year variable and irrigation treatment were non-significant (Table 4). There were no significant two-or three-way interactions for the treatments' effects. Therefore, data were discussed only for the two cultivars as the average of 3 years and four irrigation treatments. The NFB value of cv. 'Ninfa' (86.3 buds) was significantly higher than that of cv. 'Canino' (77.9 buds).

Source of
ANOVA on NVB showed a significant effect for irrigation treatment, while the year variable and cultivar were non-significant (Table 4). There were no significant two-or threeway interactions for the treatments' effects. According to ANOVA, the data were presented as the average of 3 years and two cultivars for NVB ( Figure 5B). The NVB value of CDI treatment (29.3 buds) was significantly lower than that of the RDIs or control treatments (34.0 and 36.7 buds, respectively).
ANOVA on SDF showed a significant effect for year, irrigation treatment and cultivar (Table 4). There were no significant two-or three-way interactions for the treatments' effects. Therefore, SDF data set were shown separately for years, cultivars, and irrigation treatments ( Table 5). The lowest SDF value was 66.0 in RDIm treatment for cv. 'Canino' in 2017, while the highest one was 79.0 in the RDIs treatments for cv. 'Canino' in 2016 (Table 5). Mean SDF values for cultivars showed that cv. 'Ninfa' showed significantly higher mean SDF values in 2017 compared to cv. 'Canino' (73.5 and 67.9, respectively). The significance of SDF values within irrigation treatments greatly varied among years and cultivar. Consistent significant differences (at P < 0.05) among SDF values were found between RDIm and CDI treatments in the overall year analyses for cv. 'Ninfa' but this effect was significant only in 2017 for cv. 'Canino'. Table 5. The effects of three years (2016, 2017, 2018), deficit irrigation treatments (control, moderate regular deficit irrigation-RDIm, severe regular deficit irrigation-RDIs and continuous deficit irrigation-CDI) and two cultivars ('Ninfa' and 'Canino') on starting date of flowering (SDF) in a commercial apricot orchard at Wadi El Natrun, Northern Egypt. Explanation for irrigation treatments, calculations of significant differences, and explanations for LSD 0.05 values are given in Figure 3; ns: nonsignificant. Different letters represent significant differences among DI treatments. ANOVA on NOFB showed that there was a significant (P < 0.05) effect for irrigation treatment and cultivar, while the year variable was non-significant (Table 4). There were no significant two-or three-way interactions for the treatments' effects. According to ANOVA, the data were presented as the average of 3 years for NOFB ( Figure 5C). The lowest NOFB value was 52.3 buds in CDI treatment for cv. 'Canino', while the highest one was 68.3 buds in the control treatments for cv. 'Ninfa' (Figure 5C). Cultivar 'Ninfa' showed significantly higher mean NOFB values compared to cv. 'Canino' (means were 64.1 and 56.4 buds, respectively). The two cultivars showed different effect of the DI treatments on NOFB. Treatments of RDIs and CDI resulted in significant reduction in NOFB compared to control treatments for both cultivars but these two DI treatments did not differ significantly from each other.

Cultivar/Irrigation Treatment
ANOVA on FBA showed a significant effect for the year variable and irrigation treatment, while cultivar was non-significant. Significant two-or three-way interactions were not detected for the treatments' effects (Table 4). According to ANOVA, FBA data were presented as the average of the two cultivars ( Figure 5D). The lowest FBA value was 18.9% in the control treatment in 2016, while the highest one was 33.9% in the CDI treatments in 2017 ( Figure 5D). Mean FBA values in 2017 was significantly higher than FBA values in 2016 or 2018 (means were 24.9, 28.1, and 24.5, for 2016, 2017, and 2018, respectively). In general, DI treatments increased FBA values compared to the control treatment but the treatment of CDI resulted in a significant increase in all years compared to the control treatment.

Yield Related Measures for Tree Physiological Properties
ANOVA on FS, SVG, FNT, and FY showed that there were significant (P < 0.05) effects for irrigation treatment and cultivar, while the year variable was non-significant (Table 6). There were no significant two-or three-way interactions for the treatments' effects. According to ANOVA, the data of FS, SVG, FN, and FY were presented as the average of 3 years for FS ( Figure 6A-D).
The lowest mean FS value was 26.1% in CDI treatment for cv. 'Canino', while the highest one was 39.7% in the RDIm treatments for cv. 'Ninfa' ( Figure 6A). Cultivar 'Ninfa' showed significantly higher mean FS values compared to cv. 'Canino' (means were 38.2 and 29.4%, respectively). The two cultivars showed different effect of the DI treatments on FS. Treatments of RDIm and RDIs resulted in a significant increase in FS values compared to control or CDI treatments, and FS values were significantly higher in the RDIm treatment than RDIs one for cv. 'Ninfa'. FS values in RDIs and CDI treatments were significantly lower than in the control or RDIm treatments for cv. 'Canino'. Table 6. Analyses of variance for the effects of year, deficit irrigation treatment and cultivar on yield related measures of tree physiological properties (fruit set-FS, seasonal vegetative growth-SVG, fruit number per tree-FNT, fruit weight-FW, fruit yield-FY) in a commercial apricot orchard in Egypt (2016-2018, cvs 'Ninfa' and 'Canino'). Explanations for bold figures are given in Table 3.  Seasonal vegetaive growth (SVG))))))))) Fruit number per tree (FNT) ))  The lowest mean SVG value was 23.4 cm in CDI treatment for cv. 'Canino', while the highest one was 43.0 cm in the control treatments for cv. 'Ninfa' (Figure 6B). Cultivar 'Ninfa' showed significantly higher mean SVG values compared to cv. 'Canino' (means were 37.0 and 31.5 cm, respectively). The two cultivars showed different effects of the DI treatments on SVG. Treatments of RDIm, RDIs, and control resulted in a significant increase in SVG compared to the CDI treatment for cv. 'Ninfa', but SVG values of the three treatments were not significantly different ( Figure 6B). Values of SVG in the RDIs and RDIm treatments were significantly higher compared to the CDI treatment but were significantly lower compared to the control treatment for cv. 'Canino'.

Source of
The lowest mean FNT value was 138.8 in CDI treatment for cv. 'Canino', while the highest one was 367.7 in the RDIm treatments for cv. 'Ninfa' (Figure 6C). Cultivar 'Ninfa' showed significantly higher mean FNT values compared to cv. 'Canino' (means were 265.3 and 182.3, respectively). FNT values of all DI treatments differed from each other at P = 0.05 for both cultivars. The lowest FNT values were obtained for the CDI treatments for both cultivars. All other DI treatments were significantly higher compared to CDI treatment for both cultivars.
The lowest mean FY value was 5.0 kg tree −1 in CDI treatment for cv. 'Canino', while the highest one was 13.5 in the RDIm treatments for cv. 'Ninfa' ( Figure 6D). Cultivar 'Ninfa' showed significantly higher mean FY values compared to cv. 'Canino' (means were 10.0 and 6.4, respectively). The two cultivars showed different effect of the DI treatments on FY. The RDIm treatment resulted in significantly higher FY values compared to RDIs, CDI, and control treatments for cv. 'Ninfa', but FY values of these three treatments were not significantly different. Values of FY in the fully irrigated control treatments were significantly higher compared to FY values of RDIm, RDIs, and CDI treatments for cv. 'Canino' but FY values of these three treatments were not significantly different.
ANOVA on FW showed that there was a significant (P < 0.05) effect for cultivar, while the year and irrigation treatment were non-significant. There were no significant two-or three-way interactions for the treatments' effects (Table 6). Therefore, data were discussed only for the two cultivars as the average of 3 years and 4 irrigation treatments. The mean FW value of cv. 'Ninfa' (38.3 g) was significantly higher than that of cv. 'Canino' (34.9 g).

Measures of Tree Chemical Properties
ANOVA on TSC showed a significant effect for year, irrigation treatment and cultivar (Table 7). Table 7. Analyses of variance for the effects of year, deficit irrigation treatment, and cultivar on yield related measures of tree physiological properties (fruit set-FS, seasonal vegetative growth-SVG, fruit number per tree-FNT, fruit weight-FW, fruit yield-FY) in a commercial apricot orchard at Wadi El Natrun, Northern Egypt over 2016-2018 on two apricot cultivars 'Ninfa' and 'Canino'. Explanations for bold figures are given in Table 3. There were no significant two-or three-way interactions for the treatments' effects. Therefore, TSC data set were shown separately for years, cultivars and irrigation treatments ( Table 8). The lowest TSC value was 109.4 mg g −1 in CDI treatment for cv. 'Canino' in 2017, while the highest one was 153.0 mg g −1 in the control treatment for cv. 'Ninfa' in 2018 (Table 8). Mean TSC values for cultivars showed that cv. 'Ninfa' showed significantly higher mean TSC values compared to cv. 'Canino' in all years (overall mean values were 132.5 and 124.7 mg g −1 , respectively). TSC values of RDIm, RDIs and CDI treatments were significantly lower for both cultivars compared to the control treatment in all years. TSC values of the RDIm treatment were significantly higher for cv. 'Ninfa' compared to the RDIs and CDI treatments in all years. TSC values of RDIm, RDIs, and CDI treatments were not significantly different from each other for cv. 'Canino' in the three years. Table 8. The effects of three years (2016, 2017, 2018), deficit irrigation treatments (control, moderate regular deficit irrigation-RDIm, severe regular deficit irrigation-RDIs and continuous deficit irrigation-CDI) and two cultivars ('Ninfa' and 'Canino') on the soluble carbohydrates content (TSC, mg g −1 ) in a commercial apricot orchard at Wadi El Natrun, Northern Egypt. Explanation for irrigation treatments, calculations of significant differences, and explanations for LSD 0.05 values are given in Figure 3; ns: nonsignificant. Different letters represent significant differences among DI treatments. ANOVA on TPC showed that there were significant (P < 0.05) effects for irrigation treatment and cultivar, while year was non-significant (Table 7). There were no significant two-or three-way interactions for the treatments' effects. According to ANOVA, the data of TPC were presented as the average of 3 years for TPC (Figure 7). The lowest mean TPC value was 26.3 µg g −1 in control treatment for cv. 'Ninfa', while the highest one was 67.4 µg g −1 in the CDI treatments for cv. 'Ninfa' (Figure 7). Cultivar 'Ninfa' showed significantly higher mean TPC values compared to cv. 'Canino' (means were 43.7 and 39.6 µg g −1 , respectively). Treatments of RDIm, RDIs and CDI resulted in significant increase in TPC compared to the control treatment for cv. 'Ninfa', and TPC values were significantly different in all the four irrigation treatments. TPC values in RDIs and CDI treatments were significantly higher than in control or RDIm treatments for cv. 'Canino'. TPC values in the control or RDIm treatments were not significantly different for cv. 'Canino'.

Prediction Based Model Analysis
Prediction based model analysis showed that the differences of the GLM predicted (violin plot) vs. the observed (boxplot) data ( Figure 8  ))))))))

Prediction Based Model Analysis
Prediction based model analysis showed that the differences of the GLM predicted (violin plot) vs. the observed (boxplot) data ( Figure 8) were larger than 0 for all the 17 tree measures and for all the three factors (years, irrigation treatments, and cultivars). Thus the GLM models slightly underestimated the real (observed) values for the 17 measures (median RMSEs were 23.39, 185.27, 7.77, 62.7, 9.04, 8.81, 9.14, 3.76, 6.24, 5.18, 7.03, 3.74, 43.2, 5.63, 2.49, 8.12, and 5.32, for CR, HR, DEDFH, sGGD, TNB, NFB, NVB, SDF, NOFB, FBA, FS, SVG, FNT, FW, FY, TSC, and TPC, respectively). However, the efficiency of most models was satisfactory based on the relationship between RMSE of the models and SD of the observed data. According to the calculated NSE coefficient values (ranging between 0.65 and 0.9), the best goodness-of-fit of the predicted GLM models was reached for the measures of HR, NOFB, FS, SVG, FNT, TSC, and TPC.

Pearson Correlation Analyses
When all data set for irrigation treatments, cultivars and years were analyzed together, 136 variable pairs were correlated from which 29 pair-variables correlated significantly at P = 0.05 probability level (Table 9A). When data set of 4 × 136 parameter pairs were analyzed separately for each of the four irrigation treatments, 45, 55, 53, and 43 pairvariables correlated significantly in the control, RDIm, RDIs, and CDI irrigation treatments, respectively (Table 9B,C,D,E). Table 9. Pearson's correlation coefficients (r) amongst 17 measures for four deficit irrigation treatments (control-Cont., moderate regular deficit irrigation-RDIm, severe regular deficit irrigation-RDIs and continuous deficit irrigation-CDI) in a commercial apricot orchard at Wadi El Natrun, Northern Egypt over 2016-2018, on two apricot cultivars 'Ninfa' and 'Canino'. Bold figures represent the significant (P < 0.05) correlation coefficient values in the deficit irrigation treatments. CR: chilling requirements, HR: heat requirements, DEDFH: days from end-dormancy until fruit harvest, sGDD: sum of growing degree days, TNB: total number of buds, NFB: number flower buds, NVB: number of vegetative buds, SDF: starting date of flowering, NOFB: number of opened flower buds, FBA: flower bud abscission (%), FS: fruit set, SVG: seasonal vegetative growth, FNT: fruit number per tree, FY: fruit yield, FW: fruit weight, TSC: total soluble carbohydrates and TPC: total proline content in leaf. A: Total = all data were used and combined for years, cultivars and irrigation treatments. B: Cont. = data from control deficit irrigation treatment were used for calculating correlations; data were combined for years and cultivars. C: RDIm = data from RDIm treatment were used for calculating correlations; data were combined for years and cultivars. D: RDIs = data from RDIs treatment were used for calculating correlations; data were combined for years and cultivars. E: CDI = data from CDI treatment were used for calculating correlations; data were combined for years and cultivars.

Pearson Correlation Analyses
When all data set for irrigation treatments, cultivars and years were analyzed together, 136 variable pairs were correlated from which 29 pair-variables correlated significantly at P = 0.05 probability level (Table 9A). When data set of 4 × 136 parameter pairs were analyzed separately for each of the four irrigation treatments, 45, 55, 53, and 43 pairvariables correlated significantly in the control, RDIm, RDIs, and CDI irrigation treatments, respectively (Table 9B,C,D,E). Table 9. Pearson's correlation coefficients (r) amongst 17 measures for four deficit irrigation treatments (control-Cont., moderate regular deficit irrigation-RDIm, severe regular deficit irrigation-RDIs and continuous deficit irrigation-CDI) in a commercial apricot orchard at Wadi El Natrun, Northern Egypt over 2016-2018, on two apricot cultivars 'Ninfa' and 'Canino'. Bold figures represent the significant (P < 0.05) correlation coefficient values in the deficit irrigation treatments. CR: chilling requirements, HR: heat requirements, DEDFH: days from end-dormancy until fruit harvest, sGDD: sum of growing degree days, TNB: total number of buds, NFB: number flower buds, NVB: number of vegetative buds, SDF: starting date of flowering, NOFB: number of opened flower buds, FBA: flower bud abscission (%), FS: fruit set, SVG: seasonal vegetative growth, FNT: fruit number per tree, FY: fruit yield, FW: fruit weight, TSC: total soluble carbohydrates and TPC: total proline content in leaf. A: Total = all data were used and combined for years, cultivars and irrigation treatments. B: Cont. = data from control deficit irrigation treatment were used for calculating correlations; data were combined for years and cultivars. C: RDIm = data from RDIm treatment were used for calculating correlations; data were combined for years and cultivars. D: RDIs = data from RDIs treatment were used for calculating correlations; data were combined for years and cultivars. E: CDI = data from CDI treatment were used for calculating correlations; data were combined for years and cultivars.     (Table 9).

Linear Regression Analyses
The significant 22 pair-variables were further presented by linear regression analysis, showing the data of the four DI treatments (control, RDIm, RDIs, and CDI) with separated symbols (Figure 9). All the 22 pair-variables showed significant linear relationships with r values ranging from 0.689 to 0.903, and with P values ranging from 0.049 to 0.001. However, no differences were observed among the slope parameters for all pair-variables among the control, RDIm, RDIs, and CDI treatments as t-tests showed P values ranging from 0.847 to 0.059.

Principal Component Analyses
PCA explained the 82% of the total variance for the treatments (Figure 10A (Table 9).

Linear Regression Analyses
The significant 22 pair-variables were further presented by linear regression analysis, showing the data of the four DI treatments (control, RDIm, RDIs, and CDI) with separated symbols (Figure 9). All the 22 pair-variables showed significant linear relationships with r values ranging from 0.689 to 0.903, and with P values ranging from 0.049 to 0.001. However, no differences were observed among the slope parameters for all pair-variables among the control, RDIm, RDIs, and CDI treatments as t-tests showed P values ranging from 0.847 to 0.059.

Principal Component Analyses
PCA explained the 82% of the total variance for the treatments ( Figure 10A,B). Number of PCs had been justified by the RMSR: 0.06 (5 PCs), indicating a good fit.
PC1 explained 23% of the total variance and correlated with the HR, TNB, FS, FNT, and FY ( Figure 10A,B). PC2 explained 21% of the total variance and correlated with the FBA, SVG, TSC, and TPC. PC3 explained 15% of the total variance and correlated with the NFB, NVB, and NOFB. PC4 explained 13% of the total variance and correlated with the CR, TNB, and NVB. In addition, 10% of the variance was accounted for PC5 and correlated with the sGDD and SDF.
The biplot figure for DI treatments showed that both PC1 and PC2 axes were dominant for all the four DI treatments (control, RDIm, RDIs, and CDI) ( Figure 10A). The biplot figure showed that PC2 axes were dominant for CDI treatment while PC1 axes were the most dominant for RDIm treatment ( Figure 10A).
The biplot figure for cultivars showed that both PC1 and PC2 axes were dominant for the two cultivars ('Ninfa' and 'Canino') ( Figure 10B). The biplot figure showed that PC1 axes were more dominant for cv. 'Canino' while PC2 axes were more dominant for cv. 'Ninfa' ( Figure 10B).
The decisive role of these two PC groups was also indicated by long arrows' length and proximity of the variables of HR, TSC, SVG, FNT, FY, TPC, and FBA. These results indicated strong associations among tree physiological and chemical properties ( Figure  10A,B).

Discussion
This study reported the various effects of different DI strategies on tree physiological and chemical attributes on two commonly grown apricot cultivars ('Ninfa' and 'Canino') in a semiarid Mediterranean region under sandy loam soil conditions in order to promote possible solutions to the water scarcity problem in apricot production [3].
In this study, ANOVA showed that CR were significantly affected by DI treatment and apricot cultivars under semiarid Mediterranean conditions in Egypt (Table 3) and CR increased by increasing water deficit ( Figure 3A). This result was in agreement with the PC1 explained 23% of the total variance and correlated with the HR, TNB, FS, FNT, and FY ( Figure 10A,B). PC2 explained 21% of the total variance and correlated with the FBA, SVG, TSC, and TPC. PC3 explained 15% of the total variance and correlated with the NFB, NVB, and NOFB. PC4 explained 13% of the total variance and correlated with the CR, TNB, and NVB. In addition, 10% of the variance was accounted for PC5 and correlated with the sGDD and SDF.
The biplot figure for DI treatments showed that both PC1 and PC2 axes were dominant for all the four DI treatments (control, RDIm, RDIs, and CDI) ( Figure 10A). The biplot figure showed that PC2 axes were dominant for CDI treatment while PC1 axes were the most dominant for RDIm treatment ( Figure 10A).
The biplot figure for cultivars showed that both PC1 and PC2 axes were dominant for the two cultivars ('Ninfa' and 'Canino') ( Figure 10B). The biplot figure showed that PC1 axes were more dominant for cv. 'Canino' while PC2 axes were more dominant for cv. 'Ninfa' (Figure 10B).
The decisive role of these two PC groups was also indicated by long arrows' length and proximity of the variables of HR, TSC, SVG, FNT, FY, TPC, and FBA. These results indicated strong associations among tree physiological and chemical properties ( Figure 10A,B).

Discussion
This study reported the various effects of different DI strategies on tree physiological and chemical attributes on two commonly grown apricot cultivars ('Ninfa' and 'Canino') in a semiarid Mediterranean region under sandy loam soil conditions in order to promote possible solutions to the water scarcity problem in apricot production [3].
In this study, ANOVA showed that CR were significantly affected by DI treatment and apricot cultivars under semiarid Mediterranean conditions in Egypt (Table 3) and CR increased by increasing water deficit ( Figure 3A). This result was in agreement with the study of Nasrabadi et al. [43] who showed that CR was significantly affected by DI treatments on pomegranate trees. The authors revealed that water deficit in the soil resulted in an increase of chilling requirement; but the level of response varied by cultivars [43]. This cultivar difference under DI treatments was also detected in this study for apricot ( Figure 3A). All these results showed that water deficit increases CR for apricot and pomegranates. This indicates that the required duration for winter cold temperature for completing CR has to be longer under increasing water deficit conditions in order to ensure regular flowering, fruit set, and fruit yield of the trees. This result is a strong practical indication for temperate fruit tree production of the Mediterranean regions as negative effects of the insufficient rates of chilling accumulation can be more severe in the Mediterranean regions, such as in Egypt compared to more cool regions of temperate fruit tree production [44,45]. Our results also indicated that growers in the Mediterranean regions have to choose an apricot cultivar with a lower CR if a water DI scheme is planned to be applied during fruit production. For an adequate cultivar selection, previous apricot studies [26,28,46,47] including this study, reported the CR values (between 322 and 1266 chill units) of apricot cultivars under Mediterranean climate conditions.
In this study, ANOVA showed that HR was significantly affected by DI treatment and apricot cultivars under Mediterranean climate conditions in Egypt (Table 3) but HR had no consistent increase or decrease with increasing water deficit ( Figure 3B). HR was not significantly affected by different DI treatments for the pomegranate tree in a previous study of Nasrabadi et al. [43], while the cultivars showed significantly different heat requirements, which coincided with this study too (Table 3 and Figure 3B). The differences among the HR values of the apricot cultivars may be genetically determined as a wide range of HR values was reported for various apricot genotypes in previous studies [26,28,46].
In this study, ANOVA showed that irrigation treatments had significant effect on TNB and NVB ( Table 4) which indicated that these bud measures were sensitive to DI treatments. However, increasing water deficit showed no clear decrease or increase in the bud numbers ( Figure 5A,B). For instance, values of the TNB significantly decreased in the DI treatments of RDIs compared to fully irrigated control treatment for cv. 'Ninfa', while TNB significantly increased in the RDIs treatments compared to control for cv. 'Canino' (Figure 5A). Values of NVB in the RDIs treatments increased again compared to the control one ( Figure 5B). However, the number of flower buds (NFB) was insensitive to DI treatments (Table 4). This latter result coincide with pervious results of Demirtas et al. [48] and Blanco et al. [49] in their sweet cherry studies, where authors showed that different water DI strategies did not result in significant differences for number of flower buds. However, Blanco et al. [49] also indicated that different level of water deficit resulted in an inconsistent response on the number of flower buds.
In this study, DI treatments showed significant effect on flowering and fruit set including the measures of SDF, NOFB, FBA, and FS (Tables 4-6, Figures 5 and 6). Among bud and flower measures, cultivars showed significant effects on SDF, NOFB and FS but not on FBA (Tables 4-6, Figures 5 and 6). Similarly to our results, several previous studies on fruit trees showed that periods of flowering and fruit set were sensitive to water restrictions [50][51][52][53][54][55][56][57][58][59][60][61][62]. Several studies found that water stress during flowering and fruit set stages decreased fruit set and increased the abscission of flowers, fall of young fruits; and consequently, reduced yield on citrus trees [51,53,54], sweet oranges [60], pear [57], apricot [55], and peach [61]. On the other hand, some other studies showed that the DI strategy enhanced flowering on pomegranate [62], loquat [58], pear [50], mango [52], and apricot [56]. In addition, some further studies showed no effects of DI on flowering and fruit sets on sweet cherry [49] and apricot [15,18,23,55,63]. Some earlier studies emphasized that lower fruit sets of deciduous trees occurred only under severe water deficit conditions during the postharvest period [64]. It was also reported that the use of a moderate water deficit strategy caused no effect on fruit sets if it was used in the non-sensitive phenological periods [15,18,23,49]. The various results of the studies suggest that joined factors can play a role in the effect of DI on flowering and fruit set. Such an essential combined effect can be the level of water deficit and crop load [65], as the high fruit load can decrease the fruit growth with increasing water deficit and may also delay the recovery from the water stress.
In this study, ANOVA showed that SVG was significantly affected by the DI treatments and cultivars (Table 6) and SVG significantly decreased with severe water deficit ( Figure 6B). The increasing water stresses significantly lowered SVG on cv. 'Canino' in all DI treatments; however, on cv. 'Ninfa', only the most severe water deficit treatment lowered SVG considerably ( Figure 6B). Many previous studies showed that yearly vegetative growth significantly decreased in the DI treatments on apricot [11,18,66], pomegranate [43], lemon [67,68], peach and pear [6,61,[69][70][71][72], apple [7,73], sweet cherry [49,[74][75][76], and olive [77]. Consequently, the RDI technique can be a useful tool to manipulate vegetative growth [78,79]. However, for a successful application of RDI, a clear separation is needed between the vegetative and fruit growth periods, because a detrimental effect can occur on yield if vegetative and fruit growth processes overlap.
In this study, ANOVA showed that FNT and FY were significantly affected by DI treatment and cultivars (Table 6) while FW was significantly affected by only cultivars (Table 6). In case of cv. 'Ninfa' both yield measures significantly increased in moderate RDI (RDIm) treatments compared to the fully irrigated control one ( Figure 6C,D). In previous studies, the fruit growth was not reduced during the early fruit growth stages when moderate DI managements were applied compared with the fully irrigated trees. Moreover, even a stimulated fruit growth occurred when irrigation was increased to 100% ETc during the subsequent phenological stages. This was confirmed on apricot [23,64], citrus [80,81], pear [70,82] and apple [83] although no such results have been reached in other studies [61,67]. In this study, FNT and FY were significantly decreased in the severe water deficit treatment (CDI) on both apricot cultivars ( Figure 6C,D). Our results coincided with the study of Peres-Pastor et al. [66], who showed that that RDI reduced yield when water deficits were severe. Such severe deficit induced an additional decrease in the vegetative growth of the trees which resulted in a significant reduction in the fruit numbers per tree. Peres-Pastor et al. [66] also showed that moderate water deficit treatments provided similar yields to the well irrigated treatment when the water deficit was applied before the fruit growth stage III; which also coincided with the results of this study in the RDIm treatment for cv. 'Ninfa'. Therefore, the final stages of fruit growth (stage III) have to be considered as the most sensitive period for apricot irrigation. This is in agreement with other studies as different water deficit intensities during the fruit growth stage III caused a reduction of fruit diameter and fruit number, which adversely affected yield too [54,64,66,68,81,84,85].
The overall results on vegetative and generative growth indicated relative clear separations between shoot and fruit growth periods for apricot trees, which give the option to apply an RDI strategy successfully. Yield can remain unaffected with an appropriate use of RDI, while shoot growth can be reduced and water can be saved [18,23,86,87]. The strong connections between vegetative and yield parameters were also supported by the strong correlation coefficient values between SVG and FNT, and between SVG and FY in all DI treatments of this study (Table 9). In addition, this study showed overall strong associations among measures of bud vs. vegetative tree parts vs. generative tree parts vs. yield for all DI treatments; in correlation analyses for instance between flower bud vs. opened flower bud, fruit bud vs. shoot growth, fruit bud vs. yield, and fruit set vs. yield (Table 9); and in factor analyses of PC1, for instance, among bud numbers, fruit set, and yield ( Figure 10).
This study showed that DI treatments, cultivars, and years had significant effects on the carbohydrate contents of apricot leaves (Tables 7 and 8). In agreement with this study, Nasrabadi et al. [43] also found that the effect of DIs, cultivars, and years were significant for non-structural carbohydrate content in pomegranate leaves. In this study, the carbohydrate contents were gradually decreased by the reductions of supplied water. In line with our study, Soliman et al. [7] and Nasrabadi et al. [43] reported that the leaf area and non-structural carbohydrates significantly decreased on apple and pomegranate, respectively, in water deficit treatments compared to the well irrigated ones. Wong et al. [88] showed that a reduction in carbohydrate contents of leaves under water deficit conditions is due to the reduction in the photosynthesis of the leaves. This was in line with the study of Romero et al. [89] who reported that the photosynthetic capacity significantly decreased in young branches of almond trees under DI conditions. The reduction in the photosynthesis of leaves under water deficit conditions is related to the demolition of chlorophyll structure [90]. In addition, xylem structures also changed under water deficit conditions, which caused a decrease in the water movement from the roots to the leaves, and the leaves closed their stomata to prevent water losses [91]. The closed stomata reduce the photosynthetic rates of the leaves and the plant use the non-structural carbohydrate reserves [92]. Wong et al. [88,93] also showed that the amount of stored starch is essential for seasonal growth and tree vitality. In agreement with the above statements, the results of this study showed that SVG and TSC contents decreased in the DI treatments and they have significant inter-correlations with each other ( Figure 6B, Tables 7-9). Thus, a reduced water supply caused a reduction in the photosynthetic activity of the apricot leaves, which resulted in decreased carbohydrate contents and a subsequent delay in the vegetative growth of apricot tree.
Our study showed significant negative correlations between the content on carbohydrate and chilling requirements in all DI treatments (Table 9). In agreement with this study, Nasrabadi et al. [43] also showed that pomegranate trees with high CR values had a low concentration of non-structural carbohydrates at the end of the growing season. Nasrabadi et al. [43] also noted that the correlation between the carbohydrate concentration and the heat requirement was nonsignificant, which corroborated with our non-significant relationships between TSC on HR in both correlation and regression analyses for apricot (Table 9).
This study showed that DI treatments and cultivars had significant effects on leaf proline contents (Table 7 and Figure 7). Treatments with increasing water deficit resulted in a significant increase in the leaf proline content for both apricot cultivars compared to the fully irrigated control treatment (Figure 7), which were in agreement with the study of Faghih et al. [94] on apple and Laajimi et al. [16] on apricot. Previous studies showed that proline accumulation in plants plays a protective role against the environ-mental stresses [95][96][97], such as limited availability of water due to soil water depletion or high temperature stress inducing an increased transpiration rates [98,99]. Proline is commonly increased during drought stress, which eliminates ROS and protects the plant cell membranes and antioxidant enzymes [100], and thus, the increased proline content can effectively prevent leaves from water loss [101,102]. Our study also indicated that apricot cultivars can accumulate different levels of proline, which was also in agreement with previous results on citrus trees [98,99].

Conclusions
The shortage of water resources and low chilling hours are considered as key challenges for apricot production in saline soil under semiarid climate conditions such as Egypt. Therefore, here we showed (i) overall effects of three water saving DI treatments (RDIm, RDIs, and CDI) on 17 measures of tree physiological and chemical properties, (ii) prediction based model analyses for the 17 measures and (iii) the inter-correlations of the 17 measures specified to DI treatments on two apricot cultivars in Egypt.

•
Results showed that values of CR, HR, TNB, SDF, NOFB, FS, SVG, FNT, FY, TSC, and TPC values were affected significantly by both DI treatments and cultivars. Values of FBA were significantly affected by years and DI treatments, while sGDD by years and cultivars. Values of DEDFH, NFB, and FW were significantly influenced only by cultivars, while NVB only by DI treatments.

•
Prediction based GLMs analysis demonstrated that GLMs can be predictors for the measured parameters in the DI treatments for the two cultivars and for the three years.
The best goodness-of-fit of the predicted GLMs was reached for the measures of HR, NOFB, FS, SVG, FNT, TSC, and TPC.

•
In all the four DI treatments, 22 pair-variables correlated significantly in Pearson correlation analyses, which were further confirmed by linear regression analysis. • PCA explained 82% of the total variance and PC1, PC2, and PC3 explained 23, 21, and 15% of the total variance and correlated with the HR, TNB, FS, FNT, and FY; FBA, SVG, TSC, and TPC; and NFB, NVB, and NOFB, respectively, indicating strong relationships among tree physiological and chemical properties. • Treatment effect, prediction model based analyses and inter-correlation analyses together showed that the RDIm treatment and cv. 'Ninfa' provided the most acceptable values for most measured properties compared to the fully irrigated control treatment.
Overall results showed that DI techniques using moderate water deficits (RDIm) can be managed successfully under semiarid Mediterranean climate conditions if suitable DI timings and appropriate cultivars are selected for apricot growing.