Much of the research on (ephemeral) gully erosion comprises the determination of the geometry of these eroded channels, especially their width and depth. This is not a simple task due to uncertainty generated by the wide range of variability in gully cross section shapes found in the field. However, in the literature, this uncertainty is not recognized so that no criteria for their measurement are indicated. The aim of this work is to make researchers aware of the ambiguity that arises when characterizing the geometry of an ephemeral gully and similar eroded channels. In addition, a measurement protocol is proposed with the ultimate goal of pooling criteria in future works. It is suggested that the geometry of a gully could be characterized through its mean equivalent width and mean equivalent depth, which, together with its length, define an “equivalent prismatic gully” (EPG). The latter would facilitate the comparison between different gullies.
The classic forms of water erosion are caused by non-concentrated or laminar flow and concentrated flow; in the latter, a distinction is made between rill and gully erosion (Hutchinson and Pritchard, 1976). Rill erosion occurs in the form of numerous channels of a few centimetres in depth, distributed uniformly and randomly over sloping land (Soil Science Society of America, 2015) and easily obliterated by conventional tillage (Hutchinson and Pritchard, 1976). Furthermore, permanent gullies are distinguished from ephemeral ones (Foster, 1986; Thorne et al., 1986; Casalí et al., 1999). Permanent gullies are erosion channels which are too large to be eliminated by conventional tillage (Soil Science Society of America, 2015). Ephemeral gullies – present in agricultural soils – are, like rills, small enough for it to be possible to eliminate them by traditional tillage (Soil Science Society of America, 2015), hence their being qualified as ephemeral. However, when they form again, they tend to appear in the same places, contrary to what is observed in rills. This is explained by the fact that the ephemeral gullies are formed in the thalweg, which constitutes the confluence of two opposing slopes, a fact which conditions the trajectory of the runoff. Rills, however, occur entirely on one single slope (Casalí et al., 1999); their formation is, therefore, mainly subject to the high spatial variability of intrinsic factors of the soil (structural stability, hydraulic conductivity, etc.) and of its tillage.
The objectives of a large number of works on gully erosion have been the estimation of the spatial and/or temporal evolution of a gully or a network of gullies under different conditions (e.g. climate, land use) (e.g. Casalí et al., 2006; Gabet and Bookter, 2008; Campo-Bescós et al., 2013). For that purpose, as a first step, a morphological characterization of these channels is made. The most frequent way to do so is by measuring their width and depth and the ratio between both parameters (e.g. Giménez et al., 2009); their typology is also studied (for example, whether their cross section presents a particular shape, such as a U or a V). If the measurement of the length of the gully is added to this, it may be possible to determine its volume (eroded soil).
Examples of cross sections of typical ephemeral gullies (Navarre, Spain).
Consequently, for a precise description of the geometry of a gully, the correct determination of its width is a key factor. This is not always an easy task, especially when faced with cross sections with intricate shapes and diffuse limits. However, in the numerous scientific works on the subject, no uncertainty whatever is expressed regarding this measurement, and the criteria followed in the procedure are not specified. We believe that, as a general rule, it is assumed that a gully's width is defined by the imaginary line whose ends are located at the points of the two banks where an abrupt change in slope is manifested. This criterion is followed both in direct measurements in situ and in indirect ones taken from digital elevation models and mathematic algorithms ad hoc (e.g. Evans and Lindsay, 2010; Parker et al., 2012; Castillo et al., 2014). However, this procedure, at first sight reasonable and unquestionable, raises two objections. First, there is the presence of more than one point of slope inflection in one or both banks. Second, even when only one visible inflection point is presented on the slope of each bank – with the width of the channel thus being clearly defined – this also poses a question: do the limits of this channel, defined in this way, really correspond to the transversal limits of the erosive process which gave rise to the gully? Only by knowing the topography of the land in the moments before the formation of the gully can one answer that question with any certainty.
On the other hand, the width of a gully defines the upper limit of its cross section, therefore conditioning the subsequent determination of the depth of that channel. Furthermore, in this latter measurement (depth of the gully), another important ambiguity is added, i.e. the determination of the lower limit of the cross section (channel bed). This latter limit is usually located – in our belief – at the lowest point of the cross section, which is questionable in beds with a highly irregular cross-sectional profile. Nevertheless, the difficulty inherent in measuring a gully depth is not usually emphasized in the literature either.
In short, the lack of any protocol or universal criterion in determining the geometry of gullies would then cause a particular uncertainty when comparing the experimental results obtained by different researchers, for example, erosion rate values.
In this work it is sought to alert the scientific community to the uncertainties of which researchers are often unaware and which are triggered when characterizing the geometry of an ephemeral gully; for this purpose some examples of real cases will be shown. In addition, a measurement protocol is proposed, with the ultimate aim of pooling criteria in future works and experimentation. Although the criteria are proposed for ephemeral gullies, they would equally apply to similar erosion channels.
Uncertainty in the determination of a width in a cross section of a gully (real example). Arrows show different cross section widths and their corresponding depths (same colour). See Sect. 2 for details.
Illustration of the effect that the criterion followed to determine
the cross section width exerts on the computed volume of a gully reach.
Researchers, especially those new to the subject, when confronted with the measurement of gully geometry, assume that the limits of the erosion channel will present themselves in the field as being clearly defined, and, in fact, this is often true (see Fig. 1.1–1.3). However, on many occasions this is not the case (Fig. 1.4–1.6). It is therefore possible that a clear break in the slope of one of the banks (Fig. 1.6) or in both of them (Fig. 1.5) may not be noticed. Another possible ambiguity – independent of or additional to the previous one – is that which arises when both banks of the channel are uneven (Fig. 1.4, 1.6). This means that determining a single height value to trace an imaginary horizontal line between both banks is highly subjective. The length of this line should be understood as defining the width of the cross section being measured.
By contrast, when defining the depth of a gully, the lower limit of the cross section is usually well defined by the lowest point of the bed (see Fig. 1.2). However, what usually happens is that the location of this limit is also controversial, as can be seen in the cross sections in Fig. 1.1 and 1.3, where it is not clear whether this limit would really be represented by the lower height of the bed.
An incorrect determination of the width and/or depth of a certain gully may cause (important) errors in the determination of its volume, i.e. in the estimation of the eroded soil (Figs. 2 and 3). The magnitude of this potential experimental error would be less obvious, and even underestimated, if we analysed the cross sections individually (Fig. 2). However, an overall review of all the sections making up the gully being studied would give a better assessment of this measurement error. Figure 3 aims to illustrate the effect that the criterion followed to determine the cross section width exerts on the computed volume of a gully reach. A real gully reach was selected and three cross sections were used for calculating the volume of the reach (P1, P2 and P3) (Fig. 3a), the distance between cross sections being known. First, the eroded volume was calculated considering one possible criterion for defining the gully cross section width (in blue, Fig. 3b). Then the eroded soil was calculated again but considering a different possible criterion for defining the gully cross section widths (in red, Fig. 3b). The difference in the calculated volume is remarkable, increasing by 96 % from option b to option c. Figure 3 is just one example illustrating (i) the great differences in volumes that can be obtained in fixing the gully widths arbitrarily; (ii) the error that can be generated and (iii) the necessity of establishing rigorous and objective criteria and protocols. The purpose of Fig. 3 is similar to Fig. 2, the latter depicting the effect of the uncertainty in the determination of width in a single cross section of a gully.
Let us suppose that we have a detailed digital elevation model (DEM) of a
gully whose geometry we wish to determine (Fig. 4a). Similarly, we would also
have a DEM, not more than 1 year old, of the same area, but before the
gully in question had formed. Remember that the cycle of the formation
and obliteration of an ephemeral gully is conditioned by the periodicity
(usually 1 year) of the agricultural tillage responsible for it. We shall
call the DEM prior to the appearance of the gully DEM
Let us imagine now that, at any point
If we now carry out the subtraction of both DEMs but on their entire
surface, we should obtain the volume
Furthermore, knowing
Finally, the gully could be represented as a rectangular-based prism
(
In effect, we believe that the concept of equivalent prismatic gully shows
several benefits and applications. The principal one is that it
permits the determination of the most important characteristics of a complete
gully (
The width of a gully cross section, as defined in this paper, depends on the DEM pixel size and it depends on the type and size of the studied channel. Hengl (2006) concluded that, to prevent the loss of relevant information, the maximum pixel size must be the average of the minimum distances between sampling points. In the same way, Garbrecht and Martz (1994) fixed the pixel size to the size of the smallest distinguishable object. Additionally, the new methodologies available (terrestrial or aerial lidar, 3-D photo-reconstruction, etc.) provide very detailed information, which may be more than enough, in our opinion, for the purposes of these studies. However, these thresholds should be explored in future research.
In order to progress in gully erosion research, clear criteria to define and
determine the key morphological characteristics of gullies and their related
properties (such as volumes) are needed. In this paper, a new proposal for
advancing towards that goal has been submitted. Thus, starting from a precise
definition of the width of each gully cross section, the mean equivalent
gully width and depth are defined, as is the equivalent prismatic gully. This approach permits the determination of the most important
characteristics of a complete gully (
This study was partly funded by the Spanish Ministry of Science and Innovation (project CGL2011-24336). Edited by: E. Taguas