Next Article in Journal
The Effects of Rainfall Intensities and Duration on SCS-CN Model Parameters under Simulated Rainfall
Next Article in Special Issue
B-Spline Method of Lines for Simulation of Contaminant Transport in Groundwater
Previous Article in Journal
Utilization of Steel Slag in Blind Inlets for Dissolved Phosphorus Removal
Previous Article in Special Issue
Influence of Streambed Heterogeneity on Hyporheic Flow and Sorptive Solute Transport
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Simulation of Water and Salt Dynamics in the Soil Profile in the Semi-Arid Region of Tunisia—Evaluation of the Irrigation Method for a Tomato Crop

1
INRGREF, Laboratory of Rural Engineering, University of Carthage, Ariana 2080, Tunisia
2
INAT, Tunis 1082, Tunisia
3
Department of Environmental Sciences, University of California Riverside, Riverside, CA 92521, USA
4
LHyGES, Universite de Strasbourg/EOST/ENGEES, CNRS, 67084 Strasbourg, France
5
INRAT, Laboratory of Horticulture, University of Carthage, Ariana 2080, Tunisia
*
Authors to whom correspondence should be addressed.
Water 2020, 12(6), 1594; https://doi.org/10.3390/w12061594
Submission received: 4 May 2020 / Revised: 26 May 2020 / Accepted: 28 May 2020 / Published: 3 June 2020

Abstract

:
In Tunisia, water used for irrigation is often saline, increasing the risk of salinization for soils and crops. In this study, an experiment was conducted on a tomato crop cultivated on a silty-clay soil irrigated with three different water qualities: 0, 3.5, and 7 dS·m−1. Experimental data were then used to calibrate and validate the Hydrus-1D model, which simulates water flow and salt transfer in soils. The successfully-calibrated and validated model was then used to study the combined effects of the soil osmotic and soil matrix potentials on root water uptake. The values of the root mean square error (RMSE), the coefficient of determination (CD), the modeling efficiency (EF), and the coefficient of residual mass (CRM) were close to their optimal values for both soil water content and soil electrical conductivity profiles, indicating the reliability of the model to reproduce water and salt dynamics. Relative yields (Yr), indirectly estimated using actual and potential root water uptake (transpiration), indicated that the multiplicative stress response model (using the S-shape model) satisfactorily simulates measured yields and reproduces the effects of irrigation with saline waters on crop yields. An alternative scenario using a reduction of water requirements by 50% was investigated to assess an irrigation method with considerable water savings. As the results show that relative yields, Yr, were only slightly reduced, the crop water requirements estimated by CROPWAT 8.0 must have been overestimated. The variation of the soil salinity in the root zone highlighted a high salinization risk in the short-term when water of 7 dS·m−1 is used for irrigation.

1. Introduction

In arid and semi-arid countries such as Tunisia, the use of marginal waters such as saline waters or treated wastewater is a necessity. However, the use of these waters carries a long-term risk for soils, groundwater, and crops [1,2,3]. Irrigation with saline waters generally results in salinization of the topsoil [4,5,6], and, in the long term, salts leaching under the effects of rainfall can contaminate groundwater [7,8]. On the other hand, multiple authors have shown that the salinization phenomenon in the topsoil can be cyclic, and the effect of salinity on crops can be attenuated with proper management that involves appropriate irrigation for salt leaching [9,10,11]. However, in arid regions, reduction in agricultural water use is part of the water resource management policy. This, requires lower irrigation doses and involves irrigation strategies such as deficit irrigation and partial root-zone drying. If one of these techniques is associated with the use of brackish water, the risk of salinization can become dramatic, especially for a strategic crop like tomatoes (in Tunisia).
Numerical models are important and useful tools for assessing the risk of soil salinization and its adverse effects on crop yields. The Hydrus-1D model [12], widely used for this purpose, simulates the movement of water and solutes by numerically solving the Richards and advection-dispersion equations, respectively. The model also considers the simultaneous effects of the soil osmotic and matrix potentials on root water uptake. The model does not directly simulate crop yield but calculates a relative yield as a ratio of actual and potential root water uptake (transpiration) [13]. Several authors have applied the Hydrus-1D model to study water and salt dynamics in agricultural soils and their effects on crops in semi-arid regions, as well as developed management scenarios and long-term forecasts using this model (e.g., [14,15]).
In this context, the objectives of this study are (a) to calibrate and validate the Hydrus-1D model for simulating water and salt dynamics in a soil cultivated with a tomato crop irrigated with saline waters and (b) to estimate crop water requirements using this calibrated and validated model.

2. Materials and Methods

2.1. Field Experiment and Irrigation Management

The field experiment was carried out in the city of Ariana (36°50′40.791″ N, 10°11′13.795″ E) in Tunisia, a semi-arid Mediterranean region. The Ariana region has a semi-arid Mediterranean climate with mild, wet winters and dry and hot summers. The average annual temperature is 18.7 °C, with the winter minimum of 6 °C and the summer maximum exceeding 45 °C. The annual average rainfall is 450 mm, and the average annual reference evapotranspiration is 1100 mm. The soil at the site is a Fluvisol with a loamy-clay texture. The groundwater level is beneath 8 m. A tomato crop (Rio Grande) (Petoseed Saticoy, California, USA) was cultivated with a planting density of 1750 plants/ha. The crop was planted on 17 April 2017 and harvested on 3 August 2017.
Tomato water requirements were estimated from 10-year climate data (2007–2017). Mean values of meteorological variables were used to estimate reference evapotranspiration using the Penman–Monteith equation [16]. The crop coefficient (Kc) for three crop stages were taken from [16] to calculate actual evapotranspiration. For the initial stage Kc = 0.6, for the medium stage Kc = 1.15, and for the final stage Kc = 0.8. CROPWAT 8.0 software (FAO, Rome, Italy) [17] was used to establish irrigation scheduling. CROPWAT 8.0 estimated water requirements of the tomato crop to be 720 mm.
Surface irrigation was used to deliver water to plants (30 plants for each treatment) in three qualities:
  • Freshwater (FW) with a salinity of 0 dS·m−1;
  • Saline water with a salinity of 3.5 dS·m−1;
  • Highly saline water with a salinity of 7 dS·m−1.

2.2. The Hydrus-1D Model

2.2.1. Model Inputs

1.
Soil Hydraulic Properties
The soil hydraulic properties were measured at the Laboratory of Rural Engineering (INRGREF, Ariana, Tunisia) using soil samples collected in March 2017 at depths of 20, 40, 60, and 80 cm, with three replicates for each depth. The pressure plate apparatus [18] was used to measure the volumetric water content at multiple pressure heads for each soil sample. The collected experimental data were fitted to the van Genuchten retention model [19] using the RETC software (U.S. Department of Agriculture, Riverside, USA) [20] to estimate the residual water content θr, the saturated water content θs, and the two shape parameters α and n. Saturated hydraulic conductivity (Ks) was determined by a falling head permeameter method. The values of these parameters for different soil depths are given in Table 1.
2.
Solute Transport Parameters
Soil dispersivity λ (cm) was estimated using laboratory experiments on undisturbed soil columns [21]. The distribution coefficient Kd for total dissolved solutes was measured using batch experiments [22]. The bulk density (BD) was measured for each soil layer using the cylinder method. The measured solute transport parameters for different depths are summarized in Table 2.
3.
Crop Measurements and Models
The plant height, rooting depth, and the leaf area index (LAI) were measured every ten days from the planting date to the harvest date with three replicates for each parameter (Table 3). The LAI was estimated by the analysis of the scanned leaves (5 to 10) images using Mesurim Pro [23]. The measured values of the relative crop yield were determined by a direct weighting using an accurate digital balance.
The S-shape [24] and threshold and slope models [25] were used to simulate the effect of the soil matrix and osmotic potentials on root water uptake, respectively. It was considered that the effects of saturation and salinity were multiplicative. The threshold and slope parameters were 4.59 dS·m−1 and 5%, respectively, as suggested by [25], and the two parameters of the S-shape model were P50 = −800 cm and P3 = 3.

2.2.2. Calibration and Validation

Soil sampling was carried out on the first day (initial) of the experiment, and then after 30, 50, 70, and 109 days (final) at depths of 20, 40, 60, and 80 cm. Soil water content and salt content were measured using the gravimetric method and the saturated paste extract, respectively.
Hydrus-1D calibration was performed for the experiment carried out with fresh water. The two parameters α and n were slightly adjusted during the water movement calibration process and Kd during the solute transport calibration process; all three parameters were calibrated for the two soil layers of 20–40 cm and 40–60 cm. The simulation duration was 109 days, and the output dates corresponded with sampling days. Irrigation/rainfall and evapotranspiration were used to define the atmospheric boundary condition for the soil surface (Figure 1), and a free drainage boundary condition was used at the bottom boundary. Only two rainfall events were recorded on days 49 and 50, with a total amount of 40 mm. The concentration flux boundary condition was used at both boundaries for solute transport. Hydrus-1D validation was performed using the two other experiments with saline waters. The same parameters obtained by calibration for experiments with freshwater were used as input during validation, except for the irrigation water quality.

2.2.3. Statistical Evaluation

The Hydrus-1D results were evaluated graphically by plotting measured and simulated soil water content and soil salinity as a function of soil depth and statistically by calculating the root mean square error (RMSE), the coefficient of determination (CD), the modeling efficiency (EF), and the coefficient of residual mass (CRM) [26]. RMSE values show how much the simulations under- or overestimated the measurements. CD represents the ratio between the scatter of simulated values to the average measured value. The EF value compares simulated values to the average measured value. A negative EF value indicates that the average measured value gives a better estimate than simulated values. CRM is a measure of the tendency of the model to overestimate or underestimate measurements. Positive CRM values indicate that the model underestimates measurements, and negative CRM values indicate a tendency to overestimate measurements. For a perfect fit between observed and simulated data, values of RMSE, CD, EF, and CRM should equal 0, 1, 1, and 0, respectively.
R M S E   % = i = 1 n s i m i 2 n m ¯ × 100
C D = i = 1 n m i m ¯ 2 i = 1 n s i m ¯ 2
E F = 1 i = 1 n s i m i 2 i = 1 n m i m ¯ 2
C R M = i = 1 n s i m i n · m ¯
where si are simulated values, mi are measured values, m ¯ is the average value of observed data, and n is the number of observations.

3. Results

3.1. Water and Salts Movement

Figure 2 and Figure 3 show variations in soil water content and salts profile, respectively, for the three treatments with fresh and saline waters. The soil water content profiles show that:
-
Soil water contents increased with time in all soil horizons compared to the initial soil water contents due to irrigation events.
-
The soil water content in the soil surface layer (0–40 cm) reached its maximum value after 50 days, and then it gradually decreased. In the treatments with saline water, the soil water content reached higher values, especially in the 20–40 cm layer.
-
Variations in the soil water content are lower in the lower part of the soil profile (40–80 cm). Soil water contents reached 0.45 m3·m−3 in the three treatments.
The soil salts profiles measured in the three treatments show that:
-
The salts were continuously leached from the soil profile irrigated with fresh water. Desalination was more important in the surface layer than in deeper layers.
-
The opposite phenomenon was observed in the other two treatments with saline waters, i.e., continuous accumulation of salts in the soil profile. Soil salinity reached 6.7 dS·m−1 at the soil surface and 4.6 dS·m−1 deeper in the soil profile of the treatment with irrigation water quality of 3.5 dS·m1. Maximum salinity values of 11 dS·m−1 and 6 dS·m−1 were, respectively, measured at the soil surface and the bottom of the soil profile in the treatment with irrigation water quality of 7 dS·m−1.

3.2. Hydrus-1D Calibration and Validation

3.2.1. Soil Water Profiles

The measured and simulated soil water profiles for the treatment irrigated with freshwater are plotted in Figure 4. The simulated soil water content values are close to the measured values. The Hydrus-1D model tends to slightly underestimate measured water contents, with RMSE (%) values (for the whole profile) of 5.4, 5.3, 2.6, and 4.9%, at 30, 50, 70, and 109 days, respectively. The CD, EF, and CRM indicators are close to their optimal values and demonstrate a good agreement between simulated (using calibrated parameters) and observed water contents.
The calibrated parameters were then used in model validation simulations. The measured and simulated soil water content profiles for the treatments with saline water of 3.5 dS·m−1 and 7 dS·m−1 are depicted in Figure 5 and Figure 6, respectively.
Graphically, the simulated values of soil water contents are close to those measured at each sampling depth and each sampling date in both validation experiments. Slight deviations from the measured values were observed in the layers between 20 and 40 cm on days 70 and 109. These graphical results are confirmed by the RMSE, CD, EF, and CRM values calculated for each profile (Table 4) with statistical indicators having less than optimal values on days 70 and 109, confirming model deviations from measured values.
The simulated cumulative bottom flux, calculated by Hydrus-1D, was 24 cm for the scenario irrigated with FW, 31 cm for the scenario irrigated with the 3.5 dS·m−1 water, and 38 cm for the scenario irrigated with the 7 dS·m−1 water. The leached amounts of water increased because of reduced root water uptake due to the effects of the saline water used for irrigation.

3.2.2. Soil Salts Profiles

The salinity profiles measured during the tomato crop experiment with the irrigation water quality of 0 dS·m−1 are plotted together with simulated values of soil salts contents in Figure 7. Overall, a good fitting is observed between measured and simulated (using calibrated parameters) values of soil salts contents at the four sampling dates. The soil electrical conductivity is generally slightly overestimated by the model, with RMSE values (for all depths per date) of 1.6, 4.3, 5.1, and 9.6%, respectively. The CD, EF, and CRM values indicate that Hydrus-1D was able to reproduce the salts dynamic in the soil satisfactorily.
The validation results for the salts profiles for the two treatments with saline waters are given in Figure 8 and Figure 9. The measured and simulated salts content profiles are similar at each sampling depth and each sampling date. The graphical results are confirmed by the statistical indicators (Table 5).

3.2.3. Crop Yields

The relative crop yield (Yr) cannot be estimated directly by the Hydrus-1D model. According to [8], the relative crop yield (Yr) can be estimated as the ratio between the actual and potential (or maximum) root water uptake (transpiration), and these two values are available in the model. Figure 10 displays measured and simulated relative crop yields as a function of the irrigation water quality (ECw). According to [27], Yr = Ya/Ymax, where Ya is the measured yield, and Ymax is the maximum yield. In our case, Ymax corresponds to the weight obtained in the experiment with freshwater irrigation and is equal to 2.53 kg/plant.
As shown in Figure 10, the relative yield was reduced as the salinity of the irrigation water increases. The model was able to reproduce the observed decrease in relative yield, although not to the same extent as observed. Measured relative yields, Yr, were reduced only slightly when irrigation water with ECw = 3.5 dS·m−1 was used and nearly by half with ECw = 7 dS·m−1. In contrast, simulated yields were overestimated by the model in both treatments irrigated with saline waters. The observed effects of the salinity stress were larger than those simulated. The parameters (either the threshold or the slope) of the salinity stress response function would have to be adjusted to obtain a better correspondence between simulated and measured relative yields. Note that the model predicted, due to the saturation stress, a slight reduction in relative yield even for the experiment irrigated with fresh water.

3.3. Evaluation of the Irrigation Method

The success of the calibrated and validated Hydrus-1D model in simulating water and salt dynamics for experiments with the tomato crop under field conditions, also motivated employing the model for evaluating alternative irrigation management strategies.
Figure 11 depicts relative yield, Yr, simulated using Hydrus-1D with the S-shape model for scenarios with three irrigation water qualities when 50% less irrigation water was used. Since the tomato yields decreased only slightly when freshwater was used, the irrigation doses calculated by the CROPWAT software must be grossly overestimated [28,29]. The results of Figure 12 show a substantial reduction in the relative yield when the high salinity water was used for irrigation. These results are in agreement with those of Saad et al. [30], who observed a reduction in tomato fruit yield of 40% using irrigation of 50% of ETc and water of 9 dS·m−1.
Figure 12 shows the simulated root zone soil salinity (ECe) as it increases over the season with the salinity of the irrigation water. ECe exceeds 8 dS·m−1 after three months of irrigation with the most saline water (ECw = 7 dS·m−1), which increases the risk of soil salinization. On the other hand, soil salinity in the root zone remains under 4.5 dS·m−1 when moderately saline water is used (ECw = 3.5 dS·m−1). Such root zone salinity can be tolerable to tomato crops, reducing relative yield by only about 15%.

4. Discussion

One of the major weaknesses of numerical models is the lack of calibration and validation in field conditions [31]. In this study, a field experiment involving water and solute dynamics with root water uptake allowed for the successful calibration and validation of the Hydrus-1D model. The RMSEs for measured and simulated values of water contents and salt concentrations were low and showed the ability of the model to reproduce water and salt profiles, and subsequently to evaluate the effects of irrigation with brackish waters on the soil and tomato crop [32,33,34]. Irrigation with saline water with ECw = 3.5 dS·m−1 can be tolerable to the crop and induces a low soil salinization risk. However, saline water with ECw = 7 dS·m−1 reduces Yr nearly by half and increases the soil salinity in the short-term [35] and could be disastrous for the crop after several irrigation cycles [36]. This risk is aggravated when the leaching fraction is reduced in the case of a deficit irrigation strategy [31]. The use of a drip irrigation system associated with proper management practices such as increasing irrigation frequency [37] or controlled deficit irrigation [38], can mitigate these effects on both soil and crops. Replication of this experiment over at least three years can provide a strong statistical significance of the effects of brackish water on tomato crop yields.
The used modeling approach was based on the simulation of the overall salinity. An alternative geochemical approach that considers the multicomponent solute reactions and transport (as implemented in the Unsatchem module of HYDRUS-1D [12]) for the calculation of soil electrical conductivity may be more appropriate for optimizing the management of irrigation with highly saline waters [11,39].

5. Conclusions

In arid and semi-arid regions like Tunisia, the scarcity of water resources of good quality has led to the use of brackish water in agriculture. Salts accumulated in the soil due to irrigation with these waters have negative consequences for both soil and crop yields. Experiments carried out on a tomato crop irrigated with three water qualities (0 dS·m−1, 3.5 dS·m−1, and 7 dS·m−1) allowed the study the variations of water and salt dynamics in the soil. Observed soil salt profiles indicated progressive soil salinization when the plot was irrigated with water of 7 dS·m−1, with soil salinity reaching a value of 8 dS·m−1 at the soil surface. The effects of irrigation with saline waters on the tomato crop were studied using a modeling approach with Hydrus-1D. The study of the scenario where the crop water requirements were reduced by half showed that, in all three experiments, the crop yields were not significantly affected. The irrigation amounts of tomato crop estimated by CROPWAT 8.0 were overestimated. This overestimation could have been caused by the use of crop coefficients not adapted to tomato under the local conditions of the current study.

Author Contributions

Conceptualization, S.K. and B.B.N.; data curation, S.K., R.I., and H.B.; formal analysis, S.K., J.Š., and S.B.M.; funding acquisition, R.I., B.B.N., and M.A.B.A.; investigation, S.K., I.D., and M.R.; methodology, S.K., R.I., and M.R.; project administration, S.K., M.R., and B.B.N.; resources, S.K., H.B., and M.A.B.A.; software, S.K., S.B.M., and B.B.N.; supervision, B.B.N.; validation, I.D., A.Y., M.R., H.B., and M.A.B.A.; visualization, I.D., J.Š., A.Y., R.I., S.B.M., M.R., H.B,. and M.A.B.A.; writing—original draft, S.K. and J.Š.; writing—review and editing, I.D., J.Š., and A.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors would like to appreciate the valuable comments by anonymous reviewers, which helped to improve the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jahin, H.S.; Abuzaid, A.S.; Abdellatif, A.D. Using multivariate analysis to develop irrigation water quality index for surface water in Kafr El-Sheikh Governorate, Egypt. Environ. Technol. Innov. 2020, 17, 100532. [Google Scholar] [CrossRef]
  2. Bueno, M.; Lendínez, M.L.; Calero, J.; Cordovilla, M.D.P. Salinity responses of three halophytes from inland saltmarshes of Jaén (southern Spain). Flora 2020, 266, 151589. [Google Scholar] [CrossRef]
  3. Hussain, M.I.; Muscolo, A.; Farooq, M.; Ahmad, W. Sustainable use and management of non-conventional water resources for rehabilitation of marginal lands in arid and semiarid environments. Agric. Water Manag. 2019, 221, 462–476. [Google Scholar] [CrossRef]
  4. Liu, B.; Wang, S.; Kong, X.; Liu, X.; Sun, H. Modeling and assessing feasibility of long-term brackish water irrigation in vertically homogeneous and heterogeneous cultivated lowland in the North China Plain. Agric. Water Manag. 2019, 211, 98–110. [Google Scholar] [CrossRef]
  5. Pedrero, F.; Kalavrouziotis, I.; Cabañero, J.J.A.; Koukoulakis, P.; Asano, T. Use of treated municipal wastewater in irrigated agriculture—Review of some practices in Spain and Greece. Agric. Water Manag. 2010, 97, 1233–1241. [Google Scholar] [CrossRef]
  6. Ayers, R.S.; Westcot, D.W. Water Quality for Agriculture; Irrigation and Drainage, No. 29; FAO: Rome, Italy, 1985. [Google Scholar]
  7. Golshan, M.; Colombani, N.; Mastrocicco, M. Assessing Aquifer Salinization with Multiple Techniques along the Southern Caspian Sea Shore (Iran). Water 2018, 10, 348. [Google Scholar] [CrossRef] [Green Version]
  8. Kanzari, S.; Ben Nouna, B.; Ben Mariem, S.; Rezig, M. Hydrus-1D model calibration and validation in various field conditions for simulating water flow and salts transport in a semi-arid region of Tunisia. Sustain. Environ. Res. 2018, 28, 350–356. [Google Scholar] [CrossRef]
  9. Azad, N.; Behmanesh, J.; Rezaverdinejad, V.; Abbasi, F.; Navabian, M. Developing an optimization model in drip fertigation management to consider environmental issues and supply plant requirements. Agric. Water Manag. 2018, 208, 344–356. [Google Scholar] [CrossRef]
  10. Yang, H.; Du, T.; Mao, X.; Ding, R.; Shukla, M.K. A comprehensive method of evaluating the impact of drought and salt stress on tomato growth and fruit quality based on EPIC growth model. Agric. Water Manag. 2019, 213, 116–127. [Google Scholar] [CrossRef]
  11. Phogat, V.; Mallants, D.; Cox, J.; Šimůnek, J.; Oliver, D.; Awad, J. Management of soil salinity associated with irrigation of protected crops. Agric. Water Manag. 2020, 227, 105845. [Google Scholar] [CrossRef]
  12. Šimůnek, J.; Genuchten, M.T.; Šejna, M. Recent Developments and Applications of the HYDRUS Computer Software Packages. Vadose Zone J. 2016, 15, 25. [Google Scholar] [CrossRef] [Green Version]
  13. Oster, J.; Letey, J.; Vaughan, P.; Wu, L.; Qadir, M. Comparison of transient state models that include salinity and matric stress effects on plant yield. Agric. Water Manag. 2012, 103, 167–175. [Google Scholar] [CrossRef]
  14. He, K.; Yang, Y.; Yang, Y.; Chen, S.; Hu, Q.; Liu, X.; Gao, F. HYDRUS Simulation of Sustainable Brackish Water Irrigation in a Winter Wheat-Summer Maize Rotation System in the North China Plain. Water 2017, 9, 536. [Google Scholar] [CrossRef] [Green Version]
  15. Li, H.; Yi, J.; Zhang, J.; Zhao, Y.; Si, B.; Hill, R.; Cui, L.; Liu, X. Modeling of Soil Water and Salt Dynamics and Its Effects on Root Water Uptake in Heihe Arid Wetland, Gansu, China. Water 2015, 7, 2382–2401. [Google Scholar] [CrossRef] [Green Version]
  16. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop Evapotranspiration—Guidelines for Computing Crop Water Requirements; FAO Irrigation and Drainage Paper 56; FAO: Rome, Italy, 1998. [Google Scholar]
  17. Smith, M. CROPWAT: A Computer Program for Irrigation Planning and Management; No. 46; FAO: Rome, Italy, 1992; p. 126. [Google Scholar]
  18. Richards, L.A.; Fireman, M. Pressure-plate apparatus for measuring moisture sorption and transmission by soils. Soil Sci. 1943, 56, 395–404. [Google Scholar] [CrossRef]
  19. Genuchten, M.T. A Closed-form Equation for Predicting the Hydraulic Conductivity of Unsaturated Soils. Soil Sci. Soc. Am. J. 1980, 44, 892–898. [Google Scholar] [CrossRef] [Green Version]
  20. Van Genuchten, M.T.; Leij, F.T.; Yates, S.R. The RETC Code for Quantifying the Hydraulic Functions of Unsaturated Soils; U.S. Environmental Protection Agency, Robert S. Kerr Environmental Research Laboratory: Ada, OK, USA, 1991. [Google Scholar]
  21. Kanzari, S.; Bouhlila, R. Laboratory method for estimating solute transport parameters of an unsaturated soil. Am. J. Geochem. Geophys. 2015, 1, 151–156. [Google Scholar]
  22. Kanzari, S.; Sahraoui, H. Laboratory method for characterization of soil/water adsorption coefficient. Exper. J. 2014, 29, 1952–1956. [Google Scholar]
  23. Madre, J.F. Logiciel Mesurim Pro v3.4. Available online: www.dosadi.com.2013 (accessed on 29 May 2020).
  24. Van Genuchten, M.T. A Numerical Model for Water and Solute Movement in and Below the Root Zone; Research Report No. 121; U.S. Salinity Laboratory, USDA, ARS: Riverside, CA, USA, 1987. [Google Scholar]
  25. Van Straten, G.; De Vos, A.; Rozema, J.; Bruning, B.; Van Bodegom, P. An improved methodology to evaluate crop salt tolerance from field trials. Agric. Water Manag. 2019, 213, 375–387. [Google Scholar] [CrossRef]
  26. Loague, K.; Green, R.E. Statistical and graphical methods for evaluating solute transport models: Overview and application. J. Contam. Hydrol. 1991, 7, 51–73. [Google Scholar] [CrossRef]
  27. Steduto, P.; Hsiao, T.C.; Fereres, E.; Raes, D. Crop Yield Response to Water; FAO Irrigation and Drainage Paper No. 66; FAO: Rome, Italy, 2012; p. 505. [Google Scholar]
  28. Etissa, E.; Dechassa, N.; Alemayehu, Y. Estimation of Yield Response (Ky) and Validation of CropWat for Tomato under Different Irrigation Regimes. Irrig. Drain. Syst. Eng. 2016, 5, 167. [Google Scholar] [CrossRef] [Green Version]
  29. Kifle, T. Evaluation of tomato response to deficit irrigation at Humbo Woreda, Ethiopia. J. Nat. Sc. Res. 2018, 6, 57–68. [Google Scholar]
  30. Saad, A.F.; Shalaby, A.A.; Mokhtar, A.M. Influence of Deficit Irrigation Using Saline Water on Yield of Tomato under Two Irrigation Systems. Alex. Sci. Exch. J. 2018, 39, 35–47. [Google Scholar] [CrossRef] [Green Version]
  31. Bastiaanssen, W.; Allen, R.; Droogers, P.; D’Urso, G.; Steduto, P. Twenty-five years modeling irrigated and drained soils: State of the art. Agric. Water Manag. 2007, 92, 111–125. [Google Scholar] [CrossRef]
  32. Gawad, G.A.; Arslan, A.; Gaihbe, A.; Kadouri, F. The effects of saline irrigation water management and salt tolerant tomato varieties on sustainable production of tomato in Syria (1999–2002). Agric. Water Manag. 2005, 78, 39–53. [Google Scholar] [CrossRef]
  33. Malash, N.; Flowers, T.; Ragab, R. Effect of irrigation systems and water management practices using saline and non-saline water on tomato production. Agric. Water Manag. 2005, 78, 25–38. [Google Scholar] [CrossRef]
  34. Reca, J.; Trillo, C.; Sánchez, J.; Martínez, J.; Valera, D. Optimization model for on-farm irrigation management of Mediterranean greenhouse crops using desalinated and saline water from different sources. Agric. Syst. 2018, 166, 173–183. [Google Scholar] [CrossRef]
  35. Du, Y.-D.; Cao, H.-X.; Liu, S.-Q.; Gu, X.; Cao, Y.-X. Response of yield, quality, water and nitrogen use efficiency of tomato to different levels of water and nitrogen under drip irrigation in Northwestern China. J. Integr. Agric. 2017, 16, 1153–1161. [Google Scholar] [CrossRef] [Green Version]
  36. Minhas, P.; Ramos, T.B.; Ben-Gal, A.; Pereira, L.S. Coping with salinity in irrigated agriculture: Crop evapotranspiration and water management issues. Agric. Water Manag. 2020, 227, 105832. [Google Scholar] [CrossRef]
  37. Liu, H.; Duan, A.-W.; Li, F.; Sun, J.-S.; Wang, Y.-C.; Sun, C.-T. Drip Irrigation Scheduling for Tomato Grown in Solar Greenhouse Based on Pan Evaporation in North China Plain. J. Integr. Agric. 2013, 12, 520–531. [Google Scholar] [CrossRef]
  38. Valcárcel, M.; Lahoz, I.; Campillo, C.; Martí, R.; Leiva-Brondo, M.; Roselló, S.; Cebolla, J. Controlled deficit irrigation as a water-saving strategy for processing tomato. Sci. Hortic. 2020, 261, 108972. [Google Scholar] [CrossRef]
  39. Phogat, V.; Mallants, D.; Cox, J.W.; Šimůnek, J.; Oliver, D.P.; Pitt, T.; Petrie, P.R. Impact of long-term recycled water irrigation on crop yield and soil chemical properties. Agric. Water Manag. 2020, 237, 106167. [Google Scholar] [CrossRef]
Figure 1. Daily evapotranspiration, rainfall, and irrigation during the experiments.
Figure 1. Daily evapotranspiration, rainfall, and irrigation during the experiments.
Water 12 01594 g001
Figure 2. Measured soil water content profiles for treatments with (a) fresh, (b) saline (3.5 dS·m−1), and (c) saline (7 dS·m−1) water at 0, 30 50, 70, and 109 days.
Figure 2. Measured soil water content profiles for treatments with (a) fresh, (b) saline (3.5 dS·m−1), and (c) saline (7 dS·m−1) water at 0, 30 50, 70, and 109 days.
Water 12 01594 g002
Figure 3. Measured soil salts profiles for treatments with (a) fresh, (b) saline (3.5 dS·m−1), and (c) saline (7 dS·m−1) water at 0, 30 50, 70, and 109 days. Different scales are used on the horizontal axis.
Figure 3. Measured soil salts profiles for treatments with (a) fresh, (b) saline (3.5 dS·m−1), and (c) saline (7 dS·m−1) water at 0, 30 50, 70, and 109 days. Different scales are used on the horizontal axis.
Water 12 01594 g003
Figure 4. Measured and simulated soil water content profiles for the treatment with the irrigation water quality of 0 dS·m−1 at 30, 50, 70, and 109 days.
Figure 4. Measured and simulated soil water content profiles for the treatment with the irrigation water quality of 0 dS·m−1 at 30, 50, 70, and 109 days.
Water 12 01594 g004
Figure 5. Measured and simulated soil water content profiles for the treatment with the irrigation water quality of 3.5 dS·m−1 at 30, 50, 70, and 109 days.
Figure 5. Measured and simulated soil water content profiles for the treatment with the irrigation water quality of 3.5 dS·m−1 at 30, 50, 70, and 109 days.
Water 12 01594 g005
Figure 6. Measured and simulated soil water content profiles for the treatment with the irrigation water quality of 7 dS·m−1 at 30, 50, 70, and 109 days.
Figure 6. Measured and simulated soil water content profiles for the treatment with the irrigation water quality of 7 dS·m−1 at 30, 50, 70, and 109 days.
Water 12 01594 g006
Figure 7. Measured and simulated soil salts profiles for the treatment with the irrigation water quality of 0 dS·m−1 at 30, 50, 70, and 109 days.
Figure 7. Measured and simulated soil salts profiles for the treatment with the irrigation water quality of 0 dS·m−1 at 30, 50, 70, and 109 days.
Water 12 01594 g007
Figure 8. Measured and simulated soil salts profiles for the experiment with the irrigation water quality of 3.5 dS·m−1 at 30, 50, 70, and 109 days.
Figure 8. Measured and simulated soil salts profiles for the experiment with the irrigation water quality of 3.5 dS·m−1 at 30, 50, 70, and 109 days.
Water 12 01594 g008
Figure 9. Measured and simulated soil salts profiles for the experiment with the irrigation water quality of 7 dS·m−1 at 30, 50, 70, and 109 days.
Figure 9. Measured and simulated soil salts profiles for the experiment with the irrigation water quality of 7 dS·m−1 at 30, 50, 70, and 109 days.
Water 12 01594 g009
Figure 10. The measured and simulated tomato relative yields (Yr) as a function of irrigation water electrical conductivity (ECw).
Figure 10. The measured and simulated tomato relative yields (Yr) as a function of irrigation water electrical conductivity (ECw).
Water 12 01594 g010
Figure 11. The simulated relative yields (Yr) of tomato for scenarios with different electrical conductivities of irrigation water (ECw) and with a 50% irrigation reduction.
Figure 11. The simulated relative yields (Yr) of tomato for scenarios with different electrical conductivities of irrigation water (ECw) and with a 50% irrigation reduction.
Water 12 01594 g011
Figure 12. Simulated rootzone soil salinity for scenarios with different water qualities and a 50% reduction in irrigation water.
Figure 12. Simulated rootzone soil salinity for scenarios with different water qualities and a 50% reduction in irrigation water.
Water 12 01594 g012
Table 1. Soil hydraulic parameters [19] for different soil depths.
Table 1. Soil hydraulic parameters [19] for different soil depths.
Depth (cm)θr (cm3·cm−3)θs (cm3·cm−3)α (cm−1)nKs (cm·d−1)
0–200.010.500.0121.2312.60
20–400.090.480.0151.328.90
40–600.090.500.0101.4212.60
60–800.090.470.0121.3912.50
Table 2. Solute transport parameters for different soil depths.
Table 2. Solute transport parameters for different soil depths.
Depth (cm)BD (g·cm−3)λ (cm)Kd (cm3·g−1)
0–201.435.800.14
20–401.465.400.20
40–601.394.600.30
60–801.512.800.35
Table 3. Tomato crop parameters for each water quality scenario.
Table 3. Tomato crop parameters for each water quality scenario.
Day after Planting
Water QualityParameter10203040506070
FreshwaterCrop Height (cm)56914333538
LAI0.35122.52.752.82.85
Root Depth (cm)5101520253035
3.5 dS·m−1Crop Height (cm)4.55.5810182225
LAI0.350.981.982.252.22.22.2
Root Depth (cm)481216202428
7 dS·m−1Crop Height (cm)461011131519
LAI0.350.961.962.21.81.51.5
Root Depth (cm)471114182125
Table 4. The statistical indicators for measured and simulated soil water content profiles at different sampling dates for treatments irrigated with saline waters.
Table 4. The statistical indicators for measured and simulated soil water content profiles at different sampling dates for treatments irrigated with saline waters.
Water Quality30 Days50 Days70 Days109 Days
RMSE (%)0 dS·m−15.405.302.604.90
3.5 dS·m15.003.308.107.00
7 dS·m15.104.3010.407.50
CD0 dS·m10.680.670.680.95
3.5 dS·m10.911.623.332.60
7 dS·m11.130.695.031.81
EF0 dS·m10.950.850.50.65
3.5 dS·m10.960.93−0.280.77
7 dS·m10.97−0.07−0.210.4
CRM0 dS·m10.020.010.020.04
3.5 dS·m1−0.020.010.010.02
7 dS·m10.020.04−0.020.05
Table 5. The statistical indicators for measured and simulated soil salts profiles at different sampling dates for treatments irrigated with saline waters.
Table 5. The statistical indicators for measured and simulated soil salts profiles at different sampling dates for treatments irrigated with saline waters.
Water Quality30 Days50 Days70 Days109 Days
RMSE (%)0 dS·m15.405.302.604.90
3.5 dS·m15.003.308.107.00
7 dS·m15.104.3010.407.50
CD0 dS·m10.910.870.970.91
3.5 dS·m10.980.920.980.90
7 dS·m11.051.001.000.97
EF0 dS·m10.990.990.990.99
3.5 dS·m10.990.970.990.98
7 dS·m10.990.990.990.99
CRM0 dS·m1−0.01−0.03−0.05−0.09
3.5 dS·m10.040.05−0.010.01
7 dS·m10.010.010.01−0.01

Share and Cite

MDPI and ACS Style

Kanzari, S.; Daghari, I.; Šimůnek, J.; Younes, A.; Ilahy, R.; Ben Mariem, S.; Rezig, M.; Ben Nouna, B.; Bahrouni, H.; Ben Abdallah, M.A. Simulation of Water and Salt Dynamics in the Soil Profile in the Semi-Arid Region of Tunisia—Evaluation of the Irrigation Method for a Tomato Crop. Water 2020, 12, 1594. https://doi.org/10.3390/w12061594

AMA Style

Kanzari S, Daghari I, Šimůnek J, Younes A, Ilahy R, Ben Mariem S, Rezig M, Ben Nouna B, Bahrouni H, Ben Abdallah MA. Simulation of Water and Salt Dynamics in the Soil Profile in the Semi-Arid Region of Tunisia—Evaluation of the Irrigation Method for a Tomato Crop. Water. 2020; 12(6):1594. https://doi.org/10.3390/w12061594

Chicago/Turabian Style

Kanzari, Sabri, Issam Daghari, Jiří Šimůnek, Anis Younes, Riadh Ilahy, Sana Ben Mariem, Mourad Rezig, Béchir Ben Nouna, Hassouna Bahrouni, and Mohamed Ali Ben Abdallah. 2020. "Simulation of Water and Salt Dynamics in the Soil Profile in the Semi-Arid Region of Tunisia—Evaluation of the Irrigation Method for a Tomato Crop" Water 12, no. 6: 1594. https://doi.org/10.3390/w12061594

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop