Coral reef habitat maps derived from a high-spatial-resolution multi-spectral satellite image using object based image analysis

Coral reef maps at various spatial scales and extents are needed for mapping, monitoring, modelling, and management of these environments. High spatial resolution satellite imagery, pixel <10 m, integrated with field survey data and processed with various mapping approaches, can provide these maps. These approaches have been accurately applied to single reefs (10-100 km**2), covering one high spatial resolution scene from which a single thematic layer (e.g. benthic community) is mapped. This article demonstrates how a hierarchical mapping approach can be applied to coral reefs from individual reef to reef-system scales (10-1000 km**2) using object-based image classification of high spatial resolution images guided by ecological and geomorphological principles. The approach is demonstrated for three individual reefs (10-35 km**2) in Australia, Fiji, and Palau; and for three complex reef systems (300-600 km**2) one in the Solomon Islands and two in Fiji. Archived high spatial resolution images were pre-processed and mosaics were created for the reef systems. Georeferenced benthic photo transect surveys were used to acquire cover information. Field and image data were integrated using an object-based image analysis approach that resulted in a hierarchically structured classification. Objects were assigned class labels based on the dominant benthic cover type, or location-relevant ecological and geomorphological principles, or a combination thereof. This generated a hierarchical sequence of reef maps with an increasing complexity in benthic thematic information that included: 'reef', 'reef type', 'geomorphic zone', and 'benthic community'. The overall accuracy of the 'geomorphic zone' classification for each of the six study sites was 76-82% using 6-10 mapping categories. For 'benthic community' classification, the overall accuracy was 52-75% with individual reefs having 14-17 categories and reef systems 20-30 categories. We show that an object-based classification of high spatial resolution imagery, guided by field data and ecological and geomorphological principles, can produce consistent, accurate benthic maps at four hierarchical spatial scales for coral reefs of various sizes and complexities.

Habitat map of Danajon Bank, Philippines, derived from a high-spatial-resolution multi-spectral satellite image and georeferenced point intercept transect and spot-check survey field data, using an object-based image classification method

A mosaic of two WorldView-2 high resolution multispectral images (Acquisition dates: October 2010 and April 2012), in conjunction with field survey data, was used to create a habitat map of the Danajon Bank, Philippines (10°15'0'' N, 124°08'0'' E) using an object-based approach. To create the habitat map, we conducted benthic cover (seafloor) field surveys using two methods. Firstly, we undertook georeferenced point intercept transects (English et al., 1997). For ten sites we recorded habitat cover types at 1 m intervals on 10 m long transects (n= 2,070 points). Second, we conducted geo-referenced spot check surveys, by placing a viewing bucket in the water to estimate the percent cover benthic cover types (n = 2,357 points). Survey locations were chosen to cover a diverse and representative subset of habitats found in the Danajon Bank. The combination of methods was a compromise between the higher accuracy of point intercept transects and the larger sample area achievable through spot check surveys (Roelfsema and Phinn, 2008, doi:10.1117/12.804806). Object-based image analysis, using the field data as calibration data, was used to classify the image mosaic at each of the reef, geomorphic and benthic community levels. The benthic community level segregated the image into a total of 17 pure and mixed benthic classes.

Circumarctic wetland dataset based on ENVISAT ASAR Global Monitoring Mode with links to geotiff image and NetCDF files (2005-01 to 2011-12)

A circumpolar representative and consistent wetland map is required for a range of applications ranging from upscaling of carbon fluxes and pools to climate modelling and wildlife habitat assessment. Currently available data sets lack sufficient accuracy and/or thematic detail in many regions of the Arctic. Synthetic aperture radar (SAR) data from satellites have already been shown to be suitable for wetland mapping. Envisat Advanced SAR (ASAR) provides global medium-resolution data which are examined with particular focus on spatial wetness patterns in this study. It was found that winter minimum backscatter values as well as their differences to summer minimum values reflect vegetation physiognomy units of certain wetness regimes. Low winter backscatter values are mostly found in areas vegetated by plant communities typically for wet regions in the tundra biome, due to low roughness and low volume scattering caused by the predominant vegetation. Summer to winter difference backscatter values, which in contrast to the winter values depend almost solely on soil moisture content, show expected higher values for wet regions. While the approach using difference values would seem more reasonable in order to delineate wetness patterns considering its direct link to soil moisture, it was found that a classification of winter minimum backscatter values is more applicable in tundra regions due to its better separability into wetness classes. Previous approaches for wetland detection have investigated the impact of liquid water in the soil on backscatter conditions. In this study the absence of liquid water is utilized. Owing to a lack of comparable regional to circumpolar data with respect to thematic detail, a potential wetland map cannot directly be validated; however, one might claim the validity of such a product by comparison with vegetation maps, which hold some information on the wetness status of certain classes. It was shown that the Envisat ASAR-derived classes are related to wetland classes of conventional vegetation maps, indicating its applicability; 30% of the land area north of the treeline was identified as wetland while conventional maps recorded 1-7%.

Mosaic image of the artificial light produced by Berlin on the night of September 11, 2010 observed from an altitude of 10,000 ft

Aerial observations of light pollution can fill an important gap between ground based surveys and nighttime satellite data. Terrestrially bound surveys are labor intensive and are generally limited to a small spatial extent, and while existing satellite data cover the whole world, they are limited to coarse resolution. This paper describes the production of a high resolution (1 m) mosaic image of the city of Berlin, Germany at night. The dataset is spatially analyzed to identify themajor sources of light pollution in the city based on urban land use data. An area-independent 'brightness factor' is introduced that allows direct comparison of the light emission from differently sized land use classes, and the percentage area with values above average brightness is calculated for each class. Using this methodology, lighting associated with streets has been found to be the dominant source of zenith directed light pollution (31.6%), although other land use classes have much higher average brightness. These results are compared with other urban light pollution quantification studies. The minimum resolution required for an analysis of this type is found to be near 10 m. Future applications of high resolution datasets such as this one could include: studies of the efficacy of light pollution mitigation measures, improved light pollution simulations, economic and energy use, the relationship between artificial light and ecological parameters (e.g. circadian rhythm, fitness, mate selection, species distributions, migration barriers and seasonal behavior), or the management of nightscapes. To encourage further scientific inquiry, the mosaic data is freely available at Pangaea.

Spectral reflectance library of selected biotic and abiotic coral reef features in Glovers Reef, Belize

Underwater spectral reflectance was measured for selected biotic and abiotic coral reef features of Glovers Reef, Belize from March 6 - 10, 2005. Spectral reflectance's of 63 different benthic types were obtained in-situ. An Ocean Optics USB2000 spectrometer was deployed in an custom made underwater housing with a 0.5 m fiber-optic probe mounted next to an artificial light source. Spectral readings were collected with the probe (bear fibre) about 5 cm from the target to ensure that the target would fill the field of view of the fiber optic (FOV diameter ~4.4 cm), as well as to reduce the attenuating effect of the intermediate water (Roelfsema et al., 2006). Spectral readings included for one target included: 1 reading of the covered spectral fibre to correct for instrument noise, 1 reading of spectralon panel mounted on divers wrist to measure incident ambient light, and 8 readings of the target. Spectral reflectance was calculated for each target by first subtracting the instrument noise reading from each other reading. The corrected target readings were then divided by the corrected spectralon reading resulting in spectral reflectance of each target reading. An average target spectral reflectance was calculated by averaging the eight individual spectral reflectance's of the target. If an individual target spectral reflectance was visual considered an outlier, it was not included in the average spectral reflectance calculation. See Roelfsema at al. (2006) for additional info on the methodology of underwater spectra collection.

Digital surface model, hillshade and Lambertian reflectance model of the rock glaciers Oelgrube and Aeusseres Hochebenkar (Oetztal Alps, Tyrol, Austria)

With full-waveform (FWF) lidar systems becoming increasingly available from different commercial manufacturers, the possibility for extracting physical parameters of the scanned surfaces in an area-wide sense, as addendum to their geometric representation, has risen as well. The mentioned FWF systems digitize the temporal profiles of the transmitted laser pulse and of its backscattered echoes, allowing for a reliable determination of the target distance to the instrument and of physical target quantities by means of radiometric calibration, one of such quantities being the diffuse Lambertian reflectance. The delineation of glaciers is a time-consuming task, commonly performed manually by experts and involving field trips as well as image interpretation of orthophotos, digital terrain models and shaded reliefs. In this study, the diffuse Lambertian reflectance was compared to the glacier outlines mapped by experts. We start the presentation with the workflow for analysis of FWF data, their direct georeferencing and the calculation of the diffuse Lambertian reflectance by radiometric calibration; this workflow is illustrated for a large FWF lidar campaign in the Ötztal Alps (Tyrol, Austria), operated with an Optech ALTM 3100 system. The geometric performance of the presented procedure was evaluated by means of a relative and an absolute accuracy assessment using strip differences and orthophotos, resp. The diffuse Lambertian reflectance was evaluated at two rock glaciers within the mentioned lidar campaign. This feature showed good performance for the delineation of the rock glacier boundaries, especially at their lower parts.

Benthic and substrate cover data derived from photo-transect surveys in Heron Reef

Underwater georeferenced photo-transect survey was conducted on September 23 - 27, 2007 at different sections of the reef flat, reef crest and reef slope in Heron Reef. For this survey a snorkeler or diver swam over the bottom while taking photos of the benthos at a set height using a standard digital camera and towing a surface float GPS which was logging its track every five seconds. A standard digital compact camera was placed in an underwater housing and fitted with a 16 mm lens which provided a 1.0 m x 1.0 m footprint, at 0.5 m height above the benthos. Horizontal distance between photos was estimated by three fin kicks of the survey diver/snorkeler, which corresponded to a surface distance of approximately 2.0 - 4.0 m. The GPS was placed in a dry-bag and logged its position as it floated at the surface while being towed by the photographer. A total of 3,586 benthic photos were taken. A floating GPS setup connected to the swimmer/diver by a line enabled recording of coordinates of each benthic. Approximation of coordinates of each benthic photo was done based on the photo timestamp and GPS coordinate time stamp, using GPS Photo Link Software (www.geospatialexperts.com). Coordinates of each photo were interpolated by finding the gps coordinates that were logged at a set time before and after the photo was captured. Benthic or substrate cover data was derived from each photo by randomly placing 24 points over each image using the Coral Point Count excel program (Kohler and Gill, 2006). Each point was then assigned to 1 out of 80 cover types, which represented the benthic feature beneath it. Benthic cover composition summary of each photo scores was generated automatically using CPCE program. The resulting benthic cover data of each photo was linked to gps coordinates, saved as an ArcMap point shapefile, and projected to Universal Transverse Mercator WGS84 Zone 56 South.

TanDEM-X ice thickness dataset for Wilkins Ice Shelf, Antarctic Peninsula, with link to GeoTIFF

The dataset shows the ice thickness over Wilkins Ice Shelf, Antarctic Peninsula derived from TanDEM-X Interferometry. The data has been acquired between June and August 2012. The TanDEM-X heights have been linked to CryoSAT-2 heights (V. Helm) from the respective time stamp. Elevations have been transformed from WGS84 ellipsoidal heights to the EGM2008 geoid. The ice shelf thickness was estimated assuming hydrostatic equilibrium and a mean ice density of 915 kg/m³.

Time-series of seagrass and land cover maps from 1972 to 2010, in South East Queensland, Australia, with links to GeoTIFFs

Long term global archives of high-moderate spatial resolution, multi-spectral satellite imagery are now readily accessible, but are not being fully utilised by management agencies due to the lack of appropriate methods to consistently produce accurate and timely management ready information. This work developed an object-based remote sensing approach to map land cover and seagrass distribution in an Australian coastal environment for a 38 year Landsat image time-series archive (1972-2010). Landsat Multi-Spectral Scanner (MSS), Thematic Mapper (TM) and Enhanced Thematic Mapper (ETM+) imagery were used without in situ field data input (but still using field knowledge) to produce land and seagrass cover maps every year data were available, resulting in over 60 map products over the 38 year archive. Land cover was mapped annually using vegetation, bare ground, urban and agricultural classes. Seagrass distribution was also mapped annually, and in some years monthly, via horizontal projected foliage cover classes, sand and deep water. Land cover products were validated using aerial photography and seagrass maps were validated with field survey data, producing several measures of accuracy. An average overall accuracy of 65% and 80% was reported for seagrass and land cover products respectively, which is consistent with other studies in the area. This study is the first to show moderate spatial resolution, long term annual changes in land cover and seagrass in an Australian environment, created without the use of in situ data; and only one of a few similar studies globally. The land cover products identify several long term trends; such as significant increases in South East Queensland's urban density and extent, vegetation clearing in rural and rural-residential areas, and inter-annual variation in dry vegetation types in western South East Queensland. The seagrass cover products show that there has been a minimal overall change in seagrass extent, but that seagrass cover level distribution is extremely dynamic; evidenced by large scale migrations of higher seagrass cover levels and several sudden and significant changes in cover level. These mapping products will allow management agencies to build a baseline assessment of their resources, understand past changes and help inform implementation and planning of management policy to address potential future changes.

Bathymetric map of Heron Reef, Australia, derived from airborne hyperspectral data at 1 m resolution

A simple method for efficient inversion of arbitrary radiative transfer models for image analysis is presented. The method operates by representing the shape of the function that maps model parameters to spectral reflectance by an adaptive look-up tree (ALUT) that evenly distributes the discretization error of tabulated reflectances in spectral space. A post-processing step organizes the data into a binary space partitioning tree that facilitates an efficient inversion search algorithm. In an example shallow water remote sensing application, the method performs faster than an implementation of previously published methodology and has the same accuracy in bathymetric retrievals. The method has no user configuration parameters requiring expert knowledge and minimizes the number of forward model runs required, making it highly suitable for routine operational implementation of image analysis methods. For the research community, straightforward and robust inversion allows research to focus on improving the radiative transfer models themselves without the added complication of devising an inversion strategy.

Multi-temporal mapping of seagrass cover, species and biomass of the Eastern Banks, Moreton Bay, Australia, with links to shapefiles

The spatial and temporal dynamics of seagrasses have been studied from the leaf to patch (100 m**2) scales. However, landscape scale (> 100 km**2) seagrass population dynamics are unresolved in seagrass ecology. Previous remote sensing approaches have lacked the temporal or spatial resolution, or ecologically appropriate mapping, to fully address this issue. This paper presents a robust, semi-automated object-based image analysis approach for mapping dominant seagrass species, percentage cover and above ground biomass using a time series of field data and coincident high spatial resolution satellite imagery. The study area was a 142 km**2 shallow, clear water seagrass habitat (the Eastern Banks, Moreton Bay, Australia). Nine data sets acquired between 2004 and 2013 were used to create seagrass species and percentage cover maps through the integration of seagrass photo transect field data, and atmospherically and geometrically corrected high spatial resolution satellite image data (WorldView-2, IKONOS and Quickbird-2) using an object based image analysis approach. Biomass maps were derived using empirical models trained with in-situ above ground biomass data per seagrass species. Maps and summary plots identified inter- and intra-annual variation of seagrass species composition, percentage cover level and above ground biomass. The methods provide a rigorous approach for field and image data collection and pre-processing, a semi-automated approach to extract seagrass species and cover maps and assess accuracy, and the subsequent empirical modelling of seagrass biomass. The resultant maps provide a fundamental data set for understanding landscape scale seagrass dynamics in a shallow water environment. Our findings provide proof of concept for the use of time-series analysis of remotely sensed seagrass products for use in seagrass ecology and management.

Surface temperature measurements near Ny-Ã…lesund, Svalbard, Norway in 2008 with links to data files

The ground surface temperature is one of the key parameters that determine the thermal regime of permafrost soils in arctic regions. Due to remoteness of most permafrost areas, monitoring of the land surface temperature (LST) through remote sensing is desirable. However, suitable satellite platforms such as MODIS provide spatial resolutions, that cannot resolve the considerable small-scale heterogeneity of the surface conditions characteristic for many permafrost areas. This study investigates the spatial variability of summer surface temperatures of high-arctic tundra on Svalbard, Norway. A thermal imaging system mounted on a mast facilitates continuous monitoring of approximately 100 x 100 m of tundra with a wide variability of different surface covers and soil moisture conditions over the entire summer season from the snow melt until fall. The net radiation is found to be a control parameter for the differences in surface temperature between wet and dry areas. Under clear-sky conditions in July, the differences in surface temperature between wet and dry areas reach up to 10K. The spatial differences reduce strongly in weekly averages of the surface temperature, which are relevant for the soil temperature evolution of deeper layers. Nevertheless, a considerable variability remains, with maximum differences between wet and dry areas of 3 to 4K. Furthermore, the pattern of snow patches and snow-free areas during snow melt in July causes even greater differences of more than 10K in the weekly averages. Towards the end of the summer season, the differences in surface temperature gradually diminish. Due to the pronounced spatial variability in July, the accumulated degree-day totals of the snow-free period can differ by more than 60% throughout the study area. The terrestrial observations from the thermal imaging system are compared to measurements of the land surface temperature from the MODIS sensor. During periods with frequent clear-sky conditions and thus a high density of satellite data, weekly averages calculated from the thermal imaging system and from MODIS LST agree within less than 2K. Larger deviations occur when prolonged cloudy periods prevent satellite measurements. Futhermore, the employed MODIS L2 LST data set contains a number of strongly biased measurements, which suggest an admixing of cloud top temperatures. We conclude that a reliable gap filling procedure to moderate the impact of prolonged cloudy periods would be of high value for a future LST-based permafrost monitoring scheme. The occurrence of sustained subpixel variability of the summer surface temperature is a complicating factor, whose impact needs to be assessed further in conjunction with other spatially variable parameters such as the snow cover and soil properties.

Forest disturbance map for the Carpathian ecoregion, 1984-2012

Detailed knowledge of forest cover dynamics is crucial for many applications from resource management to ecosystem service assessments. Landsat data provides the necessary spatial, temporal and spectral detail to map and analyze forest cover and forest change processes. With the opening of the Landsat archive, new opportunities arise to monitor forest dynamics on regional to continental scales. In this study we analyzed changes in forest types, forest disturbances, and forest recovery for the Carpathian ecoregion in Eastern Europe. We generated a series of image composites at five year intervals between 1985 and 2010 and utilized a hybrid analysis strategy consisting of radiometric change classification, post-classification comparison and continuous index- and segment-based post-disturbance recovery assessment. For validation of the disturbance map we used a point-based accuracy assessment, and assessed the accuracy of our forest type maps using forest inventory data and statistically sampled ground truth data for 2010. Our Carpathian-wide disturbance map achieved an overall accuracy of 86% and the forest type maps up to 73% accuracy. While our results suggested a small net forest increase in the Carpathians, almost 20% of the forests experienced stand-replacing disturbances over the past 25 years. Forest recovery seemed to only partly counterbalance the widespread natural disturbances and clear-cutting activities. Disturbances were most widespread during the late 1980s and early 1990s, but some areas also exhibited extensive forest disturbances after 2000, especially in the Polish, Czech and Romanian Carpathians. Considerable shifts in forest composition occurred in the Carpathians, with disturbances increasingly affecting coniferous forests, and a relative decrease in coniferous and mixed forests. Both aspects are likely connected to an increased vulnerability of spruce plantations to pests and pathogens in the Carpathians. Overall, our results exemplify the highly dynamic nature of forest cover during times of socio-economic and institutional change, and highlight the value of the Landsat archive for monitoring these dynamics.