Veronika Kopačková

Mineral Classification of Land Surface Using Multispectral LWIR and Hyperspectral SWIR Remote-Sensing Data. A Case Study over the Sokolov Lignite Open-Pit Mines, the Czech Republic

Remote-sensing techniques offer an efficient alternative for mapping mining environments and assessing the impacts of mining activities. Airborne multispectral data in the thermal region and hyperspectral data in the optical region, acquired with the Airborne Hyperspectral Scanner (AHS) sensor over the Sokolov lignite open-pit mines in the Czech Republic, were analyzed. The emissivity spectrum was calculated for each vegetation-free land pixel in the longwave infrared (LWIR)-region image using the surface-emitted radiation, and the reflectance spectrum was derived from the visible, near-infrared and shortwave-infrared (VNIR–SWIR)-region image using the solar radiation reflected from the surface, after applying atmospheric correction. The combination of calculated emissivity, with the ability to detect quartz, and SWIR reflectance spectra, detecting phyllosilicates and kaolinite in particular, enabled estimating the content of the dominant minerals in the exposed surface. The difference between the emissivity values at λ = 9.68 µm and 8.77 µm was found to be a useful index for estimating the relative amount of quartz in each land pixel in the LWIR image. The absorption depth at around 2.2 µm in the reflectance spectra was used to estimate the relative amount of kaolinite in each land pixel in the SWIR image. The resulting maps of the spatial distribution of quartz and kaolinite were found to be in accordance with the geological nature and origin of the exposed surfaces and demonstrated the benefit of using data from both thermal and optical spectral regions to map the abundance of the major minerals around the mines.

Application of high altitude and ground-based spectroradiometry to mapping hazardous low-pH material derived from the Sokolov open-pit mine

Mineral spectroradiometry, both from airborne/spaceborne sensors and ground measurements, represents an alternative to conventional methods and an efficient way to characterize mines and assess the potential for acid mine drainage (AMD) discharge. High-altitude spectroradiometry (advanced spaceborne thermal emission and reflection radiometer [ASTER] satellite data) together with ground- and laboratory-based spectroradiometry (ASD Filedspec spectroradiometer) were employed in order to identify the locations of the most significant sources of AMD discharge at the Sokolov lignite open-pit mines, Czech Republic. As a result, a map with delineated low-pH zones was created and validated using the ground truth data.

Normalizing reflectance from different spectrometers and protocols with an internal soil standard

ABSTRACT The internal soil standard (ISS) concept in which a soil standard sample, exhibiting stable spectral performance, is used to normalize and align all other soil spectral measurements – was further examined herein. Different spectrometers (Spectral Evolution and ASD Spectral Pro) were used to measure a set of soil samples with the soil standards sample as a reference. Two sand dune samples served as the ISS to align measurements made under different conditions and the results were compared. It was shown that the ISS method was able to correct the spectral information from one spectrometer to another; however, the differences in the results obtained when using the two different soil standards are discussed. The main conclusion of this paper is that the soil spectral user community should adopt the ISS method for the benefit of all, and the sooner the better. This will allow much more effective exploitation of all data sets acquired on a daily basis by the growing soil spectral community that still lacks standardization procedures.

Utilization of hyperspectral image optical indices to assess the Norway spruce forest health status

The work is concerned with assessing the health status of trees of the Norway spruce species using airborne hyperspectral (HS) data (HyMap). The study was conducted in the Sokolov basin in the western part of the Czech Republic. First, statistics were employed to assess and validate diverse empirical models based on spectral information using the ground truth data (biochemically determined chlorophyll content). The model attaining the greatest accuracy ( D 718 / D 704 ∶ RMSE = 0.2055     mg / g , R 2 = 0.9370 ) was selected to produce a map of foliar chlorophyll concentrations ( C a b ). The C a b values retrieved from the HS data were tested together with other nonquantitative vegetation indicators derived from the HyMap image reflectance to create a statistical method allowing assessment of the condition of Norway spruce. As a result, we integrated the following HyMap derived parameters ( C a b , REP, and SIPI) to assess the subtle changes in physiological status of the macroscopically undamaged foliage of Norway spruce within the four studied test sites. Our classification results and the previously published studies dealing with assessing the condition of Norway spruce using chlorophyll contents are in a good agreement and indicate that this method is potentially useful for general applicability after further testing and validation.

Changes in Croplands as a Result of Large Scale Mining and the Associated Impact on Food Security Studied Using Time-Series Landsat Images

Geographic information systems and satellite remote sensing information are emerging technologies in land-cover change assessment. They now provide an opportunity to gain insights into land-cover change properties through the spatio-temporal data capture over several decades. The time series of Landsat images covering the 1985–2009 period is used here to explore the impacts of surface mining and reclamation, which constitute a dominant force in land-cover changes in the northwestern regions of the Czech Republic. Advanced quantification of the extent of mining activities is important for assessing how these land-cover changes affect ecosystem services such as croplands. The images employed from 1985, 1988, 1990, 2000, 2002, 2003, 2004, 2005, 2006, 2007, 2008, and 2009 assist in mapping the extent of surface mines and mine reclamation for large surface mines in a few selected areas of interest. The image processing techniques are based on pixel-by-pixel calculation of the vegetation index, such as NDVI. The NDVI values are classified into the defined classes based on CORINE Land Cover 2000 data in a 3280 km2 strip of Landsat images. This distribution of NDVI values is used to estimate the land-cover classes in the local areas of interest (184 km2, 368 km2, 737 km2, and 1,474 km2). Thus, the approximate land-cover stability of the 3,280 km2 strip during the whole 1985–2009 period is used to explore land-cover disturbances in the local areas of surface mines. In the case of NDVI, it also includes variations, presumably caused by seasonal vegetation effects, and local meteorological conditions. However, the main trends related to mining activities during the long-term period can be clearly understood. As a result, other objectives can be explored in the 1985–2009 period, such as cropland changes to other land use classes, changes of cropland patterns, and their impacts on food security. The presented spatio-temporal modeling based on long time series from 12 satellite images provides considerable experience for processing NDVI in the framework of identification of land-cover classes and also, to a certain degree, cropland variability with its impact on food security.

Modelling Diverse Soil Attributes with Visible to Longwave Infrared Spectroscopy Using PLSR Employed by an Automatic Modelling Engine

The study tested a data mining engine (PARACUDA®) to predict various soil attributes (BC, CEC, BS, pH, Corg, Pb, Hg, As, Zn and Cu) using reflectance data acquired for both optical and thermal infrared regions. The engine was designed to utilize large data in parallel and automatic processing to build and process hundreds of diverse models in a unified manner while avoiding bias and deviations caused by the operator(s). The system is able to systematically assess the effect of diverse preprocessing techniques; additionally, it analyses other parameters, such as different spectral resolutions and spectral coverages that affect soil properties. Accordingly, the system was used to extract models across both optical and thermal infrared spectral regions, which holds significant chromophores. In total, 2880 models were evaluated where each model was generated with a different preprocessing scheme of the input spectral data. The models were assessed using statistical parameters such as coefficient of determination (R2), square error of prediction (SEP), relative percentage difference (RPD) and by physical explanation (spectral assignments). It was found that the smoothing procedure is the most beneficial preprocessing stage, especially when combined with spectral derivation (1st or 2nd derivatives). Automatically and without the need of an operator, the data mining engine enabled the best prediction models to be found from all the combinations tested. Furthermore, the data mining approach used in this study and its processing scheme proved to be efficient tools for getting a better understanding of the geochemical properties of the samples studied (e.g., mineral associations).

ASSESSING FOREST HEALTH VIA LINKING THE GEOCHEMICAL PROPERTIES OF A SOIL PROFILE WITH THE BIOCHEMICAL PARAMETERS OF VEGETATION

The transfer of chemical elements/compounds within the soil–plant chain is a part of the biochemical cycling, and this system is controlled by biotic and abiotic factors which determine the final mobility and availability of chemical variables. Heavy metal contamination and low pH are stress factors that lead to changes in the contents of important foliage compounds, which can be used as non-specific indicators of plant stress. In this study, Norway spruce forests in the Sokolov region, being a part of the “Black Triangle,” were selected to assess geochemical and biochemical interactions in the natural soil/plant system. The authors studied the relationship between soil and spruce needle contents of macronutrients and potentially toxic elements and tested whether the soil parameters and their vertical distribution within a soil profile (two organic and two mineral horizons) affect foliage biochemical parameters (contents of photosynthetic pigments, phenolic compounds and lignin). Factor analysis was used to identify underlying variables that explained the pattern of correlations within and between the biochemical and geochemical datasets. Aluminum (Al) and arsenic (As) were identified as toxic elements with high bio-availability for spruce trees, and both were taken up by trees and translocated to the foliage. The correlations between two toxic element contents in needles (Al and As) and the contents of soluble phenolic compounds and total carotenoid to chlorophyll ratio suggest that these latter two biochemical parameters, which both proved to be sensitive to the soil geochemical conditions, can serve as suitable non-specific stress markers.

Testing a Modified PCA-Based Sharpening Approach for Image Fusion

Image data sharpening is a challenging field of remote sensing science, which has become more relevant as high spatial-resolution satellites and superspectral sensors have emerged. Although the spectral property is crucial for mineral mapping, spatial resolution is also important as it allows targeted minerals/rocks to be identified/interpreted in a spatial context. Therefore, improving the spatial context, while keeping the spectral property provided by the superspectral sensor, would bring great benefits for geological/mineralogical mapping especially in arid environments. In this paper, a new concept was tested using superspectral data (ASTER) and high spatial-resolution panchromatic data (WorldView-2) for image fusion. A modified Principal Component Analysis (PCA)-based sharpening method, which implements a histogram matching workflow that takes into account the real distribution of values, was employed to test whether the substitution of Principal Components (PC1–PC4) can bring a fused image which is spectrally more accurate. The new approach was compared to those most widely used—PCA sharpening and Gram–Schmidt sharpening (GS), both available in ENVI software (Version 5.2 and lower) as well as to the standard approach—sharpening Landsat 8 multispectral bands (MUL) using its own panchromatic (PAN) band. The visual assessment and the spectral quality indicators proved that the spectral performance of the proposed sharpening approach employing PC1 and PC2 improve the performance of the PCA algorithm, moreover, comparable or better results are achieved compared to the GS method. It was shown that, when using the PC1, the visible-near infrared (VNIR) part of the spectrum was preserved better, however, if the PC2 was used, the short-wave infrared (SWIR) part was preserved better. Furthermore, this approach improved the output spectral quality when fusing image data from different sensors (e.g., ASTER and WorldView-2) while keeping the proper albedo scaling when substituting the second PC.

Spectroscopy as a tool for geochemical modeling

This study focused on testing the feasibility of up-scaling ground-spectra-derived parameters to HyMap spectral and spatial resolution and whether they could be further used for a quantitative determination of the following geochemical parameters: As, pH and Clignite content. The study was carried on the Sokolov lignite mine as it represents a site with extreme material heterogeneity and high heavy-metal gradients. A new segmentation method based on the unique spectral properties of acid materials was developed and applied to the multi-line HyMap image data corrected for BRDF and atmospheric effects. The quantitative parameters were calculated for multiple absorption features identified within the VIS/VNIR/SWIR regions (simple band ratios, absorption band depth and quantitative spectral feature parameters calculated dynamically for each spectral measurement (centre of the absorption band (λ), depth of the absorption band (D), width of the absorption band (Width), and asymmetry of the absorption band (S)). The degree of spectral similarity between the ground and image spectra was assessed. The linear models for pH, As and the Clignite content of the whole and segmented images were cross-validated on the selected homogenous areas defined in the HS images using ground truth. For the segmented images, reliable results were achieved as follows: As: R2=0.84, Clignite: R2=0.88 and R2 pH: R2= 0.57.

Integration of Absorption Feature Information from Visible to Longwave Infrared Spectral Ranges for Mineral Mapping

Merging hyperspectral data from optical and thermal ranges allows a wider variety of minerals to be mapped and thus allows lithology to be mapped in a more complex way. In contrast, in most of the studies that have taken advantage of the data from the visible (VIS), near-infrared (NIR), shortwave infrared (SWIR) and longwave infrared (LWIR) spectral ranges, these different spectral ranges were analysed and interpreted separately. This limits the complexity of the final interpretation. In this study a presentation is made of how multiple absorption features, which are directly linked to the mineral composition and are present throughout the VIS, NIR, SWIR and LWIR ranges, can be automatically derived and, moreover, how these new datasets can be successfully used for mineral/lithology mapping. The biggest advantage of this approach is that it overcomes the issue of prior definition of endmembers, which is a requested routine employed in all widely used spectral mapping techniques. In this study, two different airborne image datasets were analysed, HyMap (VIS/NIR/SWIR image data) and Airborne Hyperspectral Scanner (AHS, LWIR image data). Both datasets were acquired over the Sokolov lignite open-cast mines in the Czech Republic. It is further demonstrated that even in this case, when the absorption feature information derived from multispectral LWIR data is integrated with the absorption feature information derived from hyperspectral VIS/NIR/SWIR data, an important improvement in terms of more complex mineral mapping is achieved.