Next Article in Journal
Temporal Dynamics of Fish Assemblages as a Reflection of Policy Shift from Fishing Concession to Co-Management in One of the World’s Largest Tropical Flood Pulse Fisheries
Previous Article in Journal
Controls on Reservoir Heterogeneity of a Shallow-Marine Reservoir in Sawan Gas Field, SE Pakistan: Implications for Reservoir Quality Prediction Using Acoustic Impedance Inversion
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Impact of Land Use Changes on Soil Erosion in the River Basin of Miocki Potok, Montenegro

1
Faculty of Philosophy, Geography, University of Montenegro, D. Bojovica, 81400 Niksic, Montenegro
2
Biotechnical Faculty, University of Montenegro, 81000 Podgorica, Montenegro
3
Planning and Budgeting Organization of North Khorasan, Bojnurd 94714, Iran
4
Faculty of Agriculture, University of Belgrade, 11070 Zemun, Serbia
5
Faculty of Management, Comenius University in Bratislava, 82005 Bratislava, Slovakia
6
Faculty of Economics and Engineering Management, University Business Academy, 21102 Novi Sad, Serbia
7
International Platform for Dryland Research and Education, Arid Land Research Centre, Tottori University, Tottori 680-0001, Japan
*
Authors to whom correspondence should be addressed.
Water 2020, 12(11), 2973; https://doi.org/10.3390/w12112973
Submission received: 15 September 2020 / Revised: 12 October 2020 / Accepted: 13 October 2020 / Published: 23 October 2020

Abstract

:
Land use change in all river basins leads to changes in hydrologic response, soil erosion, and sediment dynamics characteristics. Those changes are often viewed as the main cause of accelerated erosion rates. We studied the impact of land use changes on soil erosion processes in one of the watersheds in Montenegro: the Miocki Potok, using this watershed as a pilot river basin for this area. We simulated responses of soil erosion processes by using a process-oriented soil erosion Intensity of Erosion and Outflow (IntErO) model, with different settings of land use for the years 1970, 1980, 1990, 2000, 2010, and 2020. The model provides fast, effective, and affordable insight into the effects of land use change on soil erosion processes. Testing of the applied procedures was important for the further establishment of watershed management methodologies at the national level, for the other 300 river basins of Montenegro. For the current state of land use, calculated peak discharge for the Miocki Potok was 364 m3 s−1 (2020)–372 m3 s−1 (1970) for the incidence of 100 years, and there is a possibility for large flood waves to appear in the studied basin. Real soil losses, Gyear, were calculated at 13680 m3 year−1 (2020) and specific 333 m3 km−2 year−1 (2020). A Z coefficient value of 0.439 (2020) indicated that the river basin belongs to destruction category III. The strength of the erosion process was medium, and according to the erosion type, it was mixed erosion. According to our analysis, the land use changes in the last 50 years influenced a decrease in the soil erosion intensity for 14% in the Miocki Potok River Basin. Further studies should be focused on the detailed analysis of the land use changes trends with the other river basins at the national level, closely following responses of soil erosion to the changed land use structure, and effects of plant-and-soil interaction on soil erosion and sediment dynamics.

1. Introduction

Soil is one of the most important and most complex natural resources, but current developments (climate change, soil erosion, and urbanization) increasingly threaten this valuable resource. Soil degradation caused by erosion and rapid population increase is ranked among the most important environmental problems in the world [1,2,3]. Erosion is a key driver of land degradation, heavily affecting sustainable land management in various environments worldwide [4,5,6,7,8]. Human activities may have an important impact on the Earth’s surface processes and contribute to changing the Earth’s surface significantly since soil erosion and sediment dynamics are expected to change accordingly.
The negative effects of soil erosion on agricultural productivity and ecosystem are almost invisible over a shorter period of time [9,10]. Under normal conditions, soil erosion is a constant, progressing phenomenon, and, if no accelerated land degradation process such as gullying or landsliding occurs, we commonly are aware of the consequences only when it could be too late or too expensive to solve it. Conditions of increasing soil erosion not only affect agriculture but may dramatically reduce the life expectancy of reservoirs. In the long term, such processes may lead to changes in the provision of ecosystem services and to deplete fundamental natural resources.
The analysis based on gauging data can contribute to identifying the effects of land use changes to soil erosion and sediment dynamics, but these data are available only for a small number of medium to large catchments, whereas for small basins and a variety of landscapes, this kind of information is still missing. Thus, quantification of the responses to soil erosion and sediment yield to land use variations can only be achieved by numerical simulation of soil erosion processes [11].
In their 2009 study, Volk [12] encouraged researchers to conduct such analysis highlighting the importance of achieving the aims of the European Water Framework Directive. According to this directive, land use and land management options are to be used as tools for water quantity and quality control in order to achieve and maintain ecologically stable and productive water bodies. This idea on policy issues has been well received among researchers in Montenegro and Serbia, taking into consideration the current EU accession agenda of these two countries.
Numerous semi-quantitative methods are at the moment in use for assessing sediment yield at the catchment scale. These methods take into account parameters such as geology, soil, geomorphology, topography, climate, rainfall erosivity, and land cover [13,14,15].
The objective of this research was to simulate responses of soil erosion processes with the different settings of the land use in the River Basin of Miocki Potok in Montenegro over a 50 year time interval (1970–2020) by using a process-oriented soil erosion Intensity of Erosion and Outflow (IntErO) model [13]. The model provides fast, effective, and affordable insight into the effects of land use change on soil erosion processes and sediment yield. Simulation results under different scenarios of land use change were used to analyze the response characteristics of soil erosion processes to the variations in land use with the aim of optimizing current and future land use options and to predict or optimize future land and water management options in the wider region.
It is well known that soil erosion processes can be, to some extent, controlled by appropriate land management. This requires the collection of field data and the use of predictive models for the evaluation of different management scenarios for the protection of soils. Field measurements of erosion and sedimentation using classical techniques are time-consuming and expensive; therefore, modeling of the erosion process represents a valid alternative to adequately predict both the runoff and soil loss.
According to the various reports of the European Soil Bureau Institute for Environment and Sustainability and the European Environment Agency, the Universal Soil Loss Equation model (USLE) is extensively used in the following European countries: Austria, Bosnia, and Herzegovina (including the Republic of Srpska), Bulgaria, Italy, Hungary, Norway, Romania, Slovakia, Finland, Czech Republic, Spain, and Switzerland. The Revised Universal Soil Loss Equation (RUSLE) is used in Belgium; the UK, Germany, and France are using their domestic/national models.
In the Western Balkans, the Erosion Potential Method (EPM) for mapping the intensity of water erosion at the country scale is the preferred model [12,13,14,15]. The model is based on the long-term research of the soil erosion process in the Grdelica Gorge, led by professor Gavrilovic [16]. It is important to highlight that the EPM and its successor, IntErO, has also been tested in catchment areas all over the world: Serbia [14,15,16,17], Bosnia and Herzegovina [18,19], Croatia [9,20], North Macedonia [21], Italy [22], Slovenia [23], Chile [24], Brazil [25], Nepal [4], Morocco [26], Iraq [27], and in Iran [28,29,30,31,32,33,34,35]. The researchers concluded that the results are compatible with field observations.
Global change is severely impacting the biosphere. Such impacts are expected to increase unless mitigation management actions are implemented [35]. Despite the call from the scientific and political arenas, few studies have investigated and assessed the effectiveness of actions on freshwater-related services [36].
Assessing the sustainability of water supply is of key importance, especially in Montenegro, the wider Balkans, and the Mediterranean. Montenegro has been identified as a hot spot in regards to the impact of climate change and as a region in which water demand is continuously growing for population needs, tourism development, and irrigated areas expansion.
Climate projections in some scenarios predict an increase in temperature of up to 2 °C and a decrease in daily precipitation by 20%. In this context, there is a great need for experts in the Balkans to assess future changes in water distribution capacities in the basins, taking into account the newly emerging climate change, management of reservoirs and dams. It is necessary to work in a broader context on the preparation of scenarios for the use of available water.
This research attempts to establish one of the sustainable forms of modeling that would be calibrated and validated in the Polimlje region of the River Lim Basin and further used to assess how water resources can meet water requirements for local, agricultural, and eco sectors.

2. Materials

2.1. Study Area

Northeastern Montenegro is mainly mountainous, with the presence of deep valleys incised into limestone ranges. Other parts of the northeast are rather hilly and underlain by Palaeozoic rocks. In this region, the highest peaks of Montenegro are found, including Komovi (2487 m a.s.l.) and Zla Kolata (2535 m a.s.l.) in the Prokletije Mountains. Nevertheless, this area is quite densely populated and includes the towns of Plav, Andrijevica, Berane, and Bijelo Polje. The rivers in this region drain to the Black Sea, and some of them form deep canyons crossing limestone formations. Further downstream, they form broader valleys flowing through softer Palaeozoic material [13].
This study was conducted in the area of the river basin of Miocki Potok (Figure 1).
It is a left-hand tributary of the river Lim (External coordinates: 43°10′19.1″ N 19°43′41.6″ E-North; 43°06′26.3″ N 19°42′42.6″ E-South; 43°10′07.7″ N 19°41′14.2″ E-West; 43°08′12.2″ N 19°46′32.3″ E-East).
The river basin of Miocki Potok stretches from its confluence with the Lim (lowest elevation (Hmin), 535 m a.s.l) up to the tops of the Markov Kamen and Lisa mountains, where the highest elevation (Hmax) is 1553 and 1503 m a.s.l, respectively. The lower reaches (from 561 m to 535 m a.s.l) of the Miocki Potok are called the Kanjski Potok. The drainage network includes the following rivers: Miocki Potok–Kanjski Potok, with its right-hand tributary Dulev Potok and left-hand tributaries Dobrinjski Potok and Nikolji Potok (Figure 1).
This watershed encompasses an area of 41 km2, and the natural length of the main watercourse, Lv, is 6.4 km.
A flat area occurs around the village of Kanje in the lower reach, and steep slopes make up the upper part of the watershed. The average gradient of the catchment, Isr, is 31.6% and indicates that very steep slopes are present in the river basin. The average river basin altitude, Hsr, is 917.83 m a.s.l; the average elevation difference of the basin, D, is 382.83 m. The drainage density is low indicating a rather permeable substrate.
On the basis of more than 70 years (1948–2020) of temperature and precipitation data, the study area is characterized by a mountain Mediterranean climate with rainy autumns and springs, cold winters, and a deficit of precipitation in the summer months. Basic data on the area needed for the calculation of soil erosion intensity and runoff are presented in Table 1.
The absolute maximum air temperature ever recorded was 39.2 °C. Winters are severe, with negative temperatures as low as −27.6 °C. The average annual air temperature, t0, was 8.9 °C. The average annual precipitation, Hyear, was 873 mm. The temperature coefficient for the region, T, was calculated at 0.99. The torrential rain, hb, was calculated at 84.7 mm.

2.2. The Geological Structure and Soils of the Area

Montenegro is a part of the Dinaric Alps, which are, in turn, parts of the complex thrust-and-fold system of the Mediterranean area. This area represents a proxy of long-lasting interactions between Eurasia and Gondwana, resulting in a system of fold-and-thrust belts and associated foreland and back-arc basins. The system cannot be interpreted as the end product of one single Alpine orogeny as the major suture zones result from various tectonic events that closed different oceanic basins of the former Tethyan Ocean. The Dinarides-Albanides-Hellenides orogenic belt is caused by a Tertiary collision between the Adriatic promontory and the Serbo-Macedonian-Rhodope blocks. The belt is bordered to the west by a foreland basin in the Eastern Adriatic basin filled with Eocene-Quaternary deep marine sediments.
The study area at large consists of various types of sediment, magmatic and metamorphic rocks generated in the long, Palaeozoic to Quaternary, interval. Most of the terrain is underlain by Mesozoic formations of carbonate composition, while magmatic and silico-clastic rocks are substantially less present. Paleozoic geological formations consist of sedimentary and metamorphic, silico-clastic rocks found mostly in the north-eastern parts of Montenegro, while Cainozoic rocks of carbonate and clastic composition occur almost in all regions of Montenegro.
The main rocks outcropping in the area are clastic and subordinate carbonate rocks from the Paleozoic, Triassic clastites, volcanites, tuffs, limestone and dolomites, Jurassic clastic rocks with diabasic effusions and metamorphic rocks and Quaternary, mainly alluvial and colluvial deposits.
In order to define the permeability of the rocks of the study area, we used the Geological Atlas of Serbia [37] and extracted a geological map of the Miocki Potok watershed from the Geological map of Montenegro [38]. The basin of the Miocki Potok consists of Devonian-Carboniferous (D+C) and Permian (P) phyllites, argyllo-phyllites, metasandstones, and conglomerates.
The coefficient of the region’s permeability, S1, was calculated to be 0.91. Within the study basin, the area underlain by limestone is very permeable (class fp, 10%), and the rest has poor permeability (class fo, 90%).
According to the results of the field visits and supplementary laboratory analysis, but also using the previous research data of the project Soils of Montenegro (1964–1988) carried out by the team of the Biotechnical institute of the University of Montenegro [39] and Spalevic [13], the most common soil types in the study basin were: Dystric Cambisols and Fluvisols and Colluvial Fluvisols in the lower alluvial plain.

3. Methods

Morphometric methods were used to determine slope gradient, specific lengths, aspect, and shape. The depth of the erosion base and the density of erosion rills were identified during field visits. Google Earth was used to further investigate the morphological features of the watershed with a realistic 3D view. In that manner, some minor sites or phenomena that were not reached during the fieldwork were also taken into consideration. According to Chang et al. [40], Google Earth is considered more than adequate to identify and to characterize soil erosion processes.
Some pedological profiles were excavated, and soil samples were collected for physical and chemical analysis. The grain-size distribution of the soils was determined by the pipette method [41]; the soil samples were air-dried at 105 °C and dispersed using sodium pyrophosphate. The soil reaction (pH in H2O and nKCl) was determined with a potentiometer. Total carbonates were determined by the volumetric Scheibler method [42]; the content of the total organic matter was determined by the Kotzman method [43]; easily accessible phosphorous and potassium were determined by the Al-method [44], and the adsorptive complex (y1, S, T, V) was determined by the Kappen [45] method.
Blinkov and Kostadinov [46] evaluated the applicability of various erosion risk assessment methods for engineering purposes. The factors to be considered depends on the scale, various erosion tasks, and local specific needs. According to their findings, the Erosion Potential Method—EPM [16] was the most suitable at the catchment level for watershed management in southeastern European countries. The EPM is distinguished by its high degree of reliability in calculating sediment yields, suspended sediment yield, and reservoir sedimentation [47].
In Montenegro, the EPM proved to return good predictions of runoff and soil erosion intensity, specifically in the Region of Polimlje [13,48,49,50]. In this study, the Intensity of Erosion and Outflow (IntErO) program package [13] was used to obtain estimates of soil erosion intensity in the river basin of Miocki Potok. The method has clearly defined procedures, and subjective evaluations are reduced to a minimum. The EPM is embedded in the algorithm of the IntErO model.

3.1. The IntErO Model

The IntErO model uses the Erosion Potential Method in its algorithm background [13,26]. The IntErO, an upgrading of the River Basins [48,51] and the Surface and Distance Measuring [48] programs, is simple in handling and can be used to calculate a large number of data with the processing of 22 input parameters, returning, after the processing, 26 result parameters (Coefficient of the river basin form, A; Coefficient of the watershed development, m; Average river basin width, B; (A)symmetry of the river basin, a; Density of the river network of the basin, G; Coefficient of the river basin tortuousness, K; Average river basin altitude, Hsr; Average elevation difference of the river basin, D; Average river basin decline, Isr; The height of the local erosion base of the river basin, Hleb; Coefficient of the erosion energy of the river basin’s relief, Er; Coefficient of the region’s permeability, S1; Coefficient of the vegetation cover, S2; Analytical presentation of the water retention in inflow, W; Energetic potential of water flow during torrent rains, 2×gDF^½; Maximal outflow from the river basin, Qmax; Temperature coefficient of the region, T; Coefficient of the river basin erosion, Z; Production of erosion material in the river basin, Wyear; Coefficient of the deposit retention, Ru; Real soil losses, Gsp; Real soil losses per km2. The model considers six factors related to lithology (rocks permeability by percentage: fp, permeable; fpp, semipermeable; fo, low permeability) and soil type (erodibility coefficient, Y), topographic and relief data (I coefficient), monthly mean and annual precipitation (P coefficient), temperatures annual averages (t coefficient), land cover data (X coefficient), the state of erosion patterns, and development of the watercourse network ( Φ coefficient). The IntErO model can be characterized as semi-quantitative because it is based on a combination of descriptive and quantitative procedures. Compared to other semi-quantitative methods, this is the most quantitative because it uses descriptive evaluation for three parameters only: soil erodibility, soil protection, and the extent of erosion in the catchment. All other parameters are quantitative catchment descriptors. IntErO flowchart is presented in Figure 2.
The annual volume of soil detached due to soil erosion (Wyear) was calculated in m3 km−2 y−1 by the following Equation (1):
W y e a r = T · P · π · Z 3 · F s
  • P is the annual average rainfall in mm.
  • T is the temperature coefficient, calculated from the following equation:
T = t 10 + 0.1
t represents the mean annual temperature in °C, which is extracted from the meteorological data collected from the State Hydrometeorological Institute of Montenegro and the Biotechnical faculty of the University of Montenegro.
The Z coefficient describes the intensity of the erosive process, can be classified according to the degree of erosion (Table 2); it is calculated using Equation (3).
Z = X · Y · ( Φ + I )
X represents the coefficient of soil protection (Table 3), it is a dimensionless parameter based on the vegetal cover and catchment’s land use, and it varies from 0.05 to 1; values close to 0 indicate low soil protection while values close to 1 mean high soil protection [52].
Y is the coefficient of soil erodibility, which depends on the pedological and lithological characteristics of the watershed and indicates the resistance of soils to erosion. The values of this dimensionless factor can be determined either by laboratory experiments or by field measurements [48,52] and usually ranges from 0.25 to 2. Values close to 0.2 indicate low erodibility, whereas values close to 2 represent strongly erodible formations [52,53]. In this study, we used the soil map data and the geological map of the study area to evaluate the Y factor. Subsequently, coefficient values were assigned to each soil type based on the EPM guidelines and previous studies [52,54,55,56,57].
φ is a coefficient that depends on the active erosion and on the degree of extension of linear erosion and mass movements. It is a dimensionless factor whose values range from 0.1 to 1 [58,59]. Each type of erosion form has taken on a value of Φ according to the guidelines of the EPM method [52,54,57]. In calculating Z, Y, Xa, and φ, the weighted average value was calculated all the polygons with various values [26].
Isr represents the average mean basin slope of the studied area expressed as a percentage. This parameter is extracted from the Digital Terrain Model of the catchment (DTM-ASTER-GDEM) with a horizontal resolution of 30 m and a vertical resolution of 20 m (NASA earth data), while F is a watershed area in km2.
The total volume of sediment produced in the different areas of the watershed does not fully reach downstream. A portion is redeposited in streams or other areas of the basin; therefore, it is essential to calculate the specific real sediment production (Gsp) in m3 km−2 y−1 by the following equation:
G s p = W y e a r · R u
This is possible by multiplying the average annual production of sediments W y e a r by a delivery ratio coefficient Ru, as calculated in Equation (5).
R u = O × D 0.25 × ( L v + 10 )
Lv represents the length of the main river in Km, O is the perimeter of the catchment area in Km, and D is the average height of the basin in relation to the closing section [26].

3.2. IntEro Model Verification

In order to support the idea of using the model for watershed management planning in the study region, it was necessary to perform model validation and verification. To validate the sediment yield calculations for the Polimlje region, which includes the Miocki Potok basin, bathymetric measurements were carried out in the Potpec reservoir on the river Lim in 2017 to calculate the accumulation of sediment in the 1967–1976; 1967–1981; 1967–1986; 1967–1991; 1967–1999; 1967–2006; 1967–2012; 2012–2017 intervals.
The “Potpec” reservoir is about 17 km long, with a maximum depth of about 30 m close to the river dam. Over the years, it has been receiving large amounts of sediment that have reduces its useful volume. In order to determine the quantities of deposited sediment, a survey of the bottom of the Lim riverbed was organized in the period of 11–26 April 2017.
Field measurements of the accumulation in 2012 were performed with professional hydrographic equipment. We used a GPS rover, which provides measurement accuracy of 10 mm horizontally, and 20 mm vertically, and a Trimble R6 base with an internal radio modem, while the depth in the accumulation was measured with a single-frequency portable echo sounder Oda Hydrotrac which provides measurement accuracy of 1 cm. The equipment is installed on a boat that is structurally adapted for the mentioned type of such research work. The trigonometric points closest to the accumulation on the left and right banks were used as the geodetic basis for the survey.
After the field recording, we developed a digital terrain model of the bottom of the accumulation. From the digital terrain model, transverse profiles are drawn whose position is defined by coordinates.

4. Results

4.1. Rainfall Change

Rainfall and overland flow are key factors in soil particle detachment, entrainment, and transport. Overland flow depth depends on several factors, including soil field capacity, slope gradient, and the volume and intensity of precipitation. The rainfall factor is, in fact, included in almost all the most renowned models to predict soil erosion and sediment yield, such as the Modified Universal Soil Loss Equation (MUSLE), the Unit Sediment graph (USG), the Water Erosion Prediction Project (WEPP), etc.. Rainfall is the soil erosion triggering factor and, after land use/cover, it is the most variable, especially over short timescales.
In the study area, long time series of rainfall data (1950–2016) is available only for two meteo-stations: Berane and Bijelo Polje. Soil erosion is not influenced by the rainfall amount and mostly by rainfall intensity. For this study, unfortunately, only daily precipitation, for Bijelo Polje (1966–2019), and maximum snow thickness and duration data, For Berane, Bijelo Polje, and Plav, are available. Langbein and Schumm [60] were the first to investigate the relationship between mean annual precipitation and sediment yield. Since that pioneering paper, very many other studies on the prediction of sediment yield followed, and many of them used annual precipitation as an entry parameter [61,62]. The variability of annual precipitation for Berane and Bijelo Polje is reported in Figure 3. Berane does not show any trend, whereas the Bijelo Polje diagram is characterized by a moderate increasing trend of 99 mm across a 66-year interval.
From the daily rainfall data set two time series were extracted: the mean monthly maximum daily intensity and the absolute highest daily intensity of each year. These two time series are reported in Figure 4, which shows clearly that there was no trend for the average monthly maximum daily intensity and a very moderate increase of only 9 mm in 24 hours across the 54-year interval for the annual peak values of daily intensity. The highest precipitation of 158 mm in one day was recorded in 1992; other peaks around 100 mm in 24 hours were recorded in 1972, 1986, and 1996.
Runoff generation and volume are largely influenced by rainfall intensity; however, the threshold for soil particles detachment and beginning of sediment supply to the drainage network varies for different environments according to the climatic conditions and vegetation cover type and density. Though in his plot experiment, sediment yield was not measured [63], showed that cycles of intermittent rain with an intensity of 10 mm/h for periods of 3–25 minutes are able to generate an appreciable runoff of 1–7 mm depth.
Some photos in Dunkerley’s paper [63] show evidence of soil saturation, and soil particle entrain met should be expected. Based on Dunkerley’s experiments and on the author’s experience of field observations during rainstorms, a daily intensity of 20 mm/24 h was considered as an appropriate threshold for the generation of suspended sediment transport events in the study catchment. In order to surrogate the information about rainfall intensity at an hourly scale, we built a time series of the number of days with daily rainfall higher than 20 mm and investigated the occurrence of any trend (Figure 5).
The data indicated no trend, but the diagram in Figure 5 shows that in the mid-late 1970s, more than 10 days in a year with rainfall depth of 200 mm or more were rather common. From 1990 to 2006, the highest number of days (19) and several years with 13–16 days are recorded. In the last decade (2007–2019), most of the data plotted below ten days and lower than average sediment yield values should be expected. In order to further investigate the role of annual precipitation on sediment yield of the study catchment, the Modified Fournier Index-MFI [64] was also calculated as:
MFI = p2/P
In which p is the amount of rain in a month (mm), and P is the annual precipitation (mm).
The MFI has been developed and used as a surrogate rainfall erosivity index, and it is classified into five classes in which the values less than 90 indicate low erosivity; between 90 and 16 indicate moderate erosivity, and values beyond 160 indicate high and very high erosivity.
The time series of the MFI calculated for Berane and Bijelo Polje across the 1950–2016 interval are reported in Figure 6. This diagram shows that most of the MFI values were included in the moderate rainfall erosivity class and that no significant trend was present. High erosivity peaks occurred only in 1974, 1985, and in 1992. It is worth noticing that the latter two peaks corresponded to the peaks of absolute maximum daily intensity (Figure 6).
This index is used to represent the rainfall erosivity. The dashed lines indicate the transition between low to moderate and from moderate to high erosivity. The solid black lines are the trend lines; their interpolating equations are also reported.

4.2. Demographic Changes in the River Basin of Miocki Potok

According to the 2011 census, the municipality of Bijelo Polje, covers an area of 923 km2, with 13,199 households, 137 settlements, and 46,051 inhabitants, which became 42,191 in 2020 (Figure 7). The borders of the Municipality of Bijelo Polje are also largely the borders of Montenegro to Serbia, and they extend to the tops of the mountains.
In the middle of the past century, 1267 inhabitants lived in the catchment area of the Miocki Potok that is part of the Bijelo Polje municipality (Figure 7). The average density of the population was 30.9 inhabitants per km2, which was slightly less than the density of the population of the municipality of Bijelo Polje (39.8).
Until the 1970s, population growth was registered in all settlements of Miocki Potok basin. Index rates ranged from 101.4 to 118.8 in the 1948–1953 interval, and the population continued to grow until 1961. Since the census of 1971, the number of inhabitants started to decrease in a few settlements (like in Ograde), whereas it increased in others (like in Kanje).
In 1981, all settlements had a decrease in population compared to 1971. The trend of population decline continued in all settlements. In 2011 (Source: Census 2011), the study area records a very high depopulation rate, with only 32.2% of the 1961 population. From 1970 to 2020, the decrease in population was 24.2% (Figure 8).
The population decline since the 1970s was caused by industrialization, which initially had positive effects, as the villages were relieved of latent redundancies. Workers moved primarily went to the city center of Bijelo Polje and later to other cities in Montenegro and beyond (especially to other republics of former Yugoslavia). However, as the industrialization process intensified, the countryside was increasingly neglected, and the rural exodus, instead of declining, increased over time, leading to a significant depopulation of all rural settlements, not only in the municipality of Bijelo Polje but also throughout the Northern region of Montenegro.
Migration regained intensity at the beginning of the 21st century. In many villages, the birth rate decreased to a minimum or to zero, which led to a lack of labor force and economic decline of the local communities. The emigration regarded mainly the young who left behind villages consisting of elderly households. With such a marked change in the age composition of the villages and equivalent marked change in the land use/cover was inevitable.

4.3. Vegetation and Land Use

In Montenegro, vegetation distribution is strongly influenced by (i) the distance to the Mediterranean Sea (continental or maritime climate) and (ii) the altitude (temperate or mountain climate). Due to Montenegro’s rugged topography, there are large differences in vegetation over short distances [65]. The Miocki Potok river basin is located in the Dinarid Province of the middle-south-eastern European mountainous biogeographical region. Forests dominate this river basin accounting for 60% of the total vegetation cover.
Most of the river basin is covered by degraded beech forests (Fagetum montanum). The degraded forests are located near settlements and roads because of the common practice of firewood harvesting. These forests are characterized by a terminate canopy and by a large number of species of ground flora, shrubs, and lower trees. They differ from beech forests in the inner parts of the basin, especially from the sub-association Fagetum montanum typicum, which is characterized by a dense canopy. Near the villages of Milovo and Dobrakovo, especially on southern exposures, there are forests of Sessile oak and Turkish oak (Quercetum petraeae cerridis Lak.). In the lower part of the basin, a narrow belt along the river channel is covered with hydrophilic forest (Alnetea glutinosae, Salicetea herbacea).
The State-owned forest portions have a much better structure than the private forests. We measured seven clusters within the basin. Position of clusters in the river basin study area presented in the Figure 9.
Data obtained by field measurements showed the prevalence of degraded forests. The average timber volume in the forests is 111 m3 ha−1, with an average of 327 trees per hectare. In five representative clusters, the dominant species is beech. Beech wood volume varies from 81.5 m3 ha−1 to 331 m3 ha−1. In the Miocki Potok basin, we also recorded Betula verrucosa, Quercus cerris, Carpinus betulus, Acer campestre, Quercus petraea, and Prunus avium.
In recent years, the migration of the rural population to towns resulted in an increase in forest areas with decreasing processes of forest degradation and devastation.
In order to describe the impact of land use on soil erosion in the Miocki Potok catchment, we used a series of data from The First National Forest Inventory of Montenegro [66], Forest Management Plans, Cadastre, Landsat images, and Statistical Yearbooks to quantify variations in land use across six intervals: 1970, 1980, 1990, 2000, 2010, and 2020.
The data indicate that forests cover about half of the study basin and that the forest area has been slightly increasing during the last decades (Figure 10). The findings of this study are in line with observations of Nyssen [67], who stated that there is marked forest regrowth in the former Yugoslavia. Overall, forest cover has been naturally increasing in this region.
This denser vegetation has led to higher water infiltration into the soil. The process of forest regrowth may be further accelerated by applying good silviculture measures and forest management practices. Meadows, pastures, and orchards cover around 40% of the river basin, but this land cover type has been decreasing in recent years. According to Grimes [68], some marginal lands were cultivated again in the first decade of this century, as a consequence of the economic crisis in the region.
This trend was also confirmed by our research on land use in the small Miocki Potok basin. For the period of 2000–2020, a decrease in the area of arable and cultivated land was recorded. Land use changes in the Miocki Potok basin are presented in Figure 11.
Within the basin area, degraded forests are the most widespread form of vegetation cover (29.0%, in 1990) and well-developed forests (28.3 in 2020). Other types are as follows: mountain pastures (21.5%, on average for the period 1970–2020), meadows (16.1% on average for the period 1970-2020), plow-lands (12.6% on average for the period 1970–2020), and orchards (less than 1%, on average for the period 1970–2020) (Table 4).
The coefficient of the vegetation cover, S2, was calculated to range from 0.72 to 0.74, the coefficient of the river basin planning, Xa, from 0.45 to 0.50 Figure 12 and Figure 13.

4.4. Current Erosion and the Impact of Land Use on Soil Erosion Intensity

The dominant erosion form in the study area is from surface runoff, including severe forms of erosion such as rills, gullies, and ravines. Surface erosion has affected all soils and slopes, and it is most pronounced on steep slopes with scarce or degraded vegetation cover (Figure 14).
The IntErO model was used to calculate the soil erosion intensity and the peak discharge for the Miocki Potok river basin. The results of the impact of land use changes on soil erosion intensity in the Miocki Potok river basin in six time profiles (years: 1970, 1980, 1990, 2000, 2010, and 2020) are presented in the Table 5 and Table 6 and Figure 15 and Figure 16. Calculated soil losses in 57 river basins of the Polimlje region are presented in Table 7, as well.
Asymmetry of the river basin, a, is calculated at 0.55, and, according to that result, there is a possibility for large flood waves to appear in the Miocki Potok basin. Soil loss, Gyear, was calculated as 13680 m3 year−1 and the specific soil loss is 333 m3 km−2 year−1 (240 t km−2 year−1). The value of the Z coefficient of 0.439 indicates that the river basin belongs to destruction category III. The strength of the erosion process was moderate, and, according to the erosion type, was mixed erosion.
According to our analysis, land use changes, specifically the change of the forest structure (fs), as well as the slight decrease of the plow-lands (fg), influenced the decrease of the soil erosion intensity of 14% in the Miocki Potok River Basin.
In parallel, the value of peak discharge, Qmax in the 1970–2020 interval has decreased by 3.5% (Figure 16).

4.5. Lim River Reservoir Sedimentation Survey

From the bathymetry of Potpec reservoir, it can be seen that the backfilling of the reservoir shows a declining trend compared to previous measurement periods (Table 8). The small amount of sediment in the reservoir can also be affected by the occasional flush of water through the lower dam gates.
In order to carry out the IntErO model verification for the study area, the average sediment yields of the 57 river basins of the Polimlje region in Montenegro were calculated as 347,273 m3 year−1 [13]. For the measured annual backfill for 2017 is 340,000 m3 per year (by the team of authors, Figure 17, Figure 18 and Figure 19; see also Supplementary Materials, Figure S1: Maps of the bathymetric profiles, Potpec reservoir), whereas the value obtained by the IntErO model is 330,000 m3 per year for 2020.
The Potpec reservoir undertakes the Polimje region, which includes the Miocki Potok basin.
The difference between measured accumulation and predicted sediment yield is less than 1%. This result indicates that the use of the IntErO model is acceptable for the preparation of land use scenarios on erosion intensity and runoff in the Polimlje basins and similar basins in the region; though factors such as the trap efficiency of the Potpec reservoir and the bottom gates flush routing should also be considered for more accurate predictions. The proportion of the incoming sediment that actually accumulates in a reservoir depends on many factors, including the reservoir trap efficiency, reservoir capacity, inflow discharge, and reservoir length.
Brune [74] and Churchill [75] presented curves relating to the percentage of sediment that is likely to be retained in a reservoir to the ratio of the reservoir capacity to the mean annual runoff. Roberts [76] modified Churchill’s curve to determine the trap efficiency of a reservoir with respect to a sedimentation index, in which the reservoir length is an important factor. The Potpec reservoir is very long (about 17 km) and very narrow (width ranging between 100 and 200 m). The longitudinal profiles showed that only in the very vicinity of the dam was there evidence of bottom-gate sediment flushing. On the basis of these considerations, though we do not have data of sediment transport beyond the dam, we can assume that likely 80 to 90% of the incoming sediment was retained in the reservoir. The sedimentation rate of the Potpec reservoir, obtained from bathymetric measurements, was 340.000 m3 yr−1, so the actual sediment input can be assessed to be about 390–400,000 m3 yr−1. The average sediment input predicted by the IntEro model was slightly smaller (347,273 m3 yr−1), but given the size of the Lim river catchment with respect to that of the individual 57 small catchments used for the model calibration, the model predictions could be considered accurate and reliable, and the IntErO model could be used to quantify the sediment yield variations induced by the land use changes associated with the demographic changes in the Polimje region and the Miocki Potok basin as well.

5. Discussion

Over the last six decades, anthropogenic factors have significantly increased the pressure on agricultural and forest land in Montenegro, degrading the vegetation cover, which eventually results in serious degradation and loss of fertile soil. In Montenegro, water erosion is the most important erosion type. Water erosion is primarily caused by precipitation and consecutive runoff, but also by fluvial erosion in streams [77]. Given the extreme precipitation values in some parts of the country (the highest of Europe), the influence of this erosion type on the landscape is enormous.
According to Nyssen [67], land-use changes in Montenegro show that, on average, the area covered by dense vegetation increased from about 35% in the early-20th century to 56% of today. In the mountain region, the wooded areas increased slightly. The share of agricultural land, represented by meadow and farmland, remained constant. This pattern is similar to our findings, whereby the forests cover increased from 46% in 1970 to 52% in 2020. Our results for forest land were about four percent less than those of Nyssen [67] because the Miocki Potok study area had a higher than average number of villages and population density.
Depending on available equipment and accessibility of the terrain, different verification methods can be applied. With the idea to compare the data processed by the model measurements of sediment yield on erosion plots were conducted at the Rokava (Dragonja) River basin in Slovenia by Zorn and Komac in 2011 [23]. In North Macedonia Milevski in 2008, [21] found a very good correspondence between the results obtained for the Bregalnica river basin, using the Erosion Potential Method and on-site measurements. Bagherzadeh, in 2011 [78], verified model outputs by a field survey using GPS and a visual comparison of areas characterized as areas with moderate and heavy annual sediment yields. Kouhpeima, in 2011 [33], analyzed five different catchments in Semnan Province in Iran (Amrovan, Atary, Ali Abad, Ebrahim Abad, Royan catchments) and used comparison results to measure sediment deposits in the reservoir as a verification method. Haghizadeh in 2009 [79] used a comparison of model output results with field observations and a GLASGOD (Global Assessment of Soil Degradation) map as a verification method. In the Prescudin catchment in Italy, Bemporad in 1997 [80] measured the sediment deposits, showing minimum deviation between predicted and measured sediment yield values.
Nuclear probes for suspended-load measurements were used at the Esino and Musone river basin in Italy by Tazioli [22]. The measurements exhibited some deviations in comparison with the overall sediment yield production estimated with the EPM method, but, overall, a good order-of-magnitude correspondence was obtained concerning yearly sediment yield. It was concluded that further measurements were necessary because the EPM considers total sediment load, whereas the measurements were conducted to take into account suspended load only.
Other verification methods include the use of a personal digital assistant (PDA) device and on-site observations [18], and certain verification methods remain unspecified in papers [81,82] but provide poor overall ratings for the EPM method, which is said to overestimate sediment yield [81].
Surface erosion has affected all soils and slopes, and it is most pronounced on steep slopes with scarce or degraded vegetation cover. A similar situation has been described for Brazilian watersheds, where higher slopes [83,84,85] and poorly designed rural roads [85] are the main reason for increasing surface erosion. Under the same conditions of soil, slope, and vegetation cover, Avanzi [86] reported that a more concave slope could concentrate the runoff flow leading to a greater erosion rate. These changes can significantly affect the further purpose and use of land, especially in the future agricultural and overall economic development [87].
The calculated erosion rates were rather low at international standards, and that explained the nature of the catchment Table (Table 6 and Table 7). Besides, Babic Mladenovic [88] calculated real soil losses for the total Lim River basin at 350 m3 km−2 year−1. By using the IntErO model, we found the average value was 331.78 m3 km−2 year−1 for the 57 river basins of Polimlje in Montenegro [13], and in the range of 394 m3 km−2 yr−1 (1990–2000) to 333 m3 km−2 yr−1 (2020) for the Miocki Potok river basin (Figure 15). This correspondence suggests that the results of the assessment obtained by the IntErO model were appropriate for the study area.
In parallel, the value of peak discharge, Qmax in the 1970–2020 interval, decreased by 3.5% (Figure 16). This is in line with the findings of Nyssen [67] that in Montenegro, in recent times, runoff has decreased due to an increase in vegetation cover.
The calculated erosion rates are relatively low compared to other Mediterranean regions. That is likely related to the large occurrence of forest in the Polimje region and especially in the Miocki Potok basin.
The demographic changes observed in the study area are similar and with similar effects to those recorded in other Mediterranean areas. In the northern Adriatic rivers, for instance, the sediment flux to the sea was much reduced because of the reforestation, both spontaneous and controlled by the forest authorities, following the abandonment of cultivated areas, especially at in the mountain areas for the fast industrialization of Italy and the demand for a larger workforce [89].

6. Conclusions

The scientific literature about the impact of land use/cover change on the hydrological systems in the Mediterranean is rather abundant. Typically, studies focussed on the influence of different cultivations on runoff generation and soil loss. Others investigated the role of deforestation in increasing the risk of flooding and/or land degradation. Very few papers report about the effect of vegetation regeneration on sediment yield and hydrology at the catchment scale. This study used the IntErO model to quantify the sediment yield and runoff response of a small watershed in Montenegro to past variations of land use and natural forest regeneration, driven by rural population migration. The results of this study, though obtained for a small catchment in a mountainous area of Montenegro, provide an insight on the interconnection between social issues and hydrological processes that can be used to depict realistic scenarios of soil erosion and hydrological changes in other regions with a similar bio-physical setting. The approach proposed in this study proved to be able to indicate with a good level of accuracy the sediment yield and hydrological changes that may be expected if vegetation is allowed to regrow.
In the study area, the migration of rural population to towns proved to be a crucial factor. In fact, the resulting abandonment of farmland in recent years resulted in an increase in the forested area. In the Miocki Potok watershed, whose biophysical characteristics can be considered representative of the whole region of the upper Lim River Basin, the forests increased from 46% in 1970 to 52% in 2020. Such an increase in the forested area resulted in a decrease of 14% and 3.5% for soil loss and peak discharge, respectively, which is in line with the findings of Nyssen et al. (2014) showing that in Montenegro, recent runoff decreases are due to vegetation increase. During the 1970–2020 interval, in fact, no significant change in annual precipitation was found. Yet, no relevant change in both the daily precipitation maximum intensity and the number of days with rainfall higher than 20 mm in a year was observed. These results proved that during the 1970–2020 interval, rainfall had a secondary role compared to the forest regrowth in the predicted hydrological changes. The IntErO model accuracy in predicting sediment yield was tested against the sediment accumulation within the Potpec reservoir on the Lim River. The IntErO model predicted an average annual sediment yield of 330,000 m3, whereas from repeated topographic filed surveys, an average annual sedimentation rate in the reservoir of 340,000 m3. The Potpec reservoir is very long and narrow, and the average volume of the annually accumulated sediment was assessed by Roberts [76] equation around 400,000 m3 yr−1. This result and the IntErO model prediction were very close, thus confirming the usefulness of the IntErO model to calculate runoff and sediment yield at a watershed scale in southeastern Europe, and areas with physical characteristics similar to the Polimlje Region of Montenegro. In the study area, the calculated erosion rates are relatively low compared to other Mediterranean regions. This result can be accounted for by the large proportion of forested areas in the Polimje region and especially in the Miocki Potok basin.
Seven longitudinal profiles of the Potpec reservoir bottom were surveyed in 1967, 1976, 1981, 1999, 2000, 2012, and 2017. The largest sediment accumulation rates occurred from 1967 to 1981, then a period of decreasing sedimentation followed with very low rates observed in the last two decades. This pattern, though with a decade delay, very well reflects the migration phases of the rural population and the consequent abandonment of their farmlands.
The fast industrialization of Montenegro and the growing need for a larger workforce produced deep demographic changes, especially in the rural, mountain areas. The consequent abandonment of the cultivated areas resulted in a natural process of reforestation whose hydrological implications are well predicted by the IntErO model.
Though the approach proposed in this study proved to be successful for the Balkan region, other studies in other countries are necessary to define the environmental constraints of the model application.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4441/12/11/2973/s1, Figure S1: Maps of the bathymetric profiles, Potpec reservoir.

Author Contributions

Conceptualization, V.S., G.B., M.C., N.D., B.D., and P.B.; methodology, V.S., M.C., N.D., B.D., and P.B.; software, V.S.; validation, V.S., G.B., and D.V., and P.B.; formal analysis, V.S., M.C., N.D., B.D., and P.B.; investigation, V.S., G.B., D.V., and M.C.; resources, V.S., and B.D.; data curation, V.S., D.V., and M.C.; writing—review and editing, V.S., M.C., M.B., and P.B.; writing—review and editing, V.S., G.B., D.V., M.C., M.B., N.D., B.D., and P.B.; visualization, V.S.; supervision, V.S. and P.B.; project administration, V.S.; funding acquisition, V.S. and B.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministry of Science of Montenegro.

Acknowledgments

The authors are grateful to Jan Nyssen and Amaury Frankl from the Ghent University (Belgium) for the initial support in the conceptualization of this research, but also for providing us the literature and other valuable support; to the assistance of Slavisa Lucic and Aleksandar Stijovic of the Institute of Forestry of Montenegro for providing the specific inputs on forestry management and thematic maps; to Slobodan Radusinovic from the Institute of Geology (Montenegro) for providing and commenting some thematic geological maps; to Ernad Suljevic Deputy Mayor of Bijelo Polje who provides us data on migration for the villages Ograde, Mioce, Dobrinje, Miovo, and Kanje. Authors are grateful to Will Mahoney (United States), Soil Fund Ambassador of the International Erosion Control Association (IECA), for his valuable comments and suggestions for improving the text of the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gudino-Elizondo, N.; Biggs, T.W.; Bingner, R.L.; Langendoen, E.J.; Kretzschmar, T.; Taguas, E.V.; Taniguchi-Quan, K.T.; Liden, D.; Yuan, Y. Modelling Runoff and Sediment Loads in a Developing Coastal Watershed of the US-Mexico Border. Water 2019, 11, 1024. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Dimotta, A.; Cozzi, M.; Romano, S.; Lazzari, M. Soil Loss, Productivity and Cropland Values GIS-Based Analysis and Trends in the Basilicata Region (Southern Italy) from 1980 to 2013; ICCSA 2016, Part IV, Lecture Notes in Computer Science, 9789; Gervasi, O., Ed.; Springer: Cham, Switzerland, 2016; pp. 29–45. [Google Scholar] [CrossRef]
  3. Dimotta, A.; Lazzari, M.; Cozzi, M.; Romano, S. Soil Erosion Modelling on Arable Lands and Soil Types in Basilicata, Southern Italy; ICCSA 2017, Part V, Lecture Notes in Computer Science LNCS, 10408; Gervasi, O., Ed.; Springer: Cham, Switzerland, 2017; pp. 57–72. [Google Scholar] [CrossRef]
  4. Chalise, D.; Kumar, L.; Spalevic, V.; Skataric, G. Estimation of Sediment Yield and Maximum Outflow Using the IntErO Model in the Sarada River Basin of Nepal. Water 2019, 11, 952. [Google Scholar] [CrossRef] [Green Version]
  5. Hazbavi, Z.; Sadeghi, S.H.; Gholamalifard, M. Dynamic analysis of soil erosion-based watershed health. Geogr. Environ. Sustain. 2019, 12, 43–59. [Google Scholar]
  6. Nikolic, G.; Spalevic, V.; Curovic, M.; Khaledi Darvishan, A.; Skataric, G.; Pajic, M.; Kavian, A.; Tanaskovik, V. Variability of Soil Erosion Intensity Due to Vegetation Cover Changes: Case Study of Orahovacka Rijeka, Montenegro. Not. Bot. Horti Agrobot. Cluj-Napoca 2018, 47, 237–248. [Google Scholar] [CrossRef] [Green Version]
  7. Ballesteros-Cánovas, J.A.; Czajka, B.; Janecka, K.; Lempa, M.; Kaczka, R.J.; Stoffel, M. Flash floods in the Tatra Mountain streams: Frequency and triggers. Sci. Total Environ. 2015, 511, 639–648. [Google Scholar] [CrossRef]
  8. Verheijen, F.G.A.; Jones, R.J.A.; Rickson, R.J.; Smith, C.J. Tolerable versus actual soil erosion rates in Europe. Earth Sci. Rev. 2009, 94, 23–38. [Google Scholar] [CrossRef] [Green Version]
  9. Dragicevic, N.; Karleusa, B.; Ozanic, N. Different Approaches to Estimation of Drainage Density and Their Effect on the Erosion Potential Method. Water 2019, 11, 593. [Google Scholar] [CrossRef] [Green Version]
  10. Kisic, I.; Basic, F.; Nestroy, O.; Sabolic, M. Global Change-Challenges for soil management, Advances in Geoecology. In Soil Erosion under Different Tillage and Cropping Systems in Central Croatia; Catena Verlag GMBH: Reiskirchen, Germany, 2010; pp. 141–150. [Google Scholar]
  11. Glavan, M.; Milicic, V.; Pintar, M. Finding options to improve catchment water quality Lessons learned from historical land use situations in a Mediterranean catchment in Slovenia. Ecol. Model. 2013, 261, 58–73. [Google Scholar] [CrossRef]
  12. Volk, M.; Liersch, S.; Schmidt, G. Towards the implementation of the European Water Framework Directive? Lessons learned from water quality simulations in an agricultural watershed. Land Use Policy 2009, 26, 580–588. [Google Scholar] [CrossRef]
  13. Spalevic, V. Impact of Land Use on Runoff and Soil Erosion in Polimlje. Ph.D. Thesis, Faculty of Agriculture, University of Belgrade, Belgrade, Serbia, 2011; pp. 1–260. [Google Scholar]
  14. Kostadinov, S.; Braunović, S.; Dragićević, S.; Zlatić, M.; Dragović, N.; Rakonjac, N. Effects of Erosion Control Works: Case Study—Grdelica Gorge, the South Morava River (Serbia). Water 2018, 10, 1094. [Google Scholar] [CrossRef] [Green Version]
  15. Gocić, M.; Dragićević, S.; Radivojević, A.; Martić Bursać, N.; Stričević, L.; Đorđević, M. Changes in Soil Erosion Intensity Caused by Land Use and Demographic Changes in the Jablanica River Basin, Serbia. Agriculture 2020, 10, 345. [Google Scholar] [CrossRef]
  16. Gavrilovic, S. Engineering of Torrential Flows and Erosion; Izgradnja: Beograd, Serbia, 1972; 272p. [Google Scholar]
  17. Dragicevic, S.; Pripuzic, M.; Zivkovic, N.; Novkovic, I.; Kostadinov, S.; Langovic, M.; Milojkovic, B.; Cvorovic, Z. Spatial and Temporal Variability of Bank Erosion during the Period 1930–2016: Case Study-Kolubara River Basin (Serbia). Water 2017, 9, 748. [Google Scholar] [CrossRef] [Green Version]
  18. Tosic, R.; Dragicevic, S.; Zlatic, M.; Todosijevic, M.; Kostadinov, S. The impact of socio-demographic changes on land use and soil erosion (case study: Ukrina River catchment). Geogr. Rev. 2012, 46, 69–78. [Google Scholar]
  19. Kaloper, S.E.; Čadro, S.; Uzunović, M.; Cherni-Čadro, S. Determination of erosion intensity in Brka watershed, Bosnia and Herzegovina. Agric. For. 2020, 66, 79–92. [Google Scholar] [CrossRef]
  20. Dragičević, N.; Karleuša, B.; Ožanić, N. GIS based monitoring database for Dubračina river catchment area as a tool for mitigation and prevention of flash flood and erosion. In Proceedings of the thirteenth International Symposium on Water Management and Hydraulic Engineering, Bratislava, Sovakia, 9–12 September 2013; pp. 553–565. [Google Scholar]
  21. Milevski, I.; Blinkov, I.; Trendafilov, A. Soil erosion processes and modelling in the upper Bregalnica catchment. In Proceedings of the Conference of the Danubial Countries on the Hydrological Forecasting and Hydrological Bases of Water Management (XXIV), Bled, Slovenia, 2–4 June 2008; pp. 1–10. [Google Scholar]
  22. Tazioli, A. Evaluation of erosion in equipped basins: Preliminary results of a comparison between the Gavrilovic model and direct measurements of sediment transport. Environ. Geol. 2009, 56, 825–831. [Google Scholar] [CrossRef]
  23. Zorn, M.; Komac, B. The importance of measuring erosion processes on the example of Slovenia. Hrvat. Geogr. Glas. 2011, 73, 19–34. [Google Scholar] [CrossRef] [Green Version]
  24. Kayimierski, L.D.; Irigoyen, M.; Re, M.; Menendey, A.N.; Spalletti, P.; Brea, J.D. Impact of climate change on sediment yield from the upper Plata basin. Int. J. River Basin Manag. 2013, 11, 411–421. [Google Scholar]
  25. Tavares, A.S.; Spalevic, V.; Avanzi, J.C.; Nogueira, D.A.; Silva, M.L.N.; Mincato, R.L. Modeling of water erosion by the erosion potential method in a pilot subbasin in southern Minas Gerais. Semin. Ciências Agrárias 2019, 40, 555–572. [Google Scholar] [CrossRef] [Green Version]
  26. Ouallali, A.; Aassoumi, H.; Moukhchane, M.; Moumou, A.; Houssni, M.; Spalevic, V.; Keesstra, S. Sediment mobilization study on Cretaceous, Tertiary and Quaternary lithological formations of an external Rif catchment, Morocco. Hydrol. Sci. J. 2020, 65, 1568–1582. [Google Scholar] [CrossRef]
  27. Ali, S.S.; Al-Umary, F.A.; Salar, S.G.; Al-Ansari, N.; Knutsson, S. GIS Based Soil Erosion Estimation Using EPM Method, Garmiyan Area, Kurdistan Region, Iraq. J. Civ. Eng. Archit. 2016, 10, 291–308. [Google Scholar] [CrossRef] [Green Version]
  28. Hazbavi, Z.; Azizi, E.; Sharifi, Z.; Alaei, N.; Mostafazadeh, R.; Behzadfar, M.; Spalevic, V. Comprehensive estimation of erosion and sediment components using IntErO model in the KoozehTopraghi Watershed, Ardabil Province. Environ. Eros. Res. J. 2020, 10, 92–110. [Google Scholar]
  29. Barovic, G.; Silva, M.L.N.; Batista, P.V.G.; Vujacic, D.; Soares Souza, W.; Avanzi, J.C.; Behzadfar, M.; Spalevic, V. Estimation of sediment yield using the IntErO model in the S1-5 Watershed of the Shirindareh River Basin, Iran. Agric. For. 2015, 61, 233–243. [Google Scholar] [CrossRef] [Green Version]
  30. Behzadfar, M.; Tazioli, A.; Vukleic-Shutoska, M.; Simunic, I.; Spalevic, V. Calculation of Sediment yield in the S1-1 watershed, Shirindareh watershed, Iran. Agric. For. 2014, 60, 207–216. [Google Scholar]
  31. Moradi, S.; Limaei, S.M.; Khanmohammadi, M. Calculation of sediment yield in the Zemkan River Basin of Iran using analytical methods and GIS concept. Agric. For. 2015, 61, 157–171. [Google Scholar] [CrossRef]
  32. Amini, H.; N’Jalaliyan, A.H.; Khalilizadeh, M.; Baharlouie, J. A comparison of EPM and WEPP models for estimating soil erosion of Marmeh Watershed in the South Iran. Agric. For. 2014, 60, 299–315. [Google Scholar]
  33. Kouhpeima, A.; Hashemi, S.A.A.; Feiznia, S. A study on the efficiency of erosion potential model (EPM) using reservoir sediments. Elixir Pollut. 2011, 38, 4135–4139. [Google Scholar]
  34. Sadeghi, S.H. Semi-Detailed Technique for Soil Erosion Mapping Based on BLM and Satellite Image Applications. J. Agric. Sci. Technol. 2005, 7, 133–142. [Google Scholar]
  35. Poornazari, N.; Khalilimoghadam, B.; Hazbavi, Z.; Bagheri, M. Land degradation assessment in the dust hotspot of southeastern Ahvaz, Iran. Land Degrad. Dev. 2020. Accepted Author Manuscript. [Google Scholar] [CrossRef]
  36. Jorda-Capdevila, D.; Casals, L.C. Water Management in the Media and Research: Dissemination in Catalonia and Its Capture by Private Companies. Int. Multidiscip. J. Soc. Sci. 2019, 8, 267–288. [Google Scholar] [CrossRef]
  37. Dimitrijevic, M.D. Coloured map of the Geology of Serbia, with overlay and legend. In Geological Atlas of Serbia, 1:2,000,000; Grafimex: Belgrade, Serbia, 1992. [Google Scholar]
  38. Zivaljevic, M. Tumac Geoloske Karte SR Crne Gore, 1:200 000; Posebna izdanja Geoloskog Glasnika, Knjiga VIII; Titograd, Montenegro, 1989. [Google Scholar]
  39. Fustic, B.; Djuretic, G. Soils of Montenegro [Zemljista Crne Gore]; Spalevic, V., Ed.; Biotechnical Institute, University of Montenegro: Podgorica, Montenegro, 2000; pp. 1–626. [Google Scholar]
  40. Chang, A.Y.; Parrales, M.E.; Jimenez, J.; Sobieszczyk, M.E.; Hammer, S.M.; Copenhaver, D.J.; Kulkarni, R.P. Combining Google Earth and GIS mapping technologies in a dengue surveillance system for developing countries. Int. J. Health Geogr. 2009, 8, 49. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Karkanis, P.G.; Au, K.; Schaalje, G.B. Comparison of 4 measurement schedules for determination of soil particle-size distribution by the hydrometer method. Can. Agric. Eng. 1991, 33, 211–215. [Google Scholar]
  42. Thun, R.; Herrmann, R. Die Untersuchung von Boden; Neumann Velag: Radebeul/Berlin, Germany, 1949; pp. 15–28. [Google Scholar]
  43. Jakovljevic, M.; Pantovic, M.; Blagojevic, S. Laboratory Manual of Soil and Water Chemistry (in Serbian); Faculty of Agriculture: Belgrade, Serbia, 1995. [Google Scholar]
  44. Egner, H.H.; Domingo, W.R. Untersuchungen über die chemische Bodenanalyse als Grundlage für die Beurteilung des Nährstoffzustandes der Böden. II. Chemische Extraktionsmethoden zur Phosphor-und Kaliumbes-timmung. Kungliga Lantbrukshˆgskolans Annaler 1960, 26, 199–215. [Google Scholar]
  45. Kappen, H. Die Bodenazidität; Springer: Berlin/Heidelberg, Germany, 1929; 363p. [Google Scholar]
  46. Blinkov, I.; Kostadinov, S. Applicability of various erosion risk assessment methods for engineering purposes. In Proceedings of the BALWOIS 2010 Conference, Ohrid, North Macedonia, 25–29 May 2010. [Google Scholar]
  47. Ristic, R.; Radic, B.; Vasiljevic, N.; Nikic, Z. Land use change for flood protection—A prospective study for the restoration of the river Jelasnica watershed. Bull. Fac. For. 2011, 103, 115–130. [Google Scholar] [CrossRef]
  48. Spalevic, V. Application of Computer-Graphic Methods in the Studies of Draining Out and Intensities of Ground Erosion in the Berane Valley. Master’s Thesis, Faculty of Agriculture of the University of Belgrade, Belgrade, Serbia, 1999; 135p. [Google Scholar]
  49. Vujacic, D.; Barovic, G.; Tanaskovikj, V.; Kisic, I.; Song, X.; Silva, M.L.N.; Spalevic, V. Calculation of runoff and sediment yield in the Pisevska Rijeka Watershed, Polimlje, Montenegro. Agric. For. 2015, 61, 225–234. [Google Scholar] [CrossRef]
  50. Vujacic, D.; Spalevic, V. Assessment of Runoff and Soil Erosion in the Radulicka Rijeka Watershed, Polimlje, Montenegro. Agric. For. 2016, 62, 283–292. [Google Scholar] [CrossRef] [Green Version]
  51. Spalevic, V.; Dlabac, A.; Spalevic, B.; Fustic, B.; Popovic, V. Application of computer—Graphic methods in the research of runoff and intensity of ground erosion—I program “River basins”. Agric. For. 2000, 46, 19–36. [Google Scholar]
  52. Efthimiou, N.; Lykoudi, E. Soil erosion estimation using the EPM model. In Proceedings of the 14th International Conference of the Geological Society of Greece, Thessaloniki, Greece, 25–27 May 2016; Volume L. [Google Scholar]
  53. Shahabi, H.; Salari, M.; Bin Ahmad, B.; Mohammadi, A. Soil Erosion Hazard Mapping in Central Zab Basin Using EPM Model in GIS Environment. Int. J. Geogr. Geol. 2016, 5, 224–235. [Google Scholar] [CrossRef]
  54. Noori, H.; Siadatmousavi, S.M.; Mojaradi, B. Assessment ofsediment yield using RS and GIS at two sub-basins of Dez Watershed, Iran. Int. Soil Water Conserv. Res. 2016, 4, 199–206. [Google Scholar] [CrossRef] [Green Version]
  55. Kojchevska, T. EPM for Soil Loss Estimation in Different Geomorphologic Conditions and Data Conversion by Using GIS. Ph.D. Thesis, Mediterranean Agronomic Institute, Chania, Greece, 2014. [Google Scholar]
  56. Dragicevic, N.; Karleusa, B.; Ozanic, N. Erosion Potential Method (Gavrilovic method) sensitivity analysis. Soil Water Res. 2017, 12, 51–59. [Google Scholar] [CrossRef] [Green Version]
  57. Milanesi, L.; Pilotti, M.; CleriCi, A.; GavriloviC, Z. Application of an improved version of the Erosion Potential Method in Alpine areas. Ital. J. Eng. Geol. Environ. 2015, 1, 17–30. [Google Scholar] [CrossRef]
  58. Gavrilovic, Z. Use of an empirical method (Erosion Potential Method) for calculating sediment production and transportation in unstudied or torrential streams. In International Conference on River Regime; Hydraulics Research Ltd.: Wallingford, UK, 1988; pp. 411–422. [Google Scholar]
  59. Staut, M. Recent Erosional Processes in the Catchment of the Dragonja River. Ph.D. Thesis, University of Ljubljana, Ljubljana, Slovenia, 2004. Unpublished. [Google Scholar]
  60. Langbein, W.B.; Schumm, S.A. Yield of sediment in relation to mean annual precipitation. Eos Trans. Am. Geophys. Union 1958, 39, 1076–1084. [Google Scholar] [CrossRef] [Green Version]
  61. Walling, D.E. Sediment yields and sediment budgets. In Encyclopaedia of Hydrological Sciences; Anderson, M.G., McDonnell, J.J., Eds.; Wiley: Hoboken, NJ, USA, 2005; pp. 1283–1304. [Google Scholar]
  62. White, S. Sediment yield prediction and modelling. In Encyclopaedia of Hydrological Sciences; Anderson, M.G., McDonnell, J.J., Eds.; Wiley: Hoboken, NJ, USA, 2005; pp. 1315–1326. [Google Scholar]
  63. Dunkerley, D. How is overland flow produced under intermittent rain? An analysis using plot-scale rainfall simulation on dryland soils. J. Hydrol. 2018, 556, 119–130. [Google Scholar] [CrossRef]
  64. Arnoldus, H.M.J. An approximation of the rainfall factor in the Universal Soil Loss Equation. In Assessment of Erosion; Wiley: Hoboken, NJ, USA, 1980; pp. 127–132. [Google Scholar]
  65. Wraber, T. Map of Natural Potential Vegetation of SFR Yugoslavia, Scale 1:1 000 000; Jovan Hadzi Biological Institute ZRC SAZU: Skopje, North Macedonia, 1983. [Google Scholar]
  66. Dees, M.; Andjelic, M.; Fetic, A.; Jokanovic, B.; Tepavcevic, V.; Borota, D.; Vasiljevic, A.; Weinreih, A.; Hahn, N.; Markovic, D.; et al. The First National Forest Inventory of Montenegro; Ministry of Agriculture and Rural Development, Government of Montenegro: Podgorica, Montenegro, 2013; pp. 1–347.
  67. Nyssen, J.; Van den Branden, J.; Spalevic, V.; Frankl, A.; Van de Velde, L.; Curovic, M.; Billi, P. Twentieth century land resilience in Montenegro and consequent hydrological response. Land Degrad. Dev. 2014, 25, 336–349. [Google Scholar] [CrossRef] [Green Version]
  68. Grimes, A.P.; Latif, M.A.; Vujadinovic-Colic, J.; Buskovic, V. Biodiversity Assessment Update for Serbia and Montenegro; USAID/Serbia-Montenegro: Belgrade, Serbia, 2005. [Google Scholar]
  69. Spalevic, V.; Barovic, G.; Fikfak, A.; Kosanovic, S.; Djurovic, M.; Popovic, S. Sediment yield and Land use changes in the Northern Montenegrin Watersheds: Case study of Seocki Potok of the Polimlje Region. J. Environ. Prot. Ecol. 2016, 17, 990–1002. [Google Scholar]
  70. Spalevic, V.; Curovic, M.; Tanaskovic, V.; Djurovic, N.; Lenaerts, T.; Nyssen, J. Application of the IntErO model for the assessment of the soil erosion intensity and runoff of the river basin Dragovo Vrelo, Montenegro. Turk. J. Agric. Sci. 2015, 1, 1072–1079. [Google Scholar]
  71. Spalevic, V.; Radanovic, D.; Behzadfar, M.; Djekovic, V.; Andjelkovic, A.; Milosevic, N. Calculation of the sediment yield of the Trebacka Rijeka, Polimlje, Montenegro. Agric. For. 2014, 60, 259–272. [Google Scholar]
  72. Spalevic, V.; Grbovic, K.; Gligorevic, K.; Curovic, M.; Billi, P. Calculation of runoff and soil erosion on the Tifran watershed, Polimlje, North-East of Montenegro. Agric. For. 2013, 59, 5–17. [Google Scholar]
  73. Spalevic, V.; Mahoney, W.; Djurovic, N.; Üzen, N.; Curovic, M. Calculation of soil erosion intensity and maximum outflow from the Rovacki river basin, Montenegro. Agric. For. 2012, 58, 7–21. [Google Scholar]
  74. Brune, G.M. Trap efficiency of reservoirs. Trans. Am. Geophys. Union 1953, 34, 407–418. [Google Scholar] [CrossRef]
  75. Churchill, M.A. Analysis and Use of Reservoir Sedimentation Data. In Proceedings of the Federal Interagency Sedimentation Conference (FIASC), Denver, CO, USA, 6–8 May 1947; Gottschalk, L.C., Ed.; pp. 139–140. [Google Scholar]
  76. Roberts, C.P.R. Flow Profile Calculations; Hydro 82; University of Pretoria: Pretoria, South Africa, 1982. [Google Scholar]
  77. Kostadinov, S.; Zlatic, M.; Dragovic, N.; Gavrilovic, Z. Soil Erosion in Serbia and Montenegro, Soil Erosion in Europe; Bordman, J., Poesen, J., Eds.; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2006; pp. 271–277. [Google Scholar]
  78. Bagherzadeh, A.; Daneshvae, M.R.M. Sediment yield assessment by EPM and PSIAC models using GIS data in semi-arid region. Front. Earth Sci. 2011, 5, 207. [Google Scholar] [CrossRef]
  79. Haghizadeh, A.; Shui, L.T.; Godarzi, E. Forecasting sediment with erosion potential method with emphasis on Land use changes at basin, electronic. J. Geotech. Eng. 2009, 14, 1–12. [Google Scholar]
  80. Bemporad, G.A.; Alterach, J.; Amighetti, F.F.; Peviani, M.; Saccardo, I. A distributed approach for sediment yield evaluation in Alpine regions. J. Hydrol. 1997, 197, 370–392. [Google Scholar] [CrossRef]
  81. Abadi, L.Z.; Ahmadi, H. Comparison of EPM and geomorphology methods for erosion and sediment yield assessment in Kasilian Watershed, Mazandaran Province, Iran. Desert 2011, 16, 103–109. [Google Scholar]
  82. Ghobadi, Y.; Pirasteh, S.; Pradhan, B.; Ahmad, N.B.; Shafri, H.Z.B.M.; Sayyad, G.A.; Kabiri, K. Determine of correlation coefficient between EPM and MPSIAC models and generation of erosion maps by GIS techniques in Baghmalek watershed, Khuzestan, Iran. In Proceedings of the 5th Symposium on Advances in Science and Technology SAStech, Mashhad, Iran, 12 May 2011; pp. 1–12. [Google Scholar]
  83. Silva, M.A.; Silva, M.L.N.; Curi, N.; Oliveira, A.H.; Avanzi, J.C.; Norton, L.D. Water erosion risk prediction in eucalyptus plantation. Ciência Agrotecnol. 2014, 38, 160–172. [Google Scholar] [CrossRef] [Green Version]
  84. Silva, M.A.; Silva, M.L.N.; Curi, N.; Avanzi, J.C.; Leite, F.P. Sistemas de manejo em plantios florestais de eucalipto e perdas de solo e água na região do Vale do Rio Doce, MG. Ciência Florest. 2011, 21, 765–776. [Google Scholar] [CrossRef] [Green Version]
  85. Avanzi, J.C.; Silva, M.L.N.; Curi, N.; Norton, L.D.; Beskow, S.; Martins, S.G. Spatial distributed model for assessing soil erosion risk in a small watershed. Ciência Agrotecnol. 2013, 37, 427–434. [Google Scholar] [CrossRef] [Green Version]
  86. Avanzi, J.C. Soil Properties, Conditions and Soil Losses for South and East Brazilian Forest Areas. Ph.D. Thesis, Federal University of Lavras, Lavras, Brazil, 2009; pp. 1–76. [Google Scholar]
  87. Melovic, B.; Cirovic, D.; Dudic, B.; Vulic, T.B.; Gregus, M. The Analysis of Marketing Factors Influencing Consumers’ Preferences and Acceptance of Organic Food Products—Recommendations for the Optimization of the Offer in a Developing Market. Foods 2020, 9, 259. [Google Scholar] [CrossRef] [Green Version]
  88. B Babic Mladenovic, M.; Obuskovic, Z.; Knezevic, Z. Zasipanje akumulacija u Srbiji-Problemi i pravci resavanja [Siltation in Serbia-Problems and directions of solving]. Vodoprivreda 2003, 35, 387–393. [Google Scholar]
  89. Preciso, E.; Salemi, E.; Billi, P. Land use changes, torrent control works and sediment mining: Effects on channel morphology and sediment flux, case study of the Reno River (Northern Italy). Hydrol. Process. 2012, 26, 1134–1148. [Google Scholar] [CrossRef]
Figure 1. Study area of the Miocki Potok river basin; Examples of the landscape and of the typical land use/cover in the study catchment: From the upper reaches of the Miocki Potok to the inflow of Kanjski Potok into the Lim River, with a panoramic view on the villages Milovo and Kanje.
Figure 1. Study area of the Miocki Potok river basin; Examples of the landscape and of the typical land use/cover in the study catchment: From the upper reaches of the Miocki Potok to the inflow of Kanjski Potok into the Lim River, with a panoramic view on the villages Milovo and Kanje.
Water 12 02973 g001
Figure 2. IntErO model Flow chart with the details on inputs needed (27) and basic outputs (22).
Figure 2. IntErO model Flow chart with the details on inputs needed (27) and basic outputs (22).
Water 12 02973 g002
Figure 3. Long term variability of annual precipitation in the study area.
Figure 3. Long term variability of annual precipitation in the study area.
Water 12 02973 g003
Figure 4. Long term variability of mean monthly maximum daily intensity (Mean max) and annual absolute highest daily intensity (Abs max).
Figure 4. Long term variability of mean monthly maximum daily intensity (Mean max) and annual absolute highest daily intensity (Abs max).
Water 12 02973 g004
Figure 5. Time series of the number of days with rainfall higher than 20 mm in a year. The red line discriminates the years with more than ten days with a daily intensity of 20 mm or more.
Figure 5. Time series of the number of days with rainfall higher than 20 mm in a year. The red line discriminates the years with more than ten days with a daily intensity of 20 mm or more.
Water 12 02973 g005
Figure 6. Long term variation of the Modified Fournier Index-MFI [64].
Figure 6. Long term variation of the Modified Fournier Index-MFI [64].
Water 12 02973 g006
Figure 7. Demographic changes in the Bijelo Polje Region, Montenegro.
Figure 7. Demographic changes in the Bijelo Polje Region, Montenegro.
Water 12 02973 g007
Figure 8. Demographic changes in the River Basin of Miocki Potok, Bijelo Polje Region, Montenegro.
Figure 8. Demographic changes in the River Basin of Miocki Potok, Bijelo Polje Region, Montenegro.
Water 12 02973 g008
Figure 9. Position of clusters in the river basin study area (left) and forests within the basin (right).
Figure 9. Position of clusters in the river basin study area (left) and forests within the basin (right).
Water 12 02973 g009
Figure 10. Land-use changes (period 1970–2020): fs, forests (suma); ft, grass (trava); fg, bare land (goleti).
Figure 10. Land-use changes (period 1970–2020): fs, forests (suma); ft, grass (trava); fg, bare land (goleti).
Water 12 02973 g010
Figure 11. Land use in the Miocki Potok river basin (1970, 1980, 1990, 2000, 2010, and 2020).
Figure 11. Land use in the Miocki Potok river basin (1970, 1980, 1990, 2000, 2010, and 2020).
Water 12 02973 g011
Figure 12. Coefficient of the vegetation cover, S2.
Figure 12. Coefficient of the vegetation cover, S2.
Water 12 02973 g012
Figure 13. River basin planning, coefficient of the river basin planning, Xa.
Figure 13. River basin planning, coefficient of the river basin planning, Xa.
Water 12 02973 g013
Figure 14. Collection on data on land use and on soil erosion processes from the field and using Google data; Kanje (43°08′22.0″ N 19°45′40.9″ E; 43.139430, 19.761351).
Figure 14. Collection on data on land use and on soil erosion processes from the field and using Google data; Kanje (43°08′22.0″ N 19°45′40.9″ E; 43.139430, 19.761351).
Water 12 02973 g014
Figure 15. Specific real soil losses per km2 (m3 km−2 yr−1).
Figure 15. Specific real soil losses per km2 (m3 km−2 yr−1).
Water 12 02973 g015
Figure 16. Calculated Peak discharge for the studied basin (m3 s−1).
Figure 16. Calculated Peak discharge for the studied basin (m3 s−1).
Water 12 02973 g016
Figure 17. Bathymetric survey of the Potpec reservoir and detail of the equipment used.
Figure 17. Bathymetric survey of the Potpec reservoir and detail of the equipment used.
Water 12 02973 g017
Figure 18. Bathymetric map of Potpec reservoir and representative cross-section with evidence of sediment accumulation.
Figure 18. Bathymetric map of Potpec reservoir and representative cross-section with evidence of sediment accumulation.
Water 12 02973 g018
Figure 19. Longitudinal profiles of the Potpec reservoir bottom surveyed in different years.
Figure 19. Longitudinal profiles of the Potpec reservoir bottom surveyed in different years.
Water 12 02973 g019
Table 1. Precipitation and temperature for the period 1948–2020 recorded at Bijelo Polje.
Table 1. Precipitation and temperature for the period 1948–2020 recorded at Bijelo Polje.
JanFebMarAprMayJunJulAugSepOctNovDec
Max. daily precipitation in mm
Max.68.692.87393.342.658.597.855.895.6157.6101.679.4
Aver.23.422.821.52421.920.621.721.125.3292923.5
St.D.15.618.613159.912.415.211.917.224.616.214.7
Mean monthly temperatures in °C
Max.2.95.87.812.615.818.220.820.917.712.68.64.4
Min.−5.6−5.2−0.76.19.814.116.214.311.36.2−1.6−4.7
Aver.−1.60.84.68.913.316.318.117.714.39.44.50.1
St.D.2.22.72.11.31.311.11.41.51.42.12.2
Max. daily temperatures in °C
Max.15.420.925.628.132.435.536.839.23629.52319.2
Aver.11.714.520.123.627.630.432.832.829.424.818.613.6
St.D.2.833.12.32.22.522.52.62.62.83.3
Min. daily temperatures in °C
Min.−27.6−24.5−16.5−7.5−401.22.6−4−7.2−15.4−21.7
Aver.−15.1−13−8.4−2.80.94.86.56.12.3−2.5−7.3−12.6
St.D.5.34.74.11.821.82.11.52.52.33.74.6
Source: Data from the Hydrometeorological Institute of Montenegro and the Biotechnical faculty of the University of Montenegro.
Table 2. The values of the Z coefficients classified according to the degree of erosion intensity.
Table 2. The values of the Z coefficients classified according to the degree of erosion intensity.
CategoriesErosion IntensityErosion Coefficient (Z)Average of Z
IVery SevereZ > 1.00Z = 1.25
IISevere0.71 < Z < 1.00Z = 0.85
IIIModerate0.41 < Z < 0.70Z = 0.55
IVWeak0.20 < Z < 0.40Z = 0.30
VVery WeakZ < 0.19Z = 0.10
Table 3. Values of the factors used in the EPM/IntErO model.
Table 3. Values of the factors used in the EPM/IntErO model.
Coefficient of Soil CoverX Value
Areas Without Vegetal Cover (Bare Land, Building Area, Water)0.8–0.9
Crop Fields, Meadows, Grasslands0.6–0.8
Built-up Areas and Crops, Degraded Shrublands0.4–0.6
Arboricultural Lands, Clear Shrublands0.2–0.4
Reforested Areas, Dense Forests, Dense Shrublands0.05–0.2
Coefficient of Soil ResistanceY Value
Marls, Clays, Poorly Consolidated Yellow Sands and Other Rock With Little Resistance1.3–1.7
Weak Rock, Fine Clayey Pelites With Microbereccia Beds, Recent Quaternary Scree1–1.3
Rock With Moderate Erosion Resistance, Limestone, Quaternary Fluvial Terraces0.6–1
Hard Rock, Sandstone of the Numidian Nappe0.5–0.6
Coefficient of Type and Extent of ErosionΦ Value
Deep Ravines, Landslides, Badlands Areas, and Bank Undercutting0.8–0.9
Sheet Erosion, Less than 50% of the Catchment Area With Rill and Gullies Erosion0.6–0.7
20% of the Area Attacked by Surface Erosion, Minor Slips in Stream Channels0.3–0.5
Land Surface Without Visible Erosion, Mostly Crop Fields0.1–0.2
Table 4. Land use in the Miocki Potok river Basin for the period 1970–2020.
Table 4. Land use in the Miocki Potok river Basin for the period 1970–2020.
Land Use Type197019801990200020102020
Forests-fs (sume)0.460.480.480.470.500.52
Grass-ft (trave)0.390.400.400.380.370.38
Bare Lands-fg (goleti)0.150.120.120.150.130.10
Table 5. The IntErO report of Inputs for the Miocki Potok river basin.
Table 5. The IntErO report of Inputs for the Miocki Potok river basin.
InputSymbolValueUnit
River Basin AreaF41.06km2
The Length of the WatershedO30.7km
Natural Length of the Main WatercourseLv6.33km
The Shortest Distance (the Fountainhead and Mouth)Lm6.03km
The Total Length of the Main Watercourse With Tributaries ΣL13.88km
River Basin Length Measured by A Series of Parallel LinesLb11.12km
The Area of the Bigger River Basin PartFv26.2km2
The Area of the Smaller River Basin PartFm14.86km2
Altitude of the First Contour Lineh0600m
EquidistanceΔh100m
The Lowest River Basin ElevationHmin535m
The Highest River Basin ElevationHmax1553m
Very Permeable Products From Rocks (Limestone, Gravel)fp0.12
Medium Permeable Rocks (Slates, Marls, Brownstone)fpp0.04
Poor Water Permeability (Heavy Clay, Compact Eruptive)fo0.84
The Volume of the Torrent Rainhb84.7mm
Average Annual Air Temperaturet08.9°C
Average Annual PrecipitationH yr873.7mm
Types of Soil Products and Related TypesY1.1
Numeral Equivalents of Clearly Exposed Erosion Processφ0.33
Table 6. The results of the IntErO model with the highlights on the impact of land use changes on soil erosion intensity in the Miocki Potok, Bijelo Polje Region, Polimlje, Montenegro.
Table 6. The results of the IntErO model with the highlights on the impact of land use changes on soil erosion intensity in the Miocki Potok, Bijelo Polje Region, Polimlje, Montenegro.
197019801990200020102020
A Part of the River Basin Under Forests [-]fs0.460.480.480.470.490.52
Grass, Meadows, Pastures and Orchards [-]ft0.390.400.400.380.370.38
Bare Land, Plow-Land, Without Grass [-]fg0.150.120.120.150.130.10
Coefficient of the Basin Planning [-]Xa0.490.490.500.500.470.45
Coefficient of the River Basin Form [-]A0.950.950.950.950.950.95
Coefficient of the Watershed Development [-]m0.280.280.280.280.280.28
Average River Basin Width [km]B3.693.693.693.693.693.69
(A)symmetry of the River Basin [-]a0.550.550.550.550.550.55
Density of the River Network of the Basin [-]G0.340.340.340.340.340.34
Coefficient of the River Basin Tortuousness [m]K1.051.051.051.051.051.05
Average River Basin Altitude [m]Hsr917.81917.81917.81917.81917.81917.81
Average Elevation Difference of the Basin [m]D382.81382.81382.81382.81382.81382.81
Average River Basin Decline [%]Isr31.7631.7631.7631.7631.7631.76
The Height of the Local Erosion Base of Basin [m]Hleb101810181018101810181018
Coefficient of Erosion Energy of Basin’s Relief [-]Er128.01128.01128.01128.01128.01128.01
Coefficient of the Region’s Permeability [-]S10.910.910.910.910.910.91
Coefficient of the Vegetation Cover [-]S20.740.730.730.740.730.72
The Water Retention in Inflow [m]W1.04881.04881.04881.04881.04881.0488
Energetic Potential of Flow on Torrent Rains [m km s]2 × gDF^½555.32555.32555.32555.32555.32555.32
Maximal Outflow From the River Basin [m3 s−1]Qmax372.2367.01366.35371.43366.82361.18
Coefficient of the River Basin Erosion [-]Z0.4830.4850.490.4910.4630.439
Production of Erosion Material in Basin [m3 yr−1]W yr376063789038492385633536532574
Coefficient of the Deposit Retention [-]Ru0.420.420.420.420.420.42
Real Soil Losses [m3 yr−1]G yr157931591216165161951485213680
Real Soil Losses Per [m3/km2 yr]G yr km−2384.66387.56393.72394.45361.34333.19
Table 7. Calculated soil losses in 57 river basins of Polimlje (m3 km−2 year−1).
Table 7. Calculated soil losses in 57 river basins of Polimlje (m3 km−2 year−1).
(1)645.40(12)278.33(23)324.45(34)325.19(45)212.39(56)269.25
(2)521.84(13)427.63(24)212.67(35)195.47(46)254.63(57)* 413.66
(3)104.31(14)330.12(25)385.41(36)264.43(47)200.43Average
331.78
m3 km−2 year−1
(4)288.97(15)429.10(26)492.68(37)286.07(48)514.60
(5)562.60(16)403.46(27)232.47(38)327.04(49)200.56
(6)399.52(17)370.61(28)305.76(39)452.92(50)247.93
(7)328.96(18)244.32(29)268.09(40)210.32(51)140.64
(8)180.22(19)219.39(30)266.21(41)131.23(52)315.28
(9)327.69(20)286.90(31)197.99(42)122.48(53)216.30
(10)298.19(21)88.66(32)296.45(43)194.76(54)250.39
(11)417.68(22)470.42(33)255.60(44)198.08(55)256.39
Source: [69,70,71,72,73]; * Data for the Studied Miocki Potok Basin [13].
Table 8. Data of Potpec reservoir filling and percentage of backfill since the beginning of measurements with respect to a free volume of the reservoir of 42.94 106 m3.
Table 8. Data of Potpec reservoir filling and percentage of backfill since the beginning of measurements with respect to a free volume of the reservoir of 42.94 106 m3.
Observation PeriodTotal YearsTotal Accumulation
106 m3
Total Sediments
106 m3
Percentage Backfilling
%
Annual Backfilling
106 m3
Annual Backfilling
%
1967–1976942.948.0018.630.892.07
1967–19811442.9411.0025.620.791.83
1967–19861942.9412.2028.410.641.50
1967–19912442.9411.6527.130.491.13
1967–19993242.9414.4833.720.451.05
1967–20063942.9415.3835.820.390.92
1967–20124542.9416.3738.120.360.85
2012–20175042.9416.8039.130.340.78
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Spalevic, V.; Barovic, G.; Vujacic, D.; Curovic, M.; Behzadfar, M.; Djurovic, N.; Dudic, B.; Billi, P. The Impact of Land Use Changes on Soil Erosion in the River Basin of Miocki Potok, Montenegro. Water 2020, 12, 2973. https://doi.org/10.3390/w12112973

AMA Style

Spalevic V, Barovic G, Vujacic D, Curovic M, Behzadfar M, Djurovic N, Dudic B, Billi P. The Impact of Land Use Changes on Soil Erosion in the River Basin of Miocki Potok, Montenegro. Water. 2020; 12(11):2973. https://doi.org/10.3390/w12112973

Chicago/Turabian Style

Spalevic, Velibor, Goran Barovic, Dusko Vujacic, Milic Curovic, Morteza Behzadfar, Nevenka Djurovic, Branislav Dudic, and Paolo Billi. 2020. "The Impact of Land Use Changes on Soil Erosion in the River Basin of Miocki Potok, Montenegro" Water 12, no. 11: 2973. https://doi.org/10.3390/w12112973

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