Next Article in Journal
Modelling Water Absorption in Micronized Lentil Seeds with the Use of Peleg’s Equation
Next Article in Special Issue
Spatiotemporal Evolution of Landscape Ecological Risk Based on Geomorphological Regionalization during 1980–2017: A Case Study of Shaanxi Province, China
Previous Article in Journal
Modeling the Spatial Dimensions of Warehouse Rent Determinants: A Case Study of Seoul Metropolitan Area, South Korea
Previous Article in Special Issue
Environmental Risk Assessment of Accidental Pollution Incidents in Drinking Water Source Areas: A Case Study of the Hongfeng Lake Watershed, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Potential of Gully Erosion on the Yamal Peninsula, West Siberia

by
Aleksey Sidorchuk
Lomonosov Moscow State University, Faculty of Geography, Laboratory of Soil Erosion and Fluvial Processes, Leninskiye Gory 1, 119899 Moscow, Russia
Sustainability 2020, 12(1), 260; https://doi.org/10.3390/su12010260
Submission received: 13 November 2019 / Revised: 13 December 2019 / Accepted: 24 December 2019 / Published: 28 December 2019
(This article belongs to the Collection Risk Assessment and Management)

Abstract

:
The Yamal Peninsula occupies the northern part of the West Siberian Plain in Russia. This territory has rapidly developed due to the exploitation of several gas fields. At the same time, the Yamal Peninsula is one of the most severely gullied landscapes in the Arctic. The potential risk of damage to the environment or structures and the cost of such damages are very high there. The erosion potential is the cumulative erosion by runoff above critical, calculated for each point at a catchment. Calculations take into account the geomorphic, lithological, and vegetation cover thresholds, realized in the form of critical runoff depth of erosion initiation. It also takes into account action of all flows between the critical and maximum runoff. The calculations for several gullied catchments on the Yamal Peninsula show the uneven distribution of erosion potential level with the maximum of gully erosion on the steep banks of the river valleys and on gully heads with bare soil. The area with potential erosion in these catchments varies within the range of 17–33%. The erosion on the Yamal Peninsula is mainly of natural origin. It occurs on steep slopes and at the heads of gullies. These landforms are not used for exploitation camps and settlements. Nevertheless, the linear structures, such as railways, roads and pipelines, can cross these unstable landforms with the risk of damage. Erosion potential increases at the spots with bare soil, which appear due to both construction work and natural processes, such as slumping.

1. Introduction

Gully erosion has been and remains one of the most dangerous processes of land damage [1,2]. Although it mainly occurs on agricultural lands [3], in recent years more and more gullies are formed in urbanized and industrial areas [4,5]. A typical example is the territory of the Yamal Peninsula where lands previously used for pastures only by local residents [6] are actively used at present in the development of gas and oil fields. It is important to find methods for the development of these lands that would not lead to significant changes in the vegetation cover, and in unstable landforms.
The Yamal Peninsula occupies the northern part of the West Siberian Plain (Figure 1). The peninsula is surrounded by the waters of the Baydaratskaya Bay and Kara Sea from the west and north, and of the Ob Gulf in the east. The area of the peninsula is 122,000 km2, and it is 750 km long from north to south and 140–240 km wide. The territory of the peninsula is open to the air masses from Arctic Ocean and the continent. The summer is short and cold, the winter is windy and frosty. The average annual temperature in the south of Yamal is −6.6 °C, and in its northern part −10.2 °C [7]. The peninsula is located within the permafrost zone. The depth of the seasonally thawed layer varies from 0.2 m in the northern part of Yamal to 2.0 m in the south. The entire Yamal Peninsula is a lowland plain, composed of marine loam and alluvial sands. The predominant altitudes are 40–70 m above sea level, the highest point is 102 m.
The Yamal peninsula is an area of rapid technogenic development due to gas field exploitation. There are several fields already prepared for exploitation, and more will appear in future. At the same time, the Yamal Peninsula is one of the most severely gullied landscape in the Russian Arctic (Figure 2) with the density of gullies is up to 1–2 km/km2. Therefore, natural environment, buildings and structures potentially can be damaged by gully erosion and the cost of such damage can be very high.
The methodology of evaluation of gully erosion potential is related to the ideas of geomorphic thresholds and magnitudes and frequencies of erosion events. In the classic paper “Geomorphic thresholds: the concept and its applications”. Schumm [8] argued: “geomorphic thresholds can be of two types, and they can be redefined in the following way. A geomorphic threshold is a threshold of landform stability that is exceeded either by intrinsic change of the landform itself, or by a progressive change of an external variable”. Schumm [8] presented a gully system development as an example of both intrinsic and extrinsic threshold conditions. Gully head erosion and gully length growth occur in the condition of exceeding of an intrinsic geomorphic threshold. After drainage area and/or slope steepness above a gully head decreases below a threshold values, a gully head erosion and gully length growth stops. The value of this threshold can increase or decrease due to external variable changes: surface runoff or vegetation cover transformation. In the case of threshold decrease gully can become unstable again.
In the other classic paper, Wolman and Miller [9] wrote: “Almost any specific mechanism requires that a certain threshold value of force be exceeded. However, above this threshold or critical limit, there occurs a wide range in magnitude of forces which results from variations in intensity of precipitation, wind speed, etc.”. “Available evidence indicates that evaluation of the effectiveness of a specific mechanism and of the relative importance of different geomorphic processes in molding specific forms involves the frequency of occurrence as well as the magnitude of individual events”.
According to this statement, a measure of the potential risk that the natural environment or structure (if it exists) will be damaged by gullies can be taken proportional to the combination of the values and frequencies of all erosion events that exceed a certain geomorphic threshold above which this form of relief becomes unstable [3]. This combination is the amount of cumulative erosion during various events, or the quantity of sediment transported.
Patton and Schumm [10] proposed the formula to calculate threshold conditions for gully initiation that includes “intrinsic change of the landform itself” and “a change of an external variable”.
S c r   =   a A c r b
where S is the gully catchment slope gradient, A is the contributing area, and a and b the coefficient and exponent. The regional implications of Equation (1) showed broad variability of the coefficient a, and more stable value of the exponent b. Vandaele et al. [11] compiled the parameters of Equation (1) for 25 gullied basins and obtained the coefficient a in the range 0.02–0.5 and exponent b in the range 0.15–0.47. With the increasing of the number of investigated gullied basins, these ranges increase also ([12,13,14] and many others). It is obvious that the applicability of empirical Equation (1) for gully erosion potential prediction is limited due to uncertainty of the parameters. This paper describes a novel methodology of gully erosion potential evaluation in the conditions of high Arctic and presents several examples of calculation for the gullied catchments on the Yamal Peninsula.

2. Materials and Methods

2.1. Theoretical Meanings of Empirical Equations

There is a large number of theoretical justifications of Equation (1). There are several ways to derive Equation (1) from the expression for the critical tangential stress in the flow τ, at which erosion begins [15,16,17,18], or from the expression for the critical flow velocity Ucr [19,20]. The way used by Begin and Schumm [15,16] is following. If shear stress in flow is equal or more than its threshold value
τ τ c r ,
then the basin is gullied or potentially will be gullied in future. Tangential stress in a flow is the product of water density ρ, acceleration due to gravity g, flow depth d and flow slope S:
τ   =   ρ g d S
Flow depth d is a function of discharge Q
d   =   p Q m
and discharge Q is
Q   =   k M A
where M is runoff depth, uniform for a gully catchment, and k is the coefficient of transformation of units. Note that Begin and Schumm [15] used non-uniform M, which depended on contributing area A.
These formulas transform coefficients in Equation (1) into
a   =   τ c r ρ g p ( k M ) m ;   b   =   m .
Inequality (2) takes the form
ρ g p [ ( k M A ) m S ] c r τ c r .
Inequality (7) shows the critical spatial-temporal conditions in the basin to be potentially eroded at a given point with inclination Scr and contributing area Acr at a given time when surface runoff depth of some probability is equal to Mcr.
It is also possible to use the critical velocity of erosion initiation as a threshold condition.
U U c r .
Quantitatively, expression (8) is easy to derive from Equation (2). Nevertheless, critical velocity better describes erosion initiation, as the driving forces (lift and drag) depend on flow velocity, not on shear stress [19].
Flow velocity is represented by Chezy-Manning formula
U   =   S n d 2 / 3 .
Here, n is roughness coefficient. Then coefficients in Equation (1) can be calculated as
a   =   U c r 2 n 2 p 4 3 ( k M ) 4 3 m ;   b = 4 3 m
and inequality (8) takes the form
S c r n p 2 / 3 ( k M A ) c r 2 m / 3 U c r .
There are several other ways to derive Equation (1) from the formulas of hydraulics. For example, Zorina [20] transformed Equation (9) into
Q   =   U c r 4 n 3 S 3 2 W d .
Width/depth ratio is the function of discharge
W d   =   c Q f
Therefore,
a   =   U c r 8 / 3 n 2 c 2 3 ( k M ) 2 3 ( 1 f ) ;   b   =   2 3 ( 1 f ) .
The changes in exponent b is controlled by the method used to derive Equation (1), as well as by regional variability of the exponent m or f [21]. The factors of the changes in coefficient a are more complicated. This coefficient depends on the critical shear stress or velocity and on the surface runoff depth in addition to the factors of variability of the exponent b.

2.2. The Novel Approach to Erosion Potential Estimation

One of the main problems of calculations with Equation (7) or Equation (11) is an indefinite value of surface runoff depth for estimation of critical velocity or shear stress. The proposed methodology aims to avoid this uncertainty. According to the proposed novel approach, the inequality (11) (or any other formulation of Equation (2) or Equation (8)) can be transformed in such way that the indefinite critical surface runoff depth becomes the function of known arguments:
M M c r   =   ( n U c r ) 3 / 2 m S c r 3 / 4 m p 1 / m ( k A ) c r
Inequality (15) shows that a gully catchment will be eroded at a given point with inclination Scr and contributing area Acr when the surface runoff depth M is equal or more than the critical value. The ground and vegetation cover conditions on this catchment correspond to critical velocity of erosion initiation Ucr.
The potential erosion R (depth of erosion at a given point) is the cumulative erosion during subcritical runoff.
R   =   T o T E M d t .
Here EM is erosion rate under all given runoff depths M during time interval dt, T-T0 is duration of surface runoff depth on the catchment above critical M c r for a given point. According to [22],
E M   =   k e q S .
Here, ke is erodibility coefficient, q-specific discharge Q/W. The channel width can be represented as W = pwQmw, then
E M   =   k e ( A M ) 1 m w p w S .
Gully erosion potential is calculated with the Equation (16) for those lines on the catchment, where flow velocity is more than critical, and erosion by water can be initiated. Outside these lines, erosion by water cannot be initiated at a given morphological, hydrological, lithological and vegetation cover conditions. All the processes of gully formation related to water flow can occur only within the area with initially critical trigger factors, although these factors (flow velocity and slope) can change through time. It is in this area that there is a risk of linear erosion by a concentrated flow.
The erosion by water in the Arctic conditions develops in two ways: mechanical erosion described by Equation (17) and thermal erosion, when frozen soil is thawed by more warm water, and the thawed layer is removed by flow. Investigation showed [23] that thermal erosion rate by a given flow is always less than values calculated with Equation (17).
The second important process, which is not described with Equation (17), is slumping of thawed soil from the gully walls to the bottom. This process in general decreases the rate of gully incision as the gully bottom occasionally rise due to accumulation of slumped sediments. Therefore, Equation (17) describes the upper limit of the combination of erosion rate and erosion factors.
Slumping from gully walls increases the gully top width and, to some extent, the gully length. This process is slower than erosion by water. It becomes significant at the late stages of gully development. With the general stabilization of a gully, its slopes reach the angle of repose ϕ. Then the gravitational processes stop. Nevertheless, if we take into account gravitational processes, this increases the estimate of the area with gully erosion potential, calculated with Equation (16). This increase should be taken into account at the late stages of the development of the gully. For example, a gully can become longer by Δ L   =   d H cot ( ϕ ) , where dH is gully depth at the head. For the stable gullies on the Yamal Peninsula the typical angle of repose ϕ is 0.6 radians, gully depth at the head is in the range 2–4 m [23]. Therefore, an additional length of stable gullies compares to the calculated with Equation (16) is around 3–6 m.
Potential erosion in Equation (16) for a given catchment is controlled by seven local parameters (Mmax, Ucr, n, m, p, mw, pw, ke), three variables (Scr, kAcr, Mcr) distributed in space, and one variable (M) changing through time t between T0 and T. There are different ways to estimate these variables and parameters for the Yamal Peninsula. Two of these variables (slope S and contributing area A) were taken from a digital elevation model, six parameters (Ucr, n, m, p, mw, pw, ke) were derived from field experiments, and other parameters and variables (Mmax, Mcr and M) were calculated with appropriate models.
The digital elevation models of the catchments were extracted from ArcticDEM [24] with the horizontal resolution 2 m. Modules of geographic information system QGIS with GIS SAGA [25] were used sequentially to fill sinks (SAGA Terrain Analysis-Hydrology-Fill sinks), to calculate flow accumulation (SAGA Terrain Analysis-Hydrology-Catchment area) and to calculate slope (SAGA Terrain Analysis-Morphology-Slope, aspect, curvature).
Field measurements in 1990–97 and in 2007 showed [23] that for the gullies on the Yamal Peninsula the exponent m is 0.3 and coefficient p is 0.21, the exponent mw is 0.45 and coefficient pw is 3.4 (Figure 3). A typical roughness coefficient n is 0.06 (see Table 4 in ref. [23]). Erodibility coefficient vary from 5.3 to 18.6, depending on soil type (see Figures 2 and 3 in ref. [23]).
The conditions of the erosion initiation for cohesive soils, which are typical for the gully erosion, are well investigated for a broad variety of soil types, both in laboratory and in the field [19]. It is possible to select an appropriate value for a given environment. According to field investigations on the Yamal Peninsula, the erosion under natural vegetation cover is negligible. Only bare ground is eroded by mechanical and thermal action of water. Measurements in [26] showed that erosion initiation in the loam deposits typical of the most part of Yamal Peninsula begins at the velocity 0.2 m/s. Therefore, critical daily surface runoff in mm is:
M c r   =   3.17 S c r 2.5 A
and the erosion rate in kg/(m2 s) is
E M   =   5.5 ( A M ) 0.55 S .

2.3. Hydrological Modelling

The systematic hydrological measurements are absent for the Yamal Peninsula, only some information from a few short expeditions is available [27,28]. Therefore, hydrological modelling is the main source of information on the changes in surface runoff depth through time and its frequency distribution in this region. Meteorological data for hydrological modelling are available from a few stations on the peninsula or from the reanalysis grid. In this paper, Era Interim reanalysis [29] with the grid step of 0.75° was used for the hydrological calculations.
In the conditions of Arctic climate in the region with deep permafrost, hydrological modelling bears specific characteristics. The main simplification of such modelling is the possibility to neglect water infiltration into frozen or highly saturated soil. Two main sources of water flow are typical for this environment: snow thaw during the spring and rainfall, mostly during the summer. The surface runoff depth Ms for the period of snow thaw was calculated with the models from [30,31,32,33] using the depth of snow pack (water equivalent), precipitation, air temperature and evaporation. The algorithm of calculating the runoff during the summer rains is based on precipitation depth and evaporation from Era Interim reanalysis with additional losses due to filling the depressions on the catchment surface calculated after Popov [34].
Hydrological modelling shows that the duration TD (days) of daily surface runoff depth ≥ M more than 1 mm is well approximated by an exponent:
T D   =   A N exp ( B N M ) .
Therefore, the module of time interval dt of daily surface runoff depth M is
d t   =   A N B N exp ( B N M ) d M .
Using this approximation for surface runoff frequency, the formula for calculating the gully erosion potential level (GEPL) takes the following form:
R L   =   ln { 5.5 A N B N A 0.55 S M c r M m a x M 0.55 exp ( B N M ) d M }   = ln { 5.5 A N B N A 0.55 S [ γ ( 1.55 ; B N M m a x ) γ ( 1.55 ; B N M c r ) ] } .
where γ is the lower incomplete gamma function. The logarithms in Equation (23) are to decrease the range of the calculated values.
The integral in Equation (23) is the duration of runoff above Mcr, calculated with the weight function M0.55. The logarithm of this integral (R0) decreases with increase of Mcr and BN (Figure 4). The maximum value of R0 (when Mcr = 0) is
R 0 max   =   ln [ γ ( 1.55 ; B N M m a x ) ] .
The minimum value of R0 (when Mcr = Mmax) is minus infinity. It is necessary to assign some maximum Mcr slightly lower than Mmax. If Mcr is 99% of Mmax then
R 0 min   =   ln [ γ ( 1.55 ; B N M m a x ) γ ( 1.55 ; 0.99 B N M m a x ) ] .
Value -R0min is added to Equation (23) to avoid negative erosion potential levels.
It is obvious that at the points with critical runoff depth Mcr exceeding the annual maximum of the surface runoff Mmax the erosion rate is zero. Therefore the mean gully erosion potential (MGEP) R for a given territory or a catchment can be calculated as the ratio of number of points (pixels) with Mcr < Mmax to the total number of points (pixels).

2.4. Model Validation

It is not easy to validate the results of erosion potential calculations, as this potential almost never occurs in the natural conditions, except in the badlands. Usually erosion potential is calculated to prevent dangerous results of erosion by soil conservation measures and in practice gully erosion stop at some stage of its development. The second difficulty in validation procedures is the inaccuracy of existing DEMs used for erosion potential calculations. Linear erosion can follow along such lines on the catchment as vehicle tracks, locally damaged vegetation, the local sources of water, etc., which cannot be implemented into DEM of any resolution.
The common definition of validation offered by Society for Computer Simulation Technical Committee on Model Credibility [35] is “Substantiation that a computerized model within its domain of applicability possesses a satisfactory range of accuracy consistent with the intended application of the model”. In our case, “a satisfactory range of accuracy” depends on the difficulties mentioned above. The territory of Yamal Peninsula is not badlands, although in several parts with erosion net density about 2 km/km2 it is close to such condition. ArcticDEM, the base for our calculations, depends on the resolution of the images used for DEM construction. Therefore, the model validation can be only approximate. The catchments selected for the validation procedure (point 2 in Figure 1 and Figure 5A) are situated close to the severely gullied part of the peninsula. The human activity at the KEKH exploitation camp (Bovanenkovskoye gas field) in 1986–1997 led to a nearly complete destruction of the vegetation cover there. Images from 1986–1995 and up to 2016 show that several bank gullies [23], which already became stable in virgin conditions, resumed their growth in length, mostly in the form of deep rills and ephemeral gullies (Figure 5A). At these catchments, the distribution of erosion potential level (Figure 5B) calculated with Equation (23) and Ucr = 0.2 m/s is quite similar to the erosion net, which already existed in 1995 (Figure 5C). Although there is no complete spatial overlapping of these two nets, the qualitative picture shows a consistence between the model results and measurements and “a satisfactory range of accuracy”. The mean erosion potential for this territory is 70%, which corresponds to the degree of the damage to the vegetation cover.

3. Results

The distribution of erosion potential level was calculated for seven small gullied catchments on the Yamal Peninsula (see Figure 1, Table 1). Hydrological characteristics were calculated for the nearest nodes of Era Interim net for the period 1979–2018. The maximum of the daily surface runoff Mmax vary at those points in the range of 74–118 mm. For investigated catchments the value of -R0min = 4 was used in calculations.
The distribution of the GEPL calculated with Equation (23) for the gully basin at the headwaters of a small river Nekhedeyvar-Yakha in the north of the Yamal Peninsula (72.85° N, 70.87° E) is highly uneven (Figure 6 and Figure 7).
A combination of critical surface runoff, local slope and contribution area for a given critical velocity of erosion initiation develops a complex pattern of erosion potential level. The maximum GEPL is associated with steep banks of river valleys and with gully heads. These landforms are vulnerable to shallow landslides which damage the vegetation cover [36]. In this case, we used in the calculations critical velocity for erosion initiation 0.2 m/s typical for bare ground conditions (Figure 6). The risk of gully erosion on the same catchment in the conditions of well-preserved tundra vegetation (Ucr = 0.5 m/s) is significantly lower (Figure 7). Still, the heads of gullies bear a high potential of erosion.
The calculations for seven small catchments situated in different parts of the Yamal Peninsula (Table 1) show similar, highly uneven distributions of GEPL with the same maximums on the steep banks of the river valleys and in gully heads.

4. Discussion

The gully erosion potential level RL for a given catchment is the sum of two logarithms: logarithm of product of slope S and contributing area A powered on 0.55 and logarithm of weighted duration of surface runoff depth above Mcr. The first component depends on paired distribution of slope and contributing area and is invariant for a given catchment. The second component depends on Mcr, which is also the function of slope and contributing area, as well as on parameters Mmax, Ucr, n, m, p, mw, pw, ke. There are several options of estimation Mcr. We shall discuss only three of them that use Chezy-Manning formula (9) and critical flow velocity Ucr.
In option 1, the regime equation used in Equation (9) is d   =   p Q m .
Then,
M c r   =   ( n U c r ) 3 / 2 m S c r 3 / 4 m p 1 / m ( k A ) c r
Using n = 0.06, p = 0.21, m = 0.3, k = 1/86400000 we obtain
M c r   =   n 5 U c r 5 S c r 2.5 p 3.33 ( k A ) c r   =   12140 U c r 5 S c r 2.5 A
If in Equation (9) (option 2) d   =   Q W d   then
M   =   1 k A [ U c r 5 / 3 n p w 2 / 3 S ] 3 2 ( 1 m w )
Using mw = 0.45 and pw = 3.4, we obtain
M   =   1 k A [ U c r 5 / 3 n p w 2 / 3 S ] 2.73   =   369963 U c r 4.55 S 1.37 A
According to option 3 in the approach of Zorina [20]
M   =   1 k A [ U c r 4 n 3 c S 3 / 2 ] 1 1 f
Using f = 0.15 and c = 16.2 we obtain
M   =   1 k A [ U c r 4 n 3 c S 3 / 2 ] 1.18   =   109244 U c r 4.72 S 1.77 A
Calculations for the gully basin at the headwaters of a small river Nekhedeyvar-Yakha in the north of the Yamal Peninsula (72.85° N, 70.87° E) with these three options show that the maximum differentiation of erosion potential is achieved if using option 1 (Equation (27) for Mcr) (Figure 8). Therefore, this option was used in all calculations described in the paper.
All three options are very sensitive to variations in critical velocity, so that possible error of Mcr calculation is 4.55–5-times larger than the error of Ucr estimation. The error affects both the mean gully erosion potential and distribution of its level. Sensitivity of calculations to the value of critical velocity must be taken into account in the situations when the quality of vegetation cover becomes one of the main factors of gully erosion. This is important when the restoration of territories with industrial damage is achieved by the improvement of vegetation cover. In such cases, the calculations of gully erosion potential should be accompanied with careful estimations of critical velocity of erosion initiation for a given quality of vegetation cover [37,38].

5. Conclusions

The novel measure of gully erosion potential can be taken proportional to the combination of magnitudes and frequencies of all erosion events exceeding some threshold, above which a given landform becomes unstable [9]. This combination was implemented into Equation (16) and represents the erosion depth. This equation takes into account the geomorphic, lithological and vegetation cover thresholds, realized in the form of critical runoff depth Mcr of erosion initiation (Equation (15) or (27)). It also takes into account action of all flows between the values of this threshold and maximum runoff depth Mmax. The logarithm of potential erosion depth was assigned as gully erosion potential level (GEPL) at a given point (Equation (23)), and the ratio of number of points (pixels) with Mcr < Mmax to the total number of points (pixels) was assigned as mean potential of gully erosion (MPGE) for entire territory or catchment. Model validation shows the consistence of calculations with the measurements and “a satisfactory range of accuracy”.
The calculations for seven gullied catchments on the Yamal Peninsula (Table 1) shows the highly uneven distribution of erosion potential level with the maximums on steep banks of the river valleys and in gully heads. The same landforms are vulnerable to shallow landslides, which damage the vegetation cover. Therefore, critical velocity for erosion initiation 0.2 m/s, which is suitable for bare ground, was used in calculations. The histograms of erosion potential level above zero are unimodal with probability maximums at 9.7–10.4. The mean erosion potential for these territories vary in a rather narrow range 0.17–0.33. This erosion potential is mostly natural, linked to the landforms with steep slopes and unstable vegetation cover. Such landforms are not used for exploitation camps and settlements. Nevertheless, the linear structures, such as railways, roads and pipelines, can cross these unstable landforms with high level of risk to be damaged. This erosion potential increases due to appearance of the spots of bare soil during construction works. This was shown for the territory of KEKH exploitation camp, where MPGE is 0.7. Such spots of bare soil can appear even on flat territories with an initially low level of erosion potential. Further work should be focused on calculations of erosion potential level for the territories of future development of gas exploitation to avoid possible damage to the environment and infrastructure.

Funding

This research was funded by RFBR grant 18-05-60147 “Extreme hydrometeorological phenomena in the Kara Sea and the Arctic coast”.

Acknowledgments

Era Interim data were collected by T.A. Matveeva.

Conflicts of Interest

The author declare no conflict of interest.

References

  1. Bennett, S.; Wells, R. Gully erosion processes, disciplinary fragmentation, and technological innovation. Earth Surf. Process. Landf. 2019, 44, 46–53. [Google Scholar] [CrossRef] [Green Version]
  2. Valentin, C.; Poesen, J.; Li, Y. Gully erosion: Impacts, factors and control. CATENA 2005, 63, 132–153. [Google Scholar] [CrossRef]
  3. Liu, X.; Li, H.; Zhang, S.; Cruse, R.M.; Zhang, X. Gully Erosion Control Practices in Northeast China: A Review. Sustainability 2019, 11, 5065. [Google Scholar] [CrossRef] [Green Version]
  4. Carvalho, R.F.; Lopes, P.; Leandro, J.; David, L.M. Numerical Research of Flows into Gullies with Different Outlet Locations. Water 2019, 11, 794. [Google Scholar] [CrossRef] [Green Version]
  5. Gong, C.; Lei, S.; Bian, Z.; Liu, Y.; Zhang, Z.; Cheng, W. Analysis of the Development of an Erosion Gully in an Open-Pit Coal Mine Dump During a Winter Freeze-Thaw Cycle by Using Low-Cost UAVs. Remote Sens. 2019, 11, 1356. [Google Scholar] [CrossRef] [Green Version]
  6. Forbes, B.; Stammler, F.; Kumpula, T.; Mysshtyb, M.; Pajunen, A.; Kaarlejärvi, E. High resilience in the Yamal-Nenets social-ecological system, West Siberian Arctic, Russia. Proc. Natl. Acad. Sci. USA 2009, 106, 22041–22048. [Google Scholar] [CrossRef] [Green Version]
  7. Nature of Yamal (Priroda Yamala); Dobrinski, L.N. (Ed.) Nauka: Ekaterinburg, Russia, 1995; p. 435. (In Russian) [Google Scholar]
  8. Schumm, S.A. Geomorphic thresholds: The concept and its applications. Trans. Inst. Br. Geogr. 1979, 4, 485–515. [Google Scholar] [CrossRef]
  9. Wolman, M.G.; Miller, J.P. Magnitude and frequency of forces in geomorphic processes. J. Geol. 1960, 68, 54–74. [Google Scholar] [CrossRef] [Green Version]
  10. Patton, P.C.; Schumm, S.A. Gully erosion, northwestern Colorado: A threshold phenomenon. Geology 1975, 3, 83–90. [Google Scholar] [CrossRef]
  11. Vandaele, K.; Poesen, J.; Govers, G.; van Wesemael, B. Geomorphic threshold conditions for ephemeral gully incision. Geomorphology 1996, 16, 161–173. [Google Scholar] [CrossRef]
  12. Desmet, P.J.; Poesen, J.; Govers, G.; Vandaele, K. Importance of slope gradient and contributing area for optimal prediction of the initiation and trajectory of ephemeral gullies. Catena 1999, 37, 377–392. [Google Scholar] [CrossRef]
  13. Vandekerckhove, L.; Poesen, J.; Oostwoudwijdenes, D.; Nachtergaele, J.; Kosmas, C.; Roxo, M.J.; de Figueiredo, T. Thresholds for gully initiation and sedimentation in Mediterranean Europe. Earth Surf. Process. Landf. 2000, 25, 1201–1220. [Google Scholar] [CrossRef]
  14. Garrett, K.K.; Wohl, E.E. Climate-invariant area–slope relations in channel heads initiated by surface runoff. Earth Surf. Process. Landf. 2017, 42, 1745–1751. [Google Scholar] [CrossRef]
  15. Begin, Z.B.; Schumm, S.A. Instability of alluvial valley floor: A method for its assessment. Trans. Am. Soc. Agric. Eng. 1979, 22, 347–350. [Google Scholar] [CrossRef]
  16. Begin, Z.B.; Schumm, S.A. Gradational thresholds and landscape singularity: Significance for Quaternary studies. Quat. Res. 1984, 21, 267–274. [Google Scholar] [CrossRef]
  17. Harvey, M.D.; Watson, C.C.; Schumm, S.A. Gully Erosion; Tech. Note 366; U.S. Dep. of the Inter., Bur. of Land Management: Washington, DC, USA, 1985; p. 181.
  18. Montgomery, D.R.; Dietrich, W.E. Landscape dissection and drainage area-slope thresholds. In Process Models and Theoretical Geomorphology; Kirby, M.J., Ed.; Wiley: New York, NY, USA, 1994; pp. 221–246. [Google Scholar]
  19. Mirtskhulava, T.Y. Osnovy Fiziki i Mekhaniki Erozii Rusel (Principles of Physics and Mechanics of Channel Erosion); Gidrometeoizdat: Leningrad, Russia, 1988; p. 303. (In Russian) [Google Scholar]
  20. Zorina, Y.F. Quantitative methods of estimating of gully erosion potential. In Eroziya Pochv i Ruslovyye Protsessy (Soil Erosion and Channel Processes); Chalov, R.S., Ed.; Moscow Univ. Press: Moscow, Russia, 1979; Volume 7, pp. 81–90. (In Russian) [Google Scholar]
  21. Park, C.C. World-wide variations in hydraulic geometry exponents of stream channels: An analysis and some observations. J. Hydrol. 1977, 33, 133–146. [Google Scholar] [CrossRef]
  22. Sidorchuk, A. Dynamic and static models of gully erosion. Catena 1999, 37, 401–414. [Google Scholar] [CrossRef]
  23. Sidorchuk, A. Gully erosion in the cold environment: Risks and hazards. Adv. Environ. Res. 2015, 44, 139–192. [Google Scholar]
  24. ArcticDEM, 2018. Harvard Data V1. Available online: https://doi.org/10.7910/DVN/OHHUKH (accessed on 28 February 2019).
  25. Quantum GIS Development Team, 2004–2019. V3.10. Available online: https://qgis.org (accessed on 28 October 2019).
  26. Sidorchuk, A.; Grigor’ev, V. Soil erosion on the Yamal peninsula (Russian Arctic) due to gas field exploitation. Adv. GeoEcol. 1998, 31, 805–811. [Google Scholar]
  27. Bobrovitskaya, N.N.; Baranov, A.V.; Vasilenko, N.N.; Zubkova, K.M. Hydrological conditions. In Erosion Processes at the Central Yamal (Erozionnyye Protsessy Tsentral’nogo Yamala); Sidorchuk, A., Baranov, A., Eds.; RNII KPN: St. Petersburg, Russia, 1999; pp. 90–105. (In Russian) [Google Scholar]
  28. Hydrology of the Wetlands of the Permafrost Zone of Western Siberia (Gidrologiya Zabolochennykh Territoriy Zony Mnogoletney Merzloty Zapadnoy Sibiri); Novikov, S.M. (Ed.) VVM Publ. House: St. Petersburg, Russia, 2009; p. 536. (In Russian) [Google Scholar]
  29. Berrisford, P.; Dee, D.P.; Poli, P.; Brugge, R.; Fielding, M.; Fuentes, M.; Kållberg, P.W.; Kobayashi, S.; Uppala, S.; Simmons, A. The ERA-lnterim archive, Version 2.0. ERA Rep. Ser. 2011, 1, 23. [Google Scholar]
  30. Komarov, V.D.; Makarova, T.T.; Sinegub, E.S. Calculation of the hydrograph of floods of small lowland rivers based on thaw intensity data. Proc. USSR Hydrometeorol. Cent. 1969, 37, 3–30. (In Russian) [Google Scholar]
  31. Vinogradov, Y.B. Mathematical Modeling of Flow Formation Processes (Matematicheskoye Modelirovaniye Protsessov Formirovaniya Stoka); Gidrometeoizdat: Leningrad, Russia, 1988; p. 312. (In Russian) [Google Scholar]
  32. Vinogradov, Y.B.; Vinogradova, T.A.; Zhuravlev, S.A.; Zhuravleva, A.D. Mathematical modeling of hydrographs from the unstudied river basins on the Yamal peninsula. Bull. St. Petersburg State Univ. 2014, 7, 71–81. (In Russian) [Google Scholar]
  33. Gelfan, A.N. Model of formation of river flow during snowmelt and rain. In Erosion Processes at the Central Yamal (Erozionnyye Protsessy Tsentral’nogo Yamala); Sidorchuk, A., Baranov, A., Eds.; RNII KPN: St. Petersburg, Russia, 1999; pp. 205–225. (In Russian) [Google Scholar]
  34. Popov, Y.G. Analysis of Lowland River Flow Formation (Analiz Formirovaniya Stoka Ravninnykh Rek); Gidrometeoizdat: Leningrad, Russia, 1956; p. 131. (In Russian) [Google Scholar]
  35. Schlesinger, S.; Crosbie, R.E.; Gagne, R.E.; Innis, G.S.; Lalwani, C.S.; Loch, J.; Sylvester, R.J.; Wright, R.D.; Kheir, N.; Bartos, D. Terminology for model credibility. Simulation 1979, 32, 103–104. [Google Scholar]
  36. Voskresenski, K.S.; Golovenko, S.S. Distribution of shallow landslides at the central Yamal. In Erosion Processes at the Central Yamal (Erozionnyye Protsessy Tsentral’nogo Yamala); Sidorchuk, A., Baranov, A., Eds.; RNII KPN: St. Petersburg, Russia, 1999; pp. 133–138. (In Russian) [Google Scholar]
  37. Vasilevskaya, V.D.; Grigoriev, V.Y.; Sidorchuk, A.Y. Stability of the Northern Soils and Ecosystems to Technogenical Influence. In The Development of the North and Problems of Recultivation, Proceedings of the III International Conference, Sankt-Petersburg, Russia, 27–31 May 1996; Komi NTsUrO RAN: Syktyvkar, Russia, 1997; pp. 208–212. [Google Scholar]
  38. Gyssels, G.; Poesen, J.; Bochet, E.; Li, Y. Impact of plant roots on the resistance of soils to erosion by water: A review. Progr. Phys. Geogr. 2005, 29, 189–217. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The Yamal Peninsula and surrounding areas on the Google Earth image. Stars with the numbers show the position of catchments with calculated gully erosion potential. The upper map shows the position of Yamal Peninsula in the Northern Eurasia.
Figure 1. The Yamal Peninsula and surrounding areas on the Google Earth image. Stars with the numbers show the position of catchments with calculated gully erosion potential. The upper map shows the position of Yamal Peninsula in the Northern Eurasia.
Sustainability 12 00260 g001
Figure 2. The valley of a small river on the central Yamal (71.14° N, 69.96° E) with bank gullies (photo taken by I. Krylenko in 2007) and the map of erosion net density (the fragment of unpublished map of erosion risk in the territory of the USSR, Institute of Geography, Russian Academy of Science, kindly provided by its author, A.N. Makkaveev).
Figure 2. The valley of a small river on the central Yamal (71.14° N, 69.96° E) with bank gullies (photo taken by I. Krylenko in 2007) and the map of erosion net density (the fragment of unpublished map of erosion risk in the territory of the USSR, Institute of Geography, Russian Academy of Science, kindly provided by its author, A.N. Makkaveev).
Sustainability 12 00260 g002
Figure 3. Regime equations for the gullies of the west-central Yamal, estimated by field measurements in 1991–1997 [23].
Figure 3. Regime equations for the gullies of the west-central Yamal, estimated by field measurements in 1991–1997 [23].
Sustainability 12 00260 g003
Figure 4. The weighted duration of runoff R0 via Mcr.
Figure 4. The weighted duration of runoff R0 via Mcr.
Sustainability 12 00260 g004
Figure 5. (A). The erosion net of gullies and rills in 1995, projected on the Google Earth image of 2016 within the territory of KEKH exploitation camp, Bovanenkovskoye gas field (70.37° N, 68.42° E). (B). The distribution of erosion potential level calculated with critical velocity Ucr = 0.2 m/s. (C). Comparison of calculated erosion potential level (9.5–13.5) with the erosion net in 1995.
Figure 5. (A). The erosion net of gullies and rills in 1995, projected on the Google Earth image of 2016 within the territory of KEKH exploitation camp, Bovanenkovskoye gas field (70.37° N, 68.42° E). (B). The distribution of erosion potential level calculated with critical velocity Ucr = 0.2 m/s. (C). Comparison of calculated erosion potential level (9.5–13.5) with the erosion net in 1995.
Sustainability 12 00260 g005
Figure 6. The gully erosion potential level distribution at the headwaters of Nekhedeyvar-Yakha River (72.85° N, 70.87° E) calculated with the critical velocity Ucr = 0.2 m/s. The black contours on the map show elevations in ArcticDEM system.
Figure 6. The gully erosion potential level distribution at the headwaters of Nekhedeyvar-Yakha River (72.85° N, 70.87° E) calculated with the critical velocity Ucr = 0.2 m/s. The black contours on the map show elevations in ArcticDEM system.
Sustainability 12 00260 g006
Figure 7. The distribution of gully erosion potential level at the headwaters of Nekhedeyvar-Yakha River (72.85° N, 70.87° E) calculated with critical velocity Ucr = 0.5 m/s. The black contours on the map show elevations in ArcticDEM system.
Figure 7. The distribution of gully erosion potential level at the headwaters of Nekhedeyvar-Yakha River (72.85° N, 70.87° E) calculated with critical velocity Ucr = 0.5 m/s. The black contours on the map show elevations in ArcticDEM system.
Sustainability 12 00260 g007
Figure 8. The histograms of erosion potential level RL in the headwaters of the Nekhedeyvar-Yakha River calculated with different options of Mcr estimations (see text for explanations).
Figure 8. The histograms of erosion potential level RL in the headwaters of the Nekhedeyvar-Yakha River calculated with different options of Mcr estimations (see text for explanations).
Sustainability 12 00260 g008
Table 1. The main characteristics of gully erosion potential level calculated for the gullied catchments on the Yamal Peninsula (for positions see Figure 1).
Table 1. The main characteristics of gully erosion potential level calculated for the gullied catchments on the Yamal Peninsula (for positions see Figure 1).
Catchment1234567
Amax, km21.541.141.777.321.341.5729.55
Mean erosion potential0.260.220.330.310.280.170.25
Erosion potential level of maximum probability9.710.0410.3510.1610.210.1610.03
Mmax, mm7497.697.6107.611411886.2
AN2120.620.620.819.016.919.5
BN0.0840.0740.0740.0730.0510.0580.075
1. Nekhedeyvar-Yakha river in the north of the Yamal Peninsula (72.85° N, 70.87° E); 2. Dry valley basin with PBB exploitation camp, Bovanenkovo gas field (70.39° N, 68.35° E); 3. Khangalava-Yakha river in the west of the Yamal peninsula (70.35° N, E 68.72° E); 4. Catchments with a section of the Yamal railway at 70.32° N, 68.79° E; 5. Catchments with exploitation camps of the Novoportovskoye gas field (67.81° N, 72.53° E); 6. Gully in the east of the Yamal Peninsula (68.10° N, 73.03° E); 7. Headwaters of small river Nerovarta-Yakha in the basin of Kharasavey-Yakha river (70.82° N, 68.23° E).

Share and Cite

MDPI and ACS Style

Sidorchuk, A. The Potential of Gully Erosion on the Yamal Peninsula, West Siberia. Sustainability 2020, 12, 260. https://doi.org/10.3390/su12010260

AMA Style

Sidorchuk A. The Potential of Gully Erosion on the Yamal Peninsula, West Siberia. Sustainability. 2020; 12(1):260. https://doi.org/10.3390/su12010260

Chicago/Turabian Style

Sidorchuk, Aleksey. 2020. "The Potential of Gully Erosion on the Yamal Peninsula, West Siberia" Sustainability 12, no. 1: 260. https://doi.org/10.3390/su12010260

APA Style

Sidorchuk, A. (2020). The Potential of Gully Erosion on the Yamal Peninsula, West Siberia. Sustainability, 12(1), 260. https://doi.org/10.3390/su12010260

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