Multiple Endmember Spectral Mixture Analysis Research Articles

Quantifying Qiyi Glacier Surface Dirtiness Using UAV and Sentinel-2 Imagery

The glacier surface is composed not only of ice or snow but also of a heterogeneous mixture of various materials. The presence of light-absorbing impurities darkens the glacier surface, reducing local reflectance and thereby accelerating the glacier melting process. However, our understanding of the spatial distribution of these impurities remains limited, and there is a lack of studies on quantifying the dirty degree of glacier surfaces. During the Sentinel satellite overpass on 21 August 2023, we used an ASD FieldSpec3 spectrometer to measure the reflectance spectra of glacier surfaces with varying degrees of dirtiness on the Qiyi glacier, Qinghai–Tibet Plateau. Using Multiple Endmember Spectral Mixture Analysis (MESMA), the Sentinel imagery was decomposed to generate fraction images of five primary ice surface materials as follows: coarse-grained snow, slightly dirty ice, moderately dirty ice, extremely dirty ice, and debris. Using unmanned aerial vehicle (UAV) imagery with a 0.05 m resolution, the primary ice surface was delineated and utilized as reference data to validate the fraction images. The findings revealed a strong correlation between the fraction images and the reference data (R2 ≥ 0.66, RMSE ≤ 0.21). Based on pixel-based classification from the UAV imagery, approximately 80% of the glacier surface is covered by slightly dirty ice (19.2%), moderately dirty ice (33.3%), extremely dirty ice (26.3%), and debris (1.2%), which significantly contributes to its darkening. Our study demonstrates the effectiveness of using Sentinel imagery in conjunction with MESMA to map the degree of glacier surface dirtiness accurately.

Inversion of Cotton Soil and Plant Analytical Development Based on Unmanned Aerial Vehicle Multispectral Imagery and Mixed Pixel Decomposition

In order to improve the accuracy of multispectral image inversion of soil and plant analytical development (SPAD) of the cotton canopy, image segmentation methods were utilized to remove the background interference, such as soil and shadow in UAV multispectral images. UAV multispectral images of cotton bud stage canopies at three different heights (30 m, 50 m, and 80 m) were acquired. Four methods, namely vegetation index thresholding (VIT), supervised classification by support vector machine (SVM), spectral mixture analysis (SMA), and multiple endmember spectral mixture analysis (MESMA), were used to segment cotton, soil, and shadows in the multispectral images of cotton. The segmented UAV multispectral images were used to extract the spectral information of the cotton canopy, and eight vegetation indices were calculated to construct the dataset. Partial least squares regression (PLSR), Random forest (FR), and support vector regression (SVR) algorithms were used to construct the inversion model of cotton SPAD. This study analyzed the effects of different image segmentation methods on the extraction accuracy of spectral information and the accuracy of SPAD modeling in the cotton canopy. The results showed that (1) The accuracy of spectral information extraction can be improved by removing background interference such as soil and shadows using four image segmentation methods. The correlation between the vegetation indices calculated from MESMA segmented images and the SPAD of the cotton canopy was improved the most; (2) At three different flight altitudes, the vegetation indices calculated by the MESMA segmentation method were used as the input variable, and the SVR model had the best accuracy in the inversion of cotton SPAD, with R2 of 0.810, 0.778, and 0.697, respectively; (3) At a flight altitude of 80 m, the R2 of the SVR models constructed using vegetation indices calculated from images segmented by VIT, SVM, SMA, and MESMA methods were improved by 2.2%, 5.8%, 13.7%, and 17.9%, respectively, compared to the original images. Therefore, the MESMA mixed pixel decomposition method can effectively remove soil and shadows in multispectral images, especially to provide a reference for improving the inversion accuracy of crop physiological parameters in low-resolution images with more mixed pixels.

MODIS daily cloud-gap-filled fractional snow cover dataset of the Asian Water Tower region (2000–2022)

Abstract. Accurate long-term daily cloud-gap-filled fractional snow cover products are essential for climate change and snow hydrological studies in the Asian Water Tower (AWT) region, but existing Moderate Resolution Imaging Spectroradiometer (MODIS) snow cover products are not sufficient. In this study, the multiple-endmember spectral mixture analysis algorithm based on automatic endmember extraction (MESMA-AGE) and the multistep spatiotemporal interpolation algorithm (MSTI) are used to produce the MODIS daily cloud-gap-filled fractional snow cover product over the AWT region (AWT MODIS FSC). The AWT MODIS FSC products have a spatial resolution of 0.005° and span from 2000 to 2022. The 2745 scenes of Landsat-8 images are used for the areal-scale accuracy assessment. The fractional snow cover accuracy metrics, including the coefficient of determination (R2), root mean squared error (RMSE) and mean absolute error (MAE), are 0.80, 0.16 and 0.10, respectively. The binarized identification accuracy metrics, including overall accuracy (OA), producer's accuracy (PA) and user's accuracy (UA), are 95.17 %, 97.34 % and 97.59 %, respectively. Snow depth data observed at 175 meteorological stations are used to evaluate accuracy at the point scale, yielding the following accuracy metrics: an OA of 93.26 %, a PA of 84.41 %, a UA of 82.14 % and a Cohen kappa (CK) value of 0.79. Snow depth observations from meteorological stations are also used to assess the fractional snow cover resulting from different weather conditions, with an OA of 95.36 % (88.96 %), a PA of 87.75 % (82.26 %), a UA of 86.86 % (78.86 %) and a CK of 0.84 (0.72) under the MODIS clear-sky observations (spatiotemporal reconstruction based on the MSTI algorithm). The AWT MODIS FSC product can provide quantitative spatial distribution information on snowpacks for mountain hydrological models, land surface models and numerical weather prediction in the Asian Water Tower region. This dataset is freely available from the National Tibetan Plateau Data Center at https://doi.org/10.11888/Cryos.tpdc.272503 (Jiang et al., 2022) or from the Zenodo platform at https://doi.org/10.5281/zenodo.10005826 (Jiang et al., 2023a).

Linking crown fire likelihood with post-fire spectral variability in Mediterranean fire-prone ecosystems

Background Fire behaviour assessments of past wildfire events have major implications for anticipating post-fire ecosystem responses and fuel treatments to mitigate extreme fire behaviour of subsequent wildfires. Aims This study evaluates for the first time the potential of remote sensing techniques to provide explicit estimates of fire type (surface fire, intermittent crown fire, and continuous crown fire) in Mediterranean ecosystems. Methods Random Forest classification was used to assess the capability of spectral indices and multiple endmember spectral mixture analysis (MESMA) image fractions (char, photosynthetic vegetation, non-photosynthetic vegetation) retrieved from Sentinel-2 data to predict fire type across four large wildfires Key results MESMA fraction images procured more accurate fire type estimates in broadleaf and conifer forests than spectral indices, without remarkable confusion among fire types. High crown fire likelihood in conifer and broadleaf forests was linked to a post-fire MESMA char fractional cover of about 0.8, providing a direct physical interpretation. Conclusions Intrinsic biophysical characteristics such as the fractional cover of char retrieved from sub-pixel techniques with physical basis are accurate to assess fire type given the direct physical interpretation. Implications MESMA may be leveraged by land managers to determine fire type across large areas, but further validation with field data is advised.

Automated Surface Runoff Estimation with the Spectral Unmixing of Remotely Sensed Multispectral Imagery

This work presents a methodology for the hydrological characterization of natural and urban landscapes, focusing on accurate estimations of infiltration capacity and runoff characteristics. By combining existing methods from the literature, we created a systemic process that integrates satellite-based vegetation maps, topography, and soil permeability data. This process generates a detailed vegetation classification and slope-corrected composite curve number (CNcα) map using information at the subpixel level, which is crucial for estimating excess runoff during intense precipitation events. The algorithm designed with this methodology is automated and utilizes freely accessible multispectral imagery. Leveraging the vegetation–impervious–soil (V-I-S) model, it is assumed that land cover comprises V-I-S components at each pixel. Automated Music and spectral Separability-based Endmember Selection is employed on a generic spectral library to obtain the most relevant V-I-S endmember spectra for a particular image, which is then employed in multiple endmember spectral mixture analysis to obtain V-I-S fraction maps. The derived fractions are utilized in combination with the Normalized Difference Vegetation Index and the Modified Normalized Difference Water Index to adapt the CNcα map to different seasons and climatic conditions. The methodology was applied to Esch-sur-Alzette, Luxembourg, over a four-year period to validate the methodology and quantify the increase in the impervious surface area in the commune and the relationship with the runoff dynamics. This approach provides valuable insights into infiltration and runoff dynamics across diverse temporal and geographic ranges.

Grain yield prediction using multi-temporal UAV-based multispectral vegetation indices and endmember abundance in rice

Timely and accurately predicting grain yield before harvest is of great importance for rice production management and grain trade. Numerous methods based on remote sensing (RS) technology have explored to estimate rice grain yield. However, in most cases, these methods are negatively affected by the mixed pixels of RS images due to the effect of background and panicles. To resolve such issues, the abundance information of leaves, soil and water background and rice panicles were extracted from the multispectral images of unmanned aerial vehicle (UAV) for the multiple endmember spectral mixture analysis (MESMA) method. Based on the analysis of contribution of vegetation index (VI) and abundance (ABD) to rice grain yield, a rice grain yield prediction model was developed by combining multi-stage time-series, ABD and VI. Results showed that MESMA can mitigate endmember variability in the estimation of rice yield, and achieve higher accuracy than conventional spectral mixture analysis (SMA). The NDRE and the sum of leaf and panicle abundance (ABDL+P) can be well correlated to grain yield before heading stage (R2 = 0.75) and at heading stage (R2 = 0.72), respectively. In comparison to them, the multi-stage time-series model consisted of VI and ABD [∑(VI&ABD)] produces the highest correlation (R2 = 0.80). It could also resolve the issues about the underestimation of yield using different datasets from various multi-spectral camera in comparison to the single-parameter models ∑VI and ∑ABD, and produce the good validation accuracy with R2 = 0.73, RRMSE= 0.22 and R2 = 0.75, RRMSE= 0.15, respectively. This study suggests that the combination of UAV multi-temporal VI and ABD data can achieve accurate prediction of rice grain yield. This method effectively utilizes the optical information of leaves and rice panicles and reduces background effects, which provides a new idea for accurate rice yield prediction.

First evaluation of fire severity retrieval from PRISMA hyperspectral data

The unprecedented availability of spaceborne hyperspectral data has great potential to provide fire severity estimates that align with post-fire management needs, overcoming complex logistics and data acquisition costs of airborne hyperspectral sensors, and the suboptimal sensitivity of broadband data to several post-fire ground components. We analyzed the feasibility of the PRISMA mission -one of the first spaceborne spectrometers operationally active- to assess fire severity by leveraging hyperspectral data dimensionality through the retrieval of sub-pixel components directly related to fire severity in the field. Multispectral data provided by Sentinel-2, commonly used in fire severity quantitative assessments, were used as a benchmark method. Multiple endmember spectral mixture analysis (MESMA) was used to retrieve fractional cover of char, photosynthetic vegetation (PV), as well as non-photosynthetic vegetation and soil (NPVS) from post-fire PRISMA Level 2D and Sentinel-2 Level 2A scenes in one of the largest wildfires ever recorded in the western Mediterranean Basin. Ground-truth data were obtained using the Composite Burn Index (CBI) to procure three field-measured severity metrics: vegetation, soil and site. The relationship between the CBI data on a continuum scale and the cover of char, PV and NPVS image fractions retrieved from PRISMA and Sentinel-2 was assessed through Random Forest regression (RFR). Likewise, Ordinal Forests (OF) algorithm was used to classify categorized CBI data (low, moderate and high fire severity). PRISMA-based RFR fire severity estimates at vegetation, soil and site levels (R2 = 0.64–0.79 and RMSE = 0.33–0.41) outperformed those of Sentinel-2 (R2 = 0.27–0.53 and RMSE = 0.54–0.60), and were in line with previous studies using airborne hyperspectral sensors at higher spatial resolution. Fire severity underestimation for high field CBI values was almost unnoticeable in the PRISMA estimates. Categorical fire severity, not currently estimated using hyperspectral data but with high interest in post-fire management, were accurately predicted by PRISMA-based OF classification, with consistent user's and producer's accuracy for each fire severity category. The high confusion between moderate and low/high fire severity categories, typical when unmixing broadband multispectral data, was overcome by the PRISMA-based classification scheme. Our results suggest that new spaceborne spectrometer missions can support reliable fire severity assessments equivalent to airborne spectrometers, but readily applicable to large-scale assessments of extreme wildfire events.

Regression-based surface water fraction mapping using a synthetic spectral library for monitoring small water bodies

ABSTRACT Small water bodies (SWBs), such as ponds and on-farm reservoirs, are a key part of the hydrological system and play important roles in diverse domains from agriculture to conservation. The monitoring of SWBs has been greatly facilitated by medium-spatial-resolution satellite images, but the monitoring accuracy is considerably affected by the mixed-pixel problem. Although various spectral unmixing methods have been applied to map sub-pixel surface water fractions for large water bodies, such as lakes and reservoirs, it is challenging to map SWBs that are small in size relative to the image pixel and have dissimilar spectral properties. In this study, a novel regression-based surface water fraction mapping method (RSWFM) using a random forest and a synthetic spectral library is proposed for mapping 10 m spatial resolution surface water fractions from Sentinel-2 imagery. The RSWFM inputs a few endmembers of water, vegetation, impervious surfaces, and soil to simulate a spectral library, and considers spectral variations in endmembers for different SWBs. Additionally, RSWFM applies noise-based data augmentation on pure endmembers to overcome the limitation often arising from the use of a small set of pure spectra in training the regression model. RSWFM was assessed in ten study sites and compared with the fully constrained least squares (FCLS) linear spectral mixture analysis, multiple endmember spectral mixture analysis (MESMA), and the nonlinear random forest (RF) regression without data-augmentation. The results showed that RSWFM decreases the water fraction mapping errors by ~ 30%, ~15%, and ~ 11% in root mean square error compared with the linear FCLS, MESMA unmixings, and the nonlinear RF regression without data-augmentation respectively. RSWFM has an accuracy of approximately 0.85 in R2 in estimating the area of SWBs smaller than 1 ha.

Estimation and Validation of Sub-Pixel Needleleaf Cover Fraction in the Boreal Forest of Alaska to Aid Fire Management

Wildfires, which are a natural part of the boreal ecosystem in Alaska, have recently increased in frequency and size. Environmental conditions (high temperature, low precipitation, and frequent lightning events) are becoming favorable for severe fire events. Fire releases greenhouse gasses such as carbon dioxide into the environment, creating a positive feedback loop for warming. Needleleaf species are the dominant vegetation in boreal Alaska and are highly flammable. They burn much faster due to the presence of resin, and their low-lying canopy structure facilitates the spread of fire from the ground to the canopy. Knowing the needleleaf vegetation distribution is crucial for better forest and wildfire management practices. Our study focuses on needleleaf fraction mapping using a well-documented spectral unmixing approach: multiple endmember spectral mixture analysis (MESMA). We used an AVIRIS-NG image (5 m), upscaled it to 10 m and 30 m spatial resolutions, and applied MESMA to all three images to assess the impact of spatial resolution on sub-pixel needleleaf fraction estimates. We tested a novel method to validate the fraction maps using field data and a high-resolution classified hyperspectral image. Our validation method produced needleleaf cover fraction estimates with accuracies of 73%, 79%, and 78% for 5 m, 10 m, and 30 m image data, respectively. To determine whether these accuracies varied significantly across different spatial scales, we used the McNemar statistical test and found no significant differences between the accuracies. The findings of this study enhance the toolset available to fire managers to manage wildfire and for understanding changes in forest demography in the boreal region of Alaska across the high-to-moderate resolution scale.

Mapping Alpine Grassland Fraction Coverage Using Zhuhai-1 OHS Imagery in the Three River Headwaters Region, China

The widely spread alpine grassland ecosystem in the Three River Headwaters Region (TRHR) plays an essential ecological role in carbon sequestration and soil and water conservation. In this study, we test the latest high spatial resolution hyperspectral (Zhuhai-1 OHS) remote sensing imagery to examine different alpine grassland coverage levels using Multiple Endmember Spectral Mixture Analysis (MESMA). Our results suggest that the 3-endmember (3-EM) MESMA model can provide the highest image pixel unmixing percentage, with a percentage exceeding 97% and 96% for pixel scale and landscape scale, respectively. The overall accuracy shows that Zhuhai-1 OHS imagery obtained the highest overall accuracy (83.7%, k = 0.77) in the landscape scale, but in the pixel scale, it is not as good as Landsat 8 OLI imagery. Overall, we can conclude that the hyperspectral imagery combined 3-EM MESMA model performs better in both pixel scale and landscape scale alpine grassland coverage mapping, while the multispectral imagery with the 3-EM MESMA model can satisfy requirements of alpine grassland coverage mapping at the pixel scale. The approaches and workflow to mapping alpine grassland in this study can help monitor alpine grassland degradation; not only in the Qinghai–Tibetan Plateau (QTP), but also in other grassland ecosystems.

Effects of spectral variability due to sediment and bottom characteristics on remote sensing for suspended sediment in shallow rivers

Advances in remote sensing techniques for water environments have led to acquisition of abundant data on suspended sediment concentration (SSC). However, confounding factors, such as particle sizes, mineral properties, and bottom materials, have not been fully studied, despite their substantial interference with the detection of intrinsic signals of suspended sediments. Therefore, we investigated the spectral variability arising from the sediment and bottom using laboratory and field-scale experiments. In the laboratory experiment, we focused on measuring spectral characteristics of suspended sediment according to particle size and sediment type. The laboratory experiment was conducted under conditions of completely mixed sediment and non-bottom reflectance using a specially designed rotating horizontal cylinder. To investigate the effects of different channel bottoms under sediment-laden flow conditions, we performed sediment tracer tests in field-scale channels comprising sand and vegetated bottoms. Based on experimental datasets, we performed spectral analysis and multiple endmember spectral mixture analysis (MESMA) to quantify the effect of spectral variability of sediment and bottom on the relationship between hyperspectral data and SSC. The results showed that optimal spectral bands were precisely estimated under non-bottom reflectance conditions, and the effective wavelengths depended on the sediment type. The fine sediments had a higher backscattering intensity compared to the coarse sediments, and the reflectance difference according to the particle size difference increased as the SSC increased. However, in the field-scale experiment, the bottom reflectance substantially decreased the R2 in the relationship between the hyperspectral data and SSC. Nevertheless, MESMA can quantify the contribution of suspended sediment and bottom signals as fractional images. Moreover, the suspended sediment fraction had a clear exponential relationship with SSC in all cases. We conclude that MESMA-driven sediment fractions could be an important alternative for estimating SSC in shallow rivers, as it quantifies the contributions of each factor and then minimizes the bottom effect.

Mapping Fractional Vegetation Coverage across Wetland Classes of Sub-Arctic Peatlands Using Combined Partial Least Squares Regression and Multiple Endmember Spectral Unmixing

Vegetation communities play a key role in governing the atmospheric-terrestrial fluxes of water, carbon, nutrients, and energy. The expanse and heterogeneity of vegetation in sub-arctic peatland systems makes monitoring change at meaningful spatial resolutions and extents challenging. We use a field-collected spectral endmember reference library to unmix hyperspectral imagery and map vegetation coverage at the level of plant functional type (PFT), across three wetland sites in sub-arctic Alaska. This study explores the optimization and parametrization of multiple endmember spectral mixture analysis (MESMA) models to estimate coverage of PFTs across wetland classes. We use partial least squares regression (PLSR) to identify a parsimonious set of critical bands for unmixing and compare the reference and modeled coverage. Unmixing, using a full set of 110-bands and a smaller set of 4-bands, results in maps that effectively discriminate between PFTs, indicating a small investment in fieldwork results in maps mirroring the true ground cover. Both sets of spectral bands differentiate between PFTs, but the 4-band unmixing library results in more accurate predictive mapping with lower computational cost. Reducing the unmixing reference dataset by constraining the PFT endmembers to those identified in the field-site produces only a small advantage for mapping, suggesting extensive fieldwork may not be necessary for MESMA to have a high explanatory value in these remote environments.

Spectral Mixture Modeling of an ASTER Bare Soil Synthetic Image Using a Representative Spectral Library to Map Soils in Central-Brazil

Pedological maps in suitable scales are scarce in most countries due to the high costs involved in soil surveying. Therefore, methods for surveying and mapping must be developed to overpass the cartographic material obtention. In this sense, this work aims at assessing a digital soil map (DSM) built by multispectral data extrapolation from a source area to a target area using the ASTER time series modeling technique. For that process, eight representative toposequences were established in two contiguous micro-watersheds, with a total of 42 soil profiles for analyses and classification. We found Ferralsols, Plinthosols, Regosols, and a few Cambisols, Arenosols, Gleisols, and Histosols, typical of tropical regions. In the laboratory, surface soil samples were submitted to spectral readings from 0.40 µm to 2.50 µm. The soil spectra were morphologically interpreted, identifying shapes and main features typical of tropical soils. Soil texture grouped the curves by cluster analysis, forming a spectral library (SL). In parallel, an ASTER time series (2001, 2004, and 2006) was processed, generating a bare soil synthetic soil image (SySI) covering 39.7% of the target area. Multiple Endmember Spectral Mixture Analysis modeled the SL on the SySI generating DSM with 73% of Kappa index, in which identified about 77% is covered by rhodic Ferralsols. Besides the overestimation, the DSM represented the study area’s pedodiversity. Given the discussion raised, we consider including subsoil data and other features using other sensors in operations modeled by machine learning algorithms to improve results.

Characterizing Tree Species in Northern Boreal Forests Using Multiple-Endmember Spectral Mixture Analysis and Multi-Temporal Satellite Imagery

Northern boreal forests are characterized by open stands whereby trees, understory background, and shadow are all significant components of the spectral response within a pixels’ spatial footprint. To overcome this mixed pixel problem, accurate spectral characterization of these (endmember) components is necessary for spectral mixture analysis (SMA) to generate forest classifications at the species level. Obtaining these endmember spectra in the field, however, can be difficult or impossible. This study examined whether image endmember spectra can be identified using forest inventory information to derive dominant tree species classifications. This was tested using multiple-endmember SMA (MESMA) and single- and multi-date Landsat imagery of a forested area in the Northwest Territories, Canada. Image classifications (n = 80) were generated based on 20 image-date combinations and four unmixing models. Accuracies of 80% and 82% were achieved for open and medium dense forest stands, respectively using multi-date imagery, which outperformed single-date imagery acquired at peak phenology. The overall accuracy is 72%; lower due to challenges in very open stands. The multi-date MESMA approach was robust for both compositionally pure and mixed stands. The approach merits further investigation, particularly within the context of the increasing availability of regional-scale satellite imagery enabling composite time-series and spectral-temporal image features.

Estimation of boreal forest floor lichen cover using hyperspectral airborne and field data

Lichens are sensitive to competition from vascular plants, intensive silviculture, pollution and reindeer and caribou grazing, and can therefore serve as indicators of environmental changes. Hyperspectral remote sensing data has been proved promising for estimation of plant diversity, but its potential for forest floor lichen cover estimation has not yet been studied. In this study, we investigated the use of hyperspectral data in estimating ground lichen cover in boreal forest stands in Finland. We acquired airborne and in situ hyperspectral data of lichen-covered forest plots, and applied multiple endmember spectral mixture analysis to estimate the fractional cover of ground lichens in these plots. Estimation of lichen cover based on in situ spectral data was very accurate (coefficient of determination (r2) 0.95, root mean square error (RMSE) 6.2). Estimation of lichen cover based on airborne data, on the other hand, was fairly good (r2 0.77, RMSE 11.7), but depended on the choice of spectral bands. When the hyperspectral data were resampled to the spectral resolution of Sentinel-2, slightly weaker results were obtained. Tree canopy cover near the flight plots was weakly related to the difference between estimated and measured lichen cover. The results also implied that the presence of dwarf shrubs could influence the lichen cover estimates.

Spectral mixture analysis of AVIRIS-NG hyperspectral data for material identification and classification for the part of Kolkata city

Identification and extraction of material in the urban scenario are always challenging due to the presence of different heterogeneous materials. Airborne hyperspectral sensors with the advantage of both high spectral and spatial resolution can be used in such areas for the identification and mapping of urban materials. Manual extraction of the endmembers in highly dense heterogeneous urban regions is a very difficult task to achieve, but with automated extraction techniques, it yields promising results. To counter the mixed pixels problem in heterogeneous urban areas automated endmember extraction techniques like Endmember Average RMS (EAR), Minimum Average Spectral Angle (MASA), and Count Based Endmember Selection (CoB) was used for pure endmember identification in the densely populated urban area of Kolkata region using AVIRIS NG imagery. Spectral Unmixing Multiple Endmember Spectral Mixture Analysis (MESMA), Linear Spectral Unmixing (LSU), and Mixture Tuned Matched Filtering (MTMF) was performed on AVIRIS NG data using the extracted endmembers of different urban materials. The rule-based classification was performed on abundance images obtained from MESMA, MTMF, and LSU techniques for comparative analysis. Classification using MESMA abundance images outperformed in mapping mixed urban environments with an overall accuracy (OA) of 88.6 % followed by the rule-based classified output of MTMF and LSU with OA of 75.27% and 72.45% respectively.

Mapping of understorey invasive plant species clusters of Lantana camara and Chromolaena odorata using airborne hyperspectral remote sensing

The present study highlights delineation of the distribution of understorey invasive plant species (IPS) in Mudumalai Tiger Reserve (MTR), Western Ghats, India using very high resolution airborne imaging spectroscopy images by evaluating the performance of a Multiple Endmember Spectral Mixture Analysis (MESMA). Lantana camara L. and Chromolaena odorata (L.) R.M.King & H.Rob. have imposed a serious threat to natural ecosystems in MTR by altering the ecosystem processes. The Airborne Visible-Infrared Imaging Spectrometer-Next Generation (AVIRIS-NG) based analysis delineated the distribution of understorey IPS, viz., L. camara and C. odorata in MTR, with an overall accuracy of 87 % (k = 0.84) and 84 % (k = 0.8), respectively, with minimal root-mean-square error (RMSE) of MESMA classification (0.001 to 0.0025). In MTR, both the IPS were distributed in large parts and in association with all dominant canopy tree species barring Eucalyptus species. L. camara was majorly found in the southern parts of central MTR with ‘high’ to ‘very high’ presence covering 23.09 % and 11.51 % of the total area of central MTR, respectively, in association with Terminalia anogeissiana Gere & Boatwr. and Tectona grandis L.f. Whereas, C. odorata was dominantly found with high (18.11 % area of central MTR) to very high (11.85 % area of central MTR) intensity in the south and southeastern parts within the vicinity of T. anogeissiana and Terminalia elliptica Willd. The study highlights the utilization of high spectral agility in mapping of understorey IPS through spectral differentiation of the patchy distribution of IPS and its intermixing with other local species may be helpful in ecological restoration, conservation of biodiversity, and provisioning long-term management of ecosystem services in highly invaded tropical forests.

Drought impact on cropland use monitored with AVIRIS imagery in Central Valley, California

During 2012–2016 California experienced the longest and most severe drought in the last centuries. This water scarcity led to an increase in non-cultivated croplands during this period. The objective of this study was to quantify agricultural trends in the Central Valley (California) at peak growth from 2013 to 2016 during the drought and in 2017–2018 post-drought. For this purpose, we analysed yearly official harvested area reported at county level for the main crops and compared them to visible-shortwave infrared (VSWIR) spectra acquired by the Airborne Visible/Infrared Imaging-Spectrometer (AVIRIS-classic) over 2334 km2 of the Central Valley each year. Multiple Endmember Spectral Mixture Analysis (MESMA) was applied to AVIRIS data to estimate green vegetation (GV), non-photosynthetic vegetation (NPV) and soil fractions in crop fields each year. MESMA and crop reports (R2 = 0.9) showed that soil (i.e.; non-cultivated areas) increased during the summers of the drought; with the smallest GV area in 2015, the second year classified with exceptional drought in this period. According to MESMA, 34 % of the cropland was covered by GV in 2015, and 69.5·104 ha according to the crop reports. MESMA also registered the highest value of soil area in 2015 (48 %). The year with most cultivated area was 2017, the wettest year in the studied period, with 54 % of the croplands covered by GV and 75.2·104 ha. This study verified that the non-cultivated areas increased in the Central Valley during the exceptional drought period and validated the use of AVIRIS imagery to monitor broad-scale cropland use changes in future climatic extreme events.

Extraction of rocky desertification information in karst area by using different multispectral sensor data and multiple endmember spectral mixture analysis method

As a typical form of land degradation, karst rocky desertification seriously restricts the development of the regional social economy and seriously threatens the living environment of residents. Fractional vegetation cover (fVC) and bare rock fraction (fBR) are important indicators to identify and evaluate rocky desertification. However, it is a great challenge to obtain fVC and fBR due to the complex terrain and fragmentation of karst rocky desertification areas. In this study, comparisons between Sentinel-2A Multispectral Instrument (Sentinel-2), Landsat-8 Operational Land Imager (Landsat-8), and GF-6 Wide Field View (GF-6) sensors for retrieving fVC and fBR are presented. The multiple endmember spectral mixture analysis (MESMA) and measured spectral dates were used to overcome the limitations of Spectral mixture analysis (SMA). Subsequently, fVC and fBR were validated using root mean square error (RMSE), mean absolute error (MAE), and coefficient of determination (R2). The results showed that: 1) Sentinel-2 performed best in estimating fVC and obtained the highest accuracy (R2 = 0.6259, root mean square error = 0.1568, mean absolute error = 0.1215), followed by GF-6 and Landsat 8; in the extraction of the fBR, the performance was relatively the same, and the implementation of Sentinel-2 was also the best (R2 = 0.4911, root mean square error = 0.0714, mean absolute error = 0.0539), followed by GF-6 and Landsat 8. 2) Sentinel-2 images have higher resolution, the narrowest band range, and the most significant number of bands, which can better extract information about fVC and fBR in rocky desertification areas. 3) For the three optical sensors, the spatial resolution of the images is more important to extract the information of fVC and fBR in the rocky desertification areas. 4) In general, the extraction accuracy of fBR is not as good as that of fVC. The complicated ecological and geological environment of decertified areas has more influence on the effect of extraction of the fBR. 5) The Sentinel-2 achieves high accuracy for both fVC and fBR under different-level application scenarios. It thus has great potential for application in rocky desertification information extraction.

Mapping the invasive palm species Arenga obtusifolia using multiple endmember spectral mixture analysis (MESMA) and PRISMA hyperspectral data in Ujung Kulon National Park, Indonesia

Arenga obtusifolia or “langkap” is an invasive palm that interferes with the growth of native vegetation in Ujung Kulon National Park. This study aims to map A. obtusifolia distribution in Ujung Kulon using PRISMA hyperspectral imagery and MESMA (Multiple Endmember Spectral Mixture Analysis). 38 spectra samples for the A. obtusifolia, sand, wetland vegetation and other vegetation were collected, and optimized using the IES (Iterative Endmember Selection), Multiple Signal Classification (MUSIC) and Automated MUsic and Spectral Separability based Endmember Selection (AMUSES) algorithms. Here, different schemes of MESMA were tested by using spectral weighting (SPEC) and stable zone unmixing (STABLE). Our results demonstrated that MESMA and PRISMA data was able to identify A. obtusifolia with the accuracy ranging from 55.36 (AMUSES-STABLE) to 83.93% (MUSIC-STABLE). Our results indicated the potential of MESMA and PRISMA hyperspectral data using optimized spectra from MUSIC and stable zone unmixing for mapping invasive A. obtusifolia.