Predicting the effectiveness of different mulching techniques in reducing post-fire runoff and erosion at plot scale with the RUSLE, MMF and PESERA models

Wildfires in Europe: Burned soils require attention

Annually, millions of hectares of land are affected by wildfires worldwide, disrupting ecosystems functioning by affecting on-site vegetation, soil, and above- and belowground biodiversity, but also triggering erosive off-site impacts such as water-bodies contamination or mudflows. Here, we present a soil erosion assessment following the 2017's wildfires at the European scale, including an analysis of vegetation recovery and soil erosion mitigation potential. Results indicate a sharp increase in soil losses with 19.4 million Mg additional erosion in the first post-fire year when compared to unburned conditions. Over five years, 44 million Mg additional soil losses were estimated, and 46% of the burned area presented no signs of full recovery. Post-fire mitigation could attenuate these impacts by 63-77%, reducing soil erosion to background levels by the 4th post-fire year. Our insights may help identifying target policies to reduce land degradation, as identified in the European Union Soil, Forest, and Biodiversity strategies.

Analyzing the Effects of Land Cover Change on the Water Balance for Case Study Watersheds in Different Forested Ecosystems in the USA

We analyzed impacts of interannual disturbance on the water balance of watersheds in different forested ecosystem case studies across the United States from 1985 to 2016 using a remotely sensed long-term land cover monitoring record (U.S. Geological Survey Land Change Monitoring, Assessment, and Projection (LCMAP) Collection 1.0 Science products), gridded precipitation and evaporation data, and streamgaging data using paired watersheds (high and low disturbance). LCMAP products were used to quantify the timing and degree of interannual disturbance and to gain a better understanding of how land cover change affects the water balance of disturbed watersheds. In this paper, we present how LCMAP science products can be used to improve knowledge for hydrologic modeling, climate research, and forest management. Anthropogenic influences (e.g., dams and irrigation diversions) often minimize the impacts of land cover change on water balance dynamics when compared to interannual fluctuations of hydroclimatic events (e.g., drought and flooding). Our findings show that each watershed exhibits a complex suite of influences involving climate variables and other factors that affect each of their water balances differently when land cover change occurs. In this study, forests within arid to semi-arid climates experience greater water balance effects from land cover change than watersheds where water is less limited.

Prediction, validation, and uncertainties of a nation-wide post-fire soil erosion risk assessment in Portugal

Wildfires are a recurrent and increasing threat in mainland Portugal, where over 4.5 million hectares of forests and scrublands have burned over the last 38 years. These fire-affected landscapes have suffered an intensification of soil erosion processes, which can negatively affect soil carbon storage, reduce fertility and forest productivity, and can become a source of pollutants. The main objective of the present study is to produce a post-fire soil erosion risk map for the forest and shrubland areas in mainland Portugal and assess its reliability. To this end, the semi-empirical Morgan–Morgan–Finney erosion model was used to assess the potential post-fire soil erosion according to distinct burn severity and climate scenarios, and the accuracy of the predictions was verified by an uncertainty analysis and validated against independent field datasets. The proposed approach successfully allowed mapping post-fire soil erosion in Portugal and identified the areas with higher post-fire erosion risk for past and future climate extremes. The outcomes of this study comprise a set of tools to help forest managers in their decision-making for post-fire emergency stabilization, ensuring the adequate selection of areas for mitigation to minimize the economic and environmental losses caused by fire-enhanced soil erosion.

Hydrological Response of Burned Soils in Croplands, and Pine and Oak Forests in Zagros Forest Ecosystem (Western Iran) under Rainfall Simulations at Micro-Plot Scale

The post-fire hydrological processes depend on both land use and soil condition (burned or not). This study aims at understanding the variability of the water infiltration, surface runoff and erosion in burned soils under different land uses (forestland and cropland) in comparison to unburned sitesTo this aim, infiltration, runoff and soil losses after a wildfire in two pine and oak forests, and a cropland are evaluated in Zagros forests (Western Iran) using a portable rainfall simulator. This area represents one of the lands with the highest biodiversity and naturalistic value of the entire Middle East, but no similar hydrological evaluations have been conducted so far. The difference in infiltration between the burned and unburned sites under the three land uses was not significant (on the average less than 10%). The runoff and erosion due to the wildfire noticeably increased in the forestland (+95% and 60%, respectively) and slightly decreased in the cropland (−16% and −20%) in comparison to the unburned sites. In the burned croplands erosion requires much attention, because the soil loss is on an average 30-fold compared to the values measured in the forestland. This increase may be even higher, since the rainsplash erosion could be underestimated and the rill or gully erosion was not considered due to the use of a portable rainfall simulator. Therefore, the study suggests the adoption of suitable strategies in croplands of the Zagros forests, in order to limit the negative impacts of high-intensity fires and hydrogeological events. Overall, the study has provided an insight to improve the knowledge on soil hydrology under different land uses and soil conditions. This evaluation helps landscape planners to select the most suitable anti-erosive actions against erosion in fire-affected areas without any needs of long monitoring field campaigns or model implementation.

Exploring and Modeling the Short-Term Influence of Soil Properties and Covers on Hydrology of Mediterranean Forests after Prescribed Fire and Mulching

Several studies have analyzed the changes in individual soil properties and covers and quantified the hydrological response of burned forest soils (with or without post-fire treatment). Less research exists on the influence of these changes on runoff and erosion rates immediately after a prescribed fire and post-fire treatment. Moreover, hydrological modeling of burned areas is based on complex models rather than relying on simple regression equations. This study carries out a combined analysis of the hydrological response of soil and its driving factors in three forests (pine, oak, and chestnut) of Southern Italy that were subjected to prescribed fire and post-fire treatment with mulching. Moreover, simple regression models based on a limited set of soil properties/covers are proposed to predict runoff and erosion. The Principal Component Analysis has shown that the runoff coefficients increase when the water infiltration rate and litter cover decrease and repellency, ash cover, organic carbon content, and bare soil area increase. All the analyzed variables play a secondary role in influencing the sediment concentration. Due to these properties, clear differences in soil properties and covers have been found between unburned and burned soils. The distinctions between the burned soils (mulched or not) are much lower. The proposed regression models use a very low number of soil covers and two dummy variables as input parameters. These models are very accurate in simulating the surface runoff and soil erosion in all soil conditions in the short term.

Assessing Post-Fire Effects on Soil Loss Combining Burn Severity and Advanced Erosion Modeling in Malesina, Central Greece

Earth’s ecosystems are extremely valuable to humanity, playing a key role ecologically, economically, and socially. Wildfires constitute a significant threat to the environment, especially in vulnerable ecosystems, such as those that are commonly found in the Mediterranean. Due to their strong impact on the environment, they provide a crucial factor in managing ecosystems behavior, causing dramatic modifications to land surface processes dynamics leading to land degradation. The soil erosion phenomenon downgrades soil quality in ecosystems and reduces land productivity. Thus, it is imperative to implement advanced erosion prediction models to assess fire effects on soil characteristics. This study focuses on examining the wildfire case that burned 30 km2 in Malesina of Central Greece in 2014. The added value of remote sensing today, such as the high accuracy of satellite data, has contributed to visualizing the burned area concerning the severity of the event. Additional data from local weather stations were used to quantify soil loss on a seasonal basis using RUSLE modeling before and after the wildfire. Results of this study revealed that there is a remarkable variety of high soil loss values, especially in winter periods. More particularly, there was a 30% soil loss rise one year after the wildfire, while five years after the event, an almost double reduction was observed. In specific areas with high soil erosion values, infrastructure works were carried out validating the applied methodology. The approach adopted in this study underlines the significance of using remote sensing and geoinformation techniques to assess the post-fire effects of identifying vulnerable areas based on soil erosion parameters on a local scale.

Modelling forest fire and firebreak scenarios in a mediterranean mountainous catchment: Impacts on sediment loads

Forests provide a number of ecological and hydrological services, for instance, contributing to decreased water and sediment yields through increased infiltration and reduced soil erosion. However, forest fires can turn positive forest services into drawbacks, enhancing surface runoff and soil erosion and damaging both hillslopes and downstream aquatic life in rivers. Therefore, appropriate mitigation strategies should be developed to limit these negative effects. Using a runoff and erosion model (the WaterSed model), we proposed forest fire and firebreak scenarios to analyse their respective effects on sediment loads. The model reproduced the measured discharge and sediment loads over an entire hydrological year, including 21 flood events occurring from November 2010 to May 2011 in a 72-km² Mediterranean catchment (Celone catchment, Puglia, Italy). Eight different forest fire scenarios were then proposed. While the mean burnt areas remained below 2% of the total catchment area, forest fires significantly affected the sediment yield. Indeed, the sediment yield increased over the different forest fire scenarios, from 1.97 to 2.70 t ha^-1.yr^-1, corresponding to a 37% increase. At the flood-event scale, the sediment load after fire represented up to 324% of the unburnt catchment sediment load in the worst-case scenario. By using realistic firebreaks, the sediment load could be dramatically reduced, from 324% to 165%, in the worst-case scenario. Because rural catchments, such as the Celone catchment, are currently experiencing land abandonment, forested areas are expected to replace crops and expand in the future. This change will likely increase forest ecological services, which may, however, be punctually balanced by negative fire effects. More studies addressing the global impacts of forest growth, fires and firebreaks on sediment transfers are therefore needed in similar environments.

Comparing Hydric Erosion Soil Loss Models in Rainy Mountainous and Dry Flat Regions in Portugal

Soil erosion is a severe and complex issue in the agriculture area. The main objective of this study was to assess the soil loss in two regions, testing different methodologies and combining different factors of the Revised Universal Soil Loss Equation (RUSLE) based on Geographical Information Systems (GIS). To provide the methodologies to other users, a GIS open-source application was developed. The RUSLE equation was applied with the variation of some factors that compose it, namely the slope length and slope steepness (LS) factor and practices factor (P), but also with the use of different sources of information. Eight different erosion models (M1 to M8) were applied to the two regions with different ecological conditions: Montalegre (rainy-mountainous) and Alentejo (dry-flat), both in Portugal, to compare them and to evaluate the soil loss for 3 potential erosion levels: 0–25, 25–50 and >50 ton/ha·year. Regarding the methodologies, in both regions the behavior is similar, indicating that the M5 and M6 methodologies can be more conservative than the others (M1, M2, M3, M4 and M8), which present very consistent values in all classes of soil loss and for both regions. All methodologies were implemented in a GIS application, which is free and available under QGIS software.

The Use of Unmanned Aerial Vehicles (UAVs) for Estimating Soil Volumes Retained by Check Dams after Wildfires in Mediterranean Forests

Check dams act as soil collectors during floods, thus retaining a large amount of sediments. The estimation of the soil volumes stored behind a check dam is a key activity for a proper design of these control works and for evaluation of soil delivery after restoration measures at watershed level. Several topographic techniques have been proposed for this activity, but the sediment wedge mapping tools are complex and time consuming. Conversely, the use of unmanned aerial vehicles (UAVs) has been proposed to support aerophotogrammetric techniques for several survey activities with promising results. However, surveys by UAVs have never applied to calculate the size of the sediment wedge behind check dams that are built in fire-affected watersheds, where soil loss and sediment transport may be high after a wildfire. To fill this gap, this study evaluates the efficiency and efficacy of aerophotogrammetric surveys using UAVs to estimate the volume of the sediments stored behind ten check dams, built as post-fire channel treatment in a forest watershed of Castilla La Mancha (Central Eastern Spain). The results of the aerophotogrammetric technique were compared to traditional topographic surveys using a total station and GNSS/RTK, assumed as reference. The estimation of sediment wedge volume provided by UAVs was more accurate (mean RMSE of 0.432), extensive (density of mapped points of 328 m−2) and quick (two days of fieldwork) compared to surveys using the topographic method (RMSE < 0.04 m, six days of field work and density of mapped points of 0.194 m−2) by the topographic method. The differences in the sediment volume estimated by the two methods were not significant, but the UAV method was more accurate for the larger check dams. Moreover, a significant correlation was observed between the volume estimates provided by the two methods, shown by a coefficient of determination close to 0.98. Overall, these results propose a larger use of the aerial surveys for mapping activities in channels regulated by check dams, such as those built for restoration of fire-affected forest watersheds.

A Method for Estimating the Risk of Dam Reservoir Silting in Fire-Prone Watersheds: A Study in Douro River, Portugal

Forest fires are an increasing problem over recent decades. The fires, among other consequences, lead to an increase in the soil vulnerability to water erosion and a consequent increase in sedimentation rates. When barriers are present, such as dams or weirs, there is an amplified risk of sediment and ash deposition in their reservoirs, causing siltation. Thus, there is an interest in studying in more detail the risk of siltation of barriers and reservoirs in the Douro River watershed following wildfires. A detailed barrier inventory was lacking for the Douro River, hampering the identification of siltation-prone areas. In order to fill in this gap, an extensive inventory of barriers in the Douro river basin was carried out for the present study. The result was an abundant and reliable dataset on the Douro River barriers, which allowed a prognosis on the watershed siltation risk. The method for calculating the siltation risk relied on the relationship between the frequency of forest fires, the erosion risk and the frequency of reservoirs. The sub-basins with the greater siltation risk are the Tâmega, Corgo, Sousa and Paiva river basins. Most reservoirs with the highest siltation risk were from small dams. The modelling results were compared with stream connectivity and concentrations of stream water phosphorus (associated with the sediments that flow into the rivers due to the fires). With regard to connectivity, only two reservoirs were at high risk of sedimentation due to fires, so the categories of connectivity risk and fire-based sedimentation risk are probably not related. With regard to risk of high phosphorus loadings, in 8 basins the upper classes for fire-based erosion risk coincided with the upper class for phosphorus loadings suggesting that high phosphorus loading could be associated with fire-based erosion. This study works as a simple but reliable example on the assessment and mapping of siltation risk in stream networks intersected by abundant barriers. It allowed for identifying barriers that can accumulate a large quantity of fine sediments and ashes, interfering with water quality and soil erosion as well as with the storage capacity of the respective barriers.

Estimating Soil Erosion Rate Changes in Areas Affected by Wildfires

In recent decades, wildfires have become a serious threat worldwide, producing disasters in the natural and anthropogenic environment as well as serious economic losses. One of wildfire’s major impacts is soil erosion, as it may cause major problems in both the physical and anthropogenic environment and seriously affect the landscape. This study investigates the soil erosion rate changes in areas affected by wildfires and uses, as a pilot area, the drainage basin of the Pinios earth-filled dam located in the Ilia Regional Unit, western Greece, which has suffered serious erosion changes after a wildfire event. For this purpose, the Revised Universal Soil Loss Equation (RUSLE) is applied in GIS software, and the soil erosion rate changes in the selected investigation area are estimated at different time intervals. Specifically, soil erosion rate changes are calculated by importing the factors from the RUSLE equation in the GIS software and uses as a dependent variable the cover management factor C, which is strongly influenced by large destructive fires. The models that are produced are compared with each other by collating average annual soil erosion maps and rates before the fire, immediately after the fire and for the existing conditions occurring in the pilot area.

Spatial and temporal evaluation of soil erosion in Turkey under climate change scenarios using the Pan-European Soil Erosion Risk Assessment (PESERA) model

The impacts of climate change on soil erosion are mainly caused by the changes in the amount and intensity of rainfall and rising temperature. The combination of rainfall and temperature change is likely to be accompanied by negative or positive variations in agricultural and forest management. Turkey contains vast fertile plains, high mountain chains and semi-arid lands, with a climate that ranges from marine to continental and therefore is susceptible to soil erosion under climate change, particularly on high gradients and in semi-arid areas. This study aims to model the soil erosion risk under climate change scenarios in Turkey using the Pan-European Soil Erosion Assessment (PESERA) model, predicting the likely effects of land use/cover and climate change on sediment transport and soil erosion in the country. For this purpose, PESERA was applied to estimate the monthly and annual soil loss for 12 land use/cover types in Turkey. The model inputs included 128 variables derived from soil, climate, land use/cover and topography data. The total soil loss from the land surface is speculated to be approximately 285.5 million tonnes per year. According to the IPCC 5th Assessment Report of four climate change scenarios, the total soil losses were predicted as 308.9, 323.5, 320.3 and 355.3 million tonnes for RCP2.6, RCP4.5, RCP6.0 and RCP8.5 scenarios respectively from 2060 to 2080. The predicted amounts of fertile soil loss from agricultural land in a year were predicted to be 55.5 million tonnes at present, and 62.7, 59.9, 61.7 and 58.1 under RCP2.6, RCP4.5, RCP6.0 and RCP8.5 respectively. This confirms that approximately 30% of the total erosion occurs over the agricultural lands. In this respect, degraded forests, scrub and arable lands were subjected to the highest erosion rate (68%) of the total, whereas, fruit trees and berry plantations reflected the lowest erosion rates. Low soil organic carbon, sparse vegetation cover and variable climatic conditions significantly enhanced the erosion of the cultivated lands by primarily removing the potential food for organisms. Finally, process-based models offer a valuable resource for decision-makers when improving environmental management schemes and also decrease uncertainty when considering risks.

A Method to Improve the Flood Maps Forecasted by On-Line Use of 1D Model

Forecasting floods in urban areas during a heavy rainfall is the aim of every early warning system. 2D-models produce the most accurate flood maps, but they are practically useless as quasi real-time tools, because their run times are comparable to times of propagation of floods. Run times of 1D-model are of tens of seconds, but their predictions lack accuracy and many useful indicators of flood severity. Our aim is the identification of the 2D-model map that is more similar to the actual map, chosen among those simulated off-line. To this aim, we produce a rough flood map of the occurring event, through a quasi real-time simulation of the rainfall-runoff using a 1D-model. Then we apply an original method, named “ranking approach”, to perform the best matching. This method is applied to the Corace torrent (Calabria, Southern Italy), using 17 synthetic hyetographs to simulate the same number of rainfall-runoff events, using 1D (SWMM) and 2D (MIKE) models. The method proves to be effective in 65% of the cases, while in 82% of cases (i.e., for 14 cases out 17), the event produced by the same ietograph falls within the third rank.

Fire severity and soil erosion susceptibility mapping using multi-temporal Earth Observation data: The case of Mati fatal wildfire in Eastern Attica, Greece

In recent years, forest fires have increased in terms of frequency, extent and intensity, especially in Mediterranean countries. Climate characteristics and anthropogenic disturbances lead forest environments to display high vulnerability to wildfires, with their sustainability being threatened by the loss of vegetation, changes on soil properties, and increased soil loss rates. Moreover, wildfires are a great threat to property and human life, especially in Wildland-Urban Interface (WUI) areas. In light of the impacts and trends mentioned above, this study aims to assess the impact of the Mati, Attika wildfire on soil erosion. The event caused 102 fatalities, inducing severe consequences to the local infrastructure network; economy; and natural resources. As such, the Revised Universal Soil Loss Equation (RUSLE) was implemented (pre-; post-fire) at the Rafina, Attika watershed encompassing the Mati WUI. Fire severity was evaluated based on the Normalized Burn Ratio (NBR). This index was developed utilizing innovative remotely sensed Earth Observation data (Sentinel-2). The high post-fire values indicate the fire's devastating effects on vegetation loss and soil erosion. A critical "update" was also made to the CORINE Land Cover (CLC) v. 2018, by introducing a new land use class namely "Urban Forest", in order to distinguish the WUI configuration. Post-fire erosion rates are notably higher throughout the study area (4.53-5.98 t ha^-1 y^-1), and especially within the WUI zone (3.75-18.58 t ha^-1 y^-1), while newly developed and highly vulnerable cites now occupy the greater Mati area. Furthermore, archive satellite data (Landsat-5) revealed how the repeated (historical) wildfires have ultimately impacted vegetation recovery and erosional processes. To our knowledge this is the first time that RUSLE is used to simulate soil erosion at a WUI after a fire event, at least at a Mediterranean basin. The realistic results attest that the model can perform well at such diverse conditions, providing a solid basis for soil loss estimation and identification of high-risk erosion areas.

Assessment of soil erosion risk and its response to climate change in the mid-Yarlung Tsangpo River region

Soil erosion is sensitive to climate change, especially in high mountain areas. The Tibetan Plateau has experienced dramatic land surface environment changes under the impact of climate change during the last decades. In this study, we focused on the mid-Yarlung Tsangpo River (MYZ River) located in the southern part of the Tibetan Plateau. The revised universal soil loss equation (RUSLE) was applied to assess soil erosion risk. To increase its applicability to high mountain areas with longer periods of snowfall, snowmelt runoff erosivity was considered in addition to rainfall erosivity. Results revealed that soil erosion of the MYZ River region was of a moderate grade with an average soil erosion rate of 29.1 t ha^-1 year^-1 and most serious erosion in wet and cold years. Soil erosion rate in the MYZ River region showed a decreasing trend of - 1.14% year^-1 due to the precipitation, temperature, and vegetation changes from 2001 to 2015, with decreasing precipitation being the most important factor. Increasing precipitation and temperature would lead to increasing soil erosion risk in ~ 2050 based on the Coupled Model Intercomparison Project (CMIP5) and RUSLE models. It is clear that soil erosion in high mountain areas greatly depends on climate state and attentions should be paid to address soil erosion problem in the future.

Hydrological Impacts of Large Fires and Future Climate: Modeling Approach Supported by Satellite Data

Fires have significant impacts on soil erosion and water supply that may be exacerbated by future climate. The aims of this study were: To simulate the effects of a large fire event in the SWAT (Soil and Water Assessment Tool) hydrological model previously calibrated to a medium-sized watershed in Portugal; and to predict the hydrological impacts of large fires and future climate on water supply and soil erosion. For this, post-fire recovery was parametrized in SWAT based on satellite information, namely, the fraction of vegetation cover (FVC) calculated from the normalized difference vegetation index (NDVI). The impact of future climate was based on four regional climate models under the stabilization (RCP 4.5) and high emission (RCP 8.5) scenarios, focusing on mid-century projections (2020–2049) compared to a historical period (1970–1999). Future large fire events (>3000 ha) were predicted from a multiple linear regression model, which uses the daily severity rating (DSR) fire weather index, precipitation anomaly, and burnt area in the previous three years; and subsequently simulated in SWAT under each climate model/scenario. Results suggest that time series of satellite indices are useful to inform SWAT about vegetation growth and post-fire recovery processes. Different land cover types require different time periods for returning to the pre-fire fraction of vegetation cover, ranging from 3 years for pines, eucalypts, and shrubs, to 6 years for sparsely vegetated low scrub. Future climate conditions are expected to include an increase in temperatures and a decrease in precipitation with marked uneven seasonal distribution, and this will likely trigger the growth of burnt area and an increased frequency of large fires, even considering differences across climate models. The future seasonal pattern of precipitation will have a strong influence on river discharge, with less water in the river during spring, summer, and autumn, but more discharge in winter, the latter being exacerbated under the large fire scenario. Overall, the decrease in water supply is more influenced by climate change, whereas soil erosion increase is more dependent on fire, although with a slight increase under climate change. These results emphasize the need for adaptation measures that target the combined hydrological consequences of future climate, fires, and post-fire vegetation dynamics.

Soil erosion modelling: The new challenges as the result of policy developments in Europe

New challenges and policy developments after 2015 (among others, the Common Agricultural Policy (CAP), Sustainable Development Goals (SDGs)) are opportunities for soil scientists and soil erosion modellers to respond with more accurate assessments and solutions as to how to reduce soil erosion and furthermore, how to reach Zero Net Land Degradation targets by 2030. This special issue includes papers concerning the use of fallout for estimating soil erosion, new wind erosion modelling techniques, the importance of extreme events (forest fires, intense rainfall) in accelerating soil erosion, management practices to reduce soil erosion in vineyards, the impact of wildfires in erosion, updated methods for estimating soil erodibility, comparisons between sediment distribution models, the application of the WaTEM/SEDEM model in Europe, a review of the G2 model and a proposal for a land degradation modelling approach. New data produced from field surveys such as LUCAS topsoil and the increasing availability of remote sensing data may facilitate the work of erosion modellers. Finally, better integration with other soil related disciplines (soil carbon, biodiversity, compaction and contamination) and Earth Systems modelling is the way forward for a new generation of erosion process models.

Vetiver grass hedgerows significantly reduce nitrogen and phosphorus losses from fertilized sloping lands

Runoff and over-use of fertilizers have been considered as two major factors accelerating the discharge of nitrogen (N) and phosphorus (P) from agricultural fields to surface water. The practice of vetiver grass hedgerows (VGH) can check sediments and runoff pollutants from agricultural fields. However, the efficiency of VGH in reducing N and P losses while maintaining optimum crop yields is still unclear under a recommended fertilization rate. A three-year field experiment was conducted on a 10 ^o sloping land to know how VGH can reduce the discharge of runoff nutrients to surface water bodies and maintain optimum crop yields, and to understand the relationships between changing soil properties and reduction of sediments N and P due to the adoption of VGH. Five fertilization treatments to VGH were examined under VGH plus organic fertilizer (VGH + OF), VGH plus inorganic fertilizer (VGH + IF), sole organic or inorganic fertilizer (OF or IF) and no VGH and fertilizer (Control). Runoff nutrient pollutants PO₄^-, NO₃^--N and NH₄⁺-N were significantly (P < 0.01) reduced by VGH + OF compared to OF by 97%, 94% and 95% and VGH + IF compared to IF by 95%, 88% and 89% respectively for 2012, 2013 and 2014. Sediment nutrients N and P were significantly (P < 0.01) reduced by VGH + OF compared to OF by 98% and 99%, and VGH + IF compared to IF by 94% and 99%, respectively. Improved soil properties by VGH significantly (P < 0.01) reduced runoff pollutants and consequently increased maize yields. Our results imply that runoff erosion, rather than fertilization, is a major driving force for agriculture-derived water pollution. Adoption of VGH with a recommended fertilization rate could significantly reduce N and P nutrient losses from agricultural fields and consequently improve water quality as well as maintaining optimum crop yields on sloping lands.