Optimization of the Regulated Deficit Irrigation Strategy for Greenhouse Tomato Based on the Fuzzy Borda Model

It is of great significance to explore the strategy of regulated deficit irrigation (RDI) under mulched drip irrigation to stabilize tomato yield and improve quality and efficiency. This experimental study was conducted on a drip-irrigated greenhouse in two consecutive years (2020 and 2021). Three deficit levels were set for the flowering and fruit development stage (Stage I), and three were set for the fruit-ripening stage (Stage II). As a combination evaluation method, the fuzzy Borda model was used to optimize the RDI strategy of greenhouse tomato. The results showed that the net photosynthetic rate, stomatal conductance, transpiration rate, and total shoot biomass of tomato decreased with an increase in the water deficit, while the intercellular CO2 concentration had an opposite trend. The mild and moderate water deficit at Stage I reduced tomato yield by 16–24% and 30–40% compared to full irrigation. The water deficit at Stage II was able to improve various quality parameters and the water-use efficiency of tomato; the irrigation water-use efficiency (32.8–33.9 kg/m3) and leaf water-use efficiency (3.2–3.6 μmol/mmol) were the highest when the soil water content was 70–90% θf (field capacity) at Stage I and 40–60% θf at Stage II (T3). Based on the fuzzy Borda combination evaluation model, T3 was determined as the treatment with stable yield, high quality, and efficient irrigation under the experimental conditions. The irrigation regime was as follows: irrigating 20–25 mm in the transplanting stage, no irrigation in the seedling stage, irrigating 193.2–220.8 mm at Stage I, and then irrigating 27.6 mm at Stage II.


Introduction
Water is an important factor affecting crop yield and quality [1]. With the improvement of people's living standards, it is urgent for food production to shift from simply pursuing the improvement of yield to pursuing both yield and quality [2,3]. However, water shortages, high agricultural water consumption, and low utilization rates of irrigation water are the main factors restricting sustainable agricultural development [4,5]. In addition, crops have different water requirements at different growth stages [6,7]. On the basis of this background, regulated deficit irrigation (RDI) can produce appropriate water deficits in different growth stages of crops to improve crop quality and water-use efficiency without a significant yield reduction [8,9].
Tomato is one of the most common greenhouse crops [10]. The growth and physiological indexes of tomato were affected by RDI [11]. With the decrease in irrigation amount in each growth period, the dry matter mass of tomato decreased to different degrees [12,13]. In addition, crop photosynthesis also has a certain response to water deficits [14,15]. Studies have analyzed that RDI impacts photosynthesis mainly by affecting the stomatal closure of crops [16], and the degree of water deficit in different growth stages has different effects on the photosynthesis of crops [17]. At present, there are few reports on the responses of the 2 of 16 physiological and growth indexes of greenhouse tomato to the upper and lower limits of irrigation at different growth stages. Thus, these aspects must be further investigated.
RDI can significantly improve the fruit quality and water-use efficiency of tomato [18][19][20][21]. Lu et al. [22] analyzed 561 experimental groups and 145 control groups and found that RDI could increase the WUE of tomato by 2.33 kg/m 3 on average. Studies indicated that when irrigation was sufficient in other growth stages but halved in the seedling stage, there was no significant impact on tomato yield; however, water-use efficiency was significantly improved [23]. In addition, when the total irrigation amount was reduced by 40%, soluble sugar and VC increased by 20.26% and 16.10%, respectively [24]. The improvement of quality parameters such as the total soluble solids and VC of greenhouse tomatoes by deficit irrigation at the fruit-ripening stage was greater than that at other growth stages. Moreover, the yield, quality, and water-use efficiency of different soil types and tomato varieties responded differently to the degree of water deficits at different growth stages [22]. In combination with the existing research results and the soil entropy and water consumption of plants at the seedling stage, in this study, different degrees of water deficit tests were carried out at the flowering and fruit development stages (Stage I) and the fruit-ripening stage (Stage II) to investigate the effects of RDI on the growth indexes, physiology indexes, yield, quality, and water use of tomato.
When optimizing the RDI strategy of greenhouse tomato, the influence of water on a single index such as yield, quality, and utilization efficiency of irrigation water should be considered. It is more important to realize the multi-objective evaluation, which takes into account these three factors. At present, there are two commonly used multi-objective evaluation methods in the agricultural field: the single evaluation method and the combined evaluation model based on the single evaluation method. Common single evaluation methods include principal component analysis [25,26], grey relational degree analysis [27,28], membership function analysis [29,30], and TOPSIS [31,32]. In the same experiment, the results obtained by different single evaluation methods are often not completely consistent. For example, Li et al. [33] determined that the comprehensive ranking of the 50% E-pan irrigation level was better than that of the 70% E-pan level by principal component analysis and grey correlation analysis, while the results obtained by the TOPSIS method were opposite.
In order to further optimize the evaluation results, the researchers put forward a combination evaluation model based on a single set of evaluation methods. Under the guidance of the basic principles of comprehensive evaluation, a more effective combination was sought through the integration of methods to eliminate the random errors and system bias generated by using a single method. In this way, the inconsistency of multi-method evaluation results was solved [34]. Hu et al. [35] used four different combination evaluation models to comprehensively evaluate the yield and quality indexes of tomato under different water and nitrogen treatments. The results showed that the fuzzy Borda combination evaluation model was easy to calculate and offered the most obvious advantages. On the basis of this result, the fuzzy Borda combination evaluation model was used in the present study to comprehensively evaluate the yield, quality, and water use of greenhouse tomato under RDI in order to optimize the RDI strategy of greenhouse tomato.
Based on the different sensitivity levels of tomato to water deficits in different growth stages, this study set different water deficit degrees in the two water-sensitive stages (Stage I and Stage II) and carried out nine types of regulated deficit combination experiments on greenhouse tomato in two consecutive years. The objectives of this study were to (1) analyze the effects of RDI on the physiological and growth indexes, yield, quality, and water use of greenhouse tomato and (2) optimize the RDI strategy for greenhouse tomato through a fuzzy Borda combination evaluation model in order to improve the comprehensive quality and irrigation water-use efficiency as much as possible without a significant yield reduction.

Experimental Site Description
The experiment was conducted in the solar greenhouse of Liujiapu Tomato Industrial Park (37 • 64 N, 112 • 48 E), Taiyuan, Shanxi Province, China. The experimental area features a typical temperate continental monsoon arid climate, with an average annual temperature of 11 • C, rainfall of 520 mm and evaporation of 1812.7 mm, annual average sunshine duration of 2672 h, and a frost-free period of 202 days. The textural class of the experimental soil was sandy loam, and the basic physical and chemical properties in the 0-60 cm layer are shown in Table 1. The total area of the greenhouse used for the experiment was 600 m 2 (60 m × 10 m), and the area of each plot was 45 m 2 (7.5 m × 6 m). Drip irrigation under plastic mulch was used in the experiments. Each ridge consisted of two rows of tomatoes and two rows of drip irrigation strips placed inside the tomato plants, with row spacing of 60 cm and plant spacing of 50 cm. The planting density was 33,333 plants/ha. Five to seven days after transplanting, the surface was covered with black plastic film with a width of 1.5 m and a thickness of 0.008 mm. In order to prevent the lateral exchange of soil water between plots, an impervious film with a depth of 60 cm was buried between two adjacent plots.

Experimental Design
The tomato variety "Shouyan PT326" was used in the experiment, which was transplanted at the three-leaf stage (10 May 2020; 12 May 2021). To ensure the survival rate of tomato after transplanting, the planting water was irrigated with 20 mm and 25 mm in 2020 and 2021, respectively. The growth period of tomato was divided into the seedling stage (10 May to 6 June 2020; 12 May to 15 June 2021), flowering and fruit development stage (7 June to 11 August 2020; 16 June to 23 August 2021), and fruit-ripening stage (12 August to 18 September 2020; 24 August to 22 September 2021).
According to the water sensitivity of tomato in combination with the upper and lower limits of irrigation, the RDI experiment was conducted at Stage I and Stage II. A total of nine treatments was used, and each treatment was set with three replicates. The specific irrigation scheme is shown in Table 2, where 70-90% θ f indicates that the lower limit and upper limit of irrigation water were 70% and 90% of field capacity, respectively; the other levels were similar. From Stage I, the soil moisture content was measured every 7 days in vertical layers (0-60 cm), 15 cm away from the plant. When the average soil moisture content fell below the irrigation lower limit, the irrigation amount was calculated, and the soil was irrigated by a drip irrigation system. Before transplanting, 3/5 of nitrogen fertilizer (300 kg/hm 2 ) and potassium fertilizer (300 kg/hm 2 ), as well as all phosphate fertilizer (200 kg/hm 2 ) and organic fertilizer (20,000 kg/hm 2 ), were applied evenly into the tillage layer. At the fruit development stages of the first and third ear, 1/5 of N and K fertilizer was applied, respectively. The topdressing fertilizer was dissolved in buckets and applied to the soil with water through a drip irrigation system. The nitrogen, phosphorus, and potassium fertilizer contained urea (W N ≥ 46.4% represents the content of nitrogen is greater than or equal to 46.4%, similarly hereinafter), calcium magnesium phosphate (W P2O5 ≥ 15.0%), and potassium chloride (W K2O ≥ 57.0%), respectively. The physiological and growth indexes of tomato were measured at two deficit growth periods. The net photosynthetic rate (Pn), stomatal conductance (Gs), transpiration rate (Tr), and intercellular CO 2 concentration (Ci) of three evenly growing tomatoes were measured in each treatment. WUE L = Pn/Tr represents the water-use efficiency based on leaf level. The photosynthetic indexes were measured by a Li-6400 portable photosynthetic system. The measurement time was 09:00-12:00 in the morning, and the height and orientation of tomato leaves taken during measurements remained consistent. Next, the stem, leaf, and fruit of the tested plants were separated and immediately placed in an oven at 105 • C for 30 min and dried at 75 • C until they reached a constant weight. Then, the dry weight was measured using a balance with a precision of 0.01 g.

Yield and Irrigation Water Utilization Efficiency
For each tomato plant, we left 4 ears of fruit and then picked the core. An electronic balance with a precision of 0.05 kg was used to measure the yield of tomato in each treatment. The calculation formula for irrigation water utilization efficiency is as follows: where WUE Y represents the irrigation water utilization efficiency (kg/m 3 ), Y represents the yield (kg/hm 2 ), and I represents the total irrigation amount (m 3 /hm 2 ).

Fruit Quality Parameters
When the third ear of tomato was ripe, 9 fruits were randomly selected from the top, middle, and bottom parts of each treatment plant to determine the quality. The Fruit shape index (FSI) was taken as the ratio of the longitudinal diameter to the average transverse diameter of the tomato. Fruit firmness (Fn) was measured by using a fruit hardness tester (GY-4,Aipu, Quzhou, China) 3 times on the fruit body and shoulder and then taking the average value. The soluble sugar (SS) was determined by anthrone colorimetry [36] and organic acids (OA) titrated with 0.1 mol/L NaOH [3]. The sugar/acid content ratio (SAR) Agriculture 2022, 12, 324 5 of 16 was calculated by the ratio of soluble sugar to organic acid. Vitamin C (VC) was determined by the 2, 6-dichlorophenol sodium indophenol titration [37].

Model Application and Methods
The yield, quality, irrigation water-use efficiency, and water-use efficiency based on the leaf level of the tomato were taken as evaluation indexes. The principal component analysis model (PCA) [38], grey relational degree analysis model (GRA) [39], TOPSIS model based on the analytic hierarchy process (TOPSIS-AHP) [40,41], and the membership function analysis model (MFA) [42] were used to evaluate the yield, quality, and water use of tomato. A Kendall Concorde coefficient test [43] was used to test the consistency of single evaluation methods to determine whether the four single evaluation methods were compatible. The fuzzy Borda combined evaluation model [44] was used to comprehensively evaluate the results of the various single evaluation methods. The specific calculation steps of this model are as follows: (1) Calculate the membership degree of the single evaluation method µ ij : where X ij represents the score of i treatment under the j method, and µ ij represents the membership degree of i treatment under the j evaluation method.
(2) Calculate the fuzzy frequency: Calculate the fuzzy Borda number: where Q hi-j represents the score of i treatment in position h under the j evaluation method: After being sorted according to F i , a larger value indicates a higher overall score.

Statistical Analysis
IBM SPSS Statistics 26 was used for the correlation analysis, regression analysis, PCA, and analysis of variance, while the LSD method was used for significance analysis at a p ≤ 0.05 level. Excel was used to calculate the path coefficient, GRA, MFA, TOPSIS-AHP, prior inspection, and fuzzy Borda combination evaluation model.

Photosynthetic Characteristics of Tomato
The results indicated that Pn, Gs, and Tr presented similar patterns under RDI ( Figure 1) The results indicated that Pn, Gs, and Tr presented similar patterns under RDI (Figure 1). Pn, Gs, and Tr were 34.4-77.7%, 3.0-55.7%, and 40.6-64.1% more plentiful at Stage I than Stage II. Compared with full irrigation, mild and moderate deficit irrigation at Stage I reduced Pn, Gs, and Tr by 3.2-11.1% and 12.1-25.8%, 12.7-25.0% and 25.2-41.6%, and 16.9-25.5% and 27.9-44.1%, respectively. When the irrigation level was the same at Stage I, compared to mild deficit irrigation, moderate and severe deficit irrigation at Stage II reduced Pn, Gs, and Tr by 2.1-27.8% and 5.4-47.6%, 0.4-23.4% and 2.5-32.1%, and 0.6-19.1% and 2.1-25.4%, respectively. The effect of RDI on Ci was opposite to that on Pn, Gs, and Tr. When full irrigation and mild deficit irrigation were applied at Stage I, the water deficit at Stage II had no significant effect on other indexes except for Gs and Tr treated by T1 and Pn treated by T6. While moderate deficit irrigation was applied at Stage I, the water deficit at Stage II had more significant effects on Pn and Tr.  Figure 2 shows the effects of RDI on the aboveground biomass accumulation and allocation of tomato plants in the two consecutive years. At Stage I, the total biomass of the shoot under full irrigation was 7.5-14.0% and 19.0-35.5% higher, respectively, than that under mild and moderate water deficits. The average proportion of fruit weight at Stage I was 19.8-27.5% and 27.6-35.7% higher than that of the stem and leaf, respectively. The proportion of fruit weight in T4 and T7 was 3.0-4.1% and 1.6-2.6% higher than that in T1, respectively. Under the same irrigation level at Stage I, compared with the mild water deficit, the moderate and severe water deficit reduced the total biomass of the shoot  Figure 2 shows the effects of RDI on the aboveground biomass accumulation and allocation of tomato plants in the two consecutive years. At Stage I, the total biomass of the shoot under full irrigation was 7.5-14.0% and 19.0-35.5% higher, respectively, than that under mild and moderate water deficits. The average proportion of fruit weight at Stage I was 19.8-27.5% and 27.6-35.7% higher than that of the stem and leaf, respectively. The proportion of fruit weight in T4 and T7 was 3.0-4.1% and 1.6-2.6% higher than that in T1, respectively. Under the same irrigation level at Stage I, compared with the mild water deficit, the moderate and severe water deficit reduced the total biomass of the shoot by 5.3-14.7% and 8.6-22.5%, respectively, at Stage II. Except for T7 and T1, which were the same in 2021, T2-T9 accounted for 0.9% to 6.8% more fruit weight than T1 at Stage II.

Yield and Water-Use Efficiency
The results of two consecutive years showed that the water deficit at Stage I had a significant effect on tomato yield, while the water deficit at Stage II had no significant effect on tomato yield (Table 3). Compared with full irrigation, the yield of tomato decreased by 16-24% and 30-40%, respectively, due to mild and moderate water deficits at Stage I. Reducing the irrigation amount at Stage II helped improve the WUEY. Compared with the mild water deficit, the moderate and severe water deficits increased the WUEY by 2-10% and 11-21%, respectively. In addition, when the total irrigation amount was the same, the irrigation amount was 20-26% greater at Stage I than at Stage II, and the WUEY of the tomato was increased by 9-20%. A higher yield and lower irrigation amount led to the highest WUEY under T3, which was 24% higher in 2020 and 23% higher in 2021 than the lowest value (T4 and T7), respectively. As can be seen from Table 3, WUEL in T3 was the highest, which was, respectively, 48% and 37% higher than the lowest value (T9) in the two consecutive years.

Yield and Water-Use Efficiency
The results of two consecutive years showed that the water deficit at Stage I had a significant effect on tomato yield, while the water deficit at Stage II had no significant effect on tomato yield (Table 3). Compared with full irrigation, the yield of tomato decreased by 16-24% and 30-40%, respectively, due to mild and moderate water deficits at Stage I. Reducing the irrigation amount at Stage II helped improve the WUE Y . Compared with the mild water deficit, the moderate and severe water deficits increased the WUE Y by 2-10% and 11-21%, respectively. In addition, when the total irrigation amount was the same, the irrigation amount was 20-26% greater at Stage I than at Stage II, and the WUE Y of the tomato was increased by 9-20%. A higher yield and lower irrigation amount led to the highest WUE Y under T3, which was 24% higher in 2020 and 23% higher in 2021 than the lowest value (T4 and T7), respectively. As can be seen from Table 3, WUE L in T3 was the highest, which was, respectively, 48% and 37% higher than the lowest value (T9) in the two consecutive years.  Table 4 presents the data provided from the experiments on fruit appearance quality (FSI), nutritional quality (VC, OA), flavor quality (SS, SAR), and storage and transportation quality (Fn) of greenhouse tomato in the two consecutive years. It can be seen from the table that the improvement effect of the water deficit on the tomato quality indexes was more significant at Stage II than at Stage I. At the same irrigation level at Stage I, compared with the mild water deficit at Stage II, moderate and severe water deficits increased the FSI by 2.2-4.7% and 5.4-13.3%; the VC value increased by 7.9-35.0% and 22.8-62.0%; the OA value increased by 7.0-30.2% and 12.3-43.9%; the SS increased by 9.8-30.9% and 18.9-58.0%; and the SAR increased by 0.3-25.5% and 19.2-34.9% due to the increasing effect of the deficit water on soluble sugar was greater than that on organic acid. In addition, the fruit firmness of moderate and severe water deficit was 18.4-45.0% and 32.3-63.6% higher than that of mild water deficit, which was beneficial to fruit storage and transportation. In conclusion, severe water deficit at Stage II results in the best performance of tomato fruit quality parameters in the two consecutive years.

Evaluation Results of a Single Evaluation Method
Four single evaluation methods (PCA, GRA, MFA, and TOPSIS-AHP) were used to evaluate the tomato yield, quality parameters, WUE Y , and WUE L ( Table 5). It can be seen from the table that the standard deviation for the ranking of the four evaluation methods in 2020 and 2021 ranged from 0 to 1.63 and 0 to 1.26. In the two consecutive years, 6-7 of 9 treatments had a standard deviation of less than 1, while those of the other treatments were all greater than or equal to 1. Therefore, there are still some differences in the evaluation results of the four single evaluation methods. Moreover, the fuzzy Borda combined evaluation model can be used to carry out a combined evaluation of the four single evaluation methods, which can eliminate the differences between single evaluation methods.

Evaluation Results of Fuzzy Borda Combination Evaluation
Before using the fuzzy Borda combination evaluation, the single evaluation methods were tested in advance. The Kendall correlation coefficients of the order values for the four single evaluation models [45] are shown in Table 6. Here, the mean correlation coefficients between the evaluation values of each single method and other methods are between 0.67 and 0.87, and the comprehensive correlation between GRA and the other three methods is presented to be the strongest in the two consecutive years. A Kendall-W Concorde coefficient test was further used for the significance test [35,43]. The Kendall-W Concorde coefficients were calculated to be 0.91 and 0.95 (W), respectively, in the two consecutive years. And 29.00 and 30.27 for χ 2 =M(N − 1)W, respectively, in the formula, M is the number of single evaluation methods is 4, and N is the number of treatments is 9. All χ 2 values were greater than χ 2 0.01/8 =17.54, indicating that the four methods were compatible and satisfied the prior consistency test. The fuzzy Borda combination evaluation model was used to combine the results of the above single evaluation methods ( Table 7). As can be seen from the table, there was strong consistency in the evaluation results between the two consecutive years. The comprehensive evaluation value reached its maximum under T3 (35.40 and 33.78, respectively), while the comprehensive evaluation score was the lowest under T4 and T7. The results indicated that the optimal irrigation scheme would involve a combination of soil moisture content of 70-90% θ f and 40-60% θ f at Stage I and Stage II, respectively.

Effects of RDI on Photosynthesis and Its Relationship with Yield
Photosynthesis is an important physiological process of crops. The RDI affects the photosynthesis of plants through stomatal restriction and nonstomatal restriction [46]. In this paper, Ci values increased with an increase in the water deficit, indicating that the decrease in stomatal conductance was the result of a decrease in photosynthetic rate [47]. In addition, after the moderate water deficit at Stage I, the mild water deficit at Stage II was shown to have a certain compensation effect on Pn, Gs, and Tr, which is similar to the results of Yu et al. [48] and may be the result of stomatal resistance, mesophile conductivity, biochemistry, or other factors [49].
Correlation analysis was conducted between the photosynthetic indexes and yield in two growth periods (Figure 3). Except for Ci at Stage II, the correlation between photosynthetic indexes and yield in both growth stages reached a significant level (p ≤ 0.05). In addition, the correlation coefficients of Ci, Tr, and tomato yield at Stage I were higher than those at Stage II, possibly because Stage I is the key growth period for tomato yield, during which the tomato root produces abscisic acid and water losses through xylem transport to the plant leaves, thereby decreasing stomatal opening and inhibiting crop leaf transpiration [50,51]. The low correlation coefficients between Pn, Gs, and yield at Stage I may be caused by the large interannual difference of Pn and Gs.

Effects of RDI on the Yield, Quality, and Water Use of Tomato
The effect of RDI on crop yield is related to the intensity, duration, and growth stage of deficit irrigation [1,11,22]. Adequate or mild deficit irrigation at Stage I is critical for yield formation [52]. This study also concluded that the effect of the water deficit at Stage I on tomato yield was greater than that at Stage II. Liu et al. [53] obtained the following results from potted experiments in the tomato flowering and fruiting stages with different irrigation levels: the WUEY was the highest when the irrigation level was 80-90%θf at Stage I. In this study, the irrigation water-use efficiency was highest when the irrigation level was 70-90%θf at Stage I and 40-60%θf at Stage II (T3). When the total irrigation amount was the same, a greater irrigation amount at Stage I was beneficial to improve the utilization efficiency of the irrigation water.
In this study, RDI was carried out during two water-sensitive periods for tomato plants. The effect of a water deficit on the quality parameters at Stage II was greater than that at Stage I. With an increase in the water deficit at Stage II, the VC, OA, SS, and SAR improved. These results agree with those of Chen et al. [3] and may be because an increase in the water deficit reduced the fluid flow to the phloem, as well as fruit dilution, leading to an increase in the concentration of various fruit components [15,54]. In addition, with an increase in the water deficit, tomato fruit firmness gradually increased, possibly due to a decrease in the internal expansion of the fruit under conditions with insufficient water, leading to a decrease in the pressure of the cell wall, thereby improving the elasticity of the fruit epidermis [2].

Effects of RDI on the Yield, Quality, and Water Use of Tomato
The effect of RDI on crop yield is related to the intensity, duration, and growth stage of deficit irrigation [1,11,22]. Adequate or mild deficit irrigation at Stage I is critical for yield formation [52]. This study also concluded that the effect of the water deficit at Stage I on tomato yield was greater than that at Stage II. Liu et al. [53] obtained the following results from potted experiments in the tomato flowering and fruiting stages with different irrigation levels: the WUE Y was the highest when the irrigation level was 80-90% θ f at Stage I. In this study, the irrigation water-use efficiency was highest when the irrigation level was 70-90% θ f at Stage I and 40-60% θ f at Stage II (T3). When the total irrigation amount was the same, a greater irrigation amount at Stage I was beneficial to improve the utilization efficiency of the irrigation water.
In this study, RDI was carried out during two water-sensitive periods for tomato plants. The effect of a water deficit on the quality parameters at Stage II was greater than that at Stage I. With an increase in the water deficit at Stage II, the VC, OA, SS, and SAR improved. These results agree with those of Chen et al. [3] and may be because an increase in the water deficit reduced the fluid flow to the phloem, as well as fruit dilution, leading to an increase in the concentration of various fruit components [15,54]. In addition, with an increase in the water deficit, tomato fruit firmness gradually increased, possibly due to a decrease in the internal expansion of the fruit under conditions with insufficient water, leading to a decrease in the pressure of the cell wall, thereby improving the elasticity of the fruit epidermis [2].

Relationship between Tomato Yield, Growth, and Physiology Indexes
Studies have analyzed that there is a positive correlation between the yield, biomass, and photosynthetic characteristics of tomato leaves [55][56][57][58]. With a decrease in the irrigation amount, all the photosynthetic indexes, root biomass, and yield of tomato decreased [59]. In this study, different deficit levels were set at Stage I and Stage II. The yield, biomass, and leaf photosynthetic indexes (Pn, Gs, Tr) of tomato in both years were all decreased by water deficits, which was consistent with previous research results.
Path analysis [60,61] was also applied, taking the irrigation amount, photosynthetic index (X Pn , X Gs , X Ci , X Tr ), and dry weight of shoot organs (X stem , X leaf , and X fruit ) at Stage II of greenhouse tomato as independent variables and yield (Y) as dependent variables. The optimal regression equation was obtained by the stepwise regression method: Y = 57.853 + 1.022X Pn − 0.112X Ci + 0.008X fruit . The maximum direct path coefficient of dry fruit weight on the yield was 0.571, the minimum direct path coefficient of Ci on the yield was −0.241, and the maximum indirect path coefficient of Pn through dry fruit weight on the yield was 0.512. Lu et al. [62] reached similar conclusions through a correlation analysis, showing that postanthesis dry matter accumulation and final aboveground biomass play key roles in the final yield of wheat.

Comprehensive Evaluation of Yield, Quality, and Water Use of Tomato
In this study, four single evaluation methods were used to evaluate the yield, quality, and water-use indexes of tomato. The standard deviation was greater than 0 for ranking 6-7 of the 9 treatments. This phenomenon emerged because the mechanism for each single evaluation method is different, with different angles and emphasis placed on the use of information, combined with the artificial factors that exist in the process of evaluation [63]. In addition, the grey correlation analysis method has the highest correlation with the other three evaluation methods, as well as low data requirements, a lower workload, and simpler calculations [64]. It is, therefore, recommended to use the grey correlation analysis method when selecting a single evaluation method.
In this study, the fuzzy Borda combination evaluation model was used to evaluate the yield, quality, and water use of tomato under RDI. The results of the experiments were largely consistent between the two years of the study. Further analysis of the correlation between the standardized yield, average value of each quality parameter, and irrigation quota ( Figure 4). We concluded that under an irrigation quota of 266 mm, the tomato yield and comprehensive quality were balanced, the standardized yield and quality parameters were above 0.8, and the irrigation water-use efficiency was higher. These results were consistent with those of the fuzzy Borda model, which showed that this method can be used to optimize crop RDI strategies. In addition, the fuzzy Borda combination evaluation model overcame the disadvantage of extracting only the order-value information of a single evaluation method without the evaluation value information; moreover, the calculation process is relatively simple [35]. It is, therefore, suggested to use the fuzzy Borda model for the comprehensive evaluation of a single evaluation method. consistent with those of the fuzzy Borda model, which showed that this method can be used to optimize crop RDI strategies. In addition, the fuzzy Borda combination evaluation model overcame the disadvantage of extracting only the order-value information of a single evaluation method without the evaluation value information; moreover, the calculation process is relatively simple [35]. It is, therefore, suggested to use the fuzzy Borda model for the comprehensive evaluation of a single evaluation method.

Conclusions
The leaf Pn, Gs, Tr, aboveground biomass, and yield of greenhouse tomato decreased with an increase in the water deficit degree, while the Ci value increased with a decrease in the water deficit degree. The application of an appropriate water deficit (T2 and T3) facilitated the distribution of nutrients to fruit. Compared to full irrigation, the yield of tomato decreased by 16-24% and 30-40% under mild and moderate water deficits at Stage I. Controlling water content between 40% θ f to 60% θ f at Stage II improved the fruit quality parameters. When the total irrigation amount was the same, the irrigation amount was 20-26% greater at Stage I than at Stage II, and the WUE Y of tomato increased by 9-20%. In the two consecutive years (2020 and 2021), the WUE L and WUE Y of T3 were the highest at 3.2-3.6 µmol/mmol and 32.8-33.9 kg/m 3 , respectively. Using the fuzzy Borda combination evaluation model, the yield, quality, and water use of tomato in the two consecutive years were comprehensively evaluated. The results indicated that the optimal irrigation regime was as follows: irrigation of 20-25 mm in the transplanting stage, no irrigation in the seedling stage, irrigation of 193.2-220.8 mm at Stage I, and irrigation of 27.6 mm at Stage II. The optimal RDI strategy matches with the T3. The RDI may affect the growth physiology, yield, quality, and water-use indexes of crops by first affecting hormones and enzymes. However, this result needs further analysis.