Linear Relationship of a Soil Total Water Potential Function and Relative Yield—A Technique to Control Salinity and Water Stress on Golf Courses and Other Irrigated Fields

: A simple empirical approach is proposed for the determination of crop relative yield (%) through the soil total water potential (kPa). Recurring to decimal logarithms, from analytical exponential expressions, a linear simple relationship of soil total water potential Ψ t (matric Ψ m + potential Ψ o) function and crop relative yield was studied and developed. The combination of the salinity model, the soil water retention model and the matric potential approach were used to reach this objective. The representation of turfgrass crop relative yield (%) versus a function of soil total water potential f( Ψ t) values was shown through a log-normal graph (y = a + mx); the log scale axis “y” (ordinates) deﬁnes relative yield Yr, being two the origin ordinate “a” and “m” the slope; the normal decimal scale axis “x” (abscissa) is the function of soil total water potential f( Ψ t). Hence, it is possible, using only two experimental points, to deﬁne a simple linear relation between a function of soil total water potential and crop relative yield, for a soil matric potential value lower than − 20 kPa. This approach was ﬁrst tested on golf courses (perennial turfgrass ﬁelds), but it was further decided to extend it to other annual crop ﬁelds, focused on the model generalization. The experimental plots were established, respectively, in Algarve, Alentejo and Oeiras (Portugal) and in the North Negev (Israel). Sprinkler and trickle irrigation systems, under randomized blocks and/or water and salt gradient techniques, were used for water application with a precise irrigation water and salt distribution. Results indicated that there is a high agreement between the experimental and the prediction values (R 2 = 0.92). Moreover, the precision of this very simple and easy tool applied to turfgrass ﬁelds and other irrigated soils, including their crop yields, under several different sites and climatic conditions, can contribute to its generalization.


Introduction
Lack of water and high salinity of soil and irrigation water are high-cause great problems for crop productions in arid and semi-arid regions [1][2][3][4]. The major effect of soil salinity is a reduction in plant water uptake [5]. Hence, the salinity condition in the root zone hinders moisture extraction from soil by plants because of osmotic potential development in soil water, due to the presence of salts, which ultimately decreases transpiration of plants, and thereby affects crop yield [6,7]. Therefore, the additional presence of salinity stress decreases the relative yield response due to water stress; similarly, under water stress, the relative yield response to increasing salinity is reduced [8]. On the other hand, when saline water is the only source available to the farm, it is necessary to avoid the reduction of irrigation doses to prevent excessive salt accumulation in the root zone with unavoidable effects on crop yield [9]. Hence, it is assumed that water uptake depends on matric and osmotic potentials [10], being the soil total water potential a sum of several soil potentials, especially the matric and osmotic potential [11]. Soil water retention curve plays an important role in simulating soil water movement and assessing soil water holding capacity and availability [12]. Soil water retention curves are crucial for characterizing soil moisture dynamics and are particularly relevant in the context of irrigation management [13]. The quantification of the soil water potential is necessary for a variety of applications in both agricultural and horticultural systems, such as optimization of irrigation volumes and fertilization; it has been employed in a variety of studies and applications to optimize irrigation schedule and to investigate ways in which water-saving irrigation can help to save water, increase water productivity and optimize production [14]. However, the direct measurement of soil parameters for unsaturated soils requires complicated laboratory tests and it is often expensive and time consuming, according to [15][16][17]. Linear relationships of transpiration and yield have been computed for various crops and climates under different soil water conditions [18] and under various salinity conditions [19]. Investigations into water potentials in the soil-plant system are of great relevance in environments with abiotic stresses, such as salinity and drought [20]. Later, soil salinity and matric potential interactions on yield response were studied [19]. The simplified diffusion convection equation to obtain production functions, including the effects of water, salinity and nutrition conditions, was solved by [21]. Later, a model was presented in which the wilting point is a function of the soil salt content [22] at a higher salinity, the water content at wilting point is higher than at low salinity, resulting in an insufficient amount of available water, and, therefore, a reduced yield. This model shows that the movement of salts in the soil is solely dependent on the movement of water in soil; it shows that the effect of salinity is simulated by its effect on the wilting point, thus reducing the soil available water content [23]. However, when the crop foliage is wet by sprinkling with saline water, plants are subject to additional salt damage [24]. These salt effects were studied on corn leaves [1] and lettuce [25].
An inverted logistic exponential equation for quantifying crop salt tolerance was presented by [26], showing how salt tolerance data can be used by coupling an appropriate salt balance model with a least-squares optimization method. Later, it was conducted during a thorough study on this model that was "crowned "and widely used with the initials vGG [27]. These models describe the plant as a pipeline of water, and, therefore the water uptake and transpiration are synonymous terms such that the yield, which is dependent upon the transpiration rate, is given as a unique function of soil water potential or soil osmotic potential [28]. To account for the dynamic processes and the simulation of water and salts in the soil, plant and atmosphere continuum, several wellknown models were used successfully, such as the SWAP (Sil Water Atmosphere Plant) model [29]. This model was applied to a high number of simulations, such as: (1) to account for various salinity effects in field crops irrigated with saline water, simulating soil profiles of salinity and water content comparing them with observed data; moreover, it is compared to measured and calculated transpiration from a field experiment with several salinity treatments [30], (2) to simulate saline water irrigation for seed maize [31] and (3) to model soil water-heat dynamic changes in seed-maize fields [32]. Later, other models were developed, such as HYDRUS [33], SALTMED [34] and LAWSTAC [16]. However, they did not present the linearity and the simplicity of this.
The objectives of this study were: 1.
To develop a simple new model describing the effect of total soil water potential Ψt (matric Ψm + osmotic Ψo) (matric + osmotic) on relative crop production.

2.
To replace the standard models, generally demanding a high number of laboratorial determinations due to a higher number of experimental points, and to be represented by curves not as attractive as the straight line defined only by two experimental points.

Theory
According to [2], the theoretical soil-water retention curve is given by: where: A and B are soil parameters specific for each curve and θ w is soil gravimetric content (kg water/100 kg soil). The graphical representation of θ versus (pF) 2 values in a log-normal graph, with (pF) 2 in the normal scale and θ in the log scale, shows, in fact, a straight-line pattern [2] from pF 2 (field capacity zone) to pF 4.2 (wilting point zone). It was already recently suggested by [35] that the log-logistic and lognormal distributions are more suitable to model soil's pore distribution than other tested distributions. As the dry bulk density is constant for each specific soil type, and relation between gravimetric soil water content θ w and volumetric soil water content θ is depending solely on the dry bulk density, it shall be used as the soil volumetric content θ (% = m 3 water/100 m 3 soil) for our further studies. However, the Soil Science Society of America [36] considers obsolete the pF concept, and, therefore, matrix potential will be expressed by "h" cm (water), through the pF concept, as follows: pF = ln |Ψm| (2) Combining Equations (1) and (2), and using the soil volumetric water content θ instead of the gravimetric soil water content θ w , it will be obtained which can be written as: Equation (3) was studied and developed by [3] and it was confirmed by its high rigor for |Ψm| ≥ 100 cm H 2 O (retention zone). This is due to cavitation initiated by entrapped air bubbles or the liquid's own vapor pressure [37]. This air entry point is around this value and provokes a second branch curve explained by a curve inflexion that occurs near that point. These two right-lined segments with different slopes individualized [38], one of them in the retention zone and the other in the drainage zone (|Ψm| ≤ 100 cm H 2 O). Therefore, for |Ψm| ≥ 100 cm H 2 O, the exponential curve was logarithmic, and it was obtained as a high correlated linear function for a large range of soils [3], varying the correlation coefficient from 0.956 (light soils) up to 0.999 (heavy soils), as follows: By combining Equations (3) and (4), using the soil volumetric water content θv [% (m 3 water m −3 soil)] instead of the gravimetric soil water content θw and by using Ψm expressed in kPa, it will be obtained by the expression: Recurring to logarithms, Equation (4) takes the linear form as follows, and it obtained a high correlated linear function for a large range of soils [3], as follows: For non-saline soils, it was assumed by [39] that the soil available water content readily available to plants θ ASW is the difference between the volume of soil water content at field capacity (θ fc ) and at permanent wilting point (θ wp ), and Equation (8) being Ψm fc = −33 kPa value (retention zone limits the value for mineral soils and attributed to field capacity), according to the concepts of: [38][39][40][41][42][43]. Equation (9) takes the form Under soil water stress conditions, stomatal conductance decreases and, consequently, transpiration T (mm) and CO 2 are reduced [44]. Moreover, observations have shown a linear relation between irrigated crops and climates, under conditions of water deficit [45,46], as follows: where k is a homogeneity factor, obtained by [47]: where m is a crop factor and E 0 is average seasonal free water evaporation. The crop yield response to soil total water potential will be presented on a relative basis, which offers a simple way and a uniform manner of presenting data from different crops, locations and years.
According to Equations (9)- (12), it may be written ln YrM − ln Yr = Ch {[log (10 |Ψm2|)] 2 − 6.25} (13) where Ψm2 < −33 kPa; YrM represents the maximal relative yield of the irrigated crop (100%) and Yr the relative yield (%), being Fh a specific factor of homogeneity. YrM represents 100% of the relative yield and Ψm2 is the actual soil matric potential value. Equation (13) may be expressed by 2 − ln Yr = Fh {[log (10 |Ψm2|)] 2 − 6.25} (14) Equation (14) may be reduced to The salt concentration of a relatively diluted soil solution is roughly linearly related to the electrical conductivity (EC), which in itself is linearly related to the osmotic potential (Ψo). Because EC (dS m −1 ) is easily measured, it is advantageous to express it in terms of EC [48] and to convert it to osmotic potential Ψo (cm H 2 O), The soil total water potential Ψt can be expressed as the algebraic sum of the soil component potentials, which effects are acting on soil water behavior, as follows: where: Ψm, Ψo, Ψp and Ψg are, respectively, the matrix potential, soil osmotic potential, pressure potential and the gravitational potential. All these soil component potentials are additive and among them soil matric potential Ψm is the dominant potential, followed by the osmotic potential Ψo [49], and therefore Equation (17) may be reduced to Ψt = Ψm + Ψo (18) or, combining Equations (16) and (18): where Ψm may be obtained, respectively, from Equations (4)- (6) or, more easily, through its determination by the graphical representation than by the algebraic procedure. A crop production system is characterized by the link between the crop yield Y and the factor involved in it.
The crop yield response to soil total water potential will be presented on a relative basis, which offers a simple way and uniform manner of presenting data from different crops, locations and years. This response is given as a function the soil total water potential Ψt, as follows: Linear relationships between water use and yield have been modeled for various crops and climates under conditions of water deficit [50,51] and conditions of salt stress [11,[52][53][54]. On the other hand, Hanks and co-authors [55,56] have described the equation of soil water flow within a plant root extraction, as follows where t is the time (s), z the depth (m), K the soil hydraulic conductivity (m s −1 ), Ψh is the soil hydraulic potential expressed by the sum Ψm + Ψg (Pa) and A(z) is the plant root extraction function, which depends on Ψm and Ψo (Pa). The crop yield response to total soil water potential will be presented on a relative basis, which offers a simple way and uniform manner of presenting data from different crops, locations and years. According to the influence of the salt concentration on the availability of soil water [48] Equation (15) takes the following form: where Ψt YrM represents the total soil water potential when relative yield Yr reaches 100%. Equation (22) shows a linear function for a large range of irrigated crops, and may be easily graphical represented by a straight line, in a log-normal scale, as follows: where: the log scale axis "y" (ordinates) defines relative yield Yr (%), being 2 the original ordinate "b" and "m" the slope; the normal decimal scale axis "x" (abscissa) is the f (Ψt), according to Equation (22), and it is represented by

Materials and Methods
The experimental turfgrass fields were established in several golf courses and other lawn fields in Faro, Algarve, Portugal [57].  Table 1 shows the studied sprinkle irrigated turfgrasses and their distribution on the different areas of the turfgrass fields [58]. Sprinkler and trickle irrigation systems under randomized blocks or point/line source techniques were used for water application with a precise irrigation water distribution, as follows: (1) Randomized blocks (sprinkle irrigation system) experimental design was applied to corn forage [57], according to Figure 1.  Sprinkler and trickle irrigation systems under randomized blocks source techniques were used for water application with a precise irrigation bution, as follows: (1). Randomized blocks (sprinkle irrigation system) experimental design w corn forage [57], according to Figure 1. (2). Sprinkler double line source experimental design was applied to cor cabbage [59], according to Figure 2. (2) Sprinkler double line source experimental design was applied to corn, grain and cabbage [59], according to Figure 2. (3). Sprinkler point source experimental design was applied to turfgrass and sunflower [60], according to Figure 3. (4). Double emitter source experimental design was applied to cabbage and lettuce [61], according to Figure 4. (3) Sprinkler point source experimental design was applied to turfgrass and sunflower [60], according to Figure 3. (3). Sprinkler point source experimental design was applied to turfgrass and sunflower [60], according to Figure 3. (4). Double emitter source experimental design was applied to cabbage and lettuce [61], according to Figure 4. Christiansen [62] coefficient of water distribution uniformity (CUC) was always above 80% for sprinklers [9] and above 90% for drip irrigation [63]. The plots were irrigated once a day. To control soil water along the soil profile, soil water content was monitored periodically during the experiment and gravimetrically measured for a 0.0-0.3.m depth. It was based upon the direct determination of the moisture content and dry weight of the material in the oven at 105 ° C until constant weight. Soil water retention was determined by measuring the water content at six different matric potential values, which were determined by water content retained through the Richards pressure-membrane extraction apparatus. Extraction of soil solution was conducted recurring to suction cups. Seedbed and basic fertilization were made according to regional conventional agro-techniques Tables 2 and 3 show, respectively, the soil properties of turfgrass fields and soil properties and irrigation methods of the other irrigated crop fields. Replications number of measured soil volumetric water content θv for each matric potential Ψm value was always higher than 4. Christiansen [62] coefficient of water distribution uniformity (CUC) was always above 80% for sprinklers [9] and above 90% for drip irrigation [63]. The plots were irrigated once a day. To control soil water along the soil profile, soil water content was monitored periodically during the experiment and gravimetrically measured for a 0.0-0.3 m depth. It was based upon the direct determination of the moisture content and dry weight of the material in the oven at 105 • C until constant weight. Soil water retention was determined by measuring the water content at six different matric potential values, which were determined by water content retained through the Richards pressure-membrane extraction apparatus. Extraction of soil solution was conducted recurring to suction cups. Seedbed and basic fertilization were made according to regional conventional agro-techniques. Tables 2 and 3 show, respectively, the soil properties of turfgrass fields and soil properties and irrigation methods of the other irrigated crop fields. Replications number of measured soil volumetric water content θv for each matric potential Ψm value was always higher than 4. The effects of treatments were evaluated using the analysis of variance (ANOVA), when no cause-and-effect relationship was known. The statistical test Dunnett T3 was selected, in order to identify the statistical difference among multiple mean values at the 95% significance level. Direction from [64][65][66] was used. When problems of lack of randomization were known due to the point source experimental design, a geostatistical approach was undertaken [67]; this approach shows that simple random sampling and the calculation of an average, usually used for the normal procedure of soil sampling in Agriculture, is not always the best answer. These geostatistical methods describe the spatial variability and help to produce standard deviation maps, showing the confidence of the samples taken in an area, where trend removal and direction of anisotropy of some soil properties was facilitated with kriging. Moreover, spatial variable experiments substitute, very efficiently and economically, conventional experimental designs, like randomized blocks due to use of much lower areas, less pollution of the environment, elimination of borders and research costs being saved, such as equipment, water, energy, fertilizers, plants, crops, pesticides, workers and management.

Results and Discussion
Tables 4 and 5 show the function of the soil total water potential f(Ψt) and the data needed to its calculation, according to the Equations (22) and (24), respectively, to the experimental turfgrass fields and to the other irrigated crops fields. In Table 4, it is shown that the effect of soil water content variations on the function of the soil total water potential f(Ψt), for several turfgrass fields and sites under different salinity levels.   Table 5 shows the effect of both variations-soil water and salinity different levels-on the function of the soil total water potential f(Ψt), for several other irrigated crop fields and sites. Yield response results to soil total potential (Ψt) show the general tendency of the relation between yield and soil water and salt concentration corresponds to well-known results published in other scientific papers: the lower the soil total potential (Ψt), the lower relative yield [65][66][67]. Table 6 shows the turfgrass yield response to wastewater, ground water and potable water application, for five trials (wastewater VLW-ECw = 2.1 dS m −1 , groundwater VLG-1.2 m −1 , wastewater SGW-ECw = 2.4 dS m −1 , UAW-ECw = 1.6 dS m −1 and wastewater SG-ECw = 2.4 dS m −1 ). For a low soil matric potential (Ψm close to −1500 kPa), water was the limiting factor, increasing the function of the soil total water potential f(Ψt) being relative yield Yr sharply reduced. Hence, it may be seen that grass yield was lower [higher f(Ψt)] on the wastewater treatments, especially the VLW near the dry level, and higher yield [lower f(Ψt)] was obtained on the UAP treatment due to the osmotic pressure values of the soil. For intermediate soil matric potential values (−33 kPa > Ψm > −240 kPa), wastewater application triggered slightly higher yields compared to potable water application, especially if Ψm was close to higher values (−10 kPa > Ψm > −33 kPa). This was due to the probably lower concentration of nutrients of the potable water when compared to the wastewater. On the other hand, Yr was enhanced with an increasing rate of soil water content and, therefore, with the decrease of f(Ψt), especially if Ψm is close to higher values (−10 kPa > Ψm > −33 kPa). This was due probably to the lower concentration of nutrients of the potable water when compared to the wastewater.  Table 7 shows several other irrigated crop yields (corn grain CgS, corn forage CfS, sunflower SfS, lettuce LD, cabbage CS and cabbage CD) response to water and salt application, underground water and potable (sprinkler S or drip D irrigation). The negative slope means that the relative yield (logarithmic Y axis) decreases with the enhance of the function of the soil total water potential f(Ψt), shown in decimal abscissa axis. Several major aspects may be seen: for the corn grain that, for very high content saline water (6.2 dS m −1 ) and a f(Ψt) = 3.32, there is still a relative yield at about 76%; this is due to the large amounts of irrigation water used and, therefore, the water used to increase the leaching of salts, being the yields slightly influenced by the salinity effects. For low water application, the effect of water in forage yield was more pronounced that the corn forage yields were more sensitive than the grain corn yield to the [(f(Ψt)] due the to the pronounced effects of water amounts in forage being yield lower than 70% when [(f(Ψt)] > 4. Sunflower was also very sensitive to the lack of water, decreasing sharply the yields (Yr < 15%) when [(f(Ψt)] > 8. The lettuce yields were highly influenced by the salinity effects under high salinity levels. Hence, from 2.5 dS m −1 to 3.9 dS m −1 , the yield decreased up to 61% [(f(Ψt) = 0.95], and for more than 8.5 dS m −1 yield reached only less 40% [(f(Ψt) > 2.87]. This was due probably to the use of drip instead of sprinkler irrigation (the leaves were not wet by the saline water), and soil water content was around the field capacity. In relation to cabbage, on sprinkle irrigation system (low soil osmotic pressure) cabbage relative yield was decreasing monotonically with the increase of the f(Ψ0 %t); on drip irrigation system, with the increase of soil osmotic potential from dS m −1 0.9 up to 3, dS m −1 the relative yield Yr decreased near 40% [(f(Ψt) > 2.4].  Table 8 shows the regression analysis of the relationship between turfgrass relative yield Yr (%) and the function of the soil total water potential f(Ψt). It can be seen that the coefficient of determination (0.95 < R 2 < 1.00) is very high for field conditions.
The sensitivity of the yield was higher on the wastewater treatments, namely VLW (sharper increase) due to the higher salinity level when compared to potable water treatment UAP, explained by soil leaching. On the other hand, the enhancement was lower for the UAP treatment (monotonic increase), due to its lower concentration of nutrients combined with soil leaching. It can be seen that the slope is very close to 1, the intercept is quite small, and the coefficient of determination R 2 (0.95 < R 2 < 1.00) is very high. Table 8. Regression analysis of the relationship between turfgrass relative yield Yr (%) and a function of the soil total water potential f(Ψt).  Table 9 shows the regression analysis of the relationship between the relative yield Yr (%) of several irrigated crops and a function of the soil total water potential f(Ψt).

Turfgrass
It can be seen that the coefficient of determination R 2 (0.89 < R 2 < 0.96) is high for field conditions. Table 9. Regression analysis of the relationship between the relative yield Yr (%) of several irrigated crops and the function of the soil total water potential f(Ψt).

Crop
Regression The results show that the soil moisture and soil salt concentration characteristics curves could be approximated by exponentials, for a soil matric potential (Ψm) higher than −10 kPa. The hypothesis is advanced according to which the yield Y = f(Ψt) curves (Equation (20)) may be represented by a straight line (Equation (22)) and, therefore, with a linear graphical representation in a suitable axis system.
Relationship between total simulated and observed relative yields is presented in Figure 5. The slope is very close to 1, the intercept is quite small and the coefficient of determination R 2 (0.92) is very high for field conditions; however, its value is relatively lower than the R 2 obtained by the logarithmic values of the relative yield (Tables 8 and 9), once natural values are used instead of logarithmical values.
It shows that the regression is highly significant, and, therefore, the predicting ability of this approach is very good and capable of describing the relative yield response to the function soil total water potential f(Ψt).
Given the importance of this approach, more data that are experimental should be obtained to increase the number of model simulations. Hence, in the future, it is advisable to do additional research in order to obtain more modelling results, being higher than the validation of this approach.
once natural values are used instead of logarithmical values.
It shows that the regression is highly significant, and, therefore, the predicting ability of this approach is very good and capable of describing the relative yield response to the function soil total water potential f(Ψt).
Given the importance of this approach, more data that are experimental should be obtained to increase the number of model simulations. Hence, in the future, it is advisable to do additional research in order to obtain more modelling results, being higher than the validation of this approach. Recurring to decimal logarithms, from analytical exponential expressions, a linear simple relationship of soil total water potential Ψt (matric Ψm + potential Ψo) function and crop relative yield was studied and developed, according to Equations (22)- (24). The process is displayed by a flowchart that is given below (Figure 6). Recurring to decimal logarithms, from analytical exponential expressions, a linear simple relationship of soil total water potential Ψt (matric Ψm + potential Ψo) function and crop relative yield was studied and developed, according to Equations (22)- (24). The process is displayed by a flowchart that is given below (Figure 6).  (1) and (2) by processing data collected from experiments on golf courses.

Conclusions
This work showed that a simple empirical tool described the crop response to the soil total water potential Ψt (matric Ψm + osmotic Ψo) due to the low number of the involved parameters and a rapid and accurate determination. Moreover, it is possible, with only two experimental points, to define the above relationship for soil matric potential Ψm values lower than −20 kPa. On the other hand, the establishment of conventional soil water

Conclusions
This work showed that a simple empirical tool described the crop response to the soil total water potential Ψt (matric Ψm + osmotic Ψo) due to the low number of the involved parameters and a rapid and accurate determination. Moreover, it is possible, with only two experimental points, to define the above relationship for soil matric potential Ψm values lower than −20 kPa. On the other hand, the establishment of conventional soil water retention functions is a slow process, demanding a high number of laboratory determinations due to a higher number of experimental points; also, they are generally presented by curves not as attractive as the straight lines obtained by this simple approach. The results showed a high agreement between the experimental and the predicted values (R 2 = 0.92). Moreover, the precision of this tool applied to grass fields and other different irrigated fields, under different soil types and climatic conditions, can contribute to its generalization. There are multiple numbers of applications of this empirical tool, mainly related to salinity and water stress, involving the planning and management of irrigation (water quality, amounts and frequency) and desalination projects, soil leaching, fertilization enhance and the use of clean desalination techniques (decrease of irrigation amounts and use of drought tolerant and salt removing species).