{'fraction_non_alphanumeric': 0.025683060109289616, 'fraction_numerical': 0.025075895567698847, 'mean_word_length': 4.779298245614035, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Effect of Gully Headcut Treatment on Sediment Load and Gully Expansion in the Sub Humid Ethiopian Highlands

The Ethiopian government has been implementing a land restoration program that aimed to restore degraded ecosystems and double agricultural productivity throughout the country since 2010. However, the success of the restoration program has been limited due to the lack of integrating gully erosion control measures. Consequently, many reservoirs in Ethiopia and downstream riparian countries have lost their storage capacity due to sedimentation, and studies demonstrated that gully erosion is one of the degradation hotspots within watersheds and contribute considerable proportion of the total sediment loads from a particular watershed. This study was conducted in one of large gullies in the Debre-Mawi watershed, northwestern Ethiopia to quantify the effect of gully head treatment in reducing the amount of sediment load generated from uplands and from the gully itself. We measured discharge, and sediment load and concentration in 2013 and 2014 at the upstream (inlet) and downstream ends (outlet) of the studied gully. Before the 2014 rainy phase, a gully headcut was stabilized with gabions at the bed and the gully bank was regarded to 45. The gully head retreated 12 m in 2013 but gully head retreat was stopped following the implementation of the treatment in 2014. The total sediment load and sediment concentration at the outlet was reduced by 42% and 30% respectively, in 2014 (i.e., after treatment) when compared to 2013 (i.e., before treatment). The result of this study support that controlling the upward retreat of gully head is effective in reducing sediment load and concentration as well as upward movement and expansion of gullies. However, maintenance of gully head control measures is the key to sustain the benefits.


Introduction
Soil and water conservation have been implemented in Ethiopian highlands since 1970s under the umbrella of the World Food Program 'Food for Work' [1]. However, this practice has only been sustainable as long as the program was actively engaged. One of the reasons was that farmers (the main stakeholders) were not involved in planning and evaluation of the benefit of soil and water conservation (SWC) measures and their indigenous experience and knowledge were rarely included [1]. This generally has resulted in failure of the SWC programs to bring desired results and a voluntary adoption of SWC technologies by farmers over the last four decades.
Gully erosion is one of the most damaging forms of soil erosion. In the Ethiopian highlands, gullies are particularly severe covering large tracts of areas and silting up rivers and reservoirs [2][3][4]. Extensive areas of agricultural lands are affected each year leading to irreversible changes in soil productivity and affecting the food security [5]. Most gullies in Ethiopia have been advancing uphill. Previous research confirmed that the upward head migration was the most important factor in gully expansion and sediment production [6,7].
Controlling gully erosion is more difficult and expensive than sheet and rill erosion. Once formed, it requires a huge investment to rehabilitate. It also requires knowledge of the actual gully development processes in the targeted area. This is because any gully rehabilitation practices that do not try to deal with the actual gully formation mechanics, including lack of maintenance of the installed erosion control structures are inefficient [7,8]. The Ethiopian landscape is a living example (Figure 1). Failed structures can be found in almost every gully where gully treatments were attempted. Trees planted on the bank immediately collapse in to gullies. For example, a diversion structure applied in the Debre Mawi watershed that aimed at rehabilitating formed a new gully (Figure 1c and d) following the diversion [7]. Thus, understanding patterns of gully and gully-complex development in Ethiopian highlands and its mitigation mechanisms, is of paramount importance not only from a scientific perspective, but especially to enable land managers to develop sustainable planning strategies. Moreover, quantitative information is rarely available on the effectiveness of gully headcut treatment techniques. The objective of this study is therefore, to quantify the effect of gully head protective measures on reducing sediment load, sediment concentration and headcut retreat rate. To do this, we selected a large gully in the Debre Mawi watershed which was previously studied by Tebebu et al. (2010).

Study Area Description
This study was conducted in the Debre Mawi watershed, located about 30 km south of Lake Tana at 11°20'13" N and 37°25'55 " E. The watershed drains an area of 608 ha, and its elevation ranges from 2186 to 2366 m above sea level. The slope in the watershed ranges from 1 to 2% for flat bottom lands and from 8 to 30% elsewhere. Rainfall is unimodal and the mean annual rainfall is about 1240 mm. More than 70% of the annual rainfall occurs during June, July, and August. The average monthly temperature is 20°C.
About 92% of the watershed is cultivated. The remaining portion of the watershed consists of communal grazing areas (6%), shrubs and eucalyptus woodlots (2%). The dominant soil types are Nitisol, Vertisols, and Regosol [9]. Nitisols dominate the upper parts of the watershed and is highly suitable for crop production. Vertisols are found in the lower portion of the watershed where crops such as tef (Eragrostis tef), chickpea (Cicer arietinum), and grass pea (Lathyrus sativus) are cultivated. Regosols occurs on the steep hillslopes and highly eroded parts of the watershed.
While the watershed is located at the most productive area of the country, the land is being devastated by expanding gully networks ( Figure 1) and to a lesser degree by sheet and rill erosion on the cultivated slopes [6,9,10]. The gully networks are mostly found in the periodically saturated bottomlands. One of the possible reasons for the saturation could be that lava dykes in the watershed affect the hydrology, forcing subsurface flow to the surface and causing saturated source areas for surface runoff upslope [3,7,11].

Measurements of Gully Expansion and Soil Loss
A 389 m gully ( Figure 2c) for which some information was available from an earlier study [3] was selected. The mean channel slope was 0.03. The gully was divided into 35 cross sections with the distance between two successive cross-sections varying between 6 and 13 m depending mainly on the uniformity of the gully width. This helps to minimize errors during calculating the average cross-sectional area between two cross sections.
Three gully topographic surveys were conducted: (1) before the rainy phase on 28 June 2013, (2) after the rainy phase on 18 Sep 2013 ( Figure 2), and (3) after the rainy phase on 18 Oct 2014 to quantify the extent of gully expansion and determine temporal changes in gully dimension. The gully cross-sectional geometry was surveyed by dividing the cross section into trapezoidal segments at abrupt changes in the ground profile, and measuring the width and depth of the gully at each segment. Cross-sectional area (A) and Volume (V) were then calculated: where n is the number of trapezoidal segment sides of height h and located a distance w from the left gully edge in a cross section (Figure 2c), i is a trapezoidal segment index, W is cross section width, j is cross section index, N is number of cross sections, and L j is length of the gully section between cross sections j and j+1. The total soil loss volume over the monitoring period was then obtained by the difference in V after and before the 2013 and 2014 rainy phases. To convert the volume of soil loss (obtained from gully topographic measurements) to mass, four soil bulk density (BD) samples were collected from four soil layers of the sidewalls of the gully 2 m downstream from the headcut with a 98 cm 3 cylindrical core sampler.
To evaluate the effects of headcut treatment on sediment load and concentration, two rectangular weirs, one at the upper end (inlet) and the other at the lower end (outlet) of the gully were installed prior to the rainy season of 2013. During rainfall-runoff events, stage height, and stream velocity were measured and sediment samples were collected at 10-minute intervals at both weirs from the beginning of the storm until base level was reached or flow stopped. Sediment samples were collected in one-liter bottles at two-thirds of the depth of the stream starting at the initial rise of the hydrograph until the initial depth was reached again. Stage height was measured manually by tape measure and the runoff velocity in the stream was measured by timing a float, put in at a distance 5 m upstream of the weir. The weirs followed the contour of the bank and did not restrict the flow, so that over or underestimation due to the cross-sectional change is minimal. In addition, during days without rainfall (or only base flow), discharge and sediment measurements were made and samples were collected three times daily (morning, noon, and evening) at the two weirs.
The discharge was determined in two different ways: (1) by multiplying the wetted cross-sectional area above the weir by two-thirds of the surface velocity; and (2) through rating curve for the occasions that the velocity is not available. Sediment concentrations were determined by filtering the samples over pre-weighed filter papers in the laboratory, which were then oven-dried for 24 hours at 105 °C and reweighed to the nearest 0.01 g. The suspended sediment load (Y) was then calculated using the equation Y = QC where Q is discharge and C is sediment concentration. Total amounts of runoff and suspended sediment per event were then obtained by integrating the individual measurements over time.

Headcut Treatment
We conducted three consecutive group discussions with the members of the local community and agricultural experts to identify a particular gully for the implementation gully head treatment (Figure 3a). In addition, the members of the local community were involved in the implementation and evaluation of gully head treatments. Accordingly we selected one 5 m high headcut and implemented the intervention with the support from the local community. The local community participated in collecting local materials such as stones woods and excavating equipment, while the researcher provided materials such as gabion and sacks.
The selected gully headcut and the sidewalls near the headcut was regraded to 45 o to decrease the slope and hence the bank failure (Figure 3b). Plant species [13] such as elephant grass (Pennisetum purpureum) and Sesbania susban were planted on the regraded banks (Figure 3b). In addition, to dissipate the flowing energy, one gabion check-dam underneath the headcut and others consecutive check-dams made of woods, at short distances below the headcut, were installed prior to the 2014 rainy season ( Figure  3b).

Estimation of Gully Erosion by Measuring the Changes in Gully Dimension
The total length of the gully between the two weirs was 389 m and the length from the headcut to its outlet in 2014 was 309 m with average depth 4.6 m and width 26 m. Its surface area was 9580 m 2 in 2013 and 9960 m 2 in 2014. The combined change in volume obtained from manual measurements of gully geometry between 2013 and 2014 was 4463 m 3 . The measured bulk density for this gully was 1.22 g cm -3 . Hence the total soil loss from the gully during the 2013-2014 was 5445 ton or 160 t ha -1 yr -1 .

Soil Loss as Affected by Gully Head Treatment
Runoff  Table 1). The daily sediment concentration and sediment load in 2014 were decreased by 33% and 40%, respectively (Table 1). During base flow, 7 t of soil entered the gully and 573 t was transported out of the gully (Table 1a). Thus, in the 2013 rain phase, the total sediment load entering the gully was 158 t and leaving the gully was 1,809 t. During the 2014 rain phase, after the gully headcut was arrested, the storm-generated sediment load at the inlet was a total load of 61 t and at the outlet was 764 t. During base flow, 14 t and 284 t of sediment loads were measured at the inlet and outlet, respectively. Thus, total sediment loads of 75 t at the inlet and 1,048 t at the outlet of the gully were measured during the 2014 monsoon season ( Table 1).
The gully enrichment ratio for storm-generated sediment load in 2013 was 7.2 and it was 11.5 in 2014 (Table 1b) which means, about 92% of the suspended sediment load originated from the gully ( Table 1, Figures 4 & 5). Including the base flow, over 90% of the sediment passing the outlet originated from the gully or the amount of sediment transported out of the gully was 10 to 14 times greater than at the inlet (Table 1b). In 2014, sediment concentration at the outlet of the gully slightly decreased towards the end of the rainy phase (similar to the inlet weir), whereas it increased during the 2013 rainy phase ( Figure 5).
The sediment concentration at the downstream end (outlet) of the gully is significantly greater than at the upstream end (inlet) as clearly viewed in Figure 4. The two pictures (inlet and outlet flow) in Figure 4 were taken at the same instance (August 6, 2014) and the outlet discharge color is dense brown indicating high sediment concentration whereas the inlet flow is light brown indicating a lower concentration. Quantitatively, this difference amounted to a high contribution to total sediment where more than 90% of the sediment originated from the gully itself (Table 1).

Effectiveness of Gully Head Treatment
The gully headcut treatment started prior to 2014 rain phase (May, 2014) as shown in Figure 3. This treatment was successful in controlling the uphill migration of the headcut as compared to other untreated gullies in the watershed. The sediment load and sediment concentration were respectively reduced by 40% and 30% in 2014 as compared to 2013. This is attributed to the implementation of gully headcut treatment in 2014 (Table 1, Figure 5). This was supported by the reduced upward migration of gully head following the implementation of the treatment. For example, in 2013, the headcut was linearly retreated by 12 m uphill whereas in 2014, the headcut was not retreated although there was a cross sectional growth or sidewall erosion (Figure 3c). Moreover, the rainfall in 2013 and 2014 was 917 and 1107 mm respectively, which is quite similar and this also support our argument that gully head treatment is effective in reducing sediment concentration and load.
However, the left sidewall of the headcut that was not regraded effectively was eroded by 1.5 m laterally ( Figure  3c). In addition, almost all the wooden check dams installed below the gabion check dam were totally damaged and carried away by a single high flooding event (90 mm). On the right side of the treated gully, the planted vegetation (elephant grass and sesbania) started slumping as shown in the last photo of Figure 3c. This shows that in the (sub) humid Ethiopia highlands, where the soil saturates periodically, halting gully erosion is more difficult and expensive. Thus, designing measures to treat deep gullies in cost effective ways in a humid monsoon highland should be a priority in future research in order to reduce sediment loads in rivers. These consist of lowering water table elevation and regrading of the gully head and sidewall effectively. We also suggest, treating gullies at the start of their initiation and development to reduce costs of reclamation and reclaim gullies effectively. Further, research is required to evaluate the effectiveness of integrated (physical and biological measures) gully head treatment in reducing upward migration of gully heads.

Conclusions
The results supported that gully headcut treatment is effective in reducing the upward migration of gully heads, and sediment concentration and load. We also observed that unless the gully control structures are maintained, they doomed to fail. This study is a good indicator of what areas need to be targeted to reduce sediment loads and thereby reduce damages on water resources infrastructures such as irrigation and hydropower dams. The result demonstrated that stabilizing big gullies is costly and difficult to manage. Thus, we suggest that gullies should be treated at the start of initiation and development.