Next Article in Journal
Review of Advancement in Variable Valve Actuation of Internal Combustion Engines
Previous Article in Journal
Assessment of Wide-Sense Stationarity of an Underwater Acoustic Channel Based on a Pseudo-Random Binary Sequence Probe Signal
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

HYDRUS-1D Simulation of Soil Water Dynamics for Sweet Corn under Tropical Rainfed Condition

1
Department of Biological and Agricultural Engineering, Universiti Putra Malaysia, Serdang 43400, Malaysia
2
Department of Agricultural Engineering, University of Agriculture Faisalabad, Faisalabad 38000, Pakistan
3
Malaysian Agricultural and Research Development Institute (MARDI), Serdang 43400, Malaysia
*
Author to whom correspondence should be addressed.
Appl. Sci. 2020, 10(4), 1219; https://doi.org/10.3390/app10041219
Submission received: 7 November 2019 / Revised: 10 January 2020 / Accepted: 13 January 2020 / Published: 11 February 2020
(This article belongs to the Section Environmental Sciences)

Abstract

:
Assessment of soil water balance is essential to understand water dynamics for optimal use of water and fertilizers. The study intended to simulate soil water dynamics in sweet corn production under tropical rainfed conditions. Surface runoff, subsurface leaching, and evapotranspiration are the main components of water balance, especially in tropical environments. Therefore, intensive field experiments and HYDRUS-1D numerical modeling were applied to investigate the water balance components and analyzing water dynamics. The study was carried out in a sweet corn field for two growing seasons under the rainfed conditions at the Malaysian Agricultural Research and Development Institute (MARDI), Serdang, Malaysia. The total water inputs during the first and second seasons were 75.8 cm and 79.7 cm, respectively. Simulated results of evapotranspiration (ET) accounted for 40.7% and 33.1% of total water input during the first and second seasons. Surface runoff accounted for 41% and 28.6% in the first and second season, respectively. Water leaching accounted for 10.6%–26.8% of total water input during both seasons respectively. As rainfall fulfilled the crop water requirement throughout the growing seasons no additional irrigation was required. The overall simulation results validate the HYDRUS-1D as an effective tool to simulate soil water dynamics under rainfed conditions.

1. Introduction

Around 190.6 million hectares produced 1076.18 million tons of corn at the rate of 5.65 tons/ha for the year 2017–2018 worldwide [1]. In 2017, Malaysia consumed around 4 million tons of grain corn. Most of the grain corn consumed in Malaysia is imported, while sweet corn is grown locally [2]. About 10,477.30 ha of the cultivated area produced 72,560.52 tons of sweet corn for the year 2017 in the country [3]. Sweet corn is used for both human and animal consumption, along with its use in industry as raw material [4]. In Malaysia, sweet corn grows on acidic, weathered soils with low pH and soil fertility, which results in low corn yield [5]. Imbalance in fertilizer application also results in low corn production in tropical regions [6]. The change in water balance components such as evapotranspiration, runoff, soil moisture, and rainfall might also affect this fertilizer imbalance.
Accurate estimation of water balance in agricultural fields is also key to water resources management in agriculture, which is the primary consumer of the water. Evapotranspiration and leaching have considered as leading sinks of water that affect soil water status in a soil–plant–atmosphere environment [7]. That is why devising water management strategies depend on information relating to crop water requirement (represented as evapotranspiration). Measuring other water balance components such as deep percolation, runoff in field conditions is also a big challenge [8]. The response of the water balance components also varies with climate geographically. The significant differences between regions necessitate a need to evaluate responses of water balance components in various geoclimatic regions. Rainfall is the most investigated components of water balance. A significant regional variation in rainfall has also been recorded [9]. Climate change also has a significant effect in altering rainfall patterns in humid tropics [10]. Therefore, in a tropical region like Malaysia, the change in rainfall patterns is more crucial than in temperature due to climate change.
Simulation models are recognized as a valuable tool to assess water transport and the extent to which management practices affect crop yield and the environment. However, the validation of simulation models for local conditions is crucial [11]. Among different available models, a software package HYDRUS-1D [12] has been widely used for simulating water flow and solute transport in soils, to analyze flow and transport processes in agricultural fields. Several researchers have applied HYDRUS-1D to evaluate water balance successfully. As compared to other models, HYDRUS-1D has the flexibility of accommodating different boundary conditions and the capability to consider the root uptake of water and nutrient simultaneously. The model is capable of simulating soil water and solute dynamics under different management practices [13,14,15,16,17,18,19,20,21,22]. All these researchers validated HYDRUS-1D as a reliable tool for such investigations.
It is necessary to improve agricultural water management and to develop new strategies to increase agricultural productivity [23]. The heavy rainfall results in the solute loss to increase manifold in the tropical climate. Therefore, the study was carried out to quantify water losses (surface runoff and leaching) from the field that results in fertilizer loss in Malaysia. Quantifying the water balance components under rainfed conditions can be useful to improve irrigation management and to assess the scale of fertilizer loss for sweet corn production. The study can further be used to link rainfall intensity with fertilizer loss. Besides, HYDRUS-1D has rarely been used for sweet corn production in tropical environment. The model validation through the field experiments in quantifying water losses is worthwhile for sweet corn production. Indeed, heavy rainfall in Malaysia can potentially meet crop water requirements. However, the duration and frequency of rainfall are uncertain that can affect the water availability to plant. Therefore, the objectives of the study were (1) to simulate the soil water dynamics under rainfed conditions, and (2) to evaluate the HYDRUS-1D model for the rainfall potential to meet the water requirement of sweet corn in the tropical region.

2. Material and Methods

2.1. Field Experiments and Measurements

2.1.1. Site Description

The experiment field located 2°59′ N latitude and 101°42′ E longitude, 43 m above the mean sea level at the Malaysian Agricultural Research and Development Institute (MARDI) in Selangor, Malaysia. Total annual precipitation of the area has been recorded as 2632 mm for the year 2017, whereas minimum and maximum temperatures were 21.3 °C and 35.8 °C, respectively. The rainfall occurs throughout the year, with seasonal variation. Most of the rainfall occurs in months October–December. February and March are relatively dry months. Research work was carried out from February to May 2018 and September to November 2018. The soil layer classified as clay throughout the sampling depth of 0–140 cm. The physical and chemical properties of soil are presented in Table 1 and Table 2. The variety of sweet corn cultivated was ‘Hibridmas’.

2.1.2. Experimental Design and Measurements

The research work was designed based on an intensive field experiment in real field conditions. The experiment was carried out in a 770 m2 (35 m × 22 m) field. Sweet corn ‘Hibridmas’ was planted during both seasons at 10 kg/ha where plant-plant (P × P) and row-row (R × R) space were maintained as 20 cm and 75 cm respectively. Harvesting dates were on 3rd May and 15th November for the first and second season 2018 respectively. The growing period for the first season (February–May 2018) was 76 days whereas for the second season (September–November 2018) it was 71 days. The total rainfall recorded during the first and second season is 75.8 cm and 79.7 cm, respectively.
In the experimental field, rain-gauge, RBC Parshall Flume, evapotranspiration (ET)-gauge, and 5TE sensors with EM 50 data logger were installed as shown in Figure 1; Figure 2. The sprinkler irrigation system was installed in the field for supplement irrigation. The amount of precipitation was measured using a data-logging rain-gauge. The amount of water drained from the field was measured by using Parshall Flume RBC with CR200X logger (SZ-CR200X/7070) whenever the runoff event occurs at the outlet point. ET gauge was used for evapotranspiration (ET) measurements. However, the evapotranspiration (ET) calculations were also made using the Penman–Monteith equation. The daily climate data required for ET calculations, i.e., temperature, humidity, radiation, and wind speed were obtained from the adjacent meteorological station. Soil moisture measuring sensors (5TE) with EM 50 data logger was installed at 20, 40, 60, 80, and 140 cm depths. TDR-300 was used for monitoring soil moisture of topsoil. For a few samples, soil moisture measurement was done through the oven dry method in the laboratory. Laboratory and TDR-300 measurements also helped in calibrating sensor measurements.
Two observation wells were installed to monitor groundwater table daily. The water fluctuated with rainfall occurring, otherwise, it remained around 150 cm. Three tensiometers were also installed at depths 20, 40, and 60 cm to observe pressure head within the root zone.
Soil samples were collected from the field at different depths of 0–20, 20–40, 40–60, 60–80, and 80–140 cm to determine the physical and chemical properties of the soil. The standard pipette method was used for soil texture analysis, which is considered as a reference method to classify soil texture. Table 1 presents the percentages of sand, silt, and clay of soil. The textural analysis was performed at Soil Physics Laboratory, Faculty of Agriculture, UPM. For the chemical properties, soil samples were tested at United Plantations Laboratory, Teluk Intan, Perak.
The subsurface water collection system was designed and installed to collect water percolated below the root zone (Figure 3). The concrete block of 1 m3 was buried in soil for this purpose. Soil bulk density of a block filled with soil equal to field bulk density. The blocks drain was connected to a bucket with a 72-L capacity to collect water leached in the block. The collected water was pumped out at the end of each rainfall with the help of a small-diaphragm pump.

2.1.3. Water Balance Components

The soil water balance equation used for determining available water in a sweet corn field as follows:
SS   =   WI     RO     ET     L ,
where, SS is the available soil storage (cm) determined using 5TE sensor; WI is the amount of total water input (cm) recorded using rain gauge; RO is the amount of runoff water from field (cm) measured using RBC flume; ET is the amount of evapotranspiration (cm) observed using ET gauge; L is the amount of water leaching (cm) calculated using the equation.

2.2. HYDRUS-ID Model

2.2.1. Model Description

The one-dimensional numerical model HYDRUS-1D [13] was used to simulate soil water dynamics in the sweet corn field. As the water flows from surface to groundwater is in a vertical direction under rainfall, and horizontal flow such as runoff can be covered through boundary conditions, there was no need to simulate multi-dimensional movement. The HYDRUS-1D solves the Richards equation for saturated and unsaturated water flow. Following is the governing equation for water flow, which is the modified form of Richards equation assuming water flow due to thermal gradient is neglected.
θ t   =   z K h z   +   1     S ,
where θ is the volumetric water content (cm3cm−3), h is the water pressure head (cm), t is time (day), z is the spatial coordinate (cm), K is the unsaturated hydraulic conductivity function (cm day−1), and S is the sink term in the flow equation (cm3cm−3day−1).
The soil hydraulic properties were modeled using [24] as follows:
θ h   =   θ r   +   θ s     θ r 1   +   α h n m h   <   0 θ s h     0 ,
K h   =   K s S e l 1     1     S e 1 / m m 2 .
Here
S e   =   θ     θ r θ s     θ r ,   m   =   1     1 n .
Here, θs is the saturated water content (cm3cm−3), θr is the residual water content (cm3cm−3), K s is the saturated hydraulic conductivity (cm day−1), S e is the relative saturation, and a (cm−1), and n (–) and l (–) are shape parameters.
The details of simulation process and model description have been presented in the flow chart (Figure 4).

2.2.2. HYDRUS-1D Input Parameters

Estimation of Soil Hydraulic Parameters

The soil hydraulic parameters such as θs, θr, m, a, and n (Table 3) were estimated using the ROSSETA computer program for the HYDRUS-1D input. The value of the pore connectivity parameter (l) was assumed to be (0.5) [13].

Estimation of Potential ET

HYDRUS-1D requires potential evaporation and potential transpiration fluxes as input values to simulate soil water effect on transpiration. The Penman–Monteith equation [25] was used to calculate the reference evapotranspiration ( E T 0 ) using available soil, crop, and climate data. The crop potential evapotranspiration ( E T c ) was calculated from the reference evapotranspiration ( E T 0 ) and the crop coefficient ( K c ) as given in Equation (7) under normal conditions [25]. The Kc values for four stages (early, development, mid, and late seasons) of sweet corn were 0.4, 0.8, 1.15, and 1, respectively [26]. Potential evapotranspiration (ETP) was divided into potential evaporation (EP) and Potential transpiration (TP) using the AquaCrop model.
E T 0   =   0.408 Δ R n     G   +   γ 900 T   +   273 u 2 e s     e a Δ   +   γ 1   +   0.34 u 2 ,
where, E T 0 = reference crop evapotranspiration (mmday−1), Rn = net radiation at the crop surface (MJm−2day), G = soil heat flux density (MJm−2day), T = air temperature at 2 m height (°C), u2 = wind speed at 2 m height (ms−1), e s     e a = saturation vapour pressure deficit (kPa), Δ = slope vapor pressure curve (kPa/°C), γ = psychometric constant (kPa/°C), and 900 = conversion factor
E T c   =   E T 0   ×   K c .

Root Water Uptake

The root water uptake was determined from [27].
S   =   α h S m a x ,
where S is the rate of root water uptake, S m a x is the maximum rate of root water uptake, and α h is the coefficient of root water uptake. The values of coefficients ho, hopt, h2H, h2L, and h3 were taken as −10, −25, −500, −1000, and −16,000 respectively [13]. It is assumed that the water uptake is zero near saturation and for pressure heads less than the wilting point (h3). The root water uptake is assumed optimal between pressure head hopt and h2, whereas it shows a linear increase or decrease for pressure heads between h2 and h3.

2.2.3. Initial and Boundary Conditions

The initial conditions were defined in water contents. The soil profile was subjected to rainfall and evaporation. Therefore, the upper boundary conditions were defined as the atmosphere boundary condition (BC) with runoff and lower boundary condition was opted as variable pressure head. The measured and recorded values of EP, TP, and rainfall fluxes were used to represent atmospheric boundary conditions.

2.3. Model Evaluation

To evaluate the model performance, the agreement between observed and modeled data was assessed through the statistical procedure. The coefficient of determination R2, modeling efficiency (EF) and root mean square error (RMSE) were used for this purpose, which are given as [28].
Coefficient of determination
R 2   =   i = 1 n O i     O ¯ P i     P ¯ i = 1 n O i     O ¯ 2 . 5 i = 1 n P i     P ¯ 2 . 5 2 .
Modeling efficiency
E F   =   1     i = 1 n O i     P i 2 i = 1 n O i     O ¯ 2 .
Root Mean Square Error
R M S E   =   i = 1 n O i     P i 2 n ,
where P i are predicted values; P ¯ is mean of predicted values; O are observed values; and O ¯ is mean of observed values. The optimum value of R 2 and E F is 1 while for R M S E it is 0.

2.4. Model Calibration and Validation

Data sets of the first and second seasons were used to calibrate and validate the model, respectively. HYDRUS-1D requires soil hydraulic parameters as input. The required parameters were determined through ROSETTA computer program [29] using field soil data. The determined soil hydraulic parameters θ s , θ r , α , n , and K s were further calibrated using observed soil water contents of the first season. Observed and simulated soil water content were found in good agreement. The calibrated parameters were then used to simulate soil water contents of the second season, for model validation.
The model was tested by comparing observed water content values with HYDRUS-1D simulated values using statistical indicators including the determination coefficient ( R 2 ), modelling efficiency ( E F ) and the root mean square error R M S E for both seasons as shown in Table 4. The average values of the goodness of fit indicators: R2, EF, and RMSE for depths of 20, 40, 60, and 80 cm were 0.72, 0.68, and 0.017, and 0.86, 0.73, and 0.014 for the first and second seasons, respectively. The goodness of fit indicators was within the range of values reported by different researchers [22,30,31,32,33,34]. The statistical analysis shows that HYDRUS-1D captured the change in water content at all depths except at 140 cm depth. The reason for that might be the no variation in simulated water content due to the shallow water table. Therefore, the relative change simulated water content w.r.t observed water content could not be predicted.

3. Results and Discussions

3.1. Evapotranspiration (ET)

Figure 5 presents the daily evapotranspiration data for the first and second season. The mean simulated evapotranspiration (ETm) values for the first and second seasons were 4.1 mm and 3.7 mm, respectively. The mean observed evapotranspiration (ET) values for the first and second seasons were 4.0 mm and 3.7 mm, respectively. The minimum and maximum evapotranspiration during the first season were recorded as 1.6 mm and 6.2 mm whereas for the second season the values were 1.3 mm and 5.9 mm, respectively. The evapotranspiration was low at early stages of crop due to less transpiration requirements by crop canopy, whereas it reached to the maximum between 40 and 60 days of crop plantation. The variation in the evapotranspiration (ET) values for peak days might be due to the difference in sunshine hours and duration of rainfall occurrence. The rainfall durations in the region are long, the cloudy sky also affects the evaporation. The evapotranspiration (ET) of total water use was recorded as 40.7% and 33.1% for first and second season, respectively, i.e., the rainfall amount was much higher than total plant water requirements for both seasons. The ratio of evaporation to evapotranspiration (Ea/ETa) for crop during two seasons were recorded as 0.24 and 0.26 respectively. The computed values obtained are close to reported research works [16,35].

3.2. Water Contents

The comparison between simulated and observed values of water contents in the soil profile within the root zone (0–60 cm) during two crop seasons shown in Figure 6. During the first and second season of field crop the simulated water content values at 60 cm depth ranged from 0.31 to 0.46 cm3 cm−3 and 0.31 to 0.49 cm3 cm−3 with average values of 0.41 and 0.43 cm3 cm−3 respectively. Whereas observed water content values for the first season and second season of field crop at 60 cm soil depth ranged from 0.32 to 0.45 cm3 cm−3 and 0.32 to 0.46 cm3 cm−3 with average values of 0.40 and 0.41 cm3 cm−3 respectively. The comparison between observed and simulated values of water contents in the soil profile below the root zone (60–140 cm) is illustrated in Figure 7. The values of water contents at different depths increased with rainfall events and then decreased until the next event due to evapotranspiration (ET). This variation of water content was larger in shallower depths as compared to deeper [30]. Simulated water content results at 140 cm soil depth were found constant for both seasons due to shallow water table, i.e., soil remained at field capacity (FC) or near saturation for all growing season. Whereas, a small fluctuation in observed water contents was recorded for 140 cm depth. The average simulated value at 140 cm depth for both seasons was 0.42 cm3 cm−3 whereas average observed values were 0.43 cm3 cm−3 and 0.44 cm3 cm−3 respectively. The simulated and observed water contents remained above the permanent wilting point (PWP, i.e., in the range of plant available water. The performance of HYDRUS-1D in predicting water content was found satisfactory.

3.3. Surface Runoff

Figure 8 shows the comparison of simulated and observed values of surface runoff for two consecutive seasons. The results for cumulative simulated and observed values were found in good correspondence. However, the results differ for some of the individual rainfall events. The first few low-intensity rainfall events do not produce runoff in actual field conditions while HYDRUS-1D simulation shows runoff for those events for both seasons. The reason for deviation seems to be fine soil tilth at the time of crop plantations that enhance the water holding capacity of the soil. The maximum rainfall value during the first season was recorded as 11.3 cm on the 33rd day of the plantation, which produced peak simulated and observed runoff of 9.61 and 9.12 cm, respectively. Whereas the maximum rainfall value during the second season was recorded as 6.47 cm on the 29th day of the plantation, that produced peak simulated and observed runoff of 4.92 and 6.05 cm, respectively. The difference in the ratio of runoff to total rainfall varies from event to event due to soil antecedent moisture level, rainfall amount, and intensity. Surface runoff accounted for 41% and 28.6% of total water input for first and second season respectively. Effective runoff simulation also shows HYDRUS-1D ability to handle the upper boundary condition.

3.4. Root Water Uptake

Figure 9 shows the simulated results for root water uptake (transpiration). Daily and cumulative root water uptake results for two consecutive seasons. Daily uptake mainly depends on the growth stage of the plant. It shows relatively lower values at the initial stage and maximum values at peak growth stages. For the time between sowing to emergence, the water uptake values were 0 for both seasons, while maximum values for the first and second season recorded on 58th and 52nd day after plantation were 0.56 cm/day and 0.54 cm/day, respectively. After the 60th day of the plantation, the root water uptake started to decline during crop ripening stage for both seasons. However, as the crop was harvested in fresh form and it was not left in the field (as in grain corn) the decline in uptake was less. The total amount of root water uptake recorded 30.3% of total water input during the first season and 24% in the second season. This shows the availability of water is enough to meet crop water requirements.

3.5. Water Leaching

Daily as well as cumulative water leaching results for both seasons presented in Figure 10. The total amount of water percolated in the first and second seasons was 7.64 cm and 20.9 cm, respectively. Water leaching accounts for 10.6%–26.8% of total water input during both seasons respectively. The daily leaching results show negative values during early days of the experiments, i.e., days 12–15 (first season) and days 14–17 (second season), which reflects groundwater contribution for the early days because of capillary rise due to the shallow water table. The difference in water leaching between the two seasons was quite significant. Relatively longer and frequent rainfall events in the second season caused the high antecedent moisture level, and the high antecedent moisture level resulted in the increased water leaching. Some of the results of observed leaching losses through designed subsurface water collection systems were missing due to difficulty in setting up the system in the first season and leakage detection in the system in the second season. Therefore, they were not presented for comparison with simulated results. The reason that the concrete block was sealed from the bottom while the soil profile was exposed to the shallow water table in real field conditions, the difference in observed and simulated leaching was higher. The designed subsurface system can be an excellent tool to determine water leaching but needs calibration under low water table conditions to overcome the effect of capillary rise on water leaching.

3.6. Water Balance

Water balance components such as water input (WI), evapotranspiration (ET), surface runoff (RO), and leaching (L) were monitored through intensive field observations. Simulated and observed water balance components of the soil profile are presented in Table 5.
The total water inputs during the two seasons were 75.8 cm and 79.7, respectively. Simulation results of evapotranspiration (ET) for the first and second season were 40.7% and 33.1% of total water input, which shows the rainfall amount was much higher than total plant water requirements for both seasons. Surface runoff accounted for 41% and 28.6% of total water input for the first and second season, respectively. Water leaching accounts for 10.6% to 26.8% of total water input during both seasons respectively. The total water balance error for both seasons was recorded as 2.48 cm and 3.70 cm, respectively, which accounted for 4.8% and 7.4% of the total water input. This shows the model’s effectiveness in simulating the water balance components.

4. Conclusions

The HYDRUS-1D simulation model was used to evaluate water flow and water losses during two growing seasons of sweet corn in MARDI, Malaysia under the tropical rainfed conditions in the presence of a shallow groundwater table. HYDRUS-1D (despite considerable demand for input data) has proven to be a reliable tool to evaluate water movement in agricultural fields under various irrigation schemes and different crops [17,31,32]. Similarly, HYDRUS-1D was found useful in evaluating water fluxes and to evaluate the HYDRUS-1D model for the rainfall potential to meet the water requirement of sweet corn in the tropical region. The sprinkler system was set up to deal with water-stressed conditions. As rainfall fulfilled the crop water requirement throughout the growing seasons no additional irrigation (sprinkler) was applied. The soil water content and crop uptake did not drop below the threshold level. The results suggest avoiding scheduled irrigation, which is the general practice in Malaysia. Excess irrigation triggers the runoff and leaching losses that result in fertilizer loss. However, the study was carried out in the presence of shallow water table that also contributes to soil water through the capillary rise. Therefore, it is recommended to repeat the research with deep water table scenario to establish the findings for the tropical region.

Author Contributions

M.I. carried out the experimental investigation and prepared the manuscript, including the data analysis. M.R.K. supervised the research work and improved the manuscript. M.F.M. assisted for the field experiment, H.C.M. and A.W. provided valuable comments for preparing the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding

Acknowledgments

Authors are grateful to Universiti Putra Malaysia for providing research facilities under Putra Grant (Grant No: 9557200). Authors also acknowledge the University of Agriculture Faisalabad and Higher Education Commission (HEC) Islamabad, Pakistan, for providing the scholarship in the PhD study. Authors sincerely express the gratitude to MARDI for providing experimental field.

Conflicts of Interest

All the authors declare no conflict of interest.

References

  1. USDA. United States Department of Agriculture, World Agricultural Production; Circular Series; USDA: Washington, DC, USA, 2018; Volume WAP 12–18.
  2. Wahab, A.G. GAIN REPORT: Grain and Feed Annual 2018; USDA Foreign Agricultural Service: Washington, DC, USA, 2018.
  3. DOA (Department of Agriculture). Vegetables and Cash Crops Statistic; Department of Agriculture: Putrajaya, Malaysia, 2017.
  4. Saeed, I.; Abassi, K.M.; Kazmi, M. Response of Maize (Zea mays, L.) to NP Fertilization under Agro-climatic Conditions of Rawalakot Azad Jammu and Kashmir. Park. J. Biol. Sci. 2001, 4, 949–952. [Google Scholar]
  5. Shamshuddin, J.; Sharifuddin, H.A.H.; Che Fauziah, I.; Edwards, D.G.; Bell, L.C. Temporal changes in chemical properties of acid soil profiles treated with magnesium limestone and gypsum. Pertanika J. Trop. Agric. Sci. 2010, 33, 277–295. [Google Scholar]
  6. Oad, F.C.; Buriro, U.A.; Agha, S.K. Effect of organic and inorganic fertilizer application on maize fodder production. Asian J. Plant. Sci. 2004, 3, 375–377. [Google Scholar]
  7. Shelia, V.; Šimunek, J.; Boote, K.; Hoogenbooom, G. Coupling DSSAT and HYDRUS-1D for simulations of soil water dynamics in the soil-plant-atmosphere system. J. Hydrol. Hydromech. 2018, 66, 232–245. [Google Scholar] [CrossRef] [Green Version]
  8. Igbadun, H.E. Estimation of Crop Water Use of Rain-Fed Maize and Groundnut Using Mini-Lysimeters. Pac. J. Sci. Technol. 2012, 13, 527–535. [Google Scholar]
  9. Wang, W.; Peng, S.; Yang, T.; Shao, Q.; Xu, J.; Xing, W. Spatial and temporal characteristics of reference evapotranspiration trends in the Haihe River basin, China. J. Hydrol. Eng. 2011, 16, 239–252. [Google Scholar] [CrossRef]
  10. Bruijnzeel, L.A. Predicting the hydrological impacts of land cover transformation in the humid tropics: The need for integrated research. In Amazonian Deforestation and Climate; John Wiley & Sons: New York, NY, USA, 1996; pp. 15–55. [Google Scholar]
  11. Watts, D.G.; Martin, D.L. Effects of water and nitrogen management on nitrate leaching loss from sands. Trans. ASAE 1981, 24, 911–916. [Google Scholar] [CrossRef]
  12. Šimůnek, J.; van Genuchten, M.T.; Šejna, M. Development and Applications of the HYDRUS and STANMOD Software Packages and Related Codes. Vadose Zone J. 2008, 7, 587–600. [Google Scholar] [CrossRef] [Green Version]
  13. Li, Y.; Šimůnek, J.; Zhang, Z.; Jing, L.; Ni, L. Evaluation of nitrogen balance in a direct-seeded-rice field experiment using Hydrus-1D. Agric. Water Manag. 2015, 148, 213–222. [Google Scholar] [CrossRef] [Green Version]
  14. Karandish, F.; Šimůnek, J. Two-dimensional modeling of nitrogen and water dynamics for various N-managed water-saving irrigation strategies using HYDRUS. Agric. Water Manag. 2017, 193, 174–190. [Google Scholar] [CrossRef] [Green Version]
  15. He, K.; Yang, Y.; Yang, Y.; Chen, S.; Hu, Q.; Liu, X.; Gao, F. HYDRUS simulation of sustainable brackish water irrigation in a winter wheat-summer maize rotation system in the North China Plain. Water 2017, 9, 536. [Google Scholar] [CrossRef] [Green Version]
  16. Ren, D.; Xu, X.; Hao, Y.; Huang, G. Modeling and assessing field irrigation water use in a canal system of Hetao, upper Yellow River basin: Application to maize, sunflower and watermelon. J. Hydrol. 2016, 532, 122–139. [Google Scholar] [CrossRef]
  17. Li, Y.; Šimůnek, J.; Jing, L.; Zhang, Z.; Ni, L. Evaluation of water movement and water losses in a direct-seeded-rice field experiment using Hydrus-1D. Agric. Water Manag. 2014, 142, 38–46. [Google Scholar] [CrossRef] [Green Version]
  18. Gabiri, G.; Burghof, S.; Diekkrüger, B.; Leemhuis, C.; Steinbach, S.; Näschen, K. Modeling spatial soilwater dynamics in a tropical floodplain, East Africa. Water 2018, 10, 191. [Google Scholar] [CrossRef] [Green Version]
  19. Ursulino, S.; Maria, S.; Lima, G.; Coutinho, A.P. Modelling soil water dynamics from soil hydraulic parameters estimated by an alternative method in an experimental basin located in the Brazilian Northeast region. Water 2019, 11, 1007. [Google Scholar] [CrossRef] [Green Version]
  20. Martello, M.; Dal Ferro, N.; Bortolini, L.; Morari, F. Effect of incident rainfall redistribution by maize canopy on soil moisture at the crop row scale. Water 2015, 7, 2254–2271. [Google Scholar] [CrossRef] [Green Version]
  21. Hou, L.; Zhou, Y.; Bao, H.; Wenninger, J. Simulation of maize (Zea mays L.) water use with the HYDRUS-1D model in the semi-arid Hailiutu River catchment, Northwest China. Hydrol. Sci. J. 2017, 62, 93–103. [Google Scholar]
  22. Negm, A.; Capodici, F.; Ciraolo, G.; Maltese, A.; Provenzano, G.; Rallo, G. Assessing the Performance of Thermal Inertia and Hydrus Models to Estimate Surface Soil Water Content. Appl. Sci. 2017, 7, 975. [Google Scholar] [CrossRef] [Green Version]
  23. Todorovic, M.; Lamaddalena, N.; Jovanovic, N.; Pereira, L. Agricultural water management: Priorities and challenges. Agric. Water Manag. 2014, 147, 1–3. [Google Scholar] [CrossRef] [Green Version]
  24. Van Genuchten, M.T. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci. Soc. Am. J. 1980, 44, 892–898. [Google Scholar] [CrossRef] [Green Version]
  25. Allen, R.G.; Luis, S.P.; Raes, D.; Smith, M. FAO Irrigation and Drainage Paper No. 56. Crop Evapotranspiration (guidelines for computing crop water requirements). Irrigation and Drainage 1998, 300, 300. [Google Scholar]
  26. Brouwer, C.; Heibloem, M. Irrigation Water Management: Irrigation Water Needs; Training Manual; Food and Agriculture Organization: Rome, Italy, 1986; Volume 3. [Google Scholar]
  27. Feddes, R.A. Simulation of Field Water Use and Crop Yield; Pudoc: Wageningen, The Netherlands, 1978. [Google Scholar]
  28. Nash, J.E.; Sutcliffe, J.V. River flow forecasting through conceptual models part I—A discussion of principles. J. Hydrol. 1970, 10, 282–290. [Google Scholar] [CrossRef]
  29. Schaap, M.G.; Leij, F.J.; Van Genuchten, M.T. Rosetta: A computer program for estimating soil hydraulic parameters with hierarchical pedotransfer functions. J. Hydrol. 2001, 251, 163–176. [Google Scholar] [CrossRef]
  30. González, M.G.; Ramos, T.B.; Carlesso, R.; Paredes, P.; Petry, M.T.; Martins, J.D.; Aires, N.P.; Pereira, L.S. Modelling soil water dynamics of full and deficit drip irrigated maize cultivated under a rain shelter. Biosyst. Eng. 2015, 132, 1–18. [Google Scholar] [CrossRef] [Green Version]
  31. Ramos, T.B.; Šimůnek, J.; Gonçalves, M.C.; Martins, J.C.; Prazeres, A.; Pereira, L.S. Two-dimensional modeling of water and nitrogen fate from sweet sorghum irrigated with fresh and blended saline waters. Agric. Water Manag. 2012, 111, 87–104. [Google Scholar] [CrossRef]
  32. Mo’allim, A.; Kamal, M.; Muhammed, H.; Yahaya, N.; Zawawe, M.; Man, H.; Wayayok, A. An assessment of the vertical movement of water in a flooded paddy rice field experiment using hydrus-1D. Water 2018, 10, 783. [Google Scholar] [CrossRef] [Green Version]
  33. Ramos, T.; Šimůnek, J.; Gonçalves, M.; Martins, J.; Prazeres, A.; Castanheira, N.; Pereira, L. Field evaluation of a multicomponent solute transport model in soils irrigated with saline waters. J. Hydrol. 2011, 407, 129–144. [Google Scholar] [CrossRef]
  34. Phogat, V.; Skewes, M.A.; Cox, J.W.; Alam, J.; Grigson, G.; Šimůnek, J. Evaluation of water movement and nitrate dynamics in a lysimeter planted with an orange tree. Agric. Water Manag. 2013, 127, 74–84. [Google Scholar] [CrossRef] [Green Version]
  35. López-Urrea, R.; Montoro, A.; Trout, T. Consumptive water use and crop coefficients of irrigated sunflower. Irrig. Sci. 2014, 32, 99–109. [Google Scholar] [CrossRef]
Figure 1. Field layout of equipment installed for investigation of water balance components.
Figure 1. Field layout of equipment installed for investigation of water balance components.
Applsci 10 01219 g001
Figure 2. A view of the installation of different equipment in the field.
Figure 2. A view of the installation of different equipment in the field.
Applsci 10 01219 g002
Figure 3. View of the sub-surface water collection system in the field experiment.
Figure 3. View of the sub-surface water collection system in the field experiment.
Applsci 10 01219 g003
Figure 4. Flow chart presenting model simulation process and calculation of water balance components.
Figure 4. Flow chart presenting model simulation process and calculation of water balance components.
Applsci 10 01219 g004
Figure 5. Simulated and observed daily crop evapotranspiration for sweet corn during two growing seasons.
Figure 5. Simulated and observed daily crop evapotranspiration for sweet corn during two growing seasons.
Applsci 10 01219 g005
Figure 6. Simulated and observed soil water contents within the root zone (0–60 cm) during the (a) first season and (b) second season.
Figure 6. Simulated and observed soil water contents within the root zone (0–60 cm) during the (a) first season and (b) second season.
Applsci 10 01219 g006aApplsci 10 01219 g006b
Figure 7. Simulated and observed soil water contents below the root zone (60–140 cm) during the (a) first season and (b) second season.
Figure 7. Simulated and observed soil water contents below the root zone (60–140 cm) during the (a) first season and (b) second season.
Applsci 10 01219 g007
Figure 8. Simulated and observed surface runoff during two growing seasons.
Figure 8. Simulated and observed surface runoff during two growing seasons.
Applsci 10 01219 g008
Figure 9. Simulated daily and cumulative root water uptake during two growing seasons.
Figure 9. Simulated daily and cumulative root water uptake during two growing seasons.
Applsci 10 01219 g009
Figure 10. Simulated water leaching losses during two growing seasons.
Figure 10. Simulated water leaching losses during two growing seasons.
Applsci 10 01219 g010
Table 1. Physical properties of the soil in the experimental field at Malaysian Agricultural Research and Development Institute (MARDI).
Table 1. Physical properties of the soil in the experimental field at Malaysian Agricultural Research and Development Institute (MARDI).
Depth (cm)Clay %Silt %Sand %Textural ClassBulk Density (g cm−3)Conductivity Ks (cm d−1)
0–2065.3609.7924.84Clay1.3612.86
20–4053.2815.2431.45Clay1.3912.43
40–6058.6613.1228.23Clay1.3112.78
60–8070.2020.2209.58Clay1.2218.11
80–14071.7418.0310.23Clay1.2417.56
Table 2. Chemical properties of the soil in the experimental field at MARDI.
Table 2. Chemical properties of the soil in the experimental field at MARDI.
Depth (cm)pH (Airdry)EC
umoh/cm
CEC
cmol (+)/kg
Organic Matter %NP
mg/Kg
K
mg/Kg
0–206.8762.805.880.145574.1
20–405.9564.905.680.122674.1
40–606.2045.756.490.123919.5
60–805.8041.477.470.113042.9
80–1405.9885.787.130.092219.5
Table 3. Optimized value of soil hydraulic properties.
Table 3. Optimized value of soil hydraulic properties.
Depth (cm)Soil Typeθs
(cm3cm−3)
θr (cm3cm−3)a (cm−1)nKs
(cm Day −1)
L
0–20Clay0.09730.47500.02501.162512.860.5
20–40Clay0.09300.45610.02521.201912.430.5
40–60Clay0.09520.46610.02541.182812.780.5
60–80Clay0.10200.50840.02011.190318.110.5
80–140Clay0.10170.50660.02031.182517.560.5
Table 4. Statistical analysis of the comparison between observed and simulated water contents.
Table 4. Statistical analysis of the comparison between observed and simulated water contents.
SeasonDepth
(cm)
R2EFRMSE
(cm3cm−3)
10–200.720.730.017
20–400.720.630.021
40–600.750.700.018
60–800.670.650.012
80–1400.51−8.240.014
20–200.840.700.013
20–400.850.820.014
40–600.940.720.017
60–800.820.690.013
80–1400.36−14.90.019
Table 5. Simulated and observed water balance components (cm) in the soil profile during two seasons of sweet corn in MARDI, Malaysia.
Table 5. Simulated and observed water balance components (cm) in the soil profile during two seasons of sweet corn in MARDI, Malaysia.
SeasonWIETROLSSƩ
FirstSimulated75.8−30.88−31.10−07.643.702.48
Observed −30.84−36.00
SecondSimulated79.7−26.40−22.80−20.905.903.70
Observed −26.70−22.40
WI—total input water, ET—evapotranspiration, RO—surface runoff, L—water leaching, SS—soil storage, Ʃ water balance error.

Share and Cite

MDPI and ACS Style

Iqbal, M.; Kamal, M.R.; M., M.F.; Che Man, H.; Wayayok, A. HYDRUS-1D Simulation of Soil Water Dynamics for Sweet Corn under Tropical Rainfed Condition. Appl. Sci. 2020, 10, 1219. https://doi.org/10.3390/app10041219

AMA Style

Iqbal M, Kamal MR, M. MF, Che Man H, Wayayok A. HYDRUS-1D Simulation of Soil Water Dynamics for Sweet Corn under Tropical Rainfed Condition. Applied Sciences. 2020; 10(4):1219. https://doi.org/10.3390/app10041219

Chicago/Turabian Style

Iqbal, Mazhar, Md Rowshon Kamal, Mohd Fazly M., Hasfalina Che Man, and Aimrun Wayayok. 2020. "HYDRUS-1D Simulation of Soil Water Dynamics for Sweet Corn under Tropical Rainfed Condition" Applied Sciences 10, no. 4: 1219. https://doi.org/10.3390/app10041219

APA Style

Iqbal, M., Kamal, M. R., M., M. F., Che Man, H., & Wayayok, A. (2020). HYDRUS-1D Simulation of Soil Water Dynamics for Sweet Corn under Tropical Rainfed Condition. Applied Sciences, 10(4), 1219. https://doi.org/10.3390/app10041219

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