Check dams are the most widely applied engineering structures for soil and water conservation in the erodible regions. In the gully-dominated regions, check dams have been widely constructed for retaining floodwater, intercepting soil sediment, improving gully slope stabilities, and increasing farmlands. At the beginning of construction, a water body similar to a small reservoir usually existed behind the check dam. Gradually, fertile cropland formed for grain production in the region, where a water body was filled by sediment deposition.1
Unlike reservoirs, check dams are not common around the world. It has only been reported that check dams were implemented for soil erosion control in some countries and regions, such as France, Italy, Spain, China, etc.2,3 The majority of research on check dams has been conducted to investigate their impacts on geomorphology, channel morphology, river bed stability, and sediment yield through both field measurement and numerical simulation. It has been concluded that the check dams captured sediments and led to a decrease in the longitudinal gradient upstream, but accelerated downstream erosion because of the increased flow transport capacity and scour processes.1,4,5 Bombino et al.6 addressed statistically significant differences in channel form, sediment, and vegetation among the upper, the lower reaches, and the intermediate sections of check dams in Southern Italian catchment. A reducing rate of 77% sediment yield by check dams has been reported by Boix-Fayos et al.2 in the Rogativa catchment of SE Spain. Ran et al.7 analyzed the coarse sediment retention by check dams in five typical catchments in the midstream of the Yellow River, and found that the average sediment reduction ratio could reach 60% when the check dams area accounted for more than 3% of the drainage area.
In Chinese Loess Plateau, check dams have attractive advantages because of unique environmental settings and regional food supply needs.5,8,9 For years, check dams have been constructed for soil erosion control and agricultural production improvement. The check dams have become one of the most effective strategies for soil erosion control in the Loess Plateau. Especially during the recent six decades, thousands of check dams have been constructed in the region. By 2005, check dams have been built on the Loess Plateau over the past 50 years and approximately of sediments have been captured by these dams.10 Although considerable studies have been conducted as mentioned above, the locations, storage capacities, and outflow methods have not been inventoried and documented. Thus, further understanding the mechanisms on the hydrological effects of check dams was impossible.
Traditionally, check dams and their storage capacity are estimated by using the topographic maps according to the in situ hydrographic survey with an electronic planimeter.11 However, this approach is usually laborious, costly, and time consuming. Advances in remote sensing, geographic information systems, and computer technology have made satellite images convenient and less expensive for reservoir or dam measurement and monitoring.12,13 Remote-sensing images provide invaluable and timely information on spatial distribution, spectral and temporal attributes for reservoirs, and dam operation and management.14,15 However, research on check dam identification and quantitative estimation of storage capacities by means of remote-sensing images are very rare.16,17 The few studies available mainly focused on individual reservoirs or a rather small watershed and as a result, only some tentative applications have been reported.18
To better understand the potential effects of check dams on hydrological processes and sediment yield at catchment scale, an accurate estimation of the check dam distribution and their storage volume becomes very important and necessary. The objectives of this study, therefore, were:
1. to identify the check dams in a Loess Plateau catchment by making use of multisource data;
2. to quantify the actual storage capacity of the check dams through established empirical relationships; and
3. to analyze the impacts of check dams on streamflow and sediment load in combination with recently observed dataset.
Study Area and Data
The study area, Huangfuchuan catchment, is one of the most severe soil erosion areas with the average annual sediment modulus higher than (1955 to 2000). The Huangfuchuan River is a first-order tributary of the midstream Yellow River, located in the wind-water erosion criss-cross region. The river originates from the eastern part of the Ordos Plateau, flows through the transitional zone of the Loess Plateau and Inner Mongolia grassland, and empties into the main stream of the Yellow River in Fugu County of Shaanxi Province (Fig. 1). The Huangfuchuan catchment covers an area of consisting of two subtributaries: Nalin River and Changchuan River. The main stream of the Huangfuchuan River is 137 km long with an average channel slope of 2.7%. The watershed is located in the transitional belt of warm temperate and mesothermal zones with the average annual precipitation of about 380 mm, influenced by southeast monsoon in summer, and characterized by semi-arid continental climate.19 Frequently occurring floods in July and August often cause tremendous soil erosion in the catchment. Average annual discharge and sediment load are and , respectively, between 1955 and 2010, nearly 80% of which is concentrated in the rainy season from June to September.
The catchment can be described as “hilly and gully” landscape, and the land surface material is dominated by weathering sandstone, silty loess, and desert sand with coarser particles. According to the soil erosion type, the basin can be divided into three subregions: (1) loess hilly–gully region with an area of , mainly located in the eastern and southwestern parts of the watershed; (2) the sandy loess hilly–gully region with an area of , located between Nalin and Changchuan tributaries; and (3) the sandstone weathering hilly–gully region with an area of , located in the northwestern part of the catchment. Due to the extremely high gully density and soil erosion rates, check dams have become the dominant soil and water conservation measures in the Huangfuchuan catchment. By 2010, approximately 70% of the area in the Huangfuchuan watershed was controlled by check dams.
The data used in this study were collected from topographic maps (1978), Landsat images (1990, 2007), and images acquired in 2010 from Google Earth. Table 1 lists all the available datasets and sources. In total, we collected 28 topographic maps with the resolution of 1:50,000 from Bureau of Survey and Mapping Shaanxi Province. The images were downloaded from the Landsat imagery database.20 The coverage of an individual Landsat TM image is . The resolution for bands 1 to 5 and 7 is 30 m, while it is 120 m for band 6. Cloud free ( cloud cover) images were selected. All the images were taken during the flood season, extending from August to September, to avoid various undesirable errors caused by adopting images taken in different seasons. Ran and Lu21 confirmed that the water surface area does not change significantly during the flood seasons, thus, the downloaded satellite images could capture the largest water surface area of the check dams.
Summary of the processed maps and Landsat images.
|Topography map||1∶50,000||1978||Bureau of Survey and Mapping Shaanxi Province|
|Landsat 5||30 m||127 to 33||11/09/1989||http://earthexplorer.usgs.gov|
|Landsat 5||30 m||127 to 32||13/08/1990|
|Landsat 5||30 m||127 to 32||10/09/2006|
|Landsat 5||30 m||127 to 33||25/09/2007|
|Google Earth||<15 m||23/06/2010||Google Earth images|
Check Dam Identification from Multisource Data
The topographic maps were geocoded to universal transverse mercator (UTM), zone 49N, using the WGS84 ellipsoid, and were used to delineate the check dams in ArcMap manually. Landsat images were subjected to atmospheric correction, radiometric calibration, and radiometric rectification procedures to facilitate comparability between dates.22 The Landsat images were geo-referenced to 1:50,000 scale topographic maps using second-order polynomial transformation with the root mean square error of below 0.5 pixels (), and registration was verified visually by overlaying and swiping registered images. Due to the strong absorption and reflectance region in the infrared, visible red, and near-infrared bands (bands 7, 4, and 3) provided in the Landsat images, the combination of these bands was used to distinguish water bodies and land. The check dams or reservoirs were extracted by using the maximum likelihood supervised classification method in Erdas Imagine software.14,15 During the interpretation process, each 30-m pixel is classified as check dam (open water surface) or nondam area. Ponds and drainage network were distinguished from check dams through field survey. A database of the check dams (including volumes, built time, location, dyke height, remained/unfilled storage) were collected from Hydrology Bureau of the Yellow River Water Resources Commission. The volumes, built time, and dyke height were obtained through local census, and the remained storage for each check dams was calculated by combining bathymetric survey and storage curve analysis.
Currently, Google Earth has provided images of resolution of covering the entire river basin.21 In addition, GPS field survey has similar accuracy of 10 to 20 m for ground truthing.11,15 Thus, we derived the spatial distribution and area extent of the check dams to verify the interpretation results of Landsat images. The images displayed by Google Earth were taken on June 23, 2010, coinciding with the capture time of the Landsat images. An update dataset for check dams in the Huangfuchuan catchment was obtained through images interpretation from Google Earth Pro.
Field survey was carried out in the beginning and end of the rainy season (May 2012 and October 2012). The outlines of the check dams were obtained with a Garmin 76 GPS instrument. The data stored on the GPS was converted into text formats which can be read by ArcCatalog. The points were connected together to form polygons. The surface areas and perimeters were then calculated.
To assess accuracy of the extracted check dam area from different data sources, both normalized difference area index (NDAI) and deviation area index (DAI) were applied.14,23 NDAI and DAI are defined as follows:
The NDAI value falls within the range between and 1 with values close to 0 giving the best fit between the surface areas obtained from GPS and other sources, respectively. Values increasing to both extremes indicate increasing deviation between and . Similarly, the DAI values also lie between and 1, with values close to 0 meaning the best linear fit, whereas values increasing to both extremes indicate larger deviation between and .
Storage Capacity Estimation
If reservoirs or check dams can be assumed to approximate the shape of a square pyramid diagonally cut in half,15,21,24 an empirical relationship between the storage capacity and its area can be expressed as15 However, it should be noted that the parameters may differ greatly in regions with geomorphologic heterogeneities.14,21,23 Thus, more attention should be paid when using the empirical relationship.
Correlation Between Field and Satellite Identified Check Dams
Figure 2 presents the identified water surface area of the check dams from different types of dataset. In general, a good linear relationship was detected between the results from Landsat images and Google Earth images. As shown in Fig. 2(a), most points distributed close to fit line with the relatively high correlation coefficient of 0.98. Similarly, the water surface area from satellite images fitted well with field measurements [Fig. 2(b)]. However, it should be noted that some check dams built in the late 1950s and 1960s had been destroyed by flash floods or filled by deposited sediment. For example, the check dam in Fig. 3(a) was full of water in 1978 and dried up in 1990. A lifetime of about 10 years for small check dams to trap sediment efficiently was reported by Xu et al.5 Thus, the damaged or filled check dams may lead to some errors in the extraction of check dams from the images.
To further verify the water surface area derived from Landsat images, field surveys in 11 and 14 check dams were undertaken in May and October 2012, respectively. The NDAI and DAI were applied to detect the interpretation results. Table 2 compared the water surface area in check dams derived from images and field measurements. There is a very good linear fit between these two datasets. In most cases, we found that the water surface area from the Landsat images has overestimated the measurement values. Probably, the overestimation may result from the lower resolution of the Landsat images by which the surrounding nonwater features would have been misinterpreted as water. Furthermore, the results of supervised classification are largely dependent on the accuracy of the training data and their representative. In addition, the differences in water surface area may generate great discrepancies because of the uneven distributed precipitation in different seasons. For example, a minority of the check dams were dried up in May 2012, but there was a large area of water in October due to heavy storms in July.
Comparison of water body area in check dams derived from images and field measurements.
|ID||Longitude||Latitude||AreaGPS (ha)||AreaSat (ha)||DAI (%)||NDAI (%)||Storage (104 m3)||Built date|
Storage Capacity Estimation in the Catchment
In total, datasets of 74 check dams with spatial distribution, built date, storage capacities, and dike height were provided by the Yellow River Basin Committee. By representing various check dams with different storage capacities, these check dams were employed to develop the relationship between the water surface area and the storage capacity. As shown in Fig. 4, we investigated the relationship between the storage capacity and the water surface area as well as the dam-controlled area. The results indicate relatively good power relationship between the water surface area and the storage capacity, and logarithmic relationship between check dam controlled area and storage capacity. In contrast, Fig. 4(a) confirmed that the power equation represents a better relationship between check dam water surface area and storage capacity [Fig. 4(b)]. Thus, an empirical equation for calculating the check dams’ storage is
Equation (4) was applied to estimate the storage capacity of each check dam, and the total storage capacity was estimated by summing all individual check dams’ storage together (Table 3). The total storage capacities in the Huangfuchuan catchment have increased from (1978) to (2010), and the average storage capacity for the check dams is nearly . The estimated errors for storage capacities mainly result from the spatial heterogeneities in geometric shape and extracted water surface area. Ran and Lu21 addressed that each different geometric shape requires a unique parameterization to precisely calculate its storage capacity. In addition, the over- or underestimated water surface area may lead to misestimation of storage capacity for the check dams. However, compared with traditional field survey approach, the proposed method in this study was useful, efficient, and less expensive, and the estimation accuracy is satisfactory.
Statistics of check dams in the Huangfuchuan watershed from 1978 to 2010.
|Number of check dams||390||445||504||567|
|Water surface area (km2)||10.76||14.06||16.12||17.78|
|Dams controlled area (km2)||1226.56||1531.61||1925.84||2216.47|
|Total storage capacity (108 m3)||2.65||3.44||4.06||4.11|
Historical Dynamics of Check Dams
Figure 5 shows historical dynamics of check dams derived from topographic maps (1978), Landsat images (1990 and 2007), and images of Google Earth in 2010. The results indicate that a general increase for the number of check dams can be clearly seen in the Huangfuchuan catchment. As well, the spatial distribution indicates that most of the check dams lie in the eastern part of the watershed. Taking the map of 2007, e.g., check dam density in the Wulan stream, Hujitu stream, and Tela stream show evident increasing trends. In 1978, there were only 16 check dams in the Wulan stream, while it had reached 118 in 2007 (Fig. 5).
Table 3 summarizes statistics of the check dams derived from multisource data in the Huangfuchuan catchment. By 1978, there were 390 check dams built in the Huangfuchuan catchment, with the water surface area of . The dam-controlled area was , accounting for about one third of the whole catchment. By 2010, the number of check dams had reached 567, and the dam-controlled area was . This is consistent with the field survey data reported by the Hydrology Bureau of the Yellow River Water Resources Commission. Minor discrepancies were detected between Jiao et al.25 and our estimation. They addressed that 392 major check dams had been constructed until 1992 in the Huangfuchuan catchment, and approximately sediment had been trapped behind the check dams.
Potential Impacts of Check Dams on Streamflow and Sediment Load
Figure 6 shows temporal trends of discharge and sediment load at Huangfu station. The abrupt changes were detected in 1979 for both streamflow and sediment load by the sequential Mann–Kendall test.26,27 Simple linear trend analysis indicates that the time series presents a significant downward trend within the confidence level of 95%. Average decreasing rates of and ton are estimated in streamflow and sediment load, respectively. Average annual streamflow at Huangfu station was (1980 to 2010), accounting for only 44.8% of that during 1955 to 1979. Mean annual sediment load decreased from ton (1955 to 1979) to ton (1980 to 2010). Considerable evidence exists that check dams’ construction was responsible for the significant reduction in both streamflow and sediment load.8,10 Ran28 addressed that streamflow reduced by soil and water conservation measures ranged from 0.14 to from 1970 to 1996. Average annual reducing rates for sediment load lay in the range of 5.4 to 23.6 Mt. Among these, check dams are the dominant factors retaining discharge and sediment in the watershed.29 Furthermore, sediment trapped by the check dams accounted for of the total reduced sediment load in the Huangfuchuan catchment.
To better understand the impacts of check dams on streamflow and sediment load, the hydrological time series at Huangfu station were divided into two periods (1955 to 1979 and 1980 to 2010) by the Mann–Kendall test method. A referenced period from 1955 to 1979 was assumed to be the period without intensive human influences on the hydrological process and the changing period was between 1980 and 2010, which was highly affected by human activities. Figures 6(c) and 6(f) show the daily flow duration curves within two different periods. The magnitude of the daily flow generally decreased during the changing period compared with the referenced period. It can also be clearly seen that zero discharge accounts for much more during the changing period. The daily sediment delivery ratio with 0 values is nearly 60% within 1955 to 1979 but reduced to 32.6% during the changing period. There is no discharge measured at Huangfu station during all of 2011. Even a storm in 14th of August with 58.4-mm rainfall did not produce surface flow.30 Thus, it can be concluded that numerous soil and water conservation measures, especially, the check dams trapped considerable discharge and sediment in the Huangfuchuan catchment.
Sediment Particle Sizes Associated with Check Dam Construction
In general, hyper-concentrated sediment flow can be trapped in the check dams. In most cases, the check dams consisted of the embankment, spillway, and outlet. Before 1980, most check dams had no spillway or outlet, and all the upstream incoming sediment would be trapped by the dams. In the check dams with spillway, the sediment flow will be filtered. In principle, the coarse grain sediment will deposit rapidly in the water surface area of the dam, and the sediment with fine particle sizes will be flushed away through the spillway. In such cases, most coarse-grain sediment is trapped before the dams.
In the Huangfuchuan catchment, frequently occurring storms mainly concentrated in July and August. Sediment load during these months contributed approximately 70% of the total. In recent years, zero flow and sediment lasted several months during the dry season. Thus, sediment particle sizes during July and August may represent the effects of check dams on them. As shown in Table 4, observed data at Huangfu hydrological station indicated that mean and median sediment particle diameter illustrated evident decreasing trends, presenting the trapping effects of the check dams. The percentage of sediment particle diameter increased evidently from 1966 to 2010. Both the median and mean value of sediment particle diameter displayed particularly low values during 2005 and 2010 (bold values in Table 4). This is similar to the diameter changes in the middle reaches of the Yellow River.28,29
Sediment particle sizes during July and August at Huangfu station.
|Periods||Mean sediment load (104 ton)||<0.025 mm (%)||0.025 to 0.05 mm (%)||>0.05 mm (%)||Median (Ref. 29)||Mean (Ref. 29)|
|1966 to 1969||2928.16||34.95||15.56||49.48||0.059||0.152|
|1970 to 1979||2822.06||39.50||15.05||45.45||0.059||0.132|
|1980 to 1989||1989.64||41.61||15.17||43.22||0.051||0.127|
|1990 to 1995||1027.19||52.67||12.48||34.85||0.030||0.104|
|2005 to 2010||598.3||57.63||16.85||25.52||0.022||0.048|
Driving Factors of Check Dam Construction
In most cases, the check dams were constructed for soil erosion control.31 However, other objectives such as improving agricultural productivities and building railway or highway become dominant factors. As shown in Fig. 7(a), nearly 10 check dams were constructed for protecting the highway in a region of the Huangfuchuan catchment. Road construction has become an important forcing variable for check dam building. Statistics indicate that nearly 20% of the check dams are connected with the road or used as river bridges. On the other hand, the major aim for building the check dams is to increase agricultural production [Fig. 7(b)]. In general, food production in dam cropland can be 6 to 10 times higher than that in the slope cropland. The average yield is with some even up to .32 Accordingly, dam cropland accounts for about 9% of the total cropland area in the Loess Plateau, whereas the food production amounts to 20.5% of the total food production.8
This study identified the spatial distribution, water surface area, and controlled area of the check dams by making use of topographic maps and remote-sensing images in combination with a field survey. An empirical equation was applied to estimate the storage capacity of the check dams associated with the water surface area. Potential effects of the check dams on streamflow and sediment load were additionally addressed. The conclusion of the study can be summarized as follows:
A good correlation between the measured and derived water surface area of the check dams from maps and images was detected. The check dams’ information extracted from the remote sensing images is the same as those from Google Earth images, suggesting that the topographic maps and Landsat images can be used for identifying the check dams. Historically, the number of check dams has increased from 390 in 1978 to 567 in 2010, and the dam-controlled area has reached up to , accounting for nearly two thirds of the whole catchment. The large numbers of check dams indicate their substantial trapping effects on streamflow and sediment load in the watershed. Filtering effects of check dams have been confirmed due to evident decreasing sediment particle sizes. This study proved that the proposed approach can be applied for extracting check dam information and may provide decision supports for soil and water conservation in the Loess Plateau.
The work was supported by the National Natural Sciences Foundation of China (Grant No. 41201266), the Major Programs of the Chinese Academy of Sciences (KZZD-EW-04-03), West Light Foundation of the Chinese Academy of Science (2011ZD03), and the Fundamental Research Funds for the Central Universities (QN2013071). The authors would like to express their great appreciation to the Hydrology Bureau of the Yellow River Water Resources Commission for providing valuable climatic and hydrological data as well as to the reviewers for very valuable comments that greatly improved the quality of the paper.