Modelling monthly soil losses and sediment yields in Cyprus

The aim of this study was to map soil erosion on the Mediterranean island of Cyprus. The G2 model, an empirical model for month-time step erosion assessments, was used. Soil losses in Cyprus were mapped at a 100?m cell size, while sediment yields at a sub-basin scale of 0.62?km² mean size. The results indicated a mean annual erosion rate of 11.75?t?ha^?1?y^?1, with October and November being the most erosive months. The 34% of the island's surface was found to exceed non-sustainable erosion rates (>10?t?ha^?1?y^?1), with sclerophyllous vegetation, coniferous forests, and non-irrigated arable land being the most extensive non-sustainable erosive land covers. The mean sediment delivery ratio (SDR) was found to be 0.26, while the mean annual specific sediment yield (SSY) value for Cyprus was found to be 3.32?t?ha^?1?y^?1. The annual sediment yield of the entire island was found to be 2.746?Mt?y^?1. This study was the first to provide complete and detailed erosion figures for Cyprus at a country scale. The geodatabase and all information records of the study are available at the European Soil Data Centre (ESDAC) of the Joint Research Centre (JRC).     

Timing of erosion and satellite data: A multi-resolution approach to soil erosion risk mapping

Erosion reduces soil productivity and causes negative downstream impacts. Erosion processes occur on areas with erodible soils and sloping terrain when high-intensity rainfall coincides with limited vegetation cover. Timing of erosion events has implications on the selection of satellite imagery, used to describe spatial patterns of protective vegetation cover. This study proposes a method for erosion risk mapping with multi-temporal and multi-resolution satellite data. The specific objectives of the study are: (1) to determine when during the year erosion risk is highest using coarse-resolution data, and (2) to assess the optimal timing of available medium-resolution images to spatially represent vegetation cover during the high erosion risk period. Analyses were performed for a 100-km² pasture area in the Brazilian Cerrados. The first objective was studied by qualitatively comparing three-hourly TRMM rainfall estimates with MODIS NDVI time series for one full year (August 2002–August 2003). November and December were identified as the months with highest erosion risk. The second objective was examined with a time series of six available ASTER images acquired in the same year. Persistent cloud cover limited image acquisition during high erosion risk periods. For each ASTER image the NDVI was calculated and classified into five equally sized classes. Low NDVI was related to high erosion risk and vice versa. A DEM was used to set approximately flat zones to very low erosion risk. The six resulting risk maps were compared with erosion features, visually interpreted from a fine-resolution QuickBird image. Results from the October ASTER image gave highest accuracy (84%), showing that erosion risk mapping in the Brazilian Cerrados can best be performed with images acquired shortly before the first erosion events. The presented approach that uses coarse-resolution temporal data for determining erosion periods and medium-resolution data for effective erosion risk mapping is fast and straightforward. It shows good potential for successful application in other areas with high spatial and temporal variability of vegetation cover.     

Risk assessment of soil erosion in semi-arid mountainous watershed in Saudi Arabia by RUSLE model coupled with remote sensing and GIS

Soil erosion is the most important factor in land degradation and influences desertification in semi-arid areas. A comprehensive methodology that integrates revised universal soil loss equation (RUSLE) model and GIS was adopted to determine the soil erosion risk (SER) in semi-arid Aseer region, Saudi Arabia. Geoenvironmental factors viz. rainfall (R), soil erodibility (K), slope (LS), cover management and practice factors were computed to determine their effects on average annual soil loss. The high potential soil erosion, resulting from high denuded slope, devoid of vegetation cover and high intensity rainfall, is located towards the north western part of the study area. The analysis is investigated that the SER over the vegetation cover including dense vegetation, sparse vegetation and bushes increases with the higher altitude and higher slope angle. The erosion maps generated with RUSLE integrated with GIS can serve as effective inputs in deriving strategies for land planning/management in the environmentally sensitive mountainous areas.     

Fractional vegetation cover estimation in heterogeneous areas by combining a radiative transfer model and a dynamic vegetation model

ABSTRACT

A fractional vegetation cover (FVC) estimation method incorporating a vegetation growth model and a radiative transfer model was previously developed, which was suitable for FVC estimation in homogeneous areas because the finer-resolution pixels corresponding to one coarse-resolution FVC pixel were all assumed to have the same vegetation growth model. However, this assumption does not hold over heterogeneous areas, meaning that the method cannot be applied to large regions. Therefore, this study proposes a finer spatial resolution FVC estimation method applicable to heterogeneous areas using Landsat 8 Operational Land Imager reflectance data and Global LAnd Surface Satellite (GLASS) FVC product. The FVC product was first decomposed according to the normalized difference vegetation index from the Landsat 8 OLI data. Then, independent dynamic vegetation models were built for each finer-resolution pixel. Finally, the dynamic vegetation model and a radiative transfer model were combined to estimate FVC at the Landsat 8 scale. Validation results indicated that the proposed method (R²?=?0.7757, RMSE?=?0.0881) performed better than either the previous method (R²?=?0.7038, RMSE?=?0.1125) or a commonly used method involving look-up table inversions of the PROSAIL model (R²?=?0.7457, RMSE?=?0.1249).     

Soil erosion modelled with USLE and PESERA using QuickBird derived vegetation parameters in an alpine catchment

The focus of soil erosion research in the Alps has been in two categories: (i) on-site measurements, which are rather small scale point measurements on selected plots often constrained to irrigation experiments or (ii) off-site quantification of sediment delivery at the outlet of the catchment. Results of both categories pointed towards the importance of an intact vegetation cover to prevent soil loss. With the recent availability of high-resolution satellites such as IKONOS and QuickBird options for detecting and monitoring vegetation parameters in heterogeneous terrain have increased. The aim of this study is to evaluate the usefulness of QuickBird derived vegetation parameters in soil erosion models for alpine sites by comparison to Cesium-137 (Cs-137) derived soil erosion estimates. The study site (67 km²) is located in the Central Swiss Alps (Urseren Valley) and is characterised by scarce forest cover and strong anthropogenic influences due to grassland farming for centuries. A fractional vegetation cover (FVC) map for grassland and detailed land-cover maps are available from linear spectral unmixing and supervised classification of QuickBird imagery. The maps were introduced to the Pan-European Soil Erosion Risk Assessment (PESERA) model as well as to the Universal Soil Loss Equation (USLE). Regarding the latter model, the FVC was indirectly incorporated by adapting the C factor. Both models show an increase in absolute soil erosion values when FVC is considered. In contrast to USLE and the Cs-137 soil erosion rates, PESERA estimates are low. For the USLE model also the spatial patterns improved and showed “hotspots” of high erosion of up to 16 t ha⁻¹ a⁻¹. In conclusion field measurements of Cs-137 confirmed the improvement of soil erosion estimates using the satellite-derived vegetation data.     

Estimating vegetation parameter for soil erosion assessment in an alpine catchment by means of QuickBird imagery

Soil erosion rates in alpine regions are related to high spatial variability complicating assessment of risk and damages. A crucial parameter triggering soil erosion that can be derived from satellite imagery is fractional vegetation cover (FVC). The objective of this study is to assess the applicability of normalized differenced vegetation index (NDVI), linear spectral unmixing (LSU) and mixture tuned matched filtering (MTMF) in estimating abundance of vegetation cover in alpine terrain. To account for the small scale heterogeneity of the alpine landscape we used high resolved multispectral QuickBird imagery (pixel resolution = 2.4 m) of a site in the Urseren Valley, Central Swiss Alps (67 km²). A supervised land-cover classification was applied (total accuracy 93.3%) prior to the analysis in order to stratify the image. The regression between ground truth FVC assessment and NDVI as well as MTMF-derived vegetation abundance was significant (r² = 0.64, r² = 0.71, respectively). Best results were achieved for LSU (r² = 0.85). For both spectral unmixing approaches failed to estimate bare soil abundance (r² = 0.39 for LSU, r² = 0.28 for MTMF) due to the high spectral variability of bare soil at the study site and the low spectral resolution of the QuickBird imagery. The LSU-derived FVC map successfully identified erosion features (e.g. landslides) and areas prone to soil erosion. FVC represents an important but often neglected parameter for soil erosion risk assessment in alpine grasslands.     

Temporal and spatial relationships between soil erosion and ecological restoration in semi-arid regions: a case study in northern Shaanxi,China

ABSTRACT

To assess the effects of the Grain for Green Program (GGP) on soil erosion is essential to support better land management policies in the Chinese Loess Plateau. Studies on the evaluation of the effects of the GGP on soil erosion have garnered heightened attention. However, few studies examined the efficiency of GGP on soil erosion control through spatial relationship analysis. Thus, this study focuses on analyzing the spatial variation relationship between soil erosion and GGP in northern Shaanxi, Chinese Loess Plateau, from 1988 to 2015. The Universal Soil Loss Equation was used to quantify changes in soil erosion at the regional and watershed scales, and the Geographically Weighted Regression model was used to analyze the spatial relationships between land use and land cover (LULC) and soil erosion. Our results indicated that the major characteristic of LULC change during the GGP was a rapid increase of vegetation area and a rapid decrease of cropland. Bare lands contributed to the most serious soil loss, followed by croplands and sparse grasslands. The GGP had a globally positive influence on the decrease in soil erosion over the study area, but the amount of soil erosion in western and northern regions maintained a severe level. Spatial heterogeneity in the nature of the relationships among different vegetation, croplands, and soil erosion was also observed. The change rate of wood and the change rate of soil erosion in northern sub-watershed represented a negative relationship, while the change rate of sparse grassland was negatively correlated to the change rate of soil erosion in 21 sub-watersheds, account for 72% of the study area. The GGP implemented in northern sub-watersheds were more effective for soil erosion control than southern sub-watersheds. We propose that current areas of vegetation can support soil erosion control in the whole northern Shaanxi, but local-scale ecological restoration can be considered in northern sub-watersheds.     

An assessment of gully erosion along major armoured roads in south-eastern region of South Africa: a remote sensing and GIS approach

An assessment of gully erosion along road drainage-release sites is critical for understanding the contribution of roads to soil loss and for informed land management practices. Considering that road-related gully erosion has traditionally been measured using field methods that are expensive, tedious and limited spatially as well as temporally, it is important to identify affordable, timely and robust methods that can be used to effectively map and estimate the volume of gullies along the road networks. In this study, gullies along major roads were identified from remotely sensed data sets and their volumes were estimated in a Geographic Information Systems environment. Also, the biophysical and climatic factors such as vegetation cover, the road contributing surface area, the gradient of the discharge hillslope and rainfall were derived from remotely sensed data sets using Geographic Information Systems techniques to find out whether they could explain the morphology of gullies that existed in this area. The results of this study indicate that hillslope gradient (R²?=?0.69, α = 0.00) and road contributing surface area (R²?=?0.63, α = 0.00) have a strong influence on the volume of gullies along the major roads in the south-eastern region of South Africa, as might have been expected. However, other factors such as vegetation cover (R² = 0.52, α = 0.00) and rainfall (R² = 0.41 and α = 0.58) have a moderately weaker influence on the overall volume of gullies. Overall, the findings of this study highlight the importance of using remote sensing and Geographic Information Systems technologies in investigating gully erosion occurrence along major roads where detailed field work remains a challenge.     

Detecting and mapping levels of Gonipterus scutellatus-induced vegetation defoliation and leaf area index using spatially optimized vegetation indices

Gonipterus scutellatus outbreaks may severely defoliate Eucalyptus plantations growing in South Africa. Therefore, detecting and mapping the severity and extent of G. scutellatus defoliation is essential for the deployment of suppressive measures. In this study, we tested the utility of spatially optimized vegetation indices and an artificial neural network in detecting and mapping G. scutellatus-induced vegetation defoliation, using both visual estimates of percentage defoliation and optical leaf area index (LAI) measures. We tested both field methods to determine which of the two were more superior in detecting vegetation defoliation using optimized vegetation indices. These indices were computed from a WorldView-2 pan-sharpened image, which is characterized with a 0.5-m spatial resolution and eight spectral bands. The indices were resampled to spatial resolutions that best represented levels of G. scutellatus-induced defoliation. The results showed that levels of defoliation, using visual percentage estimates, were detected with an R² of 0.83 and an RMSE of 1.55 (2.97% of the mean measured defoliation), based on an independent test data-set. Similarly, LAI subjected to defoliation was detected with an R² of 0.80 and an RMSE of 0.03 (0.06% of the mean measured LAI), based on an independent test data-set. Therefore, the results indicate that the cheaper less-complicated visual percentage estimates of defoliation was the more superior model of the two. A sensitivity analysis revealed that NDRE, MCARI2 and ARI ranked as the top three most influential indices in developing both percentage defoliation and LAI models. Furthermore, we compared the optimized model with a model developed using the original image spatial resolution. The results indicated that the optimized model performed better than the original 0.5-m spatial resolution model. Overall, the study showed that vegetation indices optimized to specific spatial resolutions can effectively detect and map levels of G. scutellatus-induced defoliation and LAI subjected to defoliation.     

Quantitative assessment of vegetation cover and soil degradation factors within terrain units for planning,monitoring and assessment of renaturation along oil and gas pipelines

This paper evaluates the renaturation activities applying the quantification of vegetation cover (VC), the site suitability analysis (SSA) based on the predefined criteria (slope steepness category (SSC), soil erodibility factor (K) and VC) and soil erosion model (SEM) results within the terrain units (TUs) along pipeline rights-of-way (RoW). Quantification of VC percentage is performed to assess the overall restored VC from 2005 to 2007. The results of the quantitative analysis in 2007 show that the total area of restored VC is 10.7 km², and 8.9 km² still needs to be restored to comply with the environmental acceptance criteria. As a result of SSA, TUs were prioritized by erosion vulnerability and this allowed to better understand the landscape behaviour in regards to erosion processes. SEM provided more detailed predictions of erosion classes falling into TUs. SEM identified 40% of erosion sites occurred from 2005 to 2010.     

Geographic information system for semi‐detailed mapping of soils in a semi‐arid region

Abstract

The paper presents a GIS model for mapping soils in the semi arid region of Israel. The model is based on a priori knowledge of the soil generating factors in the study area, namely (1) the parent‐material (lithology) which determines the origin of the soil; and (2) the relief (including the drainage patents) which is responsible for erosion, deposition, and leaching processes along the catena. A special attempt was made to represent soil belts in small subwatersheds with different parent‐material. The width of these belts can be varied as a function of the local relief.

We believe that GIS techniques, in contrast to conventional survey and mapping methods, have the potential to overcome the manpower and fiscal restrictions which limit monitoring of large areas and areas which are for different reasons difficult to access. However, systematically generalizing a spatial model for a large area may introduce some errors, either due to local variations which were not taken into account, or regional variations.     

土壤水力侵蚀的遥感信息模型研究——以江西省为例

通过利用Terra/Aqua卫星上搭载的MODIS传感器计算获取的16d合成植被指数产品(MOD13A2),进一步按照最大值合成法计算月合成光谱植被指数,按照USLE模型月模式评价江西省2005年土壤侵蚀,并与传统的USLE模型年模式计算的结果进行了比较。     

Soil erosion modeling using MMF model -A remote sensing and GIS perspective

Hardly any part of the world has remained unchanged since the arrival of the speciesHomo sapiens including the mountain ecosystems. Himalayan physiographic unit of India in due course has become populated and is tolerating all kinds of human interventions. Soil erosion in this region has been identified as a major problem due to both natural and anthropogenic factors. Remote sensing and Geographical Information system (GIS) techniques hold great promises in the assessment and conservation of natural resources including the surface soil. The major objective of the present study was to apply a process based model to quantify soil erosion and to prioritize the sub-watershed on this basis. The sub-watershed located at Jakhan rao area of Western Dun in lower Himalayan belt was taken as the test site for the study at 1: 50,000 scale. Deforestation, unscientific agricultural practices, terrace farming, cattle grazing and land degradation in the sub-watershed are some of the anthropogenic factors causing soil erosion in the area. Here, MMF model was used for estimation of soil erosion by incorporating layers derived from both remote sensing and ancillary data. IRS 1C LISS III satellite data was used for the preparation of land use map that was used to derive RD map, BD map and K map. Digital Elevation Model (DEM) provided slope map, an intermediate layer used in equation 6 to generate G map, and soil map provided MS map, BD map and K map. The above intermediate layers generated were then integrated in GIS domain to estimate the amount of soil erosion in the sub-watershed area. Results show high values 4572.333 kg/m² for G map, which depicted transport capacity of overland flow. Comparatively lower values 13.15, and 7.98 kg/m²were observed for F map, which depicted soil detachment by raindrop impact. The subtracted image of the aforesaid layers produced the real picture, where in the highest value 3.770 kg/m² was found in the midland region of the site. The crossed erosion map was then classified into different erosion classes for sub-watershed area. This study illustrates the applications of remote sensing and GIS techniques for soil erosion modeling.     

Impact of the spatial resolution of climatic data and soil physical properties on regional corn yield predictions using the STICS crop model

The assimilation of Earth observation (EO) data into crop models has proven to be an efficient way to improve yield prediction at a regional scale by estimating key unknown crop management practices. However, the efficiency of prediction depends on the uncertainty associated with the data provided to crop models, particularly climatic data and soil physical properties. In this study, the performance of the STICS (Simulateur mulTIdisciplinaire pour les Cultures Standard) crop model for predicting corn yield after assimilation of leaf area index derived from EO data was evaluated under different scenarios. The scenarios were designed to examine the impact of using fine-resolution soil physical properties, as well as the impact of using climatic data from either one or four weather stations across the region of interest. The results indicate that when only one weather station was used, the average annual yield by producer was predicted well (absolute error <5%), but the spatial variability lacked accuracy (root mean square error = 1.3 t ha⁻¹). The model root mean square error for yield prediction was highly correlated with the distance between the weather stations and the fields, for distances smaller than 10 km, and reached 0.5 t ha⁻¹ for a 5-km distance when fine-resolution soil properties were used. When four weather stations were used, no significant improvement in model performance was observed. This was because of a marginal decrease (30%) in the average distance between fields and weather stations (from 10 to 7 km). However, the yield predictions were improved by approximately 15% with fine-resolution soil properties regardless of the number of weather stations used. The impact of the uncertainty associated with the EO-derived soil textures and the impact of alterations in rainfall distribution were also evaluated. A variation of about 10% in any of the soil physical textures resulted in a change in dry yield of 0.4 t ha⁻¹. Changes in rainfall distribution between two abundant rainfalls during the growing season led to a significant change in yield (0.5 t ha⁻¹ on average). Our results highlight the importance of using fine-resolution gridded daily precipitation data to capture spatial variations of rainfall as well as using fine-resolution soil properties instead of coarse-resolution soil properties from the Canadian soil dataset, especially for regions with high pedodiversity.     

Biome-level relationships between vegetation indices and climate variables using time-series analysis of remotely-sensed data

ABSTRACT

Climatic factors such as rainfall and temperature play a vital role in the growth characteristics of vegetation. While the relationship between climate and vegetation growth can be accurately predicted in instances where vegetation is homogenous, this becomes complex to determine in heterogeneous vegetation environments. The aim of this paper was to study the relationship between remotely-sensed monthly vegetation indices (i.e. Normalized Difference Vegetation Index and Enhanced Vegetation Index) and climatic variables (temperature and precipitation) using time-series analysis at the biome-level. Specifically, the autoregressive distributed lag model (ARDL1 and ARDL2, corresponding respectively to one month and two month lags) and the Koyck-transformed distributed lag model were used to build regression models. All three models estimated NDVI and EVI fairly accurately in all biomes (Relative Root-Mean-Squared-Error (RMSE): 12.0–26.4%). Biomes characterized by relative homogeneity (Grassland, Savanna, Indian Ocean Coastal Belt and Forest Biomes) achieved the most accurate estimates due to the dominance of a few species. Comparisons of lag size (one month compared to two months) generally showed similarities (Akaike Information Criterion (AIC), Bayesian Information Criterion (BIC) and log-likelihood) with quite high comparability in certain biomes – this indicates the utility of the ARDL1 and ARDL2 model, depending on the availability of appropriate data. These findings demonstrate the variation in estimation linked to the biome, and thus the validity of biome-level correlation of climatic data and vegetation indices.     

Integrating terrain and vegetation indices for identifying potential soil erosion risk area

The present paper offers an innovative method to monitor the change in soil erosion potential by integrating terrain and vegetation indices derived from remote sensing data. Three terrain indices namely, topographic wetness index (TWI), stream power index (SPI) and slope length factor (LS), were derived from the digital elevation model. Normalized vegetation index (NDVI) was derived for the year 1988 and 2004 using remote sensing images. K-mean clustering was performed on staked indices to categorize the study area into four soil erosion potential classes. The validation of derived erosion potential map using USLE model showed a good agreement. Results indicated that there was a significant change in the erosion potential of the watershed and a gradual shifting of lower erosion potential class to next higher erosion potential class over the study period.     

Basin scale rainfall-evapotranspiration dynamics in a tropical semiarid environment during dry and wet years

In semiarid regions the occurrence of alternating long drought and heavy rainfall periods directly impacts water availability, affecting human water supply, agriculture development and the provision of ecosystem services. Because of that, research on the water input and output fluxes at the basin scale is of paramount importance. In this sense, rainfall-evapotranspiration (ET) dynamics play a critical role in water, soil and vegetation interactions, in hydrometeorological modelling and in the energy fluxes dynamics of semiarid regions. Therefore, the objective of this study was to quantify daily ET during a wet year and a dry year in a watershed located in the Brazilian Semiarid, by using remote sensing data and formulations based on the Simplified Surface Energy Balance Index (S-SEBI) and the Simplified Surface Energy Balance (SSEB) algorithms. Land surface temperature, albedo and NDVI data from MODIS sensor and solar radiation data from weather stations located in the basin were used. Rainfall analysis indicated 2009 and 2012 as being representatives of anomalously wet and dry years respectively, which were selected for the quantification of ET. The proposed algorithm was adjusted and verified with data from a flux tower equipped with eddy covariance system. Daily remote sensing ET estimates showed good agreement with observed values (RMSE = 0.79 mm.d⁻¹) and the annual ET relative error was of 7.7% (35.4 mm.year⁻¹). Results showed that the native vegetation can delay its dormant state for five months during wet years. During the wet year, ET differences between land cover classes were less noticeable due to soil saturation and the urgency of vegetated surfaces to meet their physiological needs. In dry year, however, differences were more evident, with bare soil presenting lower ET rates and vegetation classes showing higher ET values.     

WorldView-2影像林地信息提取的研究与实现

林地信息的获取具有重要意义,可广泛应用于城市规划编制、林地资源分析、水土流失治理等领域。本文基于WorldView-2影像,在研究分割与对象特征的基础上,采用面向对象方法,实现了林地信息提取,其用户精度达到94%,制图精度达到97.2%。然后,将像元二分模型应用于林地类对象,计算其植被覆盖度。试验表明,面向对象方法能高效、准确、经济地提取林地信息;基于林地类对象计算的植被覆盖度能较好地表征植被覆盖情况。     

Application of machine learning techniques in groundwater potential mapping along the west coast of India

ABSTRACT

Groundwater potential mapping (GWPM) in the coastal zone is crucial for the planning and development of society and the environment. The current study is aimed to map the groundwater potential zones of Sindhudurg coastal stretch on the west coast of India, using three machine learning models: random forest (RF), boosted regression tree (BRT), and the ensemble of RF and support vector machine (SVM). In order to achieve the objective, 15 groundwater influencing factors including elevation, slope, aspect, slope length (LS), profile curvature, plan curvature, topographical wetness index (TWI), distance from streams, distance from lineaments, lithology, geomorphology, soil, land use, normalized difference vegetation index (NDVI), and rainfall were considered for inter-thematic correlations and overlaid with spring and well occurrences in a spatial database. A total of 165 spring and well locations were identified, which had been divided into two classes: training and validation, at the ratio of 70:30, respectively. The RF, BRT, and RF-SVM ensemble models have been applied to delineate the groundwater potential zones and categorized into five classes, namely very high, high, moderate, low, and very low. RF, BRT, and ensemble model results showed that 33.3%, 35.6%, and 36.8% of the research area had a very high groundwater potential zone. These models were validated with area under the receiver operating characteristics (AUROC) curve. The accuracy of RF (94%) and hybrid model (93.4%) was more efficient than BRT (89.8%) model. In order to further evaluate and validate, four different sites were subsequently chosen, and we obtained similar results, ensuring the validity of the applied models. Additionally, ground-penetrating radar (GPR) technique was applied to predict the groundwater table and validated by measured wells. The mean difference between measured and GPR predicted groundwater table was 14 cm, which reflected the importance of GPR to guide the location of new wells in the study region. The outcomes of the study will help the decision-makers, government agencies, and private sectors for sustainable planning of groundwater in the area. Overall, the present study provides a comprehensive high-precision machine learning and GPR-based groundwater potential mapping.