Scale Effects of Water Saving on Irrigation Efficiency: Case Study of a Rice-Based Groundwater Irrigation System on the Sanjiang Plain, Northeast China

This research analyzed the scale effect of water saving in Bielahonghe (BLH) Basin, a rice-cultivating district on the Sanjiang Plain, Northeast China. Water budgets with different surface irrigation water supply ratios and water-saving measures were simulated with a semi-distributed water balance model. PFnws, representing the ratio of rice evapotranspiration to net water supply (the total amount of irrigation and precipitation minus the amount of water reused), was employed to assess the water use efficiency. Seven spatial scales (noted from S1 to S7), ranging from a single field (317.87 ha) to the whole basin (about 100,800 ha) were determined. PFnws values were quantified across scales and several water-saving measures, including water-saving irrigation regimes, canal lining, and a reduction of the surface water supply ratio (SWSR). The results indicated that PFnws increased with scale and could be calculated by a fitted power function (PFnws = 0.736Area0.033, R2 = 0.58). Furthermore, PFnws increased most prominently when the scale increased from S1 to S2. The water-saving irrigation regime (WSIR) had the most substantial water-saving effect (WSE) at S1. Specifically, PFnws improved by 21.2% at S1 when high-intensity WSIR was applied. Additionally, the WSE values of S3 and S5 were slightly higher than at other scales when the branch canal water delivery coefficient increased from 0.65 to 0.80 through canal lining. Furthermore, the PFnws at each scale varied with SWSR. Specifically, PFnws from S3 to S7 improved as SWSR decreased from 0.4 to 0.3 but remained approximately constant when SWSR decreased from 0.3 to 0.


Introduction
As a primary component of water consumption, irrigation is critical to global food production.Growing water scarcity increases the need for water savings, which relies on precise assessments of irrigation efficiency [1,2].Irrigation efficiency and water savings have become increasingly significant factors in effective agricultural water management [3].A series of irrigation efficiency indices have been proposed to demonstrate the level of agricultural water management [4][5][6][7][8].The classic concept of irrigation efficiency (PF gws ) is defined as the ratio of crop evapotranspiration (ET crop ) to the total Sustainability 2018, 10, 47; doi:10.3390/su10010047www.mdpi.com/journal/sustainabilityamount of water inputs (sum of precipitation and irrigation) [9,10].As it is well known, the applied irrigation (or used irrigation water) ends up as: (1) beneficial evapotranspiration; (2) non-beneficial evapotranspiration; (3) non-recoverable runoff/percolation (e.g., by agriculture, though useful for alternative uses such as the environment); and (4) recoverable runoff/percolation [11].The first three components constitute the consumed or depleted fraction of irrigation water used, meaning that this water is not available for further use as it is consumed as evapotranspiration, incorporated into a product, or flows to a location where it cannot be readily reused (e.g., heavily saline water).However, the fourth component of the applied irrigation is not consumed and is recoverable for further/later abstractions.Surface drainage and percolation, which flow from the fields into the drains and groundwater, account for a large proportion of the water inputs [12][13][14].Water-saving is one of the main objectives to improve irrigation efficiency [15].Water-saving measures such as drip and sprinkler irrigation, deficit irrigation, and canal lining have been developed and applied in the field to reduce seepage and percolation flows, providing considerable increases in irrigation efficiency [16].
However, irrigation fields are not isolated units but, rather, are part of a whole irrigation system.Water 'savings' in one place are likely to reduce return flows to other users (e.g., field) located downstream of the basin [17].From the perspective of a whole system, a certain amount of so-called lost water (deep percolation, drainage to ditches etc.) may not truly be lost [18].Outflows at a small scale, such as percolation and drainage, may be reused at a larger scale.'Losses' at the field level can be captured and reused downstream and do not necessarily lead to true water depletion at the irrigation system level.It is necessary to assess whether water 'saving' measures are real or only produce local 'savings'.The irrigation efficiencies in the whole system as well as within river basins are more crucial for agricultural water management [19].Therefore, water reuse and water-savings quantification at different scales are two important aspects of agricultural water management [20][21][22].
Bos and Wolters [23] found that the high recycling rate of small-scale 'water loss' (recoverable runoff/percolation) could improve the water use efficiency of an entire system.It has been argued that irrigation efficiency may increase with increasing spatial scale and may be much higher at the irrigation system level than at the individual field level [24,25].Because of the different ratios of water reuse, increasing irrigation efficiency through water-saving measures varies at different scales.In addition to the non-linear characteristics of the hydrological processes in the irrigation system, the effects of water savings on irrigation efficiency are scale dependent.Willardson [26] found that the influence of water use efficiency improvements on an individual plot is not significant for watershed-scale efficiency and that water use efficiency improvements at a small scale could lead to either improvements or reductions at a large scale.Molden [27] reported that small-scale and large-scale improvements of water use efficiency are not entirely consistent.Palanisami et al. [28] discovered that improvements in water use efficiency at a point scale do not necessarily indicate improvements in water use efficiency at a large scale.An inaccurate understanding of scale effects can cause confusion when implementing water-saving measures.It is not clear whether the implementation of water-saving irrigation measures at a field scale will produce expected water-saving effects at the basin scale.To select appropriate water-saving irrigation measures and water resource allocation strategies, the scale effects of water use efficiency should be quantitatively studied.
Two main issues must be addressed in such studies.First, conclusions based on the classic index of irrigation efficiency might be incorrect if boundaries at each scale are not correctly understood at the watershed level.Ward and Pulido-Velazquez [29] suggested that the water that is not evaporated and can be recovered should be considered in irrigation efficiency measurements.Keller and Keller [30] noted that effective outflow should be reduced from the water inputs when evaluating irrigation efficiency.However, reusable water is considered, rather than actual water reuse, limiting the quantification of the implementation effect of water-saving technologies.Second, the scale effects of water savings are complicated because the increase in irrigation efficiency should be precisely evaluated after the implementation of water-saving measures.
To evaluate irrigation efficiency and quantify the water savings at different scales, detailed water accounting at corresponding scales should be conducted in advance [31].The water balance model employs a clear concept and requires few parameters compared with hydrodynamic models such as HYDRUS [32], SWAP [33], and MIKE-SHE [34]; therefore, it is suitable for water budgets with a long time-step in a large region and has been widely used throughout the world [35][36][37][38].In an irrigation system, the hydrological processes are heavily affected by artificial measures such as crop cultivation, sluice and pump stations, artificial canals and ditches for irrigation and drainage.Due to complicated landscapes with abundant water connections among surface, unsaturated and saturated soil across scales [39], a semi-distributed hydrological model considering the aforementioned physical realities is important to account for the water budget, especially with respect to water reuse.The study area is usually divided into sub-regions to consider spatial variability.A lumped method was employed to create detailed descriptions of water movement in each sub-region which indicated same parameter values were used for specific sub-region.Hydraulic connections among sub-regions were established by identifying water exchanges in groundwater, canals and ditches [40,41].A semi-distributed physically-based model requires fewer parameters than a distributed model.Hence, it is relatively easy to add physical processes to describe artificial influences on the water cycle, and it is appropriate to analyze the scale effect and water budget in an irrigated area with a semi-distributed model.
In the context of increasing water scarcity, rice has become a main target for water-saving measures considering that 34-43% of the world's irrigation water is used for rice [42].Because of the percolation and surface seepage flows from continuously ponded fields into groundwater and drains, seepage and percolation flows account for a large proportion of the water inputs [43][44][45].The contributions of these flows are particularly high in a rice-based groundwater irrigation system, where a large proportion of groundwater may come from percolation flow.The Bielahonghe (BLH) Basin, located on the Sanjiang Plain, Northeast China, is a rice-based irrigation district that is supplied by groundwater with supplemental irrigation from drainage ditches.In this study, a semi-distributed water balance model was developed and used to account for the components of the water balance.An index of irrigation efficiency considering water reuse was implemented at various scales in the BLH Basin.The scale effects of irrigation efficiency were evaluated, and the contrast with the classic irrigation efficiency index was highlighted.Then, the water budget under different water-saving scenarios (surface irrigation water supply ratios and water-saving measures) was simulated.Water use efficiency and the water reuse fraction at different levels were analyzed, and the scale effects of different water-saving technologies were quantified.

Study Area and Data
Bielahonghe (BLH) Basin, located on the Sanjiang Plain, Northeast China, is a rice-based irrigation district with an area of 10.08 × 10 4 ha (Figure 1).The average elevation of this area is 40-50 m, with a slope of 1/5000-1/10,000 and the maximum and minimum elevation are 65 m and 39 m respectively.The study area has a temperate continental monsoon climate with a mean annual temperature of 2.2 • C; the freezing period lasts for up to 120-140 days, and the maximum frozen depth can reach 2.2 m.The average annual precipitation is 545.4 mm, 50-70% of which is concentrated from July to September.Cultivated land covers 69.8% of the study area, and rice is the only crop.Other land use types include forest, grass, bare soils and residential land.Field steeping usually starts in late April; the rice is transplanted in early May and is harvested at the end of September or early October.Groundwater is the dominant source of irrigation water, causing the groundwater table to decline by 37 cm per year.The Quaternary aquifer with a thickness of 180 m is composed of 8-14 m-thick clay loam at the top and sandy loam, sand, and gravel underneath.More exactly, thickness of clay loam of topsoil was higher in west and south when compared with east and north in our study area.Meanwhile, the change of thickness of sand aquifer in the subsoil was opposite.In addition, permeability coefficient changed from 50-60 m/d in west and south to 30-40 m/d in east and north, which was in consistent with the change of elevation.Meanwhile, the horizontal water exchange could be determined based on the monitoring of groundwater in and out of the study region.Drainage system in our study area include three levels of ditches (e.g., lateral ditch, main ditch) and the BLH River.The BLH River collects the drainage from all ditches in the area and flows from west to east to the outlet of the study area.Due to the irrigation channels are being planned, there is no irrigation channel in our study area presently and the main source for irrigation is groundwater.Although farmers might also use private pumps to lift water from the ditches for irrigation, the amount is quite small in our study.Water-saving scenarios have been planned by the authorities [46,47].
Sustainability 2018, 10, 0047 10.3390/su10010047 4 of 18 with the change of elevation.Meanwhile, the horizontal water exchange could be determined based on the monitoring of groundwater in and out of the study region.Drainage system in our study area include three levels of ditches (e.g., lateral ditch, main ditch) and the BLH River.The BLH River collects the drainage from all ditches in the area and flows from west to east to the outlet of the study area.Due to the irrigation channels are being planned, there is no irrigation channel in our study area presently and the main source for irrigation is groundwater.Although farmers might also use private pumps to lift water from the ditches for irrigation, the amount is quite small in our study.Watersaving scenarios have been planned by the authorities [46,47].The hydrological, meteorological, pumped irrigation and river flow data from 2010 to 2015 were collected from the Authority of the BLH Basin.(i) Daily meteorological data were collected from three meteorological stations (Figure 1).(ii) Groundwater table data were from 20 groundwater table observation wells in the study area (Figure 1).Fourteen (four original sets and ten new sets in August 2014) of these wells were located inside the basin, and six (original ones) were outside the basin.A manual observation method was adopted for the original observation wells with 10-day intervals.An automatic recording method was adopted for new observation wells, and the recording interval was 12 h.(iii) Data on pumped water for irrigation were collected by monitoring typical wells equipped with various motors using a flow-meter to obtain the pumped water quantity.The types and machine power of these typical wells could represent the general conditions for the whole region.Then, we recorded the pumping well distribution map and the horsepower of the equipped motors for the 7400 pumping wells in the study area.The irrigation time in different zones was obtained through the monitoring data from the four wells (Figure 1).Thus, the pumped water amount and irrigation time for all polygonal units (PUs) were determined by accumulating data of wells for specific PU and PU indicated the control area for one ditch (Details about PU were shown in Section 2.2). (iv) For data on river flow, three river flow and level observation stations were established along the BLH River (Figure 1).The water levels in the BLH River were automatically recorded with a water level recorder at 12-h intervals.The river discharge was measured with an LS45-2 current meter to obtain the stage-discharge rating curves.These curves were used to ascertain the river flow changes based on water level records.
The crop coefficients (Kc) were obtained from previous study and root depth of rice was assumed to be 30 cm [48].The positions of ditches and cross-section parameters were collected from The hydrological, meteorological, pumped irrigation and river flow data from 2010 to 2015 were collected from the Authority of the BLH Basin.(i) Daily meteorological data were collected from three meteorological stations (Figure 1).(ii) Groundwater table data were from 20 groundwater table observation wells in the study area (Figure 1).Fourteen (four original sets and ten new sets in August 2014) of these wells were located inside the basin, and six (original ones) were outside the basin.A manual observation method was adopted for the original observation wells with 10-day intervals.An automatic recording method was adopted for new observation wells, and the recording interval was 12 h.(iii) Data on pumped water for irrigation were collected by monitoring typical wells equipped with various motors using a flow-meter to obtain the pumped water quantity.The types and machine power of these typical wells could represent the general conditions for the whole region.Then, we recorded the pumping well distribution map and the horsepower of the equipped motors for the 7400 pumping wells in the study area.The irrigation time in different zones was obtained through the monitoring data from the four wells (Figure 1).Thus, the pumped water amount and irrigation time for all polygonal units (PUs) were determined by accumulating data of wells for specific PU and PU indicated the control area for one ditch (Details about PU were shown in Section 2.2). (iv) For data on river flow, three river flow and level observation stations were established along the BLH River (Figure 1).The water levels in the BLH River were automatically recorded with a water level recorder at 12-h intervals.The river discharge was measured with an LS45-2 current meter to obtain the stage-discharge rating curves.These curves were used to ascertain the river flow changes based on water level records.
The crop coefficients (Kc) were obtained from previous study and root depth of rice was assumed to be 30 cm [48].The positions of ditches and cross-section parameters were collected from the Qianjin (QJ) Water Administration Bureau.(v) The soil parameters were as follows: soil samples of different depths on 24 selected plots were obtained through drilling, and the soil particle-size distributions were measured.(vi) For hydrogeological parameters, a pumping test was conducted in four plots in October 2014 and March 2015, and the horizontal hydraulic conductivity and specific yields in the aquifer were calculated with test data.(vii) The leaf area index (LAI) at different growth stages were collected from the irrigation experiment station (30 km from the eastern boundary of the study area).In addition, due to the unified management, phenology data (e.g., sowing date, harvest date) and crop cultivar are the same for the study area.

Distributed Water Balance Model
The distributed water balance model is composed of three types of water balance units: PUs, ditch/river units (DUs), and canal units (CUs).A PU was the total area served by a branch ditch, and five land use types (Figure 2) were considered in a PU.The water budget under various land use types was first simulated; then, the amount of each water balance component in PUs was obtained through accumulation based on the land use area ratio.
Sustainability 2018, 10, 0047 10.3390/su10010047 5 of 18 the Qianjin (QJ) Water Administration Bureau.(v) The soil parameters were as follows: soil samples of different depths on 24 selected plots were obtained through drilling, and the soil particle-size distributions were measured.(vi) For hydrogeological parameters, a pumping test was conducted in four plots in October 2014 and March 2015, and the horizontal hydraulic conductivity and specific yields in the aquifer were calculated with test data.(vii) The leaf area index (LAI) at different growth stages were collected from the irrigation experiment station (30 km from the eastern boundary of the study area).In addition, due to the unified management, phenology data (e.g., sowing date, harvest date) and crop cultivar are the same for the study area.

Distributed Water Balance Model
The distributed water balance model is composed of three types of water balance units: PUs, ditch/river units (DUs), and canal units (CUs).A PU was the total area served by a branch ditch, and five land use types (Figure 2) were considered in a PU.The water budget under various land use types was first simulated; then, the amount of each water balance component in PUs was obtained through accumulation based on the land use area ratio.Water-transfer relationships among the PUs were established according the distribution of canals and ditches as well as groundwater movement direction (Figures 2 and 3).Canals (CUs) carried out stepwise distribution of water according to the water supply relationship between the upper and lower levels until the PUs were fed through the final-level CUs.After that, a part of PUs' water flowed into the nearby ditches (DUs) by overland flower, and part of PUs' water transported downward to the groundwater.Meanwhile, some of groundwater entered into DUs through subsurface drainage and some exchanged sideways with other PUs' groundwater by horizontal flow.The overland flow and groundwater drainage from PUs, as well as surplus water from canals, move to branch ditches, then to main ditches, then discharge into the BLH River, eventually flowing out of the study area.
For the water balance in CUs, the amount of rainfall, water diversion from upstream canals, water surface evaporation, water distribution to downstream canals and PUs, and seepage were considered.For the water balance in DUs and the BLH River, the amounts of rainfall, drainage source from up-level ditches and PUs, drainage sinks to down-level ditches, water surface evaporation, reused water, and water interaction with groundwater were involved [49,50].
The soil profile in PUs was divided into a ponding layer, rice root layer, soil water transmission layer, and groundwater layer.The soil water transmission layer was subdivided into a clay layer and sand layer according to the soil texture.After rainfall or irrigation, the water-ponding depth in PUs would increase and water would infiltrate into the rice root layer as the supply for crop evapotranspiration, resulting in changes in soil water content in the rice root layer.When the water Water-transfer relationships among the PUs were established according the distribution of canals and ditches as well as groundwater movement direction (Figures 2 and 3).Canals (CUs) carried out stepwise distribution of water according to the water supply relationship between the upper and lower levels until the PUs were fed through the final-level CUs.After that, a part of PUs' water flowed into the nearby ditches (DUs) by overland flower, and part of PUs' water transported downward to the groundwater.Meanwhile, some of groundwater entered into DUs through sub-surface drainage and some exchanged sideways with other PUs' groundwater by horizontal flow.The overland flow and groundwater drainage from PUs, as well as surplus water from canals, move to branch ditches, then to main ditches, then discharge into the BLH River, eventually flowing out of the study area.
For the water balance in CUs, the amount of rainfall, water diversion from upstream canals, water surface evaporation, water distribution to downstream canals and PUs, and seepage were considered.For the water balance in DUs and the BLH River, the amounts of rainfall, drainage source from up-level ditches and PUs, drainage sinks to down-level ditches, water surface evaporation, reused water, and water interaction with groundwater were involved [49,50].
The soil profile in PUs was divided into a ponding layer, rice root layer, soil water transmission layer, and groundwater layer.The soil water transmission layer was subdivided into a clay layer and sand layer according to the soil texture.After rainfall or irrigation, the water-ponding depth in PUs would increase and water would infiltrate into the rice root layer as the supply for crop evapotranspiration, resulting in changes in soil water content in the rice root layer.When the water content of each unsaturated soil layer exceeded the saturated water content of the corresponding soil layer, the excess water would enter the next layer to recharge the groundwater, causing the groundwater table to increase.Overland flow did not occur until the ponding water depth exceeded the height of the ridge (for non-paddy lands, the ridge height was set to zero).When the ponding water depth in the field (ponding water controlled condition) or the soil water content of the rice root layer (soil water content controlled condition) was below the lower limit of the control, irrigation was performed.The amount of irrigation water was determined by the gap between the upper limit and the lower limit of the ponding water depth (or of the soil water content in the rice root layer) at different growth stages in paddy rice.
Sustainability 2018, 10, 0047 10.3390/su10010047 6 of 18 content of each unsaturated soil layer exceeded the saturated water content of the corresponding soil layer, the excess water would enter the next layer to recharge the groundwater, causing the groundwater table to increase.Overland flow did not occur until the ponding water depth exceeded the height of the ridge (for non-paddy lands, the ridge height was set to zero).When the ponding water depth in the field (ponding water controlled condition) or the soil water content of the rice root layer (soil water content controlled condition) was below the lower limit of the control, irrigation was performed.The amount of irrigation water was determined by the gap between the upper limit and the lower limit of the ponding water depth (or of the soil water content in the rice root layer) at different growth stages in paddy rice.The main inputs of the model included the following: meteorological data, water diversion from each canal, irrigation by pumped wells in PUs, crop coefficients at different growth stages and the ratio of crop transpiration and evaporation distribution, control rules for the water layer in different irrigation scenarios, the surface elevation of each PU, the texture distribution and soil properties of the soil profile (saturated water content, field water holding capacity and vertical permeability coefficient), the specific yield and horizontal permeability coefficient of aquifers, the control area and land use ratio of each PU, the geometrical data (length, cross-sectional size) of the ditches and rivers at all levels, the available water utilization coefficient of each canal, the initial groundwater table and the maximum groundwater depth.
The study area was divided into 365 PUs, 365 branch drains and 365 branch canals, 28 main drains and 28 main canals (subdivided into 43 main drain segments and 43 main canal segments), and 1 river (subdivided into 6 river segments).The model was run on a daily scale.

Water Use Efficiency at Multiple Levels
Recently, the net water supply quantity (irrigation and rainfall minus the amount of reused water) was regarded as the denominator in calculating water use efficiency to reflect the real utilization efficiency of water resources.Therefore, water use efficiency was calculated as follows: The main inputs of the model included the following: meteorological data, water diversion from each canal, irrigation by pumped wells in PUs, crop coefficients at different growth stages and the ratio of crop transpiration and evaporation distribution, control rules for the water layer in different irrigation scenarios, the surface elevation of each PU, the texture distribution and soil properties of the soil profile (saturated water content, field water holding capacity and vertical permeability coefficient), the specific yield and horizontal permeability coefficient of aquifers, the control area and land use ratio of each PU, the geometrical data (length, cross-sectional size) of the ditches and rivers at all levels, the available water utilization coefficient of each canal, the initial groundwater table and the maximum groundwater depth.
The study area was divided into 365 PUs, 365 branch drains and 365 branch canals, 28 main drains and 28 main canals (subdivided into 43 main drain segments and 43 main canal segments), and 1 river (subdivided into 6 river segments).The model was run on a daily scale.

Water Use Efficiency at Multiple Levels
Recently, the net water supply quantity (irrigation and rainfall minus the amount of reused water) was regarded as the denominator in calculating water use efficiency to reflect the real utilization efficiency of water resources.Therefore, water use efficiency was calculated as follows: where λ is the amount of water reuse.In the BLH Basin, irrigation water came partly or wholly from the recycling of surface and subsurface return water (groundwater pumping and drainage reuse for irrigation).
The water reuse fraction (RF) and effective return water fraction (ERFF) were selected to quantify the magnitude of water reuse and outflow outside of the scale boundaries.
RF refers to the ratio of water reused to the theoretical return water (Q CRF ), which reflects the extent to which the return water was reused (Equation ( 2)).If there was no reused water in an irrigation system, the RF was zero.However, if all theoretical return water was reused in a closed irrigation system, the RF was 1. ERFF was defined as the fraction of actual outflow across scale boundaries to the theoretical return water (Equation ( 3)).A larger ERFF indicated more outflow and little reused water [51,52].
Irrigation and reused water differed in terms of the quantity and connotation at multiple levels (Table 1).Therefore, we should define the boundaries of scales and items of irrigation and reused water to determine the water use efficiency at different scales.In this study, seven scales, ranging from a single field to the whole basin were determined, and their horizontal boundaries are shown in Figure 4. Water reuse existed at all scales except S1.In addition, the upward water reuse in S3 and S7 were the smaller value between pump output and groundwater recharge.More exactly, when pump output was larger than groundwater recharge, we could declare that all groundwater recharge was reused.But in the opposite condition, only part of groundwater recharge (pump output) could be reused (Table 1).

Arrangements of Modeling Scenarios
Parameters including field water capacity, saturated water content, soil water specific yield, and permeability coefficient in our model were firstly calibrated by observed groundwater tables and river discharges from 14 wells and 3 rivers during the growth of paddy rice in 2015.The groundwater tables of 4 wells from 2010 to 2014 were used for model validation.To simulate the effects of watersaving measures on water use efficiency, fifteen scenarios were built by considering three measures: water-saving irrigation regimes, canal lining, and surface water supply ratio (surface water irrigation/total irrigation) reductions.Specifically, two situations were considered as canal water delivery coefficients (real water delivered by canal/water intake in canal head), namely, present and branch canal lining.Three situations of the surface water supply ratio were used: the present condition with a surface water supply ratio of 0, 0.3 and 0.4.The scenario descriptions are presented in Table 2.
The irrigation regimes were classified into the present condition, with low-intensity water saving and high-intensity water saving.The control rules for ponding water depth and soil water content in the root zone under various irrigation regimes are shown in Table 3.

Arrangements of Modeling Scenarios
Parameters including field water capacity, saturated water content, soil water specific yield, and permeability coefficient in our model were firstly calibrated by observed groundwater tables and river discharges from 14 wells and 3 rivers during the growth of paddy rice in 2015.The groundwater tables of 4 wells from 2010 to 2014 were used for model validation.To simulate the effects of water-saving measures on water use efficiency, fifteen scenarios were built by considering three measures: water-saving irrigation regimes, canal lining, and surface water supply ratio (surface water irrigation/total irrigation) reductions.Specifically, two situations were considered as canal water delivery coefficients (real water delivered by canal/water intake in canal head), namely, present and branch canal lining.Three situations of the surface water supply ratio were used: the present condition with a surface water supply ratio of 0, 0.3 and 0.4.The scenario descriptions are presented in Table 2.
The irrigation regimes were classified into the present condition, with low-intensity water saving and high-intensity water saving.The control rules for ponding water depth and soil water content in the root zone under various irrigation regimes are shown in Table 3.In the scenario simulations, we applied the measured groundwater table, river and ditch water levels on 1 January 2015, and the field capacity as the initial conditions.The daily historical meteorological data from 1956 to 2014 were regarded as the input data driving the model to calculate the water budget and water use efficiency at different scales over 59 years in the study area.

Water Balance Model Calibration and Validation
A total of 1772 observations of 14 measured groundwater tables and 492 observations of 3 sets of river discharge data during the rice growing period in 2015 were used for model calibration (Figure 5).The R 2 values of the groundwater table and river discharge between the simulated and observed data were 0.92 (Figure 5a) and 0.76 (Figure 5b), respectively.All pairs of observed and simulated groundwater and river discharge data approximated 1:1 lines (red lines in Figure 5), and the linear-fitting curves (green lines in Figure 5) were also both near 1:1 lines.However, we should point out that the simulation accuracy for river discharge was a little lower than groundwater table (Figure 4).For one thing, there might be some measurement errors for river discharge when the flow velocity is very small.For the other, field drainage after irrigation is one of the main source of river discharge and each study region contains many different fields, there were some bias between the actual irrigation time and the irrigation time inputted in the model, which could also reduce the simulation accuracy for river discharge.In addition, limited weather stations might also reduce the accuracy for our model due to the climatic variability of our experimental area.Therefore, further studies are still needed to improve our model, especially for small velocity conditions.
In the validation process, the time series measurements of 4 groundwater wells from 2010 to 2014 indicated a similar trend between observations and simulations (Figure 6).The R 2 between the observed and simulated groundwater tables of all 4 groundwater wells (QJ-9, CY-4, HW-1, and HH-2) was 0.8.Based on the model performance of calibration and validation processes, our established model captured the main characteristics of water balance in our study area.However, we also noted shortcomings in the present study.For example, a slight overestimation of groundwater tables was detected from July to September in QJ-9 and HW-1.This phenomenon might have been caused by two factors.First, the wells for monitoring groundwater tables and pumping for irrigation were usually not located in the same field.Second, due to the huge area of the study site, our model is based on PUs, which is similar with hydrologic response unit (HRU) in SWAT model and the control area is approximately 200 ha; however, the observations for groundwater tables were based on a point scale.Similar issues are found in almost all semi-distributed hydrological models, such as SWAT [53] and SAHYSMOD [54].However, our model improved the method for identifying HRU (PUs in our study) to bring simulated groundwater tables closer to the observations.Previous hydrological models (e.g., SWAT) usually define HRU based on DEM (Digital Elevation Map) data, while our model focuses on real fields and ditches in an irrigation district.Therefore, PUs in our   Based on the model performance of calibration and validation processes, our established model captured the main characteristics of water balance in our study area.However, we also noted shortcomings in the present study.For example, a slight overestimation of groundwater tables was detected from July to September in QJ-9 and HW-1.This phenomenon might have been caused by two factors.First, the wells for monitoring groundwater tables and pumping for irrigation were usually not located in the same field.Second, due to the huge area of the study site, our model is based on PUs, which is similar with hydrologic response unit (HRU) in SWAT model and the control area is approximately 200 ha; however, the observations for groundwater tables were based on a point scale.Similar issues are found in almost all semi-distributed hydrological models, such as SWAT [53] and SAHYSMOD [54].However, our model improved the method for identifying HRU (PUs in our study) to bring simulated groundwater tables closer to the observations.Previous hydrological models (e.g., SWAT) usually define HRU based on DEM (Digital Elevation Map) data, while our model focuses on real fields and ditches in an irrigation district.Therefore, PUs in our Based on the model performance of calibration and validation processes, our established model captured the main characteristics of water balance in our study area.However, we also noted shortcomings in the present study.For example, a slight overestimation of groundwater tables was detected from July to September in QJ-9 and HW-1.This phenomenon might have been caused by two factors.First, the wells for monitoring groundwater tables and pumping for irrigation were usually not located in the same field.Second, due to the huge area of the study site, our model is based on PUs, which is similar with hydrologic response unit (HRU) in SWAT model and the control area is approximately 200 ha; however, the observations for groundwater tables were based on a point scale.Similar issues are found in almost all semi-distributed hydrological models, such as SWAT [53] and SAHYSMOD [54].However, our model improved the method for identifying HRU (PUs in our study) to bring simulated groundwater tables closer to the observations.Previous hydrological models (e.g., SWAT) usually define HRU based on DEM (Digital Elevation Map) data, while our model focuses on real fields and ditches in an irrigation district.Therefore, PUs in our model were determined based on the minimum control area of each ditch, which made our model more practical and accurate.Besides, water cycle in regional scale is usually very complex and affected by many factors.Lumped model could not consider the spatial variability of these factors while distributed model need too many spatial inputs which are difficult to collect in practice.Our semi-distributed model could balance the advantages and disadvantages between lumped model and distributed model and be more suitable for water cycle simulations in regional scale [55][56][57].
Due to the lack of flow observation points in our validation period (Years: 2010-2014), river discharge data were not available for model validation.Additionally, the dynamics of soil moisture were not measured in our study, although our model could consider this process.However, previous research has indicated that it would be better, but not necessary, to calibrate and validate models for all outputs [58,59].Therefore, our model is reliable and can be used for further studies on irrigation efficiency under different scenarios at different scales.In the future study, it is necessary to explore the method for obtaining the hydraulic connection among PUs automatically based on the spatial distributions of canal, ditch, field, terrain and so on in the irrigation district.Furthermore, more data about soil, water, and crop in different irrigation district are also needed to evaluate and improve our model.

Current Irrigation Efficiency Across Different Scales
The average water balance components from 2010 to 2015 used to calculate water use efficiency (PF nws ) at different scales are shown in Table 4.Under the current conditions, groundwater was the dominant irrigation water source, accounting for more than 99% of the total irrigation.Less than 1% of the irrigation water was supplied by drainage water from ditches.Drainage water and groundwater for irrigation, which were supplied by ditches and saturated aquifers inside the scale boundary, respectively, were considered reused water.Therefore, water from saturated aquifers and ditches located outside the boundary of S1, including irrigation water from groundwater and ditches, was not regarded as reused; therefore, the reused mass was 0. At S2, although 1.622 mm/y irrigation water was derived from drainage, the amount of reused water was only 0.10 mm/y because a large amount of the drainage for irrigation was supplied by ditches outside the S2 boundary.From S4 to S7, all of the drainage water for irrigation was supplied by ditches within the scale boundary; thus, the reused water from ditches was equal to the drainage water for irrigation.From S2 to S7, because groundwater recharge was less than the irrigation pumped from groundwater, the former was considered reused water.
As the scale increased, the RF increased gradually (Figure 7).When the scale moved from S1 to S2, the reuse mass changed from 0 to 190.10 mm/y, which causing the water use efficiency to increase dramatically.As the scale further increased, the amount of reused water gradually increased but at a decelerating growth rate.Meanwhile, the change in ERFF decreased as the scale increased, which was opposite from the change in RF, indicating that the outflow across scales boundaries decreased.
decelerating growth rate.Meanwhile, the change in ERFF decreased as the scale increased, which was opposite from the change in RF, indicating that the outflow across scales boundaries decreased.The current PFnws appeared to exhibit an increasing tendency as the scale increased (Figure 8).The reason was that the outflow regarded as 'water loss' in a small scale might be reused at a larger scale, reducing the loss and improving PFnws as scale increased.Meanwhile, the scale effect of PFnws could be represented by a power function, which is similar to that reported Xie and Cui [60] but is different from the findings reported by Hafeez et al. [13], who used a linear function to indicate the relationship (Figure 8).Furthermore, water use efficiency sharply increased (26.9%) when the scale increased from the field scale (S1) to S2 (Figure 8).The reason might be that S1 only contained the surface water layer and the crop root layer, whereas S2 included all of the saturated and unsaturated zones [13].
Furthermore, the dominant reuse of the return water in the study area involved pumped groundwater recharge for irrigation.
The amount of water reuse for calculation of PFnws varied from different scales (Equation ( 1)).In S1, deep percolation from root zone was regarded as return water and no water reuse in this scale (Table 1).Nevertheless, the return water in S1 can be reused for irrigation in S2 by pumping well.Therefore, water reuse should be considered in S2.Meanwhile, the amount of water reuse from The current PF nws appeared to exhibit an increasing tendency as the scale increased (Figure 8).The reason was that the outflow regarded as 'water loss' in a small scale might be reused at a larger scale, reducing the loss and improving PF nws as scale increased.Meanwhile, the scale effect of PF nws could be represented by a power function, which is similar to that reported Xie and Cui [60] but is different from the findings reported by Hafeez et al. [13], who used a linear function to indicate the relationship (Figure 8).The current PFnws appeared to exhibit an increasing tendency as the scale increased (Figure 8).The reason was that the outflow regarded as 'water loss' in a small scale might be reused at a larger scale, reducing the loss and improving PFnws as scale increased.Meanwhile, the scale effect of PFnws could be represented by a power function, which is similar to that reported Xie and Cui [60] but is different from the findings reported by Hafeez et al. [13], who used a linear function to indicate the relationship (Figure 8).Furthermore, water use efficiency sharply increased (26.9%) when the scale increased from the field scale (S1) to S2 (Figure 8).The reason might be that S1 only contained the surface water layer and the crop root layer, whereas S2 included all of the saturated and unsaturated zones [13].
Furthermore, the dominant reuse of the return water in the study area involved pumped groundwater recharge for irrigation.
The amount of water reuse for calculation of PFnws varied from different scales (Equation ( 1)).In S1, deep percolation from root zone was regarded as return water and no water reuse in this scale (Table 1).Nevertheless, the return water in S1 can be reused for irrigation in S2 by pumping well.Therefore, water reuse should be considered in S2.Meanwhile, the amount of water reuse from Furthermore, water use efficiency sharply increased (26.9%) when the scale increased from the field scale (S1) to S2 (Figure 8).The reason might be that S1 only contained the surface water layer and the crop root layer, whereas S2 included all of the saturated and unsaturated zones [13].
Furthermore, the dominant reuse of the return water in the study area involved pumped groundwater recharge for irrigation.
The amount of water reuse for calculation of PFnws varied from different scales (Equation ( 1)).In S1, deep percolation from root zone was regarded as return water and no water reuse in this scale (Table 1).Nevertheless, the return water in S1 can be reused for irrigation in S2 by pumping well.Therefore, water reuse should be considered in S2.Meanwhile, the amount of water reuse from pumping well was very big in S2 (Table 4), which made water use efficiency exhibit the largest increase at this scale.As the scale continued to increase, the water use efficiency did not improve significantly, mainly because the reuse of surface discharge was relatively low in the study area.

Scale Effects of Simulated Water-Saving Measures on Irrigation Efficiency
Based on our simulations of different scenarios (Table 2), the influence of different irrigation regimes on water use efficiency (PF nws ) at different scales was as indicated in Figure 9. Water saving irrigation regimes (WSIRs) had the most significant positive effect (an increase by 21.2%) on S1 (field scale) among scales in this study.The water-saving effect (WSE) of WSIR on other scales was influenced by the surface water supply ratio (SWSR).First, without surface water supply (Figure 9a), WSIR provided almost no WSE at scales other than S1.For example, the WSE at S7 (irrigation area scale) was only 1.13% even if high-intensity WSIR was adopted.Second, the WSE at scales from S2 to S7 began to increase with increasing SWSR.When the SWSR reached 0.4 (Figure 9b), the WSE at S7 with the application of a high-intensity WSIR was 7.17%.Third, the difference in WSE at all scales of the low-intensity and high-intensity WSIR was not significant when SWSR was less than 0.4.
Sustainability 2018, 10, 0047 10.3390/su10010047 13 of 18 pumping well was very big in S2 (Table 4), which made water use efficiency exhibit the largest increase at this scale.As the scale continued to increase, the water use efficiency did not improve significantly, mainly because the reuse of surface discharge was relatively low in the study area.

Scale Effects of Simulated Water-Saving Measures on Irrigation Efficiency
Based on our simulations of different scenarios (Table 2), the influence of different irrigation regimes on water use efficiency (PFnws) at different scales was as indicated in Figure 9. Water saving irrigation regimes (WSIRs) had the most significant positive effect (an increase by 21.2%) on S1 (field scale) among scales in this study.The water-saving effect (WSE) of WSIR on other scales was influenced by the surface water supply ratio (SWSR).First, without surface water supply (Figure 9a), WSIR provided almost no WSE at scales other than S1.For example, the WSE at S7 (irrigation area scale) was only 1.13% even if high-intensity WSIR was adopted.Second, the WSE at scales from S2 to S7 began to increase with increasing SWSR.When the SWSR reached 0.4 (Figure 9b), the WSE at S7 with the application of a high-intensity WSIR was 7.17%.Third, the difference in WSE at all scales of the lowintensity and high-intensity WSIR was not significant when SWSR was less than 0.4.WSIR exhibited a relatively substantial WSE at the field scale (S1).The primary reason might be that the amount of reused water was not considered at this scale; thus, water use efficiency (PFnws) was influenced by the gross water supply (reused water was not considered) [61].When WSIR was adopted, the gross water supply was reduced, causing an improvement in water use efficiency.Furthermore, at scales other than S1, the reused water was considered.If we ignore the change in ET after applying WSIR at scales from S2 to S7, the water use efficiency was only influenced by the net water supply, which was the gross water supply minus the reused water amount.After WSIR was implemented, decreases in gross water supply and reused water had a positive and negative effect on the improvement of water use efficiency, respectively, and WSE obviously decreased when the effects of the two factors were combined.
When the branch canal water delivery coefficient increased from 0.65 to 0.80 through canal lining, the WSE at scales from S3 (4771.52 ha) to S5 (4.29 × 10 4 ha) were more obvious than those at other scales (Figure 10a).However, the WSE was less than 2% because of the relatively low SWSR.When the SWSR was 0.4 (Figure 10b), the WSE at scales from S4 (3.24 × 10 4 ha) to S7 (10.08 × 10 4 ha) increased slightly, reaching approximately 4.3%.However, the WSE of the canal lining decreased if WSIR was applied (Figure 10c).In this study, when the SWSR was less than 0.4, the WSE of the canal lining at all scales was almost negligible when high-intensity WSIR was adopted.
Based on the above analysis, our simulation results for different scenarios indicated that the influence of water savings and water resources regulation measures on reused water across scales exhibited nonlinear characteristics.Thus, WSE presented complex responses, which were influenced WSIR exhibited a relatively substantial WSE at the field scale (S1).The primary reason might be that the amount of reused water was not considered at this scale; thus, water use efficiency (PF nws ) was influenced by the gross water supply (reused water was not considered) [61].When WSIR was adopted, the gross water supply was reduced, causing an improvement in water use efficiency.Furthermore, at scales other than S1, the reused water was considered.If we ignore the change in ET after applying WSIR at scales from S2 to S7, the water use efficiency was only influenced by the net water supply, which was the gross water supply minus the reused water amount.After WSIR was implemented, decreases in gross water supply and reused water had a positive and negative effect on the improvement of water use efficiency, respectively, and WSE obviously decreased when the effects of the two factors were combined.
When the branch canal water delivery coefficient increased from 0.65 to 0.80 through canal lining, the WSE at scales from S3 (4771.52 ha) to S5 (4.29 × 10 4 ha) were more obvious than those at other scales (Figure 10a).However, the WSE was less than 2% because of the relatively low SWSR.When the SWSR was 0.4 (Figure 10b), the WSE at scales from S4 (3.24 × 10 4 ha) to S7 (10.08 × 10 4 ha) increased slightly, reaching approximately 4.3%.However, the WSE of the canal lining decreased if WSIR was applied (Figure 10c).In this study, when the SWSR was less than 0.4, the WSE of the canal lining at all scales was almost negligible when high-intensity WSIR was adopted.
Based on the above analysis, our simulation results for different scenarios indicated that the influence of water savings and water resources regulation measures on reused water across scales exhibited nonlinear characteristics.Thus, WSE presented complex responses, which were influenced both by recycling of the return flow as well as the responses of water balance components at different scales to water-saving measures [62,63].

Conclusions
Our study evaluated the scale effect of water saving in the BLH Basin, a rice cultivation area located on the Sanjiang Plain, Northeast China.Under the current situation, percolation and drainage were regarded as 'water losses' at a small scale but were reused by pumping for irrigation.When the scale increased, the RF increased and the ERFF decreased, and water use efficiency exhibited an increased tendency with increased scale, which could be expressed by a power function by the scale

Conclusions
Our study evaluated the scale effect of water saving in the BLH Basin, a rice cultivation area located on the Sanjiang Plain, Northeast China.Under the current situation, percolation and drainage were regarded as 'water losses' at a small scale but were reused by pumping for irrigation.When the scale increased, the RF increased and the ERFF decreased, and water use efficiency exhibited an increased tendency with increased scale, which could be expressed by a power function by the scale area.Furthermore, due to the theoretical return water being reused through pumped groundwater, the water use efficiency improved dramatically when the scale boundary expanded from a crop root layer to the whole unsaturated-saturated soil zone.
The WSE of different water-saving and water resource regulation measures proved to be scale dependent.When the SWSR was zero, the application of WSIR achieved the greatest WSE (21.2%) at S1 of all 7 scales, whereas it resulted in almost no WSE at other scales except S1.However, the WSE of WSIR at other scales, except S1, gradually increased when the SWSR increased.When the branch canal water delivery coefficient increased from 0.65 to 0.80 with canal lining and the SWSR was 0.3, without WSIR, the WSE at scales from S3 to S5 correspondingly improved (<2%).Moreover, the increase in SWSR was enhanced, whereas the WSIR implementation might reduce the WSE of branch canal lining at each scale.Furthermore, the PF nws at scales larger than S3 began to ascend when the SWSR decreased from 0.4 to 0.3 without WSIR.However, a further reduction in SWSR from 0.3 to 0.1 caused no changes in PF nws at S1 and S7.
Overall, our study established a semi-distributed water-balance model to consider the effects of crop and human activities on hydrological dynamics at different scales in an irrigation district.Our model was scientifically evaluated for both river discharge (R 2 > 0.75) and groundwater tables (R 2 > 0.8).And can be used for evaluating water cycles and water use efficiency at rice cultivation sites.However, due to the complicated mechanisms of water cycle dynamics and the huge spatial data requirements for a real irrigation district, the present model should be improved in certain respects, such as refining the simulation unit, identifying flow directions among units, and representing surface water dynamics processes.

Figure 1 .
Figure 1.Location of the study area.(a) The location of the Bielahonghe (BLH) Basin; (b) Details of the experimental sites.QLS (Qinglongshan Farm), NJ (Nongjiang Farm), HH (Honghe Farm), CY (Chuangye Farm), HW (Hongwei Farm), HC (observation site near BLH River), and QJ (Qianjin Farm) are the names of the field sites, and the numeral suffixes indicate different groundwater observation sites.

Figure 1 .
Figure 1.Location of the study area.(a) The location of the Bielahonghe (BLH) Basin; (b) Details of the experimental sites.QLS (Qinglongshan Farm), NJ (Nongjiang Farm), HH (Honghe Farm), CY (Chuangye Farm), HW (Hongwei Farm), HC (observation site near BLH River), and QJ (Qianjin Farm) are the names of the field sites, and the numeral suffixes indicate different groundwater observation sites.

Figure 2 .
Figure 2. Diagram of the water balance model (E is evaporation, and ET is evapotranspiration).

Figure 2 .
Figure 2. Diagram of the water balance model (E is evaporation, and ET is evapotranspiration).

Figure 3 .
Figure 3. Sketch diagram of the water balance in the polygon unit (PU), ditch units (DU) and canal unit (CU).

Figure 3 .
Figure 3. Sketch diagram of the water balance in the polygon unit (PU), ditch units (DU) and canal unit (CU).

Figure 4 .
Figure 4. Description of different scales (green color).Abbreviations from S1 to S7 indicate the 7 different scales in our study.

Figure 4 .
Figure 4. Description of different scales (green color).Abbreviations from S1 to S7 indicate the 7 different scales in our study.

Figure 5 .
Figure 5. Observed and simulated groundwater table (a) and river discharge (b) in the calibration period (Year: 2015).

Figure 5 .
Figure 5. Observed and simulated groundwater table (a) and river discharge (b) in the calibration period (Year: 2015).

Figure 5 .
Figure 5. Observed and simulated groundwater table (a) and river discharge (b) in the calibration period (Year: 2015).

Figure 7 .
Figure 7. Relationships among the water reuse fraction (RF), effective return water fraction (ERFF) and study area in the present context.

Figure 8 .
Figure 8. Relationship between water use efficiency (PFnws) and the study area in the current situation.

Figure 7 .
Figure 7. Relationships among the water reuse fraction (RF), effective return water fraction (ERFF) and study area in the present context.

Sustainability 2018, 10 ,
0047 10.3390/su10010047 12 of 18 decelerating growth rate.Meanwhile, the change in ERFF decreased as the scale increased, which was opposite from the change in RF, indicating that the outflow across scales boundaries decreased.

Figure 7 .
Figure 7. Relationships among the water reuse fraction (RF), effective return water fraction (ERFF) and study area in the present context.

Figure 8 .
Figure 8. Relationship between water use efficiency (PFnws) and the study area in the current situation.

Figure 8 .
Figure 8. Relationship between water use efficiency (PF nws ) and the study area in the current situation.

Figure 9 .
Figure 9.The ratios of rice evapotranspiration to net water supply (water use efficiency, PFnws) across scales under present (blue lines), low (black lines), and high-intensity water-saving (red lines) irrigation regimes.

Figure 9 .
Figure 9.The ratios of rice evapotranspiration to net water supply (water use efficiency, PF nws ) across scales under present (blue lines), low (black lines), and high-intensity water-saving (red lines) irrigation regimes.

Sustainability 2018, 10 ,
0047 10.3390/su10010047 14 of 18 both by recycling of the return flow as well as the responses of water balance components at different scales to water-saving measures[62,63].

Figure 10 .
Figure 10.The ratios of rice evapotranspiration to net water supply (PFnws) across scales with (blue lines) and without (black lines) canal lining.

Figure 10 .
Figure 10.The ratios of rice evapotranspiration to net water supply (PF nws ) across scales with (blue lines) and without (black lines) canal lining.

Table 1 .
Boundaries and water reuse categories at different scales.
Note: Present a indicates that the water delivery coefficient of branch canals and main canals was 0.65 and 0.70, respectively; Present b indicates that the surface water supply ratio was 0.

Table 3 .
Rules of ponding depth and soil water content for different irrigation regimes (θ s indicates soil water content at saturation).

Table 4 .
Simulated average water balance components from 2010 to 2015 and water use efficiency across scales.All the values in this table are the average of the study area at each scale rather than in the area planted in rice.