Mangrove Forests Mapping in the Southern Part of Japan Using Landsat ETM + with DEM

A regional map of mangrove forests was produced for six islands located in the southern part of Japan by integrating the spectral analyses of Landsat Enhanced Thematic Mapper plus (ETM+) images with a digital elevation model (DEM). Several attempts were applied to propose a reliable method, which can be used to map the distribution of mangrove forests at a regional scale. The methodology used in this study comprised of obtaining the difference between Normalized Difference Water Index (NDWI) and Normalized Difference Vegetation Index (NDVI), band ratio 5/4, and band 5, from Landsat ETM+, and integrating them with the topographic information. The integration of spectral analyses with topographic data has clearly separated the mangrove forests from other vegetation. An accuracy assessment was carried out in order to check the accuracy of the results. High overall accuracy ranging from 89.3% to 93.6% was achieved, which increased the opportunity to use this methodology in other countries rich in mangrove forests.


Introduction
Mangrove forests are of great environmental and land cover-mapping importance, they extend throughout tropical and sub-tropical ecosystems.Mangroves have important aesthetic, historical and cultural values [1].In addition, mangroves have significant economical and environmental values beside their influence on global climate change.They protect coastlines against high tidal waves and floods, host a number of animal species, provide food for many marine species, and act as filters of sediments and nutrients in coastal catchments [2].
On the other hand, mangroves have strong relationships with the surrounding environments, since the occurrence of mangrove species at a certain location is related to the surrounding ecological gradients such as elevation, tidal inundation, water salinity, and soil pH [3].Within the same environment there are other types of vegetation such as pastures, coastal shrubs and forests that grow up alongside with mangroves.The problem with these other types of vegetation is that they have spectral similarity with mangroves at a finer level [4].Thus, it is anticipated that producing maps of mangrove forests at a regional scale will usually require high spatial resolution remotely sensed data.
The importance of remote sensing for mangrove forests mapping is well established.Many papers have produced to analyze the change detection of mangrove forests using aerial photographs [5], and to monitor coastal ecosystem changes supported by multi-temporal and multi-spatial resolution remote sensing data [6].In addition, several studies have been proposed for mapping mangroves at local scale using different spatial resolution satellite data, including Landsat Thematic Mapper (TM) [7], Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) [3], Satellite Pour l'Observation de la Terre (SPOT) [1], Synthetic-aperture radar (SAR) [8], Compact Airborne Spectrographic Imager (CASI) [2], IKONOS [9] and Quickbird [10].On the other hand many methods have been proposed for identifying and mapping mangrove forests from remote sensing data.For instance, [3] proposed a Bayesian probability method incorporating ecological parameters to develop a post classifier for mangrove forests that can be used as a guideline for producing mangrove maps at a finer level.[11] used principle component analysis and band ratios in addition to other techniques to compare the relative effectiveness of multi-spectral remote sensing data in mapping mangrove forests in Turks and Caicos Islands, and [12] compared different classifiers for SPOT XS and Terra ASTER data.More recently, [9] applied an ISO-DATA unsupervised classification technique on IKONOS and QuickBird images to distinguish different species of mangrove forests in Mabala and Yélitono mangrove islands of Guinea, West Africa.
Various studies, however, are primarily based on supervised classifications.For example, [4] mapped mangrove forests in the western Waitemata Harbour, Auckland, New Zealand by applying the maximum likelihood classification method (MLC) on Landsat TM and SPOT XS imagery.This study reported that supervised classification techniques can lead to more useful results when it is used for mapping mangrove forests.[13] assessed mangrove forests in the Mekong Delta in Vietnam through applying vegetative index differencing, and supervised classification method to SPOT XS. [14] monitored mangrove forest of Sundabans, the border between Bangladesh and India, by executing the MLC method on multi-temporal Landsat data.
The World Mangrove Atlas provided a realistic assessment of existing mangroves and their evolutionary trends at both global and local scales [15].It is the most widely used reference to map the global distribution of mangroves.Moreover, two global land cover projects used to map mangrove forest independently as a separate land cover class using different sources of remotely sensed data: Global Land Cover 2000 (GLC 2000) [16], and Global Land Cover by National Mapping Organizations (GLCNMO) [17].
In this context, the aim of this study is to develop an efficient method to extract mangrove forests at a regional scale using remote sensing data.The specific objectives are: using Landsat Enhanced Thematic Mapper plus (ETM+) together with a Digital Elevation Model (DEM): 1) to discriminate between mangrove and non-mangrove forests; 2) to assess the capability of data used for mapping mangroves distribution; 3) and to map the extent of mangrove forests in Japan as a case study of a regional scale.

Study Area
The study area is located in the southern part of Japan.It extends from Kyushu Island (32˚28'N, 131˚02'E) to Iriomote Island (24˚20'N, 123˚49'E), one of the Ryukyu Islands, Japan, as shown in Figure 1.The study area comprises six islands: Kyushu, Tanega, Amami-O, Okinawa, Ishigaki and Iriomote.Climatically, these islands lie in the transition from subtropical to tropical conditions.They include a diverse array of natural subtropical forests including broadleaf evergreen and deciduous forests, needle leaf evergreen and deciduous forests, mixed forest, tree open and mangrove forest in addition to agricultural land [17].Mangrove forests can be found in some of these islands in a variety of settings as a coastal fringe, and along some rivers.

Satellite Images and Reference Data
The Landsat ETM+ and DEM acquired by the Shuttle Radar Topography Mission (SRTM) of the National Geospatial Intelligence Agency (NGA) 2010 were the main satellite data used in this study.These data are available at the Global Land Cover Facility (GLCF) [18] and the US Geological Survey (USGS) Center for Earth Resources Observation and Science (EROS) [19].
In this study, it was difficult to obtain all the Landsat ETM+ images at the same date, because some images were very cloudy in which mangrove forests were difficult to be observed.Thus, other cloud-free images at different dates were used.Moreover, two Landsat ETM+ images, acquired in 2004 and 2007, were used in this study.In these images, there was stripe noise due to an error in the Scan Line Corrector (SLC) of Landsat ETM+.However, the mangrove forest locations were not affected by the SLC error, and the images were still usable.Table 1 shows the information related to the satellite images used in this study for all islands.
Various reference data were used to validate the results of this study.The reference data mainly include: 1) aerial photographs at a nominal scale of 1:20,000 obtained in 1987; 2) vegetation maps published by the Natural Environment Biodiversity Center, Ministry of the Environment, at scales 1:50,000 and 1:25,000 acquired in 2001-2006, respectively; 3) detailed topographic maps at scale 1:25,000 published during the period from 2004 to 2006.

Image Processing
A subset of each of Landsat ETM+ and DEM image covering the study area (6 Islands) were used.The images were resampled to a pixel size of 30 m× 30 m using the nearest neighbour method.The Landsat ETM+ images were in the Universal Transverse Mercator (UTM) projection (Zones 51 and 52) and the WGS84 ellipsoid, as shown in Table 1.For each Landsat ETM+ image, a false color composite image was generated from bands 4, 5, and 7 and displayed in red, green and blue (RGB) colors, respectively.The false color composite images were used for the purpose of visual interpretations and analyses.
Since mangroves are an intensive mixture between two land cover classes that are the forest and water classes, the water was masked out using the infrared bands 4, 5 and 7 in all Landsat ETM+ images.

Spectral Analysis of Mangrove Forests
Mangroves are a typical example of an evergreen forest [20], situated in the intertidal zone along tropical coastlines of saline water and can be found along riverbanks.It is easy to differentiate among mangrove forests and other evergreen inland forests when doing a field survey.However, it maybe difficult to differentiate between these forests using remote sensing data, especially when mangrove forests are not covering a wide area, as was the case for this research study.Therefore, it is important to realize the environmental conditions in which mangrove forests live, to be able to classify this vital land cover, as a separate class.As stated earlier, mangrove forests exist usually along the coastline, which means in low elevation areas.Thus, it is very important to include the topographic condition as one factor when mapping the mangrove forests.Another important factor is that the mangrove forests exist in an aquatic environment, which indicates that water is a predominant factor.These two factors are considered common factors of all mangrove forests, which are not restricted to specific species of mangrove.These two factors are the basis of mapping mangrove forests which can be used not only at a local scale but also at a regional scale.Therefore, to develop an efficient method for mapping mangroves, the spatial distribution of mangrove forests and their spectral characteristics should be taken into consideration.Specifically, normalized difference water index (NDWI), normalized difference vegetation index (NDVI), and topographic conditions should be considered, in order to maximize the contrast between mangrove and non-mangrove forests.

NDWI and NDVI
The Normalized Difference Water Index (NDWI) proposed by [21] was used to assess the vegetation water content using remote sensing data based on physical principles.NDWI is sensitive to changes of water content of vegetation canopies.It is considered as an independent vegetation index that was developed to delineate vegetation water content features, and to enhance their presence in remotely-sensed digital imagery.NDWI is expressed using the following: where NIR is the reflectance or radiance in a near infrared wavelength channel (0.78 -0.90 μm), and SWIR is the reflectance or radiance in a short wave infrared wavelength channel (1.55 -1.75 μm).NIR and SWIR correspond to bands 4 and 5 for ETM+ images, respectively.The NIR reflectance is affected by leaf internal structure and leaf dry matter content.Therefore, this channel is located in the high reflectance plateau of vegetation canopies [21], and it is frequently used for vegetation analysis.On the other hand, the SWIR reflectance reflects changes in the vegetation water content, and it is sensitive to water conditions in leaves as well as on the topsoil [13].NDWI is one of the indices that contrast the SWIR channel with the NIR channel that are sensitive to the mass or volume of water and not to the fractional percentage of water [22].
Normalized Difference Vegetation Index (NDVI) was used widely to monitor the quality and distribution of vegetation.NDVI was, in fact, often referred to as a greenness index.This index can be computed with a simple formula [23].
where Red is the reflectance or radiance in a visible wavelength channel (0.63 -0.69 μm), and corresponds to band 3 for ETM+ images.
In this article, NDVI images were generated to enhance mangrove forests that have higher NIR reflectance, and lower red light reflectance.Also, NDVI images were produced to eliminate water bodies, those of low red light reflectance, and those of very low NIR reflectance.
Since NDWI and NDVI were used to detect the vegetation water content, and to monitor the distribution of vegetation, respectively; the difference between these two indices successfully emphasized the existence of aquatic vegetation on the images [13], which implies positive illustration of mangrove forests locations (Figure 2(a)).

Band Ratio
The band ratio was used to enhance the spectral differences between bands, and to reduce the effects of topography.Dividing one spectral band by another can produce an image that provides relative band intensities.In the resultant image of band ratio, certain ground features tend to be highlighted based on unusual or anomalous ratio values, and also eliminate the effects of dark shadow.In this present study, band ratio 5/4 was used to improve the difference between mangroves and nonmangrove forests.The selection of this band ratio was based on the physical properties and canopy spectra of mangrove forests, which had already been used for mapping mangrove forests in Belize [11].

Band 5
The Band 5 of Landsat ETM+ was frequently applied to monitor vegetation moisture content.This band has wavelength (1.55 -1.75 µm), and is located in the short wave infrared, as mentioned earlier.The spectral profile of band 5 showed a lower digital number (DN) value within mangrove forest compared with other classes [24].This makes band 5 very useful for distinguishing mangrove forests from other vegetations.In band 5, the mangrove forests appeared in dark grey (Figure 2(c)).

Digital Elevation Model
The classification of a study site can be improved with the addition of a digital elevation model (DEM).A DEM is an important factor for classifying forest vegetation distributed by elevation [25].It represents the correlation among the cover types such as forest and habitats located in the study area.Thus, the result of classification can be improved when adding the DEM rather than the original classification method, because it can increase the distinction among the forests and other vegetation species, based on elevation values.
A DEM of each island, acquired by SRTM, was employed as an important criterion for classifying mangrove forests based on their DN values that represent the elevation.A DEM was integrated with the spectral data after being resampled from 90 m to 30 m spatial resolution, to improve the visual discrimination between mangroves and other vegetation located in the study area and to exclude non-mangrove pixels that had similar spectral attributes with mangrove pixels, but were located above the elevation limiting line [26].
A DEM clearly separated the mangrove forests that were located within lowland areas (along the shoreline and rivers) from other forests and vegetation located in highland areas.

Layer Stacking
The four parameters: the difference between NDWI and NDVI, band ratio 5/4, band 5 of Landsat ETM+, and the DEM, were merged and stacked together in one image resulting in the final image containing three spectral variables and one topographic variable.This final image was used to classify and map the real extent of mangrove forests in each image for each island, separately, as shown in Figure 2(d).

Selection of Training Sites and Classification
A supervised classification approach was adopted to classify the mangrove forests.The maximum likelihood algorithm was performed using carefully selected training areas.The training areas were delineated very carefully for each image to represent the real extent of mangrove forests located on various elevations, regardless of their species.
The training samples were drawn in order to get the spectral signature of mangrove forests and other different land use/cover classes in each image, separately.Thereby, a supervised classification technique using maximum Likelihood classifier was subsequently applied on each image to map the mangrove forests and other land use/ cover classes.Post-classification editing was done to merge all the different land use/cover classes, except the mangrove forests class, into one class called the nonmangrove forests class.Thus, the final classified images contain two classes: the mangrove forest class and the non-mangrove forest class.A median convolution filter, kernel size 3  3, was used to smooth the classified images and to remove isolated pixels.

Accuracy Assessment
This study relied on 1:25,000 and 1:50,000 vegetation maps to evaluate the accuracy of the classified mangrove forests obtained from Landsat ETM+ images.However, it was necessary to confirm their reliability first before using them as references.Therefore, other reference data were used in this study such as aerial photographs and detailed topographic maps to validate the vegetation maps.After validating the vegetation maps, they were used as references to check the accuracy of the classified Landsat images.
From the histogram analysis, it was found Mangrove forests in the study area distribute within 1 km from the shoreline and lower than 20 m above sea level; therefore, the random samples were distributed within the selected threshold values mentioned above.For Okinawa, Ishigaki and Iriomote Islands, a total of 250 random sample points were selected (100 points for mangrove forests class and 150 points for non-mangrove forests class), while a total of 150 random sample points were selected for Tanega and Amami-O Islands (50 points for mangrove forests class and 100 points for non-mangrove for-ests class).The number of sample points was decided based on the area of mangrove forests in each island.Table 2 shows the calculated overall accuracy, user and producer accuracies, and kappa coefficient of the classified mangrove forests class in Tanega, Amami-O, Okinawa, Ishigaki and Iriomote islands, respectively.

Results and Discussion
A total of 7.67 km 2 of mangrove forests was calculated within the study area.Table 3 presents a comparison between the land area of mangrove forests calculated from both of the classified Landsat images, and vegetation maps for each island.
The results of this study were compared with a number of existing global land cover maps such as GLC2000, GLCNMO and the world mangrove atlas, and regional maps such as the reference vegetation maps, to evaluate the relative performance of this methodology.It is worthy of mention that the GLC2000, which used daily 1-km SPOT/VEGETATION data of 2000, did not show any mangrove forest in the study area.While in GLCNMO of  2010, the mangrove forests were mapped as a separate class.In GLCNMO data, the mangrove forests were extracted by visual interpretation of Landsat ETM+ images, and by digitizing the mangrove forests available in the world mangrove atlas.In other words, the mangrove forests in both GLCNMO and world mangrove atlas were identical, for the study area.A detailed analysis of each island is discussed, separately.

Kyushu Island
The vegetation map showed that the land area of mangrove forests in this island was about 0.05 km 2 .Whereas, the results indicated that there were no mangrove forests in Kyushu Island.The use of a Landsat image to extract this small patch of mangrove was insufficient.Higher spatial resolution data such as, Compact Airborne Spectrographic Imager (CASI), Quickbird, and IKONOS images might be able to extract this small area of mangrove forests located in Kyushu Island.However, the expense of mapping an area equivalent to the coverage of this small patch would be considerable.

Tanega Island
The mangrove forests covered about 0.33 km 2 when calculated from the classified Landsat image, while the land area of mangrove forests calculated from the vegetation map was about 0.28 km 2 .This implies that there was a consistency between the mangrove forests of the classified Landsat data and the reference vegetation data (Figure 3).

Amami-O Island
Based on the vegetation map, the land area of mangrove forests in this island was approximately 0.56 km 2 .While only 0.34 km 2 of mangrove forests was calculated from the classified Landsat image.The difference of 0.22 km 2 of land area might be due to insufficiency of Landsat spatial resolution to extract very small areas of mangrove forests, which required finer spatial resolution data.Figure 4 illustrates a comparison of mangrove forests resulting from the classified Landsat image and the vegetation map.

Okinawa Island
Mangrove forests in Okinawa Island covered an area of 0.26 km 2 according to the reference vegetation map.Whereas, the calculated area of mangroves forests from the classified Landsat image was 0.20 km 2 .This means that the results of the classification were quite close to the land area calculated from the reference vegetation map (Figure 5).Copyright © 2013 SciRes.JGIS

Ishigaki Island
The Mangrove forests on Ishigaki Island covered a larger area compared with the above mentioned islands (Figure 6).The calculated land area of mangrove forests from the classified Landsat image was about 0.8 km 2 , which was very close to the land area of 1.0 km 2 calculated from the reference vegetation maps.

Iriomote Island
The largest area covered by mangrove forests was concentrated in Iriomote Island, as presented in Figure 7.
The land area of mangrove forests mapped from the Landsat image was about 6.00 km 2 , while in the referenced vegetation map, it was about 6.2 km 2 .
The comparison results between the classified Landsat images, and the reference vegetation maps provided a great support to the methodology used to extract the mangrove forests.Despite the fact that most of the land areas were less than 1 km 2 (Table 3) with the exception of Iriomote Island, the spectral analysis of Landsat images together with the DEM data were able to detect the mangrove forests clearly.
The methodology used in this study was developed after several attempts that were applied on Iriomote Island, as an example, since it has the largest area of mangrove forests in the study area, in order to get the most accurate method to illustrate the real extent of mangrove forests in the study area.One of these methods was to use the spectral data only, without using the DEM.Thus, the result was not good when it was compared with the vegetation maps, since this method was not able to separate the mangrove forests from other forests that have similar spectral characteristics located at higher elevations.The results of this method showed an overall accuracy of 74.4% and a Kappa coefficient of 0.41 as  shown in Table 2.Then, a DEM was applied through using an elevation threshold value based on [27].However, this did not lead to accurate results, which showed a larger area of mangrove forest than the actual situation, due to the interference of other surrounding vegetation that usually grows alongside mangroves.Using this method the overall accuracy was 69.2%, and the kappa coefficient was 0.27 (Table 2).Therefore, it was decided to incorporate the DEM without assigning any threshold value.
Another method was applied by setting a threshold value to the difference between NDWI and NDVI.The results showed that whether a threshold value was set or not, the mangrove mapping results did not change.Therefore, it was determined to use the difference between NDWI and NDVI without setting a threshold value.
Other attempts were used to check which were the best Landsat ETM+ bands to increase the contrast between mangrove and non-mangrove forests.Until it was decided that band 5 was the best among the other bands.The use of band 5 of Landsat ETM+ facilitated the selection of training areas that were used after that for classification.
It is well known that serious confusion occurs when using conventional methods in mangrove classification, because some types of land cover (e.g., agricultural land and forests) have similar spectral behaviors and distribu-tion features to mangroves [27].Therefore, further processes were applied to test the capability of integrating the difference between NDWI and NDVI to map the mangrove forests.In the first process, NDWI and NDVI images were merged with the band ratio 5/4, band 5 and the DEM, representing an image composed of five parameters.In the second process, only the NDWI image was merged and integrated with the other parameters without using the NDVI image.In the third process, the NDVI image was used alone, and merged with the other parameters without using the NDWI image.These three processes were tested on Iriomote Island, as an example.The classification results of the three processes showed that the area covered by mangrove forests was less than the actual by 4 km 2 when compared with the reference vegetation map, due to the interference of the mangrove forest class with the other land cover classes when using just NDWI and NDVI.
Accordingly, the spectral characteristics that were obtained from Landsat data, based on the physical properties of mangrove forests, in combination with topographic information, increased the accuracy of mapping the mangrove forests.In particular, the combination among the difference between NDWI and NDVI, band ratio 5/4, the band 5 of Landsat ETM+ and the DEM, has increased the accuracy of mapping the mangrove forests.Using this method the mangrove forests can be separated clearly from other forests located within the same area.

Conclusions
The objective of this study was to develop a method that is capable of mapping mangrove forests at a regional scale as an independent land cover class using remote sensing techniques.In this research, the mangrove forests were accurately mapped by integrating the difference between NDWI and NDVI, band ratio 5/4, and band 5, obtained from Landsat ETM+ with the topographic information.This realistic integration clearly separated the mangrove forests from other vegetation types.High accuracy results have already been achieved, leaving little opportunity for further improvements.
Although the total area of mangrove forests located in Japan is small, this methodology was able to map this important land cover class successfully with high accuracy.The production of regional mangrove maps in Japan supported the use of ancillary topographic data over the solo use of spectral data.This can increase the possibility of employing this method presented in this research to be applied as a guideline for producing mangrove maps at regional scales.Also, it is hypothesized that this technique can be applied in other countries that are rich in mangrove forests.
This study can be important for monitoring the changes of mangrove forests over time, and for mapping regional mangrove forests accurately and cost-effectively.However, this method can not be used to discriminate between different species of mangroves.Moreover, the frequent existence of clouds in tropical and subtropical areas is considered the main difficulty of using optical remote sensing data, which have to be taken into consideration when producing regional maps of mangrove forests.

Figure 1 .
Figure 1.Location map of the study area.

Figure 2 (
b) illustrates the output image of a band ratio derived from Landsat ETM+, and it can easily recognize mangrove forests, which are shown in dark grey.

Figure 2 .
Figure 2. The band analyses used to extract the mangrove forests, an example in Iriomote Island.(a) The difference between NDWI and NDVI image showed the mangrove forest in light grey color; (b) The band ratio 5/4 image showed the mangrove forest in dark grey color; (c) The mangrove forest appeared in dark grey color using band 5 of Landsat ETM+ image; (d) A color composite image of the difference between NDWI and NDVI, band ratio 5/4 image, and band 5 of Landsat ETM+ exposed through RGB color, respectively.

Figure 3 .
Figure 3.A comparison between (a) mangrove forests obtained from classifying a Landsat image, and (b) mangrove forests obtained from reference vegetation map, in Tanega Island.

Figure 4 .
Figure 4.A comparison between (a) mangrove forests obtained from classifying a Landsat image, and (b) mangrove forests obtained from reference vegetation map, in Amami-O Island Island.

Figure 6 .
Figure 6.A comparison between (a) mangrove forests obtained from classifying a Landsat image, and (b) mangrove forests obtained from reference vegetation map, in Ishigaki Island.

Figure 7 .
Figure 7.A comparison between (a) mangrove forests obtained from classifying a Landsat image, and (b) mangrove forests obtained from reference vegetation map, in Iriomote Island.

Table 2 . Results of accuracy assessment of the classified maps produced from Landsat ETM+.
Results of accuracy assessment of the methodology applied without using the DEM; 2 Results of accuracy assessment of the methodology after applying a threshold value to the DEM.