Next Article in Journal
Groundwater Modelling in Karst Areas
Previous Article in Journal
Community Composition and Function of Bacteria in Activated Sludge of Municipal Wastewater Treatment Plants
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Water-Saving Soil Conservation Measures Should Be Used in Northern China: Evidence from Runoff Plot Data

1
Key Laboratory of Water Cycle and Related Land Surface Processes, Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing 100101, China
2
College of Resources and Environment, University of Chinese Academy of Sciences, Beijing 100049, China
Water 2021, 13(6), 853; https://doi.org/10.3390/w13060853
Submission received: 17 February 2021 / Revised: 12 March 2021 / Accepted: 18 March 2021 / Published: 20 March 2021
(This article belongs to the Section Hydrology)

Abstract

:
Most of the current studies on soil conservation measures mainly focus on their soil control effect, neglecting their impact on water quantity. In the present study, the latest seven years (2014–2020) of monitored data from 22 runoff plots in the upstream catchment of the Miyun Reservoir, Beijing were used to evaluate the effects of slope, rainfall, and soil conservation measures on soil and water loss, and some implications were given in this water-scarce region. Excluding the impact of soil conservation measures, soil loss increased with the slope gradient and slope length. Runoff and soil loss were greatly affected by the rainfall amount and maximum 30-min rainfall intensity on the bare and cultivated slopes, or by rainfall amount and rainfall duration on almost all of the plots with soil conservation measures. The results indicated that the bare soil suffered the most severe soil loss, with a mean annual soil loss rate (SLR) of 4325 t km−2 year−1, followed by the cultivated lands without any measure, with an annual SLR of above 3205 t km−2 year−1. Contour tillage cannot effectively control soil loss on steep slopes. The vegetation measures and terrace, level bench, and fish scale pits, as well as their combinations, can decrease runoff by above 86% and decrease soil loss by 95%, respectively. Water-saving measures should be implemented in the study region. The measures, such as vegetation coverage, terracing, contour tillage, etc., should be carefully implemented on slopes. Bare and cultivated lands should further be implemented with soil conservation measures in this and similar regions in the world.

1. Introduction

Water erosion is one of the biggest ecological and environmental problems worldwide [1]. It can induce both on-site fertile topsoil and nutrient losses and off-site sediment siltation in river systems, and cause water resource use problems [2,3,4,5]. Severe hazards induced by water erosion have been noticed a long time ago, both in developing and developed countries [6,7].
Soil conservation measures, such as contour tillage, terrace, and vegetation plantation, were widely used in the world to combat water erosion, and many studies have been published until very recently. For example, Zhao et al. [8] summarized the runoff plot data on the Chinese Loess Plateau, and pointed out that soil conservation measures had greatly reduced soil and water losses, but not down to background levels (i.e., 500 t km−2 year−1). Wolka et al. [9] reviewed the effects of soil conservation measures on runoff, soil loss, and crop yields in Sub-Saharan Africa that showed that the effect of crop slope barrier soil conservation techniques on crop yields varied with the rainfall slope. Maetens et al. [10] evaluated the effects of soil conservation measures in reducing plot runoff and soil loss in Europe and the Mediterranean. Xiong et al. [1] also gave a global review of the effects of soil conservation measures on soil loss control. These studies greatly improve our understanding of the effects of soil conservation measures on soil and water losses. However, most of the studies focus on the reduction effects of soil conservation measures, such as contour tillage, no tillage, horizontal terracing, level ditching, etc., on soil and water losses, considering less their effect on downstream water resource use.
The Miyun Reservoir is one of the most important drinking water sources in Northern China, and provides around 70% of the drinking water for Beijing, the capital city of China [11,12,13,14]. Therefore, the water quantity and quality of the reservoir is vitally important. Since the 1980s, especially after 2000, large-scale soil conservation measures, such as terrace, fish scale pit, level bench, contour tillage, vegetation plantation, etc., have been implemented in this region by the government and local farmers [15,16]. In this region, drinking water is the priority for local daily life and economic development [17,18,19]. Due to the important role of the Miyun Reservoir, many studies have been performed in this region. For example, Qiu et al. [13,14] evaluated the effect of land use management practices on river discharge, indicating that climate variability, especially precipitation and temperature, had great effects on runoff, sediment yield, and nutrient losses, and significantly affected the efficiency of best management practices. Yan et al. [19] also predicted future climate change impacts on streamflow and nitrogen exports into the reservoir. Tang et al. [20] evaluated land use impacts on nonpoint pollution in this catchment, and the results showed that land use change decreased TN and TP by 39.1% and 23.7%, respectively. Some other studies also reported the variations of streamflow and pollutants entering the downstream Miyun Reservoir [21,22,23,24]. Regrettably, most of these studies were done at a catchment scale. Furthermore, although the impacts of soil conservation measures, rainfall characteristics, and slope gradient and length on soil loss have thoroughly been documented in the world [8], the effects of different soil conservation measures on soil and water losses from slopes in the Miyiun Reservoir catchment were less studied, with few exceptions [25,26]. Soil conservation measures can not only control soil loss, but also intercept runoff on slopes. If more water was intercepted on slopes by the implemented measures, the water quantity in the downstream reservoir could be reduced, threatening the safety of the drinking water for Beijing. However, this kind of study was done less, although a sharp decrease in water quantity has been found in the Yellow River basin and in the Miyun Reservoir basin [15]. Therefore, did these soil conservation measures influence water resource use? Are they reasonably implemented in this water priority region? Sufficiently understanding the impact of each kind of soil conservation measure on soil and water loss is urgently required to guarantee people’s drinking water safety for downstream catchments.
Therefore, runoff and soil loss data from 22 runoff plots in a catchment upstream of the Miyun Reservoir in Beijing were used to (i) identify the characteristics of runoff and soil loss on different slopes, (ii) explore the effects of slope, rainfall characteristics, and soil conservation measures on soil and water losses, and (iii) give some suggestions to reasonably implement soil conservation measures in this region and similar regions of the world. This study can help guide the implementation of soil conservation measures and reduce the loss of drinking water in the study region.
Section 2 describes the study area, runoff plot characteristics, data collection methods, and treatment. Section 3 gives runoff and soil loss characteristics and their responses to rainfall, slope length and gradient, and soil conservation measures. Some suggestions were also given in this section. Section 4 concludes this study.

2. Materials and Methods

2.1. Description of the Study Area

The present study was done in the Shixia catchment (117°4′30″ E and 40°34′40″ N), which is located upstream of the Miyun Reservoir in Beijing, Northern China (Figure 1). This catchment covers an area of 33 km2, ranging from 130 m to 390 m a.s.l. Gneiss is the main lithology, scattered with granite and limestone. This region has a temperate territorial monsoon climate, with a mean annual rainfall amount of 475 mm, above 70% of which occurs in May–September.
In the study region, the main vegetation is Robinia pseudoacacia, Platycladus orientalis, and economic forest, and the main crops are corn and wheat. Around 60.9% of the sandy loam soil contains particles larger than 0.05 mm in diameter, 18.22% contains particles 0.05–0.005 mm in diameter, 9.22% contains particles 0.005–0.001 mm in diameter, and 14.88% contains particles less than 0.001 mm in diameter [25].

2.2. Runoff Plot and Data Collection

In the Shixia catchment, there were 22 runoff plots in total (Figure 1). The boundaries of each plot were made of bricks and cement. The height of the plot boundary is around 30 cm above ground, with a 40-cm deep insertion into the soils to prevent runoff from leaving or entering the plot. Among the 22 plots, one plot is 5 m long, two plots are 20 m long, and the other plots are 10 m long. The slope degrees of these plots vary from 3.5° to 27°. These plots include five cultivated plots, four bare plots, and thirteen plots with different soil conservation measures. The soil conservation measures implemented on these plots include contour tillage, fish scale pit, level bench, terrace, forest, grass, shrub, and their combinations. Corn was planted on the cultivated plots. Detailed information is provided in Table 1.
In 2014–2020, runoff and sediment discharges were collected with a nine-hole diversion bucket and a tank at the end of each plot. After each rainfall event, the runoff amount was measured. Evenly mixed water and sediment samples were also collected with 1000 mL flasks and transported to a laboratory where sediment concentrations were determined with a dry method [15].
Surface runoff depth (H; mm) for each rainfall event was then calculated with the total rainfall amount and the area of the runoff plot. The event soil loss rate (SLR) was also calculated using the total runoff amount, sediment concentration, and runoff plot area. Annual H and annual SLR were then obtained through summing the events for each plot.
A rain gauge and rain barrel near the plots were used to record rainfall information. In the present study, an erosive rainfall was defined as the one that induced erosion on any of the plots, as proposed by Zhu and Zhu (2014). For each erosive rainfall event, rainfall duration (RD), rainfall amount (P), mean rainfall intensity (Im), and maximum rainfall intensities at 30 min (I30) and at 60 min (I60) were obtained from rainfall process data.

2.3. Data Treatment and Statistical Analysis

In the present study, Pearson correlation analysis and regression analysis were used to determine the relationships between H, SLR, and the five rainfall eigenvalues. Fisher’s protected least significant difference test was also used to compare the means of Hs and SLRs from the plots. Treatments were considered significant if p < 0.05. All of the statistics were conducted using SPSS version 14.0 for Windows.

3. Results and Discussion

3.1. Rainfall Characteristics

The mean annual Ps in 2014–2020 ranged from 385 mm in 2020 to 532.2 mm in 2018, with an average of 475.1 mm, a coefficient of variation of 0.11, and a skewness of −0.86. The annual erosive Ps varied from 239.2 mm in 2020 to 434.7 mm in 2018, with an average of 297.1 mm, or 62.3% of the annual mean P (Figure 2 and Table 2). In 2014–2020, over 50% of the annual Ps produced runoff on the plots. These percentages were higher than those on the Chinese Loess Plateau reported by Liang et al. (2020).
In 2014–2020, there were 66 erosive rainfall events. The event Ps ranged from 4.8 to 108.1mm, with an average of 31.5 mm. The Ps in the range of 10–20 mm occupied over one third of the total number of rainfall events. All of these erosive rainfall events ranged from 0.3 to 32.3 h in RD, with an average of 6.9 h, from 1.0 to 54.0 mm h−1 in Im with an average of 9.5 mm h−1, and from 3.6 to 61.7 mm h−1 and 6.0 to 64.2 mm h−1 in I60 and I30, respectively (Figure 3). Under future climate change, more extreme rainfall events could occur [13], inducing more severe soil erosion in the study region.

3.2. Runoff and Soil Loss

Annual and event Hs varied widely among the 22 runoff plots (Figure 4a and Table 3). The mean annual Hs in 2014–2020 ranged from 0 on plots 8 and 10 to 100.61 mm on plot 22. The mean annual Hs were less than 18.03 mm on most of the plots, which were significantly less than those on the bare and cultivated plots 1, 3, 4, 6, 21, and 22. Event Hs also differed greatly. The maximum event H was 10.67 mm on plot 22, and the minimum H was zero on plots 8 and 10 (Table 3). Most of the plots had a mean event H less than 2.00 mm. The event Hs on plots 1, 3, 4, 6, 21, and 22 were significantly higher than those on other plots.
Annual SLR also varied from 0 on plots 8 and 10 to 6006.84 t km−2 year−1 on plot 21 (Figure 4b), and annual SLRs on the plots 1, 3, 4, 6, 21, and 22 were higher than 2500 t km−2 year−1, which is above twelve times the tolerate value (i.e., 200 t km−2 year−1) in the region [15]. In the Miyun Reservoir catchment, these lands are around 600 km2, occupying around 5% of the catchment area. Although this percentage is smaller, it should also be paid attention to, because it is still an important pollutant source [13,14], implying that these kinds of lands should have soil conservation measures implemented. In contrast, the annual SLRs from plots 2, 5, and 7–19 were much lower. In terms of event SLRs, the maximum SLR was up to 637 t km−2 event−1 from plot 21, and plots 1, 3, 4, 6, and 22 also had high event SLRs. Conversely, other plots had a much lower event SLR than 25 t km−2 event−1, and plots 8 and 10 did not yield sediment (Table 3).

3.3. Impact of the Slope Gradient and Slope Length

The effect of the slope gradient on runoff and SLR is a controversial issue. As the slope gradients increase, surface runoff and soil loss can increase, decrease, and/or change little, resulting from the influence of complex factors, such as a forming or disruption of soil crust [27,28], rill and/or ephemeral gully development [29,30], differential soil cracking [31], and/or ponding depth [32,33], etc. Sometimes, a threshold maximum value of H or SLR exists with increasing slope gradients, resulting from the interaction of runoff generation and its sediment carrying capacity [28,34]. In the present study, both annual H and annual SLR increased first and then decreased, with a threshold annual H of 100.1 mm and a threshold annual SLR of 4748.29 t km−2 year−1 occurring at the 14.4° slope (Figure 5). A number of studies [35,36,37,38] found that the threshold slope gradients of H and SLR varied from 25° to 50°, which are much steeper than the one used in this study. The threshold slope gradient of 14.4° of annual H and annual SLR can be explained by the impact of soil conservation measures. In the present study, the bare plots with a slope degree of 14.4° had the highest annual H and annual SLR, whereas other plots with slope degrees lower and/or higher than 14.4° were implemented with soil conservation measures, and yielded lower annual Hs and annual SLRs. For the plots with the same treatment, the changes of Hs and SLRs with the slope gradient varied depending on plot treatments. For example, higher annual Hs and SLRs occurred on contour tillage plot 6 than those of contour tillage plots 18 and 20. Similar phenomenon also occurred for the terraced orchard plots 2 and 7 (Figure 4). This phenomenon could be attributed to the decrease in depressional storage and ponding depth on the plots with higher slope degrees [33,39]. However, for the shrub plots (i.e., plots 5, 9, 11, 13, 14, and 19), a threshold slope seemed to exist, whereas Hs and SLRs on plots 21 and 22 with a lower slope gradient were much higher than that of plot 4, resulting from the effect of the slope length.
Runoff and soil loss also present different changing patterns with increasing slope length [40,41,42]. In the present study, only three slope lengths (i.e., 5 m, 10 m, and 20 m) were included (Table 1). The scattered dots in Figure 6 illustrate that annual Hs and SLRs increased or decreased when the slope lengths increased from 5 m to 20 m, resulting from the impact of soil conservation measures. Further analysis was done when their impacts were excluded. The annual H decreased when the slope length increased from 5 m to 10 m for plots 11 and 14, which had similar coverage rates and slope degrees (Table 1). This is because the time and volume of rainfall per unit area required before the runoff starts are higher in the longer plot than the smaller one [43]. However, for the bare plots 3, 4, 21, and 22, the annual H increased with increasing slope length. The difference could result from different earth surface conditions. According to study results by Mermut et al. [44] and Fang et al. [28], the soil texture in the study region made the soil crust easily form on the bare plot, which reduced the infiltration rate and increased the runoff generation capacity [33]. Although soil crust was not detected in the present study, it has been verified by previous studies through the use of natural rainfall [33] and rainfall simulation experiments [28]. Given the same land surface, annual SLRs increased with increasing slope length. For example, the mean annual SLRs of bare plots 21 and 22 were higher than those of bare plots 3 and 4, resulting from an increased surface runoff amount with increasing slope length. In contrast, a higher SLR on brush plot 14 than that on brush plot 11 could be explained by increased sediment concentration in runoff due to the increased concentrated flow on a longer plot [28].

3.4. Impact of Rainfall Characteristics

Rainfall is the external force inducing surface runoff and soil loss. In the present study, five rainfall eigenvalues were used to explore the impact of rainfall characteristics on event H and SLR. Pearson correlation coefficients indicated that event H was significantly correlated with event P and rainfall intensity (I30 and/or I60) at the 0.01 level, excluding the terraced plots 2 and 7, and the densely vegetated plot 12 (Table 4), the event H and SLR of which, however, were significantly correlated with P and RD. In respect to event SLR, event SLR was significantly correlated with P, I30, and/or both of them on most of the plots. Different responses of H and SLR from the plots to different rainfall eigenvalues can result from the protection capacities of the measures [45,46,47]. For instance, the vegetated measures can effectively protect the splash impact of rainfall intensity, resulting in insignificant correlations of H and rainfall intensity for the vegetated plots.
Similar to the studies by Fang et al. [33] and Zhu and Zhu [7], a large percent of annual runoff and soil loss occurred during several rainfall events. In the present study, we compared soil loss by the maximum annual erosive event to annual soil loss through four exampled plots, which represented different soil conditions (i.e., plots 1–3 and 14) in Figure 7. Soil loss from cultivated plot 1 by the maximum annual soil erosive event in 2014–2020 ranged from 1658.67 to 5104.98 t km–2, which contributed to a mean value of 40.50% of the annual total soil loss, ranging from 30.72% in 2019 to 61.13% in 2018. The magnitude and contributions of the maximum annual erosive event on bare plot 3 were similar to those of plot 1. However, much higher contributions of the maximum annual erosive event occurred on plots 2 and 14, with terrace and shrub measures, respectively. The maximum contributions of the maximum annual erosive events were up to 85% on plot 2 and 89% on plot 14, although their maximum soil loss rates were much less than those on the bare and cultivated plots.

3.5. Impact of Soil Conservation Measures

In order to study the impact of soil conservation measures, ten kinds of land use and soil conservation measures were obtained by using the 22 plots in Table 1 (Figure 8). Among the ten treatments, mean annual Hs differed greatly, and changed from 85.80 mm on bare plots, 72.83 mm on the cultivated plot without any measure, 67.74 mm on the cultivated plots with contour tillage, 11.49 mm on the grass plots, followed by the shrub plots, orchard plots with terrace, and other plots. Noticeably, no runoff occurred on plot 8 with forest plus grass measures and plot 10 with forest and fish scale pit measures. In comparison to the bare plots, the mean runoff reduction efficiencies were above 86% for the plots with soil conservation measures, excluding plots 1 and 6. Based on this result, the large-scale implementation of soil conservation measures since 2000 in the study region could, to some extent, explain the decreased discharge into the Miyun Reservoir, which agrees with the former studies at a catchment scale [13,14,15]. In order to ensure water quantity of the Miyun Reservoir, soil conservation measures should be water-saving, which allows more runoff downstream to the reservoir.
Annual SLRs presented the same changing patterns. The bare plots had the largest SLR, with an average of 4325.79 t km−2 year−1, followed by the cultivated plot with no soil conservation measures, with an average of 3205.23 t km−2 year−1, and the cultivated plot with contour tillage, with an average of 1067.62 t km−2 year−1. In contrast, other kinds of plots had much lower mean annual SLRs. The mean annual SLR was 109.39 t km−2 year−1 for the grass plots, 74.16 t km−2 year−1 for the orchard plot with terrace, 55.97 t km−2 year−1 for the shrub plots, and 16.3 t km−2 year−1 for the forest plus shrub plots. No sediment ran out of the plots with forest plus fish scale pit and/or grass measures. Noticeably, almost all of the plots with soil conservation measures had higher soil loss reduction efficiencies of above 95%. The effective reductions of soil and water losses by vegetated measures were consistent with a number of studies [47,48,49,50,51]. Consistent with previous studies [26,47], terrace was effective in reducing runoff and soil loss. In the present study, soil loss from cultivated plot 17 with terrace was nearly 0 (i.e., 3.47 t km−2 year−1), and no runoff and soil loss occurred on plot 10, which had terrace and forest measures implemented, although its vegetation coverage was not very high. This implies that the combination of vegetation and terrace measures was more effective than a single measure in reducing soil and water losses [8]. In the present study, the mean annual SLR on plot 6 with a slope degree of 14.4° was up to 2541 t km−2 year−1, which is twelve times above the tolerable rate, although contour tillage was also implemented on this plot (Figure 3). In contrast, contour tillage can effectively reduce soil loss from plots 18 and 20, which have lower slope degrees (Table 1). This implies that contour tillage is not reasonable when it is implemented on steep slopes in the study region.
The implemented soil conservation measures also significantly interfered with the relationships between event H and event SLR. Table 5 demonstrated that, except for the bare and cultivated plots, the SLR cannot be fitted at the 0.01 level.

4. Conclusions

The latest seven years (2014–2020) of monitored data from 22 runoff plots in the upstream catchment of the Miyun Reservoir in Beijing, Northern China were employed to study soil and water loss characteristics on slopes and their responses to slope, rainfall, and soil conservation measures, and some findings and implications were obtained. The bare land suffered the severest soil loss, with a mean annual SLR of 4325 t km−2 year−1. Cultivated lands without any measures also experienced severe soil loss, with an annual SLR above 3205 t km−2 year−1. Contour tillage cannot effectively control soil loss on the cultivated land with steep slopes. Soil conservation measures, such as forest, shrub, grass, terrace, and fish scale pits, as well as their combinations, can effectively reduce soil and water losses, with soil and water loss reduction efficiencies higher than 95% and 86%, respectively.
To protect drinking water safety, water quantity is a priority in the study region, and water-saving soil control measures should be considered to resolve the contradiction of soil loss control and water shortage in the downstream reservoir. The vegetation coverage of forest, shrub, and grass should be sparser, which allows more runoff to run downslope and simultaneously and effectively control soil loss. The size of the terrace on gentle slopes should be carefully designed. Bare and cultivated lands suffered the most severe soil erosion, with SLRs above 2732 t km−2 year−1, and soil conservation measures should be implemented. Cultivated land with contour tillage still suffered a very high SLR (i.e., 2541 t km−2 year−1), and terrace or other vegetation measures (e.g., grass strips) should be implemented on this kind of land, because these types of measures can greatly reduce SLRs in the study region (Figure 4b).
Future studies can be conducted using physically based soil erosion models, such as the Water Erosion Prediction Plan (WEPP) and Soil and Water Assessment Tool (SWAT), to evaluate the effect of soil conservation measures on different slopes on downstream water resources.

Funding

This work was financially funded by the Beijing Natural Science Foundation (grant number 8202045) and the National Natural Science Foundation of China (grant number 41977066).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in the present study are available on request from the corresponding author.

Acknowledgments

Special thanks are owed to the two reviewers for their invaluable suggestions to the improvement of the manuscript.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Xiong, M.Q.; Sun, R.H.; Chen, L.D. Effects of soil conservation techniques on water erosion control: A global analysis. Sci. Total Environ. 2018, 645, 753–760. [Google Scholar] [CrossRef] [PubMed]
  2. Damnati, B.; Ibrahimi, S.; Radakovitch, Q. Quantifying erosion using 137Cs and 210Pb in cultivated soils in three Mediterranean watershed: Synthesis study from El Hachef, Raouz and Nakhla (North West Morocco). J. Afr. Earth Sci. 2013, 79, 50–57. [Google Scholar] [CrossRef]
  3. Anache, J.A.A.; Wendland, E.C.; Oliveira, P.T.S.; Flanagan, D.C.; Nearing, M.A. Runoff and soil erosion plot-scale studies under natural rainfall: A meta-analysis of the Brazilian experience. Catena 2017, 152, 29–39. [Google Scholar] [CrossRef]
  4. Borrelli, P.; Robinson, D.A.; Fleischer, L.R.; Lugato, E.; Ballabio, C.; Alewell, C.; Meusburger, K.; Modugno, S.; Schutt, B.; Ferro, V.; et al. An assessment of the global impact of 21st century land use change on soil erosion. Nat. Commun. 2017, 8, 1–13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Wuepper, D.; Borrelli, P.; Finger, R. Countries and the global rate of soil erosion. Nat. Sustain. 2020, 3, 51–55. [Google Scholar] [CrossRef]
  6. Troeh, F.R.; Hobbs, J.A.; Donahue, R.L. Soil and Water Conservation for Productivity and Environmental Protection; Prentice-Hall, Inc.: Englewood Cliffs, NJ, USA, 1980. [Google Scholar]
  7. Zhu, T.X.; Zhu, A.X. Assessment of soil erosion and conservation on agricultural sloping lands using plot data in the semi-arid hilly loess region of China. J. Hydrol. Reg. Stud. 2014, 2, 69–83. [Google Scholar] [CrossRef] [Green Version]
  8. Zhao, J.L.; Yang, Z.Q.; Gover, G. Soil and water conservation measures reduce soil and water losses in China but not down to background levels: Evidence from erosion plot data. Geoderma 2019, 337, 729–741. [Google Scholar] [CrossRef]
  9. Wolka, K.; Mulder, J.; Biazin, B. Effects of soil and water conservation techniques on crop yield, runoff and soil loss in Sub-Saharan Africa: A review. Agric. Water Manag. 2018, 207, 67–79. [Google Scholar] [CrossRef]
  10. Maetens, W.; Poesen, J.; Vanmaercke, M. How effective are soil conservation techniques in reducing plot runoff and soil loss in Europe and the Mediterranean? Earth-Sci. Rev. 2012, 115, 21–36. [Google Scholar] [CrossRef] [Green Version]
  11. He, Y.Y.; Wang, X.Y.; Duan, S.H. Characteristics of runoff and sediment during individual rainfall in upper area of Miyun Reservoir. Transac. Chin. Soc. Agric. Eng. 2015, 31, 134–141. [Google Scholar]
  12. Han, L.F.; Gao, B.; Lu, J.; Zhou, Y.; Xu, D.Y.; Gao, L.; Sun, K. Pollution characteristics and source identification of trace metals in riparian soils of Miyun Reservoir, China. Ecotox. Environ. Saf. 2017, 144, 321–329. [Google Scholar] [CrossRef]
  13. Qiu, J.L.; Shen, Z.Y.; Hou, X.S.; Xie, H.; Leng, G.Y. Evaluating the performance of conservation practices under climate change scenarios in the Miyun Reservoir Watershed, China. Ecol. Eng. 2020, 143, 105700. [Google Scholar] [CrossRef]
  14. Qiu, J.L.; Shen, Z.Y.; Huang, M.Y.; Zhang, X.S. Exploring effective best management practices in the Miyun reservoir watershed, China. Ecol. Eng. 2018, 123, 30–42. [Google Scholar] [CrossRef]
  15. Li, Z.J. Impacts of Soil and Water Conservation Measures on Annual Runoff in the Chao River Basin. Ph.D. Thesis, Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing, China, 2007. [Google Scholar]
  16. Liu, B.Y.; Bi, X.G.; Fu, S.H. Beijing Soil Loss Equation; Science Press: Beijing, China, 2010. [Google Scholar]
  17. Gao, Y.C.; Yao, Z.J.; Liu, B.Q.; Lv, A.F. Evolution trend of Minyun Reservoir inflow and its motivating factors analysis. Prog. Geogr. 2002, 21, 546–553. [Google Scholar] [CrossRef]
  18. Ma, H.; Yang, D.W.; Tan, S.K.; Gao, B.; Hu, Q.F. Impact of climate variability and human activity on streamflow decrease in the Miyun Reservoir catchment. J. Hydrol. 2010, 389, 317–324. [Google Scholar] [CrossRef]
  19. Yan, T.Z.; Bai, J.W.; Arsenio, T.; Shen, Z.Y. Future climate change impacts on streamflow and nitrogen exports based on CMIP5 projection in the Miyun Reservoir Basin, China. Ecohydrol. Hydrobiol. 2019, 19, 266–278. [Google Scholar] [CrossRef]
  20. Tang, L.H.; Yang, D.W.; Hu, H.P.; Gao, B. Detecting the effect of land-use change on streamflow, sediment and nutrient losses by distributed hydrological simulation. J. Hydrol. 2011, 409, 172–182. [Google Scholar] [CrossRef]
  21. Xu, E.Q.; Zhang, H.Q. Aggregating land use quantity and intensity to link water quality in bupper catchment of Miyun Reservoir. Ecol. Indic. 2016, 66, 329–339. [Google Scholar] [CrossRef]
  22. Zhang, L.; Zou, Z.H.; Shan, W. Development of a method for comprehensive water quality forecasting and its application in Miyun reservoir of Beijing, China. J. Environ. Sci. 2016, 56, 240–246. [Google Scholar] [CrossRef] [PubMed]
  23. Xu, W.; Cai, Y.P.; Yin, X.N.; Yan, H.; Zhang, L. Soil nitrogen supply capacity as an indicator of sustainable watershed management in the upper basin of Miyun Reservoir. Ecol. Indic. 2018, 92, 221–227. [Google Scholar] [CrossRef]
  24. Xu, R.H.; Cai, Y.Q.; Wang, X.; Li, C.H.; Liu, Q.; Yang, Z.F. Agricultural nitrogen flow in a reservoir watershed and its implications for water pollution mitigation. J. Clean. Product. 2020, 267, 122034. [Google Scholar] [CrossRef]
  25. Wang, X.Y.; Hu, Q.J.; Zhu, F.Y.; He, W.; Wang, X.F.; Wang, G. Loss of nutrients in runoff and soil in the watershed of Miyun Reservoir-as an example of Shixia experimental area. J. Cap. Norm. Univ. 2001, 22, 79–85. [Google Scholar]
  26. Gao, Y. Study on the effect of different soil conservation measures on soil erosion in Shixia catchment. Dev. Inno. 2018, 29, 255–256. [Google Scholar]
  27. Poesen, J. The influence of slope angle on infiltration rate and Hortonian overland flow. Geomorphology 1984, 49, 117–131. [Google Scholar] [CrossRef]
  28. Fang, H.Y.; Sun, L.Y.; Tang, Z.H. Effects of rainfall and slope on runoff, soil erosion and rill development: An experimental study using two loess soils. Hydrol. Process. 2015, 29, 2649–2658. [Google Scholar] [CrossRef]
  29. Bryan, R.B.; Poesen, J. Laboratory experiments on the influence of slope length on runoff, percolation, and rill development. Earth Surf. Process. Landf. 1989, 4, 211–231. [Google Scholar] [CrossRef]
  30. Slattery, M.; Bryan, R.B. Hydraulic condition for rill incision under simulated rainfall: A laboratory experiment. Earth Surf. Process. Landf. 1992, 17, 127–146. [Google Scholar] [CrossRef]
  31. Govers, G. A field study on topographical and topsoil effects on runoff generation. Catena 1990, 18, 91–111. [Google Scholar] [CrossRef]
  32. Fox, D.M.; Bryan, R.B.; Price, A.G. The influence of slope angle on final infiltration rate for interrill conditions. Geoderma 1997, 80, 181–194. [Google Scholar] [CrossRef]
  33. Fang, H.Y.; Cai, Q.G.; Chen, H.; Li, Q.Y. Temporal changes in suspended sediment transport in a gullied loess basin: The lower Chabagou Creek on the Loess Plateau in China. Earth Surf. Process. Landf. 2008, 33, 1977–1992. [Google Scholar] [CrossRef]
  34. Zheng, Z.C.; Qin, F.; Li, T.X. Changes in soil surface microrelief of purple soil under different slope gradients and its effects on soil erosion. Trans. Chin. Soc. Agric. Eng. 2015, 31, 168–175. [Google Scholar]
  35. Guo, J.Z. Discussion on the issue of the relations of slope gradient and runoff amount and sediment production. Yel. Riv. 1958, 3, 47–49. [Google Scholar]
  36. Jin, C.X. A theoretical study on critical erosion slope gradient. Acta Geogr. Sin. 1995, 50, 234–239. [Google Scholar]
  37. Cao, W.H. Slope thresholds of soil erosion. Bull. Soil Water Conserv. 1993, 13, 1–5. [Google Scholar]
  38. Liu, X.B.; Zhang, S.L.; Zhang, X.Y.; Ding, G.W.; Cruse, R.M. Soil erosion control practices in Northeast China: A mini-review. Soil Tillage Res. 2011, 17, 44–48. [Google Scholar] [CrossRef]
  39. De Ploey, J.; Savat, J.J.; Moeyersons, J. The differential impact of some soil loss factors on flow, runoff creep and rainwash. Earth Surf. Process. Landf. 1976, 1, 151–1616. [Google Scholar] [CrossRef]
  40. Luo, L.X. Primary result of water erosion amount from gully slopes through field investigation in Lianshi Gou catchment in Gansu province. Geogr. Data 1958, 2, 111–118. [Google Scholar]
  41. Zheng, F.L.; Tang, K.L.; Zhou, P.H. Study on factors affecting rill erosion on cultivated slope land. Acta Pedol. Sin. 1989, 26, 109–116. [Google Scholar]
  42. Kong, Y.P.; Zhang, K.L.; Tang, K.L. Impacts of slope length on soil erosion process under simulated rainfall. J. Soil Water Conserv. 2001, 15, 17–24. [Google Scholar]
  43. Joel, A.; Messing, I.; Seguel, O.; Casanova, M. Measurement of surface, water runoff from plots of two different sizes. Hydrol. Process. 2002, 16, 1467–14786. [Google Scholar] [CrossRef]
  44. Mermut, A.R.; Luk, S.H.; Römkens, M.J.M.; Poesen, J.W.A. Micromorphological and mineralogical components of surface sealing in loess soils from different geographic regions. Geoderma 1995, 66, 71–84. [Google Scholar] [CrossRef]
  45. Jiao, J.Y.U.; Zou, H.Y.; Jia, Y.F.; Wang, N. Research progress on the effects of soil erosion on vegetation. Acta Ecol. Sin. 2009, 29, 85–91. [Google Scholar] [CrossRef]
  46. Chen, H.; Zhang, X.P.; Alba, M.; Yan, R.; Ren, Q.F.; Ren, Z.Y.; Yang, Y.H.; Zhao, W.H.; Lin, P.F.; Liu, B.Y.; et al. Effects of vegetation and rainfall types on surface runoff and soil erosion on steep slopes on the Loess Plateau, China. Catena 2018, 170, 141–1496. [Google Scholar] [CrossRef]
  47. Liang, Y.; Jiao, J.Y.; Tang, B.Z.; Cao, B.T.; Li, H. Response of runoff and soil erosion to erosive storm events and vegetation restoration on abandoned slope farmland in the Loess Plateau region, China. J. Hydrol. 2020, 584, 1246946. [Google Scholar] [CrossRef]
  48. Kang, S.Z.; Zhang, L.; Song, X.Y.; Zhang, S.H.; Liu, X.Z.; Liang, Y.L.; Zheng, S.Q. Runoff and sediment los response to rainfall and land use in two agricultural catchments on the Loess Plateau of China. Hydrol. Process. 2001, 15, 977–9886. [Google Scholar] [CrossRef]
  49. Cerdan, O.; Le, B.Y.; Couturier, A.; Bourennane, H.; Souchere, V. Rill erosion on cultivated hillslopes during two extreme rainfall events in Normandy, France. Soil Tillage Res. 2002, 67, 99–1086. [Google Scholar] [CrossRef]
  50. Wei, W.; Chen, L.D.; Fu, B.J.; Lu, Y.H.; Gong, J. Responses of water erosion to rainfall extremes and vegetation types in a loess semiarid hilly area, NW China. Hydrol. Process. 2009, 23, 1780–17916. [Google Scholar] [CrossRef]
  51. Wang, Z.J.; Jiao, J.Y.; Rayburg, S.; Wang, Q.L.; Su, Y. Soil erosion resistance of “Grain for Green” vegetation types under extreme rainfall conditions on the Loess Plateau. China. Catena 2016, 141, 109–116. [Google Scholar] [CrossRef]
Figure 1. Location of the 22 runoff plots. Note: the lower figure was made from a Google image.
Figure 1. Location of the 22 runoff plots. Note: the lower figure was made from a Google image.
Water 13 00853 g001
Figure 2. Annual erosive and non-erosive rainfall and the percentages of erosive rainfall amounts in 2014–2020.
Figure 2. Annual erosive and non-erosive rainfall and the percentages of erosive rainfall amounts in 2014–2020.
Water 13 00853 g002
Figure 3. Frequencies of some eigenvalues of the erosive rainfall events between 2014–2020: (a) event rainfall amount, (b) rainfall duration, (c) mean rainfall intensity, (d) maximum 60-min rainfall intensity, and (e) maximum 30-min rainfall intensity.
Figure 3. Frequencies of some eigenvalues of the erosive rainfall events between 2014–2020: (a) event rainfall amount, (b) rainfall duration, (c) mean rainfall intensity, (d) maximum 60-min rainfall intensity, and (e) maximum 30-min rainfall intensity.
Water 13 00853 g003
Figure 4. Annual runoff depth (a) and annual soil loss rate (b) and their reduction efficiencies by soil conservation measures in comparison to the mean values of bare plots 3 and 4 in 2014–2020 for the 22 runoff plots. Notes: the bars represent standard error. The average values on the columns with the same letter are not significantly different at p = 0.05, as determined by Fisher’s protected least test.
Figure 4. Annual runoff depth (a) and annual soil loss rate (b) and their reduction efficiencies by soil conservation measures in comparison to the mean values of bare plots 3 and 4 in 2014–2020 for the 22 runoff plots. Notes: the bars represent standard error. The average values on the columns with the same letter are not significantly different at p = 0.05, as determined by Fisher’s protected least test.
Water 13 00853 g004
Figure 5. Change patterns of annual runoff depth (a) and annual soil loss rate (b) with increasing slope degree in 2014–2020 for the 22 runoff plots. Note: the plots in black letters represent bare plots, and those in red letters present shrub plots in Table 1.
Figure 5. Change patterns of annual runoff depth (a) and annual soil loss rate (b) with increasing slope degree in 2014–2020 for the 22 runoff plots. Note: the plots in black letters represent bare plots, and those in red letters present shrub plots in Table 1.
Water 13 00853 g005
Figure 6. Annual runoff depth (a) and soil loss rate (b) for the 22 runoff plots in 2014–2020 with slope length.
Figure 6. Annual runoff depth (a) and soil loss rate (b) for the 22 runoff plots in 2014–2020 with slope length.
Water 13 00853 g006
Figure 7. Soil loss in the maximum annual erosive event and annual soil loss on four types of land use and soil conservation measures: plot 1 (a), plot 2 (b), plot 3 (c), and plot 14 (d). The information of the plots was given in Table 1.
Figure 7. Soil loss in the maximum annual erosive event and annual soil loss on four types of land use and soil conservation measures: plot 1 (a), plot 2 (b), plot 3 (c), and plot 14 (d). The information of the plots was given in Table 1.
Water 13 00853 g007
Figure 8. Annual runoff depth (a) and soil loss rate (b) under different land use and soi conservation measures for the study plots in 2014–2020. Note: On the x-axis, bare represents plots 3, 4, 21, and 22, cultivated represents plot 1, cultivated + contour represents plots 6, 18, and 20, and cultivated + terrace represents plot 17. Orchard represents plot 2 and 7, shrub represents plots 5, 9, 11, 13, 14, and 19, forest + grass represents plot 8, forest + fish-scale pit represents plot 10, forest + shrub represents plot 12, and grass represents plots 15 and 16 in Table 1.
Figure 8. Annual runoff depth (a) and soil loss rate (b) under different land use and soi conservation measures for the study plots in 2014–2020. Note: On the x-axis, bare represents plots 3, 4, 21, and 22, cultivated represents plot 1, cultivated + contour represents plots 6, 18, and 20, and cultivated + terrace represents plot 17. Orchard represents plot 2 and 7, shrub represents plots 5, 9, 11, 13, 14, and 19, forest + grass represents plot 8, forest + fish-scale pit represents plot 10, forest + shrub represents plot 12, and grass represents plots 15 and 16 in Table 1.
Water 13 00853 g008
Table 1. Characteristics of the 22 plots. Note: all of the treatments were applied for all seven years (i.e., 2014–2020); the vegetation coverage in the plots was obtained through visual estimation.
Table 1. Characteristics of the 22 plots. Note: all of the treatments were applied for all seven years (i.e., 2014–2020); the vegetation coverage in the plots was obtained through visual estimation.
PlotSlope
(°)
Length
(m)
Area
(m2)
Soil Depth
(cm)
Land UseSoil Conservation Measure
116.5105030Cultivated (corn)-
216.5105030Orchard (chestnut)Terrace (size: 3-m wide)
316.5105030Bare-
414.4105020Bare-
514.4105020ShrubVitex negundo (coverage: 40–60%)
614.4105020Cultivated (corn)Contour tillage
79.4105020Orchard (hawthorn)Terrace (size: 3-m wide)
827105020Forest and grassRobinia and grass (coverage: 90–100%)
927105030ShrubVitex negundo (coverage: 40–50%)
1027105020ForestRobinia pseudoacacia (coverage: 40%); fish-scale pits (size: 1 × 1.5 m2)
1119.352515ShrubVitex negundo (coverage: 45–60%)
1217.1105015Forest and shrubPlatycladus orientalis and Robinia pseudoacacia (coverage: 90%); fish-scale pits (size: 1 × 1.5 m2)
1318.6105015ShrubVitex negundo (coverage: 40–50%)
1419.7105015ShrubVitex negundo (coverage: 40–50%)
1519105015GrassAlfalfa (coverage: 50–60%)
1619105015GrassAlfalfa (coverage: 10–15%)
173.5105060Cultivated (corn)Terrace (size: 4-m wide)
183.5105060Cultivated (corn)Contour tillage
196.3105020ShrubVitex negundo (coverage: 45–60%)
206.3105020Cultivated (corn)Contour tillage
2114.42010020Bare-
2214.42010020Bare-
Table 2. Statistical characteristics of annual (erosive) rainfall and the percent of erosive to total rainfall amount in 2014–2020.
Table 2. Statistical characteristics of annual (erosive) rainfall and the percent of erosive to total rainfall amount in 2014–2020.
Annual P (mm)Annual Erosive P (mm)Percent of Mean Annual Erosive Rainfall Amount to the Annual Mean P (%)
Mean475.1297.162.3
Range385.0–532.2239.1–434.750.6–81.7
Std.52.962.68.8
CV0.110.210.14
Skewness-0.861.741.57
Table 3. Characteristics of threshold rainfall, runoff, and soil loss characteristics for the 22 plots in 2014–2020.
Table 3. Characteristics of threshold rainfall, runoff, and soil loss characteristics for the 22 plots in 2014–2020.
PlotTRA
(mm) *
MAP
(mm)
Accu. P
(mm)
RGT TimesRunoff Depth (mm)Soil Loss Rate (t km−2 year−1)
Min.Max.MeanSt.DMin.Max.MeanSt.D
14.8280.391962.7620.0029.657.73b0.770.002278.23339.95cd56.10
221.8107.31751.1140.0018.080.85c0.320.00364.1615.69e6.38
34.8274.131918.9600.0018.087.99b0.750.002429.35289.81d56.99
48.4268.041876.3570.0018.087.26b0.740.002466.32404.67c66.56
515.195.47668.3170.0011.531.22c0.310.00403.8718.35e7.64
68.4272.711909.0590.0018.087.34b0.720.002687.83269.57d53.04
755.041.36289.540.003.620.14c0.080.001.270.05e0.03
800.000.000.000.000.00c0.000.000.000.00e0.00
946.125.56178.930.0011.970.29c0.200.0028.630.76e0.53
1000.000.000.000.000.00c0.000.000.000.00e0.00
1115.196.19673.3150.0018.081.60c0.430.0037.212.13e0.78
1255.060.29422.050.0017.010.47c0.280.00107.281.73e1.63
1315.196.86678.0140.008.750.80c0.240.0023.130.99e0.44
1412.686.77607.4160.0010.550.98c0.270.00387.0612.08e6.16
1515.175.69529.8100.005.210.53c0.170.0018.780.59e0.31
1615.1134.26939.8210.0014.921.91c0.460.00568.7722.61e10.36
1728.778.30548.180.005.420.36c0.130.007.710.37e0.16
188.4255.691789.8530.0018.086.96b0.760.00130.4016.77e3.09
1915.175.97531.890.009.220.47c0.180.0073.651.30e1.12
208.4264.461851.2550.0031.457.25 b0.800.00920.3653.36e16.54
214.8280.361962.5640.0081.0410.49ab1.740.002534.04637.09a79.97
224.8278.501949.5630.0088.2310.67a1.830.003024.96503.61b75.53
Notes: * indicates the recorded minimum value to generate runoff in the 66 rainfall events. TRA, MAP, Accu. P, and RGT represent the threshold of the runoff amount, mean annual rainfall amount, accumulated rainfall amount, and runoff generation times, respectively. The mean values of runoff depth and soil loss rate on the columns with the same letter were not significantly at the 0.05 level.
Table 4. Pearson coefficient matrix between the five rainfall eigenvalues, event runoff depth (H), and event soil loss rate SLR for 66 rainfall events in 2014–2020.
Table 4. Pearson coefficient matrix between the five rainfall eigenvalues, event runoff depth (H), and event soil loss rate SLR for 66 rainfall events in 2014–2020.
RDPImI30I60
HP10.293 *0.642 **0.0580.546 **0.407 **
P20.344 **0.412 **−0.1170.2120.146
P30.379 **0.601 **0.0080.473 **0.368 **
P40.300 *0.544 **0.0400.477 **0.358 **
P50.1020.306 *0.1790.410 **0.308 *
P60.263 *0.566 **0.1170.651 **0.460 **
P70.522 **0.524 **−0.1210.2060.131
P9−0.100 *0.1680.2090.1660.300 *
P110.0690.300 *0.2110.395 **0.293 *
P120.248 *0.574 **−0.0200.1220.148
P130.0700.260 *0.1980.398 **0.259 *
P14−0.0070.257 *0.254 *0.409 **0.321 **
P150.1960.359 **0.1950.365 **0.250 *
P160.1810.537 **0.1710.514 **0.410 **
P170.399 **0.656 **−0.0190.330 **0.294 *
P180.385 **0.642 **0.0410.520 **0.369 **
P190.269 *0.366 **−0.0250.1310.106
P200.281 *0.643 **0.0730.553 **0.410 **
P210.555 **0.640 **−0.0570.394 **0.274 *
P220.567 **0.682 **−0.0390.403 **0.318 **
SLRP10.0020.310 *0.1880.431 **0.317 **
P20.0510.2290.0890.283 *0.151
P3−0.0140.298 *0.2230.410 **0.366 **
P4−0.1100.265 *0.290 *0.436 **0.388 **
P5−0.1000.0990.1760.2070.167
P60.0820.347 **0.1680.441 **0.335 **
P70.337 **0.408 **−0.0880.1260.023
P9−0.1290.0550.2210.0820.245 *
P11−0.0760.2210.1290.1280.170
P120.460 **0.420 **−0.0770.265 *0.240
P13−0.1540.0000.2100.319 **0.155
P14−0.0240.2010.0700.1990.181
P15−0.0960.0240.2250.2040.097
P16−0.0520.306 *0.1130.1790.225
P170.0040.359 **0.1260.2180.304 *
P18−0.1090.273 *0.1720.290 *0.246
P190.027−0.013−0.064−0.048−0.064
P200.295 *0.498 **0.0790.396 **0.373 **
P210.0350.334 **0.2390.525 **0.408 **
P220.0040.315 **0.284 *0.519 **0.463 **
Note: * represents significance at the 0.05 level, and ** represents significance at the 0.01 level. P, RD, Im, I30, and I60 represent rainfall amount, rainfall duration, mean rainfall intensity, and maximum 30- and 60-min rainfall intensity for each erosive rainfall event.
Table 5. Relationships between event runoff depth (H) and soil loss rate (SLR) (SLR = aHb) for the 22 runoff plots.
Table 5. Relationships between event runoff depth (H) and soil loss rate (SLR) (SLR = aHb) for the 22 runoff plots.
PlotNabF-TestR2Sig.
1627.0041.59857.4890.4890.000
21419.7230.3090.1780.0150.680
3607.7421.30028.8740.3000.000
45713.7101.31042.9280.4340.000
5171.1732.0567.3710.3290.016
6593.5791.80487.4120.6050.000
740.2790.8040.2400.1070.672
931.3851.1400.3130.2380.675
11155.900−0.0870.0190.0010.894
1250.4151.2180.7290.1950.456
13140.6851.0317.1040.3720.021
14163.3181.2814.1000.2270.062
15100.2961.6254.7110.3710.062
16211.9371.2273.9590.1800.062
1783.698−0.5000.5210.080.497
18531.9580.84214.9880.2270.000
1990.3121.1810.6810.0890.436
20552.6500.94912.8680.1950.001
216456.0820.90534.7560.3590.000
226322.1891.20360.4370.5020.000
Note: the zero values for the 66 erosive events were assigned to 0.01 for power function regression. The relationships between the H and SLR were not given for plots 8 and 10 because no runoff occurred on these two plots.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Fang, H. Water-Saving Soil Conservation Measures Should Be Used in Northern China: Evidence from Runoff Plot Data. Water 2021, 13, 853. https://doi.org/10.3390/w13060853

AMA Style

Fang H. Water-Saving Soil Conservation Measures Should Be Used in Northern China: Evidence from Runoff Plot Data. Water. 2021; 13(6):853. https://doi.org/10.3390/w13060853

Chicago/Turabian Style

Fang, Haiyan. 2021. "Water-Saving Soil Conservation Measures Should Be Used in Northern China: Evidence from Runoff Plot Data" Water 13, no. 6: 853. https://doi.org/10.3390/w13060853

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop