Modeling the Effects of Irrigation Water Salinity on Growth, Yield and Water Productivity of Barley in Three Contrasted Environments

Freshwater scarcity and other abiotic factors, such as climate and soil salinity in the Near East and North Africa (NENA) region, are affecting crop production. Therefore, farmers are looking for salt-tolerant crops that can successfully be grown in these harsh environments using poor-quality groundwater. Barley is the main staple food crop for most of the countries of this region, including Tunisia. In this study, the AquaCrop model with a salinity module was used to evaluate the performance of two barley varieties contrasted for their resistance to salinity in three contrasted agro-climatic areas in Tunisia. These zones represent sub-humid, semi-arid, and arid climates. The model was calibrated and evaluated using field data collected from two cropping seasons (2012–14), then the calibrated model was used to develop different scenarios under irrigation with saline water from 5, 10 to 15 dS m−1. The scenario results indicate that biomass and yield were reduced by 40% and 27% in the semi-arid region (KAI) by increasing the irrigation water salinity from 5 to 15 dS m−1, respectively. For the salt-sensitive variety, the reductions in biomass and grain yield were about 70%, respectively, although overall biomass and yield in the arid region (MED) were lower than in the KAI area, mainly with increasing salinity levels. Under the same environmental conditions, biomass and yield reductions for the salt-tolerant barley variety were only 16% and 8%. For the salt-sensitive variety, the biomass and grain yield reductions in the MED area were about 12% and 43%, respectively, with a similar increase in the salinity levels. Similar trends were visible in water productivities. Interestingly, biomass, grain yield, and water productivity values for both barley varieties were comparable in the sub-humid region (BEJ) that does not suffer from salt stress. However, the results confirm the interest of cultivating a variety tolerant to salinity in environments subjected to salt stress. Therefore, farmers can grow both varieties in the rainfed of BEJ; however, in KAI and MED areas where irrigation is necessary for crop growth, the salt-tolerant barley variety should be preferred. Indeed, the water cost will be reduced by 49% through growing a tolerant variety irrigated with saline water of 15 dS m−1.


Introduction
The world food supply is affected by environmental abiotic stresses, which damages up to 70% of food crop yields [1][2][3]. In the Near East and North Africa (NENA) region, physical water scarcity is already affecting food production [4]. The NENA region is characterized by an arid climate with a total annual rainfall much lower than the evapotranspiration of the field crops. In the Arab World, more than 85% of the available water resources are used for agriculture [5]. Despite this high-water allocation for the agriculture sector, about 50% of food requirements are imported [4]. Crop irrigation uses poor quality groundwater, which is saline in nature. The uninterrupted application of groundwater for irrigation is replete, which leads to a severe increase in soil salinity and reduction in crop yields. Climate changes, namely the increase in global temperatures and the decline in rainfall, exacerbate soil salinization, resulting in loss of production in arable lands [6]. According to recent estimates, one-fifth of the irrigated lands in the world are affected by salinity. Every day, on average, 2000 ha of irrigated land in arid and semi-arid areas is adversely affected by salinity problems [7]. The annual economic loss due to these increases in soil salinity is about USD 27.3 billion [8].
Cereals are the main crops in the Mediterranean and NENA regions, contributing to food security and social stability. Barley is one of these staple crops in the area. However, its production is constrained by abiotic factors, such as the arid climate, low and erratic rainfall, and soil and water salinity. The anticipated climate changes will further increase the negative impacts of these factors in the future [9]. Barley (Hordeum vulgare L.) is a drought-and salt-tolerant crop with considerable economic importance in Mediterranean and NENA regions since it is a source of stable farm income [10]. Indeed, barley is a staple food for over 106 countries in the world [11]. Barley is characterized by its high adaptability from humid to arid and even Saharan environments. Barley is grown in many areas of the world and is used for feed, food, and malt production [2,12].
To improve barley production in these regions, plant scientists have adopted a strategy to identify tolerant genotypes for maintaining reasonable yield on salt-affected soils [13]. Crops physiologists and breeders are working to assess how efficient a genotype is in converting water into biomass or yield. To do so, they use production parameters, with which measurement in field experiments is difficult and time-consuming. However, these complex parameters can be determined with the help of crop growth simulation models [6,13]. Dynamic simulation models describe the growth and development of crops based on the interaction with soil, water, and climate parameters. Models can be used to simulate soil and water salinity and crop management practices on the growth and yield of crops under different agro-climatic conditions [6].
Models were used to test the impact of salinity on crops under different environmental conditions and different fertilization practices [14,15].
AquaCrop is a water-driven dynamic model (Vanuytrecht et al., 2014). AquaCrop is a simulation model to study crops' water productivity. As crop-water-productivity is affected by climatic conditions, it is crucial to understand water productivity's response to changing rainfall and temperature patterns [9].
Among the available models, AquaCrop is preferred due to its robustness, precision, and the limited number of variables to be introduced [16]. It uses a small number of explicit and intuitive parameters that require simple calculation [16]. AquaCrop is a software system developed by the Land and Water Division of FAO to estimate water use efficiency and improve agricultural systems' irrigation management practices [17,18].
Water productivity (WP) can be described as the ratio of crops' net benefits, including both rain and irrigation.
According to [19], irrigation management organizations are interested in the yield per unit of irrigation water applied, as they have to improve the yield through human-induced irrigation processes. However, the downside is that not all irrigation water is used to generate crop production. Therefore, FAO defines water productivity as a ratio between a unit of output and a unit of input. Here, water productivity is used exclusively to indicate the amount or value of the product over the volume or value of water that is depleted or diverted [20]. This model was developed by the Food and Agriculture Organization (FAO) [16,21]. AquaCrop simulates the response of crop yield to water and is particularly suited to regions where water is the main limiting factor for agricultural production. The model is based on the concepts of crops' yield response to water developed by Doorenbos and Kassam [22]. The AquaCrop model (v4.0) published in 2012 can estimate yield under salt stress conditions.
The AquaCrop model has been used to predict crop yields under salt stress conditions in different parts of the world [23,24]. Kumar et al. [23] successfully used the AquaCrop model to predict the water productivity of winter wheat under different salinity irrigation water regimes. Mondal et al. [24] used AquaCrop to evaluate the potential impacts of water, soil salinity, and climatic parameters on rice yield in the coastal region of Bangladesh. The AquaCrop model has also been widely used to simulate yields of various crops under diverse environments. For example, barley (Hordeum vulgare L.) [5,25,26], teff (Eragrostis teff L.) [5], cotton (Gossypium hirsutum L.) [27], maize (Zea mays L.) [28] wheat (Triticum aestivum L.) [3].
In this study, the AquaCrop model (v4.0) is used to assess the performance of two barley genotypes under three contrasted agro-ecosystems (soil, salinity, and climate). In these areas, groundwater is primarily used for irrigation. The salinity of irrigation water ranges from 3 to 15 dS m −1 . Farmers do not know which barley variety is most tolerant to producing a reasonable yield under these saline environments. Furthermore, model simulations were also performed to evaluate the impact of three irrigation water salinity levels (5, 10, and 15 dS m −1 ) on the barley yield. A cost-benefit analysis was performed to determine the economic returns of each level of salinity water irrigation and genotype tolerance based on model simulation results. Those results should help recommend the farmers of saline areas to enhance barley yield and economic return.  Figure 1). Two barley varieties (Konouz from Tunisia and Batini 100/1 B from Oman) were used for field experiments. The Konouz variety is salt-sensitive [29,30], whereas Batini 100/1 B is salt-tolerant [29,31].

Description of Field Trial Sites
Agronomy 2020, 10, x FOR PEER REVIEW  3 of 19 This model was developed by the Food and Agriculture Organization (FAO) [16,21]. AquaCrop simulates the response of crop yield to water and is particularly suited to regions where water is the main limiting factor for agricultural production. The model is based on the concepts of crops' yield response to water developed by Doorenbos and Kassam [22]. The AquaCrop model (v4.0) published in 2012 can estimate yield under salt stress conditions.
The AquaCrop model has been used to predict crop yields under salt stress conditions in different parts of the world [23,24]. Kumar et al. [23] successfully used the AquaCrop model to predict the water productivity of winter wheat under different salinity irrigation water regimes. Mondal et al. [24] used AquaCrop to evaluate the potential impacts of water, soil salinity, and climatic parameters on rice yield in the coastal region of Bangladesh. The AquaCrop model has also been widely used to simulate yields of various crops under diverse environments. For example, barley (Hordeum vulgare L.) [5,25,26], teff (Eragrostis teff L.) [5], cotton (Gossypium hirsutum L.) [27], maize (Zea mays L.) [28] wheat (Triticum aestivum L.) [3].
In this study, the AquaCrop model (v4.0) is used to assess the performance of two barley genotypes under three contrasted agro-ecosystems (soil, salinity, and climate). In these areas, groundwater is primarily used for irrigation. The salinity of irrigation water ranges from 3 to 15 dS m −1 . Farmers do not know which barley variety is most tolerant to producing a reasonable yield under these saline environments. Furthermore, model simulations were also performed to evaluate the impact of three irrigation water salinity levels (5, 10, and 15 dS m −1 ) on the barley yield. A cost-benefit analysis was performed to determine the economic returns of each level of salinity water irrigation and genotype tolerance based on model simulation results. Those results should help recommend the farmers of saline areas to enhance barley yield and economic return.  Figure 1). Two barley varieties (Konouz from Tunisia and Batini 100/1 B from Oman) were used for field experiments. The Konouz variety is salt-sensitive [29,30], whereas Batini 100/1 B is salt-tolerant [29,31].  In Kairouan (KAI) and Medenine (MED) field trial sites were divided into two sub-plots. Each subplot was irrigated by one water salinity treatment (EC = 2 and 13 dS m −1 ). Three blocks were defined perpendicularly to the sub-plots so that both treatments were observed in each block. As Beja is located in the rainfed cereal growing area of Tunisia, no irrigation was applied.

Description of Field Trial Sites
The weather data characterize the trials sites related to temperature, and rain was described by [29]. The irrigation water applied and reference evapotranspiration (ET o ) registered in the trials during the two growing seasons are presented in Table 1. The collected data from each site were used to estimate the reference evapotranspiration (ET) according to the Penman-Monteith Evapotranspiration FAO-56 Method, and then the total water supplied was determined for each site to obtain the water barley requirement. Irrigation was applied using a drip system. To ensure water supply homogeneity, line source emitters were installed at each planting row and 33-cm spacing between emitters on the same row. Soil samples were taken from the trial sites, and physico-chemical analyses were performed. The site's soil characteristics are diverse, from soil rich in clay and organic matter in BEJ to sandy soil with impoverished organic matter continent in MED (Table 2). Crops were sown during the last week of November. Seeds were hand sown at the rate of 200 viable grains per m 2 . Nitrogen, potassium and phosphorus were applied separately at 85, 50, and 50 kg/ha rates, respectively.
At the five different stages, plants for each genotype, from three small areas (25 × 25 cm) were taken from each experimental unit and used to determine the biomass. At a final harvest stage, plot (1 × 2 m) was used for biomass and grain yield assessment. Water productivity (WP) was calculated as the ration between the collected yield expressed in kg ha −1 and the daily transpiration simulated by the model.

Description of the AquaCrop Model
The model describes soil, water, crop, and atmosphere interactions through four sub-model components: (i) the soil with its water balance; (ii) the crop (development, growth, and yield); (iii) the atmosphere (temperature, evapotranspiration, and rainfall), and carbon dioxide (CO 2 ) concentration; and (iv) the management, such as irrigation and crop fertilization soil fertility.
The AquaCrop model is based on the relationship between the relative yield and the relative evapotranspiration [22] as follows where Y x is the maximum yield, Y a is the actual yield, ET x is the maximum evapotranspiration, ET a is the actual evapotranspiration, and K y is the yield response factor between the decrease in the relative yield and the relative reduction in evapotranspiration. The AquaCrop model does not take into account the non-productive use of water for separating evapotranspiration (ET) into crop transpiration (T) and soil evaporation (E) where ET = actual evapotranspiration, E = soil evaporation and Tr = the sweating of crop. At a daily time step, the model successively simulates the following processes: (i) groundwater balance; (ii) development of green canopy (CC); (iii) crop transpiration; (iv) biomass (B); and (v) conversion of biomass (B) to crop yield (Y). Therefore, through the daily potential evapotranspiration (ET o ) and productivity of water (WP*), the daily transpiration (Tr) is converted into vegetal biomass as follows where WP* is the normalized water productivity [32,33] relative to Tr. After the normalization of water productivity for different climatic conditions, its value can be converted into a fixed parameter [34].
The estimation and prediction of performance are based on the final biomass (B) and harvest index (HI). This allows a clear distinction between impact of stress on B and HI, in response to the environmental conditions Y = HI * (B) (4) where: Y = final yield; B = biomass; HI = harvest index. During the calibration and testing of the model, we calculated water productivity (WP) as presented by Araya et al. [5] where Y is the yield expressed in kg ha −1 and Tr is the daily transpiration simulated by the model.

Crop Response to Soil Salinity Stress
The electrical conductivity of saturation soil-past extracts from the root zone (ECe) is commonly used as an indicator of the soil salinity stress to determine the total reduction in biomass production, determines the value for soil salinity stress coefficient (K s, salt ).
The coefficient of soil salinity stress (Ks salt ) varied between 0 (full effect of stress of soil salinity) and 1 (no effect). The following equation determined the reduction in biomass B rel represents the expected biomass production under given salinity stress relative to the biomass produced in the absence of salt stress. The coefficient is adjusted daily to the average ECe in the root zone [35].
Then, the thresholds values are given for the sensitive and tolerant barley genotype and expressed in dS m −1 . This allows the estimation of the lower limit (EC en ) to which the soil salinity stress begins to affect the production of biomass and the upper threshold (EC ex ), in which soil salinity stress has reached its maximum effect.

Soil Salinity Calculation
AquaCrop adopts the calculation procedure presented in BUDGET [36] to simulate the movement and retention of salt in the soil profile. The salts enter the soil profile as solutes after irrigation with saline water or through capillary rise from a shallow groundwater table (vertical downward and upward salt movement). The average ECe in the compartments of the effective rooting depth determines the effects of soil salinity on biomass production.
To explain the movement and retention of soil water and salts in the soil profile, AquaCrop divides the soil profile into 2 to 11 soil compartments called "cells", depending on the type of soil in each horizon (clay, sandy horizon) and its saturated hydraulic conductivity (Ksat in mm/day). The salt diffusion between two adjacent cells (cell j and cell j+1) is determined by the differences in salt concentration and expressed by the electrical conductivity (EC) of soil water.
AquaCrop determines the vertical salt movement in response to soil evaporation, considering the amount of water extracted from the soil profile by evaporation and the wetness of the upper soil layer. The relative soil water content of the topsoil layer determines the fraction of the dissolved salts that moves with the evaporating water.
AquaCrop determines the vertical salt movement because of the capillary rise. Finally, the salt content of a cell is determined by Salt cell is the salt content expressed in grams salts per m 2 soil surface, Wcell its volume expressed in liter per m 2 (1 mm = 1 L/m 2 ), and 0.64 a global conversion factor used in AquaCrop to convert dS/m to g/L. The electrical conductivity of the soil water (ECsw) and of the electrical conductivity of saturation soil-past extract (ECe) at a particular soil depth (soil compartment) is calculated as where n is the number of cells in each soil compartment; θ is the soil water content (m 3 /m 3 ); θ sat is the soil water content (m 3 /m 3 ) at saturation; ∆z (m) is the thickness of the soil compartment and Vol% gravel is the volume percentage of the gravel in the soil horizon of each compartment.

Input and Output Variables of the Model
The model was calibrated using data from the growing season of 2012-2013 and evaluated using data from 2013-2014. Determining parameters for crop development and production, as well as water and salinity stress, was fundamental for calibrating the AquaCrop model. The parameters of climate, soil, and crop management used for the model calibration are presented in Table 3.

Statistical Parameters Used for the Calibration and Evaluation of Model
Several statistical indices were used to evaluate the performance of the model on the field measured data. These include Percentage Error (PE), Root Mean Square Error (RMSE), Model Efficiency (ME) and Coefficient of Determination (R 2 ).
Percentage Error (PE) was determined using the following equation where S i and O i are simulated and observed values, respectively. The root means square error (RMSE) [37] is presented by the following equation with the values of RMSE close to zero indicate the best model fit.
The model efficiency (ME) [38] was applied to assess the effectiveness of the model. The ME indicator compares the variability of prediction errors by the model to those of collected data from the field. If the prediction errors are greater than the data error, then the indicator becomes negative. The upper ME bound is at 1.
The coefficient of determination (R 2 ), as a result of regression analysis, is the proportion of the variance in the dependent variable (predict value) that is predictable from the independent variable (observed value) and is computed according to [35] R 2 is between 0 and 1.

Parameters Used for Model Calibration
In total, 26 input parameters were used for the model calibration ( Table 4). Out of these, 14 parameters were considered as "conservative" because they do not change with salinity and are independent of limiting or non-limiting conditions. These parameters include normalized crop water productivity and crop transpiration coefficient. The remaining 12 are site-specific (climate, water, and soil salinity) and crop-specific (tolerant or sensitive). These input parameters were adjusted during the calibration process to obtain better adequacy between the measured and simulation values. Table 4. Final values of different model input parameters obtained after calibration for two genotypes under different salinity levels (S1 = 2 dS m −1 ; S2 = 13 dS m −1 ).

Development of Different Scenarios
After calibration and evaluation, the model was used to assess the performance of two barley varieties under three water salinity conditions scenarios i.e., 5, 10, and 15 dS m −1 , using the weather data for the growing season 2013-2014.

The Economic Gain from the Use of a Unit of Water Consumed in the Tow Barley Varieties under Different Climatic Conditions
The economic productivity of two barley varieties was estimated using the average unit cost of one water cubic meter in Tunisia and the water use predicted by AquaCrop. The crop water economic productivity of the tolerant and the sensitive barley varieties as the measure of the biophysical and then economic gain from the use of a unit of water consumed were estimated by AquaCrop model in grain yield production [20]. This is expressed in productive crop units of kg/m 3 and money unit/m 3 .

Biophysical Environments Variability of Experimental Sites
The experiments are conducted in adaptability trials set up in three contrasting biophysical environments (from the sub-humid to the arid interior). These sites, namely Beja, Kairouan and Medenine, were selected on a North-South transect ( Figure 1). The soils of the trial sites are very diverse, from soil rich in clay and poor in organic matter in BEJ to sandy soil with poor organic matter continent in MED (Table 2). Beja's sub-humid site received annual rainfall of 472 and 413 mm respectively during the two cropping seasons. However, in the semi-arid and arid sites, low rainfall was registered. The arid site of MED received an annual rainfall of 81 mm during the first cropping season and 156 mm during the second. At Kairouan, the rainfall for the 2012/2013 and 2013/2014 seasons was 152 and 180 mm, respectively ( Table 1). As Beja is located in the rainfed cereal-growing area of Tunisia, no irrigation was applied. KAI and MED field trial sites, two different salinities (EC = 2 and 13 dS m −1 ) of water were used for irrigation.
Soil calcium and potassium content was higher in KAI as compared to MED. Soil sodium content changes during the different experimentation period following irrigation with saline water in KAI and MED, where sodium is the dominant element present in the saline irrigation water. The variation between sites might be explained by the variation in the cationic exchange capacity of the sandy soil and torrential character of the rainfall in this area (Table 2).

Biomass, Grain Yield, and Water Use Efficiency
The correlation between grain yield, biomass, and water productivity values for two barley genotypes showed that the observed and simulated values are closely co-related, as evidenced by the high R 2 values, i.e., 0.91, 0.93, and 0.89 for grain yield, biomass, and water productivity, respectively ( Figure 2).

Development of Different Scenarios
After calibration and evaluation, the model was used to assess the performance of two barley varieties under three water salinity conditions scenarios i.e., 5, 10, and 15 dS m −1 , using the weather data for the growing season 2013-2014.

The Economic Gain from the Use of a Unit of Water Consumed in the Tow Barley Varieties under Different Climatic Conditions
The economic productivity of two barley varieties was estimated using the average unit cost of one water cubic meter in Tunisia and the water use predicted by AquaCrop. The crop water economic productivity of the tolerant and the sensitive barley varieties as the measure of the biophysical and then economic gain from the use of a unit of water consumed were estimated by AquaCrop model in grain yield production [20]. This is expressed in productive crop units of kg/m 3 and money unit/m 3 .

Biophysical Environments Variability of Experimental Sites
The experiments are conducted in adaptability trials set up in three contrasting biophysical environments (from the sub-humid to the arid interior). These sites, namely Beja, Kairouan and Medenine, were selected on a North-South transect ( Figure 1). The soils of the trial sites are very diverse, from soil rich in clay and poor in organic matter in BEJ to sandy soil with poor organic matter continent in MED (Table 2). Beja's sub-humid site received annual rainfall of 472 and 413 mm respectively during the two cropping seasons. However, in the semi-arid and arid sites, low rainfall was registered. The arid site of MED received an annual rainfall of 81 mm during the first cropping season and 156 mm during the second. At Kairouan, the rainfall for the 2012/2013 and 2013/2014 seasons was 152 and 180 mm, respectively ( Table 1). As Beja is located in the rainfed cereal-growing area of Tunisia, no irrigation was applied. KAI and MED field trial sites, two different salinities (EC = 2 and 13 dS m −1 ) of water were used for irrigation.
Soil calcium and potassium content was higher in KAI as compared to MED. Soil sodium content changes during the different experimentation period following irrigation with saline water in KAI and MED, where sodium is the dominant element present in the saline irrigation water. The variation between sites might be explained by the variation in the cationic exchange capacity of the sandy soil and torrential character of the rainfall in this area (Table 2).

Biomass, Grain Yield, and Water Use Efficiency
The correlation between grain yield, biomass, and water productivity values for two barley genotypes showed that the observed and simulated values are closely co-related, as evidenced by the high R 2 values, i.e., 0.91, 0.93, and 0.89 for grain yield, biomass, and water productivity, respectively ( Figure 2 The correlation between observed and simulated values of biomass yield for two barley genotypes at three locations showed proximity (Figure 3), which indicates the excellent ability of the AquaCrop model to predict biomass yield under different agro-climatic conditions. The results also show that the sensitive barley variety at MED produces the lowest biomass for both irrigation water qualities. Similar trends were observed for grain yield, where the tolerant barley variety performed better than the sensitive variety regardless of the location and the quality of irrigation water.

Canopy Cover (CC)
The maximum and minimum CC were 85% and 30% in the sub-humid and arid areas, respectively. The salinity induces a 10% reduction in the CC in the sub-humid environment and 5-30% in the dry climate of MED. CC reduction under saline irrigation water is less noticeable in the tolerant variety than the sensitive variety for both salinity levels. However, in the rainfed area of Beja, the growth of both varieties was comparable. Figure 4 shows a strong correlation between measured and simulated CC values for both varieties of barley (R 2 = 0.91 and R 2 = 0.93). In general, a good match between the observed and the simulated CC was observed in all three locations. However, the model somewhat over-estimated CC in the rainfed environment of Beja and slightly under-estimated it in the other two situations.

Canopy Cover (CC)
The maximum and minimum CC were 85% and 30% in the sub-humid and arid areas, respectively. The salinity induces a 10% reduction in the CC in the sub-humid environment and 5-30% in the dry climate of MED. CC reduction under saline irrigation water is less noticeable in the tolerant variety than the sensitive variety for both salinity levels. However, in the rainfed area of Beja, the growth of both varieties was comparable. Figure 4 shows a strong correlation between measured and simulated CC values for both varieties of barley (R 2 = 0.91 and R 2 = 0.93). In general, a good match between the observed and the simulated CC was observed in all three locations. However, the model somewhat over-estimated CC in the rainfed environment of Beja and slightly under-estimated it in the other two situations.

Effects of Soil Salinity
The maximum soil salinity was in the arid and semi-arid areas irrigated with saline water, respectively. The soil salinisation dynamic depends on the salinity of irrigation water. However, in the rainfed area of Beja, we noted the absence of any salty issue. Figure 5 shows that the simulated soil salinity trend in the root zone (up to a depth of 0.7 m) corresponds very well with the measured values under different saline water regimes across different environments throughout the growing season. The observed and modeled soil salinity correlated well, with an R 2 of 0.96. Figure 5 shows that the model reliably simulated average root zone salinity when the crop is irrigated with low-salinity water (2 dS m −1 ). However, it slightly underestimated soil salinity under higher saline water conditions (13 dS m −1 ), particularly for the late growing season.

Effects of Soil Salinity
The maximum soil salinity was in the arid and semi-arid areas irrigated with saline water, respectively. The soil salinisation dynamic depends on the salinity of irrigation water. However, in the rainfed area of Beja, we noted the absence of any salty issue. Figure 5 shows that the simulated soil salinity trend in the root zone (up to a depth of 0.7 m) corresponds very well with the measured values under different saline water regimes across different environments throughout the growing season. The observed and modeled soil salinity correlated well, with an R 2 of 0.96. Figure 5 shows that the model reliably simulated average root zone salinity when the crop is irrigated with low-salinity water (2 dS m −1 ). However, it slightly underestimated soil salinity under higher saline water conditions (13 dS m −1 ), particularly for the late growing season.

Statistical Indices for AquaCrop Model Evaluation
The statistical indices derived for evaluating the AquaCrop model's performance in predicting soil water content, yield, canopy cover percent, biomass, and water productivity (WP) of barley genotypes under different saline water regimes across different environments are given in Table 5. All statistical parameters depict a strong correlation between simulated and observed values for model calibration and evaluation periods. The correlation between all statistical parameters remained almost the same for the calibration and evaluation period, which indicates the robustness of the

Statistical Indices for AquaCrop Model Evaluation
The statistical indices derived for evaluating the AquaCrop model's performance in predicting soil water content, yield, canopy cover percent, biomass, and water productivity (WP) of barley genotypes under different saline water regimes across different environments are given in Table 5. All statistical parameters depict a strong correlation between simulated and observed values for model calibration and evaluation periods. The correlation between all statistical parameters remained almost the same for the calibration and evaluation period, which indicates the robustness of the model prediction. Based on the model calibration and evaluation results, the model was found robust enough to calculate different scenarios.

Development of Different Scenarios
Due to a shortage of surface water, farmers of KAI and MED regions have no option than to use groundwater for irrigation. The quality of groundwater ranges from 4 to15 dS m −1 in these two regions. Farmers are interested to know which barley varieties would be most suitable to grow under these groundwater quality conditions. The calibrated and evaluated model was used to assess the performance of two barley varieties under three water salinity conditions i.e., 5, 10, and 15 dS m −1 , and the results are presented in Table 6. The performance of both barley varieties in the KAI area is predicted to be much higher than MED area under all salinity levels due to prevailing climatic conditions. In the KAI area, biomass and grain yield reductions are much higher with the increasing water salinity for both varieties. For example, the biomass and yield reductions in the KAI area were about 40%with an increase in salinity from 5 to 10 and 15 dS m −1 . For the sensitive genotype, the biomass and yield reductions in the KAI area would be above 72% with a similar increase in the salinity levels. Although overall biomass and grain yields in the MED area were lower than in the KAI area, biomass and yield reductions for the salt-tolerant barley variety were only 16% and 8%, with an increase in salinity from 5 to 15 dSm −1 , respectively. However, for the sensitive genotype, reductions in biomass and yield were 12% and 43%, respectively, with a similar increase in salinity levels. Similar trends are obtained for water productivities.
Without salt stress, both varieties have the same performance. However, the tolerant variety performs better than the sensitive variety under salt stress. This is because it has better potential. Therefore, farmers can grow both varieties in the rainfed areas of BEJ, while, in KAI and MED areas where irrigation is necessary for crop growth, the salt-tolerant barley variety should be preferred. The cultivation of the salt-sensitive barley variety in the MED area will be risky, as the yields will be low, and the development of soil salinity over time will remain a challenge. This situation will be very critical for long-term sustainable crop production in the area.

Economic Productivity of Barley Varieties under Different Climatic Conditions
The economic productivity of two barley varieties was estimated using the average unit cost of one water cubic meter in Tunisia and the water use predicted by AquaCrop. The results show that the production cost of 1 kg of barley is lowest in the BEJ area compared to those areas where it is irrigated with saline water.
In the KAI region, the cost will be reduced by 13.28% 28.72% and 47.19% by growing the tolerant variety irrigated with saline water of 5, 10, and 15 dS m −1 , respectively. In the arid region of MED, the benefit will be reduced by 40%, 38%, and 49% by growing the tolerant barley variety by irrigating with saline water of 5, 10 and 15 dS m −1 , respectively ( Figure 6). However, in the sub-humid region of BEJ, there is no significant difference between susceptible and tolerant genotypes. The results show the economic interest for arid region farmers to grow the tolerant barley variety. This stresses the need for appropriate breeding programs for the saline environments for optimizing crop production instead of targeting potential yields.

Discussion
We evaluated the AquaCrop model for two barley varieties under contrasting environments and different water salinity levels. The simulated model values were close to the field measurements concerning biomass, yield and soil salinity. ME and R 2 parameters were close to 0.9, showing the model's ability to simulate the behavior of sensitive and resistant cultivars in contrasting environments and irrigation practices. Araya et al. [5] reported R 2 values of 0.80 when simulating barley biomass and grain yield using AquaCrop. El Mokh et al. [25] reported R 2 values of 0.88 when simulating barley yield under different irrigation regimes in a dry environment using AquaCrop. Mondal et al. [24] reported a 0.12 t ha −1 root mean square error after simulating the yield response of rice to salinity stress with the AquaCrop model. Our results also show a correct prediction with an RMSE of 0.45 t ha −1 (Table 5). This shows that the AquaCrop model simulates biomass production for all environments with an acceptable accuracy level.
AquaCrop model produces consistent simulation results for CC with an R 2 of 0.89 and RMSE of 2.25 ( Table 5). The model also simulated soil salinity satisfactorily for all environments (R 2 = 0.96) for all

Discussion
We evaluated the AquaCrop model for two barley varieties under contrasting environments and different water salinity levels. The simulated model values were close to the field measurements concerning biomass, yield and soil salinity. ME and R 2 parameters were close to 0.9, showing the model's ability to simulate the behavior of sensitive and resistant cultivars in contrasting environments and irrigation practices. Araya et al. [5] reported R 2 values of 0.80 when simulating barley biomass and grain yield using AquaCrop. El Mokh et al. [25] reported R 2 values of 0.88 when simulating barley yield under different irrigation regimes in a dry environment using AquaCrop. Mondal et al. [24] reported a 0.12 t ha −1 root mean square error after simulating the yield response of rice to salinity stress with the AquaCrop model. Our results also show a correct prediction with an RMSE of 0.45 t ha −1 (Table 5). This shows that the AquaCrop model simulates biomass production for all environments with an acceptable accuracy level.
AquaCrop model produces consistent simulation results for CC with an R 2 of 0.89 and RMSE of 2.25 ( Table 5). The model also simulated soil salinity satisfactorily for all environments (R 2 = 0.96) for all situations. The R 2 values exceeding 0.8 are considered excellent for model performance [39]. The ability of AquaCrop to predict yield depends on the appropriate calibration of the canopy cover curve [1,40]. Indeed, after simulation of soil water balance at a daily time step, the model simulates CC and then simulates the transpiration of a crop, biomass above the soil, and converts biomass into yield. Therefore, it is essential to make accurate predictions of the canopy cover by the proper calibration of crop traits. Therefore, through proper calibration, models can be used for additional solutions for the quantification of salinity build-up in the root zone [41].
We also noted the overestimation of the soil salinity at the end of the growing season when saline water is used for irrigation ( Figure 5). This could be due to the excessive leaching of salts from the soil profile through irrigation, as reported by Mohammadi et al. [42]. Over-or underestimation at the end of the season could be the simplification of soil salt transport calculations in the model based on some empirical functions, including the parameters of Ks and the drainage coefficient for vertical downward salt movement. Furthermore, the occasional leaching of salts from the root zone using relatively better-quality water is also recommended. Changing cropping patterns is also a useful strategy for the rehabilitation and management of saline soils, especially when only saline water is available for irrigation.
The AquaCrop model was also capable of predicting water productivity under sub-humid, semi-arid, and arid environments and the effect of salinity. Plants subjected to salinity stress show a varying response in WP. The sensitive genotype was more exposed to varying responses in WP. Besides, heat stress induced by increased temperatures and the water deficit also decreases productivity, as demonstrated by Hatfield [43]. The observed and predicted water productivities were directly affected by climate aridity and the salinity of the irrigation water. However, the tolerant barley variety was less affected by these factors. These results are in agreement with the earlier studies [16,44].
Water scarcity is already hampering agricultural production in the MENA region. Therefore, the adoption of integrated management strategies will be useful for growing tolerant genotypes under saline water conditions and increasing the water use efficiency. For the sustainable management of crop growth in saline environments, soil-crop-water management interventions consistent with site-specific conditions need to be adopted [41]. These may include cyclic or conjunctive saline water use and freshwater through proper irrigation scheduling to avoid salinity development.
There are several traits available for screening genetic material for enhanced production and WP under different climate scenarios. This study shows that, under different water salinity conditions, sensitive barley genotype is more affected by the increasing water salinity than the tolerant barley genotype. The crop yields for both genotypes under all water salinity levels were higher in KAI area compared to the MED area. Therefore, this study recommends that farmers with higher salinity water for irrigation should grow tolerant barley genotypes, allowing them to reduce the cost, on average, by 30% ( Figure 6). However, from a sustainability point of view, irrigation amounts should be kept to a minimum to optimize crop yields instead of targeting potential yields [45]. This exercise will help there be less accumulation of salts in the root zone. Besides, the occasional leaching of salts from the root zone using relatively better-quality water is also recommended. Changing cropping patterns is also regarded as a useful strategy for the rehabilitation and management of saline soils, especially when only saline water is available for irrigation [46,47].

Conclusions
The AquaCrop model with a salinity module was used to evaluate the agronomic performance of two barley varieties for the three different agro-climatic zones in Tunisia. These zones represent sub-humid, semi-arid, and arid climates. The model was calibrated and evaluated using field data from two years (2012 and 2014). The excellent correlation between the simulated and measured data of biomass, yield, and soil salinity confirms the ability of AquaCrop model to simulate crop growth under different climatic conditions. The scenario results using the calibrated model indicate that farmers with higher salinity water for irrigation should grow tolerant barley genotypes. However, from a sustainability point of view, irrigation amounts should be kept to a minimum to optimize crop yields instead of targeting potential yields.