Non-linear spectral unmixing for monitoring rapidly salinizing coastal landscapes

Coastal farmlands in the eastern United States of America (USA) are increasingly suffering from rising soil salinity, rendering them unsuitable for economically productive agriculture. Saltwater intrusion (SWI) into the groundwater reservoir or soil salinization can result in land cover modification (e.g. reduced plant growth) or land cover conversion. Two primary examples of such land cover conversion are farmland to marsh or farmland to salt patches with no vegetation growth. However, due to varying spatial granularity of these conversions, it is challenging to quantify these land covers over a large geographic scale. To address this challenge, we evaluated a non-linear spectral unmixing approach with a Random Forest (RF) algorithm to quantify fractional abundance of salt patch and marshes. Using Sentinel-2 imagery from 2022, we generated gridded datasets for salt patches and marshes across the Delmarva Peninsula, and the associated uncertainty. Moreover, we developed two new spectral indices to enhance the spectral unmixing accuracy: the Normalized Difference Salt Patch Index (NDSPI) and the Modified Salt Patch Index (MSPI). We constructed two sets of ten RF models: one for salt patches and the other for marshes, achieving high (>99 %) training and testing accuracies for classification. The consistently high accuracy and low error values across different model runs demonstrate the method's reliability for classifying spectrally similar land cover classes in the mid-Atlantic region and beyond. Validation metrics for sub-pixel fractional abundances in the salt model revealed a moderate R-squared value of 0.50, and a high R-squared value of 0.90 for the marsh model. Our method complements labor-intensive field-based salinity measurements by offering a reproducible method that can be repeated annually and scaled up to cover large geographic regions.