Next Article in Journal
Effect of Biochar Amendment in Woodchip Denitrifying Bioreactors for Nitrate and Phosphate Removal in Tile Drainage Flow
Next Article in Special Issue
Accumulation of Heavy Metal Ions from Urban Soil in Spontaneous Flora
Previous Article in Journal
Bioaugmentation-Enhanced Remediation of Crude Oil Polluted Water in Pilot-Scale Floating Treatment Wetlands
Previous Article in Special Issue
Development and Verification of the Available Number of Water Intake Days in Ungauged Local Water Source Using the SWAT Model and Flow Recession Curves
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Land Degradation and Soil Conservation Measures in the Moldavian Plateau, Eastern Romania: A Case Study from the Racova Catchment

1
Department of Geography, Faculty of Geography & Geology, Alexandru Ioan Cuza University of Iași, Carol I Blvd. 11, 700506 Iasi, Romania
2
Department of Life and Environmental Physics, Horia, Horia Hulubei National Institute for Physics and Nuclear Engineering, Reactorului Str. 30, P.O. BOX MG-6, 077125 Magurele, Romania
*
Author to whom correspondence should be addressed.
Water 2021, 13(20), 2877; https://doi.org/10.3390/w13202877
Submission received: 16 July 2021 / Revised: 5 October 2021 / Accepted: 6 October 2021 / Published: 14 October 2021
(This article belongs to the Special Issue River Basin Management and River Evolution Research)

Abstract

:
Land degradation by soil erosion, gullying and landslides and reservoir sedimentation is a major environmental threat in the Moldavian Plateau of eastern Romania. The widespread development of these processes in the last two centuries was favored mainly by traditional agriculture focused on ‘up-and-down slope’ farming on small plots. However, soil conservation measures were actively undertaken between 1970 and 1989. More recent legislation (No. 18/1991 Agricultural Real Estate Act) includes two provisions that discourage maintaining and extending soil conservation practices. Hence, the former contour farming system has been abandoned in favor of the traditional, inadequate farming methods. Thus, this paper reviews the impact of land degradation and soil conservation measures in a representative 32,908 ha catchment located in the Central Moldavian Plateau. Based on field measurements, the results show that the estimated mean long-term (1973–2017) sedimentation rate reaches 4.7 cm y−1 in the Puscasi Reservoir at the catchment outlet, resulting in an associated sediment delivery ratio of 0.28. The initial area of the Puscasi Reservoir at normal retention level has decreased by 32% and the water storage capacity has decreased by 39%. Consequently, land degradation remains a serious problem in the study area and effective soil conservation is urgently needed.

1. Introduction

Land degradation by soil erosion, gully erosion and mass movements is an important environmental threat throughout the world and poses major challenges to soil conservation measures. Over recent decades, significant progress has been made in understanding land degradation, its controlling factors and associated processes.
Soil erosion has been investigated in connection with its evolution over time, from past [1] to present [2] and even future scenarios [3]. The research scale varies largely, from planetary [4] to continental [5], regional [6], national [7] or local scales [8]. A lot of papers focus on quantitative assessment, using direct field measurements related to specific events [9], field plots [10] or modelling using GIS techniques [11]. Most frequently, research papers focus on specific issues regarding control factors such as climate change, land use, ecosystem pattern, [12,13] or the effects of erosion control measures [14,15,16]. Furthermore, a lot of research works developed in former communist countries from central-eastern Europe proves a direct impact of agriculture land use change and demographic change on sediment source and delivery over the last three post-communist decades [17,18,19,20].
Nevertheless, previous studies have focused on soil erosion in larger areas within the Moldavian Plateau, providing insufficient information at the smaller catchment scale. Many publications discuss environmental characteristics [21,22], whereas others focused on specific issues. Based on various sources, Motoc [23] estimated total erosion on agricultural lands in Romania, and suggested rates were especially high on the Barlad Plateau, with an average value of 20–30 t ha−1 y−1.
In a similar way, investigations of gully erosion include identifying magnitude (distribution and density) at different scales all over the world [24,25,26,27], controlling factors [28,29,30,31], estimating some rates of gullying [32,33,34], applying modern GIS techniques and statistical models [35,36,37], and gully control practices [38,39,40].
Two main areas of gullying were distinguished in the Moldavian Plateau [41,42]. In order to obtain a clear overview of the development of continuous gullies, 13 gullies were first sampled near the town of Barlad. Linear gully head advance, areal gully growth and sediment quantities were measured for three timeframes (1961–70, 1971–80 and 1981–90). The results indicate that gully erosion rates have decreased since 1960, due to changes in rainfall distribution patterns and the increased influence of soil conservation measures [43,44,45]. The mean denudation rate by landslides was estimated at 36.0 mm yr1 between 1968–1992 within the Barlad Plateau [46,47]. The large extent of landslides within three catchments (Upper Barlad upstream of Bacesti, Sacovat and Crasna) on the Central Moldavian Plateau was emphasized by Ionita et al. [48].
Reservoir sedimentation is a serious problem worldwide, with studies focusing on the specific sedimentation rates in small reservoirs [49], sedimentation rates and factors at different research scales [50,51], or the mitigation measures to be implemented [52]. Referring to the central and southern Moldavian Plateau, Ionita et al. [53] conducted a preliminary study using the caesium-137 (137Cs) technique to estimate sedimentation rates at two sites on the floor of the Puscasi and Pungesti-Garceana reservoirs. The 137Cs technique was also effectively used in areas of recent deposition of gully sediments, to provide chronological measures of gully development around the town of Barlad [54,55].
About three-quarters of degraded agricultural land on the Moldavian Plateau were adequately treated during the short-lived period (1970–1989) of proper soil conservation measures. Then, the traditional ‘up-and-down slope’ farming under small plots became prevalent again. In this context, the aim of the present study is to present both the current state of land degradation and soil conservation and their impact on the reservoir sedimentation in the representative Racova Catchment in the Central Moldavian Plateau of eastern Romania.

2. Study Area

Extending over ~27,000 km2, the Moldavian Plateau (MP) is the broadest and most typical plateau in Romania. Generally, it appears as an aggregate of plateaus, hills and rolling hills (collines) whose surface altitudes decrease towards the south–southeast.
The major units of the MP are the Suceava Plateau (SP), the Jijia Rolling Plain (JRP) in the northern part, and the Barlad Plateau (BP) and Covurlui High Plain (CHP) in the central–southern area. The BP is the most extensive high unit of the MP and covers > 8000 km2. The plateau comprises three major subunits: the Central Moldavian Plateau (CMP) in the north, the Tutova Rolling Hills (TRH) west of the Barlad Valley, and the Falciu Hills (FH) east of the Barlad Valley.
The Racova Catchment is located in the southern frame of the CMP at the border with the TRH (Figure 1). The Racova Catchment covers 32,908 ha, is 49 km long and is orientated in a west–east direction.
The Racova Catchment is made of clay-sandy and sandy-clay formations, almost exclusively deposited in deltaic facies. The sedimentary strata lie in a general homocline structure, with a gentle dip of 7–8 m km−1 towards the southeast [56,57].
The local relief is typically hilly, with altitudes varying from 485.5 m a.s.l. on the Mangalaria Hill in the southern watershed and 89 m a.s.l. on the Barlad floodplain.
Overall, the Racova valley shows a cross asymmetry. Generally, the left side of the Racova catchment covers 61% of the catchment, representing a broad cuesta back-slope. On the other side, the remaining 39% of the area represents a large north-facing cuesta.
The overall mean slope of the Racova Catchment is estimated to be 18.7%, which indicates that the study area has a high erosion potential (Figure 2) in most of the catchment (78.3%) [58].
The continental temperate climate is characterized by a mean annual temperature varying between 7.5 °C and 9.9 °C. The mean annual precipitation for the period 1963–2015 varied between 533 mm at 105 m a.s.l and ~700 mm in the area > 400 m a.s.l. Usually, two thirds of precipitation fall during the warm season, reaching a monthly maximum in June, with typical totals of 79–88 mm.
Bio-pedo-geographically, the higher areas are covered with deciduous forest. The sylvo-steppe is advancing in low areas. Accordingly, the zonal soils in the higher districts are Luvisols. Lower areas are predominantly Phaeozems and Chernozems. A large proportion of slopes is mantled by less productive soils such as Erodosols and Regosols, depending on the stage of land degradation processes.
Agricultural land covers 68% of the catchment area (Table A1). Some 35.4% is arable and 26% is pasture. The proportion of non-agricultural land is 32%, of which woodland covers 26.6%. This region was severely deforested during the 18th and especially the 19th century [59,60]. The area under the native forest in the Racova Catchment comprises only 17.8% of the area and has remained fairly constant since the late 19th Century. At present, woodland also includes the afforestation (silvic plantations) area, on 5.2% of the catchment area (1704 ha).

3. Methods

Several methods were deployed to estimate soil erosion losses, gully distribution, landslide inventory and reservoir sedimentation rates. Firstly, a digital elevation model (DEM) was created by digitizing the national 1:5000 topographic maps using TNT Mips version 7.1 software. The Revised Universal Soil Loss Equation (RUSLE) [61,62] based on the Universal Soil Loss Equation (USLE) was used to estimate long-term average annual soil loss from the Racova Catchment. The six factors [63,64], as adapted for Romanian conditions [65,66], are represented by: rainfall erosivity (R), soil erodibility (K), slope length & steepness (L&S), crop cover and management factors (C), and conservation practices (P).
Information for present-day land use was abstracted from the 2009 aerial orthophotos and 1:5000 topographic maps. These sources, together with the 2012 LiDAR images, have been successfully used to draw gully outlines and especially land covered by landslides. Useful information about land use and the gully network at the end of the 19th Century was drawn from both the topographic map of Moldavia (scale 1:20,000) and the 1894 Atlas of Moldavia (scale 1:50,000).
Levelling (topographic) surveys, using a Leica 407 TCR and GPS South 82V-Trimble, were conducted to obtain information about the behavior of the check-dams constructed during the early 1970s along gullies within the Ivanesti sub-catchment. The GPS South 82V-Trimble was used to obtain 2081 points on the former submerged floor of the Puscasi Reservoir, including the extensive but temporarily emerged area during spring 2017. Six bathymetrical cross-sections and a longitudinal profile were surveyed using the Midas Valeport Eco-sounder, type Bathy-500DF, in the permanently submerged area. A very detailed topographic map, consisting of a grid of 749 points covering the floor of the future Puscasi Reservoir, was completed in December 1969 by ISPIF (Institute for Land Treatments Studies and Designs), Bucharest [67]. We used the map to calculate the thickness and volume of deposited sediment in the reservoir over a 44-year period. A mean bulk density of 1.45 t m−3 has been frequently used in the study area to convert volumes to sediment yield (SY) and sediment delivery ratio (SDR) at the catchment outlet.
The 137Cs technique was used along gully floors in the Ivanesti sub-catchment to estimate the impact of soil conservation measures (check dams and afforestation). Gamma spectroscopy, associated with the Canberra MCA S100 system equipped with a Ge (Li) detector, was used to determine 137Cs concentrations in sediments. Soil surveys at a 1:10,000 scale were retrieved from OSPA Vaslui (Office for Pedological and Agrochemical Surveys) to distinguish the main soil classes and types [68]. Data processing was performed using Microsoft Office 2010 and particular attention was given to ‘ground truthing’ cartographic information.
In brief, the flowchart of the methodology starts from 1) the inventory of soil conservation measures carried out 44 years ago in a representative catchment located in the Central Moldavian Plateau and continues with 2) the assessment of their current state, being completed with 3) the analysis of the anti-erosional effect over time. By comparing the sedimentation volume within the Puscasi reservoir with the total gross erosion estimated for the entire tributary catchment, we try to demonstrate that land degradation remains a serious problem in the study area and extra soil conservation measures are urgently needed.

4. Results and Discussion

4.1. Land Degradation

Land degradation has been recognized as the major cause of environmental degradation worldwide [69,70] and, in particular, in the MP of eastern Romania. This area is highly susceptible to soil erosion [71], gullying [72] and landslides [73], which damage the local landscape by depleting soil resources and decreasing agricultural productivity.

4.1.1. Erosion

Soils in the study area are highly eroded (Figure 3A,B). A map of soil losses from agricultural land, with five erosion classes, was created by using the RUSLE as adapted to Romanian conditions (Figure 3C). The classes between 7–25 t ha−1 y−1 contributed an estimated 46% of the total, and one-third of soil loss was assigned to classes > 25 t ha−1 y−1. The mean estimated soil loss by water erosion (rill and inter-rill) on agricultural land was 21.6 t ha−1 y−1 and this area delivered an estimated 455 103 t y−1. Adjusting for the sediment contribution from woodlands, the mean specific sediment yield (SSY) decreased to 15.6 t ha−1 y−1. Finally, by adding gully erosion input, the SSY mean was 22.7 t ha−1 y−1 at the catchment scale. Therefore, rill and inter-rill erosion from agricultural land accounted for an estimated 61 or 69% of gross erosion and therefore supplies most sediment.
Gully erosion was much more limited on the CMP, due to more erosion-resistant substrata and forest cover, compared with other subunits on the BP. In the Racova Catchment, the present total gully length was 367 km, distributed in the tributary sub-catchments, and thus, gully density was 1.12 km km−2 (Figure 4A,B). These values are double those of the late 19th Century. However, at that time, the heads of the main gullies were located close to the watershed and today they enter as historical gullies [74,75]. That means former road gullies developed very rapidly after deforestation (Figure 4C). Although the area covered by gullies is small (2.7% of the total, 871 ha), they play important roles both in the triggering or reactivation of landslides and sediment detachment and transport. Using the relationships established by Ionita [44,45], based on medium-term field monitoring of some representative gullies on the BP, the SSY associated with gullying in the Racova Catchment was estimated at 7.1 t ha−1y−1 and accounted for (31.2) 31% of the sediment mass eroded by water. This relative contribution was half that of the Falciu Hills (FH), where gullies have incised and developed in a blanket of loess and loamy sands [55].

4.1.2. Landslides

Landslides are of particular concern on the Central Moldavian Plateau (CMP), both in terms of damage and affected areas. Despite the valuable research on landslides, there is a need to assess both their spatial distribution and the very recent temporal development of landslides (Figure 5A,B). The landslide inventory and map (Figure 5C) shows that landslides were highly variable in size, age, shape and form, and occurred on a total 56.2% (18,510.4 ha) of the catchment area. This was the highest identified proportion in the entire Moldavian Plateau (MP).
Based on sustained field campaigns, our landslide inventory shows that most landslides were stable (dormant) and the active ones formed only ~3% of the total landslide area (TLA), which is a representative value for the CMP. However, after the rainier 1973–1986 period, active landslides occupied 21.4% of the TLA [47]. The decline in landslide activity resulted from a decreasing tendency in precipitation totals since 1982. Most new landslides occurred by local reactivation of areas that had previously experienced landslide activity. Our conclusions are in accordance with those obtained by Pujină [47] based on field measurements or by Mărgărint et. al. [73] using modern GIS techniques.
The Landslide-Hypsometry Index (LHI) (i.e., the ratio of landslide area (LA) to total catchment area (CA) on hypsometric classes) showed a slightly asymmetrical distribution. The LHI peak value of 0.73 was typical of the 300–400 m contour interval. By cross-checking the landslide map with the slope map, it was evident that the three main slope classes (5–18, 18–27 and > 27%) amounted to one-third each of the LA. This finding underlines the even distribution of slope values within both sides of the Racova Catchment. Three-quarters of landslides developed on cuesta fronts with the remainder occurring on degraded cuesta back-slopes, especially in the upper sub-catchments. The deep-seated landslides were more frequently initiated in the Kersonian strata, because of weaknesses caused by the inter-bedding of sand and clay. These were one of the most characteristic landslide types, that is landslide amphitheatres, locally called “hartoape.”

4.2. Evolution of Soil Conservation Measures

By 1960, the traditional agricultural system on the hills of the MP consisted of ‘up-and-down-slope’ farming, with ~90% of agricultural land divided into small (< 1 ha) plots. Except in local areas, there was no concern about soil erosion and little awareness of conservation practices. After 1960, these areas were turned into co-operative farms. The remaining larger plots (~10% of the area) were transferred from farmers to State farms.
After several decades of quiescence, many new, innovative research studies on soil erosion control were initiated [66,76]. The first priority consisted of implementing one or more conservation practices, starting with contour ploughing. By late 1989, 75% (0.9 × 106 ha) of agricultural land at risk of erosion on the MP was adequately treated with conservation measures.
The new legislation (No. 18/1991 of the Agricultural Real Estate Act) includes two provisions that discourage soil conservation measures [55,77]. One of these stipulates that land reallocation must conform to the old locations; that is, plots must be orientated up-and-down slopes. The second refers to the successors’ land rights, which apply up to the fourth degree of kinship. Under these circumstances, land division increased and is now higher than before World War II. The major effect of the new law is the revival of the traditional agricultural system of ‘up-and-down slope’ farming. Another problem over recent decades is that the State ceased funding soil conservation, and investment in soil conservation has low priority among landowners.
Land changes within the Racova Catchment mirrored the general changes within the entire MP. During the 20-year period (1970–1989), much soil conservation work was accomplished, especially by IEELIF Vaslui (Enterprise for Performing and Exploiting the Land Improvement Works) [78], namely:
  • The design and construction of dams and reservoirs, such as Puscasi in 1973 (lake area 257 ha at normal retention level/NRL) and Trohan in 1982 (21 ha at NRL) on the Racova floodplain, and Pungesti-Garceni on the Garceneanca floodplain in 1976 (61 ha at NRL).
  • Design and construction of check-dams to control gully erosion in the tributaries of the River Racova.
  • Design and implementation of soil conservation practices on slopes in large farms, namely: strip-cropping, buffer strip cropping and especially bench terraces.
  • Design and building of drainage systems.
  • Filling small gullies, land reshaping using topsoil, and improving pastures.
  • Large-scale afforestation on 1704 ha on landslides and gullies, especially using black locust (Robinia pseudoacacia) and populus. Some 1001 ha of the afforested area were established by the Vaslui Silvic Enterprise and 703 ha by IEELIF Vaslui.
After implementing the provisions of Act No. 18/1991, the former short-lived contour farming system almost disappeared and the ‘up-and-down slope’ farming under small plots is prevalent again. The case of the 956 ha Ivanesti sub-catchment, located on the main Racova cuesta, illustrates this surprising evolution. Accordingly, a combination of contour strip cropping systems and bench terraces were implemented on 130.2 ha of arable land, especially northwest of the Ivanesti village. Here, on a field of 49.1 ha with an average slope of 13.6%, six 663–1024 m long bench terraces, spaced at 80–145 m intervals, were combined with six strip crops, each covering 4.3–11.5 ha (Figure 6).
Currently, the same field comprises 135 small individual plots, most of them orientated up-and-down slopes. About half of them are between two former bench terraces and their mean area is 0.37 ha. Others cross two, three, or four former strip crops and cover 0.6–1.6 ha each. Only nine plots are still on the contour, and these occupy 0.5–4.1 ha. Similar examples can be readily identified throughout the Racova Catchment (Figure 7A–F). Being one of the most important agricultural areas of Romania, during the 1990s and the first part of the 2000s, the Moldavian Plateau was characterized by the repopulation of rural areas, in the context of the collapse of the communist urban industry. This led to an increased pressure on agricultural land [79], having consequences in the further fragmentation of the agricultural land as well as in the continuous land degradation by erosion. In this part of Romania, compared to other central-eastern European regions [51], the most important land use changes related to the depopulation of rural areas and abandonment of agricultural land appear only in the second half of the 2000s, especially after 2007, the year of Romania’s accession to the EU.
In 1974, 14 check-dams (CDs) were constructed along valley bottom gullies, namely: six in the Canepa gully, five in the Balica gully and three in the Gologofta gully. They were set apart on reaches with similar original slopes of 3.0–3.2% but with various lengths, namely: 805 m in the Canepa, 811 m in the Balica and 310 m in the Gologofta gully. The bottom width (b) of their trapezoidal spillway varied between 2.8–6.6 m and the vertical drop (the difference between spillway crest and stilling basin/energy dissipator) ranged from 1.2–4.2 m Table A2. Simultaneously, afforestation on 15.1% of the total area (144.3 ha) has been deployed along the gully network and especially in areas with landslides.
Under these conditions, most streams are no longer competent to scour the bed or to undermine gully walls. The backwater effect of check-dams and the progressive impact of vegetation cover on stream flow decreased flow velocity and accelerate sediment deposition on the gully floor. In these cases, reducing bed gradient and the creation of trapezoidal cross-sections resulted in major changes in gully morphology (Figure 8).
Overall, gully depth decreased by 0.7–3.5 m and gully bottom width increased between 3.1–11.1 m relative to the original (restored) situation in 1974. Based on values abstracted from 18 gully cross-sections in the Ivanesti sub-catchment, it was possible to identify strong associations between gully depth and gully width for time periods, underlining the present-day lower values of gully depth (Figure 8A). Consequently, the lower width/depth (W/D) ratios, between 1.93–7.50 for the original situation, versus the higher values ranging from 3.39–10.39 at present, exemplifies gully infilling over 42 years (1974–2016). Figure 8B illustrates the strong correlation (r = 0.64, p < 0.01, n = 18) of the difference between the present and original W/D ratios versus the distance upstream of CDs Table A3. The difference between the original and present-day gully cross-section areas, ranging from 1.1–33.0 m2, emphasizes the growth of the infilled gully cross-section area. Figure 8C shows the trend of decreasing gully infilling with distance upstream of CDs in the Gologofta gully. This gully was particularly suitable for survey, as it had a uniform and symmetric cross-section, in which the gully banks were not disturbed by landslides and the cross-section was largely intact. Figure 8D summarizes changes of the original, design (desired) and present-day gradients in the same reach in the Gologofta gully.
Eleven CDs had a blanket of sediment in their stilling basins (SBs) that ranges from 0.5–1.8 m depth (Figure 9A). The recommended design (desired) gradient in accordance with the particle-size distribution (PSD) was 0.5% for fine sediment (clay, silt) and 1.0% for middle and coarse sand [77]. However, the sandy PSD within the Ivanesti sub-catchment, along with the mature vegetation cover, can maintain higher gradients, usually of 1.8–2.4%.
As a general conclusion related to the efficiency of these CDs, based on field measurements, all the data prove that the distance between those CDs was underestimated and their number overestimated. The SBs of two CDs were clean; consequently, those CDs were appropriately spaced. However, the SB of one CD was destroyed by stream scouring. The efficiency of these works is proven, but it could be obtained with a lower investment consisting of a smaller number of dams.
The progressively combined influence of conservation measures on sedimentation rates is emphasized by the depth distribution of 137Cs along the bottom of the Balica gully. Figure 9B,C shows the gully cross-section and the site of a 235 cm deep alluvial profile located 114 m upstream of CD1. The upper 220 cm consisted of recent sediments accumulated after 1974 as a result of the construction of the CD. The 137C peak value of 167 Bq kg−1, associated with the Chernobil accident of April 1986, occurred at 110–115 cm (Figure 9D). This indicates a mean sedimentation rate of 3.8 cm y−1 over a drier 30-year timespan (1986–2016). That rate was double between 1963 (peak year of nuclear weapons tests) and 1986, due to more precipitation [53,55]. Therefore, almost the entire column of sediment in the gully bottom was deposited after implementing conservation measures, especially during the early 1970s. Thus, the mean sedimentation rate over 42 years (1974–2016) was estimated to be 5.2 cm y−1.

4.3. Reservoir Sedimentation

Reservoir sedimentation data provides complementary information for evaluations of land degradation. As already mentioned, three dams were built in the study area, namely: Puscasi in the lower Racova Catchment, Pungesti-Garceni on the Garceneanca floodplain and Trohan in the upper Racova Catchment. They became operational in 1973, 1976, and 1982, respectively.
The Puscasi Reservoir is the most important and its initial area at the normal retention level (NRL) of 257 ha has decreased by 32.3% to its current area of 174 ha, while the water storage capacity has decreased by 38.6%, from 9.33 106 m3 to 5.73 106 m3. An accurate map of recent sediment deposition on the reservoir floor was created based on field measurements, using GPS (Figure 10). Deposition was assessed in eight classes and sediment thickness was uneven and adopted a deltaic shape.
The mean sediment thickness (STH) deposited in the Puscasi Reservoir was 206 cm and the mean sedimentation rate was 4.7 cm y−1. STH was greater in the upper-half of the reservoir and varied from 1.5–3.90 m, with a peak sedimentation rate of 9 cm y−1 over 44 years. A 137Cs profile in the area of high sedimentation showed a peak sedimentation rate of 11.5 cm y−1 during 1986–1998 [49]. The area exceeding 2.5 m STH highlighted the alluvial fan of the Racova River, which is the main contributor to siltation in the Puscasi Reservoir. However, the lateral input of close tributaries must be considered, since sediment discharge from upstream tributaries (such as the Tulbure, Ivanesti and Oprisita from the Racova cuesta front) also represents a major sediment source. An exception is the low alluvia input from densely forested sub-catchments, such as the Chelaru located upstream of the right shoulder of the dam, where forest covers 61% of the sub-catchment area.
The estimated volume of sediment within the Puscasi Reservoir was 5.3 × 106 m3, which represents a SY of 7.7 × 106 t. Estimated gross erosion (GE) of 24.8 t ha−1y−1 was associated with the 25,056 ha catchment area (without 5461 ha of the Pungesti-Garceni and Trohan reservoir catchment areas). Over 44 years, this totaled 27.3 106 t, and thus, the estimated SDR was 0.28. Trapping efficiency was assumed to be 100%, since the proportion of river sediment load captured by dams approaches this value in large reservoirs [80].
The smaller reservoirs (Pungesti-Garceni and Trohan) showed an even distribution of deposited sediment, as they are temporarily drained every 2–3 years to facilitate fish harvesting. Here, the mean sedimentation rate was ~2.7 cm y−1, as a consequence of the smaller catchment areas (3363 and 2098 ha, respectively) and the higher proportion of forest (32 and 44%, respectively).

5. Conclusions

The novelty of the current study resides in the application of an innovative and complex methodology using different field sampling and measurement methods applied in the field. Among others, the map of sediment thickness, the estimated volume of sediment within the Puscasi reservoir, as well as the estimated SDR of 0.28 are regarded as important issues in the present research. Other specific issues are:
(1)
The 32,098 ha Racova Catchment on the Moldavian Plateau of eastern Romania is highly susceptible to land degradation, due to both natural conditions and human impacts.
(2)
The mean value of soil losses by water erosion (rill and inter-rill) on agricultural land was estimated to be 21.6 t ha−1y−1. By adding the woodland contribution, this value decreased to 15.6 t ha−1y−1 and accounted for 68.7% of gross (total) erosion. The specific sediment yield by gully erosion was estimated at 7.1 t ha−1y−1 and averaged (31.2) 31.3% of the sediment mass eroded by water.
(3)
The estimated mean long-term (1973–2017) sedimentation rate reached 4.7 cm y−1 in the Puscasi Reservoir at the catchment outlet and the estimated associated sediment delivery ratio (SDR) was 0.28. The initial area of the Puscasi Reservoir at normal retention level (NRL) had decreased by 32% and the water storage capacity had decreased by 39% over 44 years.
(4)
Despite proper conservation measures that were designed and applied over a 20-year time-span from 1970–1990, after implementing the provisions of Act No. 18/1991, the contour farming system collapsed and returned to the traditional ‘up-and-down slope’ farming system on very small plots.
Furthermore, it is recommended that we intensify efforts to raise the awareness of citizens regarding the impact of land degradation and the societal importance of soil conservation. Farmers also need to revise their land use strategies, by focusing on better farming, improving pastures, and extending afforestation.
The need to conserve topsoil for future generations might be considered sufficient reason to justify establishing a Soil Conservation District covering the Racova Catchment. The recommended policy is that landowners and farmers must comply with soil control regulations provided by the District Commissioner. The Conservation District should be administered by a small board of directors, appointed by the County Council for five-year terms. The Board members should be local freeholders, residents of the Racova Catchment, and appointed on the basis of qualifications, without regard to political affiliation.

Author Contributions

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

Funding

This research was funded by the Department of Geography, Faculty of Geography & Geology, “Alexandru Ioan Cuza” University of Iași, Internal Research Program.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

Acknowledgment is given by L.N. to the infrastructure support from the Operational Program Competitiveness 2014–2020, Axis 1, under POC/448/1/1 research infrastructure projects for public R&D institutions/Sections F 2018, through the Research Center with Integrated Techniques for Atmospheric Aerosol Investigation in Romania (RECENT AIR) project, under grant agreement MySMIS no. 127324. We thank ABA Prut-Barlad (Prut-Barlad Catchment Administration) for the 2012 LiDAR images. We are also grateful to SGA Vaslui (Vaslui Water Management System) both for field access to the Puscasi Reservoir area and for the detailed topographic map made in December 1969 by ISPIF Bucharest, Romania. The assistance of Cristian Stoleriu and Alin Mihu Pintilie in performing bathymetric measurements is gratefully acknowledged.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Appendix A

Table A1. Land use in the Racova Catchment as abstracted from the 2009 aerial orthophotos.
Table A1. Land use in the Racova Catchment as abstracted from the 2009 aerial orthophotos.
CategoriesSub-CategoriesArea
ha% of Total
ArableTotal13,203.9540.13
Proper arable11,648.5835.40
Complex arable1555.374.73
Pastures 8588.9826.10
Vineyards 556.851.69
Orchards 47.330.14
Total agricultural land 22,397.1168.06
WoodlandTotal7431.5622.58
Forest5727.9817.40
Silvic plantations1703.585.18
Bushes & thorn bushes 1181.913.59
Water landLakes270.910.82
Roads 319.680.97
Land with buildingsTotal821.032.49
Yards. buildings713.652.17
Buildings & industrial yards105.800.32
Dams1.580.003
Wasteland 485.801.48
Total non-agricultural land 105,10.8931.94
General total 32,908.00100.00
Table A2. Check dams from the Ivanesti sub-catchment, Moldavian Plateau, Romania.
Table A2. Check dams from the Ivanesti sub-catchment, Moldavian Plateau, Romania.
GullyCheck Dam
(CD)
Location Upstream of Junction with (UJW) or Upstream of Check Dam (UCD)Trapezoidal SpillwayObs.
B
(m)
b
(m)
H
(m)
Useful Height (Hu) as the Difference between Spillway Crest and Stilling Basin
(m)
Sediment Thickness on Stilling Basin (SB)
(m)
CanepaCD1295.4 m UJW Ivanesti5.73.21.14.21.8
CD2182.5 m UCD16.63.81.22.21.3
CD3100.1 m UCD26.84.11.02.01
CD4241.6 m UCD37.84.81.552.951
CD5193.0 m UCD45.954.80.51.20.95
CD688.1 m UCD56.54.01.23.21.1
BalicaCD189.0 m UJW Gologofta 85.51.23.50.5
CD2277.3 m UCD14.03.20.752.10.7
CD321.2 m UCD26.14.21.21.5Destroyed SB
CD4228.9 m UCD35.452.81.02.8SB without baffle
piers
CD5283.8 m UCD45.43.01.12.71.30
GologoftaCD150.4 m UJW Balica5.902.951.33.2Clean SB
CD285.5 m UCD18.656.621.02.20.7
CD3225 m UCD26.724.251.01.71.5
Table A3. Some parameters of 18 gully cross-sections from the Ivanesti sub-catchment, Moldavian Plateau, Romania.
Table A3. Some parameters of 18 gully cross-sections from the Ivanesti sub-catchment, Moldavian Plateau, Romania.
GullyCross-Section
No.
Location
Upstream of
Check Dam
(UCD)
Ratio between the Present and Original Depth
(Dp/Do)
Difference between the Present and Original W/D Ratios
(Wp/Dp—Wo/Do)
Ratio between the Present and Original Area of Cross-Section
(CSp/CSo)
Gologofta1187 m UCD 20.621.350.80
2162 m UCD 20.591.690.71
398.6 m UCD 20.473.660.59
430.8 m UCD 20.572.760.70
545.6 m UCD 30.731.430.88
619.4 m UCD 30.641.870.80
Balica180.7 m UCD 10.722.880.77
27.6 m UCD 10.493.960.63
37.5 m UCD 20.464.710.63
4114.0 m UCD 50.751.160.86
Canepa18.8 m UCD 10.353.580.58
242.9 m UCD 20.593.570.75
310.6 m UCD 30.474.400.58
4230.4 m UCD 40.810.650.92
515.0 m UCD 40.621.760.81
681.3 m UCD 50.701.350.86
7165.6 m UCD 60.762.110.98
864.8 m UCD 60.542.640.68

References

  1. Dotterweich, M. The history of soil erosion and fluvial deposits in small catchments of central Europe: Deciphering the long-term interaction between humans and the environment—A review. Geomorphology 2008, 101, 192–208. [Google Scholar] [CrossRef]
  2. Poesen, J. Soil erosion in the Anthropocene: Research needs. Earth Surf. Process. Landf. 2018, 43, 64–84. [Google Scholar] [CrossRef]
  3. Panagos, P.; Ballabio, C.; Himics, M.; Scarpa, S.; Matthews, F.; Bogonos, M.; Poesen, J.; Borrelli, P. Projections of soil loss by water erosion in Europe by 2050. Environ. Sci. Policy 2021, 124, 380–392. [Google Scholar] [CrossRef]
  4. Garcia-Ruiz, J.M.; Begueria, S.; Nadal-Romero, E.; Gonzalez-Hidalgo, J.C.; Lana-Renault, N.; Sanjuan, Y. A meta-analysis of soil erosion rates across the world. Geomorphology 2015, 239, 160–173. [Google Scholar] [CrossRef] [Green Version]
  5. Panagos, P.; Katsoyiannis, A. Soil erosion modelling: The new challenges as the result of policy developments in Europe. Environ. Res. 2019, 172, 470–474. [Google Scholar] [CrossRef]
  6. Golosov, V.; Gusarov, A.; Litvin, L.; Yermolaev, O.; Chizhikova, N.; Safina, G.; Kiryukhina, Z. Evaluation of soil erosion rates in the southern half of the Russian Plain: Methodology and initial results. Proc. Int. Assoc. Hydrol. Sci. 2017, 375, 23. [Google Scholar] [CrossRef]
  7. Patriche, C.V. Quantitative assessment of rill and interrill soil erosion in Romania. Soil Use Manag. 2019, 35, 257–272. [Google Scholar] [CrossRef]
  8. Navas, A.; Machín, J.; Soto, J. Assessing soil erosion in a Pyrenean mountain catchment using GIS and fallout 137Cs. Agric. Ecosyst. Environ. 2005, 105, 493–506. [Google Scholar] [CrossRef] [Green Version]
  9. Niacșu, L.; Sfîcă, L.; Ursu, A.; Ichim, P.; Bobric, D.E.; Breabăn, I.G. Wind erosion on arable lands, associated with extreme blizzard conditions within the hilly area of Eastern Romania. Environ. Res. 2019, 169, 86–101. [Google Scholar] [CrossRef] [PubMed]
  10. Maetens, W.; Vanmaercke, M.; Poesen, J.; Jankauskas, B.; Jankauskiene, G.; Ionita, I. Effects of land use on annual runoff and soil loss in Europe and the Mediterranean: A meta-analysis of plot data. Prog. Phys. Geogr. Earth Environ. 2012, 36, 599–653. [Google Scholar] [CrossRef]
  11. Panagos, P.; Borrelli, P.; Poesen, J.; Ballabio, C.; Lugato, E.; Meusburger, K.; Montanarella, L.; Alewell, C. The new assessment of soil loss by water erosion in Europe. Environ. Sci. Policy 2015, 54, 438–447. [Google Scholar] [CrossRef]
  12. Borrelli, P.; Robinson, D.A.; Panagos, P.; Lugato, E.; Yang, J.E.; Alewell, C.; Ballabio, C. Land use and climate change impacts on global soil erosion by water (2015–2070). Proc. Natl. Acad. Sci. USA 2020, 117, 21994–22001. [Google Scholar] [CrossRef]
  13. Stefanidis, S.; Alexandridis, V.; Chatzichristaki, C.; Stefanidis, P. Assessing soil loss by water erosion in a typical mediterranean ecosystem of northern Greece under current and future rainfall erosivity. Water 2021, 13, 2002. [Google Scholar] [CrossRef]
  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. Maetens, W.; Poesen, J.; Vanmaercke, M. How effective are soil conservation techniques in reducing plot runoff and soil loss in Europe and the Mediterranean? Earth Sci. Rev. 2012, 115, 21–36. [Google Scholar] [CrossRef] [Green Version]
  16. Xiong, M.; Sun, R.; Chen, L. Effects of soil conservation techniques on water erosion control: A global analysis. Sci. Total Environ. 2018, 645, 753–760. [Google Scholar] [CrossRef]
  17. Gusarov, A.V. The impact of contemporary changes in climate and land use/cover on tendencies in water flow, suspended sediment yield and erosion intensity in the northeastern part of the Don River basin, SW European Russia. Environ. Res. 2019, 175, 468–488. [Google Scholar] [CrossRef]
  18. Tošić, R.; Dragićević, S.; Zlatić, M.; Todosijević, 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. Kijowska-Strugata, M.; Bucata-Hrabia, A.; Demczuk, P. Long-term impact of land use changes on soil erosion in an agricultural catchment (in the Western Polish Carpathians). Land Degrad. Dev. 2018, 29, 1871–1884. [Google Scholar] [CrossRef]
  20. Spalević, V.; Barović, G.; Vujacić, D.; Curović, M.; Behzadfar, M.; Ðurović, N.; Dudić, B.; Billi, P. The impact of land use changes on soil erosion in the river basin of Miocki Potok, Montenegro. Water 2020, 12, 2973. [Google Scholar] [CrossRef]
  21. Bacauanu, V.; Barbu, N.; Pantazica, M.; Ungureanu, A.L.; Chiriac, D. The Moldavian Plateau-Nature, Man, Society; Scientific and Encyclopedic Publishing House: Bucharest, Romania, 1980. [Google Scholar]
  22. Ungureanu, A. Geography of the Romanian Plateaus and Plains; “Al. I. Cuza” University Publishing House: Iasi, Romania, 1993. [Google Scholar]
  23. Motoc, M. Average Rate of Soil Erosional Degradation in R.S. Romania; Buletin Informativ ASAS, 2: Bucharest, Romania, 1983; pp. 67–73. [Google Scholar]
  24. Brooks, A.P.; Shellberg, J.G.; Knight, J.; Spencer, J. Alluvial gully erosion: An example from the Mitchell fluvial megafan, Queensland, Australia. Earth Surf. Process. Landf. 2009, 34, 1951–1969. [Google Scholar] [CrossRef] [Green Version]
  25. Guerra, A.J.T.; Bezerra, J.F.R.; Fullen, M.A.; Mendoça, J.K.S.; Sathler, R.; Lima, F.S.; Mendes, S.P.; Guerra, T.T. Urban Gullies in Sao Luis City, Maranhao State, Brazil, Progress in Gully Erosion Research. In Proceedings of the IVth International Symposium on Gully Erosion, Pamplona, Spain, 17–19 September 2007; Casali, J., Giménez, R., Eds.; Universidad Publica de Navarra: Pamplona, Spain, 2007; pp. 58–59. [Google Scholar]
  26. Rădoane, M.; Rădoane, N. Gully erosion. In Landform Dynamics and Evolution in Romania; Rădoane, M., Vespremeanu-Stroe, A., Eds.; Springer: Cham, Switzerland, 2016; pp. 371–396. [Google Scholar]
  27. Yermolaev, O.P. Erosion processes of the forest and forest-steppe zones in the eastern part of the Russian Plain. World Appl. Sci. J. 2014, 29, 453–459. [Google Scholar]
  28. Valentin, C.; Poesen, J.; Li, Y. Gully Erosion: Impacts, factors and control. Catena 2005, 63, 132–153. [Google Scholar] [CrossRef]
  29. Imwangana, F.M.; Vandecasteele, I.; Trefois, P.; Ozer, P.; Moeyersons, J. The origin and control of mega-gullies in Kinshasa (DR Congo). Catena 2015, 125, 38–49. [Google Scholar] [CrossRef]
  30. Wilkinson, S.N.; Kinsey-Henderson, A.E.; Hawdon, A.A.; Hairsine, P.B.; Bartley, R.; Baker, B. Grazing impacts on gully dynamics indicate approaches for gully erosion control in northeast Australia. Earth Surf. Process. Landf. 2018, 43, 1711–1725. [Google Scholar] [CrossRef]
  31. Xu, J.; Li, H.; Liu, X.; Hu, W.; Yang, Q.; Hao, Y.; Zhen, H.; Zhang, X. Gully erosion induced by snowmelt in Northeast China: A case study. Sustainability 2019, 11, 2088. [Google Scholar] [CrossRef] [Green Version]
  32. Vanmaercke, M.; Poesen, J.; Van Mele, B.; Demuzer, M.; Bruynseels, A.; Golosov, V.; Bezerra, J.F.R.; Bolysov, S.; Dvinskih, A.; Frankl, A.; et al. How fast do gully headcuts retreat? Earth Sci. Rev. 2016, 154, 336–355. [Google Scholar] [CrossRef]
  33. Nyssen, J.; Poesen, J.; Veyret-Picot, M.; Moeyersons, J.; Haile, M.; Deckers, J.; Dewit, J.; Naudts, K.; Govers, G. Assessment of gully erosion rates through interviews and measurements: A case study from northern Ethiopia. Earth Surf. Process. Landf. 2006, 31, 167–185. [Google Scholar] [CrossRef]
  34. Nachtergaele, J.; Poesen, J.; Oostwoud, D.W.; Vandekerckhove, L. Medium-term evolution of a gully developed in a loess-derived soil. Geomorphology 2002, 46, 223–239. [Google Scholar] [CrossRef]
  35. Gomez-Gutierrez, A.; Schnabel, S.; Berenguer-Sempere, F.; Lavado-Contador, F.; Rubio-Delgado, J. Using 3D photo-reconstruction methods to estimate gully headcut erosion. Catena 2014, 120, 91–101. [Google Scholar] [CrossRef]
  36. Niculiță, M.; Mărgărint, M.C.; Tarolli, P. Chapter 10-Using UAV and LiDAR data for gully geomorphic changes monitoring. In Developments in Earth Surface Processes; Tarolli, P., Mudd, S.M., Eds.; Elsevier: Amsterdam, The Netherlands, 2020; Volume 23, pp. 271–315. [Google Scholar]
  37. Poesen, J.; Torri, D.; Vanwalleghem, T. Gully erosion: Procedures to adopt when modelling soil erosion in landscapes affected by gullying. In Handbook of Erosion Modelling; Morgan, R., Nearing, M., Eds.; Blackwell Publishing Ltd.: Hoboken, NJ, USA, 2011; pp. 360–386. [Google Scholar]
  38. Liu, X.B.; Li, H.; Zhang, S.M.; Cruse, R.M.; Zhang, X.Y. Gully erosion control practices in Northeast China: A review. Sustainability 2019, 18, 5065. [Google Scholar] [CrossRef] [Green Version]
  39. Guyassa, E.; Frankl, A.; Zenebe, A.; Poesen, J.; Nyssen, J. Effects of check dams on runoff characteristics along gully reaches, the case of Northern Ethiopia. J. Hydrol. 2017, 545, 299–309. [Google Scholar] [CrossRef] [Green Version]
  40. Otlacan-Nedelcu, L. The uselfulness of a new model for the gully control structures effect prediction. In Sustaining the Global Farm; Stott, D.E., Mohtar, R.H., Steinhart, G.C., Eds.; Purdue University: West Lafayette, IN, USA, 2001; pp. 1000–1007. [Google Scholar]
  41. Ichim, I.; Mihaiu, G.; Surdeanu, V.; Radoane, M.; Radoane, N. Gully erosion in agricultural lands in Romania. In Soil Erosion on Agricultural Land; Boardman, J., Foster, I.D.L., Dearing, J.A., Eds.; Wiley: Hoboken, NJ, USA, 1990; pp. 55–68. [Google Scholar]
  42. Radoane, M.; Ichim, I.; Radoane, N. Gully distribution and development in Moldavia, Romania. Catena 1995, 24, 127–146. [Google Scholar] [CrossRef]
  43. Ionita, I. Geomorphological Study of Land Degradations in the Middle Barlad Catchment. Ph.D. Thesis, “Al. I. Cuza” University of Iasi, Iasi, Romania, 1998. [Google Scholar]
  44. Ionita, I. Gully Formation and Evolution in the Barlad Plateau; Corson Publishing House: Iasi, Romania, 2000. [Google Scholar]
  45. Ioniţă, I. Gully development in the Moldavian Plateau of Romania. Catena 2006, 68, 133–140. [Google Scholar] [CrossRef]
  46. Pujină, D. Research Concerning Landsliding Processes on Agricultural Lands in the Barlad Plateau and Contributions to Land Conservation Techniques. Ph.D. Thesis, “Gh. Asachi” Technical University of Iasi, Iasi, Romania, 1997. [Google Scholar]
  47. Pujină, D. Landslides in the Moldavian Plateau; Performantica Publishing House: Iasi, Romania, 2008. [Google Scholar]
  48. Ionita, I.; Chelaru, P.; Niacsu, L.; Butelca, D.; Andrei, A. Landslide distribution and their recent development within the Central Moldavian Plateau of Romania. Carpathian J. Earth Environ. Sci. 2014, 9, 241–252. [Google Scholar]
  49. Honek, D.; Michalková, M.Š.; Smetanová, A.; Sočuvka, V.; Velísková, Y.; Karásek, P.; Konecna, J.; Nemetova, Z.; Danáčová, M. Estimating sedimentation rates in small reservoirs-Suitable approaches for local municipalities in central Europe. J. Environ. Manag. 2020, 261, 109958. [Google Scholar] [CrossRef]
  50. Verstraeten, G.; Bazzoffi, P.; Lajczak, A.; Rãdoane, M.; Rey, F.; Poesen, J.; de Vente, J. Reservoir and pond sedimentation in Europe. In Soil Erosion in Europe; Boardman, J., Poesen, J., Eds.; Wiley: Hoboken, NJ, USA, 2006; pp. 759–774. [Google Scholar]
  51. Manojlović, S.; Sibinović, M.; Srejić, T.; Hadud, A.; Sabri, I. Agriculture land use change and demographic change in response to decline suspended sediment in Južna Morava River basin (Serbia). Sustainability 2021, 13, 3130. [Google Scholar] [CrossRef]
  52. Stefanidis, P.; Stefanidis, S. Reservoir sedimentation and mitigation measures. Lakes Reserv. Res. Manag. 2012, 17, 113–117. [Google Scholar] [CrossRef]
  53. Ionita, I.; Margineanu, R.M.; Hurjui, C. Assessment of the reservoir sedimentation rates from 137-Cs measurements in the Moldavian Plateau. Acta Geol. Hisp. 2000, 35, 357–367. [Google Scholar]
  54. Ionita, I.; Margineanu, R.M. Application of the 137-Cs for measuring soil erosion/deposition rates in Romania. Acta Geol. Hisp. 2000, 35, 311–319. [Google Scholar]
  55. Ionita, I.; Niacsu, L.; Petrovici, G.; Blebea-Apostu, A.M. Gully development in eastern Romania: A case study from Falciu Hills. Nat. Hazards 2015, 79 (Suppl. 1), 113–138. [Google Scholar] [CrossRef]
  56. Jeanrenand, P. Geology of Central Moldova between Siret and Prut. Ph.D. Thesis, “Al. I. Cuza” University of Iasi, Iasi, Romania, 1971. [Google Scholar]
  57. Ionesi, L. Geology of the Platform Units and the North-Dobrogea Orogen; Technical Publishing House: Bucharest, Romania, 1994. [Google Scholar]
  58. Samoila, C.; Ionita, I. Racova catchment, Geomorphological peculiarities. Lucr. Semin. Geogr. Dimitrie Cantemir 2017, 45, 109–124. [Google Scholar] [CrossRef] [Green Version]
  59. Poghirc, P. Village from Tutova Hills. Geographical Study; Scientific and Encyclopedic Publishing House: Bucharest, Romania, 1972. [Google Scholar]
  60. Stanga, I.C.; Niacsu, L. Using old maps and soil properties to reconstruct the forest spatial pattern in the late 18th century. Environ. Eng. Manag. J. 2016, 15, 1369–1378. [Google Scholar]
  61. Renard, K.G.; Foster, G.R.; Yoder, D.C.; McCool, D.K. RUSLE revisited: Status, questions, answers, and the future. J. Soil Water Conserv. 1994, 49, 213–220. [Google Scholar]
  62. Renard, K.G.; Foster, G.R.; Weesies, G.A.; McCool, D.K.; Yoder, D.C. Predicting Soil Erosion by Water: A Guide Toconservation Planning with the Revised Universal Soil Loss Equation (RUSLE). Agriculture Handbook, vol. 703; US Department of Agriculture: Washington, DC, USA, 1997; p. 384. [Google Scholar]
  63. Wischmeier, W.H.; Smith, D.D. A Universal Soil-Loss Equation to guide conservation farm planning. Trans. Int. Congr. Soil Sci. 1960, 1, 418–425. [Google Scholar]
  64. Wischmeier, W.H.; Smith, D.D. Predicting Rainfall Erosion Losses: A Guide to Conservation Planning. USDA Agriculture Handbook No. 537; Department of Agriculture, Science and Education Administration: Washington DC, USA, 1978. [Google Scholar]
  65. Motoc, M. Estimation de l’influence des facteurs d’érosion. In Proceedings of the International Water Erosion Symposium-Proceedings II, Prague, Czech Republic, 1970. [Google Scholar]
  66. Motoc, M.; Munteanu, S.; Baloiu, V.; Stanescu, P.; Mihaiu, G. Soil Erosion and Control Methods; Ceres Publishing House: Bucharest, Romania, 1975. [Google Scholar]
  67. Ionescu, A. Topographic Map of the Floor of Puscasi Reservoir; ISPIF (Institute for Land Treatments Studies and Designs): Bucharest, Romania, 1969. [Google Scholar]
  68. Pedological Study of Laza Territory. Includes Soils Map of Laza Territory-Scale 1:10,000; Office for Pedological and Agrochemical Surveys: Vaslui, Romania, 1997; p. 88. [Google Scholar]
  69. Eswaran, H.; Lal, R.; Reich, P.F. Land degradation: An overview. Response to Land Degradation. In Proceedings of the 2nd International Conference on Land Degradation (ICLD), Khon Kaen, Thailand, date of conference; Bridges, E.M., Hannam, I.D., Oldeman, L.R., DeVries, W.T.P., Scherr, S.J., Sombatpanit, S., Eds.; CRC Press: Boca Raton, FL, USA, 2001; pp. 20–35. [Google Scholar]
  70. Pravalie, R.; Patriche, C.; Borrelli, P.; Panagos, P.; Rosca, B.; Dumitrascu, M.; Nita, I.A.; Savulescu, I.; Birsan, M.V.; Bandoc, G. Arable lands under the pressure of multiple land degradation processes. A global perspective. Environ. Res. 2021, 194, 110697. [Google Scholar] [CrossRef]
  71. Ioniţă, I.; Rădoane, M.; Mircea, S. 1.13 Romania. In Soil Erosion in Europe; Boardman, J., Poesen, J., Eds.; Wiley: Hoboken, NJ, USA, 2006; pp. 155–166. [Google Scholar]
  72. Ioniță, I.; Niacșu, L.; Poesen, J.; Fullen, M.A. Controls on the development of continuous gullies: A 60 year monitoring study in the Moldavian Plateau of Romania. Earth Surf. Process. Landf. 2021, 46, 2746–2763. [Google Scholar] [CrossRef]
  73. Margarint, M.C.; Niculita, M. Landslide type and pattern in Moldavian Plateau, NE Romania. In Landform Dynamics and Evolution in Romania; Radoane, M., Vespremeanu-Stroe, A., Eds.; Springer: Cham, Switzerland, 2017; pp. 271–304. [Google Scholar]
  74. Poesen, J.; Nachtergaele, J.; Verstraeten, G.; Valentin, C. Gully erosion and environmental change: Importance and research needs. Catena 2003, 50, 91–133. [Google Scholar] [CrossRef]
  75. Poesen, J. Challenges in gully erosion research. Landf. Anal. 2011, 17, 5–9. [Google Scholar]
  76. Motoc, M.; Ionita, I.; Nistor, D.; Vatau, A. Soil erosion control in Romania. State of the Art. In U.S. Central and Eastern European Agro-Environmental Program, Soil Erosion Prevention and Remediation Workshop; Academy of Agricultural and Forestry Sciences: Bucharest, Romania, 1992; pp. 111–133. [Google Scholar]
  77. Gaspar, R.; Traci, C.; Apostol, A.; Necula, F.; Mesina, P. Normative for Design of the Correction Works of Torrents and Forest Improvement of Eroded Lands; Institute for Land Improvement Studies and Projects: Bucharest, Romania, 1967. [Google Scholar]
  78. Enterprise for Performing and Exploiting the Land Improvement Works Vaslui. In Soil Conservation Work Plans; IEELIF: Vaslui, Romania, 1970–1989.
  79. Rusu, A.; Ursu, A.; Stoleriu, C.C.; Groza, O.; Niacșu, L.; Sfîcă, L.; Minea, I.; Stoleriu, O.M. Structural changes in the Romanian economy reflected through Corine Land Cover datasets. Remote Sens. 2020, 12, 1323. [Google Scholar] [CrossRef] [Green Version]
  80. Verstraeten, G.; Poesen, J. Estimating trap efficiency of small reservoirs and ponds: Methods and implications for the assessment of sediment yield. Prog. Phys. Geogr. 2000, 24, 219–251. [Google Scholar] [CrossRef]
Figure 1. Location of the study area.
Figure 1. Location of the study area.
Water 13 02877 g001
Figure 2. Middle Racova catchment as seen next to an outcrop in lower Maeotian layers at 350 m on Poienesti Hill, towards Magura Hill (444 m) in the background.
Figure 2. Middle Racova catchment as seen next to an outcrop in lower Maeotian layers at 350 m on Poienesti Hill, towards Magura Hill (444 m) in the background.
Water 13 02877 g002
Figure 3. Excessive soil erosion under cambic Chernozems in the left side of Harsova valley on 4 April 2016 (A), severe erosion under Luvisols from upper Racovitza on 20 October 2013 (B), and map of annual soil losses on agricultural land in the Racova Catchment (C).
Figure 3. Excessive soil erosion under cambic Chernozems in the left side of Harsova valley on 4 April 2016 (A), severe erosion under Luvisols from upper Racovitza on 20 October 2013 (B), and map of annual soil losses on agricultural land in the Racova Catchment (C).
Water 13 02877 g003
Figure 4. Continuous gullies incised in the upper Poienesti sub-catchment on 20 October 2013 (A), new discontinuous incisions in an old V-shaped gully from Oprisita forest on 4 April 2016 (B), and map of gully distribution in the Racova Catchment (C).
Figure 4. Continuous gullies incised in the upper Poienesti sub-catchment on 20 October 2013 (A), new discontinuous incisions in an old V-shaped gully from Oprisita forest on 4 April 2016 (B), and map of gully distribution in the Racova Catchment (C).
Water 13 02877 g004
Figure 5. Reactivation of old landslides after deforestation in the upper catchment of the Racova upstream of Slobozia on 20 October 2013 (A), ample Canepa landslide amphitheatre (hartop), underlain mainly by Maeotian layers in the second step of the Racova cuesta front on 20 October 2013 (B), and landslide inventory map from the Racova Catchment, Moldavian Plateau, Romania (C).
Figure 5. Reactivation of old landslides after deforestation in the upper catchment of the Racova upstream of Slobozia on 20 October 2013 (A), ample Canepa landslide amphitheatre (hartop), underlain mainly by Maeotian layers in the second step of the Racova cuesta front on 20 October 2013 (B), and landslide inventory map from the Racova Catchment, Moldavian Plateau, Romania (C).
Water 13 02877 g005
Figure 6. Map of the current land use & land cover in the Ivanesti sub-catchment. Among the conservation measures implemented during 1970–1989, those presented on the map are still functional.
Figure 6. Map of the current land use & land cover in the Ivanesti sub-catchment. Among the conservation measures implemented during 1970–1989, those presented on the map are still functional.
Water 13 02877 g006
Figure 7. Shifting from contour conservation practices to up-and-down hill farming under small plots: (A) in the left side of the Ivanesti valley on 21 June 2017, (B) in the left side of the Oprisita valley on 4 April 2016, (C) in the middle Toporasti sub-catcment on 8 June 2016, (D) in the right side of the Poienesti valley on 4 April 2016, (E) in the Cosesti valley with remnants of former conservation practices on the right side (east looking cuesta back-slope) on 29 May 2017, and (F) in the left side (west facing cuesta front) of the Racovitza valley on 20 October 2013.
Figure 7. Shifting from contour conservation practices to up-and-down hill farming under small plots: (A) in the left side of the Ivanesti valley on 21 June 2017, (B) in the left side of the Oprisita valley on 4 April 2016, (C) in the middle Toporasti sub-catcment on 8 June 2016, (D) in the right side of the Poienesti valley on 4 April 2016, (E) in the Cosesti valley with remnants of former conservation practices on the right side (east looking cuesta back-slope) on 29 May 2017, and (F) in the left side (west facing cuesta front) of the Racovitza valley on 20 October 2013.
Water 13 02877 g007
Figure 8. Changes in gully morphology after implementing conservation measures: (A) relationship between GD and GW at present and in 1974, (B) relationship of the difference between present and original W/D ratios versus the distance upstream of check dams, (C) rate of gully infilling expressed by growth of the infilled gully cross section area (GIGCSA) to distance upstream of check dams (CDs) in the Gologofta gully, (D) longitudinal profile along a reach of 361 m from the Gologofta gully, upstream of the junction with the Balica gully.
Figure 8. Changes in gully morphology after implementing conservation measures: (A) relationship between GD and GW at present and in 1974, (B) relationship of the difference between present and original W/D ratios versus the distance upstream of check dams, (C) rate of gully infilling expressed by growth of the infilled gully cross section area (GIGCSA) to distance upstream of check dams (CDs) in the Gologofta gully, (D) longitudinal profile along a reach of 361 m from the Gologofta gully, upstream of the junction with the Balica gully.
Water 13 02877 g008
Figure 9. Gully infilling with recent alluvia after implementing conservation measures, (A) check dam no. 2 from the Gologofta gully on Nov. 29, 2016 (spillway bottom width = 6.6 m; vertical drop = 2.2 m from which 0.7 m is under sediments on stilling basin), (B) aggradation on the Balica gully floor mainly induced by afforestation on 23 November 2016, (C) cross section of the the Balica gully at the site of sediment profile, (D) depth distribution of 137Cs in the alluvia deposited on the Balica gully floor.
Figure 9. Gully infilling with recent alluvia after implementing conservation measures, (A) check dam no. 2 from the Gologofta gully on Nov. 29, 2016 (spillway bottom width = 6.6 m; vertical drop = 2.2 m from which 0.7 m is under sediments on stilling basin), (B) aggradation on the Balica gully floor mainly induced by afforestation on 23 November 2016, (C) cross section of the the Balica gully at the site of sediment profile, (D) depth distribution of 137Cs in the alluvia deposited on the Balica gully floor.
Water 13 02877 g009
Figure 10. Map of sediment thickness within the Puscasi reservoir from the lower Racova catchment, Moldavian Plateau, Romania.
Figure 10. Map of sediment thickness within the Puscasi reservoir from the lower Racova catchment, Moldavian Plateau, Romania.
Water 13 02877 g010
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Niacsu, L.; Ionita, I.; Samoila, C.; Grigoras, G.; Blebea-Apostu, A.M. Land Degradation and Soil Conservation Measures in the Moldavian Plateau, Eastern Romania: A Case Study from the Racova Catchment. Water 2021, 13, 2877. https://doi.org/10.3390/w13202877

AMA Style

Niacsu L, Ionita I, Samoila C, Grigoras G, Blebea-Apostu AM. Land Degradation and Soil Conservation Measures in the Moldavian Plateau, Eastern Romania: A Case Study from the Racova Catchment. Water. 2021; 13(20):2877. https://doi.org/10.3390/w13202877

Chicago/Turabian Style

Niacsu, Lilian, Ion Ionita, Claudia Samoila, Georgel Grigoras, and Ana Maria Blebea-Apostu. 2021. "Land Degradation and Soil Conservation Measures in the Moldavian Plateau, Eastern Romania: A Case Study from the Racova Catchment" Water 13, no. 20: 2877. https://doi.org/10.3390/w13202877

APA Style

Niacsu, L., Ionita, I., Samoila, C., Grigoras, G., & Blebea-Apostu, A. M. (2021). Land Degradation and Soil Conservation Measures in the Moldavian Plateau, Eastern Romania: A Case Study from the Racova Catchment. Water, 13(20), 2877. https://doi.org/10.3390/w13202877

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