Next Article in Journal
Application of Urea and Ammonium Nitrate Solution with Potassium Thiosulfate as a Factor Determining Macroelement Contents in Plants
Next Article in Special Issue
Simulation of Maize Growth Under the Applications of Brackish Water in Northwest China
Previous Article in Journal
Impacts of Farming Activities on Nitrogen Degradability under a Temperate Continental Monsoon Climate
Previous Article in Special Issue
Decision of Straw Deep Burial and Aluminum Sulfate Drip Irrigation in Soda Saline Soil Based on Grey Relation Analysis and TOPSIS Coupling
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Simulation of Drainage Volume and Nitrogen Loss Load in Paddy Fields under Different Irrigation and Drainage Modes and Hydrological Years

1
Heilongjiang Province Hydraulic Research Institute, Harbin 150080, China
2
College of Agricultural Science and Engineering, Hohai University, Nanjing 210098, China
3
College of Hydrology and Water Resources, Hohai University, Nanjing 210098, China
4
School of Water Conservancy & Civil Engineering, Northeast Agricultural University, Harbin 150030, China
*
Authors to whom correspondence should be addressed.
Agronomy 2024, 14(6), 1095; https://doi.org/10.3390/agronomy14061095
Submission received: 9 April 2024 / Revised: 16 May 2024 / Accepted: 19 May 2024 / Published: 22 May 2024

Abstract

:
Controlled irrigation and drainage technology for rice is crucial water management that has been widely promoted in northeastern China. It is of great significance to clarify the response mechanism of the drainage volume and nitrogen loss load in the paddy field for realizing water saving, emission reduction, pollution control and high yield in rice-planting areas. In this study, we conducted field experiments and simulations on drainage volume and nitrogen loss load regulations of paddy fields in a cold black soil region under different hydrological years and irrigation and drainage modes. The key parameters for simulating drainage volume and nitrogen loss load in paddy fields using DRAINMOD-NII were determined by combining field experiments, data analysis, and numerical simulation. The results showed that the simulated drainage volume and nitrogen loss load showed a high coefficient of determination with the observed results, which were all above 0.83. The Nash–Sutcliffe efficiency coefficient ranged from 0.72 to 0.97 in model calibration and verification, indicating that the model effectively simulated drainage volume and nitrogen loss load in paddy fields under controlled irrigation and drainage in the cold black soil region. The paddy field drainage volume was not only influenced by rainfall frequency but also by the distribution of rainfall. Compared with traditional irrigation and drainage, the controlled irrigation and drainage significantly reduced the irrigation amount by 39.07% and increased rainwater utilization efficiency by 13.07%. It also reduced the drainage volume by 44.71% and NO3-N and NH4+-N loss load by 59.38% and 44.96%. The controlled irrigation and drainage mode optimized natural rainfall resources and increased irrigation water productivity and rice yield by 97.85% and 16.88%, respectively. Controlled irrigation and drainage outperformed the traditional mode in different hydrological years, with more pronounced effects in dry years, which highlights its significant value in practical agricultural production.

1. Introduction

Water and nitrogen in paddy fields are key factors affecting rice growth, and the rational regulation of water and nitrogen can promote rice growth and improve water use efficiency [1,2]. When the soil moisture in the paddy field is sufficient, the water and nutrients are coordinated with each other to promote the release of nutrients to be absorbed and utilized by plants in time [3,4]. A field plot experiment is an important means to quantitatively study water and nitrogen transport in paddy fields, and the water transport model is an effective method to characterize the dynamic changes of water and nitrogen in paddy fields [5,6,7]. Most of the nitrogen fertilizer applied in paddy fields is partially absorbed by plants, and the ammonia volatilization accounted for about 14.0% [8]. Analysis of nitrogen loss load in paddy fields through field experiments showed that the average concentration of total nitrogen in paddy field drainage was about 6.60 mg·L−1, with an average nitrogen loss of about 16.70 kg·hm−2 [9]. Although the annual nitrogen loss load in agricultural ecosystems in China is constantly changing, the situation of nutrient loss load is not optimistic [10,11,12]. Controlled drainage improves rainwater utilization and reduces nitrogen loss in drainage by reducing paddy field drainage [13]. Lu and Xiao found that extending drainage time after rainfall effectively reduced nitrogen loss in paddy fields and increased fertilizer use efficiency [14,15].
At present, many mathematical equations directly explain the process of surface runoff and nitrogen migration in soil. The EPIC model proposed by Li et al. [16] accurately simulated soil nitrogen and phosphorus transport and crop uptake processes. The CREAMS model [17] can simulate and evaluate the effects of different tillage measures on farmland surface runoff migration and the quantitative description of non-point source pollution load. The CREAMS-PADDY coupling model proposed by Chung et al. [18] was used to simulate and predict the loss of nitrogen and phosphorus in paddy field surface water. Hao et al. [19] simulated the process of soil water and nitrogen transport and transformation in wheat fields by the HYDRUS-1D model. Li et al. [20] simulated the process of nitrogen migration and transformation in paddy fields by the HYDRUS-1D model under the premise of considering the influencing factors such as irrigation, rainfall, surface runoff and infiltration in paddy fields. The Agricultural Research Center of the United States Department of Agriculture proposed the SWAT model [21], which can use the spatial data provided by the geographic information system to simulate the migration and transformation of water quantity, water quality and solute in different hydrological processes [22]. Li et al. [23] used the DNDC model to quantitatively analyze the main sources of nitrogen pollutants in farmland. Moursi et al. [24] improved the drainage model and developed the DRAINMOD-DWR model to simulate the hydrology and crop yield of the water cycle system.
The DRAINMOD-NII model improved by Skaggs et al. [25,26,27] is a field-scale model based on process, which is used to simulate the dynamic migration process of nitrogen under different field management measures and different environmental and soil conditions. The model is based on a simple water balance principle and a complete nitrogen transport process. The number of parameter inputs is small and easy to obtain, and the model is easy to understand and operate. In recent years, many researchers have simulated farmland irrigation and drainage, nutrient and salt transport using the DRAINMOD-NII model, which provides a reference for optimizing agricultural management [28,29]. The research results of Salazar et al. [30] found that the model effectively simulated the actual drainage and nitrogen loss of farmland in southeastern Sweden by comparing the measured values of drainage and nitrogen loss with the simulation results of the model. Li et al. [31] simulated shallow groundwater utilization and salt accumulation using the model. Jia et al. [32] studied the impact of different drainage intensities on wetland hydrological processes using the mentioned model. Luo et al. [33,34] used the model to simulate the drainage in the upper reaches of the Yellow River and the hydrological processes of maize–soybean rotation in agricultural fields. The research conducted by Wang et al. [35] showed that they applied the model to simulate changes in groundwater level and drainage volume; Zhang et al. [36] used the calibrated model to simulate groundwater depth and salt content in soil profiles, and the simulation results were used for prediction of improvement in saline–alkali land. Abduljaleel, Y. et al. [37] simulated excess water in poorly drained irrigated areas through experiments to protect them from induced waterlogging and salinization problems.
In contrast to the extensive results available for water-saving irrigation or controlled drainage in paddy fields, information on drainage and nitrogen loss under controlled irrigation and drainage in the cold black soil region is scarce. In particular, the features of drainage and nitrogen loss in different hydrological years have not been reported. In this study, water depth in paddy fields was used as an indicator for irrigation and drainage regulation. By combining field plot experiments, theoretical analysis and numerical simulation, the study explores and simulates the drainage and nitrogen loss load during the whole growth period of rice under different irrigation and drainage modes, which is expected to provide a basis for water saving and emission reduction theory and technology development in the cold black soil region.

2. Materials and Methods

2.1. Experimental Site and Soil Properties

The experiments were conducted at the Station of National Agriculture Irrigation Experiment in Qingan County, which is in the southwestern part of Heilongjiang Province, northeastern China (127°30′04″ E, 46°52′41″ N). The effective accumulated perennial temperature (≥10 °C) is 2518 °C. The mean annual precipitation is 545.3 mm, and the variation range is 450~700 mm. The average annual surface evaporation is 664.5 mm, and the variation range is 462~873 mm. The soil in the experimental site is a sandy loam with heavy texture and good structural characteristics. The soil porosity was 53.4%, the organic matter mass ratio was 37.5 g·kg−1, and the total nitrogen amount was 13.2 g·kg−1.

2.2. Experimental Design and Field Management

2.2.1. Experimental Design

A field experiment was conducted to study the paddy field drainage volume and nitrogen loss load in the paddy. The rice variety selected in this experiment was Suijing 18, which has a growth period of about 120 days and is widely cultivated in the local area. It was a fragrant rice variety that was mainly distributed in the northeast of China. This Japonica rice in the cold region demonstrated strong adaptability and resistance to lodging and is famous for its high quality and great taste [38,39,40]. The water depth in the paddy field served as the indicator for regulating irrigation and drainage. The test plots (10 m × 10 m) were separated by a 30 cm wide ridge to prevent water and nitrogen permeation between adjacent plots. A 50 cm deep film was embedded in the soil ridge to eliminate the effects of water and nitrogen leakage.
The experiment comprised two irrigation treatments, namely, controlled irrigation and drainage (CID: Controlled irrigation and drainage is a new irrigation technology which determines irrigation time, irrigation times and irrigation quota according to soil moisture content in the root layer or water depth in paddy fields in each growth period after seedling transplanting) and traditional irrigation and drainage (TID: Compared with controlled irrigation and drainage, traditional irrigation and drainage refers to the irrigation of crops and drainage of paddy fields by traditional agricultural irrigation methods. Irrigation can be divided into border irrigation, furrow irrigation, inundation irrigation and flood irrigation according to different ways of wetting soil.). Each experiment treatment had three repeated test plots, and the observed results of different monitoring indicators were the average values of the repeated plots. Each plot had its own independent irrigation and drainage system. Irrigation was carried out when the water depth in the test plot dropped to the lower limit of the paddy field (CID: 1 cm at replanting stage, −10 cm at tillering stage, −20 cm from jointing and booting stage to milking ripening stage. TID: 1 cm from replanting stage to milking ripening stage), and irrigation continued until the water layer depth reached the upper limit (CID: 3 cm from replanting stage to milking ripening stage. TID: 5 cm from replanting stage to milking ripening stage.). After rainfall, if the water depth in the test plot exceeded the maximum allowable water depth for the paddy field (CID: 10 cm from replanting stag to tillering stage, 10 cm from jointing and booting stage to milk ripening stage. TID: 5 cm from replanting stage to milk ripening stage), the excess water would be drained to reach the maximum allowable water depth after rain. The paddy fields of both TID and CID received the same fertilization. The water depth designs in the paddy field at different growth stages are shown in Table 1.

2.2.2. Field Management

The management of the experimental field is consistent with that of the local high-yield field. Seeds were sown in the greenhouse. The field test plots were flooded to a depth of 5 cm for approximately 1 week, and then the seedlings were transplanted to the paddy field at the 6-leaf stage with a transplanting size of 30 cm × 40 cm. The paddy field was fertilized three times: with base fertilizer (urea: 8.7 g·m−2; P2O5: 7.5 g·m−2; KCl: 6.0 g·m−2), tillering fertilizer (urea: 6.5 g·m−2) and ear fertilizer (urea: 6.5 g·m−2; KCl: 6.0 g·m−2), respectively. Please refer to Table 2 for specific agronomic measures implemented and the corresponding time in the paddy field.

2.3. Measurement Items and Methods

2.3.1. Meteorological Indicators

The meteorological data were automatically collected from a small automatic meteorological observation station (Watchdog 2000, Spectrum, Aurora, IL, USA) and nearby weather stations. The data included rainfall, evaporation, air temperature, relative humidity, wind speed, atmospheric pressure, total radiation, and net radiation. The temperature input format of the model took in the daily maximum and minimum temperatures. The rainfall data used in the model were measured by the meteorological monitoring equipment in the experimental site. The time step for rainfall was set to 1 h.
The daily evapotranspiration of DRAINMOD refers to the standard reference crop under high uniformity and sufficient water supply, rather than the actual crop. The potential evapotranspiration of the reference crop was calculated using the Thorthwaite formula in DRAINMOD and calibrated using monthly check coefficients. The input parameters of the model were revised to improve the accuracy of the simulation results. The potential evapotranspiration of reference crops was calculated using the Penman Monteith formula (Formula (1)), which is as follows [41]:
E T 0 = 0.408 × Δ × ( R n G ) + γ × 900 273 + T × u 2 × ( e a e d ) Δ + γ × ( 1 + 0.34 × u 2 )
where ET0 is the reference crop evapotranspiration, mm·d−1; ∆ is the slope of the saturated water vapor pressure curve, kPa/°C; Rn is the net radiation at the reference crop canopy surface, MJ·m−2·d−1; G is the soil heat flux, MJ·m−2·d−1; γ is the psychrometric constant, kPa/°C; T is the daily average temperature at 2 m, °C; u2 is the Wind speed at 2 m, m·s−1; ea is the saturated vapor pressure, kPa; and ed is the actual vapor pressure, kPa.
According to the rice crop coefficient in Table 3, the potential evapotranspiration of reference crops was modified, and the actual evapotranspiration of crops was calculated according to Formula (2).
E T = K c × E T 0
where ET is the actual crop evapotranspiration, mm·d−1; Kc is the crop coefficient; and ET0 is the reference crop evapotranspiration, mm·d−1.

2.3.2. Soil Physicochemical Parameters

A total of 72 soil samples from different depths (0~20, 20~40, 40~60, 60~80 cm) were taken with a ring knife at three points in the six test plots. The saturated water content and bulk density were measured using the ring knife method. The soil water characteristic curve was observed using a pressure membrane meter (NSSP-15, SEC, Goleta, CA, USA). The saturated water conductivity was measured using a soil water conductivity meter (KSAT, METER, MUC, Germany). The soil particle gradation was measured using sieving analysis and the siphon gravity bottle. The content of organic matter was determined by the volumetric method. The total nitrogen was analyzed by the alkaline potassium persulfate digestion method [42] (MEPC, 2002). The soil physical properties and hydraulic parameters of different soil layers in the test area are shown in Table 4.

2.3.3. Irrigation and Drainage

(1)
Water Depth
A paddy water level observation well is set up in each experimental plot to measure the water depth of the paddy field. The observation well is composed of a PVC pipe with a diameter of 4 cm and a length of 200 cm, and small holes are staggered at the end of the PVC pipe, for 50 cm, and then wrapped with geotextile and buried in the field. The outer wall of the PVC pipe is in close contact with the soil, and the pipe mouth is about 30 cm from the ground. The water level is observed regularly at 8:00 am every day. The soil surface in the field plot serves as the reference point, and the water level is positive above the ground line and negative below the ground line. The vertical distance from surface to water level in each field plot was measured, accurate to 1 mm. According to the irrigation test specification, the observation well is located in the middle of the field test plot and 1 m away from the ridge.
(2)
Irrigation and Drainage
The irrigation amounts were calculated using the water meter. Drainage occurred according to the test requirements when the water depth of the paddy field exceeded the allowed depth after rainfall. The water depths before and after drainage were recorded to calculate drainage volume.
(3)
Nitrogen in Water Samples
The water samples were manually collected in 500 mL plastic bottles every 10 min at the outlet of the plot when drainage occurred, and they were stored in a container and transported to the laboratory for analysis. Each water sample had three replicates. The concentrations of ammonia nitrogen (NH4+-N) and nitrate nitrogen (NO3-N) in the water samples were determined using the indophenol blue method and the disulfonic acid–phenol method, respectively (MEPC 2002). A UV spectrophotometer (Cary 50 Spectrophotometer, Varian, Palo Alto, CA, USA) was used for measurement [42]. The nitrogen concentrations were calculated based on the three replicates.

2.3.4. Yield

The theoretical yield of rice was calculated according to Formula (3) [43].
Ty = P × G × S × W × 0.85
Notes: Ty means the theoretical yield, kg·hm−2; P means effective panicle number, 10,000·hm−2; G means the number of grains per panicle; S means the average setting rate, %; and W means thousand grain weight, g.

2.4. Drainage Design Parameter of Model

DRAINMOD was initially developed as a field hydrological model for dry farmland in humid coastal alluvial plains. The parameter needed to be adjusted when the model was applied to paddy fields. The effective diameter of the drainage pipe was replaced with the hydraulic radius of the drainage ditch, and the irrigation data were included in the corresponding rainfall data. The design parameters for the drainage system in the model include both surface drainage and underground drainage. Surface drainage parameters consist of maximum and minimum surface water depths. Underground drainage parameters included drainage spacing, drainage depth, relative impervious depth, drainage coefficient, and initial groundwater level. The relative impervious depth and drainage coefficient were calculated according to the actual situation of the paddy field combined with the model calibration. Maximum surface water storage depth (Sm) and Krikham water depth (SI) were determined based on the surface roughness and actual drainage conditions of the paddy field and relevant literature [27]. The drainage depth and drainage spacing were measured in the test area (Table 5). The maximum nitration and denitrification rate and optimum temperature were derived from model calibration, and the NH4+-N and NO3-N concentration in the rain were measured in the test area (Table 6).

2.5. Model Evaluation Parameters

The model’s applicability was assessed by comparing simulated values with observed values using various indicators. Statistical parameters and a combination of numerical and graphical methods were utilized to evaluate the model’s suitability in simulating paddy field drainage and nitrogen loss load [44]. The absolute root mean square error (RMSEa), relative error (ε), Nash–Sutcliffe efficiency coefficient (NSE), and coefficient of determination (R2) were used to assess the approximation between the simulated and observed values of the drainage and nitrogen loss load [45]. The calculation formulas are as follows:
R M S E a = 1 n i = 1 n P i Q i 2
Q ¯ = 1 n i = 1 n Q i
ε = i = 1 n Q i i = 1 n P i / i = 1 n P i × 100 %
N S E = 1 i = 1 n P i Q i 2 / i = 1 n P i P ¯ 2
R 2 = i = 1 n P i P ¯ Q i Q ¯ 2 / i = 1 n P i P ¯ 2 i = 1 n Q i Q ¯ 2
where n is the number of samples, P represents the ith simulated value, Q represents the ith observed value, P ¯ represents the average of the simulated value sequence, and Q ¯ represents the average of the observed value sequence. RMSEa represents the degree of closeness between the simulated and observed values. R2 indicates the regression effect of the model. RMSEa and ε reflect the difference between the simulated and observed values. NSE reflects the simulation efficiency of the model.

2.6. Data Statistical Analysis Methods

The test data were processed using Microsoft Excel 2016 (Microsoft, Redmond, WA, USA). SPSS 17.0 software (IBM SPSS Statistics, Chicago, IL, USA) was used to perform standard deviation and other statistical analyses. Microsoft Excel 2010 and Origin 9.1 (OriginLab, Ampton, Kissimmee, FL, USA) were employed for data visualization and chart production.

3. Results

3.1. Model Calibration

The calibration process of model parameters is to adjust the parameter values, compare the model simulation values with the measured data in the same period, and adjust the parameters to make the simulation results infinitely close to the measured values to the greatest extent. The soil parameters, drainage volume, nitrogen loss, crop parameters and irrigation systems were obtained from a field test. The soil water characteristic curve was determined by a pressure membrane meter in the laboratory after sampling, the hydraulic conductivity was determined using a soil water conductivity meter in the field, and the other soil input parameters were deduced by a theoretical method according to the soil water characteristic curve. The horizontal saturated hydraulic conductivity has an obvious influence on the simulation results. In this paper, PEST was used to focus on the adjustment of horizontal saturated hydraulic conductivity. PEST (Parameter Estimation) is comprehensive software for parameter estimation and uncertainty analysis independent of the model. The control file (pest.pst), description file (dhp.ins) and template file (dhp.tpl) were established in the installation directory of PEST while running DRAINMOD-NII, and the check command of the relevant file was run in the DOS interface to determine whether the established file was normal. PEST uses a nonlinear evaluation method, which can not only converge quickly and effectively towards the objective function but also estimate the parameters with as few model runs as possible. After 51 iterations of the DRAINMOD-NII model, the parameters were optimized (Figure 1 and Table 7, Table 8 and Table 9).

3.1.1. Paddy Field Drainage

The rainfall during the whole crop growth period was 547.9 mm in 2020. The irrigation times were two, due to the low rainfall in the early growth period. The paddy field drainage occurred in the later period of rice growth because rainfall was increased and concentrated at this time. The drainage times were four, and the overall trend of the observed and simulated field drainage values was consistent, and the time of drainage peak simulation was basically consistent with the actual time of drainage under CID mode during different growth stages in the paddy field. However, the simulated drainage values are higher than the observed values for the four drainage instances in the paddy field; the observed value of the total drainage volume was 54.0 mm, and the simulated value was 68.0 mm. The maximum drainage occurred at the heading and flowering stage, with the observed value of 28.0 mm and the simulated value of 32.0 mm (Figure 2). This may be because the simulation process considers irrigation as equivalent rainfall and does not fully account for the actual drainage of the test field. It could also be related to the accuracy of manual measurement of paddy field drainage, as the human factors of the drainage process are not considered in the model.
The scatter plot in the graph represents the simulated drainage values corresponding to the observed values. The simulation results of the model overestimate the measured values, but the dispersion falls near the equality line (Y = X), indicating a reasonable alignment between the simulated and observed drainage volumes (Figure 3). The measured values of paddy field drainage are basically consistent with the simulated values, which are 54.0 mm and 68.0 mm, respectively. The R2 and NSE are 0.88 and 0.94, and the RMSEa and ε are 2.08 and 0.14, suggesting that the errors are within an acceptable range (Table 10). These indicate that the model performs well and effectively simulates the paddy field drainage volume under CID in the cold black soil region.

3.1.2. Nitrogen Loss Load

DRAINMOD-NII introduced the nitrogen cycle module, and the water balance result was used as the input value of nitrogen migration. The multi-phase one-dimensional convection–dispersion equation is used to model the nitrogen migration, and the first-order finite difference method is used to approximate the equation; then, the migration process of nitrogen in the soil unit is simulated. The simulated values of the NO3-N load in drainage were consistently higher than the observed ones, with the simulated and observed values being 0.75 kg·hm−2 and 0.60 kg·hm−2, respectively. The simulated results of the NH4+-N loss load were higher than the observed values, except for the heading and flowering stage, when the cumulative simulated value was lower than the observed value, with a simulated value and measured value of 4.10 kg·hm−2 and 3.60 kg·hm−2, respectively (Figure 4). The paddy field drainage volume was larger and the drainage time was longer due to the heavy rainfall at the heading and flowering stage. It was possible that the nitrogen concentration varied greatly in different periods due to digestion and denitrification during the drainage process. In order to improve the accuracy of the monitoring value, it is recommended to take more samples to calculate the average value. In addition, the inaccuracy may also be caused by the interference of human factors, leading to the increase in the error of the test index.
The scatter plots illustrate the relationship between the simulated and observed values of paddy field drainage. In the model calibration period, the scatter falls around the equality line (Y = X), indicating that the simulated values tend to be higher than the observed values (Figure 5). The ε and NSE for NO3-N and NH4+-N loss load cumulative values were 0.23, 0.88 and 0.80, 0.92, respectively. The RMSEa and R2 were 0.22, 1.09 and 0.83, 0.85, respectively. These values suggest that the model performs well, with errors falling within an acceptable range. The results indicate that the simulated values of NO3-N and NH4+-N loss load in paddy field drainage are consistent with the observed ones under the CID in 2020, indicating that the model performed well in simulating nitrogen loss load (Table 11).

3.2. Model Verification

In this paper, the climate data, crop data, paddy field drainage and nitrogen loss load data in 2021 were used to verify whether the soil parameters, drainage parameters and nitrogen transport module parameters of the model were reasonable. The verification effect of the model was evaluated by the graphical display method and statistical parameter index. In the graphical display method, the observed values of paddy field drainage and nitrogen loss load data were compared with the simulated values, and the simulation effect was judged intuitively from the graphical trend and the degree of coincidence. It shows that the model parameter value can represent the actual situation in this area, and the model parameter calibration is reasonable, if the simulated value is close to the measured value in the same period, and the error is within the acceptable range.

3.2.1. Paddy Field Drainage

The rainfall during the whole growth period of rice was 455.4 mm in 2021. The irrigation times were 4, which mainly occurred at the jointing and booting stage with less rainfall. The paddy field drainage occurred at the tillering stage because the water requirement in the paddy fields was less and rainfall was more. The results indicate a close consistency between the simulated and observed values, exhibiting similar peak values and occurrence times. The drainage times were 5, and the simulated value was 70.0 mm, while the observed value was 59.0 mm, resulting in a difference of 11.0 mm (Figure 6).
The scatter points were uniformly distributed around the equality line (Y = X), indicating a good fit between the simulated and observed values (Figure 7). The ε was 0.19, the R2 was 0.90, the NSE was 0.19, and the RMSEa was 2.72 (Table 12). Thus, results indicate a satisfactory calibration of the model with the simulated paddy field drainage processes being reasonably well modeled. Additionally, these results suggest that the model effectively represents the dynamics of paddy field drainage in the cold black soil region with accuracy.

3.2.2. Nitrogen Loss Load

In the model verification period, the simulated NO3-N loss load exceeded the observed values at all growth stages except for the milking ripening stage. Specifically, the cumulative simulated value for NO3-N loss load was 0.68 kg·hm−2, slightly lower than the measured value of 0.70 kg·hm−2. On the other hand, the simulated NH4+-N loss load was consistently higher than the observed values throughout the whole growth period. The total cumulative simulated value for NH4+-N loss load in drainage was 3.60 kg·hm−2, surpassing the observed result of 3.50 kg·hm−2 (Figure 8).
During the model verification period, the simulated values exhibit good agreement with the observed ones, as evident from the evenly dispersed data points around the equality line (Y = X) (Figure 9). These results indicate that the model effectively simulates nitrogen loss load under CID in the cold black soil region. The ε for NO3-N and NH4+-N loss load were 0.08 and 0.12, respectively. The NSE were 0.97 and 0.93, the R2 were 0.88 and 0.90, and the RMSEa were 0.04 and 0.29, respectively (Table 13). Overall, the simulated values align well with the observed values, indicating a satisfactory model performance in simulating nitrogen loss in paddy fields. However, the nitrogen loss load exhibited deviations at certain time periods for the following reasons. The accuracy of experimental observation needs to be improved. The drainage volume in the experimental area was observed manually, and the observation error of drainage volume was caused by the reason that led to the error of nitrogen loss. The accuracy of ammonium nitrogen simulation depends on the dispersion coefficient, nitrification reaction parameters and the appropriate decomposition temperature of organic matter, while the accuracy of nitrate nitrogen simulation depends on denitrification parameters and nitrification parameters [46]. Experimental means should be strengthened to ensure the accuracy of field measurements. Most of the parameters in the model calculation are determined by experiments, but some parameters are obtained indirectly by other means. The selection of parameters will affect the simulation results, so the experimental means should be expanded in the subsequent work to obtain as many measured parameters as possible. These aspects introduce uncertainties in the model simulation, and it is recommended to improve the accuracy of model simulation by improving the above aspects.

3.3. Simulation of Drainage Volume and Nitrogen Loss Load in Different Hydrological Years

3.3.1. Selection of Representative Hydrological Year

According to the rainfall data of 61 years from 1961 to 2021, the sum of the rainfall from May to September during the main growth period of rice in each year was calculated and the frequency was calculated. Based on the optimum line method of the Pilson-Ⅲ type, the rainfall of different years during the growth period was obtained, and the rainfall frequency of p ≤ 37.5%, 37.5% < p ≤ 62.5% and p > 62.5% were divided into the wet year group, normal year group and dry year group, respectively. According to the years corresponding to the rainfall frequency of 25%, 50% and 75%, the typical representative years of wet, normal and dry year are determined, and the typical representative years and rainfall of different rainfall frequencies are obtained (Table 14). The model was subsequently employed to simulate the drainage volume and nitrogen loss load in paddy fields in different representative hydrological years, and the variations in drainage volume and nitrogen loss load in different hydrological years were analyzed.

3.3.2. Simulation of Drainage and Nitrogen Loss Load in Different Hydrological Years

A comparison of simulation results from different hydrological years simulating paddy field drainage volume is depicted in Figure 10. Rainfall had a significant impact on paddy irrigation and drainage. Paddy field drainage varies significantly across different hydrological years. The greater the rainfall during the rice growth period, the smaller the irrigation quota and the greater the drainage volume. The paddy field drainage of CID mode in the normal representative year decreased by 74.7%, and the drainage times decreased by 50.0% compared with the wet years. Furthermore, no drainage occurred in the representative year of 75% hydrology. Because the paddy field makes full use of natural rainfall, the drainage volume and drainage times were reduced, so as to achieve the effect of controlling irrigation and emission reduction and saving manpower. During the rice growing period, the paddy field drainage is 57.4%, 82.9%, and 296.0% higher than that in the 5%, 25%, and 50% hydrological years, respectively, but 7.7% lower than that in the 10% hydrological year. This discrepancy is attributed to the more concentrated drainage-producing rainfall in the 10% hydrological year, resulting in higher drainage compared with the 1% hydrological year. Furthermore, no drainage occurred in the representative year of 75% hydrology, indicating that paddy field drainage is influenced not only by the total annual rainfall amount but also by the distribution pattern of rainfall throughout the year.

4. Discussion

4.1. The Simulation of Drainage Volume and Nitrogen Loss Load in the Paddy Field

It was highlighted that the model was calibrated and validated using two years of data from field experiments. The simulated paddy field drainage showed a high level of agreement with the observed results, as evidenced by the R2 above 0.88, the ε within 0.19, the RMSEa within 2.72, and the NSE above 0.72. These results demonstrated that the model was capable of accurately simulating the paddy field drainage under controlled irrigation and drainage in the cold black soil region. Some scholars applied the DRAINMOD-NII to the other irrigation areas and believed that it simulated the dynamic changes of paddy field drainage well, and these research results are basically consistent with those described in this study [47,48].
Simulation analysis of paddy drainage in different hydrological years revealed that the paddy field drainage volume was not only determined by annual precipitation, but was also significantly affected by precipitation distribution during the crop growth period. Specifically, the drainage volume tended to increase when the rainfall was more concentrated during the growth period. In addition, the risk of nitrogen loss in rice fields was higher under extreme climate conditions with 1%, 5% and 10% rainfall frequencies.
The model effectively simulated the nitrogen loss load of the paddy field. The simulated cumulative loss load of NO3-N and NH4+-N closely matched the observed values, as indicated by the R2 above 0.83, the ε within 0.80, the RMSEa within 1.09, and the NSE above 0.88. These findings substantiated the ability of the model to successfully simulate the nitrogen loss load under CID mode in the cold black soil region. Some scholars use the DRAINMOD-NII to explore the nitrogen changes during drainage [49]. Gao et al. [46] simulated the nitrogen loss load of paddy fields using DRAINMOD-NII and found that the model predicted the nitrogen loss characteristics of paddy fields effectively in typical irrigation areas in southern China, and the relative error ranged from 2.8% to 16.0%. Hashemi et al. [50] found that the DRAINMOD-NII simulated the NO3-N loss load well in paddy fields, and the relative error ranged from 1.0% to 35.0%. Morris [51] found that the simulated drainage and NH4+-N loss load in the whole growth period of rice were in good agreement with the measured values, with relative errors of 4.6% and 2.8%, respectively. However, in terms of the different growing stages of paddy rice, there was a large deviation between simulated and observed values. Preliminary analysis indicated that the parameters used in the model were not calibrated well due to limited observed data. On the whole, DRAINMOD-NII can efficiently describe and predict drainage and nitrogen loss processes in paddy fields, which can play an important role in improving the efficiency of water and nitrogen use in the paddy field as well as in reducing agricultural pollution. Those findings are consistent with the results of this study conducted in the cold black soil region.

4.2. Optimizing Water Management of Paddy Fields

The CID mode is crucial for optimizing resource management in paddy fields. The results showed that the CID mode reduced paddy field drainage compared to the TID mode, with decreases of 60.93% and 24.48% in 2020 and 2021, respectively. The loss load of NO3-N and NH4+-N in the CID paddy field decreased by 59.38% and 44.96%, respectively, compared with TID. Additionally, the yield of the CID mode surpassed that of the TID mode in different years, with increases of 14.15% and 19.62% in 2020 and 2021, respectively (Table 15). The rainfall utilization efficiency of a paddy field is determined by the ratio of the difference between rainfall and paddy field drainage to the total rainfall. Meanwhile, the irrigation water productivity is measured by the ratio of yield to irrigation, and the crop water productivity is evaluated using the ratio of yield to water consumption. In the CID paddy field, the rainfall utilization efficiency is 89.12% and 90.50% in 2020 and 2021, respectively. It is noteworthy that the CID mode exhibits a higher rainfall utilization efficiency compared to the TID mode, indicating an enhanced capability of the paddy field to store rainwater (Figure 11). Furthermore, the CID mode demonstrates higher crop water productivity compared to the TID mode, with increases of 17.04% and 22.57% in 2020 and 2021, respectively. The irrigation water productivity of the CID mode reached 162.65 and 77.12 kg·hm−2·mm−1 in 2020 and 2021, respectively, with the TID mode of 70.67 and 46.58 kg·hm−2·mm−1. The irrigation water productivity showed significant improvements, with increases of 130.13% and 65.56% in 2020 and 2021, respectively, compared to the TID mode. This indicated that the CID mode effectively reduced paddy field drainage and improved the irrigation water productivity of rice. Xie et al. [52] simulated the drainage volume and nitrogen loss load under different irrigation and drainage modes and found that the CID mode utilized nutrients and water efficiently and gave full play to the wetland function of paddy fields compared with the TID mode. The CID mode regulated the water status of the paddy field to achieve a high yield of rice, and the yield was increased by 10% compared with the TID mode; the difference was not significant [53]. The results are basically consistent with those of this experiment.
The CID mode proves to be an effective method of utilizing natural rainfall resources, and it not only improves rice yield and rainwater utilization efficiency but also reduces the drainage volume and nitrogen loss load. As a result, it achieves the goals of water saving, emission reduction, and high yield of the paddy field. During the main period of rice growth and development in the experimental area, from July to August, the CID mode effectively captured and utilized rainfall resources. Water saving was achieved by reducing paddy field drainage and improving rainfall utilization efficiency. Moreover, the CID mode optimally utilized the water storage function of the paddy field during flood seasons, reducing downstream flood control pressure [47]. Simultaneously, the CID mode minimized the loss of nitrogen in the paddy field and reduced the risk of non-point source pollution [48]. China’s rice planting area exceeds 30 million ha, which is widely distributed in the cold black soil region, where the growth period is synchronized with the rainy period. The CID mode improved the effective utilization efficiency of rainwater and reduced the irrigation quota, while the lower irrigation limit improved the water storage capacity of the paddy field and achieved the goal of water saving. Compared with TID, the irrigation amount of the paddy field in the CID mode decreased by 11.89~29.88% [48]. Therefore, the CID mode can effectively reduce runoff, which is of practical significance for reducing downstream drainage and flood control pressure.
However, there are challenges and considerations when implementing the CID mode, which should be adapted to suit local soil types, climate, and water availability. Tailoring practices to specific environmental conditions ensures the effectiveness of water resource optimization strategies. In conclusion, controlled irrigation and drainage play a pivotal role in optimizing water resources in paddy fields, contributing to sustainable and efficient rice cultivation. By implementing efficient and sustainable practices, it is possible to conserve water, improve nutrient utilization, and ultimately enhance the overall productivity and sustainability of rice cultivation.

5. Conclusions

These results demonstrate that the model effectively simulated drainage volume and nitrogen loss load in paddy fields under controlled irrigation and drainage in the cold black soil region; the coefficient of determination between the simulated paddy field drainage volume and nitrogen loss load results and the observed values were all above 0.83, and the Nash–Sutcliffe efficiency coefficient ranged from 0.72 to 0.97 in model calibration and verification. The CID mode effectively reduced both irrigation amount and frequency, enhanced the utilization of natural rainfall resources, significantly reduced the irrigation amount by 39.07% and increased rainwater utilization efficiency by 13.07% compared with the TID mode. The CID paddy field decreased the drainage volume by 44.71% and NO3-N and NH4+-N loss load by 59.38% and 44.96%, respectively, compared with TID mode. The CID mode can optimize natural rainfall resources and increase irrigation water productivity by 97.85%. At the same time, the yield of rice increased by 16.88% compared with the TID mode, thereby achieving the goal of water saving, emission reduction, and high-yield goals in paddy fields. The results of this study can serve as a valuable reference for rice production under similar environmental conditions in other regions. Field studies were conducted for only two years on a sandy loam with a heavy texture and good structural properties. It is necessary to test and promote research in subsequent years and on the other types of soil, different field scales, different hydrological years and geographical divisions to improve the model.

Author Contributions

Conceptualization, Y.M. and Y.Y.; data curation, Y.M.; investigation, Y.M. and S.Y.; methodology, Y.M. and S.Y. resources, Y.M. and Y.Y.; software, Y.M.; writing original draft preparation, Y.M.; writing-review and editing, Y.Y., P.C. and K.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, grant: 51879074; the Research Business Expenses of Heilongjiang Provincial Research Institutes: CZKYF2023-1-C031; Natural Science Foundation of Heilongjiang Province: LH2019E122; Postdoctoral Fellowship Program of CPSF: GZC20230668.

Data Availability Statement

Data supporting the findings of this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Tan, X.Z.; Shao, D.G.; Gu, W.Q.; Liu, H.H. Field analysis of water and nitrogen fate in lowland paddy fields under different water managements using HYDRUS-1D. Agric. Water Manag. 2015, 150, 67–80. [Google Scholar] [CrossRef]
  2. Yang, S.H.; Peng, S.Z.; Xu, J.Z.; Hou, H.J.; Gao, X.L. Nitrogen loss from paddy fields with different water and nitrogen managements in Taihu Lake region of China. Commun. Soil Sci. Plant Anal. 2013, 44, 2393–2407. [Google Scholar] [CrossRef]
  3. Liu, X.F.; Fei, L.J.; Meng, Z.J.; Zhang, J.Y.; Liu, X.L.; Zhang, Y.Y. Effects of water and nutrient synergies on dry matter translocation and nitrogen uptake and utilization in winter wheat. J. Plant Nutr. Fertil. 2018, 24, 905–914. [Google Scholar]
  4. Ye, Y.S.; Liang, X.Q.; Chen, Y.X.; Liu, J.; Gu, J.T.; Guo, R.; Li, L. Alternate wetting and drying irrigation and controlled-release nitrogen fertilizer in late-season rice. Effects on dry matter accumulation, yield, water and nitrogen use. Field Crops Res. 2013, 144, 212–224. [Google Scholar] [CrossRef]
  5. Yang, J.Z.; Wang, L.Y.; LU, X.Y.; Jayawardane, N. Experiment and numerical simulation of nitrogen transport in soils irrigated with treated sewage. Irrig. Drain. 2008, 57, 203–217. [Google Scholar] [CrossRef]
  6. Antonopoulos, V.Z. Modelling of water and nitrogen balances in the ponded water and soil profile of paddy fields in Northern Greece. Agric. Water Manag. 2010, 98, 321–330. [Google Scholar] [CrossRef]
  7. Tan, X.Z.; Shao, D.G.; Liu, H.H.; Sun, C.M. Experimental and simulated water-nitrogen balance in rice paddy fields under water-saving irrigation with controlled drainage. Trans. Chin. Soc. Agric. Eng. 2011, 27, 193–198. [Google Scholar]
  8. Ferm, M. Atmospheric ammonia and ammonium transport in Europe and critical loads, a review. Nutr. Cycl. Agroecosystems 1998, 51, 5–17. [Google Scholar] [CrossRef]
  9. Huang, S.F.; Shen, G.X.; Tang, H.; Lu, Y.T. Study on nitrogen loss in paddy fields in suburban Shanghai. Environ. Pollut. Control 2005, 27, 651–654. [Google Scholar]
  10. Wesstrom, I.; Joel, A.; Messing, I. Controlled drainage and subirrigation—A water management option to reduce non-point source pollution from agricultural land. Agric. Ecosyst. Environ. 2014, 198, 74–82. [Google Scholar] [CrossRef]
  11. Wei, Q.; Xu, J.Z.; Liu, Y.Z.; Wang, D.; Chen, S.Y.; Qian, W.H.; He, M.; Chen, P.; Zhou, X.Y.; Qi, Z.M. Nitrogen losses from soil as affected by water and fertilizer management under drip irrigation: Development, hotspots and future perspectives. Agric. Water Manag. 2024, 296, 108791. [Google Scholar] [CrossRef]
  12. Youssef, M.A.; Abdelbaki, A.M.; Negm, L.M.; Skaggs, R.W.; Thorp, K.R.; Jaynes, D.B. DRAINMOD-simulated performance of controlled drainage across the US Midwest. Agric. Water Manag. 2018, 197, 54–66. [Google Scholar] [CrossRef]
  13. He, X.C.; Shao, D.G.; Liu, W.Y.; Dai, T. Research progress and prospect of resource utilization of agricultural drainage. Trans. Chin. Soc. Agric. Eng. 2006, 22, 176–179. [Google Scholar]
  14. Lu, W.X.; Cheng, W.G.; Zhang, Z.; Xin, X.; Wang, X.H. Differences in rice water consumption and yield under four irrigation schedules in central Jilin Province, China. Paddy Water Environ. 2016, 14, 473–480. [Google Scholar] [CrossRef]
  15. Xiao, M.H.; Yu, S.E.; She, D.L.; Hu, X.J.; Chu, L.L. Nitrogen and phosphorus loss and optimal drainage time of paddy field under controlled drainage condition. Arab. J. Geosci. 2015, 8, 4411–4420. [Google Scholar] [CrossRef]
  16. Li, J.; Shao, M.N.; Zhang, X.C. Mathematical simulation of soil nitrogen and phosphorus transport and crop nutrition in EPIC model. J. Plant Nutr. Fertil. 2005, 11, 166–173. [Google Scholar]
  17. Beasley, D.; Huggins, L.; Monke, E. ANSWERS: A model for watershed planning. Trans. ASAE 1980, 23, 938–944. [Google Scholar] [CrossRef]
  18. Chung, S.; Kim, H.; Kim, J. Model development for nutrient loading from paddy rice fields. Agric. Water Manag. 2003, 62, 1–17. [Google Scholar] [CrossRef]
  19. Hao, F.H.; Sun, W.; Zeng, A.Y.; Li, P.; Zhang, J.X.; Yue, Y. Simulation of nitrogen migration by HYDRUS-1D model under different irrigation scenarios in Hetao irrigation district. J. Environ. Sci. 2008, 28, 853–858. [Google Scholar]
  20. Li, Y.; Simunek, J.; Jin, L.F.; Zhang, Z.T.; Ni, L.X. Evaluation of nitrogen balance in a direct-seeded-rice field experiment using Hydrus-1D. Agric. Water Manag. 2015, 148, 213–222. [Google Scholar] [CrossRef]
  21. Li, S. Simulation of Non-Point Source Nitrogen and Phosphorus Pollution and Response of Lake Sediments in Nansi Lake Basin Based on SWAT Model. Ph.D. Thesis, Shandong Normal University, Jinan, China, 2012. [Google Scholar]
  22. Hoang, L.; Griensven, A.; Mynett, A. Enhancing the SWAT model for simulating denitrification in riparian zones at the river basin scale. Environ. Model. Softw. 2017, 93, 163–179. [Google Scholar] [CrossRef]
  23. Li, H.; Qiu, J.J.; Gao, C.Y.; Wang, L.G. Estimation of nitrogen leaching potential of farmland in a typical basin around Bohai Sea based on DNDC model. Trans. Chin. Soc. Agric. Eng. 2012, 28, 127–134. [Google Scholar]
  24. Moursi, H.; Youssef, M.A.; Chescheir, G.M. Development and application of DRAINMOD model for simulating crop yield and water conservation benefits of drainage water recycling. Agric. Water Manag. 2022, 266, 107592. [Google Scholar] [CrossRef]
  25. Skaggs, R.W. A Water Management Model for Shallow Water Table Soils; Water Resources Research Institute of the University of North Carolina: Chapel Hill, NC, USA, 1978. [Google Scholar]
  26. Breve, M.A. Modeling the Movement and Fate of Nitrogen in Artificially Drained Soils; North Carolina State University: Raleigh, NC, USA, 1994. [Google Scholar]
  27. Youssef, M.A.; Skaggs, R.W. The nitrogen simulation model, DRAINMOD-NII: Field testing and model. Am. Soc. Civ. Eng. 2006, 1–11. [Google Scholar] [CrossRef]
  28. Hong, L.; Luo, W.B. Dynamic simulation of nitrogen loss in surface runoff from farmland based on DRAINMOD. Adv. Water Sci. 2011, 22, 703–709. [Google Scholar]
  29. Qi, Z.M.; Du, X.; Feng, H.; Matthew, J.H. Influence of winter long-term covering rye on subsurface drainage and nitrate nitrogen loss based on DRAINMOD-NII. Trans. Chin. Soc. Agric. Eng. 2017, 33, 153–161. [Google Scholar]
  30. Salazar, O.; Wesstrom, I.; Youssef, M.A.; Skaggs, R.W.; Joel, A. Evaluation of the DRAINMOD-N II model for predicting nitrogen losses in a loamy sand under cultivation in south-east Sweden. Agric. Water Manag. 2009, 96, 267–281. [Google Scholar] [CrossRef]
  31. Li, S.; Luo, W.; Jia, Z.H.; Pan, Y.X.; Wu, D.; Zhang, D.K. Estimation of shallow groundwater withdrawal and salt accumulation in irrigated area based on DRAINMOD model. Trans. Chin. Soc. Agric. Eng. 2015, 31, 89–97. [Google Scholar]
  32. Jia, Z.H.; Luo, W.; Mo, F.; Cheng, H.Y. Effects of drainage and inflow on wetland hydrology using DRAINMOD model under different climatic conditions. J. Soil Water Conserv. 2003, 5, 54–58. [Google Scholar]
  33. Luo, W.; Jia, Z.; Skaggs, R.W.; Xi, W.Y.; Zhang, Y.F. Modeling of drainage process in yinnan irrigation district with DRAINMOD. Trans. Chin. Soc. Agric. Eng. 2006, 22, 53–57. [Google Scholar]
  34. Luo, W.; Gary, R.; Jing, W.H.; Jia, Z.H. Hydrological effects of planting perennial grasses in the "corn belt" based on simulation with DRAINMOD. Trans. Chin. Soc. Agric. Eng. 2010, 26, 89–94. [Google Scholar]
  35. Wang, S.L.; Wang, X.K.; Prasher, S.O.; Tan, C.S. Simulation of groundwater table and surface drainage using DRAINMOD model. Trans. Chin. Soc. Agric. Eng. 2006, 22, 54–59. [Google Scholar]
  36. Zhang, Z.; Zhang, Y.; Zhang, J.; Ma, H.Y.; Zhu, W.Y. Simulation of subsurface drainage in coastal saline-alkali farmland based on DRAINMOD-S Model. Adv. Water Sci. 2012, 23, 782–788. [Google Scholar]
  37. Abduljaleel, Y.; Awad, A.; Al-Ansari, N.; Salem, A.; Negm, A.; Gabr, M.E. Assessment of subsurface drainage strategies using DRAINMOD model for sustainable agriculture: A review. Sustainability 2023, 15, 1355. [Google Scholar] [CrossRef]
  38. Gao, Y.; Xu, X.; Sun, C.; Ding, S.B.; Huo, Z.L.; Huang, G.H. Parameterization and modeling of paddy rice growth and water use in cold regions, yield and water-saving analysis. Agric. Water Manag. 2021, 250, 106864. [Google Scholar] [CrossRef]
  39. Hansen, J.; Fuller, F.; Gale, F.; Crook, F.; Wailes, E.; Moore, M. China’s japonica rice market, growth and competitiveness. In Rice Situation and Outlook Yearbook; Economic Research Service: Washington, DC, USA, 2002; pp. 32–37. [Google Scholar]
  40. Zade, A. Technologies of quality, the role of the Chinese state in guiding the market for rice. East Asian Sci. Technol. 2020, 5, 461–477. [Google Scholar] [CrossRef]
  41. Liu, S.; Chen, X.F.; Liu, L. Evapotranspiration algorithm and application based on the Penman-Monteith equation. Heilongjiang Sci. Technol. Inf. 2013, 40–41. [Google Scholar]
  42. The State Environmental Protection Administration. Water and Wastewater Monitoring and Analysis Methods, 4th ed.; Environmental Science Press: Beijing, China, 2002. [Google Scholar]
  43. SL 13-2015; Water Industry Standard of the People’s Republic of China, Specifications for Irrigation Experiment. Ministry of Water Resources of Pelple’s Repubic of China: Beijing, China, 2015.
  44. Bouman, B.A.M.; Laar, H.H. Description and evaluation of the rice growth model ORYZA2000 under nitrogen-limited conditions. Agric. Syst. 2006, 87, 249–273. [Google Scholar] [CrossRef]
  45. Nash, J.; Sutcliffe, J.V. River flow forecasting through conceptual models part I—A discussion of principles. J. Hydrol. 1970, 10, 282–290. [Google Scholar] [CrossRef]
  46. Gao, X.R.; Dong, B.; Qin, D.Y.; Sun, N.N.; Zhang, Z.Y. Simulation of drainage process and nitrogen loss in paddy field using DrainMOD. Trans. Chin. Soc. Agric. Eng. 2011, 27, 52–58. [Google Scholar]
  47. Wu, D.; Luo, W.; Zhou, M.Y.; Jia, Z.H.; Li, S.; Li, L. Paddy field drainage processes and its influence by lateral ditches in thepolder area of southeastern China. J. Irrig. Drain. 2015, 34, 28–33. [Google Scholar]
  48. Golmohammadi, G.; Rudra, R.P.; Prasher, S.O.; Madani, A.; Goel, P.K.; Mohammadi, K. Modeling the impacts of tillage practices on water table depth, drain outflow and nitrogen losses using DRAINMOD. Comput. Electron. Agric. 2016, 124, 73–83. [Google Scholar] [CrossRef]
  49. Wang, J.; Yu, S.E.; Wang, M.; Zhang, Y.D.; Gao, S.K. Effect of water depth over paddy field on nitrogen and phosphorus in the water after storms at tillering stage. J. Irrig. Drain. 2018, 37, 71–75, 82. [Google Scholar]
  50. Hashemi, S.Z.; Darzi-Naftchali, A.; Qi, Z. Assessing water and nitrate-N losses from subsurface-drained paddy lands by DRAINMOD-NII. Irrig. Drain. 2020, 69, 776–787. [Google Scholar] [CrossRef]
  51. Morris, M.D. Factorial sampling plans for preliminary computational experiments. Technometrics 1991, 33, 161–174. [Google Scholar] [CrossRef]
  52. Xie, Y.C.; Xu, M.; Gao, S.K.; Shen, L.Q. Simulation of water and nitrogen transport in rice field with different irrigation and drainage models based on DRAINMOD model. J. Irrig. Drain. 2021, 40, 37–44. [Google Scholar]
  53. Shao, G.C.; Cui, J.T.; Yu, S.E.; Lu, B.; Brian, B.J.; Ding, J.H.; She, D.L. Impacts of controlled irrigation and drainage on the yield and physiological attributes of rice. Agric. Water Manag. 2015, 149, 156–165. [Google Scholar] [CrossRef]
Figure 1. Soil water characteristic curve of field plot.
Figure 1. Soil water characteristic curve of field plot.
Agronomy 14 01095 g001
Figure 2. Drainage simulation of paddy field in 2020 (model calibration). (SIM: simulated value; OBS: observed value).
Figure 2. Drainage simulation of paddy field in 2020 (model calibration). (SIM: simulated value; OBS: observed value).
Agronomy 14 01095 g002
Figure 3. Comparison of observed and simulated values of drainage in paddy field in 2020 (model calibration).
Figure 3. Comparison of observed and simulated values of drainage in paddy field in 2020 (model calibration).
Agronomy 14 01095 g003
Figure 4. Simulation of NH4+-N and NO3-N loss load in paddy field in 2020 (model calibration).
Figure 4. Simulation of NH4+-N and NO3-N loss load in paddy field in 2020 (model calibration).
Agronomy 14 01095 g004
Figure 5. The observed and simulated values of NO3-N and NH4+-N loss load in paddy field in 2020 (model calibration).
Figure 5. The observed and simulated values of NO3-N and NH4+-N loss load in paddy field in 2020 (model calibration).
Agronomy 14 01095 g005
Figure 6. Drainage simulation of paddy field in 2021 (model verification).
Figure 6. Drainage simulation of paddy field in 2021 (model verification).
Agronomy 14 01095 g006
Figure 7. Comparison of observed and simulated values of drainage in 2021 (model verification).
Figure 7. Comparison of observed and simulated values of drainage in 2021 (model verification).
Agronomy 14 01095 g007
Figure 8. Simulation of NH4+-N and NO3-N loss in paddy field in 2021 (model verification).
Figure 8. Simulation of NH4+-N and NO3-N loss in paddy field in 2021 (model verification).
Agronomy 14 01095 g008
Figure 9. The observed and simulated values of NO3-N and NH4+-N loss in 2021 (model verification).
Figure 9. The observed and simulated values of NO3-N and NH4+-N loss in 2021 (model verification).
Agronomy 14 01095 g009
Figure 10. Simulation of paddy field drainage in different hydrological years.
Figure 10. Simulation of paddy field drainage in different hydrological years.
Agronomy 14 01095 g010aAgronomy 14 01095 g010b
Figure 11. Rainfall utilization efficiency and water productivity under different modes.
Figure 11. Rainfall utilization efficiency and water productivity under different modes.
Agronomy 14 01095 g011
Table 1. Test design.
Table 1. Test design.
ModeGrowth PeriodWater Depth of Paddy Field
CIDReplanting
Tillering
Jointing and booting
Heading and flowering
Milking ripening
Yellow ripening
1~3~5
−10~3~5
−20~3~10
−20~3~10
−20~3~10
Dry naturally after drainage
TIDReplanting~Milking ripening
Yellow ripening
1~5~5
Dry naturally after drainage
Notes: The regulation value of the water depth in the paddy field was based on the field surface. The water depth refers to the vertical distance from the paddy field water level to the field surface. For example, in −20~3~10, the left value represents the lower limit of paddy water depth, the middle value represents the upper limit of paddy water depth in the paddy field, and the right value represents the allowable maximum water depth after rainfall.
Table 2. Agronomic measures and time of rice cultivation.
Table 2. Agronomic measures and time of rice cultivation.
YearRaise SeedlingsSteeping FieldTransplantingBase FertilizerTillering FertilizerEar FertilizerHarvest
20204.55.85.155.125.267.59.25
20214.15.15.185.125.287.189.28
Table 3. Kc of rice.
Table 3. Kc of rice.
MonthMayJuneJulyAugustSeptember
Kc0.380.781.341.060.45
Table 4. Soil physical properties and hydraulic parameters of different soil depths.
Table 4. Soil physical properties and hydraulic parameters of different soil depths.
Soil Depth
/cm
Distribution of Soil Particle Size/%Bulk Density
/g·cm−3
Field Capacity
/cm3·cm−3
Wilting Coefficient
/cm3·cm−3
Saturated Moisture Content
/cm3·cm−3
Lateral Saturation Conductivity
/cm·d−1
2.0~0.05
/mm
0.05~0.002
/mm
<0.002
/mm
0~2020.067.212.81.470.290.0930.491.40
20~409.176.914.01.480.280.0680.431.12
40~6011.374.913.81.530.260.0580.440.95
60~8017.871.710.51.290.230.0510.410.75
Table 5. Main input parameters of drainage design parameters.
Table 5. Main input parameters of drainage design parameters.
Relative Impervious Depth
/cm
Maximum Surface Water Storage Depth (Sm)
/cm
Drainage Coefficient
/cm·d−1
Drainage Depth
/cm
Drainage Spacing
/cm
Krikham Water Depth (SI)
/cm
28010.03.05010000.8
Table 6. Main input parameters of nitrogen parameters.
Table 6. Main input parameters of nitrogen parameters.
Maximum Nitration Rate
/µg·ng−1·d−1
Optimum Nitration Temperature
/°C
Maximum Denitrification Rate
/µg·ng−1·d−1
Optimum Denitrification Temperature
/°C
Nitrogen Dissolution Rate
/d−1
NO3-N Concentration
in Rain
/mg·L−1
NH4+-N Concentration
in Rain
/mg·L−1
25253.5201.51.00.75
Table 7. Drainage and evaporation capacity of soil.
Table 7. Drainage and evaporation capacity of soil.
Groundwater depth/cm15253575120150200500
Drainage/cm0.070.220.401.804.256.8511.5055.10
Submersible rising flux/cm0.300.280.250.190.180.150.100.05
Table 8. Green–Ampt infiltration parameters of the model.
Table 8. Green–Ampt infiltration parameters of the model.
Water table/cm010406080100200
A/cm·h−10.850.150.721.201.651.903.60
B/cm·h−10.900.900.900.900.900.900.90
Table 9. Saturated water conductivity.
Table 9. Saturated water conductivity.
Soil depth/cm0~2020~4040~60
Vertical saturated hydraulic conductivity/cm·h−11.740.920.78
Horizontal saturated hydraulic conductivity/cm·h−11.251.020.45
Table 10. Evaluation value of paddy field drainage volume in 2020 (model calibration).
Table 10. Evaluation value of paddy field drainage volume in 2020 (model calibration).
Evaluation ParameterRMSEaεNSER2
Drainage volume2.080.140.940.88
Table 11. Evaluation value of nitrogen loss load in 2020 (model calibration).
Table 11. Evaluation value of nitrogen loss load in 2020 (model calibration).
Evaluation ParameterRMSEaεNSER2
NO3-N loss load0.220.230.880.83
NH4+-N loss load1.090.800.920.85
Table 12. Evaluation value of paddy field drainage volume in 2021 (model verification).
Table 12. Evaluation value of paddy field drainage volume in 2021 (model verification).
Evaluation ParameterRMSEaεNSER2
Drainage volume2.720.190.720.90
Table 13. Evaluation value of nitrogen loss load in paddy field in 2021 (model verification).
Table 13. Evaluation value of nitrogen loss load in paddy field in 2021 (model verification).
Evaluation ParameterRMSEaεNSER2
NO3-N loss load0.040.080.970.88
NH4+-N loss load0.290.120.930.90
Table 14. Rainfall of different typical years in the test area.
Table 14. Rainfall of different typical years in the test area.
Rainfall Frequency
/%
Rainfall
/mm
Hydrological YearRepresentative Year
1742.4Wet year2012
5700.4Wet year1985
10633.3Wet year1988
25551.3Wet year1998
50483.5Normal year1966
75416.4Dry year2006
Notes: % is the average probability of one occurrence per hundred years.
Table 15. Evaluation and comparison of different irrigation and drainage modes.
Table 15. Evaluation and comparison of different irrigation and drainage modes.
ModeYearWater SavingEmission ReductionHigh Yield
Irrigation Amount
/mm
Rainwater Utilization Rate
/%
Drainage
/mm
NH4+-N
Loss Load
/kg·hm−2
NO3-N
Loss Load
/kg·hm−2
Yield
/kg·hm−2
CID202074.589.1254.03.60.612,117.15
2021140.690.5059.03.50.710,843.51
TID2020150.275.64138.28.52.110,615.35
2021194.683.6582.54.41.19065.12
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Meng, Y.; Yu, Y.; Yu, S.; Chen, P.; Cao, K. Simulation of Drainage Volume and Nitrogen Loss Load in Paddy Fields under Different Irrigation and Drainage Modes and Hydrological Years. Agronomy 2024, 14, 1095. https://doi.org/10.3390/agronomy14061095

AMA Style

Meng Y, Yu Y, Yu S, Chen P, Cao K. Simulation of Drainage Volume and Nitrogen Loss Load in Paddy Fields under Different Irrigation and Drainage Modes and Hydrological Years. Agronomy. 2024; 14(6):1095. https://doi.org/10.3390/agronomy14061095

Chicago/Turabian Style

Meng, Yan, Yanmei Yu, Shuangen Yu, Peng Chen, and Kaihua Cao. 2024. "Simulation of Drainage Volume and Nitrogen Loss Load in Paddy Fields under Different Irrigation and Drainage Modes and Hydrological Years" Agronomy 14, no. 6: 1095. https://doi.org/10.3390/agronomy14061095

APA Style

Meng, Y., Yu, Y., Yu, S., Chen, P., & Cao, K. (2024). Simulation of Drainage Volume and Nitrogen Loss Load in Paddy Fields under Different Irrigation and Drainage Modes and Hydrological Years. Agronomy, 14(6), 1095. https://doi.org/10.3390/agronomy14061095

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