Next Article in Journal
Manual but Not Instrument-Assisted Cervical Manipulation Reduces Pain and Disability in Subjects with Nonspecific Neck Pain: Double-Blinded, Randomized Clinical Trial
Previous Article in Journal
Influence of Moisture in Concrete on the Photothermal Response—A New Approach for a Measurement Method?
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Appropriate Method to Estimate Farmland Drainage Coefficient in the Wanyan River Surface Waterlogged Area in Suibin County of the Sanjiang Plain, China

1
State Key Laboratory of Simulation and Regulation of Water Cycle in River Basin, China Institute of Water Resource and Hydropower Research, Beijing 100038, China
2
Department of Irrigation and Drainage, China Institute of Water Resources and Hydropower Research, Beijing 100038, China
3
Heilongjiang Province Service Center for Rural Water Conservancy and Hydropower, Harbin 150040, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2023, 13(5), 2769; https://doi.org/10.3390/app13052769
Submission received: 5 January 2023 / Revised: 17 February 2023 / Accepted: 19 February 2023 / Published: 21 February 2023

Abstract

:
Drainage coefficient has an important influence on the determination of the scale of drainage engineering and agricultural production. This paper calculated the drainage coefficient based on the calculation results of theoretical runoff using the empirical formula (QE) and the average draining method (QA), and discussed the drainage duration (TE) corresponding to the empirical formula and the area range (FA) corresponding to the average draining formula. The results showed that the observed runoff was less than the theoretical runoff. The design flow of the Wanyan River pumping station does not meet the drainage demand of the 5-year rainfall return period. For the 1-day rainfall with 3-year return period and 3-day rainfall, QE was higher than QA, and the design drainage duration was higher than TE, so it was safer to calculate the drainage coefficient using the empirical formula in the drainage engineering design. For 1-day rainfall greater than the 3-year return period, QA was higher than QE, and the design drainage duration was less than TE, and the results of the average draining formula were safer in the drainage engineering design. FA increased with the increase in rainfall duration and decreased with the increase in return period. When the rainfall duration and the drainage duration were long, the calculation results of the average draining formula were much smaller.

1. Introduction

IPCC’s Assessment Report 6 (AR6) noted that human activities have contributed to an increase of 1.09 °C in global temperature since industrialization [1]. A global concern has been raised due to the climate change issues caused by temperature rise. In addition to affecting rainfall, runoff, infiltration, evaporation, and other components of the water cycle [2], climate change will also change the distribution of regional water resources [3] and increase the frequency and scale of rainstorms and floods [4,5]. Lots of researches have demonstrated that there have been more global floods than ever before since 2000 [6], which has aggravated farmland waterlogging and adversely affected human production activities and food security [7].
One of the most common agricultural disasters is farmland waterlogging caused by excessive rain. On the one hand, farmland waterlogging affects crop growth, lowers the activity of the soil’s microbial population [8], and influences the productivity of the land [9,10]. It is reported that farmland waterlogging decreases crop yield by 32.9% on average [11]. On the other hand, farmland waterlogging causes pesticides and fertilizers to enter the water cycle, which increases the risk of water pollution and salinization [12], causes agricultural non-point source pollution, and damages the agricultural environment [13].
Sustainable crop production depends on effective farmland drainage, which is also a key solution to farmland waterlogging. A reasonable scale of drainage system can ensure crop growth and food security. Therefore, it is very important to determine the design parameters of the surface drainage system. Usually, a surface drainage system is designed based on the removal rate of excess surface water, which is expressed as the water depth drained off from a given area in a day or the design discharge per unit area, termed as the surface drainage coefficient [14]. The drainage coefficient is the main index for determining the scale and capacity of the surface drainage system. The drainage coefficient for any region varies with the geographical locations; land use; size of the area; rainfall intensity, frequency, and duration; and other climatic factors. The key to determine the drainage coefficient is the estimation of surface runoff.
At present, the commonly used surface runoff calculation methods include the runoff coefficient method, the infiltration curve method, the excess water storage method, the SCS-CN, etc. Internationally, the runoff coefficient method was proposed in the 1930s [15], and was often used to calculate the runoff of urban impervious surfaces where evaporation and infiltration are negligible. China began to study runoff coefficient in the 1950s and has established a complete set of runoff coefficient values so far [16]. The principle of the infiltration curve model is that the rainwater after deducting evaporation, vegetation interception, infiltration, and groundwater storage in the basin is surface runoff. Infiltration calculation includes the Horton formula [17], Green-Ampt formula [18] and the Philip formula [19]. In the 1960s, China first proposed the excess water storage model, called the Xin’anjiang model [20], assuming the rainwater would be fully converted into runoff when the soil water content in the aeration zone reached the field water holding capacity. The Soil Conservation Service-Curve Number (SCS-CN) method, developed by the National Resources Conservation Service (NRSC), United States Department of Agriculture in 1969, was used for estimating direct runoff depth based on rainstorm depth [21]. This method was relatively easy to use, simple in structure, suitable for small catchment, and has been widely used in runoff estimation all over the world. Italy, Spain, Greece, and France are countries that widely use the continuation method and the movement method to calculate surface runoff [22,23]. In a word, the water balance principle is the essence of all runoff calculation methods [24]. In a practical application, the impact of climate, topography, underlying surfaces, and human activities on runoff should be comprehensively considered and the most suitable method should be chosen to calculate surface runoff.
Some researchers have established the relationship among the drainage coefficient, catchment area and surface runoff, such as the Cypress-Creak formula in the United States and the empirical formula in China and the former Soviet Union [25,26,27]. These formulas are suitable for the drainage design of river channels and ditches in plain areas with large watershed areas. The Rational Formula [28] was proposed on the base of the empirical relationship between the peak runoff rate and rainfall intensity, runoff coefficient, as well as catchment area, which is mainly suitable for small catchments in hilly areas. The average draining formula, which is a suitable method for the drainage design of small plain basins, is to drain off the design runoff in the drainage area within the designed drainage duration, and the drainage coefficient can be obtained. Paul et al. [29] analyzed the daily rainfall data of 21 years and calculated the drainage coefficient by using the simplified hydrological method, SCS-CN method, Cypress-Creak formula, and Rational formula. The results not only obtained the drainage flow of the surface drainage system, but also proved that the Rational formula was suitable for the design of river channels and outfalls. Luo et al. [26] computed the drainage coefficient in the Luoshan drainage area by using the average draining formula and the empirical formula. The results showed that the drainage coefficient calculated by the average draining formula was greater than that calculated by the empirical formula for 1 day rainfall duration and was less than that calculated by the empirical formula for 3 day rainfall duration under the same return period. Wang et al. [27] calculated and compared the drainage coefficient of the Huaibei plain calculated by the empirical formula and the average draining formula and adjusted the peak value relation index (m) in the empirical formula. It was found that if the drainage duration in the average draining formula was set as 1 day and m in the empirical formula was set as one, the calculation results of the two formulas were close and the adjusted m range was taken as 0.84–0.87, meaning that it was more reasonable to calculate the drainage flow of the Huaibei Plain using the empirical formula. In addition, Indian scholars calculated the drainage coefficient by subtracting the basic infiltration rate from the consecutive days maximum rainfall. Dabral et al. calculated the drainage coefficients of North Lakhimpur (Assam) [30] and Doimukh (Arunachal Pradesh), India, [14] for different return periods in 2008 and 2016, respectively. The former was primarily made up of loam with an infiltration rate ranging from 1 to 5 mm/h, whereas the latter was primarily made up of sandy soil with a basic infiltration rate of 10 mm/h. The results showed that the drainage coefficient with a 25-year return period was the largest, while the drainage coefficient with a 5-year return period was very small, and certain drainage measures were required for both regions. At the same time, the author did not consider evapotranspiration, surface retention, and raindrop interception.
The average draining formula or the empirical formula is often used in China, but both methods have some drawbacks. The parameters in the empirical formula were determined based on hydrological observation data in the 1970s; whether the parameters meet the requirements of the current conditions needs further study. The maximum drainage flow in the drainage area is the design basis of the empirical formula, which considers the influence of runoff, area, and shape of the waterlogging area on the drainage coefficients, but does not consider the influence of drainage duration. The average draining formula is applied well in small catchments and the calculation results are small in plain areas with an area of more than 10 square kilometers. However, in the existing studies, the range of area in which this method can be applied is from 0 to 150 km2 [26,27], and the applicable area of the method needs to be further studied.
The Wanyan River waterlogged area in Suibin County of the Sanjiang Plain, China, was chosen as the study area in this paper, and the theoretical runoff and observed runoff were calculated based on rainfall data, pumping station data, and land use type data over the years. The drainage coefficient was calculated by theoretical runoff. This paper analyzed the drainage duration corresponding to the drainage coefficient estimated by the empirical formula and the area range corresponding to the average draining formula.

2. Materials and Methods

2.1. Study Area

Suibin County (131°08′121″–132°31′00″ E, 47°11′55″–47°45′23″ N) is located in the northeast of Hegang City, Heilongjiang Province, China, in the confluence zone of the Heilongjiang River and Songhua River, covering an area of 3344 km2. It is an important part of the Sanjiang Plain. There is the Heilongjiang River between Suibin County and Russia in the north, Luobei County in the west, and Huachuan County and Fujin City in the south across the Songhua River. Suibin County belongs to the mid-temperate continental monsoon climate zone, with a large temperature difference between day and night in spring and autumn, a hot and short summer, abundant rainfall, and severe cold and a long winter. The average annual rainfall is 567.9 mm, which is unevenly distributed throughout the year, mostly concentrated from July to September. The average annual evaporation is 650.8 mm, which is the largest in spring and easily causes spring drought. The average annual temperature is 1.8–2.4 °C, with a maximum of 21.9 °C in July and a minimum of minus 19.6 °C in January.
There are 25 rivers in Suibin County, including 5 rivers in the Heilongjiang River Basin and 20 rivers in the Songhua River Basin. Due to the uneven distribution of annual rainfall, the distribution of annual runoff is also uneven, and the runoff from July to September accounts for more than half of the annual runoff. The terrain in the area is flat, with an average gradient of 1/6000–1/10,000, and the altitude range is between 55 and 77 m. The landforms belong to terraces, high floodplains, and low floodplains. Most of the soil in Suibin County is loamy soil with strong water permeability, and the infiltrability is 0.144–0.47 m/d. Because there is an aquifer layer at a depth of 50~100 cm below the surface, the upper soil has little water retention capacity, and the soil is vulnerable to drought and waterlogging. At present, Suibin County has become one of the largest grain producing counties in Heilongjiang Province, mainly planting wheat, soybeans, rice, and corn.
The Wanyan River surface waterlogged area, namely the study area, is located in the central part of Suibin County and is one of the major waterlogged areas in the Sanjiang Plain. It is bounded by the watershed between the Heilongjiang River basin and the Wanyan River basin in the north; the watershed between the Songhua River basin, the Wanyan River basin, and the Songhua River dam in the south; and the No. 1 road of Suibin Farm in the west and the Songhua River dam in the east, as shown in Figure 1. The study area is 76.88 km long and 16 km wide with a total area of 1230 km2, including 528.8 km2 of Suibin County, 263.67 km2 of Suibin Farm, and 437.53 km2 of 290 Farm.
The Wanyan River is the main drainage ditch with a length of 92.40 km. The Wanyan River pumping station is located at the end outlet, which was built in 1982 and renovated in 2011, with a design flow of 57.97 m3/s. There are also 18 sub-main ditches with a total length of 163.89 km.

2.2. Research Data and Analysis

The daily rainfall of 63 years from 1957 to 2019 was collected, the annual maximum value method and the Pearson type three distribution (P-Ⅲ) [31,32] were used to obtain 1-day and 3-day design maximum rainfall [33]. The results showed that the maximum 1-day rainfall was 56.36 mm and the maximum 3-day rainfall was 70.96 mm at the 3-year return period.
The Wanyan River surface waterlogged area land use data were obtained from the China Land Use/Cover Dataset (CLUD) in 1990, 1995, 2000, 2005, 2010, 2015, 2018, and 2020. The data set was provided by the Data Center for Resources and Environmental Sciences, Chinese Academy of Sciences (RESDC) (http://www.resdc.cn (accessed on 15 November 2022)). From 1990 to 2010, Liu et al. [34,35,36] compared the actual field survey data with cartographic results to evaluate the accuracy. CLUD first-level accuracy was greater than 80%. The accuracy of the second-level CLUD, which started in 2010, was 91.2%. Based on the actual field investigation records from 2010 to 2020, Liu et al. [37,38,39] used a confusion matrix to analyze the accuracy of CLUD first- and second-level type mapping, both of which were above 90%. The land use of the study area can be divided into five categories, that is, dry land, paddy field, wood meadows, water bodies, and construction land [40]. The land use for other years was obtained by linear interpolation. The drainage flow of the Wanyan River pumping station was collected and calculated to represent the actual runoff, which is hereinafter referred to as the observed runoff.

2.3. Selection of Drainage Criteria

The basic requirement of farmland waterlogging control is to discharge the excess surface water generated by rainstorms with a certain frequency below a certain depth within the specified time in order to protect crops from being flooded. The drainage criteria, which is the key index for determining the surface drainage coefficient, is defined by the design rainfall return period, rainfall duration, and drainage duration. Based on current drainage criteria and local conditions, the design criteria of discharging 1-day rainfall in 2 days and 3-day rainfall in 4 days are selected separately to compute the drainage coefficient.

2.4. Methods

2.4.1. Theoretical Runoff Calculation Method

Different land use types used different methods to calculate theoretical runoff. The runoff of paddy fields (Rpf, in mm) was calculated by the infiltration curve method. The runoff of dry fields (Rdf, in mm) was calculated by the saturation-excess runoff formula [32] using the 1-layer evaporation model. The runoff of wood meadows (Rwm, in mm) was regarded as 40% of dry field runoff. The runoff of construction land (Rcl, in mm) was calculated by the runoff coefficient. Combining the area and runoff results of different land use types in different years, the yearly runoff for 1-day and 3-day rainfall at 3-, 5-, 10-, and 20-year return period was calculated according to the area weighting method. Therefore, theoretical runoff was calculated as follows [41]:
R pf = P h ( E + D ) t R df = P ( W M W 0 ) E 0 W 0 W M R w m = 40 % R d f R cl = α P R = f p f R p f + f d f R d f + f w m R w m + f c l R c l
where R is theoretical runoff (mm); f p f , f d f , f w m , and f c l are the area proportion of paddy fields, dry lands, wood meadows, and construction land, respectively (%); P is rainfall (mm); h is the depth of paddy submergence tolerance (mm); E is the evaporation rate of paddy fields (mm/d); D is the infiltration rate of the surface water (mm/d); W0 is the average initial water storage (mm); WM is the average maximum water storage (mm); E0 is water surface evaporation (mm); and α is the runoff coefficient. According to the relevant references [16,41,42,43], h = 30, E = 4.29, D = 3.59, W0 = 40, WM = 80, E0 = 3.1, and α = 0.3.

2.4.2. Empirical Formula

Based on the relationship between drainage coefficient, design runoff depth, and basin area, the formula can be expressed as follows [26,32]:
Q E = 86.4 K R m F n
where QE is the drainage coefficient (mm/d); F is the drainage area (km2); R is the design runoff (mm); K is the comprehensive coefficient, which is related to the river network matching level, the slope of drainage ditches, rainfall duration, the basin shape, and other factors; m is the index reflecting the relationship between the flood peak and volume; and n is the descending exponent reflecting the relationship between the drainage coefficient and the drainage area. The related parameters used in the study area refer to those used in the central plain area of Liaoning Province, K = 0.0127, m = 0.93, and n = −0.176 [44,45].

2.4.3. Average Draining Formula

The average draining formula under pumping drainage conditions is as follows:
Q A = 86.4 R φ 3.6 T t
where QA is the drainage coefficient (mm/d); R is the design runoff (mm); φ is the channel storage coefficient, which is related to the basin area; T is the drainage duration (d); and t is set as 24 h under gravity drainage conditions and 23 h under pumping drainage conditions, indicating the pumping station runtime in a day.

2.4.4. Drainage Duration (TE) Corresponding to Empirical Formula

When the catchment area in the plain area is large, the process of runoff generation will be slow. The catchment area’s regulating function during the drainage process makes the drainage flow average and the peak flow reduction. The duration of the regulating time will have an impact on the drainage flow, while the empirical formula does not consider the influence of the change in the drainage duration on the drainage coefficient.
In order to estimate the drainage duration (TE, in d) corresponding to the empirical formula, it is assumed that the drainage flow calculated by the empirical formula is discharged averagely and the drainage duration can be deduced as follows:
T E = 86.4 R φ 3.6 23 Q E

2.4.5. Area (FA) Corresponding to the Average Draining Formula

The average draining formula was mostly used in the Heilongjiang Province to calculate the drainage coefficient. The researchers did not consider the applicable area range of this method as well as the impact of the area size on the drainage coefficient. To study the area (FA, in km2) corresponding to the average draining formula, it is assumed that the area has an impact on the drainage coefficient calculated by the average draining formula. In other words, QA is equal to QE, and the following equation can be deduced:
F A = ( Q A 86.4 K R m ) 1 n

3. Results

3.1. Runoff Analysis

3.1.1. Observed Runoff Estimation

Seven sets of rainfall and drainage data with 1-day rainfall and 3-day rainfall were selected separately in the study area, as shown in Table 1.
As the 3-day rainfall of the number seven event was 72.6 mm in Table 1, which was close to the maximum 3-day rainfall of 70.96 mm, its return period was regarded as 3-year. Apart from the number seven event, the rainfall of the other six events was less than the 3-year return period. In this paper, they were all defaulted to be less than the 3-year return period, and no specific division was made.

3.1.2. Changes of Cultivated Land in the Past 30 Years

Figure 2 describes the change of land use area in the study area in the last 30 years. The area of farmland increased continually from 1990 to 2020 and accounted for more than 80% of the whole area, with a minimum ratio of 83.7% in 1990 and a maximum ratio of 93.43% in 2020. In the meantime, the area of paddy fields increased rapidly from 0.73% in 1990 to 72.1% in 2020 with an annual increase of 2.37%, of which the growth trend was slow from 1990 to 2010 and much faster from 2010 to 2018. On the contrary, the area of dry land decreased continually from 82.97% in 1990 to 21.33% in 2020 with an annual decrease rate of 2.05%. The ratio of paddy area to dry area showed a continuous increasing trend, which was the same as that in the paddy field. The water area remained steady before 2015 and reduced rapidly to 0.5% in 2020. The construction land showed a fluctuating increasing trend, but the trend was not obvious. The area of wood meadows decreased at first, reaching its minimum value in 2015, and then increased at a slow growth rate after 2015.

3.1.3. Theoretical Runoff Estimation

The years with observed runoff data and the years with runoff turning points (where the runoff changed from increasing to decreasing or from decreasing to increasing) were chosen in the theoretical runoff as typical years, and the results are shown in Table 2.
The runoff increased when the return period was longer under the same rainfall duration or the rainfall duration was longer under the same return period in Table 2. There was a positive correlation between the return period and the runoff and between the rainfall duration and the runoff.
The annual runoff change rate of 1-day rainfall increased with the increase in the return period. The runoff of 3-day rainfall decreased slowly with time, so the annual runoff change rates were negative values. According to the absolute value of the runoff change rate under 3-day rainfall, the annual runoff change trend decreased with the increase in the return period, and the annual runoff change trend of 1-day rainfall was greater than that of 3-day rainfall.

3.1.4. Comparison of Runoff

Both the observed runoff and the theoretical runoff increased when the rainfall was higher and the rainfall duration was longer. Although the observed runoff was not strongly correlated with rainfall, it also showed a positive correlation with rainfall.
The observed runoff was lower than the theoretical runoff. Under 3-day rainfall with a 3-year return period, the observed runoff was 8.17 mm, while the theoretical runoff was 18.93 mm. The actual rainfall in other years was lower than the design rainfall, so the observed runoff was lower than the theoretical runoff. It was observed that the Heilongjiang River and Songhua River were at high water levels when the rainfall was heavy in summer. Limited by the capacity of the pumping station and the high river water level at the boundary of the study area, the runoff generated by a rain event could not be fully drained by the pumping station, which could cause excess water in the farmland. To solve this problem, the workers temporarily add some small pumps to pump water from the drainage ditch. The observed runoff in this study did not include these temporary pumping volumes, which was one of the reasons for the small observed runoff.
Comparing number four and number five events under the 1-day rainfall duration in Table 1, it was noticed that the runoff generated by a light rain was more than that generated by a heavy rain, which is inconsistent with the general law. This phenomenon showed that some treatments might lead to the deviation of the results when calculating the observed runoff. For example, the larger the basin area, the longer the rainfall confluence time; especially when the groundwater confluence time was long and the drainage process of this rainfall was not completed before the next rainfall came, it was difficult to accurately divide two rainfall-runoff processes. Compared with the observed runoff, the theoretical runoff affected by land use change was more regular, rational, and reliable. Therefore, the theoretical runoff was selected for the calculation of the drainage coefficient.

3.2. Analysis of Drainage Coefficient

According to the results of theoretical runoff, the drainage coefficients calculated by the average draining method and the empirical formula are shown in Table 3 and Table 4. Considering the current actual farmland drainage demand, the table only lists the drainage coefficient results for the 3-, 5-, and 10-year return periods.
As shown in Table 3, QA was greater than QE when the return period was greater than 3-year, and less than QE when the return period was 3-year under 1-day rainfall duration. Table 4 shows that QE was always greater than QA for any return period under 3-day rainfall duration. The longer the drainage duration, the greater the regulating effect of farmland on runoff. The average draining formula could effectively reflect the regulating influence of the drainage duration on runoff. Thus, the QA was small when discharging 3-day rainfall in 4 days.
It can be seen from Table 3 and Table 4 that QA and QE both decreased with the decrease in the rainfall return period under the same rainfall duration in the same year. QE increased with the increase in rainfall duration under the same rainfall return period in the same year. The effect of rainfall duration on QA had no obvious regularity. For 1-day rainfall duration, the drainage coefficient increased with time, and its annual change rate decreased with the decrease in the return period. The annual change trend of QA was greater than that of QE when the return period was greater than 3-year. The drainage coefficient decreased with time for the 3-day rainfall duration. According to the absolute value of the drainage coefficient change rate under 3-day rainfall duration, the change trend of the annual drainage coefficient increased with the decrease in the return period, and that of 3-day rainfall duration was less than that of 1-day rainfall duration. The change trend of the annual drainage coefficient was the same as that of theoretical runoff.
As shown in Table 1, the observed runoff of the Wanyan River pumping station in 2018 was 8.17 mm under the 3-day rainfall duration with 3-year return period. Therefore, the observed QE was 2.21 mm/d, and QA was 0.98 mm/d when the drainage duration was 4 days, so the observed drainage flow of the Wanyan River pumping station was 13.97–31.49 m3/s. As can be seen from Table 4, under 3-day rainfall duration with 3-year return period in 2018, the theoretical QE was 4.83 mm/d, QA was 2.47 mm/d, and the corresponding drainage flow of the pumping station was 35.11–68.73 m3/s. The theoretical flow was greater than the actual flow. In the same way, under 5-year return period in 2018, the theoretical drainage flow was 79.74–87.32 m3/s for 1-day rainfall duration and 67.97–105.6 m3/s for 3-day rainfall duration. The design flow of the Wanyan River pumping station is 57.97 m3/s, which cannot meet the drainage demand of the 5-year return period and is also an important reason for the local waterlogging disaster.

3.3. Analysis on Drainage Duration of Empirical Formula

By substituting the drainage coefficient calculated by the empirical formula into Formula (3), the drainage duration required when the drainage flow calculated by the empirical formula was discharged at an average rate can be obtained. The results are shown in Table 5.
As shown in Table 5, TE varied with rainfall and increased with the increase in rainfall duration and return period. According to the average value, under 1-day rainfall duration, TE of the 5-year return period was 15.14% higher than that of the 3-year return period, and TE of the 10-year return period was 30.05% higher than that of the 5-year return period. Under 3-day rainfall duration, TE of the 5-year return period was 19.46% higher than that of the 3-year return period, and TE of the 10-year return period was 16.29% higher than that of the 5-year return period. Under 1-day rainfall duration, TE increased with time and the annual change of TE was small. Under 3-day rainfall duration, TE decreased with time, and its annual change trend increased with the decrease in return period.
When the design drainage duration was equal to TE, QA was equal to QE. TE ranged from 2.11 d to 2.82 d when the 1-day rainfall greater than the 3-year return period was discharged in 2 days, TE was greater than the design drainage duration, and QA was larger than QE, so it was safer to use the average draining formula to calculate the drainage coefficient during the construction of drainage system. TE ranged from 1.83 d to 1.87 d when discharging 1-day rainfall in 2 days with 3-year return period and from 2.04d to 3.09 d when discharging 3-day rainfall in 4 days. In the above two cases, the design drainage durations were greater than the TE, respectively, so QE was greater than QA, and it was safer to apply the empirical formula to calculate the drainage coefficient.
The results of the 1-day rainfall duration with the 3-year return period are different from the existing literature, but the other results are the same as those in the existing literature [24]. The reason might be that the existing literature only considers the influence of the design drainage duration on TE, and the effects of rainfall duration and return period are also considered in our research.
Therefore, the longer drainage duration can be used for places with good regulating effects. In order to avoid idle pumps and waste human resources and funds, the empirical formula can be selected for engineering design; at the same time, the economic feasibility of engineering also needs to be considered.

3.4. Analysis on the Area of Average Draining Formula

Substituting the parameters of the empirical formula into Formula (4), the FA under different rainfall durations and different return periods was calculated. The results are shown in Table 6.
As shown in Table 6, FA ranged from 173 to 2029 km2 for 1-day rainfall duration and from 5353 to 55,883 km2 for 3-day rainfall duration. For the same rainfall duration, FA kept rising as the return period decreased. FA increased with the increase in rainfall duration under the same return period. Therefore, the smaller the rainfall duration and the greater the return period, the smaller FA. On the contrary, the larger the rainfall duration and the smaller the return period, the larger FA. FA decreased with time under the 1-day rainfall duration and increased with time under the 3-day rainfall duration. According to the absolute value of the annual FA change rate, the annual FA change rate of 1-day rainfall duration was smaller than that of the 3-day rainfall duration, and both of them increased with the decrease in the return period.
FA was small under 1-day rainfall duration and relatively large under the 3-day rainfall duration. Under 1-day rainfall duration with a 10-year return period, FA was the lowest, varying from 173 to 203 km2. Under 3-day rainfall duration with a 3-year return period, FA was the highest, varying from 24,407 to 55,883 km2. It was evident that the average draining formula gave the small results for long duration of rainfall and drainage, making it potentially unsafe for application in engineering design.

4. Discussion

Runoff has an important influence on the drainage coefficient, whether the drainage coefficient is calculated using the average draining formula or the empirical formula. The difference is that the drainage coefficient calculated by the empirical formula varies with the change of catchment area. The empirical formula considers the impact of runoff, topography, and catchment shape on the waterlogged area, so it requires higher timeliness of parameters and the calculation results are more realistic. The drainage duration has a significant impact on the drainage coefficient calculated by the average draining formula. In our research, the channel storage coefficient was added to the average draining formula, which reflected the effect of area on the drainage coefficient.
Whether the empirical formula or the average draining formula is used, the related parameters should be modified, such as the parameters of the empirical formula, the paddy field stagnant water depth in the average draining formula, cultivated land infiltration coefficient, etc. Further study is needed in the future because our research used the regional empirical parameters directly. These two methods do not systematically evaluate the impact of underlying surface changes on drainage flow, nor do they consider the flow change during a drainage process. A model will be used to link the discharge coefficient with the peak of the hydrograph.
In our research, it was assumed that the size of the waterlogged area has an impact on the drainage coefficient calculated by the average draining formula. This assumption was used to calculate FA for different rainfall durations and return periods. In reality, the area size has an impact on the catchment’s runoff, which will change the drainage coefficient. Therefore, when the results in this paper are applied to other areas, the FA needs to be further confirmed.

5. Conclusions

The observed runoff was obtained by drainage data from the Wanyan River pumping station, and the theoretical runoff was calculated using the surface runoff calculation method based on the rainfall and land use data. Based on the theoretical runoff, the drainage coefficient was calculated using the empirical formula and the average draining formula. The main conclusions were as follows.
(1)
As rainfall and rainfall duration increased, both the observed runoff and the theoretical runoff increased. The observed runoff was typically less than the theoretical runoff due to the limitation of the pumping capacity and the high water levels in the Heilongjiang River and Songhua River. The drainage coefficient calculated by theoretical runoff showed that the design flow of the Wanyan River pumping station could not meet the drainage demand under 5-year return period rainfall.
(2)
When the 1-day rainfall greater than the 3-year return period was discharged in 2 days, the results of the average draining formula were greater than those of the empirical formula, and the TE was greater than the design drainage duration. In this condition, it was safer to use the average draining formula to calculate the drainage coefficient for drainage engineering design. Under the 1-day rainfall duration with a 3-year return period and discharging 3-day rainfall in 4 days, the results of the empirical formula were greater than those of the average draining formula, and the design drainage duration was greater than TE. In this condition, it was safer to use the empirical formula to calculate the drainage coefficient. Therefore, when calculating the drainage coefficient, a longer drainage duration could be selected for places with good water storage conditions.
(3)
FA ranged from 173 to 2029 km2 for 1-day rainfall duration and from 5353 to 55,883 km2 for 3-day rainfall duration. FA increased with the increase in rainfall duration and decreased with the increase in return period. When the rainfall duration and drainage duration were long, the calculation results of the average draining formula were much smaller, which was unsafe for engineering design.

Author Contributions

Conceptualization, S.W.; methodology, L.Y., S.W., Y.T. and Y.L.; validation, S.W., Y.T. and Y.L.; formal analysis, S.W., Y.T. and Y.L.; investigation, L.Y., Y.L. and Y.T.; data curation, L.Y. and Y.L.; writing—original draft preparation, L.Y., Y.T. and Y.L.; writing—review and editing, L.Y., Y.T.; project administration, S.W.; funding acquisition, S.W. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the National Natural Science Foundation Program of China (No. 52079145, No. 51779274).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The datasets used and analyzed during the current study are available from the corresponding authors upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. IPCC. Summary for policymakers. In Climate Change 2021: The Physical Science Basis; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2021. [Google Scholar]
  2. Liu, S.M.; Wang, H.; Yan, D.H.; Qin, T.L.; Wang, Z.L.; Wang, F.X. Crop Growth Characteristics and Waterlogging Risk Analysis of Huaibei Plain in Anhui Province, China. J. Irrig. Drain. Eng. 2017, 143, 04017042. [Google Scholar] [CrossRef]
  3. Milly, P.C.D.; Betancount, J.; Falkenmark, M.; Hirsch, R.M.; Kundzewicz, Z.W.; Lettenmaier, D.P.; Souffer, R.J. Stationarity Is Dead: Whither Water Management? Science 2008, 319, 573–574. [Google Scholar] [CrossRef] [PubMed]
  4. Beniston, M.; Stephenson, D.B.; Christensen, O.B.; Ferro, A.T.F.; Frei, C.; Goyette, S.; Halsnaes, K.; Holt, T.; Jylhä, K.; Koffi, B.; et al. Future extreme events in European climate: An exploration of regional climate model projections. Clim. Change 2007, 81, 71–95. [Google Scholar] [CrossRef] [Green Version]
  5. IPCC. Climate Change 2013: The Physical Science Basis; An Overview of the Working Group 1 Contribution to the Fifth Assessment Report; IPCC: Geneva, Switzerland, 2013; p. 3544. [Google Scholar]
  6. Mwangi, H.; Julich, S.; Patil, S.; Mcdonald, M.; Feger, K.H. Modelling the impact of agroforestry on hydrology of Mara River basin in East Africa using a distributed model. In Proceedings of the Egu General Assembly Conference, Vienna, Austria, 17–22 April 2016. [Google Scholar]
  7. Liu, S.M.; Wang, H.; Yan, D.H.; Ren, Q.F.; Wang, D.; Gong, B.Y. Analysis of Spatiotemporal Evolution of Isolated Rainstorm Events in Huai River Basin, China. Adv. Meteorol. 2017, 2017, 3010295. [Google Scholar] [CrossRef] [Green Version]
  8. Norton, U.; Mosier, A.R.; Morgan, J.A.; Derner, J.D.; Ingramd, L.J.; Stahld, P.D. Moisture pulses, trace gas emissions and soil C and N in cheatgrass and native grass-dominated sagebrush-steppe in Wyoming, USA. Soil Biol. Biochem. 2008, 40, 1421–1431. [Google Scholar] [CrossRef]
  9. Greenway, H.; Gibbs, J. Review: Mechanisms of anoxia tolerance in plants. II. Energy requirements for maintenance and energy distribution to essential processes. Funct. Plant Biol. 2003, 30, 999–1036. [Google Scholar] [CrossRef]
  10. Dat, J.F.; Capelli, N.; Folzer, H.; Bourgeade, P.; Badot, P.-M. Sensing and signalling during plant flooding. Plant Physiol. Biochem. 2004, 42, 273–282. [Google Scholar] [CrossRef]
  11. Tian, L.X.; Zhang, Y.C.; Chen, P.L.; Zhang, F.F.; Feng, B.L. How Does the Waterlogging Regime Affect Crop Yield? A Global Meta-Analysis. Front. Plant Sci. 2021, 12, 634898. [Google Scholar] [CrossRef]
  12. Chen, H.R.; Zeng, W.Z.; Jin, Y.L.; Mi, B.Y.; Zha, Y.Y.; Zhang, S.H. Development of a waterlogging analysis system for paddy fields in irrigation districts. J. Hydrol. 2020, 591, 125325. [Google Scholar] [CrossRef]
  13. Brody, S.D.; Peacock, W.G.; Gunn, J. Ecological Indicators of Flood Risk along the Gulf of Mexico. Ecol. Indic. 2012, 18, 493–500. [Google Scholar] [CrossRef]
  14. Dabral, P.P.; Kumar, A.; Tana, G. Determination of surface drainage coefficient—A case study of Doimukh (Arunachal Pradesh), India. Agric. Eng. Int. CIGR E-J. 2016, 18, 1–10. [Google Scholar]
  15. Merz, R.; Blöschl, G.; Parajka, J. Spatio-temporal variability of event runoff coefficients. J. Hydrol. 2006, 331, 591–604. [Google Scholar] [CrossRef]
  16. Ministy of Housing and Urban-Rural Develpoment of the People’s Republic of China. GB 50014-2021; Code for Design of Outdoor Wastewater Engineering. China Planning Press: Beijing, China, 2021; pp. 13–14.
  17. Grimaldi, S.; Petroselli, A.; Romano, N. Green-Ampt Curve Number mixed procedure as an empirical tool for rainfall-runoff modelling in small and ungauged basins. Hydrol. Process. 2013, 27, 1253–1264. [Google Scholar] [CrossRef]
  18. Li, J.; Wang, Z.G.; Liu, C.M. A combined rainfall infiltration model based on Green-Ampt and SCS-curve number. Hydrol. Process. 2015, 29, 2628–2634. [Google Scholar] [CrossRef]
  19. Wang, Q.J.; Lai, J.B.; Yi, L. Comparison of Green-Ampt Model with Philip Infiltration Model. Trans. Chin. Soc. Agric. Eng. 2002, 18, 13–16. [Google Scholar]
  20. Song, X.M.; Kong, F.Z.; Zhan, C.S.; Han, J.W. Hybrid Optimization Rainfall-Runoff Simulation Based on Xinanjiang Model and Artificial Neural Network. J. Hydrol. Eng. 2012, 17, 1033–1041. [Google Scholar] [CrossRef]
  21. Kumar, A.; Kanga, S.; Taloor, A.K.; Singh, S.K.; Đurin, B. Surface runoff estimation of Sind river basin using integrated SCS-CN and GIS techniques. Hydro Res. 2021, 4, 61–74. [Google Scholar] [CrossRef]
  22. Ruggiero, C. Bnifiche Idauliche; Milani: Bologna, Italy, 1930. [Google Scholar]
  23. Pasin, P. Coefficient Udometrici Desunti dal Lavoro Delle Machine Nelle Bonifiche Mecaniche, Agosto, Italy, 1940.
  24. Praveen, K.; Pradeep, H.; Lokesh, A.; Akarshraj, K.H.; Surendra, H.J.; Avinash, S.D. Estimation of Runoff using Empirical Equations and Fuzzy Logic method: A case study. Int. J. Sci. Eng. Res. 2016, 7, 5. [Google Scholar]
  25. Smedema, L.K.; Mooren, J.H.M.; Veraa, C.F.J.M. Area Reduction Factors for Land Drainage Design in Flat Basins in Moderate Temperate Climates. In Hydraulic Design in Water Resources Engineering: Land Drainage; Rycroft, D.W., Ed.; Springer: Berlin/Heidelberg, Germany, 1986; pp. 87–99. [Google Scholar] [CrossRef]
  26. Luo, W.B.; Wang, X.G.; Luo, Q. Comparison of methods for calculating farmland drainage modulus. Trans. Chin. Soc. Agric. Eng. 2013, 29, 85–91. [Google Scholar] [CrossRef]
  27. Wang, J.F.; Zhu, M.; Wang, Z.l.; Zhang, T.T.; Qi, Y.F.; Zhu, L.L. Improvement of farmland drainage calculation method under environmental change in Huaibei plain. Irrig. Drain. 2021, 70, 1347–1358. [Google Scholar] [CrossRef]
  28. Wang, G.A.; He, S.D.; Li, R.R.; Li, C.Q.; Gai, Y.G.; Wang, P. On Basin Principle and Applicable Conditions of Reasoning Formula. Yellow River 2010, 32, 4. [Google Scholar] [CrossRef]
  29. Paul, J.C.; Panigrahi, B.; Behera, B.P.; Patnaik, P.K. Design of surface drainage system for managing water congestion in mahanadi delta: A case study. J. Agric. Eng. 2014, 51, 44–50. [Google Scholar]
  30. Dabral, P.P.; Baithuri, N. Determination of surface drainage coefficient through rainfall analysis for North Lakhimpur (Assam). J. Soil Water Conserv. 2008, 7, 38–42. [Google Scholar]
  31. Kong, X.M.; Hao, Z.D.; Zhu, Y. Entropy Theory and Pearson Type-3 Distribution for Rainfall Frequency Analysis in Semi-arid Region. IOP Conf. Ser. Earth Environ. Sci. 2020, 495, 012042. [Google Scholar] [CrossRef]
  32. Zhan, D.J.; Xu, X.Y.; Chen, Y.F. Engeering Hydrology, 4th ed.; China Water & Power Press: Beijing, China, 2010. [Google Scholar]
  33. You, L.J.; Wang, S.L.; Liu, C.R.; Zhang, X.B.; Fan, L.S.; Tao, Y.; Peng, Z. Characteristics analysis of rainfall-induced flood and waterlogging in the Sanjiang Plain from 1957 to 2019: A case study of Suibin County. J. Water Resour. Water Eng. 2022, 33, 40–49. [Google Scholar] [CrossRef]
  34. Liu, J.Y.; Liu, M.l.; Deng, X.Z.; Zhuang, D.F.; Zhang, Z.X.; Luo, D. The land use and land cover change database and its relative studies in China. J. Geogr. Sci. 2002, 12, 275–282. [Google Scholar]
  35. Liu, J.Y.; Liu, M.L.; Zhuang, D.F.; Deng, X.Z.; Zhang, Z.X. Recent Spatial pattern analysis of land use change in China. Sci. Sin. 2002, 1031–1040+1058–1060. [Google Scholar]
  36. Liu, J.Y.; Zhang, Z.X.; Xu, X.L.; Kuang, W.H.; Zhou, W.C.; Zhang, S.W.; Li, R.D.; Yan, C.Z.; Yu, D.S.; Wu, S.X.; et al. Spatial Patterns and Driving Forces of Land Use Changein China in the Early 21st Century. Acta Geogr. Sin. 2009, 64, 1411–1420. [Google Scholar]
  37. Liu, J.Y.; Kuang, W.H.; Zhang, Z.X.; Xu, X.L.; Qin, Y.W.; Ning, J.; Zhou, W.C.; Zhang, S.W.; Li, R.D.; Yan, C.Z.; et al. Spatiotemporal characteristics, patterns, and causes of land-use changes in China since the late 1980s. J. Geogr. Sci. 2014, 24, 195–210. [Google Scholar] [CrossRef]
  38. Ning, J.; Liu, J.Y.; Kuang, W.H.; Xu, X.L.; Zhang, S.W.; Yan, C.Z.; Li, R.D.; Wu, S.X.; Hu, Y.F.; Du, G.M.; et al. Spatiotemporal patterns and characteristics of land-use change in China during 2010–2015. J. Geogr. Sci. 2018, 28, 547–562. [Google Scholar] [CrossRef] [Green Version]
  39. Kuang, W.H.; Zhang, S.W.; Du, G.M.; Yan, C.Z.; Wu, S.X.; Li, R.D.; Lu, D.S.; Pan, T.; Ning, J.; Guo, C.Q.; et al. Remotely sensed mapping and analysis of spatio-temporal patterns of land use change across China in 2015–2020. Acta Geogr. Sin. 2022, 77, 1056–1071. [Google Scholar]
  40. Ministry of Land and Resources of the People’s Republic of China. GB/T 21010-2017; Current Land Use Classification. Standards Press of China: Beijing, China, 2017; pp. 1–6.
  41. Luo, W.B. Effects of Underlying Surface Changes on Farmland Drainage Discharge; Wuhan University: Wuhan, China, 2014. [Google Scholar]
  42. Zhou, R.Y. A complete analysis of the scale review process of Suibin Pump Station renovation project. Heilongjiang Sci. Technol. Water Conserv. 2011, 39, 42–43. [Google Scholar]
  43. Zhu, Y.C.; Liu, Q.; Zhou, G. Analysis of drainage coefficient in Naoli River surface waterlogged area of Raohe County. Water Conserv. Sci. Technol. Econ. 2009, 15, 629–630. [Google Scholar]
  44. Ministry of Water Resources of the People’s Republic of China. SL 723-2016; Standard for Waterlogging Control. China Water & Power Press: Beijing, China, 2016; p. 54.
  45. Guo, Y.Y. Farmland Water Conservancy, 3rd ed.; China Water & Power Press: Beijing, China, 2007. [Google Scholar]
Figure 1. Location of the study area.
Figure 1. Location of the study area.
Applsci 13 02769 g001
Figure 2. Change of land use area in the study area in recent 30 years.
Figure 2. Change of land use area in the study area in recent 30 years.
Applsci 13 02769 g002
Table 1. Observed runoff under different rainfall durations.
Table 1. Observed runoff under different rainfall durations.
FID1-Day Rainfall3-Day Rainfall
DateRainfall
(mm)
Drainage
(104 m3)
Runoff (mm)DateRainfall (mm)Drainage
(104 m3)
Runoff (mm)
111 September 199114.744.050.3528 August–1 September 198814.2197.271.60
2 26 August 2009 18.2 52.670.4316 August–18 August 199148.3344.022.80
3 28 August 2009 21.1 90.970.743 October–5 October 1994 23.3 201.11.60
4 21 July 2014 8.9 38.300.3119 May–21 May 2010 33.3 121.620.99
5 23 July 2014 31.8 21.070.1714 August–16 August 2010 28.5 118.260.96
6 5 August 2016 7.5 7.660.0611 August–13 August 2014 9.5 50.750.41
722 September 201817.815.030.1224 July–26 July 201872.61004.348.17
Table 2. Theoretical runoff under 1-day and 3-day rainfall with different return periods (mm).
Table 2. Theoretical runoff under 1-day and 3-day rainfall with different return periods (mm).
Year1-Day Rainfall3-Day Rainfall
20-Year10-Year5-Year3-Year20-Year10-Year5-Year3-Year
199139.2630.5017.6013.0155.9444.8532.8222.81
199540.3331.3318.0913.3757.4046.0133.6523.38
200941.1032.0718.8014.0755.8744.4732.0721.78
201041.1232.1018.8214.0955.8244.4232.0221.72
201441.9332.8719.5414.8054.3942.9530.5020.16
201542.1333.0619.7214.9754.0342.5830.1219.77
201643.3034.0320.4115.5654.3242.6329.9119.35
201846.1636.4622.1817.1055.5843.3430.0018.93
Annual change rate25.56%22.07%16.96%15.15%−1.33%−5.59%−10.44%−14.37%
Table 3. Drainage coefficient under discharging 1-day rainfall in 2 days with different return periods (mm/d).
Table 3. Drainage coefficient under discharging 1-day rainfall in 2 days with different return periods (mm/d).
Year10-Year5-Year3-Year
QEQAQEQAQEQA
19917.5310.344.524.783.413.12
19957.7210.624.634.913.503.21
20097.8910.884.805.103.673.38
20107.9010.884.815.113.673.38
20148.0711.154.985.303.843.55
20158.1211.215.025.353.893.59
20168.3411.545.185.544.033.74
20188.8912.555.606.134.404.10
Annual change rate5.04%8.19%4.01%5.03%3.65%3.63%
Table 4. Drainage coefficient under discharging 3-day rainfall in 4 days with different return periods (mm/d).
Table 4. Drainage coefficient under discharging 3-day rainfall in 4 days with different return periods (mm/d).
Year10-Year5-Year3-Year
QEQAQEQAQEQA
199110.788.318.065.395.753.39
199511.048.528.255.535.883.48
200910.698.247.895.275.513.18
201010.688.237.885.265.493.17
201410.357.967.534.935.132.74
201510.277.897.444.875.032.63
201610.287.907.404.844.932.52
201810.448.037.424.774.832.47
Annual change rate−1.25%−1.04%−2.39%−2.29%−3.39%−3.42%
Table 5. Drainage duration corresponding to the empirical formula under 1-day and 3-day rainfall with different return periods (TE, in d).
Table 5. Drainage duration corresponding to the empirical formula under 1-day and 3-day rainfall with different return periods (TE, in d).
Year1-Day Rainfall3-Day Rainfall
10-Year5-Year3-Year10-Year5-Year3-Year
19912.752.111.833.082.682.36
19952.752.121.833.092.682.36
20092.762.121.843.082.672.31
20102.762.121.843.082.672.31
20142.762.131.853.072.622.13
20152.762.131.853.072.622.09
20162.772.141.853.072.622.05
20182.822.191.873.072.572.04
Average2.772.131.853.072.642.21
Annual change rate0.26%0.30%0.15%−0.04%−0.41%−4.52%
Table 6. The area corresponding to the average draining formula under 1-day and 3-day rainfall with different return periods (FA, in km2).
Table 6. The area corresponding to the average draining formula under 1-day and 3-day rainfall with different return periods (FA, in km2).
Year1-Day Rainfall3-Day Rainfall
10-Year5-Year3-Year10-Year5-Year3-Year
19912038972029540812,07724,646
19952018872008535311,95724,407
20091998741967542612,18827,762
20101998731966542912,19627,789
20141978601928550213,61743,612
20151968571919552113,68649,080
20161948451890551813,72455,405
20181737341821548215,03455,883
Annual change rate−1.11−6.04−7.702.74109.521156.93
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

You, L.; Wang, S.; Tao, Y.; Liu, Y. Appropriate Method to Estimate Farmland Drainage Coefficient in the Wanyan River Surface Waterlogged Area in Suibin County of the Sanjiang Plain, China. Appl. Sci. 2023, 13, 2769. https://doi.org/10.3390/app13052769

AMA Style

You L, Wang S, Tao Y, Liu Y. Appropriate Method to Estimate Farmland Drainage Coefficient in the Wanyan River Surface Waterlogged Area in Suibin County of the Sanjiang Plain, China. Applied Sciences. 2023; 13(5):2769. https://doi.org/10.3390/app13052769

Chicago/Turabian Style

You, Lijun, Shaoli Wang, Yuan Tao, and Yongji Liu. 2023. "Appropriate Method to Estimate Farmland Drainage Coefficient in the Wanyan River Surface Waterlogged Area in Suibin County of the Sanjiang Plain, China" Applied Sciences 13, no. 5: 2769. https://doi.org/10.3390/app13052769

APA Style

You, L., Wang, S., Tao, Y., & Liu, Y. (2023). Appropriate Method to Estimate Farmland Drainage Coefficient in the Wanyan River Surface Waterlogged Area in Suibin County of the Sanjiang Plain, China. Applied Sciences, 13(5), 2769. https://doi.org/10.3390/app13052769

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