Lithological mapping of ophiolite complex with emphasis on chromite and magnesite exploration using remote sensing techniques

This research employed the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) and Landsat Thematic Mapper (TM) data for lithological mapping and delineating of high potential chromite zone mineralization in ophiolite complexes. Abdasht, Soghan and Sikhoran chromite mining areas located in Sanandaj-Sirjan technically a part of the Esphandagheh ophiolite complex zone in Kerman province, southeastern of Iran have been selected for this research. In order to discriminate and to demarcate of the high potential chromite and magnesite rock zone, ASTER and Landsat TM bands properties have been utilized for running principal components analysis (PCA), band ratio (BR), minimum noise fraction (MNF), de-correlation stretch, log residual, spectral mapping methods and feature level fusion. A comparison between the image processing results with field investigation and primary geological map confirmed the concentration of chromite and magnesite mineralized zone associated with serpentinized dunite and hurzburgite. A new geological map showing high potential chromite zones and the boundary of lithological units was produced based on the interpretation of remote sensing data. The map can be used for geological exploration and mine engineering purposes. The data and methods used have emphasized high ability of the ASTER data to provide geological information for detecting chromite host rock such as serpentinized dunites and hurzburgite as well as lithological mapping at both district and regional scales. Additionally, Landsat TM data have also produced suitable results for lithological purposes on a regional scale. The approach used in this study is broadly applicable for exploring new chromite prospects and lithological mapping of the ophiolitic complexes especially in the arid and semi-arid regions of the earth

Satellite remote sensing for hydrothermal alteration minerals mapping of subtle geothermal system in unexplored aseismic environment

Mapping prospective geothermal (GT) resources and monitoring associated surface manifestations can be challenging and prohibitively expensive in subtle systems especially when using conventional survey methods. Remote sensing offers a synoptic and costeffective capability for identification of GT systems. The objective of this research is to refine and develop methods of identifying unconventional GT systems by evaluating the applicability of the ASTER, Landsat 8 and Hyperion satellite data for mapping hydrothermal alteration indicator minerals as proxy for detecting subtle GT targets in unexplored aseismic settings. The study area is Yankari Park in North Eastern Nigeria, characterized by the thermal springs; Wikki, Mawulgo, Gwana and Dimmil. Spectral Angle Mapper (SAM), Linear spectral Unmixing (LSU) and Mixture Tuned Matched Filtering (MTMF) were comparatively evaluated by using image derived spectra and corresponding library spectra for mapping pixel abundance of GT indicator minerals in a novel and efficient manner. The results indicated that employing image derived spectra from field validated and laboratory verified regions of interest as reference, gives more accurate results than using library spectra around known alteration zones remotely detectable on the imagery. The MTMF provided high performance subpixel target detection with an accuracy of 50-100% and 70-100% subpixel abundance for argillicphyllic- silicic and propylitic alteration mineral assemblages respectively, as compared to less than 10% for the same endmembers when using library spectra. The MTMF is thus best suited for mapping alterations associated with subtle GT systems than the less selective LSU. The per-pixel SAM was unsuitable for target detection of alteration indicators of interest with poor overall accuracy of 33.81% and 0.24 Kappa coefficient at 0.02 radian angle. Results of mapping thermally anomalous pixels do not conform to known locations of the thermal springs signifying the limitations of the current thermal sensors in mapping low temperature GT systems even at 60m spatial resolution. However, examining the spatial correlation of the anomaly areas with the major geologic structure systems from geological map of the study area indicates a close affinity between them and with previously reported thermal gradients within heat insulating sedimentary formations. This study establishes the integrative applicability of Multispectral and Hyperspectral data for mapping subtle GT targets in unexplored regions using in-situ validated alteration mineral mapping and thermal anomaly detection. This has significant implication for the GT green energy industry as the developed methods and GT prospect map could aid the prefeasibility stage narrowing of targets for in-depth geophysical, geochemical, geothermometric and related surveys

Classification of North Africa for Use as an Extended Pseudo Invariant Calibration Sites (Epics) for Radiometric Calibration and Stability Monitoring of Optical Satellite Sensors

An increasing number of Earth-observing satellite sensors are being launched to meet the insatiable demand for timely and accurate data to help the understanding of the Earth’s complex systems and to monitor significant changes to them. The quality of data recorded by these sensors is a primary concern, as it critically depends on accurate radiometric calibration for each sensor. Pseudo Invariant Calibration Sites (PICS) have been extensively used for radiometric calibration and temporal stability monitoring of optical satellite sensors. Due to limited knowledge about the radiometric stability of North Africa, only a limited number of sites in the region are used for this purpose. This work presents an automated approach to classify North Africa for its potential use as an extended PICS (EPICS) covering vast portions of the continent. An unsupervised classification algorithm identified 19 “clusters” representing distinct land surface types; three clusters were identified with spatial uncertainties within approximately 5% in the shorter wavelength bands and 3% in the longer wavelength bands. A key advantage of the cluster approach is that large numbers of pixels are aggregated into contiguous homogeneous regions sufficiently distributed across the continent to allow multiple imaging opportunities per day, as opposed to imaging a typical PICS once during the sensor’s revisit period. In addition, this work proposes a technique to generate a representative hyperspectral profile for these clusters, as the hyperspectral profile of these identified clusters are mandatory in order to utilize them for performing cross-calibration of optical satellite sensors. The technique was used to generate the profile for the cluster containing the largest number of aggregated pixels. The resulting profile was found to have temporal uncertainties within 5% across all the spectral regions. Overall, this technique shows great potential for generation of representative hyperspectral profiles for any North African cluster, which could allow the use of the entire North Africa Saharan region as an extended PICS (EPICS) dataset for sensor cross-calibration. Furthermore, this work investigates the performance of extended pseudo-invariant calibration sites (EPICS) in cross-calibration for one of Shrestha’s clusters, Cluster 13, by comparing its results to those obtained from a traditional PICS-based cross-calibration. The use of EPICS clusters can significantly increase the number of cross-calibration opportunities within a much shorter time period. The cross-calibration gain ratio estimated using a cluster-based approach had a similar accuracy to the cross-calibration gain derived from region of interest (ROI)-based approaches. The cluster-based cross-calibration gain ratio is consistent within approximately 2% of the ROI-based cross-calibration gain ratio for all bands except for the coastal and shortwave-infrared (SWIR) 2 bands. These results show that image data from any region within Cluster 13 can be used for sensor crosscalibration. Eventually, North Africa can be used a continental scale PICS

Matched filter stochastic background characterization for hyperspectral target detection

Algorithms exploiting hyperspectral imagery for target detection have continually evolved to provide improved detection results. Adaptive matched filters, which may be derived in many different scientific fields, can be used to locate spectral targets by modeling scene background as either structured geometric) with a set of endmembers (basis vectors) or as unstructured stochastic) with a covariance matrix. In unstructured background research, various methods of calculating the background covariance matrix have been developed, each involving either the removal of target signatures from the background model or the segmenting of image data into spatial or spectral subsets. The objective of these methods is to derive a background which matches the source of mixture interference for the detection of sub pixel targets, or matches the source of false alarms in the scene for the detection of fully resolved targets. In addition, these techniques increase the multivariate normality of the data from which the background is characterized, thus increasing adherence to the normality assumption inherent in the matched filter and ultimately improving target detection results. Such techniques for improved background characterization are widely practiced but not well documented or compared. This thesis will establish a strong theoretical foundation, describing the necessary preprocessing of hyperspectral imagery, deriving the spectral matched filter, and capturing current methods of unstructured background characterization. The extensive experimentation will allow for a comparative evaluation of several current unstructured background characterization methods as well as some new methods which improve stochastic modeling of the background. The results will show that consistent improvements over the scene-wide statistics can be achieved through spatial or spectral subsetting, and analysis of the results provides insight into the tradespaces of matching the interference, background multivariate normality and target exclusion for these techniques

Coupled canopy-atmosphere modelling for radiance-based estimation of vegetation properties

Vegetation is an important component of the Earth’s biosphere and therefore plays a crucial role in the carbon exchange of terrestrial ecosystems. Vegetation variables, such as leaf area index (LAI) and leaf chlorophyll content (Cab), can be monitored at global scale using remote sensing (RS). There are two main categories of approaches for estimating the vegetation variables from RS data: empirical and physically-based approaches. Physically-based approaches are more widely applicable because they rely on radiative transfer (RT) models, which can be adapted to the observation conditions and to the observed vegetation. For estimating the vegetation variables, however, the RT model has to be inverted, and this inversion is usually an ill-posed and under-determined problem. Several regularization methods have been proposed to allow finding stable and unique solutions: model coupling, using multi-angular data, using a priori information, as well as applying spatial or temporal constraints. Traditionally, radiance data measured at top-of the atmosphere (TOA) are pre-processed to top-of-canopy (TOC) reflectances. Corrections for atmospheric effects, and, if needed, for adjacency, directional, or topographic effects are usually applied sequentially and independently. Physically, however, these effects are inter-related, and each correction introduces errors. These errors propagate to the TOC reflectance data, which are used to invert the canopy RT model. The performance of the TOC approach is therefore limited by the errors introduced in the data during the pre-processing steps. This thesis proposes to minimize these errors by directly using measured TOA radiance data. In such a TOA approach, the atmospheric RT model, which is normally inverted to perform the atmospheric correction, is coupled to the canopy RT model. The coupled canopy-atmosphere model is inverted directly using the measured radiance data. Adjacency, directional and topographic effects can then be included in the coupled RT model. The same regularization methods as used for TOC approaches can be applied to obtain stable and unique estimates. The TOA approach was tested using four case studies based on mono-temporal data. A) The performance of the TOA approach was compared to a TOC approach for three Norway spruce stands in the Czech Republic, using near-nadir Compact High Resolution Imaging Spectrometer (CHRIS) data. The coupled model included canopy directional effects and simulated the CHRIS radiance data with similar accuracy as the canopy model simulated the atmospherically-corrected CHRIS data. Local sensitivity analyses showed that the atmospheric parameters had much less influence on the simulations than the vegetation parameters, and that the sensitivity profiles of the latter were very similar for both TOC and TOA approaches. The dimensionality of the estimation problem was evaluated to be 3 for both approaches. Canopy cover (Cv), fraction of bark material (fB), Cab, and leaf dry matter content (Cdm) were estimated using look-up tables (LUT) with similar accuracy with both approaches. B) Regularization using multi-angular data was tested for the TOA approach, using four angular CHRIS datasets, for the same three stands as used in A). The coupled model provided good simulations for all angles. The dimensionality increased from 3 to 6 when using all four angles. Two LUTs were built for each stand: a 4-variable LUT with fB, Cv, Cdm, and Cab, and a 7-variable LUT where leaf brown pigment concentration (Cs), dissociation factor (D), and tree shape factor (Zeta) were added. The results did not fully match the expectation that the more angles used, the more accurate the estimates become. Although their exploitation remains challenging, multi-angular data have higher potential than mono-angular data at TOA level. C) A Bayesian object-based approach was developed and tested on at-sensor Airborne Prism Experiment (APEX) radiance data for an agricultural area in Switzerland. This approach consists of two steps. First, up to six variables were estimated for each crop field object using a Bayesian optimization algorithm, using a priori information. Second, a LUT was built for each object with only LAI and Cab as free variables, thus spatially constraining the values of all other variables to the values obtained in the first step. The Bayesian object-based approach estimated LAI more accurately than a LUT with a Bayesian cost function approach. This case study relied on extensive field data allowing defining the objects and a priori data. D) The Bayesian object-based approach proposed in C) was applied to a simulated TOA Sentinel-2 scene, covering the area around Zurich, Switzerland. The simulated scene was mosaicked using seven APEX flight lines, which allowed including all spatial and spectral characteristics of Sentinel-2. Automatic multi-resolution segmentation and classification of the vegetated objects in four levels of brightness in the visible domain enabled defining the objects and a priori data without field data, allowing successful implementation of the Bayesian object-based approach. The research conducted in this thesis contributes to the improvement of the use of regularization methods in ill-posed RT model inversions. Three major areas were identified for further research: 1) inclusion of adjacency and topography effects in the coupled model, 2) addition of temporal constraints in the inversion, and 3) better inclusion of observation and model uncertainties in the cost function. The TOA approach proposed here will facilitate the exploitation of multi-angular, multi-temporal and multi-sensor data, leading to more accurate RS vegetation products. These higher quality products will support many vegetation-related applications.</p

An Integrated physics-based approach to demonstrate the potential of the Landsat Data Continuity Mission (LDCM) for monitoring coastal/inland waters

Monitoring coastal or inland waters, recognized as case II waters, using the existing Landsat technology is somewhat restricted because of its low Signal-to-Noise ratio (SNR) as well as its relatively poor radiometric resolution. As a primary task, we introduce a novel technique, which integrates the Landsat-7 data as a surrogate for LDCM with a 3D hydrodynamic model to monitor the dynamics of coastal waters near river discharges as well as in a small lake environment. The proposed approach leverages both the thermal and the reflective Landsat-7 imagery to calibrate the model and to retrieve the concentrations of optically active components of the water. To do so, the model is first calibrated by optimizing its thermal outputs with the surface temperature maps derived from the Landsat-7 data. The constituent retrieval is conducted in the second phase where multiple simulated concentration maps are provided to an in-water radiative transfer code (Hydrolight) to generate modeled surface reflectance maps. Prior to any remote sensing task, one has to ensure that a dataset comes from a well-calibrated imaging system. Although the calibration status of Landsat-7 has been regularly monitored over multiple desert sites, it was desired to evaluate its performance over dark waters relative to a well-calibrated instrument designed specifically for water studies. In the light of this, several Landsat- 7 images were cross-calibrated against the Terra-MODIS data over deep, dark waters whose optical properties remain relatively stable. This study is intended to lay the groundwork and provide a reference point for similar studies planned for the new Landsat. In an independent case study, the potential of the new Landsat sensor was examined using an EO-1 dataset and applying a spectral optimization approach over case II waters. The water constituent maps generated from the EO-1 imagery were compared against those derived from Landsat-7 to fully analyze the improvement levels pertaining to the new Landsat\u27s enhanced features in a water constituent retrieval framework

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

Advances in Data Mining Knowledge Discovery and Applications

Advances in Data Mining Knowledge Discovery and Applications aims to help data miners, researchers, scholars, and PhD students who wish to apply data mining techniques. The primary contribution of this book is highlighting frontier fields and implementations of the knowledge discovery and data mining. It seems to be same things are repeated again. But in general, same approach and techniques may help us in different fields and expertise areas. This book presents knowledge discovery and data mining applications in two different sections. As known that, data mining covers areas of statistics, machine learning, data management and databases, pattern recognition, artificial intelligence, and other areas. In this book, most of the areas are covered with different data mining applications. The eighteen chapters have been classified in two parts: Knowledge Discovery and Data Mining Applications