Modelling Salinity and Sodicity Risks of Long-Term Use of Recycled Water for Irrigation of Horticultural Crops

: Long-term use of recycled water (RW) for irrigation in arid and semiarid regions usually changes the soil solution composition and soil exchange characteristics, enhancing the risk for salinity and sodicity hazards in soils. This modelling study focuses on developing alternative management options that can reduce the potentially harmful impacts of RW use on the irrigation of wine grapes and almonds. The multicomponent UNSATCHEM add-on module for HYDRUS-1D was used to evaluate the impact of long-term (2018–2050) use of irrigation waters of different compositions: good-quality low-salinity (175 mg/L) water (GW), recycled water with 1200 mg/L salinity (RW), blended water of GW and RW in the 1:1 proportion (B), and monthly (Alt1) and half-yearly (Alt6) alternate use of GW and RW. The management options include different levels of annual gypsum applications (0, 1.7, 4.3, and 8.6 t/ha soil) to the calcareous (Cal) and hard red-brown (HRB) soils occurring in the Northern Adelaide Plain (NAP) region, South Australia. Additional management scenarios involve considering different leaching fractions (LF) (0.2, 0.3, 0.4, and 0.5) to reduce the salinity build-up in the soil. A new routine in UNSATCHEM to simulate annual gypsum applications was developed and tested for its applicability for ameliorating irrigation-induced soil sodicity. The 1970–2017 period with GW irrigation was used as a warmup period for the model. The water quality was switched from 2018 onwards to reﬂect different irrigation water qualities, gypsum applications, and LF levels. The data showed that the GW, B, Alt1, and Alt6 irrigation scenarios resulted in lower soil solution salinity (EC sw ) than the RW irrigation scenario, which led to increased EC sw values (4.1–6.6 dS/m) in the soil. Annual gypsum applications of 1.7, 4.3, and 8.6 t/ha reduced pH, SAR, and ESP in both soils and reduced the adverse impacts of irrigation, especially in surface soils. A combination of water blending or cyclic water use with 3.8 t/ha annual gypsum applications showed promise for the SAR and ESP control. Additionally, irrigation with RW, a 0.2 LF, and annual gypsum applications limited the harmful salinity impacts in the soils. However, in the RW irrigation scenario, EC sw and ESP at the bottom of the crop root zone (90–120 cm depth) in the HRB soil were still higher than the wine grape and almond salinity thresholds. Thus, annual amendment applications, combined with the long-term use of blended water or cyclic use of RW and GW, represent a sustainable management option for crop production at the calcareous and hard red-brown soils.


Introduction
Availability of freshwater for irrigation is diminishing worldwide due to a rapid increase in demand by other sectors (e.g., industrial and urban), erratic rainfall distribution, and dwindling surface and groundwater water resources. Therefore, other sources of irrigation water are increasingly being explored to sustain irrigated agriculture. These waters include wastewaters emanating from public utilities and industries. These waters may contain numerous disease-causing organisms, hazardous chemicals, metal ions, and other impurities. Although direct use of these waters is prohibited, post-treated use for irrigation is not uncommon, especially in arid and semiarid regions, including Australia. There is wide-ranging experience using recycled water for irrigation of different crops involving varied soil textures, climate, and management conditions [1,2].
Generally, the use of recycled water (RW) for irrigation can have both positive and negative impacts on crops and soils, depending on the water quality and the extent and duration of exposure. The benefits include reducing freshwater demand, recycling nutrients, minimizing the discharge of pollutants into waterways [3,4], and improving soil fertility, nutrient uptake, plant growth, and crop yield [4][5][6][7][8][9]. Despite these beneficial impacts, negative influences of long-term use of RW are also common. Therefore, appropriate management of these waters is essential to realize their promise for sustainable use for crop production [10][11][12].
Sustainable use of treated/recycled water (RW) for irrigation refers to developing strategies and management tools to help reduce harmful impacts at the farm and regional scales while maintaining the soil's productive capacity. The main factor affecting the management of RW irrigation is the interaction between the total mass of salts in the soil, salinity of irrigation water (characterized by electrical conductivity, EC), the ratio of sodium and other cations in the solution (represented by the sodium absorption ratio, SAR), and the fraction of the solid phase occupied by sodium cations (characterized by the exchangeable sodium percentage, ESP) in the soil. In recent years, the role of high concentration of K and Mg on the soil exchange and soil solution, typical of recycled municipal wastewater [13], is being reassessed by many researchers. Numerous studies have shown that Mg can have adverse impacts on the permeability of irrigated soils [14][15][16]. The magnitude of impacts depends on soil clay mineralogy, as well as on the relative amounts of K and Mg over Na and Ca in the soil exchange and in irrigation water [14,[17][18][19]. Recognizing the role of K and Mg, Rengasamy and Marchuk [14] proposed a cation ratio of structural stability (CROSS) alternate to SAR. This average value for the Mg coefficient is consistent with the value 0.60 utilized by Rengasamy and Sumner [20]. Zhu et al. [18] found that the numerical coefficient for the K concentration in CROSS was in the range of 0.309-0.777, with an average value of 0.605. These factors control soil salinity, soil structure degradation, water flow, and leaching of soluble salts from the root zone.
Soil salinity impacts plant growth through osmotic and toxicity effects and causes nutrient imbalances in the soil [21,22]. High SARSAR or ESP causes soil sodicity (resulting in soil dispersion/swelling), thus increasing soil resistance and reducing root growth and water movement through the soil [23]. High sodicity reduces the leaching of salts from the root zone and may cause oxygen depletion through waterlogging [24] and restricted drainage conditions in the soil.
The adverse impacts of RW irrigation depend on factors such as the chemical composition of irrigation water, climate parameters, soil texture, vertical soil heterogeneity (e.g., the presence of clay-rich horizons), and the water needs and tolerance levels of crops. However, there is a strong tendency for RW waters to have higher salinity and higher concentrations of sodium (Na) relative to other cations compared to other irrigation waters [25]. To better understand salinity problems, one must realize that every irrigation water adds salts into the soil system. For example, irrigation with 1 ML of good-quality water (0.5 dS/m, similar to rainwater) can add approximately 320 kg of salts into the soil. As plants usually contain no more than 3.5% of salts in the dry matter [26], no more than 350 kg of salts applied with irrigation water can be exported in 10 t of harvestable products. The salt mass in the Soil Syst. 2021, 5, 49 3 of 30 soil increases with each irrigation unless the salts are leached below the crop rooting zone by an excess of water above what the soil can store [26]. Therefore, a minimum leaching fraction (LF) is required in all irrigation systems, regardless of water quality, to remove the inevitable build-up of salts from the rooting zone [27]. If salts accumulated due to RW irrigation are not properly leached, increased osmotic pressures can potentially lead to a severe reduction in the sap flow and water absorption by the plants [28,29] and ultimately have a drastic impact on crop growth and yield [2,24].
On the other hand, high SAR of irrigation water can have harmful impacts on the soil by increasing the soil's ESPESP, which induces severe soil structural (e.g., crusting) and hydraulic (pore blockage and permeability reduction) problems. RWs typically have SAR values ranging from 7.5 to 8.5, which can potentially lead to the development of high ESP in the soil [25,30]. The soil analysis of the NAP region in the current project also indicated that most soils have high ESP, especially in the subsoil [10]. Other studies have also observed significant ill-effects of long-term irrigation with RW on soils' physicochemical properties [11,12,24,[31][32][33]. Increased ESP causes deterioration of soils' physical properties, such as clay dispersion, with a subsequent collapse of soil structure, pore blockage, and eventually, a decrease in soil permeability, leading to waterlogging, poor plant performance, decreased leaching, and salinization [34][35][36]. Intensive irrigation with RW of loam and clay soils has been shown to result in a significant increase in clay dispersion and its eluviation from the upper soil layers [37]. Notably, clay dispersion may be greatly increased by high concentrations of Na + in the soil solution due to Na + dominating the soil exchange complex. In addition, the reduced hydraulic conductivity can result in the increased surface runoff of irrigation water, which may then contaminate other water bodies with both dissolved and suspended materials.
Although salt accumulation and sodicity will inevitably increase with irrigation, this development can be slowed down and even reversed within the crop root zone by applying careful management. However, judicious use for irrigation requires developing site-specific management options to mitigate the harmful impacts [38]. To address RW-induced problems and devise site-specific practical guidelines, there is a need to evaluate multiple long-term strategies to reduce expected harmful impacts on the soil while sustaining crop production. If irrigation water salinity is too high for crop growth, dilution with lower salinity water could be a good approach [39]. Grattan and Rhoades [40] advocated alternate irrigation with fresh and saline waters in preference to blending, suggesting that this improves soil conditions by leaching salts with lower salinity water. Experiments with grapevines and saline waters [41,42] showed that irrigation with non-saline water (rainfall) for two months usually returned soil EC e to pre-saline irrigation conditions. Irrigation with fresh or blended water prior to planting is also a useful option since it provides the seedlings, which are more sensitive to salts [43], with a less saline environment during their establishment. This would also reduce leaching requirements during subsequent crop growth [44]. Balancing fresh, saline, and/or reclaimed water irrigation may be valuable for maintaining soil quality under irrigated crops. However, this should be practiced with caution as this can lead to increased soil degradation problems.
On the other hand, to address high SAR-and ESP-related problems, frequent application of Ca (e.g., gypsum/organics) is essential for reducing the harmful impacts of RW with high exchangeable Na + . Gypsum directly affects the flocculation and aggregation of the subsoil and indirectly improves rooting, which leads to greater aggregation [45]. However, the extent and frequency of such amendments depend on continuous changes in the soil solution due to RW irrigation and inherent sodicity problems. Other options include evaluating the impact of blending RW with good-quality water (GW), harvesting stormwater, or using RW and GW alternatively. Thus, appropriate irrigation water use is essential to reduce the adverse impacts of soil salinity and sodicity in the soil system. This modelling study used the HYDRUS 1D-UNSATCHEM simulator to evaluate the impact of different management scenarios involving long-term irrigation with waters of Soil Syst. 2021, 5, 49 4 of 30 different qualities applied to two texturally different soils and crops with different salt tolerance to sustainably manage water and soil resources. To the best of our knowledge, this is the first time that such a comprehensive set of modelling scenarios has been evaluated, focusing on the impact of annual gypsum applications and leaching fraction to ameliorate the adverse impact of the use of RW water for irrigation. Results will help devise better management strategies for using recycled water in combination with other available water resources.

Description of the Study Area
The study region covers a part of the Northern Adelaide Plains (NAP), South Australia, where a considerable expansion of irrigated horticulture has been initiated, predominantly utilizing available RW from the adjoining Bolivar treatment plant [46]. The wastewater is treated using the "Dissolved Air Flotation and Filtration" (DAFF) technique through a tertiary filtration and disinfection process, making it suitable for irrigating horticultural crops. The NAP has a Mediterranean-type climate characterized by hot, dry summers and cool to cold winters. Long-term (1900-2016) average rainfall in the region amounts to 475 mm [47], and annual evapotranspiration amounts to 1308 mm. Therefore, horticultural orchards require intensive irrigation for profitable production. A more detailed description of the study area can be found in [10].
Average seasonal RW quality data were estimated from raw data from 2002-2017 [48]. The final values used for the modelling study are taken from Phogat et al. [12]. The water quality falls into the medium category of the salinity rating in Australia [49], which means that moderately tolerant crops can be grown using RW as an irrigation source.
Fourteen representative soil profiles were excavated in the study area to determine the relevant physico-chemical properties of the main soil groups. The sampling strategy was focused on ensuring a reasonable representation of major soil groups occurring in the region. More details about the soil sampling sites are described in [50]. The prominent soil groups occurring in the region are hard red-brown (HRB), deep uniform to gradational (DuG), calcareous (Cal), and sand over clay (SOC) textures covering 52, 13, 12, and 10% of the study area, respectively [51]. A number of soil profiles were sampled from each soil group: 6 from HRB, 3 from Cal, 3 from SOC, and 2 from DuG textural groups. Undisturbed (single) and disturbed (triplicate) soil samples were collected from each profile at depths of 0-10, 10-30, and 30-60 cm. Additional samples were collected from greater depths where sampling was not impeded by hard soil layers (60-90 and 90-120 cm). This investigation has put most effort on the calcareous and hard red-brown soil groups, which are the predominant soil groups in the NAP region [51]. Disturbed soil samples were used for particle size analysis, while undisturbed cores were used to measure the bulk density, the water content-matric suction relationship, and the saturated hydraulic conductivity [52]. More details on measuring soil hydraulic and soil chemical properties can be found in Appendix A.

Modelling Scenarios
We used the UNSATCHEM module of HYDRUS-1D [53] to evaluate the long-term dynamics of coupled water flow, solute transport, and major ion geochemistry in the predominant soils in the NAP region under different applications of irrigation water quality and soil amendments [10]. Few studies have evaluated the impact of the long-term use of RW on soil properties and crop yield [12,54]. This study mainly focuses on managing the long-term use of RW, utilizing scarce GW, and their various alternate and blending combinations. The blending mode consists of 1:1 mixing of RW and GW (B), and the RW and GW alternate use comprises monthly (Alt1) and half-yearly cycles (Alt6). The scenarios with different annual gypsum applications of 0, 1.7, 4.3, 8.6, and 12.9 t/ha are designated as G0, G1, G2, G3, and G4, respectively. The various leaching fraction scenarios (LF = 0.2, 0.3, 0.4, and 0.5) are listed as LF0.2, LF0.3, LF0.4, and LF0.5, respectively. The LF0.2 scenario Soil Syst. 2021, 5, 49 5 of 30 was adopted as a base scenario for different irrigation application modes. Wine grapes and almonds were considered as test crops in the calcareous and hard red-brown soils, respectively. Detailed information on individual scenarios for different NAP soils is given in Tables 1 and 2.  Alt1 " 0 * RW = recycled water, # GW = good-quality water, + B (1:1) = blending of RW and GW in a 1:1 ratio, " Alt1 = monthly alternate use of RW and GW, † Alt6 = six-monthly alternate use of RW and GW, HRB = hard red-brown, LF = leaching fraction.
Apart from the scenarios listed in Table 1, five more scenarios involving a much higher annual gypsum application (12.9 t/ha) were performed for hard red-brown (HRB) soils. Three leaching scenarios were carried out to evaluate the impact of different leaching fractions (LF; 0.3, 0.4, 0.5) for recycled water and blending irrigation options.
Gypsum applications are usually estimated based on soil texture, cation exchange capacity, ESP of the soil, soil depth, crop type, water quality, etc. [55][56][57]. The estimated Soil Syst. 2021, 5, 49 6 of 30 amount of gypsum is mixed in the surface soil layer, usually in the 0-20 cm depth. However, the effectiveness of gypsum depends on many factors, such as the purity of gypsum, its granule size, application method, mixing uniformity, etc. [33]. In addition to a control scenario (zero gypsum application), we have considered three typical gypsum application rates (1.72, 4.3, and 8.6 t/ha), which fall within the general recommended range of gypsum applications based on the sodicity (ESP) classes occurring in Australian soils [58]. These classes are slightly sodic (ESP < 6), sodic (ESP = 6-15), and highly sodic (ESP >15) for alkaline soils [59]. However, farmers sometimes resort to much higher amounts of gypsum. Therefore, we have also tested a higher dose of 12.9 t/ha as an annual application, especially in hard red-brown soils ( Table 2), where inherent subsoil constraints, including high sodicity, are a common occurrence.

Model Description
A detailed description of the HYDRUS-1D model is available in Šimůnek et al. [53,60]. The one-dimensional Richards equation, which assumes that the air phase plays an insignificant role in liquid flow processes and that water flow due to thermal gradients can be neglected, is solved in HYDRUS-1D using a Galerkin-type linear finite element scheme [51]. The governing one-dimensional water flow equation is described as: where θ is the soil water content (L 3 L −3 ), t is the time (T); h is the soil water pressure head (L), z is the vertical coordinate (positive upwards), K(h) is the hydraulic conductivity function (LT −1 ), α is the angle between the flow direction and the vertical axis (e.g., 1 for vertical flow and 0 for horizontal flow), and S(h, z, t) is a sink term accounting for water uptake by plant roots (L 3 L −3 T −1 ).
The HYDRUS-1D model is run with the best available soil, crop, climate, and irrigation management parameters. The soil parameters were mostly directly measured on representative soil samples or estimated using auxiliary soil data. Most of the crop parameters are from the literature, although some calibrated parameters were available for wine grapes. Climate data were partly measured on-site (historical climate) and partly derived from downscaled climate projections. Irrigation management parameters were obtained from surveys of irrigators [48]. The model was, thus, not specifically calibrated, as preference was given to the use of the measured input parameters. Note that estimating model parameters through calibration using an inadequate conceptual model could lead to the so-called parameter bias [61].

Soil Hydraulic Characteristics
The soil hydraulic characteristics are modelled using the water retention and hydraulic conductivity functions described by the van Genuchten-Mualem constitutive relationships [62]: where θ s is the saturated water content (L 3 L −3 ), θ r is the residual water content (L 3 L −3 ), K s is the saturated hydraulic conductivity (LT −1 ), l is a shape factor, and m, α, and n are empirical shape parameters, where m = 1 − 1/n. S e is the relative saturation (dimensionless), which is defined as The parameters (θ r , θ s , K S , α, and n) for different soils were derived from the measurements on soil cores collected from soil profiles from the NAP, as described in Section 2.1. Average parameter values were then derived for the two most important textural soil groups (Appendix A, Table A1).

Solute Transport Properties for UNSATCHEM
The major ion chemistry module UNSATCHEM enables the simulation of equilibrium geochemical reactions involving Ca, Mg, Na, K, SO 4 , Cl, NO 3 , H 4 SiO 4 , HCO 3 , and CO 2 [60,63]. The model accounts for equilibrium chemical reactions between these components, such as aqueous complexation, cation exchange, and precipitation-dissolution. Initial average concentrations of soil solution species, soil exchange cations, and estimated Gapon selectivity coefficients for soils collected from the NAP area are given in [9,10]. Note that the initial soil solution is brought into equilibrium with the cation exchange complex during the model warming up period (1970-2018).

Root Water Uptake Parameters
Water extraction S(h, h s , x, z, t) from the soil is computed according to the Feddes macroscopic approach [64]. In this method, the potential transpiration rate, T p , is distributed over the root zone using a normalized root-density distribution function {β(x, z, t)} and multiplied by dimensionless water α 1 (h) and salinity α 1 (h s ) stress response functions as: This model reduces potential plant root water uptake rates according to the local soil water pressure head, h, and osmotic head, h s , at any point in the root zone. This model defines how potential transpiration (T p ) is reduced when the soil is no longer capable of supplying the amount of water required by plants under prevailing climatic and soil conditions. The multiplicative model for the uptake reduction due to the osmotic stress is considered in this study: The reduction of root water uptake due to the water stress, α 1 (h), is described as: where h 1 , h 2 , h 3 , and h 4 are the threshold model parameters. Water uptake is at the potential rate when the pressure head is between h 2 and h 3 , decreases linearly when h > h 2 or h < h 3 , and becomes zero when h < h 4 or h > h 1 . These critical values of the pressure head for viticulture and almonds were taken from previous investigations in South Australia [65][66][67].
The threshold model uses two variables to simulate the osmotic stress α 1 (h s ): the osmotic head, below which water is extracted at the maximum rate, and the slope, i.e., a fractional reduction of water uptake per unit increase in the osmotic head above the threshold. These parameters were obtained from Zhang et al. [68] for wine grapes and from Ayers and Westcott [69] for almonds.

Estimation of Potential Evaporation (E s ) and Potential Transpiration (T p )
The HYDRUS-1D/UNSATCHEM [44] model requires daily inputs of rainfall as well as potential evaporation (E s ) and potential transpiration (T p ). Therefore, the daily E s and T p values for all wine grapes and almonds for current and future climates were estimated following the FAO-56 dual crop coefficient (DCC) approach [70]. The FAO-56 DCC methodology requires a considerable amount of data for crops, soils, and climate. All relevant information and required data for this approach for the NAP soils and crops are given in Mallants et al. [10] and Phogat et al. [12]. Daily climate data for the historical climate (1970-2018) were obtained from the Bureau of Meteorology (BOM), the Edinburgh RAAF station, while future climate (2018-2050) data were taken from the Goyder climate change median climate projections for the same station [71].
The median data are based on the downscaled series obtained from the GFDL-ESM2M Global Climate Model (GCM), one of the six better-performing GCMs. A median decrease in annual rainfall by 2050 is 6.8% (relative to the 1986-2005 baseline), the 10th percentile decrease is 8.8%, and the 95th percentile decrease is 3.5% (for the intermediateemission Representative Concentration Pathway RCP 4.5). RCP 4.5 is a scenario representing long-term global emissions of greenhouse gases, short-lived species, and land-use-cover, which stabilizes radiative forcing at 4.5 W/m 2 (approximately 650 ppm CO 2 -equivalent) in the year 2100 without ever exceeding that value [72]. A single future climate scenario is used rather than a range of future climate scenarios to keep the overall number of modelling scenarios manageable.
Apart from daily E s and T p , the irrigation requirements for wine grape and almonds were also estimated following the FAO-56 dual crop coefficient (DCC) approach [70], as discussed in detail by Mallants et al. [10] and Phogat et al. [12,73]. For wine grapes, measured LAI values were taken from similar studies [67,74] in the study area. In addition, values for canopy geometry and fractional cover for wine grape previously calibrated by Phogat et al. [65] have been used here. The functional relation followed in the FAO dual crop coefficient approach [70] is given as: where ET C is crop evapotranspiration (LT −1 ), ET 0 is reference crop evapotranspiration (LT −1 ), K cb is the basal crop coefficient, representing the plant transpiration component, and K e is the soil evaporation coefficient. Some of the basic information, such as Kcb, crop duration, plant height, rooting depth, and depletion factor, is available from FAO-56 [70]. The details on data requirements and other relevant information for each crop, soil, and climate conditions can be found in Mallants et al. [10]. In this approach, standard K cb values [70] of the crops were adjusted for the local climate, taking into consideration crop height, wind speed, and minimum relative humidity averages. Estimated K cb 's in the current study compare well with previously calibrated values for wine grape [58]. Soil-specific information (soil texture, field capacity, θfc; permanent wilting point, θwp; readily available water, RAW; and total available water, TAW) has been drawn from the soil analyses [50]. The values of daily potential transpiration (T p ) and soil evaporation (E s ) were used as time-variable boundary conditions in the HYDRUS model, along with the irrigation schedule for different crops and precipitation received at the site during the simulation period. The amount and timing of irrigation were also imposed as a time-variable flux boundary.

Modelling Domain, Initial and Boundary Conditions
Simulations for vine grapes and almonds were performed for a 200 cm deep soil profile divided into 100 finite elements. The finite element nodes were distributed so that elements at the soil surface were half the size of the elements at the bottom of the soil profile to ensure fine discretization at the soil surface where major water and solute dynamic processes occur. The measured textural heterogeneity in the Cal and HRB soils Soil Syst. 2021, 5, 49 9 of 30 in the NAP area was maintained by dividing the domain into five soil layers (0-15, 15-30, 30-60, 60-100, and 100-200 cm). At the soil surface, an atmospheric boundary condition with surface runoff was imposed. A free drainage boundary condition was applied at the bottom boundary. For solute transport, a concentration flux boundary condition was used at the soil surface, and a zero-concentration gradient boundary condition was imposed at the bottom to allow gravitational outflow with drainage water. A constant pressure head of-100 cm was considered to initialize the model runs. The measured daily rainfall and calculated daily E s and E p for different crops were applied during the warmup period (1970-2018). The measured soil solution and exchange parameters were assumed as initial conditions for the multicomponent solute transport module. Initially, the model was warmed up for 48 years from July 1970 to June 2018 using the historical climate parameters from the Bureau of Meteorology (BOM) Station No. 023,083 at the Edinburgh RAAF (latitude-34.7111, longitude 138.6222, elevation 17 m) to initialize the soil water contents and to attain equilibrium conditions for chemical species for both soils, crops, and irrigation quality and management scenarios. Subsequently, the model was run using the future climate (2018-2050) projection data for the same station [71]. Soil amelioration doses (Go, G1, G2, G3, and G4) and leaching fractions (LF0.3, LF0.4, and LF0.5) were introduced during the future climate simulations.
More details about the input parameters required for executing various management scenarios using HYDRUS-1D (UNSATCHEM) have been described in our previous publications [11,12] and project reports [10,48,50]. However, it is worth mentioning that we have tested a new subroutine implemented in HYDRUS-1D/UNSATCHEM specifically for this study, i.e., yearly applications of gypsum to manage adverse impacts of RW on the soil properties, especially soils' SAR and ESP. Therefore, simulations were conducted with different annual applications of gypsum to optimize the application dose for different NAP soils for sustainable horticultural production.

Model Calibration
This study was focused on assessing the feasibility of the expansion of irrigated agriculture in the Northern Adelaide Plain, South Australia, by utilizing the increased capacity of recycled water (RW) from the water treatment plant. The RW water delivery and the irrigation system in the primary production area selected for the irrigation expansion are not yet operational. Therefore, conducting a field experiment with RW on various soil groups (or even a single soil group) was not possible. HYDRUS-2D (the two-dimensional version of HYDRUS) has been sufficiently calibrated and validated in the adjoining area (Riverland) for water balance and salinity dynamics under a range of horticultural industries, including almonds, grapevine, and citrus [65][66][67]. These studies provide convincing evidence that the HYDRUS model can generate reliable estimates when field-measured input data are used for simulations. Similarly, HYDRUS-UNSATCHEM has been successfully calibrated and validated in other studies involving a variety of water qualities, crops, and environmental conditions (e.g., [75][76][77][78]). Numerous other published studies [79][80][81][82] involving UNSATCHEM simulations similar to the current study have been able to unravel varied intricacies of the soil-plant-atmosphere system and generate valuable information. Therefore, we strongly believe that in spite of the lack of site-specific model calibration, numerical simulations with measured and reliable input data can produce reliable outcomes.

Results and Discussion
Different soils responded to recycled water irrigation differently, depending on the dynamic interactions between soil solution and exchange properties of the soil. Therefore, management options are discussed below for different soil types.

Changes in Soil pH in the Calcareous Soil
In the calcareous soil, average pH values in the 0-15 cm soil depth ranged from 7.3 to 8.2 under different irrigation quality scenarios (RW, GW, B, Alt1, and Alt6) with no gypsum additions ( Figure 1). The pH values decreased with depth except under good-quality (GW) irrigation. Modelled values are comparable to those estimated with the 0.01 M CaCl 2 solution in the calcareous soil profiles collected from the study area [50]. It is worth noting that the pH values stabilized for the scenarios with 8.6 t/ha annual gypsum (G3) applications. The overall average reduction in pH in various depths, irrespective of water quality and gypsum doses, ranged from 4 to 6.2%. Simulated average pH values were 7.6, 7.5, and 7.3 in the 0-15, 30-60, and 90-120 cm soil depths, respectively.

Salinity Distribution in the Calcareous Soil
Average ECsw in the 0-15, 30-60, and 90-120 cm soil depths in the G0 scenario were lower than the salinity threshold of wine grapes (i.e., 4.2 dS/m; [68]), except for the RW, blending, and cyclic water use treatments in the 90-120 cm depth ( Figure 1). However, the addition of gypsum increased ECsw in all depths for all water quality scenarios. Average annual ECsw for scenarios with annual gypsum additions of 1.7, 4.3, and 8.6 t/ha increased by 43, 55.8, and 55.6%, respectively ( Figure 1). Gypsum is a moderately soluble salt, which can increase soluble salts in the soil. As a soil amendment, it improves the S1Gw S1B S1Alt1 S1Alt6 S1Rw

Salinity Distribution in the Calcareous Soil
Average EC sw in the 0-15, 30-60, and 90-120 cm soil depths in the G0 scenario were lower than the salinity threshold of wine grapes (i.e., 4.2 dS/m; [68]), except for the RW, blending, and cyclic water use treatments in the 90-120 cm depth ( Figure 1). However, the addition of gypsum increased EC sw in all depths for all water quality scenarios. Average annual EC sw for scenarios with annual gypsum additions of 1.7, 4.3, and 8.6 t/ha increased by 43, 55.8, and 55.6%, respectively ( Figure 1). Gypsum is a moderately soluble salt, which can increase soluble salts in the soil. As a soil amendment, it improves the cation composition on the soil exchanger by replacing Na ions with Ca ions, which reduces the ESP of the soil. Overall, there was a general pattern of increasing salt storage with depth in the calcareous soil, irrespective of irrigation water quality. However, the extent of salinity build-up depends on the ionic composition of the irrigation water.
Similarly, annual EC sw values were also below the salinity threshold of sensitive crops (4 dS/m = EC e of 2 dS/m), such as vegetables and wine grapes in the 0-15 cm soil depth under all types of irrigation, including with recycled water (RW), blending (B), and cyclic use (Alt1 and Alt6) with an 8.6 t/ha gypsum application ( Figure 2). Stevens et al. [83] estimated the EC e values to be 2.8 dS/m (or an equivalent EC sw of 5.6 dS/m) in the 20 cm depth of the NAP soil irrigated with recycled water for 10-28 years. However, their virgin soils had EC e values comparable with EC e values in scenarios irrigated with good-quality water. On the other hand, EC sw in the 30-60 cm depth showed an increasing trend for all irrigation water qualities, with average annual values ranging from 1.5 to 4.49 dS/m in different scenarios. An increase of 53-86% was comparable to an increase in the 0-15 cm depth ( Figure 1). EC sw in the RW scenario compares closely with the Stevens et al. [25] values in the 20 cm depth. Similarly, in the 90-120 cm depth, average EC sw varied from 3.11 to 7.35 dS/m, a fourfold increase compared to the values in the 0-15 cm soil depth ( Figure 1). These values are lower than the salinity tolerance threshold for moderately tolerant and tolerant crops [48]. In contrast, Stevens et al. [83] observed much higher salinities in deeper depths in virgin soils (EC e = 6.9 dS/m or approximately equivalent EC sw = 13.8 dS/m) in the NAP region. It was also observed that when wastewater was used to irrigate grapevines, inferior growth was associated with inadequate accumulation of nutrients in vine biomass [84]. Additionally, a progressive increase in soil salinity throughout several years of irrigation was observed [85,86]. However, a reduction in hydraulic conductivity, textural differences, and inherent subsoil constraints (e.g., sodicity) in deeper depths may also play a vital role in the salt accumulation in these soils. Apparently, the occurrence of such soil characteristics may reduce the salts' leaching, thus favoring storage in deeper depths in the calcareous soil.

Dynamics of the Sodium Adsorption Ratio (SAR) in the Calcareous Soil
Annual gypsum additions at a rate of 1.7, 4.3, and 8.6 t/ha showed a remarkable reduction in average SAR values ( Figure 1). Initial average SAR values (4.6-13.1) decreased by almost 50% with annual applications of 1.7 t/ha of gypsum (G1). These values were further reduced by 1-25% with increased annual gypsum additions of 4.3 (G2) and 8.6 t/ha (G3). However, this impact was observed only in the surface soil layer (0-15 cm). The resulting SAR values fell below the critical threshold (SAR < 4) for calcareous soils [12] under irrigation with all water qualities, except in the case of continuous use of recycled water. For example, in the RW scenario, SAR (10-13) values decreased to 4.6-8.4 with the annual gypsum application of 8.6 t/ha. An almost similar reduction in SAR was observed for other water quality scenarios. irrigate grapevines, inferior growth was associated with inadequate accumulation of nutrients in vine biomass [84]. Additionally, a progressive increase in soil salinity throughout several years of irrigation was observed [85,86]. However, a reduction in hydraulic conductivity, textural differences, and inherent subsoil constraints (e.g., sodicity) in deeper depths may also play a vital role in the salt accumulation in these soils. Apparently, the occurrence of such soil characteristics may reduce the salts' leaching, thus favoring storage in deeper depths in the calcareous soil. Figure 2. Impact of different water qualities (recycled water, RW; good water, GW; blending of RW and GW in a 1:1 ratio, B; RW and GW used in monthly cycles, Alt1; and RW and GW used in half-yearly cycles, Alt6) and annual gypsum applications of 8.6 t/ha (G3) on ECsw for calcareous soil (S1) under wine grape cultivation.

S1RwG3
S1GwG3 S1BG3 S1Alt1G3 S1Alt6G3 a) 0-15 cm b) 30-60 cm c) 90-120 cm Years Figure 2. Impact of different water qualities (recycled water, RW; good water, GW; blending of RW and GW in a 1:1 ratio, B; RW and GW used in monthly cycles, Alt1; and RW and GW used in half-yearly cycles, Alt6) and annual gypsum applications of 8.6 t/ha (G3) on EC sw for calcareous soil (S1) under wine grape cultivation.
Average annual SARs in different depths obtained with the annual gypsum application of 8.6 t/ha are shown in Figure 3. SAR in the surface soil layer (0-15 cm) was effectively maintained below the calcareous soil threshold (SAR = 4) under all irrigation treatments, except for RW during some years, although these values were close to the critical threshold as well. In the 30-60 cm soil depth, average annual SAR values decreased significantly compared to G0, but remained above the critical value in all water quality scenarios except for GW. Average annual SAR remained constant during the initial 2-3 years, then increased before stabilizing after 10-12 years of irrigation at values larger than 6 for RW and between 4 and 6 for blending and a cyclic use of RW and GW. Similarly, in the 90-120 cm depth, SAR was maintained near the critical value for 5-6 years of irrigation, then increased and attained a constant value after 20 years of irrigation. Subsequent small changes in the SAR at different soil depths were induced by the rainfall variability rather than the water quality as the soil profile attains a quasi-steady state. The final value was higher than 10 for RW and between 7 and 9 for blending and cyclic treatments. A gradual increase in SAR in deeper depths of the calcareous soil seems to be influenced by RW irrigation. Therefore, annual gypsum applications of 8.6 t/ha effectively reduced SAR below the threshold in the surface soil (0-15 cm) but could not fully reclaim the soil in deeper depths.
Average annual SAR values for RW irrigation were higher than the threshold (SAR = 4) for calcareous soils estimated in the previous [10,12] and current studies. However, it is appropriate to note that reclaimed SAR values are below the US Salinity Lab classification [87]. It is expected that marginally higher SAR values may not dramatically impact the overall crop performance, even though they may lead to high ESP in the soils, which is visible in average ESP values (Figure 1). Apart from RW, blending (B) and cyclic (Alt1 and Alt6) irrigation scenarios also showed ESP values corresponding to the sodic class (ESP 6-15%), as per the Australian sodic soil categorization [58]. However, these values are well below the US standards (ESP > 15%; [87]) defining sodic soils. changes in the SAR at different soil depths were induced by the rainfall variability rather than the water quality as the soil profile attains a quasi-steady state. The final value was higher than 10 for RW and between 7 and 9 for blending and cyclic treatments. A gradual increase in SAR in deeper depths of the calcareous soil seems to be influenced by RW irrigation. Therefore, annual gypsum applications of 8.6 t/ha effectively reduced SAR below the threshold in the surface soil (0-15 cm) but could not fully reclaim the soil in deeper depths. Figure 3. Impact of different water qualities (recycled water, RW; good water, GW; blending RW and GW in a 1:1 ratio, B; RW and GW used in monthly cycles, Alt1; and RW and GW used in half-yearly cycles, Alt6) and annual gypsum applications of 8.6 t/ha (G3) on annual SAR for calcareous soil (S1) under wine grape cultivation.

S1RwG3
S1GwG3 S1BG3 S1Alt1G3 S1Alt6G3 a) 0-15 cm b) 30-60 cm c) 90-120 cm Years Figure 3. Impact of different water qualities (recycled water, RW; good water, GW; blending RW and GW in a 1:1 ratio, B; RW and GW used in monthly cycles, Alt1; and RW and GW used in half-yearly cycles, Alt6) and annual gypsum applications of 8.6 t/ha (G3) on annual SAR for calcareous soil (S1) under wine grape cultivation.

Exchangeable Sodium Percentage (ESP) in the Calcareous Soil
Average annual ESP for different treatments without gypsum ranged from 14.1-21.9, 24.6-40, and 35.7-40.6% in the 0-15, 30-60, and 90-120 cm soil depths, respectively ( Figure 1). However, annual gypsum applications of 1.7 t/ha reduced ESP to between one-third and one-half of the corresponding values without gypsum applications (Figure 1). Further additions of gypsum at 4.3 and 8.6 t/ha reduced ESP further, but the impact was much smaller and only observed in the surface depth (0-15 cm). Average ESP values in the 30-60 and 90-120 cm depths under different water use scenarios remained higher than the threshold, with most values falling into the highly sodic class (ESP > 15; [58]), which is equivalent to the normal threshold in the US standards (ESP > 15%; [87]) defining sodic soils. In fact, irrigation with good-quality water also showed ESP values above the threshold (ESP > 6%) in the 90-120 cm depth. The combination of other water qualities and gypsum applications had only a small impact on reducing SAR and ESP in the soils in deeper depths. Application and mixing of gypsum in the upper soil layer (0-20 cm) had a less positive impact in deeper depths.
The impact of gypsum applications of 8.6 t/ha on annual ESP values in different depths is less prominent as the ESP values remained higher than the threshold (Figure 4). It is interesting to note that ESP values under GW irrigation remained almost constant over the entire simulation period (2018-2050). ESP values lower than critical values were observed only in the surface depths. Annual ESP in the 30-60 and 90-120 cm depths showed a similar trend as observed for SAR. Especially in deeper depths, ESP remained constant for 6-7 years, then gradually increased and stabilized at a value larger than 27 and 33 under the RW treatment after 10 (30-60 cm) and 20 years (90-120 cm) years of irrigation, respectively ( Figure 4). Subsequent small changes in the ESP were associated with the rainfall variability. This further suggests that observed constant ESP values probably represent inherent sodicity that typically occurs in the NAP soils in deeper depths and is more prevalent in calcareous soils. Assouline and Narkis [24] also found that high SAR and high salinity of treated effluent irrigation water are often accompanied by high ESP values that increase with depth. Nonetheless, high sodicity in the 90-120 cm soil depth may not impact plant growth, including wine grapes, as more than 90% of the roots of agricultural crops remain within the first-meter soil depth [70]. Likewise, 50% of the roots of cereals, pulses, and oilseed crops are accumulated within the 8-20 cm soil depth [88]. Nevertheless, such high ESP conditions can lead to restricted drainage or perched water table conditions, usually prevalent in the NAP region. Many studies have shown that ESP generally increases linearly with SAR of irrigation water [89][90][91].
Soil Syst. 2021, 5, 49 15 of 32 Figure 4. Impact of different water qualities (recycled water, RW; good water, GW; blending RW and GW in a 1:1 ratio, B; RW and GW used in monthly cycles, Alt1; and RW and GW used in half-yearly cycles, Alt6) and annual gypsum applications of 8.6 t/ha (G3) on annual ESP for calcareous soil (S1) under wine grape cultivation.

Role of Leaching on Salinity Control
It is evident from the discussion in the previous section that annual gypsum applications at a rate of 8.6 t/ha and blending or cyclic use irrigation options are unable to lower ECsw, SAR, and ESP of calcareous soils as per Australian standards, especially in deeper depths. Therefore, simulations were conducted to evaluate the impact of different leaching fractions (0.2, 0.3, 0.4, and 0.5) on reducing the harmful impact of RW irrigation. It can be seen that increasing the leaching fraction (LF) from 0.2 to 0.5 has a moderate effect on reclamation of the 0-15 cm soil depth because ECsw, SAR, and ESP were reduced only by 6.3, 5.1, and 9.9% ( Figure 5). Nevertheless, ECsw, SAR, and ESP were maintained near or below the respective threshold levels with RW irrigation. Blending and cyclic treatments had much lower values compared to RW irrigation (results not shown).

S1RwG3
S1GwG3 S1BG3 S1Alt1G3 S1Alt6G3 a) 0-15 cm b) 30-60 cm c) 90-120 cm Years Figure 4. Impact of different water qualities (recycled water, RW; good water, GW; blending RW and GW in a 1:1 ratio, B; RW and GW used in monthly cycles, Alt1; and RW and GW used in half-yearly cycles, Alt6) and annual gypsum applications of 8.6 t/ha (G3) on annual ESP for calcareous soil (S1) under wine grape cultivation.

Role of Leaching on Salinity Control
It is evident from the discussion in the previous section that annual gypsum applications at a rate of 8.6 t/ha and blending or cyclic use irrigation options are unable to lower EC sw , SAR, and ESP of calcareous soils as per Australian standards, especially in deeper depths. Therefore, simulations were conducted to evaluate the impact of different leaching fractions (0.2, 0.3, 0.4, and 0.5) on reducing the harmful impact of RW irrigation. It can be seen that increasing the leaching fraction (LF) from 0.2 to 0.5 has a moderate effect on reclamation of the 0-15 cm soil depth because EC sw , SAR, and ESP were reduced only by 6.3, 5.1, and 9.9% ( Figure 5). Nevertheless, EC sw , SAR, and ESP were maintained near or below the respective threshold levels with RW irrigation. Blending and cyclic treatments had much lower values compared to RW irrigation (results not shown).
Average EC sw reductions in the 30-60 cm soil depth by the 0.3, 0.4, and 0.5 LFs with B water irrigation were observed to be 7.3, 12.5, and 16.4%, respectively, compared to the 0.2 LF (Figure 6). Similarly, average SAR and ESP reductions by the 0.5 LF were 18.4 and 15.5%, respectively. This has reduced EC sw and SAR for B water below or near the threshold values of 4.2 dS/m and 4, respectively. However, ESP values (average ESP = 15.7) were still higher and fall into the highly sodic class [58]. This may result from the precipitation of added Ca as gypsum in this depth of calcareous soils. As expected, EC sw , SAR, and ESP values were relatively higher for RW irrigation than for blending. were still higher and fall into the highly sodic class [58]. This may result from the precipitation of added Ca as gypsum in this depth of calcareous soils. As expected, ECsw, SAR, and ESP values were relatively higher for RW irrigation than for blending. Similarly, in the 90-120 cm soil depth, ECsw, SAR, and ESP reductions for the 0.5 LF with B water irrigation were 21.2, 14.3, and 11%, respectively (Figure 7). In this scenario, average ECsw and SAR were reduced to 5.5 dS/m and 5.4, respectively. However, ESP was still maintained at an average value of 20.1%, which falls into the highly sodic class [59] for medium-and high-tolerance crop groups [49]. ECsw and SAR values are marginally higher than the wine grape threshold, but other more tolerant crops such as pistachios and barley can be grown very well under such conditions. Organic matter additions, such Similarly, in the 90-120 cm soil depth, EC sw , SAR, and ESP reductions for the 0.5 LF with B water irrigation were 21.2, 14.3, and 11%, respectively (Figure 7). In this scenario, average EC sw and SAR were reduced to 5.5 dS/m and 5.4, respectively. However, ESP was still maintained at an average value of 20.1%, which falls into the highly sodic class [59] for medium-and high-tolerance crop groups [49]. EC sw and SAR values are marginally higher than the wine grape threshold, but other more tolerant crops such as pistachios and barley can be grown very well under such conditions. Organic matter additions, such as compost, could further lower harmful impacts of high ESP, enabling the use of blending or cycling of RW and GW as sustainable management options for crop production in calcareous soils. However, continuous and indiscriminate use of recycled water, especially for growing salinity-sensitive crops, should be avoided in these soils.

Changes in Soil pH
In hard red-brown soils, the average soil pH remained higher than 8 for all irrigation treatments without gypsum additions (Figure 8). Under RW irrigation, pH values increased with depth and ranged from 8.0 to 8.4. Stevens et al. [25] also reported similar pH values and trends in the NAP soils irrigated with recycled water from the Bolivar plant. However, in all other treatments, reverse trends were observed, and average pH values varied from 8.0 to 8.3, similar to those observed for RW irrigation. Stevens et al. [25] also reported insignificant impacts of different irrigation water qualities on soil pH. However, gypsum applications were found to significantly impact soil pH [92,93], which is also visible in the present study. Annual additions of 1.7 t/ha of gypsum reduced pH below 8, irrespective of irrigation water quality. The pH reduction (1.9-3.9%) was much higher for gypsum additions of 1.72 t/ha. Further gypsum additions resulted in gradually smaller pH reductions. pH values for annual gypsum applications of 8.6 t/ha were almost similar to those for 12.9 t/ha applications. The reduction in pH varied from 5% to 6.2% as compared to the zero gypsum application.

Changes in Soil pH
In hard red-brown soils, the average soil pH remained higher than 8 for all irrigation treatments without gypsum additions (Figure 8). Under RW irrigation, pH values increased with depth and ranged from 8.0 to 8.4. Stevens et al. [25] also reported similar pH values and trends in the NAP soils irrigated with recycled water from the Bolivar plant. However, in all other treatments, reverse trends were observed, and average pH values varied from 8.0 to 8.3, similar to those observed for RW irrigation. Stevens et al. [25] also reported insignificant impacts of different irrigation water qualities on soil pH. However, gypsum applications were found to significantly impact soil pH [92,93], which is also visible in the present study. Annual additions of 1.7 t/ha of gypsum reduced pH below 8, irrespective of irrigation water quality. The pH reduction (1.9-3.9%) was much higher for gypsum additions of 1.72 t/ha. Further gypsum additions resulted in gradually smaller pH reductions. pH values for annual gypsum applications of 8.6 t/ha were almost similar to those for 12.9 t/ha applications. The reduction in pH varied from 5% to 6.2% as compared to the zero gypsum application.

Salinity Distribution in the Soil
Average ECsw values in scenarios without gypsum applications were observed to remain below the almond threshold salinity (ECsw = 3 dS/m) in the 0-15 cm soil depth (Figure 8). However, ECsw increased in deeper depths, especially for RW irrigation, for which average values reached 4.2 and 5.4 dS/m in the 30-60 and 90-120 cm depths, respectively. These values are comparable to the values for sensitive crops as reported in the Australian water quality guidelines [49]. Similarly, average ECsw in deeper depths were higher than the threshold, except for irrigation with GW water. However, annual additions of gypsum in the 0-15 cm soil depth rapidly increased soluble salt contents in this layer. It must be noted that gypsum is a salt with moderate solubility (2-2.5 g/L) and, therefore, its addition increased the total soluble salt contents in the soils. It can be seen that ECsw for the scenario

Salinity Distribution in the Soil
Average EC sw values in scenarios without gypsum applications were observed to remain below the almond threshold salinity (EC sw = 3 dS/m) in the 0-15 cm soil depth ( Figure 8). However, EC sw increased in deeper depths, especially for RW irrigation, for which average values reached 4.2 and 5.4 dS/m in the 30-60 and 90-120 cm depths, respectively. These values are comparable to the values for sensitive crops as reported in the Australian water quality guidelines [49]. Similarly, average EC sw in deeper depths were higher than the threshold, except for irrigation with GW water. However, annual additions of gypsum in the 0-15 cm soil depth rapidly increased soluble salt contents in this layer. It must be noted that gypsum is a salt with moderate solubility (2-2.5 g/L) and, therefore, its addition increased the total soluble salt contents in the soils. It can be seen that EC sw for the scenario with 1.7 t/ha of gypsum additions and GW irrigation increased three times (2.97 dS/m) in the 0-15 cm layer. Similarly, in other treatments, there were considerable increases in the soluble salt content in the upper soil layer. However, the impact of annual gypsum additions of 1.7 and 4.3 t/ha was almost negligible in deeper depths.
A further increase in annual gypsum additions to rates of 8.6 and 12.9 t/ha resulted in an increase in EC sw in deeper depths as well. This suggests that most gypsum added up to a rate of 4.3 t/ha gradually solubilized and contributed to sodic soil reclamation, even though its effects on conditions in deeper soil layers were limited. However, a subsequent increase in annual gypsum additions (8.6 and 12.9 t/ha) provided extra amounts of dissolved calcium, which gradually leached to deeper soil layers and increased the salt content there as well. Nonetheless, EC sw at the shallowest depth (0-15 cm) showed a decreasing trend with higher gypsum additions under almost all water quality scenarios. This indicates that the soil has reached its full capacity for gypsum solubilization in this layer.
Yearly fluctuations in EC sw in the 0-15, 30-60, and 90-120 cm soil depths for different water quality scenarios with annual gypsum additions of 8.6 t/ha are shown in Figure 9. Except for GW irrigation in the 0-15 and 30-60 cm depths, EC sw was higher than the almond (a sensitive crop) tolerance threshold (EC sw = 3 dS/m; [94]) under other water quality scenarios. Gypsum applications increased soluble salt contents in the soils. However, EC sw values were less than 6 dS/m in all soil depths, which is approximately equivalent to EC e of 3 dS/m. This soil salinity level is lower than the normally accepted EC e value for saline soils [87]. EC sw also increased in response to the relative composition of irrigation water, being higher for RW irrigation in all depths. Increases in soil salinity are often associated with wastewater irrigation [23,25,86] and can directly affect crop performance and many biological, chemical, and physical properties of the soil [25,95]. Soil profile normally attains a quasi-steady state in 3-5 years in response to different qualities of irrigation water and subsequent changes are associated with annual rainfall. A further increase in annual gypsum additions to rates of 8.6 and 12.9 t/ha resulted in an increase in ECsw in deeper depths as well. This suggests that most gypsum added up to a rate of 4.3 t/ha gradually solubilized and contributed to sodic soil reclamation, even though its effects on conditions in deeper soil layers were limited. However, a subsequent increase in annual gypsum additions (8.6 and 12.9 t/ha) provided extra amounts of dissolved calcium, which gradually leached to deeper soil layers and increased the salt content there as well. Nonetheless, ECsw at the shallowest depth (0-15 cm) showed a decreasing trend with higher gypsum additions under almost all water quality scenarios. This indicates that the soil has reached its full capacity for gypsum solubilization in this layer.
Yearly fluctuations in ECsw in the 0-15, 30-60, and 90-120 cm soil depths for different water quality scenarios with annual gypsum additions of 8.6 t/ha are shown in Figure 9. Except for GW irrigation in the 0-15 and 30-60 cm depths, ECsw was higher than the almond (a sensitive crop) tolerance threshold (ECsw = 3 dS/m; [94]) under other water quality scenarios. Gypsum applications increased soluble salt contents in the soils. However, ECsw values were less than 6 dS/m in all soil depths, which is approximately equivalent to ECe of 3 dS/m. This soil salinity level is lower than the normally accepted ECe value for saline soils [87]. ECsw also increased in response to the relative composition of irrigation water, being higher for RW irrigation in all depths. Increases in soil salinity are often associated with wastewater irrigation [23,25,86] and can directly affect crop performance and many biological, chemical, and physical properties of the soil [25,95]. Soil profile normally attains a quasi-steady state in 3-5 years in response to different qualities of irrigation water and subsequent changes are associated with annual rainfall. Figure 9. Impact of different water qualities (recycled water, RW; good water, GW; blending of RW and GW in a 1:1 ratio, B; RW and GW used in monthly cycles, Alt1; and RW and GW used half-yearly cycles, Alt6) and annual gypsum applications of 8.6 t/ha (G3) on ECsw for hard red-brown soil (S2) under almond cultivation.  Figure 9. Impact of different water qualities (recycled water, RW; good water, GW; blending of RW and GW in a 1:1 ratio, B; RW and GW used in monthly cycles, Alt1; and RW and GW used half-yearly cycles, Alt6) and annual gypsum applications of 8.6 t/ha (G3) on EC sw for hard red-brown soil (S2) under almond cultivation.

Dynamics of the Sodium Adsorption Ratio (SAR) in the HRB Soil
Average SAR in the hard red-brown soil in the 0-15 cm soil depth under RW irrigation without gypsum additions was almost three times higher than under GW irrigation ( Figure 8). However, SAR increased under GW irrigation as well, by 44% in the 30-60 cm depth and threefold in the 90-120 cm soil depth. The corresponding SAR values for these depths under RW irrigation were 12.9 and 14.8, respectively. Similarly, in other irrigation scenarios, including blending (B) and cycled use (Alt1 and Alt6) of GW and RW, SAR values were in between the two extremes discussed above and showed patterns with depth similar to RW and GW. However, SAR values for all water quality scenarios were higher than the threshold (SAR = 3.5) for the hard red-brown soil [12]. On the other hand, gypsum additions produced a considerable reduction in SAR values in all depths. A gypsum application rate of 1.7 t/ha (G1) reduced SAR by 14.7-33.5%, 18-30.1%, and 23.2-58.1% in the 0-15, 30-60, and 90-120 cm depths in different scenarios, respectively ( Figure 8). Overall, there was a SAR reduction of 28.3, 40.7, 52.0, and 54.3% for annual gypsum applications of 1.7, 4.3, 8.6, and 12.9 t/ha, respectively, irrespective of the irrigation water. This shows that the reduction in SAR stabilized at annual gypsum applications of 8.6 t/ha, and higher annual applications may not be economical. In addition, annual gypsum applications of 12.9 t/ha produced only a moderate reduction (0.2-5.8%) in SAR in all the scenarios and in all soil depths.
Annual SAR values remained higher than the threshold (3.5) for the hard red-brown soil in all depths for annual gypsum applications of 8.6 t/ha in all water quality scenarios, except for GW irrigation ( Figure 10). For RW irrigation, SAR was reduced to 6 in the 30-60 cm depth, while it attained an annual value of 4.5 in blending and alternate use of GW and RW scenarios. Similarly, annual SAR was reduced to 9, 6, and 6 under RW, B, and Alt treatments in the 90-120 cm depth, respectively. Constant SAR values (quasisteady state) in the 30-60 and 90-120 cm depths were attained after 6 and 10 years of irrigation, respectively ( Figure 10). However, these values can still pose a serious threat to developing high ESP in the hard red-brown soil. Conversely, these values are comparable to the values for sensitive crops reported in the Australian water quality guidelines [49]. However, Suarez et al. [95] demonstrated that an increase in irrigation water SAR from 2 to 4 significantly decreased infiltration rates in clay soils. Lower infiltration rates would result in higher surface runoff during rainfall events, leading to less water availability for crops and ultimately lower crop yields. Almond is extremely sensitive to SAR in the range of 2-8 [49] under non-saline conditions. These conditions can lead to leaf tip burn, and leaf scorch can have an immense impact on the tree growth and kernel yield. However, increased levels of soluble salts can Almond is extremely sensitive to SAR in the range of 2-8 [49] under non-saline conditions. These conditions can lead to leaf tip burn, and leaf scorch can have an immense impact on the tree growth and kernel yield. However, increased levels of soluble salts can help in reducing the harmful impact of high SAR. SAR of 6-10 and EC sw of 1.87 dS/m (RW) can lead to moderate to high sodicity hazards, depending on soil texture and rainfall trends. Heavy textured soils, such as the hard red-brown soil, can develop high sodicity hazards and create considerable management issues for sustainable production.

Exchangeable Sodium Percentage (ESP) in the HRB Soil
Similarly, average ESP values in the 0-15 cm soil depth for scenarios without gypsum additions ranged from 17.5% to 36.1%, being higher under RW irrigation and lower under GW irrigation (Figure 8). Average ESP values were within this range for blending and alternate use of RW and GW scenarios. Overall, in the 30-60 and 90-120 cm depths, ESP values increased by 24 and 70.5%, respectively, as compared to the surface layer. Usually, ESP in the 90-120 cm depth was more than 45% higher under all irrigation scenarios. Hulugalle et al. [30] found that wastewater irrigation increased exchangeable Na in all depths. This indicates significant sodicity and associated inherent subsoil hydraulic constraints in deeper depths in the hard red-brown soil. However, average ESP values were reduced by 27, 37, 45, and 47% with annual gypsum additions at a rate of 1.7, 4.3, 8.6, and 12.9 t/ha, respectively, irrespective of water quality and soil depth, indicating a gradual decrease in the ESP reduction with gypsum dose (Figure 8). Nevertheless, ESP values indicated highly sodic conditions (ESP > 15%), especially in the middle (30-60 cm) and deeper (90-120 cm) depths. The exception is under GW irrigation, where the maximum effectiveness of gypsum applications in reducing ESP occurred in all depths. This suggests that along with gypsum, water quality plays a critical role in reducing sodic conditions.
Annual ESP values, though reduced by annual gypsum applications (8.6 t/ha), are still higher than the critical threshold value of 6 for Australian soils in all irrigation scenarios and in all soil depths ( Figure 11). Constant ESP values of 9, 9, and 14 occur in the 0-15, 30-60, and 90-120 cm depths, respectively, under GW irrigation. However, corresponding ESP values doubled and tripled, respectively, for scenarios with blending (B) and cyclic use (Alt1 and Alt6) and RW in all depths compared to GW irrigation. Therefore, only some shallowrooted crops can be grown with appropriate annual gypsum applications (8.6 t/ha soil) and good-quality water. All other water qualities are likely to create salinity and sodicity problems at the soil surface, as well as in deeper soil layers. Soil profile normally attains a quasi-steady state ESP in 5-10 years in response to different qualities of irrigation water and subsequent changes are associated with the variability in the annual rainfall. In fact, RW irrigation could induce severe salinity and sodicity conditions for almonds and other sensitive crops. These issues are not fully addressed by annual gypsum applications for ameliorating sodic soil conditions. However, annual gypsum applications of 8.6 t/ha reduced the soil ESP by more than half. High ESP leads to a decline of soil physical properties, in the form of clay dispersion, with the subsequent collapse of soil structure, pore blockage, and eventually a decrease in soil permeability, leading to waterlogging, poor plant performance, decreased leaching, and salinization. Brady and Weil [96] reported that for successful irrigation with reclaimed water, the saturated hydraulic conductivity (K s ) should ideally be in the range of 1-15 cm/h. Intensive irrigation with treated wastewater in loam and clay soils was shown to result in a significant increase in clay dispersion and eluviation from the upper soil layers [37]. These results corroborate the usefulness of UNSATCHEM for simulating salinity and sodicity hazards of long-term use of RW water for irrigation and predicting the amelioration by gypsum use as shown in other studies [81,97].
draulic conductivity (Ks) should ideally be in the range of 1-15 cm/h. Intensive irrigation with treated wastewater in loam and clay soils was shown to result in a significant increase in clay dispersion and eluviation from the upper soil layers [37]. These results corroborate the usefulness of UNSATCHEM for simulating salinity and sodicity hazards of long-term use of RW water for irrigation and predicting the amelioration by gypsum use as shown in other studies [81,97]. Figure 11. Impact of different water qualities (recycled water, RW; good water, GW; blending of RW and GW in a 1:1 ratio, B; RW and GW used in monthly cycles, Alt1; and RW and GW used in half-yearly cycles, Alt6) and annual gypsum applications of 8.6 t/ha (G3) on annual ESP for hard red-brown soil (S2) under almond cultivation.

Role of Leaching for Salinity Control
Long-term modelling simulations were conducted to evaluate the impact of different leaching fractions (0.2, 0.3, 0.4, and 0.5) on reducing salinity and sodicity hazards in the HRB soil ( Figure 12). It can be observed that relative average reductions in ECsw, SAR, and ESP under RW irrigation were significantly higher in all depths compared to calcareous soils because the amount of water used for irrigation of almonds grown in the HRB soil was three times higher compared to irrigation of wine grapes in the Cal soil. In the 0-15  Figure 11. Impact of different water qualities (recycled water, RW; good water, GW; blending of RW and GW in a 1:1 ratio, B; RW and GW used in monthly cycles, Alt1; and RW and GW used in half-yearly cycles, Alt6) and annual gypsum applications of 8.6 t/ha (G3) on annual ESP for hard red-brown soil (S2) under almond cultivation.

Role of Leaching for Salinity Control
Long-term modelling simulations were conducted to evaluate the impact of different leaching fractions (0.2, 0.3, 0.4, and 0.5) on reducing salinity and sodicity hazards in the HRB soil ( Figure 12). It can be observed that relative average reductions in EC sw , SAR, and ESP under RW irrigation were significantly higher in all depths compared to calcareous soils because the amount of water used for irrigation of almonds grown in the HRB soil was three times higher compared to irrigation of wine grapes in the Cal soil. In the 0-15 cm soil depth, average reductions of EC sw , SAR, and ESP at the 0.5 LF of RW irrigation were 16.9, 19.3, and 12.1%, respectively ( Figure 12). Therefore, average EC sw , SAR, and ESP attained values of 4.7 dS/m, 5.5, and 22.8%, respectively, under RW irrigation. This extent of salinity and SAR may not create any problems for tolerant crops but can reduce the yield for sensitive crops.
On the other hand, the blending option with the 0.5 LF reduced EC sw to 3.55 dS/m and SAR to 3.3 in the 0-15 cm depth, which are comparable values to the almond thresholds ( Figure 13). ESP was reduced to an average value of 16.3%, which is higher than the lower limit for highly sodic (ESP > 15%) conditions. Meanwhile, in the 30-60 cm depth, the blending (B) option with the 0.5 LF reduced the average EC sw and SAR to 1.96 dS/m and 3.4, respectively, showing a positive impact of leaching ( Figure 14). It is interesting to note that these values are well below the sensitive crop thresholds, including almonds. However, the average ESP was maintained at 15.7%, slightly lower than the corresponding value at the surface depth. Reductions in EC sw , SAR, and ESP in the 30-60 cm soil depth ranged from 3.4-60.4, 15.5-30.4, and 12.7-26.9%, respectively, with corresponding mean values of 49.0, 22.0, and 17.9%, respectively, compared to the 0.2 LF (Figure 14). cm soil depth, average reductions of ECsw, SAR, and ESP at the 0.5 LF of RW irrigation were 16.9, 19.3, and 12.1%, respectively ( Figure 12). Therefore, average ECsw, SAR, and ESP attained values of 4.7 dS/m, 5.5, and 22.8%, respectively, under RW irrigation. This extent of salinity and SAR may not create any problems for tolerant crops but can reduce the yield for sensitive crops. On the other hand, the blending option with the 0.5 LF reduced ECsw to 3.55 dS/m and SAR to 3.3 in the 0-15 cm depth, which are comparable values to the almond thresholds ( Figure 13). ESP was reduced to an average value of 16.3%, which is higher than the lower limit for highly sodic (ESP >15%) conditions. Meanwhile, in the 30-60 cm depth, the blending (B) option with the 0.5 LF reduced the average ECsw and SAR to 1.96 dS/m and 3.4, respectively, showing a positive impact of leaching ( Figure 14). It is interesting to note that these values are well below the sensitive crop thresholds, including almonds. However, the average ESP was maintained at 15.7%, slightly lower than the corresponding value at the surface depth. Reductions in ECsw, SAR, and ESP in the 30-60 cm soil depth ranged from 3.     Similarly, average values of EC sw , SAR, and ESP in the 90-120 cm depth for the 0.5 LF with blended water were observed to be reduced by 33.1, 34.8, and 28.9%, respectively, compared to the corresponding values for the 0.2 LF (Figure 15). Reading et al. [98] also reported that long-term leaching involving significantly higher pore volumes than is commonly reported in the literature is required to obtain a stable Ks in sodic clay soils. Average values of EC sw , SAR, and ESP were similar to the respective values obtained in the 30-60 cm soil depth, and ESP was still higher than the Australian classification of sodic soils. However, values are close to the ESP threshold in the US Salinity Lab classification [87]. Overall, blending RW and GW with annual gypsum applications of 8.6 t/ha coupled with the 0.5 LF reduced the EC sw and SAR in the HRB soil to the thresholds for sensitive crops. These results are similar to Shaygan et al. [78]. They found UNSATCHEM an effective tool to simulate the salt leaching in amended sodic soils. However, ESP was maintained at a relatively higher level, creating restricted drainage/perched water table conditions if not managed appropriately. However, the harmful impacts of recycled water irrigation can be further managed with the addition of organic manures/composts, which is a common practice among growers in the NAP region.
On the other hand, a further increase in the LF (>0.5) could increase the irrigation water amount to a level which would not be economically viable and which would insert immense amounts of salts into the soil, which would leach down to deeper soil layers. Ultimately, these salts can pose a serious threat to underground aquifers. Thus, immense trade-offs exist between the extent of irrigation, salts leached from the root zone, and increased risks of groundwater pollution. Apart from adopting annual gypsum applications and appropriate leaching fractions, organic matter/compost additions, and appropriate nutrient management would further help reduce the harmful impacts of constrained use of RW for sustainable crop production in the NAP region.
These results are similar to Shaygan et al. [78]. They found UNSATCHEM an effective tool to simulate the salt leaching in amended sodic soils. However, ESP was maintained at a relatively higher level, creating restricted drainage/perched water table conditions if not managed appropriately. However, the harmful impacts of recycled water irrigation can be further managed with the addition of organic manures/composts, which is a common practice among growers in the NAP region. On the other hand, a further increase in the LF (>0.5) could increase the irrigation water amount to a level which would not be economically viable and which would insert immense amounts of salts into the soil, which would leach down to deeper soil layers. Ultimately, these salts can pose a serious threat to underground aquifers. Thus, immense trade-offs exist between the extent of irrigation, salts leached from the root zone, and increased risks of groundwater pollution. Apart from adopting annual gypsum applications and appropriate leaching fractions, organic matter/compost additions, and appropriate nutrient management would further help reduce the harmful impacts of constrained use of RW for sustainable crop production in the NAP region.
While it is important for any modelling study to validate the results with field observations to increase confidence in the modelling outcome, conducting long-term multifactor field studies involving varied water qualities, soil types, crops, and management op- While it is important for any modelling study to validate the results with field observations to increase confidence in the modelling outcome, conducting long-term multifactor field studies involving varied water qualities, soil types, crops, and management options is very complex and expensive to undertake. Furthermore, the verification of predictions made into the future is not normally possible, except by commencing a long-term field trial and waiting for several years, after which a so-called retrospective study (hindcasting) can be done. If validation data are not available, model parameters and boundary conditions should be constrained as much as possible using site-specific data. This approach was adopted in the current study.

Conclusions
This modelling study evaluated the impact of long-term (2018-2050) management options involving the use of either only recycled water (RW) or RW in various combinations with good-quality water (GW), such as blending and alternate use. Management options included different leaching fractions (0.2, 0.3, 0.4, and 0.5) and the annual application of varied doses of gypsum (0, 1.72, 4.3, 8.6, and 12.9 t/ha) to reduce salinity (EC sw ) and sodicity hazards (SAR and ESP).
Results showed that the degree of salinity and sodicity problems increased manifolds in the calcareous and hard red-brown soils with the continuous use of RW irrigation. Therefore, continuous soil amelioration should be aligned with the use of RW for irrigation. Results further revealed that applying 50% higher irrigation than leaching (LF 0.5) reduced the salinity levels below crop thresholds. Incorporating annual gypsum applications of 8.6 t/ha in calcareous and hard red-brown soils decreased SAR below the threshold (<3), especially in the blending scenario and in other similar combinations (Alt1 and Alt6) of RW and GW, but not under RW alone. Combined applications of gypsum (8.6 t/ha) and LF of 0.5 reduced the ESP fourfold, which demonstrated the important role of the long-term use of such ameliorants. However, reliance on a regular soil amendment with gypsum should be accompanied by a better understanding of soil water interactions, clay mineralogy, and the fate of salts, as well as the effectiveness of drainage in the soils.
Author Contributions: Conceptualization, methodology, resources, supervision, project administration, and funding acquisition: J.W.C., D.M. and V.P.; software: J.Š. and V.P.; formal analysis: V.P.; investigation: V.P., D.M. and J.W.C.; data curation: V.P.; writing-original draft preparation: V.P.; writing-review and editing: V.P., D.M., J.Š., J.W.C., P.R.P. and T.P.; visualization: V.P. All authors have read and agreed to the published version of the manuscript. is described in [10]. Values of soil exchangeable cations and Gapon constants for both soils used in the modelling simulations are shown in Table A2. Soil chemical properties (Table A2) measured in the 0-30 cm soil layer were assumed to extend to 200 cm (the domain depth) as the soils had almost uniform texture below 20 cm. Running the model during the warm-up period (48 years) allowed the chemical species to acquire a stable chemical composition in the solution and at the soil exchanger. θ r is the residual water content; θ s represents the saturated water content; α represents the inverse of an air entry value; n and L are pore distribution parameters, and K s is the soil saturated hydraulic conductivity.