SciELO - Scientific Electronic Library Online

 
vol.48 issue4Assessment of satellite-derived rainfall and its use in the ACRU agro-hydrological modelIrrigation scheduling for green beans grown in clay loam soil under a drip irrigation system author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Services on Demand

Article

Indicators

Related links

  • On index processCited by Google
  • On index processSimilars in Google

Share


Water SA

On-line version ISSN 1816-7950
Print version ISSN 0378-4738

Water SA vol.48 n.4 Pretoria Oct. 2022

http://dx.doi.org/10.17159/wsa/2022.v48.i4.3857 

RESEARCH PAPER

 

Calibration, validation and application of the SWAT model to determine the hydrological benefit of wetland rehabilitation in KwaZulu-Natal, South Africa

 

 

BC Scott-ShawI; R Lechmere-OertelII; TR HillIII

IDiscipline of Hydrology, School of Agricultural, Earth and Environmental Sciences, University of KwaZulu-Natal, Pietermaritzburg Campus, South Africa
IIBiodiversity Planning and Management, Box 1362, Hilton 3245, South Africa
IIIDiscipline of Geography, School of Agricultural, Earth and Environmental Sciences, University of KwaZulu-Natal, Pietermaritzburg Campus, South Africa

Correspondence

 

 


ABSTRACT

In South Africa, with highly variable and intense land-use practices, coupled with limited soil fertility and water resources, there has been a long history of encroachment of arable lands (sugarcane and timber plantations) nto surrounding wetlands. Although wetland delineation within the timber and sugar sectors is well-defined in policy, and existing and proposed legislation, there are significant areas of non-compliance. The spatially-explicit Soil Water Assessment Tool (SWAT) was adopted to investigate the interactions of climate, land-use and soil on the water-use of natural and encroached wetlands. This paper documents the calibration, validation and application of the SWAT model on Quaternary Catchment (QC) U20G, which is a 498 km2 catchment that forms part of the uMngeni River basin. The SWAT-CUP parameter sensitivity and optimization model was tested with daily observed streamflow data for this catchment. Parameters were modified using the sequential uncertainty fitting (SUFI-2) analysis routine to calibrate the model. The simulated flow had a close fit to the observed flow with a regression coefficient (r2) of 0.87 and a Nash-Sutcliffe (NS) coefficient of 0.8. Through the buffer scenario analysis, the model showed that if the wetland and a 20-m buffer were to be returned to a natural state, there could be a 16% increase in the annual streamflow contribution, with an upper limit of a 60% increase in some hydrologic response units (HRUs). Thus there would be a hydrological gain if wetlands and sensitive buffer areas were to be cleared of commercial timber species and sugarcane.

Keywords: SWAT, modelling, calibration, stream flow reduction, commercial plantations, sugarcane


 

 

INTRODUCTION

Commercial forestry has been shown to increase the green water (water lost by total evaporation) and decrease the 'blue water (water in rivers and dams) across South Africa (Jewitt, 2006; Everson et al, 2008; Gush, 2011). The National Water Act (Act No. 36 of 1998) declares commercial forestry to be a streamflow reduction activity (SFRA), whilst sugarcane, although not legislated as an SFRA, has shown to be a high water user (Bezuidenhout et al., 2006; Jarmain et al., 2014). All timber and most sugarcane in South Africa is rainfed or non-irrigated (Olivier and Singels, 2015), with soil water availability being the key limiting factor in their production. Historically, wetland areas, particularly in KwaZulu-Natal, have proved to be attractive for timber and sugarcane farmers, leading to the complete degradation of these wetland systems and, inter alia, an increased evaporation loss as a consequence of these 'thirsty' crops.

Although wetland delineation, within the timber and sugarcane sectors, has become well-defined in policy and existing (and proposed) legislation, and despite considerable pull-back in many timber estates, there are significant areas of non-compliance (Edwards and Roberts, 2006). The full implications of complete wetland delineation for production, regional economics, hydrological security and biodiversity have not been quantified and are poorly understood. Current drivers of policy and legislation are based on assumed, but untested, hydrological and biodiversity benefits. This study aims to calibrate the model selected to quantify the potential hydrological benefits of achieving complete delineation and buffering of wetlands within a significant forestry and sugarcane area of KwaZulu-Natal.

Hydrological modelling is an appropriate approach to investigate the complex interactions of climate, land-use and soil on the water-use of agricultural systems, in particular where spatial heterogeneity exists. KwaZulu-Natal has a spatially variable climate and complex land-use patterns, resulting in many of the commonly used models being unsuitable for spatial scenario testing. The objective of this study is to develop a methodology to apply a hydrological model to forestry and sugarcane delineation and buffering scenarios that could guide legislation, policy and management decisions. The Soil and Water Assessment Tool (SWAT) has emerged as one of the most widely used water quality, watershed and river basin-scale models, applied for a diverse range of hydrologic and/or environmental problems (Gassman et al., 2007; Gassman et al., 2014). The recent development of SWAT+ has allowed for the implementation of landscape units, improving flow and pollutant routing across the landscape (Bieger et al., 2017).

South Africa has been sub-divided into primary, secondary, tertiary and quaternary catchments (QCs). These interlinked and hydrologically cascading quaternary catchments vary in size from 48 to 18 100 km2 (Schulze and Horan, 2007). In this research, the calibrated baseline state, with the current patterns of land-use for QC U20G, was compared to two scenarios: (i) where all encroached wetlands were delineated and restored to a natural state, and (ii) an additional 20 m buffer surrounding the wetland was restored to natural vegetation. The model, using newly derived calibrated input values, was simulated for a time period beyond the calibration period. In addition, a highly afforested catchment with higher rainfall was simulated and compared with observed flow data. This catchment (QC U40A) was selected as it is one of a few QCs with good quality observed flow and no upstream catchments. For the calibration period, the model provided a close fit to the observed data, allowing for the model to be validated.

The SWAT model

In South Africa, models such as the Agricultural Catchments Research Unit (ACRU), Systéme Hydrologique Européen (SHE) model group and WAVES have been extensively applied (Lewarne, 2009). However, given the spatially complex nature of the data, the recent development of the ArcSWAT GIS interface, and recent studies adopting this model, SWAT (ArcSWAT ver. 2012.10_5.21) was selected as the most appropriate model. SWAT is a conceptual continuous time model developed in the early 1990s to assist in water resource management, to assess the impact of management and climate on water supplies and non-point source pollution in watersheds and large river basins (Arnold and Fohrer, 2005). Recently SWAT has been applied in tropical and sub-tropical regions of Africa (Easton et al., 2010; Schuol et al., 2008; Everson et al., 2018; Akoko et al, 2021), Asia (Thampi et al., 2010; Wagner et al., 2011), and Latin America (Strauch et al., 2013). It is physically based, uses readily available inputs and is computationally efficient to operate on large catchments over extended time periods (Everson et al., 2007). The required data inputs for SWAT to drive flows and direct sub-basin routing include; soil distributions, land-use patterns, management plans, elevation and daily/monthly climate (Arnold and Fohrer, 2005). SWAT integrates these parameters into hydrologic response units (HRUs), effectively replacing the underlying spatial data; these units are grouped according to topography, soils (type/structure/ depth/chemical properties), land use and slope. The SWAT model uses the water balance equation (Eq. 1) in its simulation of the hydrological cycle (Neitsch et al., 2011).

where SW( is the final soil water content (mm); SW0 the initial soil water content on day i (mm); Rday the precipitation on day i (mm); Qsurf the surface runoff on day i (mm); Ea the total evaporation on day i (mm); Wseep the water entering the vadose zone on day i (mm) and the return flow on day i (mm). The wetland module in SWAT is simple, where wetland area is allocated per sub-basin. The simulation uses unidirectional seepage from the wetland to groundwater, surface runoff and lateral subsurface generated from surrounding uplands and the wetland water surface area, with the discharge of water from the wetland to the river (Raham et al., 2016).

The study area

Quaternary Catchment U20G drains into the Mkhabela River, a key tributary in the uMngeni River that is managed by a state-owned entity, Umgeni Water (Fig. 1). U20G is heavily transformed primarily by commercial forestry and agriculture, with some built-up areas such as the town of Wartburg. Rainfall in the region occurs mostly in summer (December to February), with a mean annual precipitation of 863 mm, ranging from 650 mm towards the lowland southern areas to 1 060 mm in the upland northern areas. The reference potential evaporation (ET0) is approximately 1 675 mm (A-pan equivalent, after Schulze, 2011) and the mean annual evaporation is between 1 200 and 1 300 mm, which exceeds the annual rainfall, suggesting the catchment is water limited. Summers are warm to hot and winters are cool. The mean annual temperature ranges between 20.1°Cin summer and 11.7°C in winter. This 498 km2 catchment ranges from 1 120 m asl to 293 m asl at the catchment outlet. The underlying geology of the site is Natal Group Shale that leads to the formation of deep, well-drained soils such as the Avalon, Westleigh and Hutton forms.

The catchment area is dominated by commercial forestry (hybrids of Eucalyptus grandis, Pinus patula I P. elliotii, Acacia mearnsii) and sugarcane (Saccharum officinarum). The areas covered by these crops are extensive (Table 1) and there are numerous scientific findings indicating their high water use in contrast to the natural vegetation that they replaced (Olbrich et al., 1996; Everson et al., 2007; Dye et al., 2008; Gush and Dye, 2009; Gush and Dye, 2015). Many of these planted lands have historically, and in some instances currently, encroached into wetland areas or their buffers, sometimes in their entirety. This is largely due to the attractive arable conditions associated with wetlands: high fertility, easy access to groundwater, and amenable topography. Indeed, government policies and subsidies in past decades encouraged the conversion of wetlands to arable lands, particularly sugarcane and timber. The few remnant natural areas are a combination of KwaZulu-Natal Sandstone Sourveld, KwaZulu-Natal Hinterland Thornveld and Ngongoni Veld of Valley Thicket Biome and Natal Central Bushveld of the sub-escarpment savanna bioregion (Mucina and Rutherford, 2006).

 

METHODS

Elevation and topography

The output of SWAT is largely dependent on the resolution of the input data, in particular the digital elevation model (DEM). Much of the model development time was spent translating data into suitable input data. Figure 2 shows the elevation model, soils layer and land-use layer with the five climate stations indicated in red. These layers were used to derive a total of 3 745 hydrological response units (HRUs). This allowed for the spatial complexity of the model outputs to be maintained.

The 30 m Shuttle Radar Topography Mission (SRTM) 1 Arc-Second Global DEM was used to delineate the watershed (improved using verified point and contour data to correct inaccuracies associated with tall vegetation, and interpolated into a higher resolution DEM). Soils and hydropedological units were defined according to research published by Neitsch et al. (2011) and through personal communication (Van Tol and Lorentz, 2019). Daily data from five spatially variable weather stations were obtained and modified to drive the model.

Daily observed flow from the stream gauges at the two inlets (U2H014 below Albert Falls Dam and U2H012 on the Mpolweni River) and a sub-catchment outlet (U2H059 upstream of Nagle Dam) were used in the model. This allowed for the contributing catchment above U20G (U20A, U20B, U20C, U20D, U20E and U20F, totalling 2 605 km2), where less information was available, to be excluded within the model using accurate observations.

Land use and vegetation

A combination of provincial and industry-based GIS layers were integrated to create a new land-use layer. User-defined vegetation growth input parameters were constructed for the study area based on site observations, scientific publications, available databases (EKZNW, 2014) and expert opinion. The accuracy of the timber and sugarcane boundaries of the provincial landcover layer (EKZNW, 2014) was improved by cutting in the relatively accurate commercial forestry and sugarcane layers obtained from the various industry GIS databases. The accuracy of the wetland delineations was improved by cutting in a detailed wetland layer digitised from aerial imagery (Lechmere-Oertel, 2017), including a 20 m buffer around the delineated wetland edge.

The wetland and 20 m buffer layers retained the detail of the underlying landcover classes, essentially creating a series of new classes, such as 'natural wetland, 'timber in wetland', and 'sugarcane in buffer. The original wetland areas were still simulated with wetland characteristics but with a modified land use. Any areas that would not realistically change (e.g. built up areas and roads) were left unchanged within the wetland and buffer areas. The SWAT Land-use Definition Tool (ArcSWAT interface) was used to clip the land-use layer to the catchment boundary and to reclassify the land-use layer to match the attributes contained in the SWAT database. The resulting landcover layer indicates the significance of timber and sugarcane in U20G, and shows the degree to which they encroach into wetlands and their 20 m buffer (Table 1).

Specific SWAT input parameters were obtained for the key vegetation types in this area (Scott-Shaw, 2018, Table 2). The Eucalyptus, wattle and sugarcane parameters were modified to match species and hybrids grown in KwaZulu-Natal. New parameters for 'Wattle in Wetlands', 'Eucalyptus in Wetlands' and 'Sugarcane in Wetlands' were generated. From a land management perspective, it is important to determine which areas should be recognised as priority areas for management changes (i.e. changes in future planting areas). These are areas that provide the greatest hydrological benefit due to their current level of encroachment and location in the catchment. To identify these areas, three scenarios were applied:

1. Current/baseline state defined by the latest landcover intersected with detailed timber, sugar and wetlands outlines (including probable wetlands encroached by commercial forestry and sugarcane).

2. Baseline landcover with all encroached or invaded wetlands rehabilitated to a natural state (either hygrophilous/sedge dominated or riparian forest without management).

3. Baseline landcover with all encroached or invaded wetlands rehabilitated to a natural state and a 20 m buffer applied to the wetland edge. The buffer was delineated as a natural grassland (excluding natural forest areas).

Soil physical and chemical properties

No SWAT soil database currently exists in South Africa. As such, soils information was taken from cxistingsites within the catchment boundary (described by Le Roux et al., 2015 - Hydrology of South African Soils and Hillslopes) and extrapolated, where possible, using terrain models. Some soil attribute data were obtained from ongoing research projects in the region. The structure, depth, number of layers and texture were used to construct a detailed soil layer with up to five variable soil horizons. Where available, the South African Soil Classification system (Soil Classification Working Group, 1991) and the soil hydrologic group (NRCS, 1996) were used to determine the soil form and family and to estimate the key physical and chemical soil properties for the typical soil forms required by SWAT (Table 3). Additional soils data were used to create a more comprehensive soil database, and are not included in Table 3. The use of the Soil-Landscape Estimation and Evaluation Program (SLEEP) was considered but not utilised, as a suitable soil layer was manually derived.

Climate

SWAT weather data definitions were modified to include data from five climate stations in U20G, all of which had daily rainfall records. However, only two had daily temperature, relative humidity, solar radialion and wind speed records. The climate period (from July 1971 to December 2018) accounted for wet, dry and average climate years. The averaged monthly data can be seen in Fig. 3.

Management

Within a land-use type, management is crucial for hydrological simulations as it has a significant impact on the hydrological partitions of the lands. SWAT allows for the adoption of management operations through potential heat units (PHU) or through fixed dates of specific operations. The management operations were modified in SWAT at an HRU level to specify the initial growing state and periods during harvest, fallow lands and planting. Local management practices were accounted for where information could be obtained. Depending on the vegetation type, a management period was applied allowing for the complete removal of the vegetation at specified intervals. A 10-year rotation was applied for Eucalyptus and Acacia species, while a 15-year rotation was applied for pulp Pinus species. The planting operation was distributed according to satellite observations following periods where significant areas were clearfelled. However, fixed 'kill' and 'plant' dates do not adequately represent evergreen species as the LAI decline is linear towards zero, rather than sigmoidal towards the LAIMIN (Strauch and Volk, 2013). Although considered in this study, a process-based soil moisture approach (modified plant growth module - Strauch and Volk, 2013) was not used but is recommended for future studies.

Within the management routine of the model, for rainfed sugarcane in KwaZulu-Natal, the land use was ratooned up to a maximum of 6 times, taking approximately 10 to 12 years. During the replanting period, tillage was initiated in the form of a deep disk and rip. Further information can be obtained from Abdalla etal. (2019).

Model calibration and validation

It is commonly accepted that deterministic calibration approaches are outdated and fundamentally flawed (Abbaspour et al., 2015). On the contrary, stochastic calibration approaches recognize the errors and uncertainties in our models and attempt to capture, to some degree, the lack of understanding of the processes in natural systems (Abbaspour et al., 2015). However, as SWAT input parameters are process-based, they must be held within a realistic uncertainty range during the calibration process (Arnold et al., 2015). Automatic calibration may lead to substantial errors if limitations exist in measured data, there is a lack of knowledge on physical processes and operational procedures, and there is uncertainty on mathematical equations and model sensitivity (Moriasi et al., 2007). Thus stricter performance ratings should be adopted during model calibration than during validation, with the inclusion of uncertainty analysis (Moriasi et al., 2007).

In sequential uncertainty fitting (SUFI-2), uncertainty in parameters, expressed as ranges (uniform distributions), accounts for all sources of uncertainty, such as in driving input variables (e.g., rainfall), conceptual model, parameters, and measured data (Abbaspour et al., 2015). SWAT-CUP was run using observed streamflow data extracted for the outlet of the catchment using verifieddata obtainedfrom the Department ofWater andSanitation. A total of 300 simulations were performed with the relative parameters activated. Model sensitivity is defined as the change in model output per unit change in parameter input (Tesfahunegn et al., 2012). After pre-processing of the required input for the SWAT model, monthly flow simulations were performed for 46 years of recording periods (1972 to 2018). Three years were used as a 'warm-up' period, following which the simulation was used for a sensitivity analysis of hydrologic parameters and calibration of the model. The sensitivity analysis was performed using the standalone SWAT-CUP sensitivity analysis tool that uses the Latin hypercube one-factor-at-a-time sampling (LH-OAT).

A validation period for a higher rainfall catchment with intensive commercial forestry was selected. The model (using parameterized inputs) was used to simulate QC U40A (Mistley catchment) over a 31-year period (1985 to 2018). This catchment was selected as it is relatively close to QC U20G, has good quality observed streamflow data and is an isolated catchment, limiting the uncertainty present where upstream catchments exist.

 

RESULTS AND DISCUSSION

Sensitivity analysis

Results from the sensitivity analysis using SWAT-CUP allowed for the mean relative sensitivity of each input parameter to be determined and ranked (Table 4). The ranking was determined from the resultant P-value (significance of the sensitivity where the parameter becomes significant if the P-values are close to zero). The SCS runoff curve number (CN2) and maximum canopy storage (CANMX) were the most sensitive parameters. The authors took consideration where fitted values were produced which were outside of realistic known bounds to ensure that unrealistic input values were not used. The fitted value indicates the best value of the objective function for the simulation with the best fit. A minimum and maximum value range was used when parameterizing the model, where relevant.

Model calibration

Although the pre-calibration results provided a reasonable fit with the observed data, the model generally over-simulated streamflow Possible causes include: inaccurate spatially distributed climate data, under-estimation of abstraction for irrigation or, more likely, incorrect input parameters (particularly for lesser studied land types). The post-calibration simulated monthly flow had a close fit to the observed flow (Fig. 4) with a regression coefficient (r2) of 0.87 and a Nash-Sutcliffe (NS) coefficient of 0.8. The model uncertainty was predicted using the computed percent prediction uncertainty (95 PPU) and dotty plots for each parameter (plots of parameter values or relative changes versus objective function).

Model validation

Observed streamflow was obtained from U4H002 at the Mistley station, representing QC U40A. The model simulated monthly volumes well with a r2 of 0.67, Nash-Sutcliffe of 0.8 and a root mean square error (RMSE) of 0.70 (Fig. 5). The p-factor (percentage of observations covered by the 95 PPU) was 67% and the r-factor (average thickness of the 95PPU band divided by the standard deviation of the measured data) was 0.59, indicating an acceptable goodness of fit.

Low flows were well simulated in this catchment. However, some peak events were over-simulated. Through further investigation, it is likely that the large wetland situated above the outlet of the catchment attenuates these peak flows in reality and is not accommodated for in this model simulation. Although wetlands were included as a land use, the detailed functioning of this wetland was not included, which would release water when the volume exceeds the storage of the wetland. Thus, if the wetland were not at full capacity, flows would be attenuated, much like it would in a pond or reservoir. However, it does illustrate that the output volumes from the model are realistic and could be applied to similar catchments within KwaZulu-Natal. The inclusion of significant wetland areas is recommended for future modelling studies.

Model application

The calibrated model was simulated for a larger catchment area comprising of seven QCs. The results show that if wetlands and a 20 m buffer were to be returned to a natural state, as much as a 16% increase in streamflow could be gained (with an upper limit of 60% streamflow contribution per HRU). This is a significant increase and highlights key priority areas that not only are heavily encroached by commercial forestry or sugarcane, but are important water generating areas.

The water yield output is defined as the total amount of water leaving the HRU and contributing to the main channel for the given time step. The model clearly shows that water yield is not uniform across the study area (Fig. 6) but is rather correlated to bio-climatic variations. Of widest applicability is the relationship between land-use and water yield, with greatest yields over areas of natural grassland. The impact of land-use can be clearly seen throughout the upper uMngeni (particularly in QC U20F) where the plantation and sugarcane areas have a low water yield compared to the surrounding natural land.

Evapotranspiration (ET) is an important component as it comprises almost half of the water balance in most catchments in the area and is the primary determinant of streamflow reduction. ET, throughout the study area, is highly variable and dependent on the landcover input (Fig. 6). Water bodies, wetlands, plantations and sugarcane use the highest amount of water. The ET is further increased in wetlands and riparian areas where encroachment by timber and sugarcane has occurred. In addition to ET, commercial forestry and sugarcane can significantly impact upon the groundwater component. This is accounted for in SWAT through the 'revap' coefficient, where water may move from the shallow aquifer into the overlying unsaturated zone, allowing deep-rooted plants to be able to take up water directly (Neitsch et al., 2011). Areas of low recharge can correspond to soils with higher clay contents and vegetation with a high biomass and deep rooting systems. Areas of higher rainfall have a greater recharge to the shallow and deep aquifers whereas areas of low rainfall, coupled with geology and soils that do not promote infiltration, have a low recharge (Fig. 6).

The results highlight the spatial variability in response to clearing scenarios, at scales ranging from farm units to quaternary catchments, when compared to a historic baseline time period. Clearly there are priority areas where the co sts of removing timber and/or sugarcane and rehabilitating the natural vegetation are more likely to have a greater return on investment. The importance of these spatial outputs is the relative difference between scenarios which could indicate key areas for clearing at a very fine scale. There is also potential for climate scenarios to be considered, ranging from simple temperature and rainfall corrections to more complex climate models and vegetation parameter modifications.

Priority areas were identified by converting the water yield output from each scenario to a high-resolution raster grid (5 x 5 m) calculating the difference in water yield between the different scenarios. The distribution of priority areas (where the % increase between scenarios is > 50%) is spatially highly uneven, suggesting that there are distinct areas where it would be beneficial to enact the scenarios (Fig. 7).

 

CONCLUSION

The overall objective of the study was to calibrate the SWAT model for a well-researched and data-rich catchment in KwaZulu-Natal, with the aim of comparing hydrological gains due to the removal of exotic commercial species from wetlands and surrounding buffer areas. Previous modelling studies in this area have not captured the spatial complexity that exists in these catchments, where there are significant hydrological differences throughout each catchment. In addition, studies have not considered the existing high water-use of wetland and riparian systems. This study has modelled, calibrated and validated a highly complex catchment system using the SWAT model. The SWAT-CUP parameter sensitivity and optimization model (Abbaspour, 2015) was tested with monthly observed streamflow data. Parameters were modified using the SUFI-2 analysis routine to calibrate the model. This study has allowed for the calibration and validation of the SWAT model and the subsequent application of the model to a broader study area comprising of 17 QCs. The model calibration showed that with slight adjustment to sensitive input parameters, while keeping within realistic bounds, a strong correlation could be obtained between simulated and observed streamflow. The simulation was within the upper- and lower percent prediction uncertainty bounds. However, 49% of the observations were within the uncertainty band (p-factor of 0.49) and the r-factor (average thickness of the uncertainty band divided by standard deviation of observed data) was 0.36. This suggests that the model uncertainty is within an acceptable range, although there is potential for improvement by increasing the number of iterations and including additional observations of additional parameters. Thus one potential direction for future research is to extend the scale (catchment level), undertake similar modelling exercises for neighbouring and/or paired catchments and increase diversity of land-use types and management practices, to increase the regional-specific variables and diversity. This will provide us with greater confidence in model output at the regional scale and begin to allow comparison between biophysical characteristics and land-use practices between catchments and crops.

The spatial distribution of total evaporation (ET), water yield and groundwater recharge provide a valuable output for scenario testing. This study went beyond determining an area-weighted increase due to clearing, as it allowed for a soil-, slope- and climate-specific simulation for each land unit. The changes in hydrological partitions were output based on landcover, climate, soils and the proximity of the landcover to wetlands. This allowed for a highly detailed output of the spatial distribution of the hydrological partitions.

The results show that the SWAT model is a suitable model to be used for a range of land uses and soil types in KwaZulu-Natal. The calibration, which allowed for 200 behavioural simulations, resulted in the percentage of data bracketed by the 95% prediction uncertainty (95PPU) being 0.49 (p-factor of 0.49 and r-factor of 0.36). The regression coefficient (r2) was 0.87 with a Nash-Sutcliffe (NS) coefficient of 0.8, indicating a good fit and confidence in the input parameters used. The model, using newly derived calibrated input values, was run for a time period beyond the calibration period. In addition, a nearby catchment was simulated and compared with observed flow data. Much like the calibration period, the model provided a close fit to the observed data, allowing for the model to be validated. The model, through a scenario comparison, allowed for the compilation of annual water balances, time-series data and spatially explicit data that were linked back to the derived CIS data to provide a useful spatial priority-benefit product. Some of the cleared hydrological response units had an increase of streamflow contribution exceeding 60%. Within the greater catchment area (including upstream cleared catchments), an average annual increase of 18.2 million cubic meters was observed under the cleared scenario. This was as a result of 152 km2 of commercial forestry and 26.8 km2 of sugarcane being cleared from the wetland and the surrounding 20 m buffer.

The management component in SWAT is detailed and suitable for catchments in KwaZulu-Natal that are afforested and cultivated with sugarcane. If this model were to be calibrated for additional spatially variable catchment areas, the model could be adopted for a broad range of applications in South Africa. However, a major limitation is the lack of SWAT-ready input data, such as soils, land use and climate. Although climate change projections have not been considered in this study, the methods displayed, and the subsequent results, provide a suitable platform for the application and testing of likely climatic extremes.

 

ACKNOWLEDGEMENTS

The authors wish to acknowledge the financial contribution of the SANBI-GEF 5 'Biodiversity and Landuse project' (UNDP Grant). The Water Research Commission (Project No. K5/2402: Assessing the impact of erosion and sediment yield from farming and forestry systems in selected catchment of South Africa) and Fountain Hill Estate, where much of the observed data were obtained and analysed, are thanked for their time and funding.

 

REFERENCES

ABBASPOUR KC (2015) SWAT-CUP: SWAT Calibration and Uncertainty Programs - A User Manual. Department of Systems Analysis, Integrated Assessment and Modelling (SIAM), Eawag - Swiss Federal Institute of Aquatic Science and Technology, Duebendorf. 100 pp.         [ Links ]

ABDALLA K, DICKEY M, HILL TH and SCOTT-SHAW BC (2019) Assessment of soil erosion under rainfed sugarcane in KwaZulu-Natal, South Africa. Nat. Resour. Forum, https://doi.org/10.1111/1477-8947.12179        [ Links ]

AKOKO G, LE TU, GOMI T and KATO T (2021) A review of SWAL model application in Africa. Water. 13 (9) 1313. https://doi.org/10.3390/W13091313        [ Links ]

ALMENDINGER JE, MURPHY MS and ULRICH JS (2014) Use of the Soil and Water Assessment Tool to scale sediment delivery from field to watershed in an agricultural landscape with topographic depressions. J. Environ. Qual. 43 9-17. https://doi.org/10.2134/jeq2011.0340        [ Links ]

ARNOLD JG and FOHRER N (2005) SWAT2000: Current capabilities and research opportunities in applied watershed modelling. Hydrol. Process. 19 (3) 563-572. https://doi.org/10.1002/hyp.5611        [ Links ]

ARNOLD JG, WILLIAMS JR, SRINIVASEN R and KING KW (1999) SWAT: Soil and Water Assessment Tool-Model Documentation. USDA, Agricultural Research Service, Texas, USA.         [ Links ]

ARNOLD JG, MORIASI DN, GASSMAN PW, ABBASPOUR KC, WHITE MJ, SRINIVASAN R, SANTHI C, HARMEL RD, VAN GRIENSVEN A, VAN LIEW MW, KANNAN N and JHA MK (2012) SWAT: Model use, calibration, and validation. Trans. ASABE. 55 (4) 1491-1508. https://doi.org/10.13031/2013.42256        [ Links ]

ARNOLD JG, YOUSSEF MA, YEN H, WHITE MJ, SHESHUKOV AY, SADEGHI AM, MORIASI DN, STEINER JL, AMATYA DM, SKAGGS RW, HANEY EB, JEONG J, ARABI M and GOWDA PH (2015) Hydrological processes and model representation: impact of soft data on calibration. Trans. ASABE. 58 (6) 1637-1660. https://doi.org/10.13031/trans.58.10726        [ Links ]

BEESON PC, SADEGHI AM, LANG MW, TOMER MD and DAUGHTRY CST (2014) Sediment delivery estimates in water quality models altered by resolution and source of topographic data. J. Environ. Qual. 43 26-36. https://doi.org/10.2134/jeq2012.0148        [ Links ]

BETRI GD, MOHAMED YA, VAN GRIENSVEN A and SRINIVASAN R (2011) Sediment management modelling in the Blue Nile Basin using SWAT model. Hydrol. Earth Syst. Sci. 15 807-818. https://doi.org/10.5194/hess-15-807-2011        [ Links ]

BIEGER K, ARNOLD JG, RATHJENS H, WHITE MJ, BOSCH DD, ALLEN PM, VOLK M and SRINIVASAN R (2017) Introduction to SWAT+, a completely restructured version of the Soil and Water Assessment Tool. /. Am. Water Resour. Assoc. 53 (1) 115-130. https://doi.org/10.1111/1752-1688.12482        [ Links ]

BORAH DK and BERA M (2003) Watershed-scale hydrologic and nonpoint-source pollution models: review of mathematical bases. Trans. ASABE.46(6) 1553-1566. https://doi.org/10.13031/2013.15644        [ Links ]

CLULOW AD, EVERSON CS and GUSH MB (2011) The long-term impact of Acacia mearnsii trees on evaporation, streamflow and groundwater resources. WRC Report No. TT 505/11. Water Research Commission, Pretoria.         [ Links ]

DYE PJ, GUSH MB, EVERSON CS, JARMAIN C, CLULOW A, MENGISTU M, GELDENHUYS CJ, WISE R, SCHOLES RJ, ARCHIBALD S and SAVAGE MJ (2008) Water-use in relation to biomass of indigenous tree species in woodland, forest and/or plantation conditions. WRC Report No. TT 361/08. Water Research Commission, Pretoria. 156 pp.         [ Links ]

EDWARDS MBP and ROBERTS PJT (2006) Managing forests for water: the South African experience. Int. For. Rev. 8 (1) 65. https://doi.org/10.1505/ifor.8.1.65        [ Links ]

EASTON ZM, FUKA DR, WALTER MT, COWAN DM, SCHNEIDERMAN EM and STEENHUIS TS (2008) Re-conceptualizing the Soil and Water Assessment Tool (SWAT) model to predict runoff from variable source areas. J. Hydrol. 348 (3-4) 279-291. https://doi.Org/10.1016/j.jhydrol.2007.10.008        [ Links ]

EVERSON CS, CLULOW AD, BECKER M, WATSON, A, NGUBO C, BULCOCK H, MENGISTU M, LORENTZ S and DEMLIE M (2014) Lhe long term impact of Acacia mearnsii trees on evaporation, streamflow, low flows and ground water resources. Phase II: Understanding the controlling environmental variables and soil water processes over a full crop rotation. WRC Report No. 2022/1/13. Water Research Commission, Pretoria.         [ Links ]

EVERSON CS, CLULOW AD and MENGITSU M (2009) Feasibility study on the determination of riparian evaporation in non-perennial systems. WRC Report No. TT 424/09. Water Research Commission, Pretoria.         [ Links ]

EVERSON CS, GUSH M, MOODLEY M, JARMAIN C, GOVENDER M and DYE P (2007) Can effective management of the riparian zone vegetation to significantly reduce the cost of catchment management and enable greater productivity of land resources. WRC Report No. 1284/1/07. Water Research Commission, Pretoria.         [ Links ]

GASSMAN PW, REYERS MR, GREEN CH and ARNOLD JG (2007) Lhe Soil and Water Assessment Tool: Historical development, applications, and future research directions. Trans. ASABE. 50 (4) 1211-1250. https://doi.org/10.13031/2013.23637        [ Links ]

GASSMAN PW, SADEGHI AM and SRINIVASAN R (2014) Applications of the SWAT model special section: overview and insights. J. Environ. Qual. 43 1-8.         [ Links ]

GUSH MB, DE LANGE WJ, DYE PJ and GELDENHUYS CJ (2015) Water use and socio-economic benefit of the biomass of indigenous trees. Volume 1. WRC Report No. 1876/1/15, Water Research Commission, Pretoria.         [ Links ]

LE MAITRE DC, VAN WILGEN BW, GELDERBLOM CM, BAILEY C, CHAPMAN R and NEL J (2002) Invasive alien trees and water resources in South Africa: case studies of the costs and benefits of management. For. Ecol. Manage. 160 143-159. https://doi.org/10.1016/S0378-1127(01)00474-1        [ Links ]

LEWARNE M (2009) Setting up ArcSWAT hydrological model for the Verlorenvlei Catchment. MSc Thesis, Stellenbosch University.         [ Links ]

MORIASI DN, ARNOLD JG, VAN LIEW MW, BINGNER RL, HARMEL RD and VEITH TL (2007) Model evaluation guidelines for systematic quantification of accuracy in watershed simulations. Trans. ASABE. 50 (3) 885-900. https://doi.org/10.13031/2013.23153        [ Links ]

MUCINA L and RUTHERFORD MC (2006) Lhe vegetation of South Africa, Lesotho and Swaziland. Strelitzia 19. South African National Biodiversity Institute, Pretoria.         [ Links ]

NATIONAL RESOURCES CONSERVATION SERVICE SOIL SURVEY STAFF (1996) National Soil Survey Handbook. 430-VI. US Government Printing Office, Washington, DC        [ Links ]

NEITSCH SL, ARNOLD JG, KINIRY JR and WILLIAMS JR (2011) Soil and Water Assessment Tool Theoretical Documentation Version 2009. Texas Water Resources Institute.         [ Links ]

NABI G, HUSSAIN F, WU R, NANGIA, V, BIBIR and MAJID A (2017) Optimizing micro watershed management for soil erosion control under various slope gradient and vegetation cover conditions using SWAT modelling. Hydrol. Earth Syst. Sci. Discuss, https://doi.org/10.5194/hess-2017-532        [ Links ]

OLBRICH B, OLBRICH K, DYE PJ and SOKO SA (1996) Year-long comparison of water-use efficiency of stressed and non-stressed E. grandis and P. patula: findings and management recommendations. CSIR report FOR-DEA 958. CSIR, Pretoria.         [ Links ]

RAHMAN MM, THOMPSON JR and FLOWER RJ (2016) An enhanced SWAT wetland module to quantify hydraulic interactions between riparian depressional wetlands, rivers and aquifers. Environ. Model. Softw. 84 (2016) 263e289. https://doi.Org/10.1016/j.envsoft.2016.07.003        [ Links ]

SCHULZE RE (1995) Hydrology and agrohydrology: a text to accompany the ACRU-3.00 Agro-Hydrological Modelling System. Water Research Commission Report No. TT 69/95. Water Research Commission, Pretoria. 552 pp.         [ Links ]

SCHULZE RE and HORAN MJC (2007) Hydrological modelling as a tool for ecosystem services trading: case studies from the Drakensberg Region of South Africa. ACRUcons Report, 56. School of Bioresources Engineering and Environmental Hydrology, University of KwaZulu-Natal, Pietermaritzburg.         [ Links ]

SCHUOL J, ABBASPOUR KC, SRINIVASAN R and YANG H (2008) Estimation of fresh water availability in the West African subcontinent using the SWAT hydrologic model. J. Hydrol. 352 30-49. https://doi.Org/10.1016/j.jhydrol.2007.12.025        [ Links ]

SCOTT-SHAW BC (2019) Water-use dynamics of alien plant invaded riparian forests in South Africa. PhD dissertation. School of Agricultural, Earth and Environmental Sciences, College of Agriculture, Engineering and Science, University of KwaZulu-Natal.         [ Links ]

STRAUCH M and VOLK M (2013) SWAT plant growth modification for improved modelling of perennial vegetation in the tropics. J. Ecol. Model. 269 98-112. https://doi.Org/10.1016/j.ecolmodel.2013.08.013        [ Links ]

TAN ML, GASSMAN PW, SRINIVASAN R, ARNOLD JG and YANG XY (2019) A Review of SWAT Studies in Southeast Asia: applications, challenges and future directions. Water. 11 (5) 914. https://doi.org/10.3390/wl1050914        [ Links ]

THAMPI SG, RANEESH SG and SURYA KY (2010) Influence of scale on SWAT model calibration for streamflow in a river basin in the humid tropics. Water Resour. Manage. 24 4567-4578. https://doi.org/10.1007/s11269-010-9676-y        [ Links ]

TOLOSA D (2015) Modelling sediment yield using ArcSWAT and analyzing the most prominent remedial measures: cases of Melka Wakena reservoir, Ethiopia. Addis Ababa University.         [ Links ]

VAN GRIENSVEN A, NDOMBA P, YALEW S and KILONZO T (2012) Critical review of SWAT applications in the upper Nile basin countries. Hydrol. Earth Syst. Sci. 16 3371-3381. https://doi.org/10.5194/hess-16-3371-2012        [ Links ]

WAGNER PD, KUMAR S, FIENER P and SCHNEIDER K (2011) Hydrological modelling with SWAT in a monsoon-driven environment: experience from the Western Ghats, India. Irani. ASABE. 54 1783-1790. https://doi.org/10.13031/2013.39846        [ Links ]

ZABALETA A, MEAURIO M, RUIZ E and ANTIGUEDAD I (2014) Simulation climate change impact on runoff and sediment yield in a small watershed in the Basque Country, northern Spain. J. Environ. Qual. 43 235-245. https://doi.org/10.2134/jeq2012.0209        [ Links ]

 

 

Correspondence:
BC Scott-Shaw
Email: brucecharlesss@gmail.com

Received: 8 October 2020
Accepted: 5 October 2022

Creative Commons License All the contents of this journal, except where otherwise noted, is licensed under a Creative Commons Attribution License