Increasing human populations around the global coastline have caused extensive loss, degradation and fragmentation of coastal ecosystems, threatening the delivery of important ecosystem services
. As ...a result, alarming losses of mangrove, coral reef, seagrass, kelp forest and coastal marsh ecosystems have occurred
. However, owing to the difficulty of mapping intertidal areas globally, the distribution and status of tidal flats-one of the most extensive coastal ecosystems-remain unknown
. Here we present an analysis of over 700,000 satellite images that maps the global extent of and change in tidal flats over the course of 33 years (1984-2016). We find that tidal flats, defined as sand, rock or mud flats that undergo regular tidal inundation
, occupy at least 127,921 km
(124,286-131,821 km
, 95% confidence interval). About 70% of the global extent of tidal flats is found in three continents (Asia (44% of total), North America (15.5% of total) and South America (11% of total)), with 49.2% being concentrated in just eight countries (Indonesia, China, Australia, the United States, Canada, India, Brazil and Myanmar). For regions with sufficient data to develop a consistent multi-decadal time series-which included East Asia, the Middle East and North America-we estimate that 16.02% (15.62-16.47%, 95% confidence interval) of tidal flats were lost between 1984 and 2016. Extensive degradation from coastal development
, reduced sediment delivery from major rivers
, sinking of riverine deltas
, increased coastal erosion and sea-level rise
signal a continuing negative trajectory for tidal flat ecosystems around the world. Our high-spatial-resolution dataset delivers global maps of tidal flats, which substantially advances our understanding of the distribution, trajectory and status of these poorly known coastal ecosystems.
Maps that categorise the landscape into discrete units are a cornerstone of many scientific, management and conservation activities. The accuracy of these maps is often the primary piece of ...information used to make decisions about the mapping process or judge the quality of the final map. Variance is critical information when considering map accuracy, yet commonly reported accuracy metrics often do not provide that information. Various resampling frameworks have been proposed and shown to reconcile this issue, but have had limited uptake. In this paper, we compare the traditional approach of a single split of data into a training set (for classification) and test set (for accuracy assessment), to a resampling framework where the classification and accuracy assessment are repeated many times. Using a relatively simple vegetation mapping example and two common classifiers (maximum likelihood and random forest), we compare variance in mapped area estimates and accuracy assessment metrics (overall accuracy, kappa, user, producer, entropy, purity, quantity/allocation disagreement). Input field data points were repeatedly split into training and test sets via bootstrapping, Monte Carlo cross-validation (67:33 and 80:20 split ratios) and k-fold (5-fold) cross-validation. Additionally, within the cross-validation, we tested four designs: simple random, block hold-out, stratification by class, and stratification by both class and space. A classification was performed for every split of every methodological combination (100’s iterations each), creating sampling distributions for the mapped area of each class and the accuracy metrics. We found that regardless of resampling design, a single split of data into training and test sets results in a large variance in estimates of accuracy and mapped area. In the worst case, overall accuracy varied between ~40–80% in one resampling design, due only to random variation in partitioning into training and test sets. On the other hand, we found that all resampling procedures provided accurate estimates of error, and that they can also provide confidence intervals that are informative about the performance and uncertainty of the classifier. Importantly, we show that these confidence intervals commonly encompassed the magnitudes of increase or decrease in accuracy that are often cited in literature as justification for methodological or sampling design choices. We also show how a resampling approach enables generation of spatially continuous maps of classification uncertainty. Based on our results, we make recommendations about which resampling design to use and how it could be implemented. We also provide a fully worked mapping example, which includes traditional inference of uncertainty from the error matrix and provides examples for presenting the final map and its accuracy.
•Resampling designs are compared for image classification and accuracy assessment.•Resampling provides robust accuracy and area estimates with confidence intervals.•A single split into training/test data often gives inaccurate or misleading results.•Recommendations, examples and code are given for implementing resampling.
The utility of land cover maps for natural resources management relies on knowing the uncertainty associated with each map. The continuous advances typical of remote sensing, including the increasing ...availability of higher spatial and temporal resolution satellite data and data analysis capabilities, have created both opportunities and challenges for improving the application of accuracy assessment. There are well established accuracy assessment methods, but their underlying assumptions have not changed much in the last couple decades. Consequently, revisiting how map error and accuracy have been performed and reported over the last two decades is timely, to highlight areas where there is scope for better utilization of emerging opportunities. We conducted a quantitative literature review on accuracy assessment practices for mapping via remote sensing classification methods, in both terrestrial and marine environments. We performed a structured search for land and benthic cover mapping, limiting our search to journals within the remote sensing field, and papers published between 1998–2017. After an initial screening process, we assembled a database of 282 papers, and extracted and standardized information on various components of their reported accuracy assessments. We discovered that only 56% of the papers explicitly included an error matrix, and a very limited number (14%) reported overall accuracy with confidence intervals. The use of kappa continues to be standard practice, being reported in 50.4% of the literature published on or after 2012. Reference datasets used for validation were collected using a probability sampling design in 54% of the papers. For approximately 11% of the studies, the sampling design used could not be determined. No association was found between classification complexity (i.e. number of classes) and measured accuracy, independent from the size of the study area. Overall, only 32% of papers included an accuracy assessment that could be considered reproducible; that is, they included a probability-based sampling scheme to collect the reference dataset, a complete error matrix, and provided sufficient characterization of the reference datasets and sampling unit. Our findings indicate that considerable work remains to identify and adopt more statistically rigorous accuracy assessment practices to achieve transparent and comparable land and benthic cover maps.
In the Yellow Sea region of East Asia, tidal wetlands are the frontline ecosystem protecting a coastal population of more than 60 million people from storms and sea-level rise. However, unprecedented ...coastal development has led to growing concern about the status of these ecosystems. We developed a remote-sensing method to assess change over ~4000 km of the Yellow Sea coastline and discovered extensive losses of the region's principal coastal ecosystem - tidal flats - associated with urban, industrial, and agricultural land reclamations. Our analysis revealed that 28% of tidal flats existing in the 1980s had disappeared by the late 2000s (1.2% annually). Moreover, reference to historical maps suggests that up to 65% of tidal flats were lost over the past five decades. With the region forecast to be a global hotspot of urban expansion, development of the Yellow Sea coastline should pursue a course that minimizes the loss of remaining coastal ecosystems.
Our understanding of Earth surface processes is rapidly advancing as new remote sensing technologies such as LiDAR and close-range digital photogrammetry become more accessible and affordable. A ...very-high spatial resolution digital terrain model (DTM) and orthophoto mosaic (mm scale) were produced using close-range digital photogrammetry based on ‘Structure-from-Motion’ (SfM) algorithms for a 250m transect along a shallow coral reef flat on Heron Reef, Great Barrier Reef. The precise terrain data were used to characterise surface roughness, a critical factor affecting ecological and physical processes on the reef. Three roughness parameters, namely the root mean square height, tortuosity (or rugosity) and fractal dimension, were derived and compared in order to asses which one better characterises reef flat roughness. The typical relief across the shallow reef flat was 0.1m with a maximum value of 0.42m. Coral reef terrain roughness, as characterised by the three chosen parameters, generally increased towards the middle of the transect where live coral covers most of the reef flat and decreases towards the edges of the transect. The fractal dimension (values ranging from 2.2 to 2.59) best characterised reef roughness, as evidenced by a closer agreement with the distribution of known coral benthic substrates. This is the first study quantifying scale-independent roughness of a coral reef at benthic and biotope/patch levels (cm-m). The readily available and cost-effective methods presented are highly appropriate for data collection, processing and analysis to generate very-high spatial resolution DTMs and orthophoto mosaics of shallow and energetic coral reefs.
•A very-high spatial resolution DTM was produced for an intertidal coral reef flat.•SfM is a cost-effective approach for mapping the terrain of shallow reef flats.•Methods presented are readily available and require minimum training.•Coral reef roughness was better characterised using the fractal dimension parameter.•This is the first study quantifying scale-independent roughness at cm resolution.
Big data reveals new, stark pictures of the state of our environments. It also reveals ‘bright spots’ amongst the broad pattern of decline and—crucially—the key conditions for these cases. Big data ...analyses could benefit the planet if tightly coupled with ongoing sustainability efforts.
In the Australian summer season of 2022, exceptional rainfall events occurred in Southeast Queensland and parts of New South Wales, leading to extensive flooding of rural and urban areas. Here, we ...map the extent of flooding in the city of Brisbane and evaluate the change in electricity usage as a proxy for flood impact using VIIRS nighttime brightness imagery. Scanning a wide range of possible sensors, we used pre-flood and peak-flood PlanetScope imagery to map the inundated areas, using a new spectral index we developed, the Normalized Difference Inundation Index (NDII), which is based on changes in the NIR reflectance due to sediment-laden flood waters. We compared the Capella-Space X-band/HH imaging radar data captured at peak-flood date to the PlanetScope-derived mapping of the inundated areas. We found that in the Capella-Space image, significant flooded areas identified in PlanetScope imagery were omitted. These omission errors may be partly explained by the use of a single-date radar image, by the X-band, which is partly scattered by tree canopy, and by the SAR look angle under which flooded streets may be blocked from the view of the satellite. Using VIIRS nightly imagery, we were able to identify grid cells where electricity usage was impacted due to the floods. These changes in nighttime brightness matched both the inundated areas mapped via PlanetScope data as well as areas corresponding with decreased electricity loads reported by the regional electricity supplier. Altogether we demonstrate that using a variety of optical and radar sensors, as well as nighttime and daytime sensors, enable us to overcome data gaps and better understand the impact of flood events. We also emphasize the importance of high temporal revisit times (at least twice daily) to more accurately monitor flood events.
Providing accurate maps of coral reefs where the spatial scale and labels of the mapped features correspond to map units appropriate for examining biological and geomorphic structures and processes ...is a major challenge for remote sensing. The objective of this work is to assess the accuracy and relevance of the process used to derive geomorphic zone and benthic community zone maps for three western Pacific coral reefs produced from multi-scale, object-based image analysis (OBIA) of high-spatial-resolution multi-spectral images, guided by field survey data. Three Quickbird-2 multi-spectral data sets from reefs in Australia, Palau and Fiji and georeferenced field photographs were used in a multi-scale segmentation and object-based image classification to map geomorphic zones and benthic community zones. A per-pixel approach was also tested for mapping benthic community zones. Validation of the maps and comparison to past approaches indicated the multi-scale OBIA process enabled field data, operator field experience and a conceptual hierarchical model of the coral reef environment to be linked to provide output maps at geomorphic zone and benthic community scales on coral reefs. The OBIA mapping accuracies were comparable with previously published work using other methods; however, the classes mapped were matched to a predetermined set of features on the reef.
Providing accurate maps of mangroves, where the spatial scales of the mapped features correspond to the ecological structures and processes, as opposed to pixel sizes and mapping approaches, is a ...major challenge for remote sensing. This study developed and evaluated an object-based approach to understand what types of mangrove information can be mapped using different image datasets (Landsat TM, ALOS AVNIR-2, WorldView-2, and LiDAR). We compared and contrasted the ability of these images to map five levels of mangrove features, including vegetation boundary, mangrove stands, mangrove zonations, individual tree crowns, and species communities. We used the Moreton Bay site in Australia as the primary site to develop the classification rule sets and Karimunjawa Island in Indonesia to test the applicability of the rule sets. The results demonstrated the effectiveness of a conceptual hierarchical model for mapping specific mangrove features at discrete spatial scales. However, the rule sets developed in this study require modification to map similar mangrove features at different locations or when using image data acquired by different sensors. Across the hierarchical levels, smaller object sizes (i.e., tree crowns) required more complex classification rule sets. Incorporation of contextual information (e.g., distance and elevation) increased the overall mapping accuracy at the mangrove stand level (from 85% to 94%) and mangrove zonation level (from 53% to 59%). We found that higher image spatial resolution, larger object size, and fewer land-cover classes result in higher mapping accuracies. This study highlights the potential of selected images and mapping techniques to map mangrove features, and provides guidance for how to do this effectively through multi-scale mangrove composition mapping.
Savannas comprise a major component of the Earth system and contribute ecosystem services and functions essential to human livelihoods. Monitoring spatial and temporal trends in savanna vegetation ...and understanding change drivers is therefore crucial. Widespread greening has been identified across southern Africa; yet its drivers and manifestations on the ground remain ambiguous. This study removes the effects of precipitation on an NDVI time-series, thereby identifying trends not driven by rainfall. It utilizes the significant correlation between vegetation and precipitation as captured using MODIS and rainfall estimates. A linear regression between variables was used to derive its residual (corrected) time-series, and the rate and spatial extent of trends were evaluated in relation to biomes. A random sample-based qualitative interpretation of high spatial resolution imagery was then used to evaluate the nature of the trend on the ground. 23.25% of the country, including all biomes exhibited positive trends. We propose that greening may be related to a reduction in woody species richness, loss of the large trees and a shift towards drought tolerant shrub species, as has been shown in other sub-Saharan environments. 3.23% of the country exhibited negative trends, which were mostly associated with more humid (forested) regions pointing to deforestation as a cause; these manifested as vegetation clearing, identifiable using high resolution multi-temporal imagery. Greening trends could not be identified using this approach; instead, they point to the occurrence of gradual vegetation change caused by indirect drivers.
Display omitted
•The effects of precipitation were removed from a MODIS NDVI time-series.•Rate, intensity and spatial extent of trends were computed for each biome.•Significant greening covered 27.14% of Namibia.•Negative trends occurred mostly in the desert biome, implying land degradation.•High resolution imagery could identify negative trends as vegetation clearing.