The data concept behind the data: From metadata models and labelling schemes towards a generic spectral library

Spectral libraries play a major role in imaging spectroscopy. They are commonly used to store end-member and spectrally pure material spectra, which are primarily used for mapping or unmixing purposes. However, the development of spectral libraries is time consuming and usually sensor and site dependent. Spectral libraries are therefore often developed, used and tailored only for a specific case study and only for one sensor. Multi-sensor and multi-site use of spectral libraries is difficult and requires technical effort for adaptation, transformation, and data harmonization steps. Especially the huge amount of urban material specifications and its spectral variations hamper the setup of a complete spectral library consisting of all available urban material spectra. By a combined use of different urban spectral libraries, besides the improvement of spectral inter- and intra-class variability, missing material spectra could be considered with respect to a multi-sensor/ -site use. Publicly available spectral libraries mostly lack the metadata information that is essential for describing spectra acquisition and sampling background, and can serve to some extent as a measure of quality and reliability of the spectra and the entire library itself. In the GenLib project, a concept for a generic, multi-site and multi-sensor usable spectral library for image spectra on the urban focus was developed. This presentation will introduce a 1) unified, easy-to-understand hierarchical labeling scheme combined with 2) a comprehensive metadata concept that is 3) implemented in the SPECCHIO spectral information system to promote the setup and usability of a generic urban spectral library (GUSL). The labelling scheme was developed to ensure the translation of individual spectral libraries with their own labelling schemes and their usually varying level of details into the GUSL framework. It is based on a modified version of the EAGLE classification concept by combining land use, land cover, land characteristics and spectral characteristics. The metadata concept consists of 59 mandatory and optional attributes that are intended to specify the spatial context, spectral library information, references, accessibility, calibration, preprocessing steps, and spectra specific information describing library spectra implemented in the GUSL. It was developed on the basis of existing metadata concepts and was subject of an expert survey. The metadata concept and the labelling scheme are implemented in the spectral information system SPECCHIO, which is used for sharing and holding GUSL spectra. It allows easy implementation of spectra as well as their specification with the proposed metadata information to extend the GUSL. Therefore, the proposed data model represents a first fundamental step towards a generic usable and continuously expandable spectral library for urban areas. The metadata concept and the labelling scheme also build the basis for the necessary adaptation and transformation steps of the GUSL in order to use it entirely or in excerpts for further multi-site and multi-sensor applications

Spectroscopy-supported digital soil mapping

Global environmental changes have resulted in changes in key ecosystem services that soils provide. It is necessary to have up to date soil information on regional and global scales to ensure that these services continue to be provided. As a result, Digital Soil Mapping (DSM) research priorities are among others, advancing methods for data collection and analyses tailored towards large-scale mapping of soil properties. Scientifically, this thesis contributed to the development of methodologies, which aim to optimally use remote and proximal sensing (RS and PS) for DSM to facilitate regional soil mapping. The main contributions of this work with respect to the latter are (I) the critical evaluation of recent research achievements and identification of knowledge gaps for large-scale DSM using RS and PS data, (II) the development of a sparse RS-based sampling approach to represent major soil variability at regional scale, (III) the evaluation and development of different state-of-the-art methods to retrieve soil mineral information from PS, (IV) the improvement of spatially explicit soil prediction models and (V) the integration of RS and PS methods with geostatistical and DSM methods. A review on existing literature about the use of RS and PS for soil and terrain mapping was presented in Chapter 2. Recent work indicated the large potential of using RS and PS methods for DSM. However, for large-scale mapping, current methods will need to be extended beyond the plot. Improvements may be expected in the fields of developing more quantitative methods, enhanced geostatistical analysis and improved transferability to other areas. From these findings, three major research interests were selected: (I) soil sampling strategies, (II) retrieval of soil information from PS and (III) spatially continuous mapping of soil properties at larger scales using RS. Budgetary constraints, limited time and available soil legacy data restricted the soil data acquisition, presented in Chapter 3. A 15.000 km2 area located in Northern Morocco served as test case. Here, a sample was collected using constrained Latin Hypercube Sampling (cLHS) of RS and elevation data. The RS data served as proxy for soil variability, as alternative for the required soil legacy data supporting the sampling strategy. The sampling aim was to optimally sample the variability in the RS data while minimizing the acquisition efforts. This sample resulted in a dataset representing major soil variability. The cLHS sample failed to express spatial correlation; constraining the LHS by a distance criterion favoured large spatial variability over short distances. The absence of spatial correlation in the sampled soil variability precludes the use of additional geostatistical analyses to spatially predict soil properties. Predicting soil properties using the cLHS sample is thus restricted to a modelled statistical relation between the sample and exhaustive predictor variables. For this, the RS data provided the necessary spatial information because of the strong spatial correlation while the spectral information provided the variability of the environment (Chapter 3 and 6). Concluding, the RS-based cLHS approach is considered a time and cost efficient method for acquiring information on soil resources over extended areas. This sample was further used for developing methods to derive soil mineral information from PS, and to characterize regional soil mineralogy using RS. In Chapter 4, the influences of complex scattering within the mixture and overlapping absorption features were investigated. This was done by comparing the success of PRISM’s MICA in determining mineralogy of natural samples and modelled spectra. The modelled spectra were developed by a linearly forward model of reflectance spectra, using the fraction of known constituents within the sample. The modelled spectra accounted for the co-occurrence of absorption features but eluded the complex interaction between the components. It was found that more minerals could be determined with higher accuracy using modelled reflectance. The absorption features in the natural samples were less distinct or even absent, which hampered the classification routine. Nevertheless, grouping the individual minerals into mineral categories significantly improved the classification accuracy. These mineral categories are particularly useful for regional scale studies, as key soil property for parent material characterization and soil formation. Characterizing regional soil mineralogy by mineral categories was further described in Chapter 6. Retrieval of refined information from natural samples, such as mineral abundances, is more complex; estimating abundances requires a method that accounts for the interaction between minerals within the intimate mixture. This can be done by addressing the interaction with a non-linear model (Chapter 5). Chapter 5 showed that mineral abundances in complex mixtures could be estimated using absorption features in the 2.1–2.4 µm wavelength region. First, the absorption behaviour of mineral mixtures was parameterized by exponential Gaussian optimization (EGO). Next, mineral abundances were successfully predicted by regression tree analysis, using these parameters as inputs. Estimating mineral abundances using prepared mixes of calcite, kaolinite, montmorillonite and dioctahedral mica or field samples proved the validity of the proposed method. Estimating mineral abundances of field samples showed the necessity to deconvolve spectra by EGO. Due to the nature of the field samples, the simple representation of the complex scattering behaviour by a few Gaussian bands required the parameters asymmetry and saturation to accurately deconvolve the spectra. Also, asymmetry of the EGO profiles showed to be an important parameter for estimating the abundances of the field samples. The robustness of the method in handling the omission of minerals during the training phase was tested by replacing part of the quartz with chlorite. It was found that the accuracy of the predicted mineral content was hardly affected. Concluding, the proposed method allowed for estimating more than two minerals within a mixture. This approach advances existing PS methods and has the potential to quantify a wider set of soil properties. With this method the soil science community was provided an improved inference method to derive and quantify soil properties The final challenge of this thesis was to spatially explicit model regional soil mineralogy using the sparse sample from Chapter 3. Prediction models have especially difficulties relating predictor variables to sampled properties having high spatial correlation. Chapter 6 presented a methodology that improved prediction models by using scale-dependent spatial variability observed in RS data. Mineral predictions were made using the abundances from X-ray diffraction analysis and mineral categories determined by PRISM. The models indicated that using the original RS data resulted in lower model performance than those models using scaled RS data. Key to the improved predictions was representing the variability of the RS data at the same scale as the sampled soil variability. This was realized by considering the medium and long-range spatial variability in the RS data. Using Fixed Rank Kriging allowed smoothing the massive RS datasets to these ranges. The resulting images resembled more closely the regional spatial variability of soil and environmental properties. Further improvements resulted from using multi-scale soil-landscape relationships to predict mineralogy. The maps of predicted mineralogy showed agreement between the mineral categories and abundances. Using a geostatistical approach in combination with a small sample, substantially improves the feasibility to quantitatively map regional mineralogy. Moreover, the spectroscopic method appeared sufficiently detailed to map major mineral variability. Finally, this approach has the potential for modelling various natural resources and thereby enhances the perspective of a global system for inventorying and monitoring the earth’s soil resources. With this thesis it is demonstrated that RS and PS methods are an important but also an essential source for regional-scale DSM. Following the main findings from this thesis, it can be concluded that: Improvements in regional-scale DSM result from the integrated use of RS and PS with geostatistical methods. In every step of the soil mapping process, spectroscopy can play a key role and can deliver data in a time and cost efficient manner. Nevertheless, there are issues that need to be resolved in the near future. Research priorities involve the development of operational tools to quantify soil properties, sensor integration, spatiotemporal modelling and the use of geostatistical methods that allow working with massive RS datasets. This will allow us in the near future to deliver more accurate and comprehensive information about soils, soil resources and ecosystem services provided by soils at regional and, ultimately, global scale.</p

Hyperspectral Imaging for Fine to Medium Scale Applications in Environmental Sciences

The aim of the Special Issue “Hyperspectral Imaging for Fine to Medium Scale Applications in Environmental Sciences” was to present a selection of innovative studies using hyperspectral imaging (HSI) in different thematic fields. This intention reflects the technical developments in the last three decades, which have brought the capacity of HSI to provide spectrally, spatially and temporally detailed data, favoured by e.g., hyperspectral snapshot technologies, miniaturized hyperspectral sensors and hyperspectral microscopy imaging. The present book comprises a suite of papers in various fields of environmental sciences—geology/mineral exploration, digital soil mapping, mapping and characterization of vegetation, and sensing of water bodies (including under-ice and underwater applications). In addition, there are two rather methodically/technically-oriented contributions dealing with the optimized processing of UAV data and on the design and test of a multi-channel optical receiver for ground-based applications. All in all, this compilation documents that HSI is a multi-faceted research topic and will remain so in the future

Combining Field and Imaging Spectroscopy to Map Soil Organic Carbon in a Semiarid Environment

Semiarid regions are especially vulnerable to climate change and human-induced land-use changes and are of major importance in the context of necessary carbon sequestration and ongoing land degradation. Topsoil properties, such as soil carbon content, provide valuable indicators to these processes, and can be mapped using imaging spectroscopy (IS). In semiarid regions, this poses difficulties because models are needed that can cope with varying land surface and soil conditions, consider a partial vegetation coverage, and deal with usually low soil organic carbon (SOC) contents. We present an approach that aims at addressing these difficulties by using a combination of field and IS to map SOC in an extensively used semiarid ecosystem. In hyperspectral imagery of the HyMap sensor, the influence of nonsoil materials, i.e., vegetation, on the spectral signature of soil dominated image pixels was reduced and a residual soil signature was calculated. The proposed approach allowed this procedure up to a vegetation coverage of 40% clearly extending the mapping capability. SOC quantities are predicted by applying a spectral feature-based SOC prediction model to image data of residual soil spectra. With this approach, we could significantly increase the spatial extent for which SOC could be predicted with a minimal influence of a vegetation signal compared to previous approaches where the considered area was limited to a maximum of, e.g., 10% vegetation coverage. As a regional example, the approach was applied to a 320 km2 area in the Albany Thicket Biome, South Africa, where land cover and landuse changes have occurred due to decades of unsustainable land management. In the generated maps, spatial SOC patterns were interpreted and linked to geomorphic features and land surface processes, i.e., areas of soil erosion. It was found that the chosen approach supported the extraction of soil-related spectral image information in the semiarid region with highly varying land cover. However, the quantitative prediction of SOC contents revealed a lack in absolute accuracy

Quantitative Mapping of Soil Property Based on Laboratory and Airborne Hyperspectral Data Using Machine Learning

Soil visible and near-infrared spectroscopy provides a non-destructive, rapid and low-cost approach to quantify various soil physical and chemical properties based on their reflectance in the spectral range of 400–2500 nm. With an increasing number of large-scale soil spectral libraries established across the world and new space-borne hyperspectral sensors, there is a need to explore methods to extract informative features from reflectance spectra and produce accurate soil spectroscopic models using machine learning. Features generated from regional or large-scale soil spectral data play a key role in the quantitative spectroscopic model for soil properties. The Land Use/Land Cover Area Frame Survey (LUCAS) soil library was used to explore PLS-derived components and fractal features generated from soil spectra in this study. The gradient-boosting method performed well when coupled with extracted features on the estimation of several soil properties. Transfer learning based on convolutional neural networks (CNNs) was proposed to make the model developed from laboratory data transferable for airborne hyperspectral data. The soil clay map was successfully derived using HyMap imagery and the fine-tuned CNN model developed from LUCAS mineral soils, as deep learning has the potential to learn transferable features that generalise from the source domain to target domain. The external environmental factors like the presence of vegetation restrain the application of imaging spectroscopy. The reflectance data can be transformed into a vegetation suppressed domain with a force invariance approach, the performance of which was evaluated in an agricultural area using CASI airborne hyperspectral data. However, the relationship between vegetation and acquired spectra is complicated, and more efforts should put on removing the effects of external factors to make the model transferable from one sensor to another.:Abstract I Kurzfassung III Table of Contents V List of Figures IX List of Tables XIII List of Abbreviations XV 1 Introduction 1 1.1 Motivation 1 1.2 Soil spectra from different platforms 2 1.3 Soil property quantification using spectral data 4 1.4 Feature representation of soil spectra 5 1.5 Objectives 6 1.6 Thesis structure 7 2 Combining Partial Least Squares and the Gradient-Boosting Method for Soil Property Retrieval Using Visible Near-Infrared Shortwave Infrared Spectra 9 2.1 Abstract 10 2.2 Introduction 10 2.3 Materials and methods 13 2.3.1 The LUCAS soil spectral library 13 2.3.2 Partial least squares algorithm 15 2.3.3 Gradient-Boosted Decision Trees 15 2.3.4 Calculation of relative variable importance 16 2.3.5 Assessment 17 2.4 Results 17 2.4.1 Overview of the spectral measurement 17 2.4.2 Results of PLS regression for the estimation of soil properties 19 2.4.3 Results of PLS-GBDT for the estimation of soil properties 21 2.4.4 Relative important variables derived from PLS regression and the gradient-boosting method 24 2.5 Discussion 28 2.5.1 Dimension reduction for high-dimensional soil spectra 28 2.5.2 GBDT for quantitative soil spectroscopic modelling 29 2.6 Conclusions 30 3 Quantitative Retrieval of Organic Soil Properties from Visible Near-Infrared Shortwave Infrared Spectroscopy Using Fractal-Based Feature Extraction 31 3.1 Abstract 32 3.2 Introduction 32 3.3 Materials and Methods 35 3.3.1 The LUCAS topsoil dataset 35 3.3.2 Fractal feature extraction method 37 3.3.3 Gradient-boosting regression model 37 3.3.4 Evaluation 41 3.4 Results 42 3.4.1 Fractal features for soil spectroscopy 42 3.4.2 Effects of different step and window size on extracted fractal features 45 3.4.3 Modelling soil properties with fractal features 47 3.4.3 Comparison with PLS regression 49 3.5 Discussion 51 3.5.1 The importance of fractal dimension for soil spectra 51 3.5.2 Modelling soil properties with fractal features 52 3.6 Conclusions 53 4 Transfer Learning for Soil Spectroscopy Based on Convolutional Neural Networks and Its Application in Soil Clay Content Mapping Using Hyperspectral Imagery 55 4.1 Abstract 55 4.2 Introduction 56 4.3 Materials and Methods 59 4.3.1 Datasets 59 4.3.2 Methods 62 4.3.3 Assessment 67 4.4 Results and Discussion 67 4.4.1 Interpretation of mineral and organic soils from LUCAS dataset 67 4.4.2 1D-CNN and spectral index for LUCAS soil clay content estimation 69 4.4.3 Application of transfer learning for soil clay content mapping using the pre-trained 1D-CNN model 72 4.4.4 Comparison between spectral index and transfer learning 74 4.4.5 Large-scale soil spectral library for digital soil mapping at the local scale using hyperspectral imagery 75 4.5 Conclusions 75 5 A Case Study of Forced Invariance Approach for Soil Salinity Estimation in Vegetation-Covered Terrain Using Airborne Hyperspectral Imagery 77 5.1 Abstract 78 5.2 Introduction 78 5.3 Materials and Methods 81 5.3.1 Study area of Zhangye Oasis 81 5.3.2 Data description 82 5.3.3 Methods 83 5.3.3 Model performance assessment 85 5.4 Results and Discussion 86 5.4.1 The correlation between NDVI and soil salinity 86 5.4.2 Vegetation suppression performance using the Forced Invariance Approach 86 5.4.3 Estimation of soil properties using airborne hyperspectral data 88 5.5 Conclusions 90 6 Conclusions and Outlook 93 Bibliography 97 Acknowledgements 11

Caracterização e estudo comparativo de exsudações de hidrocarbonetos e plays petrolíferos em bacias terrestres das regiões central do Irã e sudeste do Brasil usando sensoriamento remoto espectral

Orientador: Carlos Roberto de Souza FilhoTese (doutorado) - Universidade Estadual de Campinas, Instituto de GeociênciasResumo: O objetivo desta pesquisa foi explorar as assinaturas de exsudações de hidrocarbonetos na superfície usando a tecnologia de detecção remota espectral. Isso foi alcançado primeiro, realizando uma revisão abrangente das capacidades e potenciais técnicas de detecção direta e indireta. Em seguida, a técnica foi aplicada para investigar dois locais de teste localizados no Irã e no Brasil, conhecidos por hospedar sistemas ativos de micro-exsudações e afloramentos betuminosos, respectivamente. A primeira área de estudo está localizada perto da cidade de Qom (Irã), e está inserida no campo petrolífero Alborz, enterrado sob sedimentos datados do Oligoceno da Formação Upper Red. O segundo local está localizado perto da cidade de Anhembi (SP), na margem oriental da bacia do Paraná, no Brasil, e inclui acumulações de betume em arenitos triássicos da Formação Pirambóia. O trabalho na área de Qom integrou evidências de (i) estudos petrográficos e geoquímicos em laboratório, (ii) investigações de afloramentos em campo, e (iii) mapeamento de anomalia em larga escala através de conjuntos de dados multi-espectrais ASTER e Sentinel-2. O resultado deste estudo se trata de novos indicadores mineralógicos e geoquímicos para a exploração de micro-exsudações e um modelo de micro-exsudações atualizado. Durante este trabalho, conseguimos desenvolver novas metodologias para análise de dados espectroscópicos. Através da utilização de dados simulados, indicamos que o instrumento de satélite WorldView-3 tem potencial para detecção direta de hidrocarbonetos. Na sequência do estudo, dados reais sobre afloramentos de arenitos e óleo na área de Anhembi foram investigados. A área foi fotografada novamente no chão e usando o sistema de imagem hiperespectral AisaFENIX. Seguiu-se estudos e amostragem no campo,incluindo espectroscopia de alcance fechado das amostras no laboratório usando instrumentos de imagem (ou seja, sisuCHEMA) e não-imagem (ou seja, FieldSpec-4). O estudo demonstrou que uma abordagem espectroscópica multi-escala poderia fornecer uma imagem completa das variações no conteúdo e composição do betume e minerais de alteração que acompanham. A assinatura de hidrocarbonetos, especialmente a centrada em 2300 nm, mostrou-se consistente e comparável entre as escalas e capaz de estimar o teor de betume de areias de petróleo em todas as escalas de imagemAbstract: The objective of this research was to explore for the signatures of seeping hydrocarbons on the surface using spectral remote sensing technology. It was achieved firstly by conducting a comprehensive review of the capacities and potentials of the technique for direct and indirect seepage detection. Next, the technique was applied to investigate two distinctive test sites located in Iran and Brazil known to retain active microseepage systems and bituminous outcrops, respectively. The first study area is located near the city of Qom in Iran, and consists of Alborz oilfield buried under Oligocene sediments of the Upper-Red Formation. The second site is located near the town of Anhembi on the eastern edge of the Paraná Basin in Brazil and includes bitumen accumulations in the Triassic sandstones of the Pirambóia Formation. Our work in Qom area integrated evidence from (i) petrographic, spectroscopic, and geochemical studies in the laboratory, (ii) outcrop investigations in the field, and (iii) broad-scale anomaly mapping via orbital remote sensing data. The outcomes of this study was novel mineralogical and geochemical indicators for microseepage characterization and a classification scheme for the microseepage-induced alterations. Our study indicated that active microseepage systems occur in large parts of the lithofacies in Qom area, implying that the extent of the petroleum reservoir is much larger than previously thought. During this work, we also developed new methodologies for spectroscopic data analysis and processing. On the other side, by using simulated data, we indicated that WorldView-3 satellite instrument has the potential for direct hydrocarbon detection. Following this demonstration, real datasets were acquired over oil-sand outcrops of the Anhembi area. The area was further imaged on the ground and from the air by using an AisaFENIX hyperspectral imaging system. This was followed by outcrop studies and sampling in the field and close-range spectroscopy in the laboratory using both imaging (i.e. sisuCHEMA) and nonimaging instruments. The study demonstrated that a multi-scale spectroscopic approach could provide a complete picture of the variations in the content and composition of bitumen and associated alteration mineralogy. The oil signature, especially the one centered at 2300 nm, was shown to be consistent and comparable among scales, and capable of estimating the bitumen content of oil-sands at all imaging scalesDoutoradoGeologia e Recursos NaturaisDoutor em Geociências2015/06663-7FAPES

Mapping landscape function with hyperspectral remote sensing of natural grasslands on gold mines

Thesis submitted in fulfilment of the requirements for the degree of Doctor of Philosophy. School of Animal, Plant and Environmental Science, University of the Witwatersrand, Johannesburg, South Africa. October 2016.Mining has negative impacts on the environment in many different ways. One method developed to quantify some of these impacts is Landscape Function Analysis (LFA) and this has been accepted by some mining companies and regulators. In brief, LFA aims at quantifying the organization of vegetative and landscape components in a landscape into patches along a transect and quantifying, in a relative manner, three basic processes important to landscape functioning, namely: soil stability or susceptibility to erosion, infiltration or runoff, and nutrient cycling or organic matter decomposition. However, LFA is limited in large heterogeneous environments, such as those around mining operations, due to its localized nature, and the man hours required to collect a representative set of measurements for such large and complex environments. Remote sensing using satellite-acquired data can overcome these limitations by sampling the entire environment in a rapid and objective manner. What is required is a method of connecting these satellite-based measurements to LFA measurements and then being able to extrapolate these measurements across the entire mine surface. The aim of this research was to develop a method to use satellite-based hyperspectral imagery to predict landscape function analysis (LFA) using partial least squares regression (PLSR). This was broken down into three objectives: (1) Collection of the LFA data in the field and validation of the LFA indices against other environmental variables collected at the same time, (2) validation of PLSR models predicting LFA indices and various environmental variables from ground-based spectra, and (3) production of risk maps based on predicting LFA indices and above-ground biomass using PLSR models and Hyperion satellite-based hyperspectral imagery. Although the study was based in grasslands at two mining regions, West Wits and Vaal River, a suitable Hyperion image was only available for Vaal River. A minimum of 374 points were sampled for LFA indices, ground-based spectra, above-ground biomass and soil cores along 2880 m of LFA transect from both mine sites. Soil cores were weighed fresh before sieving with a 2 mm sieve to separate root and stone fractions. The sieved soil fraction was tested for pH, EC, SOM, and for the West Wits samples, organic nitrogen and total extractable inorganic nitrogen. There was one modification to the LFA method where grass patches were collapsed into homogenous units as it was deemed not feasible to sample 180 m transects at grass tuft scales of 10 – 30 cm, but other patch definitions followed the LFA manual (Tongway and Hindley, 2004). Evidence suggested that some of the different patch types, in particular the bare/biological soil crust – bare grass – sparse grass patch types, represented successional stages in a continuum although this was not conclusive. There also was evidence that the presence or absence of cattle play a role in some processes active in these grasslands and erosion is mainly through deflation, rain splash and sheet wash. Generally the environmental variables supported the LFA indices although the nutrient cycling index was representative of above-ground nutrient cycling but not below-ground nutrient cycling. Models derived with PLSR to predict the LFA indices from ground-based spectral measurements were strong at both mine sites (West Wits: LFA stability r2 = 0.63, P < 0.0001; LFA infiltration r2 = 0.75, P < 0.0001; LFA nutrient cycling r2 = 0.73, P < 0.0001; Vaal River: LFA stability r2 = 0.39, P < 0.0001, LFA infiltration r2 = 0.72, P < 0.0001, LFA nutrient cycling r2 = 0.54, P < 0.0001), as were PLSR models predicting above-ground biomass (West Wits above-ground biomass r2 = 0.55, P = 0.0003; Vaal River above-ground biomass r2 = 0.79, P < 0.0001) and soil moisture (West Wits soil moisture r2 = 0.45, P = 0.0017; Vaal River soil moisture r2 = 0.68, P < 0.0001). However, for soil organic matter (r2 = 0.50, P < 0.0001) and EC (r2 = 0.63, P < 0.0001), Vaal River had strong prediction models while West Wits had weak models for these variables (r2 = 0.31, P = 0.019 and r2 = 0.10 and P < 0.18, respectively). For EC, the wide range of soil values at Vaal River in association with gypsum crusts, and low values throughout West Wits explained these model results but for soil organic matter, no clear explanation for these site differences was identified. Patch-based models could accurately discriminate between spectrally well-defined patch types such S. plumosum patches but were less successful with patch types that were spectrally similar such as the bare/biological soil crust – bare grass – sparse grass patch continuum. Clustering similar patch types together before PLSR modelling did improve these patch-based spectral models. To test the method proposed to predict LFA indices from satellite-based hyperspectral imagery, a Hyperion image matching 6 transects at Vaal River was acquired by NASA’s EO-1 satellite and downloaded from the USGS Glovis website. LFA transects were partitioned to match and extract pixel spectra from the Hyperion data cube. Thirty-one spectra were separated into calibration (20) and validation (11) data. PLSR models were derived from the calibration data, tested with validation data to select the optimum model, and then applied to the entire Hyperion data cube to produce prediction maps for five LFA indices and above-ground biomass. The patch area index (PAI) produced particularly strong models (r2 = 0.79, P = 0.0003, n =11) with validation data, whereas the landscape organization index (LOI) produced weak models. It is argued that this difference between these two essentially similar indices is related to the fact that the PAI is a 2-dimensional index and the LOI is a 1-dimensional index. This difference in these two indices allowed the PAI to compensate for some burned pixels on the transects by “seeing” the density pattern of grass tufts and patches whereas the linear nature of the LOI was more susceptible to the changing dimensions of patch structure due to the effects of fire. Although validation models for the three LFA indices of soil stability, infiltration and nutrient cycling were strong (r2 = 0.72, P = 0.004; r2 = 0.66, P = 0.008; r2 = 0.70, P = 0.005, n = 9 respectively), prediction maps were confounded by the presence of fire on some transects. The poor quality of the Hyperion imagery also meant great care had to be taken in the selection of models to avoid poor quality prediction maps. The 31 bands from the VNIR (478 – 885 nm) portion of the Hyperion spectra were generally the best for PLSR modelling and prediction maps, presumably because of better signal-to-noise ratios due to higher energy in the shorter wavelengths. With two satellite-based hyperspectral sensors already operational, namely the US Hyperion and the Chinese HJ-1A HSI, and a number expected to be launched by various space agencies in the next few years, this research presents a method to use the strengths of LFA and hyperspectral imagery to model and predict LFA index values and thereby produce risk maps of large, heterogeneous landscapes such as mining environments. As this research documents a method of partitioning the landscape rather than the pixel spectra into pure endmembers, it makes a valuable contribution to the fields of landscape ecology and hyperspectral remote sensing.LG201

Satellite Imagery to Map Topsoil Organic Carbon Content over Cultivated Areas: An Overview

There is a need to update soil maps and monitor soil organic carbon (SOC) in the upper horizons or plough layer for enabling decision support and land management, while complying with several policies, especially those favoring soil carbon storage. This review paper is dedicated to the satellite-based spectral approaches for SOC assessment that have been achieved from several satellite sensors, study scales and geographical contexts in the past decade. Most approaches relying on pure spectral models have been carried out since 2019 and have dealt with temperate croplands in Europe, China and North America at the scale of small regions, of some hundreds of km(2): dry combustion and wet oxidation were the analytical determination methods used for 50% and 35% of the satellite-derived SOC studies, for which measured topsoil SOC contents mainly referred to mineral soils, typically cambisols and luvisols and to a lesser extent, regosols, leptosols, stagnosols and chernozems, with annual cropping systems with a SOC value of similar to 15 g.kg(-1) and a range of 30 g.kg(-1) in median. Most satellite-derived SOC spectral prediction models used limited preprocessing and were based on bare soil pixel retrieval after Normalized Difference Vegetation Index (NDVI) thresholding. About one third of these models used partial least squares regression (PLSR), while another third used random forest (RF), and the remaining included machine learning methods such as support vector machine (SVM). We did not find any studies either on deep learning methods or on all-performance evaluations and uncertainty analysis of spatial model predictions. Nevertheless, the literature examined here identifies satellite-based spectral information, especially derived under bare soil conditions, as an interesting approach that deserves further investigations. Future research includes considering the simultaneous analysis of imagery acquired at several dates i.e., temporal mosaicking, testing the influence of possible disturbing factors and mitigating their effects fusing mixed models incorporating non-spectral ancillary information