Temporal stability of soil moisture and radar backscatter observed by the advanced Synthetic Aperture Radar (ASAR)

The high spatio-temporal variability of soil moisture is the result of atmospheric forcing and redistribution processes related to terrain, soil, and vegetation characteristics. Despite this high variability, many field studies have shown that in the temporal domain soil moisture measured at specific locations is correlated to the mean soil moisture content over an area. Since the measurements taken by Synthetic Aperture Radar (SAR) instruments are very sensitive to soil moisture it is hypothesized that the temporally stable soil moisture patterns are reflected in the radar backscatter measurements. To verify this hypothesis 73 Wide Swath (WS) images have been acquired by the ENVISAT Advanced Synthetic Aperture Radar (ASAR) over the REMEDHUS soil moisture network located in the Duero basin, Spain. It is found that a time-invariant linear relationship is well suited for relating local scale (pixel) and regional scale (50 km) backscatter. The observed linear model coefficients can be estimated by considering the scattering properties of the terrain and vegetation and the soil moisture scaling properties. For both linear model coefficients, the relative error between observed and modelled values is less than 5 % and the coefficient of determination (R-2) is 86 %. The results are of relevance for interpreting and downscaling coarse resolution soil moisture data retrieved from active (METOP ASCAT) and passive (SMOS, AMSR-E) instruments

Soil erosion in the Alps : causes and risk assessment

The issue of soil erosion in the Alps has long been neglected due to the low economic value of the agricultural land. However, soil stability is a key parameter which affects ecosystem services like slope stability, water budgets (drinking water reservoirs as well as flood prevention), vegetation productivity, ecosystem biodiversity and nutrient production. In alpine regions, spatial estimates on soil erosion are difficult to derive because the highly heterogeneous biogeophysical structure impedes measurement of soil erosion and the applicability of soil erosion models. However, remote sensing and geographic information system (GIS) methods allow for spatial estimation of soil erosion by direct detection of erosion features and supply of input data for soil erosion models. Thus, the main objective of this work is to address the problem of soil erosion risk assessment in the Alps on catchment scale with remote sensing and GIS tools. Regarding soil erosion processes the focus is on soil erosion by water (here sheet erosion) and gravity (here landslides). For these two processes we address i) the monitoring and mapping of the erosion features and related causal factors ii) soil erosion risk assessment with special emphasis on iii) the validation of existing models for alpine areas. All investigations were accomplished in the Urseren Valley (Central Swiss Alps) where the valley slopes are dramatically affected by sheet erosion and landslides. For landslides, a natural susceptibility of the catchment has been indicated by bivariate and multivariate statistical analysis. Geology, slope and stream density are the most significant static landslide causal factors. Static factors are here defined as factors that do not change their attributes during the considered time span of the study (45 years), e.g. geology, stream network. The occurrence of landslides might be significantly increased by the combined effects of global climate and land use change. Thus, our hypothesis is that more recent changes in land use and climate affected the spatial and temporal occurrence of landslides. The increase of the landslide area of 92% within 45 years in the study site confirmed our hypothesis. In order to identify the cause for the trend in landslide occurrence time-series of landslide causal factors were analysed. The analysis revealed increasing trends in the frequency and intensity of extreme rainfall events and stocking of pasture animals. These developments presumably enhanced landslide hazard. Moreover, changes in land-cover and land use were shown to have affected landslide occurrence. For instance, abandoned areas and areas with recently emerging shrub vegetation show very low landslide densities. Detailed spatial analysis of the land use with GIS and interviews with farmers confirmed the strong influence of the land use management practises on slope stability. The definite identification and quantification of the impact of these non-stationary landslide causal factors (dynamic factors) on the landslide trend was not possible due to the simultaneous change of several factors. The consideration of dynamic factors in statistical landslide susceptibility assessments is still unsolved. The latter may lead to erroneous model predictions, especially in times of dramatic environmental change. Thus, we evaluated the effect of dynamic landslide causal factors on the validity of landslide susceptibility maps for spatial and temporal predictions. For this purpose, a logistic regression model based on data of the year 2000 was set up. The resulting landslide susceptibility map was valid for spatial predictions. However, the model failed to predict the landslides that occurred in a subsequent event. In order to handle this weakness of statistic landslide modelling a multitemporal approach was developed. It is based on establishing logistic regression models for two points in time (here 1959 and 2000). Both models could correctly classify >70% of the independent spatial validation dataset. By subtracting the 1959 susceptibility map from the 2000 susceptibility map a deviation susceptibility map was obtained. Our interpretation was that these susceptibility deviations indicate the effect of dynamic causal factors on the landslide probability. The deviation map explained 85% of new independent landslides occurring after 2000. Thus, we believe it to be a suitable tool to add a time element to a susceptibility map pointing to areas with changing susceptibility due to recently changing environmental conditions or human interactions. In contrast to landslides that are a direct threat to buildings and infrastructure, sheet erosion attracts less attention because it is often an unseen process. Nonetheless, sheet erosion may account for a major proportion of soil loss. Soil loss by sheet erosion is related to high spatial variability, however, in contrast to arable fields for alpine grasslands erosion damages are long lasting and visible over longer time periods. A crucial erosion triggering parameter that can be derived from satellite imagery is fractional vegetation cover (FVC). Measurements of the radiogenic isotope Cs-137, which is a common tracer for soil erosion, confirm the importance of FVC for soil erosion yield in alpine areas. Linear spectral unmixing (LSU), mixture tuned matched filtering (MTMF) and the spectral index NDVI are applied for estimating fractional abundance of vegetation and bare soil. To account for the small scale heterogeneity of the alpine landscape very high resolved multispectral QuickBird imagery is used. The performance of LSU and MTMF for estimating percent vegetation cover is good (r²=0.85, r²=0.71 respectively). A poorer performance is achieved for bare soil (r²=0.28, r²=0.39 respectively) because compared to vegetation, bare soil has a less characteristic spectral signature in the wavelength domain detected by the QuickBird sensor. Apart from monitoring erosion controlling factors, quantification of soil erosion by applying soil erosion risk models is done. The performance of the two established models Universal Soil Loss Equation (USLE) and Pan-European Soil Erosion Risk Assessment (PESERA) for their suitability to model erosion for mountain environments is tested. Cs-137 is used to verify the resulting erosion rates from USLE and PESERA. PESERA yields no correlation to measured Cs-137 long term erosion rates and shows lower sensitivity to FVC. Thus, USLE is used to model the entire study site. The LSU-derived FVC map is used to adapt the C factor of the USLE. Compared to the low erosion rates computed with the former available low resolution dataset (1:25000) the satellite supported USLE map shows “hotspots” of soil erosion of up to 16 t ha-1 a-1. In general, Cs-137 in combination with the USLE is a very suitable method to assess soil erosion for larger areas, as both give estimates on long-term soil erosion. Especially for inaccessible alpine areas, GIS and remote sensing proved to be powerful tools that can be used for repetitive measurements of erosion features and causal factors. In times of global change it is of crucial importance to account for temporal developments. However, the evaluation of the applied soil erosion risk models revealed that the implementation of temporal aspects, such as varying climate, land use and vegetation cover is still insufficient. Thus, the proposed validation strategies (spatial, temporal and via Cs-137) are essential. Further case studies in alpine regions are needed to test the methods elaborated for the Urseren Valley. However, the presented approaches are promising with respect to improve the monitoring and identification of soil erosion risk areas in alpine regions

Imaging multi-age construction settlement behaviour by advanced SAR interferometry

This paper focuses on the application of Advanced Satellite Synthetic Aperture Radar Interferometry (A-DInSAR) to subsidence-related issues, with particular reference to ground settlements due to external loads. Beyond the stratigraphic setting and the geotechnical properties of the subsoil, other relevant boundary conditions strongly influence the reliability of remotely sensed data for quantitative analyses and risk mitigation purposes. Because most of the Persistent Scatterer Interferometry (PSI) measurement points (Persistent Scatterers, PSs) lie on structures and infrastructures, the foundation type and the age of a construction are key factors for a proper interpretation of the time series of ground displacements. To exemplify a methodological approach to evaluate these issues, this paper refers to an analysis carried out in the coastal/deltaic plain west of Rome (Rome and Fiumicino municipalities) affected by subsidence and related damages to structures. This region is characterized by a complex geological setting (alternation of recent deposits with low and high compressibilities) and has been subjected to different urbanisation phases starting in the late 1800s, with a strong acceleration in the last few decades. The results of A-DInSAR analyses conducted from 1992 to 2015 have been interpreted in light of high-resolution geological/geotechnical models, the age of the construction, and the types of foundations of the buildings on which the PSs are located. Collection, interpretation, and processing of geo-thematic data were fundamental to obtain high-resolution models; change detection analyses of the land cover allowed us to classify structures/infrastructures in terms of the construction period. Additional information was collected to define the types of foundations, i.e., shallow versus deep foundations. As a result, we found that only by filtering and partitioning the A-DInSAR datasets on the basis of the above-mentioned boundary conditions can the related time series be considered a proxy of the consolidation process governing the subsidence related to external loads as confirmed by a comparison with results from a physically based back analysis based on Terzaghi's theory. Therefore, if properly managed, the A-DInSAR data represents a powerful tool for capturing the evolutionary stage of the process for a single building and has potential for forecasting the behaviour of the terrain-foundation-structure combination

From imagery to ecology: Leveraging time series of all available Landsat observations to map and monitor ecosystem state and dynamics

https://zslpublications.onlinelibrary.wiley.com/doi/full/10.1002/rse2.24Published versio

A multi-temporal phenology based classification approach for Crop Monitoring in Kenya

The SBAM (Satellite Based Agricultural Monitoring) project, funded by the Italian Space Agency aims at: developing a validated satellite imagery based method for estimating and updating the agricultural areas in the region of Central-Africa; implementing an automated process chain capable of providing periodical agricultural land cover maps of the area of interest and, possibly, an estimate of the crop yield. The project aims at filling the gap existing in the availability of high spatial resolution maps of the agricultural areas of Kenya. A high spatial resolution land cover map of Central-Eastern Africa including Kenya was compiled in the year 2000 in the framework of the Africover project using Landsat images acquired, mostly, in 1995. We investigated the use of phenological information in supporting the use of remotely sensed images for crop classification and monitoring based on Landsat 8 and, in the near future, Sentinel 2 imagery. Phenological information on crop condition was collected using time series of NDVI (Normalized Difference Vegetation Index) based on Landsat 8 images. Kenyan countryside is mainly characterized by a high number of fragmented small and medium size farmlands that dramatically increase the difficulty in classification; 30 m spatial resolution images are not enough for a proper classification of such areas. So, a pan-sharpening FIHS (Fast Intensity Hue Saturation) technique was implemented to increase image resolution from 30 m to 15 m. Ground test sites were selected, searching for agricultural vegetated areas from which phenological information was extracted. Therefore, the classification of agricultural areas is based on crop phenology, vegetation index behaviour retrieved from a time series of satellite images and on AEZ (Agro Ecological Zones) information made available by FAO (FAO, 1996) for the area of interest. This paper presents the results of the proposed classification procedure in comparison with land cover maps produced in the past years by other projects. The results refer to the Nakuru County and they were validated using field campaigns data. It showed a satisfactory overall accuracy of 92.66 % which is a significant improvement with respect to previous land cover maps

Urban energy exchanges monitoring from space

One important challenge facing the urbanization and global environmental change community is to understand the relation between urban form, energy use and carbon emissions. Missing from the current literature are scientific assessments that evaluate the impacts of different urban spatial units on energy fluxes; yet, this type of analysis is needed by urban planners, who recognize that local scale zoning affects energy consumption and local climate. However, satellite-based estimation of urban energy fluxes at neighbourhood scale is still a challenge. Here we show the potential of the current satellite missions to retrieve urban energy budget, supported by meteorological observations and evaluated by direct flux measurements. We found an agreement within 5% between satellite and in-situ derived net all-wave radiation; and identified that wall facet fraction and urban materials type are the most important parameters for estimating heat storage of the urban canopy. The satellite approaches were found to underestimate measured turbulent heat fluxes, with sensible heat flux being most sensitive to surface temperature variation (-64.1, +69.3 W m-2 for ±2 K perturbation); and also underestimate anthropogenic heat flux. However, reasonable spatial patterns are obtained for the latter allowing hot-spots to be identified, therefore supporting both urban planning and urban climate modelling

Evaluation of a global soil moisture product from finer spatial resolution sar data and ground measurements at Irish sites

In the framework of the European Space Agency Climate Change Initiative, a global, almost daily, soil moisture (SM) product is being developed from passive and active satellite microwave sensors, at a coarse spatial resolution. This study contributes to its validation by using finer spatial resolution ASAR Wide Swath and in situ soil moisture data taken over three sites in Ireland, from 2007 to 2009. This is the first time a comparison has been carried out between three sets of independent observations from different sensors at very different spatial resolutions for such a long time series. Furthermore, the SM spatial distribution has been investigated at the ASAR scale within each Essential Climate Variable (ECV) pixel, without adopting any particular model or using a densely distributed network of in situ stations. This approach facilitated an understanding of the extent to which geophysical factors, such as soil texture, terrain composition and altitude, affect the retrieved ECV SM product values in temperate grasslands. Temporal and spatial variability analysis provided high levels of correlation (p < 0.025) and low errors between the three datasets, leading to confidence in the new ECV SM global product, despite limitations in its ability to track the driest and wettest conditions

Satellite based synthetic aperture radar and optical spatial-temporal information as aid for operational and environmental mine monitoring

A sustainable society is a society that satisfies its resource requirements without endangering the sustainability of these resources. The mineral endowment on the African continent is estimated to be the first or second largest of world reserves. Therefore, it is recognised that the African continent still heavily depends on mineral exports as a key contributor to the gross domestic product (GDP) of various countries. These mining activities, however, do introduce primary and secondary environmental degradation factors. They attract communities to these mining areas, light and heavy industrial establishments occur, giving rise to artisanal activities. This study focussed on satellite RS products as an aid to a mine’s operations and the monitoring of its environment. Effective operational mine management and control ensures a more sustainable and profitable lifecycle for mines. Satellite based RS holds the potential to observe the mine and its surrounding areas at high temporal intervals, different spectral wavelengths and spatial resolutions. The combination of SAR and optical information creates a spatial platform to observe and measure the mine’s operations and the behaviour of specific land cover and land use classes over time and contributes to a better understanding of the mining activities and their influence on the environment within a specific geographical area. This study will introduce an integrated methodology to collect, process and analyse spatial information over a specific targeted mine. This methodology utilises a medium resolution land cover base map, derived from Landsat 8, to understand the predominant land cover types of the surrounding area. Using very high resolution mono- and stereoscopic satellite imagery provides a finer scale analysis and identifies changes in features at a smaller scale. Combining these technologies with the synthetic aperture radar (SAR) applications for precise measurement of surface subsidence or upliftment becomes a spatial toolbox for mine management. This study examines a combination of satellite remote sensing products guided by a systematic workflow methodology to integrate spatial results as an aid for mining operations and environmental monitoring. Some of the results that can be highlighted is the successful land cover classification using the Landsat 8 satellite. The land cover that dominated the Kolomela mine area was the “SHRUBLAND/GRASS” class with a 94% coverage and “MINE” class of 2.6%. Sishen mine had a similar dominated land cover characteristic with a “SHRUBLAND/GRASS” class of 90% and “MINE” class of 4.8%. The Pléiades time-series classification analysis was done using three scenes each acquired at a different time interval. The Sishen and Kolomela mine showed especially changes from the bare soil class to the asphalt or mine class. The Pléiades stereoscopic analysis provided volumetric change detection over small, medium, large and recessed areas. Both the Sishen and Kolomela mines demonstrated height profile changes in each selected category. The last category of results focused on the SAR technology to measure within millimetre accuracy the subsidence and upliftment behaviour of surface areas over time. The Royal Bafokeng Platinum tailings pond area was measured using 74 TerraSAR-X scenes. The tailings wall area was confirmed as stable with natural subsidence that occurred in its surrounding area due to seasonal changes of the soil during rainy and dry periods. The Chuquicamata mine as a large open pit copper mine area was analysed using 52 TerraSAR-X scenes. The analysis demonstrated significant vertical surface movement over some of the dumping sites. It is the wish of the researcher that this dissertation and future research scholars will continue to contribute in this scientific field. These contributions can only assist the mining sector to continuously improve its mining operations as well as its monitoring of the primary as well as the secondary environmental impacts to ensure improved sustainability for the next generation.Environmental SciencesM. Sc. (Environmental Science