Geomorphology of the upper Mkhomazi River basin, KwaZulu-Natal, South Africa, with emphasis on late Pleistocene colluvial deposits

ABSTRACT We present a 1:50 000 scale geomorphological map of the upper Mkhomazi River basin, located in the foothills of the Drakensberg mountains in KwaZulu-Natal Province, South Africa. The sub-horizontal strata of the Permo-Triassic Beaufort Group forms plateau interfluves with a concave valley slope morphology. Locally, thick sequences of late Pleistocene colluvial deposits and associated buried paleosols (Masotcheni Formation) infill first-order tributary stream valleys and extend across the adjacent lower slopes. Surface runoff processes preferentially incise into the poorly consolidated, highly erodible sediments causing severe gully erosion that is responsible for widespread land degradation and desertification phenomena. The main purpose of this work is to derive a geomorphological map of the study area focussing on the erosional landforms to understand their spatial distribution and their relation to the colluvial deposits. Finally, a local and regional stratigraphic correlation of colluvial deposits and associated buried palaeosol profiles is proposed.


Introduction
In this work, we present a 1:50,000 scale geomorphological map of the upper Mkhomazi River basin in the central Drakensberg Mountains, Kwazulu-Natal, South Africa ( Figure 1). The area is characterized by a complex terrain morphology and is deeply incised by runoff-related soil erosion phenomena which dissect large areas of agricultural land.
The bedrock geology is composed of Permian sedimentary rocks of the Beaufort Group (Adelaide and Tarkastad Subgroups), which are sub-horizontally bedded, structurally influencing the typical deeply incised valley and high plateau terrain morphology. Other significant morphological controls are the Jurassic dolerite sills and dykes which crosscut the above-mentioned rock sequences. It was in this region that King and Fair (1944) and King (1949) associated sheet flow sediment transport and gully formation with long-term erosion, parallel slope retreat and pediment formation. However, they did not describe the detritus or transported waste that commonly overlies the cut bedrock pediment surface in this region. The association of episodic accumulation of poorly sorted sediment through the episodic accumulation of colluvial deposits on hillslopes in the Middleveld region of Swaziland was described by Price Williams et al. (1982) who noted the association of archaeological artefacts that was deposited during colluviation and the subsequent phases of soil development during hillslope stabilization. Colluvium is defined as poorly sorted sediment comprising an admixture of clay-, silt-, sand-and gravel sized particles that is eroded from upper slopes and deposited on lower hillslopes or pediments by sheet-wash processes (Watson et al., 1984). These hillslope cover deposits occur in a wide range of terrain morphological and bedrock contexts across parts of Southern Africa.
In the KwaZulu-Natal study area of South Africa these colluvial deposits are grouped as the Masotcheni Formation, a lithostratigraphic grouping of sheetwash and gully channel colluvial and buried palaeosols that represent several cycles of colluvium deposition, pedogenesis and erosion during the Quaternary (Botha, 1996;Botha et al., 1990Botha et al., , 1994Temme et al., 2008). The geological map of the Drakensberg foothills study area (Geological Survey, 1981) depicts only a simplistic, incomplete distribution of the colluvial deposits that delineates the sediments exposed by deep erosion gullies (dongas) rather than the broader deposits associated with infilled bedrock depressions that extend over wide areas between interfluves. In this study we present a more accurate and detailed assessment of the spatial distribution of the Masotcheni Formation in the upper Mkhomazi River basin based on field measurements, remote sensing, and geophysical surveys. The topographical range and lithological characteristics of the dispersive and highly erodible colluvium/palaeosol pediment cover deposits (i.e. fine clay and clay loam deposits on gentle slope (∼5°)), support rill-interrill erosion (Figure 2(a)) and gully erosion phenomena (Figure 2(b)). Gully erosion is the removal of soil by concentrated turbulent water flow through enlargement of rills and channels mainly due to intermittent stream flow after heavy rainfall events (e.g. Poesen, 1996;Sidorchuk et al., 2003). Gully erosion represents a geo-environmental problem and a serious land degradation issue in many parts of the world causing soil loss and sediment mobilization (e.g. Conforti et al., 2011Conforti et al., , 2014Kropáček et al., 2016;Pelacani et al., 2009;Sidorchuk, 1999;Sidorchuk et al., 2003;Zakerinejad et al., 2018). In their study of the physiochemical properties of the colluvial sediments and buried palaeosols typical of the Masotcheni Formation, Cama et al., 2020, Rienks et al. (2000 found that the vulnerability of the colluvium/palaeosols to gully erosion was related to the variable dispersivity and erodibility properties of the material, subsurface piping and turbulent surface runoff hydrology. The geomorphological mapdepicts four principal landforms; (i) fluvial and runoff associated landforms, (ii) slope landforms and deposits, (iii) litho-structural landforms, and (iv) anthropogenic structures. The landforms were mapped through field survey, orthophotos and Google Earth image analyses. Furthermore, gully forms and features have been characterized from a geological point of view using 20 stratigraphic columns of 14 well exposed and accessible gullies. A detailed terrain analysis based on a 12 m DEM allowed for a quantitative morphometric description of the gully systems. The map covers a regional area of about 500 km 2 and yields insights into the spatial distribution and lithostratigraphic nature of the colluvial deposits, highlighting the vulnerability of the area in terms of gully and rill-interrill erosion phenomena.
Moreover, the colluvial deposits of southern Africa are also relevant since the time elapsed between phases of colluviation permitted archaeological material to accumulate both during colluviation and the subsequent phases of stabilization and soil formation (Price Williams et al., 1982). The archaeological value of the Middle Stone Age (MSA) artefacts in the Masotcheni Formation is due to the occurrences in 'open sites' rather than the confines of cave deposits where MSA are usually excavated. The luminescence dated colluvial succession where Early-and Middle Stone Age stone artefacts were described by Botha (1996) also represents a much longer late Pleistocene temporal framework (∼130,000 years) than that preserved in most cave deposits. At the Sibudu Shelter near the coast in KwaZulu-Natal, Jacobs et al. (2008) correlated three periods of late Pleistocene hillslope colluviation in the region with periods when MSA inhabitants abandoned the shelter whereas the intervening periods bracketing palaeosol formation (Botha & Partridge, 2000) correlate with periods of occupation at Sibudu.
Finally, a detailed knowledge of the geomorphological processes and related forms allows for better land use management especially in the context of subsistence agriculture as practiced in the area today.

Geological setting
The study area (Figure 1) covers about 500 km 2 and is dominated by a mountainous landscape. The elevation ranges between 2267 m a.s.l. at the peak of the Drakensberg Mountains to ca. 1055 m a.s.l. at the Mkhomazi River. The region is characterized by a Subtropical Highland climate (Cwb), following the Köppen climate classification (Kottek et al., 2006), with significant precipitation of about 920 mm/year based on the precipitation time series from 1970 to 2005 (Nel & Sumner, 2008). From the base, the upper part of the Karoo Supergroup sedimentary succession exposed in the study area is represented by the Permo-Triassic Beaufort Group (Rutherford et al., 2015;Turner, 1981), specifically by the Adelaide Subgroup ( Figure  2(c)) and the upper Tarkastad Subgroup (Figure 2 (d)) (Catuneanu & Elango, 2001;Norman & Whitfield, 2006;Smith et al., 1993). The Karoo Supergroup succession is intruded by dolerite sills and dykes (Lyons et al., 2013) associated with the igneous event that produced the Drakensberg basaltic lavas. A long period of erosion and valley incision followed Gondwana breakup during the late Cretaceous and the study area is now covered by widespread, thin veneer of late Pleistocene colluvial deposits represented by the Masotcheni Formation. Pleistocene colluvial mantles occur in northern and central KwaZulu-Natal, the eastern Free State and some areas in central Swaziland, Botswana and Zimbabwe (Botha, 1996;Botha et al., 2016;Watson et al., 1984). Although similar colluvial palaeosol successions are widespread north and west of the Maluti mountains in the eastern Free State province and Lesotho, these occurrences are in a different 'geomorphic province' context (Partridge et al., 2010) to the Masotcheni Formation type area which extends across the Ladysmith Basin, South-eastern Coastal Hinterland and Queenstown Basin geomorphic provinces of central KwaZulu-Natal and the Eastern Cape province.
The Masotcheni Formation colluvial deposits are commonly concentrated where the hillslope morphology forces overland flow and sheetwash transported sediments to accumulate in bedrock depressions or colluvial hollows (e.g. Dietrich & Dorn, 1984;Mills, 1987;Reneau et al., 1989) and multiple episodes of gully cut-and-fill have been documented in this region (Botha, 1996).
The Masotcheni Formation includes fine clay, silty and sandy, poorly sorted, stratified colluvial sediments (Botha & Fedoroff, 1995;Lyons et al., 2013;Temme et al., 2008;Watson et al., 1984) (Figure 2(b)) generated by the erosion of weathered regolith and soils from upslope on the Drakensberg foothill interfluve ridges and deposited along their lower slopes (Botha et al., 2016). In many areas, the Masotcheni Formation is eroded by recently forming gullies, locally named 'dongas' (Lyons et al., 2013). Following the FAO soil Group Classification (IUSS Working Group WRB, 2014) the paleosols are classified as Solonetz and are characterized by silty-clay soils with a high-concentration of sodic clay in a Btn horizon. Moreover, some these paleosols are characterized by typical columnar, prismatic shaped peds with a polygonal structure on top of the Btn-horizon (Figure 2(b, e)). The colluvial deposits were deposited under semi-arid conditions (Watson et al., 1984), followed by long periods of hillslope stability represented by paleosol profiles (Botha, 1996). Concentration of runoff towards the colluviums-filled bedrock depressions on middle and lower hillslopes results in preferential erosion of the colluvium/paleosol succession which incorporates some dispersive and highly erodible sediment (Rienks et al., 2000).

Geomorphological setting
The steep, high relief topography is structurally controlled by the sub-horizontal strata of the Tarkastad and Adelaide Subgroups, and intrusive dolerite sills, leading to a classic plateau morphology (Figure 3). The dolerite intrusions mainly occur as sills and dykes that dominate the landscape due to their relative resistance to weathering and the related selective erosion processes. Severe weathering processes cause a typical 'onion skin' spheroidal weathering of dolerite deposits and a characteristic oxidation of soils giving a red colour to the latter (Figure 2(d)). The grassland and patches of sparse shrubby vegetation, combined with the lithological and topographic characteristics of the area as well as intense overland flow runoff processes, lead to widespread rill-interrill and gully processes and their related forms.
In the study area, the formation of gullies is also influenced by the presence of man-made contour terraces, where gullies can develop due to breaching of the contour banks by the concentration of surface runoff. Gully erosion exposes thick sequences of interbedded colluvial sediments and associated paleosols that define the stratigraphy of the Masotcheni Formation.
The Mkhomazi and Mkhomazana rivers represent the most important rivers of the area investigated. Long term erosion of the river system has preserved fluvial terraces at levels of up to 60 m above the current riverbed. These river terraces form flat surfaces that represent paleo-floodplains and valley bottoms. River terrace deposits are constituted by heterogeneous pebbles scattered in a fine red matrix. Shallow rotational and translational landslides and rockfall deposits occur on the steep slopes below the plateau margins ( Figure 2(f)).

Material and methods
The geomorphological forms and features were delineated, and the geomorphological map was generated following the methodology reported in Figure 4.

Geomorphological mapping and map assemblage
The geomorphological map was generated following several working steps: (i) detailed field survey in 2017/2018/2019, (ii) analysis of available orthophotos (2009-1: 10 000 scale) and Digital Elevation Model, and (iii) GIS based map compilation.
The field survey was conducted based on a 1:50,000 topographic map 2929DA Himeville (National Geo-Spatial Information, 2016). Once the principal landforms were mapped and features were checked, we validated the identified forms and features using a set of orthophotos provided by the Chief Directorate: National Geo-Spatial Information, 2009 as well as Google Earth images (© 2016 Google). The field mapping supported by aerial photo interpretation (API) was enhanced by a detailed terrain analysis using a DEM with 12 m resolution, from the TanDEM-X satellite platform provided by Deutsches Zentrum für Luftund Raumfahrt (DLR). The elevation model is an innovative spaceborne radar interferometer that is based on two TerraSAR-X radar satellites flying in close formation (Krieger et al., 2007). The terrain analysis was performed using SAGA GIS version 7.6.2 (Conrad et al., 2015). A set of different topographic indices (Table 1) were calculated following Wilson and Gallant (2000); Zakerinejad and Maerker (2014); and Bachofer et al. (2015). The topographic indices facilitate the geomorphological interpretation of landforms and features Maerker et al., 2019). Initially, the DEM was corrected to remove sinks using the algorithm proposed by Wang and Liu (2006) and subsequently we derived the indices represented in Table 1.
Furthermore, we integrated our field investigation with the existing gully inventory map provided by the Department of Environmental Affairs, Republic of South Africa, available at http://carbon-atlas.dirisa.org/download-data. The assembled database includes 279 gullies and 133 forms of sheet erosion that were digitized from Google Earth images (2016) and orthophotos (2009). Each landform mapped in the field was digitized using QGIS (version 3.2) with a unique symbol following the Italian guidelines for geomorphological cartography (after Campobasso et al., 2018). The symbols representing the processes, forms and features were differentiated into tectonic (brown), litho-structural (ochre), slope landforms and deposits (red), fluvial and runoff associated forms (green), and anthropogenic (black). The geological base map was derived from the lithostratigraphic subdivision of the existing 1:250 000 geological series map, 2928 Drakensberg (Geological Survey, 1981). The geological boundaries were digitized using QGIS and the colluvial deposits of the Masotcheni Formation partially mapped with DGPS and field surveys. The lithotypes present in the map were grouped on the basis of geomorphological characteristics and on their attitude in relation to degradationfeatures (Campobasso et al., 2018). The principal lithotypes are: (i) Sandstone and Quaternary deposits (yellow colour) represented by fluvial deposits   (FVL) and by the Masotcheni Formation (Qm), (ii) Intrusive rocks (bordeaux colour) represented by the dolerite dykes and sills (Jd) that intrude the bedrock formations, (iii) the Elliot Formation (grey colour) is characterized by red mudstone (Tre), and (iv) interstratified rocks consisting of alternating sandstone, shale and mudstone (light blue colour) representing the Tarkastad Subgroup (TRt) where sandstone and maroon coloured mudstone are the dominant lithotype. Finally, in blue colour the Adelaide Subgroup (Pa) in which the mudstone is the dominant lithotype. Additionally, three geological sections are depicted on the map to highlight the general stratigraphy and the main terrain morphology of the area.

Description of landforms and features
In the upper Mkhomazi River basin the principal geomorphological processes are related to the fluvial channel, as well as surface and subsurface erosion processes. Runoff from upslope and the sides of colluvial depressions is concentrated, causing severe erosion on the weakly consolidated colluvial deposits/paleosols that mantle the footslopes. Consequently, rill-interrill and gully erosion landforms and features are widespread in the area (see Main map). During the field investigation we conducted a sedimentological, pedological and stratigraphic description of the colluvial deposits/paleosols exposed by the gully erosion. In order to correlate the thickness and the stratigraphic distribution of the colluvial deposits and associated paleosols at a local and regional scale, we mapped and described profiles of 14 well exposed and accessible gully systems (see Main map). In each a detailed stratigraphic section was measured and the thickness, lithology, colour, and structure of the paleosols were described (see section 4.1). Furthermore, five morphometric indices derived from the DEM (elevation, slope, aspect, general curvature, and vertical distance to channel network) were statistically assessed and plotted using the software R (R Development Core Team, 2008). The indices were classified as follows; (i) the slopes ranging between 0°and 90°were split into six classes, (ii) aspect ranging between 0°to 360°was divided into eight classes, (iii) the general curvature which ranges between −0.5 and +0.5 was divided into three classes (−0.5 to −0.002 concave slopes, −0.002 to 0.0002 plan slopes, and 0.0002-0.5 convex slopes), (iv) the elevation ranging between 1000 m and 1600 m was split into six classes, and (v) the vertical distance to channel network which was divided into fourteen classes from 0 to 600 m. The other fluvial and runoff associated features occurring in fluvial channels as well as raised fluvial terrace remnants were mapped during the field survey with the support of Google Earth imagery and orthophoto maps.
River terrace deposits were identified and mapped on the field from 2017 to 2019 in the Mkhomazi and Mkhomazana basins (Figure 5(a)). A vertical distance to channel network was applied to extract the river terrace surfaces using a 5th order stream network. The QGIS plugin, QProf, was used to draw several river terrace profiles in order to correlate the river terrace deposits at the regional scale.
The study area is also characterized by dryland agricultural fields where contour terraces have been created to control and divert runoff from the ploughed lands. Poor maintenance of the contour banks often results in the development of rill-interrill and gully erosion ( Figure 5(b)). Fluvial and runoff landforms and features were mapped using a green colour. The gravitational mass movement landforms were symbolized in red polygons. The principal river systems, drainages and lakes were mapped from the orthophotos as polygonal shapes and are shown in blue.
The anthropogenic terrain forms are depicted using black lines. Quarries, landfill sites and terraced slopes were identified with vector polygons and the villages with a point symbol. Roads were mapped with white lines. Archeological areas and the locations of the stratigraphic columns were also reported in the map, with black point dots.

Geomorphological map
The landforms were identified, grouped, and described as follows: (i) Fluvial and runoff associated landforms, (ii) Slope landforms and deposit, (iii) Litho-structural landforms, and (iv) Anthropogenic forms.

Fluvial and runoff associated landforms
The fluvial landforms include river terraces, abandoned meanders (oxbow lakes), fluvial saddles, V-shaped valleys and crevasse channels.
The Mkhomazi and Mkhomazana rivers are characterized by incised meandering channels. During the field survey several raised river terrace deposits were identified. In order to extract the paleo-river terrace surfaces preserved in the Mkhomazi River basin we applied the 'Vertical Distance to Channel Network' module on a 12 m DEM. The identification of five terrace levels confirms the long-term incision by the river. The 'level 0' is the oldest terrace at a height of about 40 m above the present channel, the 'level 1' terrace is elevated 35-38 m, the 'level 2' terrace at 22-25 m, the 'level 3' terrace occurs at a height of about 18 m, and the 'level 4' terrace is the most recently incised, situated at a height of 5 m above the present day river bed. The profiles were extracted with QProf, where the elevation scale was exaggerated in order to better visualize the height variations.
Along the rivers the riverbanks are stable and vegetated. Furthermore, along the slopes, several V-shaped valleys were mapped. From the orthophoto analysis, flood crevasse channels and abandoned meanders were also identified in the area.
Finally, the resulting map depicts the rill-interrill and gully erosion features that were observed on both data sources (2016 Google Earth and 2009 orthophotos).
In the study area, gullies develop preferentially in the fine colluvial sediments and buried palaeosols of the Masotcheni Formation. The gully incision is up to 9 m deep and the bedrock and the overlying colluvial deposits can be readily distinguished. Gullies cover an area of 3.9 km 2 while rill-interrill erosion forms and sheet erosion features cover an area of just 1.3 km 2 , for a total area of 5.2 km 2 . Even if the area covered by the erosional features results small respect to the extent of the study area, these landforms manifest on the gentle slopes which are often used for agriculture purposes causing severe loss in agricultural and grazing land. Figure 6(a) shows a general and complete vertical section of a gully system (KwaThunzi gully) located in the Mkhomazi basin. In general, the Btn horizon is characterized by clay-loam soil and the A horizon is depicted by clay soil (Figure 6(b)). The burial of a sequence of paleosols is the result of successive periods of landscape stability (pedogenesis), followed by instability (erosion) and sediment accretion (gully infill) (Figure 6(c)). Btn 5 and A5 horizon were dated using Optically Stimulated Luminescence (OSL), a method used to estimate the burial ages of sediments (Adamiec et al., 2008). Preliminary OSL dating results (Figure 6(c)) indicate a the accumulation of this colluvium occurred at 66.62 ± 7.68 ka and 39.42 ± 5.27 ka which brackets the earliest colluvial deposition in the Okhombe valley which is in a similar geological and terrain morphological context towards the northwest (Temme et al., 2008). The Mkhomazi area dates do not include the earliest colluvial accretion and erosion events dated by Botha et al. (1994) and Botha (1996) in the Thukela basin to the northeast. Interestingly, these dates overlap a period during Oxygen Isotope Stage 3 when cooler and drier conditions prevailed at Sibudu shelter (Jacobs et al., 2008) when widespread erosion in the catchments made occupation of the region inhospitable to Middle Stone Age culture.
Btn5 horizon is clay-rich with thin beds of angular sandstone pebbles and does not show the typical columnar/prismatic structures of a solonetz soil. According to the Munsell soil chart the colour of this layer can be described as 10YR7/6. A5 horizon (2.5Y5/2) is rich in clayey components and includes sub-angular granule-sized siltstone clasts.
The following Btn4 is dominated by greenish clay (5YR8/2), silt and rounded sand grains. Btn3 horizon is dominated by clays (10YR8/1) and is characterized by fine, black material, which might be traces of organic matter. Btn 2 horizon and A1 horizon consist of brown clay and clay-loam, which includes spherical brown oxide particles. According to the Munsell colour charts, Btn 2 has a colour of 10YR7/2, A1 has a colour of 10YR6/1. The highest part of the succession is the current soil profile in which Btn1 horizon is characterized by a columnar structure and colour 7.5YR6/1. The topsoil is not present, having been eroded by surface runoff.
Twenty stratigraphic profiles were measured in 14 gully sidewalls in order to spatially correlate the Masotcheni paleosol horizons in the study area ( Figure  7). As highlighted in Figure 7 the most consolidated Btn5 horizon is always preserved and its thickness varies between 40 and 300 cm. Btn4 horizon ranges from brown or reddish colour depending on the influence of weathered dolerite components in the colluvial sediments. Other Btn and A layers are widely eroded in the area and are characterized by a range of thickness. This is due to the footslope inclination that leads to a differential erosion of the colluvial layers.
In order to characterize the gullies of the study area from a morphometric point of view we took into consideration the five morphometric indices derived from the 12 m DEM. The mean value from the point cloud of each gully was calculated and plotted in R where 88% are in an elevation range between 1200 m and 1400 m altitude (Figure 8(a)).

Slope landforms and deposits
The area steep upper hillslopes are presently stable and not affected by large landslides. However, localized, shallow landslides were identified in the field ( Figure  2(f)), caused by detachment of the upper soil cover, creating shallow depressions that concentrate runoff. The sub-horizontal attitude of the well-bedded Tarkastad Subgroup bedrock associated with steep slopes and high relief leads also to widespread rockfall debris below steep slope segments in the study area. Small rockfalls and topples occur in parts of the gully systems ( Figure 6(a)).

Litho-structural landforms
The litho-structural forms mapped include dykes, plateaus, structural scarps, ridges and saddles. The principal litho-structural landforms are directly related to the characteristics of the weathered Adelaide and Tarkastad Subgroups bedrock strata, i.e. alternating sandstone/mudrock beds, thickness of resistant strata as well as the sub-horizontal attitude of the bedding. Selective erosion acts on the friable clayey strata, resulting in undermining of the low scarps associated with the overlying, more resistant sandstone beds, forming the typical plateau topography.
The bedrock sedimentary succession is traversed by sub-vertical dolerite dykes which can be easily recognized in the field and from orthophoto analysis due to their positive weathering. The dolerite dykes define a preferential NW-SE aligned trend and are usually not covered by colluvial deposits.

Anthropogenic forms
In the study area the rural community is located in widespread homestead clusters and a few small villages. Dryland crop agriculture and cattle breeding are the principal sources of income in the area. Consequently, the principal evidence of human-induced landscape modifications is represented by agricultural contour bank terrace systems ( Figure 5(b)). These terraces in ploughed lands follow the contour and were established in the late 1950s (Botha, 1996) as a means of reducing rill erosion from crop lands. These structures are seldom adequately maintained and thus encourage soil erosion which starts with rill-interrill erosion and can develop into gully erosion features where contour banks are breached. Gully systems are often used as landfill or dumping area in the region. Finally, in the upper Mkhomazi basin there are small dams that have been constructed for irrigation purposes.

Conclusions
We present a 1:50 000 scale geomorphological map of the upper Mkhomazi basin, KwaZulu-Natal, South Africa. Thick sequences of fine colluvial deposits on gentle slopes have been preferentially eroded by surface runoff phenomena, yielding rill-interrill and gully erosion landforms and features. The colluvial deposits, previously only partly mapped and inadequately represented on geological maps in the study area (Geological Survey, 1981), were characterized and analysed in detail through an extensive field survey. A consistent sequence of buried, partially truncated paleosol horizons was identified. The Btn5-horizon is always preserved in the area, whereas the younger Btn and A horizons are more or less eroded. The high silt and clay component of the paleosols often includes sodic smectitic clay that imparts dispersive soil characteristics that enhance the highly erodible nature of the pedogenically altered colluvial hillslope cover deposits. High erosion rates and formation of rill and gully erosion reduce soil fertility and strip soil cover from agricultural fields. Soil erosion leads to the loss of agricultural productivity, environmental degradation and structural damage to buildings and road/bridge infrastructure. The geomorphological map represents an accurate inventory of erosional landforms, especially including gully and rill-interrill features and constitutes the basis for a quantitative assessment of erosion susceptibility in the upper Mkhomazi basin and can be used to gauge expansion of erosion features in the future.

Software
The maps were generated and realized using open source GIS (QGIS 3.0). The DEM was derived with SAGA GIS version 7.6.1. The stratigraphic columns were compiled with the free software Inkscape.