Next Article in Journal
Experimental Assessment of the Use of a Novel Superabsorbent polymer (SAP) for the Optimization ofWater Consumption in Agricultural Irrigation Process
Previous Article in Journal
Economic Assessment of Opportunities for Managed Aquifer Recharge Techniques in Spain Using an Advanced Geographic Information System (GIS)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Characteristics of Water Budget Components in Paddy Rice Field under the Asian Monsoon Climate: Application of HSPF-Paddy Model

1
Civil Engineering Division, Samsung C&T Corp., Seoul 137-858, Korea
2
Research Group on Agricultural Engineering, Rural Research Institute, Ansan 426-908, Korea
3
Department of Environmental Engineering, Andong National University, Andong 760-749, Korea
*
Author to whom correspondence should be addressed.
Submission received: 14 May 2014 / Revised: 26 June 2014 / Accepted: 9 July 2014 / Published: 17 July 2014

Abstract

: The HSPF-Paddy model was applied to the Bochung watershed in Korea to compare water budget components by the land use types under the Asian monsoon climate. The calibration of HSPF-Paddy during 1992–2001 with PEST, a package program to optimize HSPF, and validation during 1985–1991 were carried out. The model efficiencies for monthly stream flow are 0.85 for calibration and 0.84 for validation. The simulation of annual mean runoff met the criteria of water budget analysis with the acceptable error level (less than 10 percent mean error). The simulation of the movement of water from paddy rice field to watershed was successful, and application of HSPF-Paddy coupled with PEST was able to improve accuracy of model simulation with reduced time and efforts for model calibration. The results of water budget analysis show that most of the outflow (86%) for the urban area occurred through surface runoff, showing the highest rate among the land use types compared. Significant amounts of water are irrigated to paddy rice fields, and the runoff depth as well as evapotranspiration from paddy rice field is higher than other land use types. Hydrological characteristic of paddy rice field is that most of water movement occurred at the surface area, resulting from the low infiltration rate and manning’s coefficient, as well as ponded water throughout the growing season. Major impact on input and output of water were precipitation and runoff, respectively, influenced by an Asian monsoon climate.

1. Introduction

Agricultural in Asian monsoon region feeds about 60% of the world’s population by using 30% of the arable land, and by approximately 40 years after World War II, most Asia countries had achieved self-sufficiency in rice [1]. The rainfall in an Asian monsoon region is concentrated on crop-growing season especially summer, so most of the nonpoint source pollution from crop land occurred during summer season. In paddy rice fields, dikes and forced drain alter the hydrological characteristics, as well as the effects from excessive nutrient supply and fertilization to ponded water quality [2]. These characteristics of paddy rice fields make simulation of pollutant load and water movement more cumbersome.

Many researchers have tried to modify existing models or develop a new model to simulate the water budget, pollutant behavior in paddy rice fields. Most of them are field scale models such as PADDIMOD [3], GLEAMS-PADDY [4], RICEWQ [5], and the nitrogen balance model [6]. Some researchers developed watershed scale models to simulate water and pollutant loading from mixed watershed including paddy rice field. Takeuchi et al. [7] developed a cell-based distributed hydro-environmental watershed model for water budget analysis with three-zone cell profiling which consists of a surface water zone with the modified tank model, a surface soil zone with the soil moisture model, and a groundwater zone with the unconfined shallow groundwater model. Soil and Water Assessment Tool (SWAT) was modified and has been applied for simulation in paddy rice fields for the watershed scale [8,9,10]. Jeon et al. [11] developed the HSPF-Paddy model to evaluate the effect of water and nutrient loading from paddy rice fields within complex watershed on the river environment by modifying Hydrological Simulation Program-Fortran (HSPF), and evaluated the developed model for the watershed scale as well as the field scale model with high model efficiency.

The important steps in modeling are calibration and validation processes. Calibration involves minimization of the difference between observed and simulated values by adjusting model parameters [12], while validation involves the use of the second set of independent information to ensure the calibration accuracy [13]. Although newer versions of models continue to be released, much of recent research activity in model community is focused on automatic calibration by optimization techniques. The major advantage of optimization technique is to get better results without consuming significant efforts and time. Soil and Water Assessment Tool (SWAT) is currently linked with optimization technique, and has been used for automatic calibration [14,15,16]. Model-Independent Parameter Estimation and Uncertainty Analysis (PEST) is a nonlinear parameter estimation package, and widely used for automatic calibration tool of HSPF model [17,18,19,20].

The application of the model, developed to simulate a paddy rice field, to the observed data has been low due to the limited access for the model for the model to the public. In addition, there are a few study cases on characterizing water budget component of paddy rice field by comparing that of other land use types. In this study, HSPF-Paddy linked with PEST was calibrated and validated during the period between 1985 and 2001. The results were utilized to analyze the water budget components by land use type, and to characterize the water budget in paddy rice fields under an Asian monsoon climate compared with that under a Mediterranean climate.

2. Materials and Methods

2.1. Overview of HSPF-Paddy

The schematic diagram of water budget components in HSPF-Paddy model is shown in Figure 1 and the equation is as follows [11]:

STRi = STRi-1 + PREi + IRRi + EXTi - EVPi - RNFi - PRCi
where, PRE, IRR, EXT are the inflow from precipitation, irrigation, and external inflow, respectively; EVP is total actual evapotranspiration from interception, storage, upper and low zones, and groundwater; RNF is the outflow from surface, interflow, and groundwater; STR is the amount of storage in interception and surface by the dike height in paddy rice field; PRC is the deep percolation; The subscript i denotes day.

The input and output time-series data such as meteorological data are stored in Watershed Data Management Utility (WDMUtil). Dike height and irrigation was added to the HSPF-Paddy model as WDM time-series. Surface runoff occurred when the water depth in paddy rice fields is greater than the dike height. The infiltration rate can be applied separately during growing and non-growing seasons using Special Action block in the HSPF-Paddy input file.

Figure 1. Schematic diagram of the water budget components in Hydrological Simulation Program-Fortran (HSPF)-Paddy model [11].
Figure 1. Schematic diagram of the water budget components in Hydrological Simulation Program-Fortran (HSPF)-Paddy model [11].
Water 06 02041 g001 1024

2.2. Overview of PEST

Model-Independent Parameters Estimator (PEST) is the software package for parameter estimation of complex environmental computer models. The time series processor (TSPROC), which is the surface utility of PEST, was developed to automatically calibrate HSPF. The TSPROC can manage the WDM time-series of observed and simulated values. Figure 2 shows the interrelation between the HSPF, WDM time-series, and TSPROC. HSPF reads meteorological data in WDM file, run, and store the output in WDM file. The TSPROC reads simulated and observed values from WDM file, parameter ranges, and generates PEST input file (*.pst). The PEST generates random numbers for calibration parameters, run HSPF, and calculates objective function.

Figure 2. Overview interaction of The time series processor (TSPROC), Parameter Estimation and Uncertainty Analysis (PEST), HSPF, and Watershed Data Management (WDM) file.
Figure 2. Overview interaction of The time series processor (TSPROC), Parameter Estimation and Uncertainty Analysis (PEST), HSPF, and Watershed Data Management (WDM) file.
Water 06 02041 g002 1024

The PEST employs the Gauss-Marquardt-Levenberg method which requires a continuous relationship between model parameters and model output. This method helps find the optimized parameter in fewer model runs than any other parameter estimation methods. The PEST minimizes a single weighted least-squares objective function, Φ:

Water 06 02041 i001
where, B is a vector containing values of the calibrated parameters; n is the number of observation; xi is the ith observation; xi′ (B) is the simulated value corresponding to the ith observation; yi is the jth prior estimate; yi′(B) is the jth simulated value; wi is the weight for the ith observation; and wj is the weight for the jth prior estimate.

2.3. Study Area

The study area is the Bochung watershed, Korea, where Bochung river crosses the watershed (Figure 3). The total area is 553.56 km2 and average elevation and slope of watershed are 263.93 m and 32.09%, respectively. The length of total streams and main stream are 1689.70 and 68.05 km, respectively. The last stream order is 7th and total number of 1st stream order is 2592. The watershed is a traditional rural area. The major land use type is forest (65%) and the urban area accounts for about 3% of the total area (Table 1). The cropped area occupies 28%, including 16% of paddy rice fields.

Figure 3. Study area.
Figure 3. Study area.
Water 06 02041 g003 1024
Table 1. Land use distribution in study area.
Table 1. Land use distribution in study area.
Land UseUrbanCropPaddyForestPastureWaterBarren
Area (km2)16.164.384.6356.16.810.86.5
Percentage (%)3121665121

2.4. Modeling Approach

The land cover data which identifies sixteen possible land use types and digital elevation model (DEM) with 30 m cell size were achieved from Environmental Geographic Information System as GIS file [21]. Hourly meteorological data, which contains air temperature, dew point temperature, cloud cover, wind speed, precipitation, and solar radiation, were obtained from Korea Meteorological Administration (KMA). Evapotranspiration and evaporation were calculated by the WDM Utility. Daily stream flow data at the mouth of watershed were collected from the Water Management Information System (WAMIS) [22].

Pre-processors of HSPF-Paddy including watershed segmentation, land cover overlay, and input file generation were performed using Better Assessment Science Integrating Point and Nonpoint Source (BASINS). Using DEM and BASINS Delination Utility, the Bochung watershed was divided into 21 sub-watersheds automatically. Using land cover and BASINS Land Use Definition Utility, a land use division was performed for each sub-watershed. The pervious area of the urban was set 10% of total urban area. Jeon et al. [23] studied regionalized CN parameters at Nakdong River Basin, Korea using eight watersheds, and reported that the hydrologic characteristics of the urban area in Korea represented 90% of the impervious area.

BASINS technical note 6 documented by USEPA guides researchers in estimating hydrology and hydraulic parameters for HSPF [24]. The document provides “typical” and “possible” ranges for the HSPF parameters based on experiences with HSPF over the past two decades on watershed across the USA and abroad to be realistic parameters and reflect conditions on the watershed. The ranges of calibrated parameters are listed in Table 2 by referring BASINS technical note 6 [24]. The calibration period is for 10-year (1992~2001), and the validation period is for 7 years (1985~1991). The PEST program searched optimized parameters within the range for each land use type.

Table 2. Calibration parameters and possible range of value.
Table 2. Calibration parameters and possible range of value.
ParameterDescriptionPossible Range
LZSNLower zone nominal soil moisture storage5–38 cm
INFILTRelated to infiltration capacity of the soil0.25–25 mm/h
KVARYGroundwater recession flow parameter0–7.62 cm–1
AGWRCGroundwater recession rate0.833–0.999 day–1
DEEPFRFraction of groundwater inflow to deep recharge0.0–0.5
BASETPFraction of remaining ET from baseflow0.0–0.05
CEPSCInterception storage capacity0.25–10.1 mm
UZSNUpper zone nominal soil moisture storage0.1–5.0 cm
NSURManning’s n for overland flow0.05–0.5
INTFWInterflow inflow parameter1.0–10.0
IRCInterflow recession parameter0.3–0.85
LZETPLower zone ET parameter0.1–0.9

The average values of parameters varied with the land cover, such as LZSN, INFILT, CEPSC, UZSN, NSUR, and LZETP, were calculated separately according to the land use type referred by BASINS technical note [24]. Once the average values were determined for each land use type, maximum and minimum values are generated with the consideration of possible ranges shown in Table 2. The PEST program searched optimized values for calibration parameters within maximum and minimum values, and the ranges are shown in Table 3.

Table 3. Parameter ranges for land use type.
Table 3. Parameter ranges for land use type.
Land useLZSN (cm)INFILT (mm/h)CEPSC (mm)UZSN (cm)NSUR (none)LZETP (none)
Urban*3.8–18.82.7–18.00.00.1–3.00.010–0.140.3–0.8
Paddy4.1–20.30.127–0.1780.0–0.50.2–4.1<0.0100.8–2.5
Upland4.1–20.32.0–13.00.0–0.80.2–4.10.025–0.340.6–2.0
Pasture3.8–18.82.0–13.60.0–0.50.2–4.10.025–0.340.5–1.7
Forest4.1–20.31.4–9.20.0–1.00.2–5.10.037–0.500.7–2.3
Barren3.8–18.83.5–23.40.00.1–3.00.02–0.270.3–0.8

Note: * apply at pervious area of urban.

Jeon et al. [25] evaluated calibration methods of stream flow using HSPF coupled with PEST program, and proposed the calibration based on monthly flow is more effective than that on daily flow for the yearly basis water budget analysis. Therefore, objective function for automatic calibration was performed by monthly stream flow in this study. The Nash-Sutcliff Efficiency (NSE) [26] has been widely used as a “goodness-of-fit” statistics indicator to access the goodness of model to simulate the flow [27]. Ahmadi et al. [28] evaluate various statistical indices and concluded that root mean square (RMS) error, residual mass (RM) coefficient, and Nash-Sutcliff (NS) coefficient were better for the model calibration. In this study, RMS, NS, R2, and NS coefficients were used for calibration performance and these are as follows:

Water 06 02041 i002
Water 06 02041 i003
Water 06 02041 i004
Water 06 02041 i005
where, O is the mean of observed values; Pi is the predicted value; Oi is the observed value; and n is the number of data. The values for RMS, RM, NS, and R2 are 0.0, 0.0, 1.0, and 1.0, respectively, when observed and simulated vales are same.

3. Results and Discussion

3.1. Calibration and Validation

Graphical comparison of yearly runoff between simulated and observed values is displayed in Figure 4. The mean annual runoffs, simulated and observed, are 654.5 and 719.2 mm for calibration and 685.6 and 700.5 mm for validation, respectively (Table 4). Donigian [29] listed general calibration/validation tolerances or target and the ranges should be applied to annual mean values as shown in Table 5. The mean annual volumes are within the 10% difference for the calibration and validation results, showing a very good agreement.

The graphical comparison of simulated and observed monthly flow and statistical analysis are depicted in Figure 5 and Figure 6, respectively, and listed in Table 6. The hydrology results reflected the observed data well (Figure 5). The monthly NS values for calibration and validation are greater than 0.80 and monthly residual mass (RM) coefficient are consistently less than 0.10. Donigian [30] provided value ranges for coefficient of determination (R2) for assessing model performance for both daily and monthly flows. As the criteria of model performance and Table 5, the HSPF-Paddy model performances illustrated very good agreement (R2 > 0.85 and RM < 0.1) for monthly calibration and validation. Overall, the hydrology results show a very good agreement based on annual and monthly comparisons by using statistical and graphical analysis, and can be used for water budget analysis.

Figure 4. Comparison between observed and simulated yearly runoff.
Figure 4. Comparison between observed and simulated yearly runoff.
Water 06 02041 g004 1024
Table 4. Annual mean values and percent mean error for calibration and validation.
Table 4. Annual mean values and percent mean error for calibration and validation.
Calibration/validationObserved (mm)Simulated (mm)Percent Mean Error (%)
Calibration654.5719.29.9
Validation685.6700.52.2
Table 5. General calibration/validation target or tolerances for HSPF application.
Table 5. General calibration/validation target or tolerances for HSPF application.
ItemVery GoodGoodFair
Hydrology/flow<0.100.10–0.150.15–0.25
Sediment<0.200.20–0.300.30–0.45
Water quality/nutrients<0.150.15–0.250.25–0.35
Figure 5. Monthly simulated and observed runoff.
Figure 5. Monthly simulated and observed runoff.
Water 06 02041 g005 1024
Figure 6. Scatter plots of monthly simulated and observed runoff for calibration (a) and validation (b).
Figure 6. Scatter plots of monthly simulated and observed runoff for calibration (a) and validation (b).
Water 06 02041 g006 1024
Table 6. Statistical analysis for model performance.
Table 6. Statistical analysis for model performance.
Calibration/validation RMSRMNSR2
Calibration0.670.0890.810.89
Validation0.550.0050.840.87

3.2. Optimized Parameters

Optimized HSPF parameters for each land use type are listed in Table 7. Vegetated land cover which includes paddy, crop, and forest areas can retain rainfall with high CEPSCS, UZSN, and LZSN values compared to non-vegetated area. Most of parameters are optimized by PEST within the possible range recommended by USEPA [24]. The BASETP and INTFW values are optimized as upper and lower limit values, respectively. However, DEEPFR was not optimized within the possible range, showing 0.80, as the percent mean error for annual mean values couldn’t meet the criteria (Table 3) when PEST searched optimal DEEPFR value within the possible range. Some parameters of paddy rice field are outlier. The values for INFILT and NSUR are optimized under the lower limit values, indicating that, although the dike of paddy rice field can effectively reduce surface runoff by retaining rainfall, peak flow is higher and time for peak is shorter than other land use types once the surface runoff occurred.

3.3. Comparison of Water Budget Components

The results of water budget analysis during 1985–2001 are shown in Table 8. The influx for all land use types except paddy rice field is 1305 mm, contributed solely by rainfall, while that for paddy rice field is about 2695 mm with the addition of 1390 mm of irrigation. Total runoffs are ranged from 424 to 1697 mm, and are significantly influenced by land use types. Whereas runoff values for the pervious areas, such as crop, forest, pasture, and barren, show similar levels, the value for paddy rice field is higher than that for the urban area. This resulted from the low infiltration rate and manning’s coefficient despite the large irrigated volume.

Table 7. Optimized hydrological parameters of HSPF-Paddy model.
Table 7. Optimized hydrological parameters of HSPF-Paddy model.
ParameterUrbanCropPaddyForestPastureBarren
LZSN (cm)5.0020.3215.7515.855.005.00
INFILT (mm/h)1.783.430.132.272.2511.68
KVARY (1/cm)0.210.210.210.210.210.21
AGWRC (1/day)1.001.001.001.001.001.00
DEEPFR0.800.800.800.800.800.80
BASETP0.050.050.050.050.050.05
CEPSC (mm)0.003.962.415.112.740.00
UZSN (cm)0.130.190.200.250.210.22
NSUR0.040.210.010.380.340.10
INTFW0.840.840.840.840.840.84
IRC (1/day)0.700.700.700.700.700.70
LZETP0.330.690.750.690.660.04
Table 8. Comparison of water budget components by land use for simulation period (mm).
Table 8. Comparison of water budget components by land use for simulation period (mm).
ComponentUrbanCropPaddyForestPastureBarren
Influx
Rainfall130513051305130513051305
Irrigation001390000
Total130513052695130513051305
Runoff
Surface10972451251302354159
Interflow22132401127146242
Baseflow446453436114
Total1123 (86%)424 (32%)1697 (63%)463 (35%)536 (41%)515 (39%)
Groundwater
Inactive GW23 (2%)267 (20%)272 (10%)206 (16%)218 (17%)578 (44%)
Evapotranspiration
Intercept01971472271600
Upper zone16574406781141
Lower zone233312632429136
Baseflow22123181932
Total159 (12%)608 (47%)736 (27%)635 (49%)550 (42%)208 (16%)

The rates of each runoff component are shown in Figure 7. Most of the runoff is generated from the surface except the barren area. The surface runoff rate of the urban area is about 98% due to the large impervious surface area. The rate of surface runoff from the vegetation area such as crop, paddy, forest, and pasture ranged from 58% to 74%, and the highest rate of 74% is from paddy rice field. The interflow runoff from the urban area is lowest accounting only 2% of total runoff, and that from the vegetation area ranged from 24% to 31% of total runoff. The highest rate of interflow runoff is from barren areas showing about 47% of total runoff. The runoff from base flow is not significant compared to other runoff sources. The rate of base flow runoff from the vegetation area is less than 11% and that from barren is highest with 22% of total runoff. Inflow to deep ground water from the urban area is just 23 mm and those from vegetated area are ranged from 206 to 272 mm. Inflow to deep ground water from barren is about 578 mm.

Figure 7. Comparison of rates of runoff components by land use type.
Figure 7. Comparison of rates of runoff components by land use type.
Water 06 02041 g007 1024

Hydrologic characteristics of barren are somewhat different with other pervious area and are summarized as less surface runoff and more interflow and base flow runoff. Those results are obviously influenced by the hydrologic soil condition as shown in Table 9. Soil can be divided into four groups according to water transmitting soil layers with the lowest saturated hydraulic conductivity. Hydrologic soil groups A, B, C, and D are characterized by low, moderately low, moderately high, and high runoff potential, respectively [31]. Most of hydrologic soil group of pervious area is between B and C group however the area rate of hydrologic soil A of barren area is about 60% at the Bochung watershed. The calibrated infiltration rate of the barren area is highest with 11.68 mm/h, the high infiltration rate might reduce surface runoff increase interflow and base flow runoff.

Table 9. The comparison of hydrologic soil conditions with land use type at study area.
Table 9. The comparison of hydrologic soil conditions with land use type at study area.
Land useCalibrated infiltration rate (mm/h)Hydrologic soil group (HSG)
ABCD
Urban1.783149173
Crop3.431261225
Paddy0.132638315
Forest2.272592118
Pasture2.251365147
Barren11.685916719

The rates of each evapotranspiration component are shown in Figure 8. The evapotranspiration is significantly influenced by land use type and ranged from 159 to 736 mm. The high impervious surface of the urban area makes to less evapotransrate compared with other land use type. The amounts of evapotranspiration from vegetated area are more than 500 mm/yr and those from paddy rice field are 736 mm/yr showing the highest value among the land use type compared. More than 50% of total evapotranspiration from vegetated area except paddy rice field occurred from lower zone storage within soil, while about 60% of total evaporanspiration from paddy rice fields are generated from the upper zone resulting from ponded water remained during growing season (Figure 8).

Figure 8. Comparison of rates of evapotranspiration components by land use type.
Figure 8. Comparison of rates of evapotranspiration components by land use type.
Water 06 02041 g008 1024

The water budget analysis of paddy rice fields in various countries was shown in Table 10. In Asian monsoon regions including Taiwan [32], Korea [33], and in this study, around 50% of total input water came from precipitation. More than 90% of total input water was feed from irrigation and precipitation was less than 10% in Mediterranean regions including Spain [34] and Greece [6] because most precipitation is concentrated during the winter~spring season. Most of the input water in paddy rice fields was lost by runoff in an Asian monsoon area whereas deep percolation was observed in a Mediterranean climate area. Overall, much of the water entered paddy rice fields by precipitation and was lost by runoff in paddy rice fields in an Asian monsoon region.

Table 10. Comparison of water budget components in paddy rice fields by climate (mm/yr).
Table 10. Comparison of water budget components in paddy rice fields by climate (mm/yr).
ComponentMediterranean climate regionAsian monsoon region
SpainGreeceTaiwanKoreaThis study
Input20243896138020502695
Irrigation1874 (93%)3647 (94%) 624 (45%) 695 (34%)1390 (52%)
Precipitation150 (7%)249 (6%)756 (55%)1355 (66%)1305 (48%)
Output2024 3915138021432705
Evapotranspotation731 (36%)930 (24%)503 (36%)663 (31%)736 (27%)
Runoff372 (18%)1118 (29%)541 (39%)1352 (63%)1652 (61%)
Deep percolation830 (41%)1867 (48%)336 (24%)128 (6%)317 (12%)
Other91 (4%)

The calibration results represented in this study at one particular site may have more uncertainties [35]; thus, further verification at various sites and water budget components may be necessary.

4. Conclusions

The HSPF-Paddy model coupled with PEST was applied to the Bochung watershed in Korea to analyze the characteristics of water budget components in paddy rice field under an Asian monsoon climate during 1985–2001. Water movement within the watershed was well simulated by HSPF-Paddy showing high model efficiency, and simulation results met the criteria of calibration and validation with the acceptable error level. The hydrological characteristics of paddy rice fields were somewhat different to those of other vegetated areas. Runoff from paddy rice fields was much higher than that from other land use types compared due to the low infiltration rate and manning’s coefficient as well as the irrigation, implying that most runoff is generated from the surface zone in paddy rice fields. Another characteristic of paddy rice fields is a significant amount of evapotranspiration compared with other vegetated areas resulting from the ponded water maintained during the growing season. The characteristics of climate influenced the water budget in paddy rice fileds. In an Asian monsoon area, much precipitation entered the paddy rice field and significant water was lost by runoff because most precipitation is concentrated on the summer season in contrast to a Mediterranean climate area.

Acknowledgments

This subject is supported by Korea Ministry of Environment as “Global Top Project” (Project No.: GT-11-B-02-004-1).

Author Contributions

Young-Jin Kim wrote the paper and applied HSPF-Paddy linked with Model-Independent Parameters Estimator (PEST) program. Hae-Do Kim collaborated on calibration. Ji-Hong Jeon reviewed the model application and directed this research.

Conflicts of Interests

The authors declare no conflict of interest.

References

  1. Takase, K. Japan, Monsoon Asia and the global agricultural revolution in the 21st century. Paddy Water Environ. 2003, 1, 13–18. [Google Scholar] [CrossRef]
  2. Yoon, G.G.; Ham, J.-H.; Jeon, J.-H. Mass balance analysis in Korean paddy rice culture. Paddy Water Environ. 2003, 1, 99–106. [Google Scholar] [CrossRef]
  3. Jeon, J.-H.; Yoon, C.G.; Han, J.-H.; Jung, K.-W. Model development for surface drainage loading estimates from paddy rice fields. Paddy Water Environ. 2005, 3, 93–101. [Google Scholar] [CrossRef]
  4. Chung, S.O.; Kim, H.S.; Kim, J.S. Model development for nutrient loading from paddy rice fields. Agric. Water Manag. 2003, 62, 1–17. [Google Scholar] [CrossRef]
  5. Williams, M.W.; Ritter, A.M.; Cheplick, M.J.; Zdinak, C.E. RICEWQ: Pesticide runoff model for rice crops. In User Manual and Program Documentation; Version 1.6.1; Waterborne Environmental, Inc.: Leesburg, VA, USA, 1999. [Google Scholar]
  6. Antonopoulos, V.Z. Modeling of water and nitrogen balance in the ponded water of rice fields. Paddy Water Environ. 2008, 6, 387–395. [Google Scholar] [CrossRef]
  7. Takeuchi, J.; Kawachi, T.; Unami, K.; Meada, S.; Izumi, T. A distributed hydro-environmental watershed model with three-zoned cell. Paddy Water Environ. 2009, 7, 33–43. [Google Scholar] [CrossRef]
  8. Sakaguchi, A.; Eguchi, S.; Kato, T.; Kasuya, M.; Ono, K.; Miyata, A.; Tase, N. Development and evaluation of paddy module for improving hydrological simulation in SWAT. Agric. Water Manag. 2014, 137, 116–122. [Google Scholar] [CrossRef]
  9. Jung, C.-G.; Park, J.-Y.; Kim, S.-J.; Park, G.-A. The SRI (system of rice intensification) water management evaluation by SWAPP (SWAT-APEX Program) modeling in an agricultural watershed of South Korea. Paddy Water Environ. 2014, 12, 251–261. [Google Scholar] [CrossRef]
  10. Xie, X.; Cui, Y. Development and test of SWAT for modeling hydrological processes in irrigation districts with paddy rice. J. Hydrol. 2011, 1–2, 61–71. [Google Scholar] [CrossRef]
  11. Jeon, J.-H.; Yoon, C.G.; Donigian, A.S.; Jung, K.-W. Development of the HSPF-Paddy model to estimate watershed pollutant loads in paddy farming regions. Agric. Water Manag. 2007, 90, 75–86. [Google Scholar] [CrossRef]
  12. Jewell, T.K.; Nunno, T.J.; Adrian, D.D. Methodology for calibrating stormwater model. J. Environ. Eng. Div. 1978, 104, 485–501. [Google Scholar]
  13. U.S. Environmental Protection Agency (USEPA). Compendium of Tool for Watershed Assessment and TMDL Development; EPA841-B-97-006; USEPA, Office of Water: Washington, DC, USA, 1997. [Google Scholar]
  14. Shawul, A.A.; Alamirew, T.; Dinka, M.O. Calibration and validation of SWAT model and estimation of water balance components of Shaya mountainous watershed, Southeastern Ethiopia. Hydrol. Earth Syst. Sci. Discuss. 2013, 10, 13955–13978. [Google Scholar] [CrossRef]
  15. Wu, Y.; Liu, S.; Li, Z.; Dahal, D.; Young, C.J.; Schmidt, G.L.; Liu, J.; Davis, B.; Sohl, T.L.; Wemer, J.M.; et al. Development of a generic auto-calibration package for regional ecological modeling and application in the Central Plains of the United States. Ecol. Inform. 2014, 19, 35–46. [Google Scholar] [CrossRef]
  16. Wu, Y.; Liu, S. Improvement of the R-SWAT-FME framework to support multiple variables and multi-objective functions. Sci. Total Environ. 2014, 466–467, 455–466. [Google Scholar] [CrossRef]
  17. Kim, J.; Ryu, J. A threshold of basin discretization level for HSPF simulation with NEXRAD input. J. Hydrol. Eng. 2013. [Google Scholar] [CrossRef]
  18. Diaz-Ramirez, J.N.; Johnson, B.E.; McAnally, W.H.; Martin, J.L.; Alarcon, V.J. Estimation and propagation of parameter uncertainty in lumped hydrological models: A case study of HSPF model applied to Luxapallila Creek Watershed in Southeast USA. J. Hydrogeol. Hydrol. Eng. 2013. [Google Scholar] [CrossRef]
  19. Kim, S.M.; Benham, B.L.; Brannan, K.M.; Zeckoski, R.W.; Doherty, J. Comparison of hydrological calibration of HSPF using automatic and manual method. Water Resour. Res. 2007, 43. [Google Scholar] [CrossRef]
  20. Jeon, J.-H.; Lee, S.-B. Automatic calibration for noncontinuous observed data using HSPF-PEST. J. KoreanSoc. Agric. Eng. 2012, 54, 111–119. (In Korean). [Google Scholar]
  21. Environmental Geographic Information System. Available online: http://egis.me.go.kr/main.do (accessed on 13 May 2014).
  22. Water Management Information System. Available online: http://www.wamis.go.kr/ (accessed on 13 May 2014). (In Korean).
  23. Jeon, J.-H.; Choi, D.H.; Kim, J.J.; Kim, T.-D. Regionalization of CN parameters for Nakdong Basin using SCE-UA algorithm. J. Korean Soc. Water Qual. 2009, 25, 245–255. (In Korean). [Google Scholar]
  24. U.S. Environmental Protection Agency (USEPA). BASINS Technical Note 6: Estimating Hydrology and Hydraulic Parameters for HSPF; EPA-823-R00-012; USEPA, Office of Water: Washington, DC, USA, 2000. [Google Scholar]
  25. Jeon, J.-H.; Choi, D.; Kim, J.-J.; Kim, T.-D. Evaluating calibration methods of stream flow for water quality management. J. Korean Soc. Agric. Eng. 2009, 25, 432–440. (In Korean). [Google Scholar]
  26. Nash, J.E.; Sutcliffe, J.V. River flow forecasting through conceptual model, Part I: A discussion of principles. J. Hydrol. 1970, 10, 282–290. [Google Scholar] [CrossRef]
  27. Shrestha, N.K.; Shakti, P.C.; Gurung, P. Calibration and validation of SWAT model for Low Lying Watersheds: A case study on the Kliene Nete Watershed, Belgium. Hydrol. Nepal. J. Water Energy Environ. 2010, 6, 47–51. [Google Scholar]
  28. Ahmadi, S.H.; Amin, S.; Keshavarzi, A.R.; Mirzamostafa, N. Simulating watershed outlet sediment concentration using the ANSWERS model by applying two sediment transport capacity equation. Biosyst. Eng. 2006, 94, 615–626. [Google Scholar] [CrossRef]
  29. Donigian, A.S., Jr. Slide #L19–22; Washington, DC, USA, 10–14 January 2000, U.S. EPA Office of Water, Office of Science and Technology: Washington, DC, USA, 2000.
  30. Donigian, A.S. Watershed Model Calibration and Validation: The HSPF Experience, National TMDL Science and Policy Specialty Conference 2002; Water Environment Federation: Phoenix, AZ, USA, 2000. [Google Scholar]
  31. United States Department of Agriculture, Soil Conservation Service (USDA-SCS). National Engineering Handbook Chapter: Chapter 7: Hydrologic Soil Groups; USDA-SCS: Washington, DC, USA, 2004. [Google Scholar]
  32. Chien, C.-P.; Fang, W.-T. Modeling irrigation return flow for the return flow reuse system in paddy fields. Paddy Water Environ. 2012, 10, 187–196. [Google Scholar] [CrossRef]
  33. Song, J.-H.; Kang, M.-S.; Song, I.-H.; Jang, J.-R. Comparing farming methods in pollutant runoff loads from paddy fields using the CREAMS-PADDY model. Korean J. Environ. Agric. 2012, 31, 318–327. [Google Scholar] [CrossRef]
  34. Playán, E.; Pérez-Coveta, O.; Martínez-Cob, A.; Herrero, J.; García-Navarro, P.; Latorre, B.; Brufau, P.; Garcés, J. Overland water and salt flows in a set of rice paddies. Agric. Water Manag. 2008, 95, 645–658. [Google Scholar] [CrossRef]
  35. Wang, S.; Huang, J.; Li, J.; Rivera, A.; McKenney, D.W.; Sheffield, J. Assessment of water budget for sixteen large drainage basins in Canada. J. Hydrol. 2014, 512, 1–15. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Kim, Y.-J.; Kim, H.-D.; Jeon, J.-H. Characteristics of Water Budget Components in Paddy Rice Field under the Asian Monsoon Climate: Application of HSPF-Paddy Model. Water 2014, 6, 2041-2055. https://0-doi-org.brum.beds.ac.uk/10.3390/w6072041

AMA Style

Kim Y-J, Kim H-D, Jeon J-H. Characteristics of Water Budget Components in Paddy Rice Field under the Asian Monsoon Climate: Application of HSPF-Paddy Model. Water. 2014; 6(7):2041-2055. https://0-doi-org.brum.beds.ac.uk/10.3390/w6072041

Chicago/Turabian Style

Kim, Young-Jin, Hae-Do Kim, and Ji-Hong Jeon. 2014. "Characteristics of Water Budget Components in Paddy Rice Field under the Asian Monsoon Climate: Application of HSPF-Paddy Model" Water 6, no. 7: 2041-2055. https://0-doi-org.brum.beds.ac.uk/10.3390/w6072041

Article Metrics

Back to TopTop