Next Article in Journal
Network Analysis of Seed Flow, a Traditional Method for Conserving Tartary Buckwheat (Fagopyrum tataricum) Landraces in Liangshan, Southwest China
Next Article in Special Issue
Sustaining Yield of Winter Wheat under Alternate Irrigation Using Saline Water at Different Growth Stages: A Case Study in the North China Plain
Previous Article in Journal
Impact of Low-E Window Films on Energy Consumption and CO2 Emissions of an Existing UK Hotel Building
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Simulation of Saline Water Irrigation for Seed Maize in Arid Northwest China Based on SWAP Model

1
College of Hydraulic, Energy and Power Engineering, Yangzhou University, Yangzhou 225009, China
2
Jiangxi Water Resources Institute, Nanchang 330013, China
3
Centre for Agricultural Water Research in China, China Agricultural University, Beijing 100083, China
*
Author to whom correspondence should be addressed.
Sustainability 2019, 11(16), 4264; https://doi.org/10.3390/su11164264
Submission received: 1 July 2019 / Revised: 1 August 2019 / Accepted: 5 August 2019 / Published: 7 August 2019
(This article belongs to the Special Issue Water Resources and Green Growth)

Abstract

:
Water resource shortages restrict the economic and societal development of China’s arid northwest. Drawing on groundwater resources for irrigation, field experiments growing seed maize (Zea mays L.) were conducted in 2013 and 2014 in the region’s Shiyang River Basin. The Soil–Water–Atmosphere–Plant (SWAP) model simulated soil water content, salinity, and water–salt transport, along with seed maize yield, in close agreement with measured values after calibration and validation. The model could accordingly serve to simulate different saline water irrigation scenarios for maize production in the study area. Waters with a salinity exceeding 6.0 mg/cm3 were not suitable for irrigation, whereas those between 3.0 and 5.0 mg/cm3 could be acceptable over a short period of time. Brackish water (0.71–2.0 mg/cm3) could be used with few restrictions. Long-term (five years) simulation of irrigation with saline water (3.0–5.0 mg/cm3) showed soil salinity to increase by over 9.5 mg/cm3 compared to initial levels, while seed maize yield declined by 25.0% compared with irrigation with brackish water (0.71 mg/cm3). An irrigation water salinity of 3.0–5.0 mg/cm3 was, therefore, not suitable for long-term irrigation in the study area. This study addressed significance issues related to saline water irrigation and serves as a guide for future agricultural production practices.

1. Introduction

The shortage of water resources is a prominent problem restraining social and economic development in China’s arid northwest. Accounting for 24.5% of the nation’s land area, northwest China’s arid regions span about 2.5 × 106 km2. The total quantity of water resources in China’s arid northwest is estimated at 197.9 billion m3, only 5.84% of the nation’s overall supply [1]. In the region, the mean quantities of water resources per capita and per hectare are 68% and 27%, respectively, of the national average [2]. The region’s scarcity of precipitation and strong evaporative loss result in the extreme scarcity of surface water resources; however, shallow saline groundwater resources offer large reserves, with a salinity of 2–5 g/L amounting to 88.6 billion m3 [3]. The irrigation area is about 3 million hectares in the arid area of northwest China [4]. In arid areas where freshwater resources are scarce, the exploitation of underground water resources and the use of saline water irrigation technologies in agricultural production are among the important measures implemented to solve water shortage issues. With a century-old history in foreign countries, saline water irrigation was implemented successfully on a number of occasions. Several studies undertaken on the north China plain and in China’s northwest and coastal regions reported on good results with saline water irrigation practices. Studies indicated that saline water irrigation can alleviate soil drought, and the yield of some salt-tolerant crops under saline water irrigation is similar to that under freshwater irrigation, increasing production compared with no irrigation [5,6]. However, should saline water be used improperly, it can lead to the accumulation of soil salinity, as well as changes in soil physical and chemical properties, causing soil secondary salinization and threatening crop growth [7]. How to use saline water resources to irrigate scientifically, rationally, safely, and effectively was always a key concern to researchers. Most studies on saline water irrigation employed indoor or fixed field experiments to analyze the relatively short-term influence of saline water irrigation on the soil environment and crop growth [8,9,10,11]. However, the effects of saline water irrigation represent a long-term process. When employing conventional long-term field experiments, these can be affected by many external factors, in addition to placing high demands on manpower, finances and experimental facilities. Accordingly, simulating the effects of long-term saline water irrigation on the soil environment and crop growth using mathematical models was proven an effective and economical research method. The Soil–Water–Atmosphere–Plant (SWAP) model, based on the Richards and advection dispersion equations, being a comprehensive model for water, heat, and solute transport [12], can be used to simulate the water, heat, and solute fluxes in the SPAC (Soil-Plant-Atmosphere Continuum), as well as the growth process of crops. The SWAP model was widely applied to simulate soil water–salt transport and its short and long-term impact on crops and the environment around the world [13,14,15,16]. Yang et al. [17] simulated the effects of brackish water irrigation under different irrigation quotas on soil salt accumulation and spring wheat yield in China’s Hetao irrigation area. They also predicted soil root layer salt distribution and balance under long-term brackish water irrigation. Using SWAP, Kumar et al. [18] simulated crop root zone soil salt dynamics and relative wheat yield under different saline water irrigation conditions for the region of New Dehli, India. Studying spring wheat grown under a deficit of irrigation with saline water in China’s Shiyang River basin, Jiang et al. [19] used SWAP to simulate and predict soil water–salt transport under different conditions. The salinization process led to the soils reaching equilibrium after a few years’ use of saline water for irrigation. The SWAP model could be used to simulate the soil salt dynamics and wheat yield with acceptable accuracy under irrigated saline environment.
Northwest China’s unique geographical and climatic conditions, along with its abundant solar and thermal resources, led the region’s seed maize (hereafter simply maize) planting area to rise sharply and the crop to become economically important. However, few studies used the SWAP model to simulate and predict the effects of long-term saline water irrigation on soil water–salt transport and maize growth in China’s arid northwest.
The objectives of this study were (1) to calibrate and validate the applicability of SWAP model by comparing simulated results with measured data drawn from field experiments; (2) to simulate and predict the long-term effects of saline water irrigation with waters of different salinity on soil water–salt transport and maize yield.

2. Materials and Methods

2.1. Field Experiments

The field experiments were conducted in 2013 and 2014 at the Shiyang River Experimental Station of China Agricultural University (37°52′ north (N) latitude, 102°52′ east (E) longitude, elevation 1581 m). The station is located on the upper–middle reaches of the Shiyang River, near Wuwei city, Gansu Province, China (Figure 1). The region has a typical continental temperate dry climate, with an annual average (1951–2013) precipitation of 164.5 mm and an annual average evaporation of 2000 mm [20]. The mean underground water table exceeds 40 m in depth below the ground surface.
Nine non-weighing lysimeters (Figure 2) were used for the experiments, each with an area of 6.66 m2 (3.33 m × 2 m) and a depth of 3 m. Adjacent lysimeters were separated by concrete to avoid lateral seepage and surface runoff transfer. The bottom of each lysimeter was a cement floor. The physical and chemical properties of the soil prior to the onset of the experiments are presented in Table 1.
Saline water irrigation was performed with water bearing 0.71 mg/cm3 (freshwater, s0), 3.0 mg/cm3 (brackish water, s3), or 6.0 mg/cm3 (saline water, s6). These three treatments were replicated three times, with nine test plots laid out in a split-plot arrangement (Figure 3). Freshwater bearing 0.71 mg/cm3 was pumped directly from the experimental station well. Saline water bearing 3.0 mg/cm3 or 6.0 mg/cm3 was prepared artificially by dissolving NaCl, MgSO4, and CaSO4 in freshwater in mass ratio of 2:2:1. Plastic pipes equipped with water meters allowed the quantity of irrigation water to be controlled and monitored in each test pit. Maize variety “Funong 340” was sown with one line of male plants and seven lines of female plants, with 56 plants in each pit. Intra-row plant spacing was 0.25 m and inter-row spacing was 0.35 m. Sown on 20 April 2013 and 19 April 2014, maize was harvested on 13 September 2013 and 19 September 2014. The irrigation water quota refers to actual local conditions pertinent to maize production and irrigation scheduling (Table 2). Crop and field management measures followed local methods.
Using a soil auger, soil samples were taken prior to seeding and after harvest, and before and after each irrigation during the maize growth period in 2013 and 2014. Samples were obtained at depths of 0–10, 10–20, 20–40, 40–60, 60–80, and 80–100 cm. The soil moisture was determined gravimetrically by oven-drying soil at 105 °C for eight hours followed by weighing of the dry soil. Soil extractions were prepared at a soil-to-water ratio of 1:5. Electrical conductivity EC1:5 was measured using an SG-3 conductivity meter (SG3-ELK742) and translated into soil salinity through the equation, S = 0.0275EC1:5 + 0.1366 [21]. Maize plant height, leaf length (L), and leaf width (W) were measured every 15–20 days after seedling. The LAI (leaf area index) was obtained through the equation LAI = (K × L × W)/A, where K is a fitting coefficient (0.75 for maize) [22], and A is the area covered by plant leaves. The root length density was obtained through a drilling method at different maize growth stages, followed by scanning and analysis with a root software (WinRHIZO PRO 2007). The saturated hydraulic conductivity was measured using a soil permeability meter (TST-55, China) equipped with a constant head permeameter. The soil water retention curve was measured using a high-speed refrigerated centrifuge. Soil hydraulic parameters of the VG (van Genuchten) model were analyzed through the RETC (RETention Curve) program and are presented in Table 3. The maize yield of each plot was obtained by weighting method to determine yield per hectare. The meteorological data were obtained by an automatic meteorological station (Weather Hark, Campbell Scientific, Logan, UT, USA) at the experimental station. The daily rainfall, minimum temperature, and maximum temperature for the maize growth stage are shown Figure 4. Total effective precipitation during the maize growth stage in 2013 and 2014 was 64.6 mm and 141.2 mm, respectively.

2.2. SWAP Model

The SWAP model, developed by the Water Resource Group of Wageningen University, was used to simulate one-dimensional vertical soil water flow, solute transport, and crop growth [12]. The model was widely applied to simulate soil water–salt transport, crop evapotranspiration, crop growth, farmland drainage, and ground water table changing in semi-arid and arid areas around the world. The vertical soil water flow in saturated and unsaturated zone is described by Richards’ equation as follows:
θ t = C ( h ) h t = z [ K ( h ) ( h z + 1 ) ] S ( h ) ,
where θ is soil water content (cm3/cm3), t is time (days), C is differential water capacity(/cm), h is soil water pressure head (cm), z is the vertical coordinate (cm, positive uptake), k(h) is the hydraulic conductivity (cm/day), and S(h) is the soil water extraction rate by plant roots (cm3/(cm3/day)).
Solute transport is computed by the advection dispersion equation as follows:
J = q c θ ( D d i f + D d i s ) c z ,
where J is total solute flux density (g/(cm2/day)), q is vertical flow at the bottom (cm/day), c is solute concentration in soil water (g/cm3), Ddif and Ddis are the diffusion coefficient (cm2/d) and the dispersion coefficient (cm2/d), respectively, and ∂c/∂z is solute concentration gradient.
The SWAP model simulates crop growth process using the WOFOST (world food study) crop growth model, which includes a detailed crop model and a simple crop model. The simple crop growth module of Doorenbos and Kassam [23] included in SWAP was used in this study. Each growing stage, the actual yield relative to the potential yield during this growing stage is calculated by
1 Y a , k Y p , k = K y , k ( 1 T a , k T p , k ) ,
where Ya,k and Yp,k are the actual yield (kg/hm2) and the potential yield (kg/hm2), respectively, Ta,k and Tp,k are the actual transpiration (cm) and potential transpiration (cm), respectively, Ky,k is the yield response factor of growing stage k, and index k is the growing stage.
The relative yield of the full growing season is a product of the relative yield of each growing stage and is calculated as follows:
Y a Y p = k = 1 n ( Y a , k Y p , k ) ,
where Ya and Yp are the cumulative actual yield (kg/hm2) and the cumulative potential yield (kg/hm2) of the entire growing season, respectively, and n is the number of defined crop growing stages.
The input parameters required by the SWAP model are meteorological data, irrigation data, soil physical and chemical properties, soil hydraulic parameters, initial pressure head, salt data, crop data, and initial and boundary conditions. These parameters were obtained in the experiments. The measured dispersion length and molecular diffusion coefficient were 10 cm and 0.5 cm2/day, respectively. The relative uptake of solutes by roots and the solute concentration in precipitation were both zero. The soil water pressure head and the soil salinity before sowing were regarded as initial conditions. The bottom boundary condition was set as free drainage because of the large groundwater depth. A detailed introduction to the SWAP model can be found in the SWAP model theory book [12].
The RMSE (root mean square error) and MRE (mean relative error) were used to quantify the deviation of the simulated and observed data. The RMSE and MRE were calculated according to the following equations:
R M S E = 1 N i = 1 N ( P i O i ) 2 ,
M R E = 1 N i = 1 N | P i O i O i | × 100 % ,
where N is the total number of observations in the experiments, and Pi and Oi are the ith model predicted and observed values (i = 1, 2, …, N), respectively.

3. Results and Discussion

3.1. Calibration and Validation of SWAP Model

The model was calibrated and validated with field experiment data. The soil water content at depths of 10, 40, and 80 cm of s3 treatment in 2013 served to calibrate the SWAP model, while the soil water content of the same soil layers of s3 treatment in 2014 served to validate the model. The simulated soil water content agreed reasonably well with the measured values at different soil layers (Figure 5). The simulated values effectively reflected tendencies in variation of the measured data. The RMSE and MRE values were less than 0.05 cm3/cm3 and 20% in both model calibration and validation. These model accuracy statistics showed the model performance for soil water content simulation to be good. The soil hydraulic parameters after calibration are shown in Table 4.
The soil salinity at different times during the s3 treatment in 2013 served to calibrate the model, while the equivalent 2014 data served to validate the model. The simulated soil salinity agreed reasonably well with the measured data and reflected the trends in variation of the measured values at different times (Figure 6). The simulated values were mainly in especially good agreement with the measured values at the early growth stages of maize, worsening slightly at later growth stages. The RMSE and MRE values were less than 5.0 mg/cm3 and 25% in both model calibration and validation. These model accuracy statistics showed the model performance for soil salinity simulation to be good. The dispersion length and molecular diffusion coefficient after calibration were 8.5 cm and 0.5 cm2/day, respectively.
The yield of maize in 2013 and 2014 served, respectively, to calibrate and validate the SWAP model. Default crop parameters were adjusted within the range suggested by SWAP. The calibrated parameters were consistent with values in other maize studies [24]. The SWAP model’s output was the yield relative to the maximum actual yield (6303.36 kg/hm2) achieved for the s0 treatment in 2013 (set at 1.0); accordingly, the actual yield in different treatments was calculated by conversion of the simulated value. Simulated maize yield agreed reasonably well with the measured yield for the different treatments (Figure 7). The simulated maize yield mainly reflected the trends in measured yield. The RMSE and MRE values were less than 380 kg/hm2 and 10% in both model calibration and validation. Model performance for maize yield simulation was good. The minimum canopy resistance, critical level ECmax, and decline per unit ECslope of maize after calibration were 60 s/m, 1.7 dS/m, and 4.0%, respectively.
These results demonstrated that the calibrated and validated SWAP model was able to accurately simulate soil water–salt transport and maize yield, and predict these under long-term saline water irrigation with waters of different salinity in this study area.

3.2. Simulation of Saline Water Irrigation with Different Water Salinity on Soil Water–Salt Transport and Maize Yield

Because field experiments were restricted by various factors, it was difficult to comprehensively carry out saline water irrigation experiments with waters of different salinities. The calibrated SWAP model can be used to simulate saline water irrigation with different water salinity. The simulated irrigation water salinities tested were 0.71, 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, and 8.0 mg/cm3. As the Shiyang River Basin is located in an arid desert area where rainfall is scarce, there is little difference between normal flow years, and dry or wet years. During the simulation, the meteorological data were based on the average annual method (p = 50%). Drawing on 1960–2014 rainfall data, rainfall frequency analysis showed that the typical year with p = 50% was 2011, when rainfall during the maize growth period was 118.8 mm. Based on local irrigation experience and previous studies, maize is generally irrigated four or five times over its growth period. Irrigating for the fifth time is mainly to maintain maize grain plumpness, thereby obtaining a greater yield and economic income. The added irrigation occurs in the mature stage of the maize crop cycle [25]. The present simulation did not consider a fifth irrigation, but rather irrigated four times during the maize growth period, according to the crop’s water requirements. The overall irrigation water quota was 450 mm, including 120 mm at the seedling stage (5 June), 120 mm at the jointing stage (30 June), 105 mm at the heading stage (20 July), and 105 mm at the filling stage (10 August). The SWAP was simulated based on initial soil water content, soil salinity, and crop growth data of the s0 treatment in 2013. The simulated soil depth range was 0–100 cm. Under these conditions, the calibrated SWAP model was used to simulate saline water irrigation with nine different levels of irrigation water salinity. Figure 8a shows that, although the irrigation water amount was unchanged, irrigation water of different salinities altered soil water transport. The soil water content increased gradually with an increase in irrigation water salinity. At the higher irrigation water salinities, salt stress worsened, inhibiting root water absorption by the crop, and allowing more water to remain in the soil [26]. The soil water content under irrigation water salinities of 0.71, 1.0, and 2.0 mg/cm3 was basically the same, indicating that soil water content was almost the same under low irrigation water salinities. Soil salinity in the 0–100-cm soil layer under saline water irrigation increased with an increase in irrigation water salinity (Figure 8b). At the end of the maize growth period, the soil salinity under irrigation with waters of 0.71, 1.0, and 2.0 mg/cm3 salinities was basically the same; however, for waters with salinities of 3.0, 4.0, 5.0, 6.0, 7.0, and 8.0 mg/cm3, it increased by 3.73, 5.03, 6.20, 7.23, 8.18, and 9.05 mg/cm3, respectively, compared to the 0.71 mg/cm3 irrigation water. Simulated soil water flux and salt flux changed in a similar manner in the 100-cm soil layer (Figure 9); the downward soil water flux and salt flux gradually increased with an increase in irrigation water salinity, reflecting the fact that soil water content could promote soil salinity movement to a certain degree.
Soil water–salt balance under saline water irrigation with different water salinities is presented in Table 5. A negative sign for soil water change indicates that soil water content was diminished through consumption. A negative sign for bottom soil water flux indicated that soil water was leaking downward. The bottom soil salt flux was negative, indicating that soil salinity moved downward. The increase of soil salinity was negative, indicating that soil salinity was leached. Based on the principle that a lesser increase in soil salinity would result in a greater yield and water use efficiency for maize, suitable irrigation water salinity was sought for the study area. It can be seen from Table 5 that the soil salinity of the 0.71, 1.0, and 2.0 mg/cm3 irrigations increased below 50 mg/cm2, the maize yield was above 5168 kg/hm2, and the maize water use efficiency was above 1.0 kg/m3. Therefore, brackish water of 0.71–2.0 mg/cm3 can be used for irrigation in the study area. The soil salinity of 6.0, 7.0, and 8.0 mg/cm3 irrigations increased more than 180 mg/cm2. Compared with the simulated maize yield of 0.71 mg/cm3, the maize yield decreased by more than 28% and the water use efficiency was less than 0.95 kg/m3. Therefore, irrigation water salinity of 6.0, 7.0, and 8.0 mg/cm3 was not suitable for irrigation. The soil salinity of 3.0, 4.0, and 5.0 mg/cm3 irrigations increased between 80 and 160 mg/cm2. The yield reduction of maize was between 12.5% and 23.0% and the water use efficiency was between 0.98 and 0.95 kg/m3, which were slightly worse than brackish water of 0.71–2.0 mg/cm3. Therefore, saline water of 3.0–5.0 mg/cm3 can be used for irrigation over a short period of time. Whether saline water of 3.0–5.0 mg/cm3 can be used for irrigation over a long period of time still needs further study.

3.3. Long-Term Simulation for Soil Salinity Transport and Maize Yield

Whether saline water of 3.0, 4.0, and 5.0 mg/cm3 can be used for a long time in the study area was determined by allowing the SWAP model to predict the long-term effects of saline water irrigation bearing 3.0, 4.0, and 5.0 mg/cm3 on soil salt transport and maize yield. In the simulation process, the initial soil water content, initial soil salinity, irrigation scheduling, and meteorological data were unchanged. Taking the simulation results of soil water content and soil salinity at the end of each year as the initial conditions for the next year, saline water irrigation with 3.0, 4.0, and 5.0 mg/cm3 was continuously simulated for five years. The influence of irrigation water salinity on soil salinity transport in the 0–100-cm soil layer over the five-year period showed a year-to-year increase in soil salinity for all three irrigation water salinities (Figure 10a). The soil salinity increased sharply in the early stages of the simulation, and more slowly in the latter stages. At the end of the fifth year, the soil salinity of 3.0, 4.0, and 5.0 mg/cm3 irrigations reached a steady state. After five years of simulation, the soil salinity of 3.0, 4.0, and 5.0 mg/cm3 irrigations increased by 9.76, 10.95, and 11.94 mg/cm3, respectively, compared with the initial soil salinity. Simulated maize yields over the five years indicated that maize yield declined with advancing years (Figure 10b). The maize yield decreased greatly in the early stages of simulation, and more slowly thereafter. After five years of simulation, maize yield reached a steady state. After five years of simulation, the maize yields under saline water irrigation with 3.0, 4.0, and 5.0 mg/cm3 were 4097.18, 3718.98, and 3403.81 kg/hm2, decreasing by 26.14%, 32.95%, and 38.64%, respectively, compared with the simulated yields under saline water irrigation with 0.71 mg/cm3 in the first year. The simulated annual declines in maize yield for saline water irrigation at a salinity of 3.0 mg/cm3 were 12.5%, 23.86%, 25.0%, 26.14%, and 26.14% from the first to the fifth year, respectively, compared with the simulated yield for maize in the first year for irrigation water bearing a salinity of 0.71 mg/cm3. Over a long period of time, an irrigation water salinity of 3.0–5.0 mg/cm3 would gradually increase soil salinity and have a significant negative impact on maize yield. The maize yields under saline water irrigation with 3.0–5.0 mg/cm3 decreased by over 25% after five years. Therefore, saline water of 3.0–5.0 mg/cm3 was not suitable for irrigation over a long period of time.
This study’s results were consistent with those of Wu et al. [27], who found that many years’ irrigation of winter wheat with water bearing a salinity of 3.0 mg/cm3 in northern China’s Huanghuai-Hai Plain led to salt accumulation in the upper soil layer. They concluded that, if saline water with a degree of mineralization of 3.0 g/L was utilized continuously, plants would be subjected to salt stress. However, the salt in the upper soil layer could be kept in the normal range for crop growth under sufficient rainfall and leaching. After five years of saline water irrigation with 3.0 mg/cm3, the soil salinity increased slightly in the later stages and reached a steady state in this study. Crescimanno and Garofalo [28] used the SWAP model to explore different saline water irrigation scenarios, which aimed at optimizing irrigation in a Sicilian vineyard. They found that irrigation water salinities of 2.1 dS/m and 6.2 dS/m applied at a fixed annual volume of 112 mm proved to be the best strategy, reducing soil salinization and enhancing crop transpiration. This is similar to the finding that saline water of 3.0–5.0 mg/cm3 could be used for irrigation in the short term in this study. However, the difference was that saline water of 3.0–5.0 mg/cm3 was not suitable for irrigation in the long term in this study. Martínez et al. [29], simulating the water balance components in an irrigated citrus orchard in Andalucía (southern Spain) with the SWAP model, found it to be useful in supporting the development of appropriated irrigation strategies for sustainable water management in Andalucía. The model could be used as a simulation tool of saline water irrigation to simulate soil water–salt transport and maize yield in this study. Therefore, the SWAP model could be used to vet saline water irrigation strategies for sustainable water management in semi-arid and arid areas [30,31].

4. Conclusions

The RMSE and MRE values were less than 0.05 cm3/cm3 and 20%, respectively, for soil water content in the SWAP model calibration and validation. For soil salinity, the RMSE values were all lower than 5.0 mg/cm3, and the MRE values were lower than 25%. The RMSE and MRE values were less than 380 kg/hm2 and 10%, respectively, for maize yield. The simulated soil water content, soil salinity, and maize yield agreed well with the measured values. The SWAP model was able to simulate soil water–salt transport and maize yield for different saline water irrigations in the study area. The simulation results of saline water irrigation with different water salinities showed that brackish water of 0.71–2.0 mg/cm3 can be used unrestrictedly for irrigation. The yield of maize was above 5168 kg/hm2 and the water use efficiency of maize was above 1.0 kg/m3. Irrigation water salinity above 6.0 mg/cm3 was not suitable for irrigation. The yield of maize was reduced by above 28% and the water use efficiency of maize was below 0.95 kg/m3. The use of irrigation water with a salinity of 3.0–5.0 mg/cm3 could be acceptable for short-term irrigation as the yield of maize was reduced by 12.5–23.0% and the water use efficiency was 0.95–0.98 kg/m3. Long-term simulation results of saline water irrigation with 3.0–5.0 mg/cm3 indicated that soil salinity would gradually increase over the years. After five years, the soil salinity of 3.0–5.0 mg/cm3 increased by above 9.5 mg/cm3 compared with the initial soil salinity. At this irrigation water salinity, maize yield decreased over the years, dropping by over 25% after five years. An Irrigation water salinity of 3.0–5.0 mg/cm3 was deemed unsuitable for long-term irrigation in the study area. This study is important in terms of saline water irrigation and agricultural production practices in the Shiyang River Basin of northwest China.

Author Contributions

C.Y. and Q.J. collected the related data; C.Y. and Z.H. analyzed and built the SWAP model; C.Y. and S.F. analyzed the results and wrote the paper; Z.H. and S.F. modified this article.

Funding

This research was financially supported by the National Natural Science Foundation of China (51179166), Specialized Research Fund for the Doctoral Program of Higher Education of China (20123250110004), and a project funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD).

Acknowledgments

We are very grateful to Z.Q. (James H. Brace Associate Professor of McGill University, Canada) for his revision of the English language of our manuscript. We greatly appreciate to the editors and anonymous reviewers for their valuable comments.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Chen, Y.N.; Yang, Q.; Luo, Y.; Shen, Y.J.; Pan, X.L.; Li, L.H.; Li, Z.Q. Ponder on the issues of water resources in the region of northwest China. Arid Land Geogr. 2012, 35, 1–9. [Google Scholar]
  2. Tian, F.Q.; Wen, J.; Hu, H.C.; Ni, G.H. Review on water and salt transport and regulation in drip irrigated fields in arid regions. J. Hydraul. Eng. 2018, 49, 126–135. [Google Scholar]
  3. Huo, H.X.; Zhang, J.G. Advance and prospects of research on soil salt movement under saline irrigation. Water Sav. Irrig. 2015, 4, 41–45. [Google Scholar]
  4. Wang, Y.J.; Qin, D.H. Influence of climate change and human activity on water resources in arid region of northwest China: An Overview. Clim. Chang. Res. 2017, 13, 483–493. [Google Scholar] [CrossRef]
  5. Wan, S.Q.; Kang, Y.H.; Wang, D.; Liu, S.P. Effect of saline water on cucumber (Cucumis sativus L.) yield and water use under drip irrigation in north china. Agric. Water Manag. 2010, 98, 105–113. [Google Scholar] [CrossRef]
  6. Singh, A.; Panda, S.N. Effect of saline irrigation water on mustard (Brassica juncea) crop yield and soil salinity in a semi-arid area of North India. Exp. Agric. 2012, 48, 99–110. [Google Scholar] [CrossRef]
  7. Mondal, M.K.; Bhuiyan, S.I.; Franco, D.T. Soil salinity reduction and prediction of salt dynamics in the coastal rice lands of Bangladesh. Agric. Water Manag. 2001, 47, 9–23. [Google Scholar] [CrossRef]
  8. Ma, W.J.; Mao, Z.Q.; Yu, Z.R.; van Mensvoort, M.E.F.; Driessen, P.M. Effects of saline water irrigation on soil salinity and yield of winter wheat-maize in North China Plain. Irrig. Drain. Syst. 2008, 22, 3–18. [Google Scholar] [CrossRef]
  9. Wu, Z.D.; Wang, Q.J. Effect on characteristics of soil water and salt movement about saline water infiltration with different infiltration amount. Trans. Chin. Soc. Agric. Mach. 2010, 41, 67–78. [Google Scholar]
  10. Feikema, P.M.; Morris, J.D.; Connell, L.D. The water balance and water sources of a Eucalyptus plantation over shallow saline groundwater. Plant Soil 2010, 332, 429–449. [Google Scholar] [CrossRef]
  11. Li, C.; Lei, J.; Zhao, Y.; Xu, X.; Li, S. Effect of saline water irrigation on soil development and plant growth in the Taklimakan desert highway shelterbelt. Soil Tillage Res. 2015, 146, 99–107. [Google Scholar] [CrossRef]
  12. Van Dam, J.C.; Huygen, J.; Wesseling, J.G.; Feddes, R.A.; Kabat, P.; Walsum, P.E.V.; Groenendijk, P.; van Diepen, C.A. Theory of SWAP Version 2.0; Simulation of Water Flow, Solute Transport and Plant Growth in the Soil-Water-Atmosphere-Plant Environment; DLO Winand Staring Centre: Wageningen, The Netherlands, 1997; pp. 1–167. [Google Scholar]
  13. Schahbazian, N.; Eitzinger, J.; Montazar, A.; Akbari, G.; Allahdadi, I. Using simulation models to improve irrigation water application. Pak. J. Water Resour. 2007, 11, 43–57. [Google Scholar]
  14. Singh, U.K.; Ren, L.; Kang, S. Simulation of soil water in space and time using an agro-hydrological model and remote sensing techniques. Agric. Water Manag. 2010, 97, 1210–1220. [Google Scholar] [CrossRef]
  15. Verma, A.K.; Gupta, S.K.; Isaaca, R.K. Use of saline water for irrigation in monsoon climate and deep water table regions: Simulation modeling with SWAP. Agric. Water Manag. 2012, 115, 186–193. [Google Scholar] [CrossRef]
  16. Verma, A.K.; Gupta, S.K.; Isaac, R.K. Calibration and validation of SWAP to simulate conjunctive use of fresh and saline irrigation waters in semi-arid regions. Environ. Model. Assess. 2014, 19, 45–55. [Google Scholar] [CrossRef]
  17. Yang, S.Q.; Ye, Z.G.; Shi, H.B.; Lan, L.T. Simulation and prediction of rotational irrigation with salty and fresh water in the Hetao irrigation area of Inner Mongolia. Trans. Chin. Soc. Agric. Eng. 2010, 26, 8–17. [Google Scholar]
  18. Kumar, P.; Sarangi, A.; Singh, D.K.; Parihar, S.S.; Sahoo, R.N. Simulation of salt dynamics in the root zone and yield of wheat crop under irrigated saline regimes using SWAP model. Agric. Water Manag. 2015, 148, 72–83. [Google Scholar] [CrossRef]
  19. Jiang, J.; Feng, S.; Huo, Z.; Zhao, Z.; Jia, B. Application of the swap model to simulate water-salt transport under deficit irrigation with saline water. Math. Comput. Model. 2011, 54, 902–911. [Google Scholar] [CrossRef]
  20. Jiang, J.; Feng, S.; Ma, J.; Huo, Z.; Zhang, C. Irrigation management for spring maize grown on saline soil based on SWAP model. Field Crops Res. 2016, 196, 85–97. [Google Scholar] [CrossRef]
  21. Wang, Q.M.; Huo, Z.L.; Zhang, L.D.; Wang, J.H.; Zhao, Y. Impact of saline water irrigation on water use efficiency and soil salt accumulation for spring maize in arid region of China. Agric. Water Manag. 2016, 163, 125–138. [Google Scholar] [CrossRef]
  22. Li, S.E.; Kang, S.Z.; Li, F.S.; Zhang, L. Evapotranspiration and crop coefficient of spring maize with plastic mulch using eddy covariance in Northwest China. Agric. Water Manag. 2008, 97, 1214–1222. [Google Scholar] [CrossRef]
  23. Doorenbos, J.; Kassam, A.H. Yield Response to Water. FAO Irrigation and Drainage Paper 33; FAO: Rome, Italy, 1979. [Google Scholar]
  24. Ma, Y.; Feng, S.Y.; Song, X.F. Evaluation of optimal irrigation scheduling and groundwater recharge at representative sites in the North China Plain with SWAP model and field experiments. Comput. Electron. Agric. 2015, 116, 125–136. [Google Scholar] [CrossRef]
  25. Jiang, X.; Kang, S.; Tong, L.; Li, F. Modification of evapotranspiration model based on effective resistance to estimate evapotranspiration of maize for seed production in an arid region of Northwest China. J. Hydrol. 2016, 538, 194–207. [Google Scholar] [CrossRef]
  26. Ben-Asher, J.; van Dam, J.; Feddes, R.A.; Jhorar, R.K. Irrigation of grapevines with saline water: II. Mathematical simulation of vine growth and yield. Agric. Water Manag. 2006, 83, 22–29. [Google Scholar] [CrossRef]
  27. Wu, Z.D.; Wang, Q.J. Effect of saline water continuous irrigation on winter wheat yield and soil physicochemical property. Trans. Chin. Soc. Agric. Mach. 2010, 41, 36–43. [Google Scholar]
  28. Crescimanno, G.; Garofalo, P. Management of irrigation with saline water in cracking clay soils. Soil Sci. Soc. Am. J. 2006, 70, 1774–1787. [Google Scholar] [CrossRef]
  29. Martínez-Ferri, E.; Muriel-Fernández, J.L.; Rodríguez Díaz, J.A. Soil water balance modelling using SWAP: An application for irrigation water management and climate change adaptation in citrus. Outlook Agric. 2013, 42, 93–102. [Google Scholar] [CrossRef]
  30. Sarwar, A.; Bastiaanssen, W.G.M.; Feddes, R.A. Irrigation water distribution and long-term effects on crop and environment. Agric. Water Manag. 2001, 50, 125–140. [Google Scholar] [CrossRef]
  31. Tedeschi, A.; Menenti, M. Simulation studies of long-term saline water use: Model validation and evaluation of schedules. Agric. Water Manag. 2002, 54, 123–157. [Google Scholar] [CrossRef]
Figure 1. Geographical localization map of study area.
Figure 1. Geographical localization map of study area.
Sustainability 11 04264 g001
Figure 2. Non-weighing lysimeters used in the experiments.
Figure 2. Non-weighing lysimeters used in the experiments.
Sustainability 11 04264 g002
Figure 3. Arrangement diagram of test plot.
Figure 3. Arrangement diagram of test plot.
Sustainability 11 04264 g003
Figure 4. Daily meteorological data for maize growth stage in 2013 (a) and 2014 (b).
Figure 4. Daily meteorological data for maize growth stage in 2013 (a) and 2014 (b).
Sustainability 11 04264 g004
Figure 5. Comparison of the simulated and measured soil water content at different soil layers in calibration 2013 (a) and validation 2014 (b).
Figure 5. Comparison of the simulated and measured soil water content at different soil layers in calibration 2013 (a) and validation 2014 (b).
Sustainability 11 04264 g005
Figure 6. Comparison of the simulated and measured soil salinity at different times: (a) 28 June 2013; (b) 19 July 2013; (c) 31 July 2013; (d) 27 August 2013; (e) 11 June 2014; (f) 29 June 2014; (g) 24 July 2014; (h) 15 August 2014.
Figure 6. Comparison of the simulated and measured soil salinity at different times: (a) 28 June 2013; (b) 19 July 2013; (c) 31 July 2013; (d) 27 August 2013; (e) 11 June 2014; (f) 29 June 2014; (g) 24 July 2014; (h) 15 August 2014.
Sustainability 11 04264 g006
Figure 7. Comparison of the simulated and measured maize yield in calibration 2013 (a) and validation 2014 (b).
Figure 7. Comparison of the simulated and measured maize yield in calibration 2013 (a) and validation 2014 (b).
Sustainability 11 04264 g007
Figure 8. Simulation of soil water content (a) and soil salinity (b) with different saline water irrigation for the 0–100-cm soil layer.
Figure 8. Simulation of soil water content (a) and soil salinity (b) with different saline water irrigation for the 0–100-cm soil layer.
Sustainability 11 04264 g008
Figure 9. Simulation of soil water flux (a) and salt flux (b) with different saline water irrigation for the 100-cm soil layer.
Figure 9. Simulation of soil water flux (a) and salt flux (b) with different saline water irrigation for the 100-cm soil layer.
Sustainability 11 04264 g009
Figure 10. Simulation of soil salinity transport (a) and maize yield (b) over five years.
Figure 10. Simulation of soil salinity transport (a) and maize yield (b) over five years.
Sustainability 11 04264 g010
Table 1. Soil physical and chemical properties.
Table 1. Soil physical and chemical properties.
Soil Depth (cm)Sand (%)Silt (%)Clay (%)Organic Content (g/kg)Soil bulk Density (g/cm3)Field Capacity (cm3/cm3)Saturated Water Content (cm3/cm3)International Soil Texture
0–2061.0328.4310.542.601.490.270.36Sandy loam
20–4058.3330.4511.222.641.540.300.40Sandy loam
40–10053.1134.1412.755.691.550.320.42Clay loam
Table 2. Irrigation scheduling in the experiments.
Table 2. Irrigation scheduling in the experiments.
TreatmentIrrigation Water Salinity (mg/cm3)Irrigation Water Quota (mm)Overall Irrigation Water Quota (mm)
5 June30 June20 July10 August29 August
s00.71120120105105105555
s33.0120120105105105555
s66.0120120105105105555
Table 3. Soil hydraulic parameters of VG (van Genuchten) model.
Table 3. Soil hydraulic parameters of VG (van Genuchten) model.
Soil Depth (cm)Residual Water Content (cm3/cm3)Saturated Water Content (cm3/cm3)Saturated Hydraulic Conductivity (cm/d)Shape Factor α (/cm)Shape Factor nShape Factor γ
0–200.0440.3632.570.0241.4340.5
20–400.0430.3829.850.0241.4170.5
40–1000.0490.4013.710.0111.4800.5
Table 4. Soil hydraulic parameters for different soil layers after calibration.
Table 4. Soil hydraulic parameters for different soil layers after calibration.
Soil Depth (cm)Residual Water Content (cm3/cm3)Saturated Water Content (cm3/cm3)Saturated Hydraulic Conductivity (cm/day)Shape Factor α (/cm)Shape Factor nShape Factor γ
0–200.0550.3430.810.0241.4020.5
20–400.0440.3831.660.0221.4080.5
40–1000.0950.3912.000.0201.3100.5
Table 5. Simulation of soil water and salt balance for different saline water irrigation.
Table 5. Simulation of soil water and salt balance for different saline water irrigation.
Irrigation Water Salinity (mg/cm3)Water BalanceSalt BalanceYield (kg/hm2)Water Use Efficiency (kg/m3)
Irrigation (mm)Rainfall and Interception (mm)Soil Water Change (mm)Bottom Flux (mm)Evapotranspiration (mm)From Irrigation (mg/cm2)Bottom Flux (mg/cm2)Increase in Soil (mg/cm2)
0.71450108.7−107.8−118.1548.431.95−35.95−4.005546.961.01
1.0450108.7−103.9−120.3542.345.00−36.758.255483.921.01
2.0450108.7−91.0−129.9519.890.00−40.6749.335168.761.00
3.0450108.7−80.7−142.4497.0135.00−46.6188.394853.590.98
4.0450108.7−73.0−157.1474.6180.00−55.19124.814538.420.96
5.0450108.7−67.3−173.6452.4225.00−66.86158.144286.280.95
6.0450108.7−62.8−190.7430.8270.00−81.49188.513971.120.92
7.0450108.7−59.3−207.8410.2315.00−98.57216.433718.980.90
8.0450108.7−56.4−224.4390.7360.00−117.80242.203466.850.89

Share and Cite

MDPI and ACS Style

Yuan, C.; Feng, S.; Huo, Z.; Ji, Q. Simulation of Saline Water Irrigation for Seed Maize in Arid Northwest China Based on SWAP Model. Sustainability 2019, 11, 4264. https://doi.org/10.3390/su11164264

AMA Style

Yuan C, Feng S, Huo Z, Ji Q. Simulation of Saline Water Irrigation for Seed Maize in Arid Northwest China Based on SWAP Model. Sustainability. 2019; 11(16):4264. https://doi.org/10.3390/su11164264

Chicago/Turabian Style

Yuan, Chengfu, Shaoyuan Feng, Zailin Huo, and Quanyi Ji. 2019. "Simulation of Saline Water Irrigation for Seed Maize in Arid Northwest China Based on SWAP Model" Sustainability 11, no. 16: 4264. https://doi.org/10.3390/su11164264

APA Style

Yuan, C., Feng, S., Huo, Z., & Ji, Q. (2019). Simulation of Saline Water Irrigation for Seed Maize in Arid Northwest China Based on SWAP Model. Sustainability, 11(16), 4264. https://doi.org/10.3390/su11164264

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