Abstract

Coastal wetlands of eastern Georgian Bay provide critical habitat for a variety of wildlife, especially spawning and nursery habitat for Great Lakes fishes. Although the eastern shoreline has been designated a World Biosphere Reserve by UNESCO, a complete inventory is lacking. Prior effort by the Great Lakes Coastal Wetland Consortium (GLCWC) was unable to fully identify coastal wetland habitat in eastern Georgian Bay due to limited data coverage. Here we outline the methodology, analyses, and applications of the McMaster Coastal Wetland Inventory (MCWI) created from a comprehensive collection of satellite imagery from 2002–2008. Wetlands were manually delineated in a GIS as two broad habitat types: coastal marsh and upstream wetland. Coastal marsh was further subdivided into low marsh (LM; permanently inundated) and high marsh (HM; seasonally inundated) habitat. Within the coastal zone of eastern and northern Georgian Bay there are 12629 distinct wetland units comprised of 5376 ha of LM, 3298 ha of HM and 8676 ha of upstream habitat. The MCWI identifies greater total wetland area within the coastal zone than does the GLCWC inventory (17350 ha versus 3659 ha resp.). The MCWI provides the most current and comprehensive inventory of coastal wetlands in eastern Georgian Bay.

1. Introduction

Wetlands represent some of the most biologically diverse ecosystems on the planet, yet globally, estimates of wetland loss due to human development range from 50% to 90% [1]. Canada has approximately 25% of the world’s wetlands [2, 3]. As a signatory in 1981 of the Ramsar convention, Canada has an obligation to identify and protect ecologically important wetlands (http://www.ramsar.org/). To date, Canada has protected more wetland area than any other country, but the prevalence of wetlands in the Canadian landscape means that there are many wetlands that have not yet been delineated. In many regions of Canada, there is still an urgent need to catalogue and monitor wetlands [3].

The Laurentian Great Lakes, shared by Canada, and the United States, represent the largest freshwater resource in the world. A review of wetland research conducted in this region by Herendorf [4] identified over 1500 large coastal wetlands with a total surface area of 1700 km2. These marshes provide many important ecosystem services including water purification, nutrient sequestration, and shoreline buffering [1, 3], as well as important feeding and nursery habitat for a wide variety of organisms including fish, birds, invertebrates, mammals, and reptiles [57]. In Ontario, majority of these wetlands have been lost or degraded as a result of human disturbance [8], except in the eastern and northern shore of Georgian Bay, where they are mostly in relatively pristine condition [9].

During the early 2000s, the Great Lakes Coastal Wetland Consortium (GLCWC) [10], consisting of both US and Canadian scientists and policy makers, worked together to create a Geographic Information Systems (GISs) inventory for all five Great Lakes. The GLCWC aimed to delineate all coastal wetlands along the Great Lakes shoreline that were >2 ha in size. The Consortium included all wetlands in the Ontario Great Lakes Coastal Wetland Atlas (OGLCWA) [11], which is a GIS atlas of all wetlands that had been evaluated by the Ontario Ministry of Natural Resources using the Ontario Wetland Evaluation System (OWES) [12]; information used in this atlas dated back to 1983 and were updated with information current to 1999. The OGLCWA was complete for Lakes Superior, Erie, Ontario, as well as all connecting channels, but was incomplete for Lake Huron because it was missing some large wetlands occurring in eastern and northern Georgian Bay. Therefore, the Consortium employed experts to manually identify wetlands using available aerial photographs for this region. For eastern Georgian Bay, aerial photographs (taken in the summers between 1984 to 2000) were available at a scale of 1 : 10000 from Severn Sound to Parry Sound, and this allowed for a pixel resolution of 0.4 m. From Parry Sound to Key River, however, photos were only available at a scale of 1 : 20000 (taken in the summers between 1986 to 1996) and allowed for a more coarse pixel resolution of 0.8 m. There was only limited coverage of aerial photos for northern Georgian Bay from the Key River to MacGregor Bay.

Lack of information along the shore and in the surrounding islands of eastern Georgian Bay is a serious concern because this region holds some of the most pristine coastal marshes in the entire Great Lakes basin [9, 13]. This area is unique due to the low levels of agriculture and urban development that have allowed these wetlands to maintain the highest proportions of disturbance-intolerant fish and plant taxa within the Great Lakes coastal system [9, 1315]. Although much of the shoreline was logged prior to the 1900s, easier access to inland logging sites and slow successional growth have prevented continuous logging along the shoreline in the past century. Hence, most of the wetlands have been able to persist in relatively natural condition, free of human disturbance. Furthermore, watersheds of eastern and northern Georgian Bay have thin, nutrient-poor soils on top of Precambrian Shield, which has created naturally oligotrophic coastal wetlands with very soft water mixed with more alkaline water of Georgian Bay [16]; this creates unique geochemical characteristics that support regionally high biodiversity of aquatic plants [17]. Lastly, the complex shoreline of eastern Georgian Bay is composed of both large riverine wetlands as well as thousands of small (<2 ha) shallow rocky embayments which are protected from the strong wind and wave action that characterize the region. Therefore, Georgian Bay has an assembly of coastal wetlands that are unique in the Great Lakes basin in terms of geochemistry, biodiversity, areal cover, and abundance and can be considered reference conditions for the Great Lakes [18].

Besides the ubiquitous potential for indiscriminate human development [16], coastal wetlands of Georgian Bay have been strongly influenced by the sustained low water levels that have prevailed over the last decade (Figure 1). This trend is expected to continue, with climate change scenarios predicting a further decline in water level of Lakes Huron-Michigan by >1 m during the next 25 years [19], accompanied by reduced interannual variation [20]. Coastal wetlands are dynamic systems, where diversity of habitat and biota are maintained by a natural disturbance in the form of fluctuating water levels; in years of high water, terrestrial vegetation dies back, and, in years of low water levels, aquatic vegetation dies back [21]. Without interannual water-level variation, either the aquatic or the terrestrial vegetation would dominate at the expense of the other. The current episode of sustained low water levels would favour terrestrial vegetation at the expense of aquatic vegetation [22], which is an important component that provides critical spawning and nursery habitat for the Lake Huron fish community [23, 24].

Although water levels encountered presently are not the lowest in recorded history for Lakes Huron-Michigan, the sustained low water levels that began in 2001 and that have persisted through to 2011 have not occurred in the last 100 years, and it is difficult to predict how the wetland community will adapt to these extremes, and how the fish community, in particular, will be able to adjust to losses and gains in aquatic habitat. Hence, there is an urgent need to conduct research to determine how this trend towards lower water levels will affect the quantity and quality of wetland habitat so that they can be monitored and protected from further human activities. To aid environmental agencies and municipal planners and to enable valid extrapolation, it is important that the research be conducted on a set of randomly chosen wetlands. To date, however, there is no such inventory, because the one created by the GLCWC is incomplete for this region. There are other limitations of the GLCWC inventory (herein GLCWCI) that make it unsuitable for research on fish habitat. First, since wetlands have been delineated from photos taken in different years (1983 to 1999) and at different water levels, size of habitat zones in wetlands within the inventory cannot be directly compared because they are not standardized to one water level. Secondly, the GLCWCI excludes most of the wetlands that occur in the rocky coastal region and island archipelagos of northeastern Georgian Bay, where there could be extensive fish habitat. The last and perhaps the greatest limitation for fish ecologists is that the inventory does not distinguish between terrestrial and aquatic habitat types and it is impossible to determine the distribution of fish habitat across the region.

In this paper, we show how a comprehensive coastal wetland inventory can be created for eastern and northern Georgian Bay that is both cost-effective and suitable for use in studies of fish habitat at the scale of the entire Georgian Bay (over 4500 km of shoreline). We propose to use high-resolution IKONOS satellite imagery acquired during a 5-year period with similar water levels to ensure that wetland habitat can be directly comparable across the region and to limit the inconsistencies that can result when imagery under different water-level scenarios are used to delineate wetlands boundaries. We will also apply a simple rule to delineate coastal wetland habitat into low-marsh zone (fish habitat) and high-marsh zone (meadow habitat). The approach we develop here can be used in coastal projects of other large lakes where there is a need to monitor changes in fish habitat at a scale of an entire lake basin.

2. Methods

2.1. Satellite Imagery

IKONOS satellite images covering all of eastern Georgian Bay and parts of the North Channel were acquired by Georgian Bay Forever (GBF; http://www.georgianbayforever.org/) and licensed to McMaster University. Images covering the regions between Severn Sound and Parry Sound were acquired in July 2002, images covering up to Key River were acquired in July 2003, images covering the McGregor Bay and Bay of Islands regions were collected in July 2005, an image covering Matchedash Bay was collected in September 2005, and an image covering Beaverstone Bay was acquired in August 2008 (Figure 2). The IKONOS satellite images used for this inventory have a pan-sharpened resolution of 1 m and provide spectral information in the red, green, blue, and near-infrared wavelengths. IKONOS satellite images were not available for two regions of northern Georgian Bay, but through the Ontario Provincial QuickBird Project (2007), we were able to obtain QuickBird satellite coverage for all of the remaining gaps except for only a 10 km stretch between the French River and Beaverstone Bay (Figure 2). Like IKONOS imagery, QuickBird images provide high-resolution (60 cm) multispectral data (visible and near-infrared spectrum). In total, 8 QuickBird images were used, all of which were acquired in September 2006. In a comparison of IKONOS and QuickBird images for the purpose mapping mangroves, Wang et al. [25] found little difference in their ability to map this habitat type. We therefore concluded that the use of different image types would have a negligible impact on the quality of the final inventory.

2.2. Manual Delineation

The McMaster Coastal Wetland Inventory (MCWI) was created by manually delineating wetlands from satellite images in a GIS. The IKONOS images were initially stacked for easier use. The three visible bands (red, green, blue) of the IKONOS images were used to create a true colour image. A second image was then created through the substitution of the near-infrared (NIR) band in place of the red band (i.e., NIR, green, blue). QuickBird images were already combined into both a true colour and a near-infrared form. The NIR wavelength is a good indicator of vegetation, especially in aquatic systems [26], and, by switching between the true colour and the near-infrared images, technicians are better able to discriminate wetland vegetation from surrounding land cover. Stacked images were then imported into ArcMap 9.2 (ESRI Inc., Redlands, California, USA, 2006) in the working projection Universal Transverse Mercator (UTM), zone 17. For each wetland, the technicians traced the boundary of wetland habitats (i.e., low marsh, high marsh, upstream wetland) following specific rules regarding the upper and lower limits of each category (outlined below). The result was a single polygon for each applicable wetland category for each wetland.

2.3. Rules for Delineations/Accuracy Assessment

The Ontario Ministry of Natural Resources (OMNR) in its Ontario Wetland Evaluation Systems (OWES) [12] defines coastal wetlands as wetlands that are influenced by large water bodies and generally found within 2 km of the high water mark. The coastal zone is therefore operationally defined as land within 2 km of the shoreline and within this zone, only wetlands that are hydrologically connected via surface water to Georgian Bay are considered coastal wetlands. This 2 km coastal zone equates to roughly 177000 ha along the shores of eastern and northern Georgian Bay. Since the major focus of this inventory is to quantify fish habitat, surface hydrologic connectivity is an essential criterion. Despite the existence of many hydrologically isolated wetlands upstream of the shoreline, disconnected wetlands were not included. These wetlands occur above the high-water mark, and therefore they do not serve as current or potential Great Lakes coastal fisheries habitat. They also can be difficult to identify visually because they exist along a continuum of succession from open water to areas that are fully forested, and their delineation would have greatly prolonged the time to completion of the MCWI inventory. They were thus omitted from this inventory.

Each wetland found to be within 2 km and hydrologically connected to Georgian Bay was delineated into three habitat categories (i.e., high marsh, low marsh, upstream wetland), and the area of each polygon was calculated in ArcMap 9.2 (ESRI Inc., Redlands, California, USA, 2006). Rules for delineating boundaries of each habitat category are as follows. (i) The high-marsh (HM) category represented wetland habitat that was inundated on a seasonal basis; this area is often referred to as “wet meadow” habitat. Wet meadows provide important habitat for a variety of species including birds, reptiles, and amphibians [27]. The lower limit of the HM habitat was defined by the shoreline and was the upper limit of LM. The upper limit of the HM was the forest boundary and/or when there was change from HM to upstream habitat (swamp, bog, or fen). (ii) The low-marsh (LM) category represented portions of the wetland that were permanently inundated and essential areas for fish spawning and foraging. Upper limits of LM habitat were defined by the water’s edge and exclude meadow vegetation. The lower limit of the wetland includes submerged aquatic vegetation (SAV), which is a critical component of fish habitat. Unfortunately, the dystrophic conditions of the water in most regions of eastern Georgian Bay did not allow us to map SAV using satellite images. Therefore, the lower limit was approximated by calculating a distance that is 2.5 times the width of the emergent and/or floating vegetation zone (visible in the image) and applying this from the water’s edge along the longest axis of the wetland. This distance was reduced if the lower limit extended beyond the opening of an embayment. This lake-ward boundary is considered a conservative estimate of the maximum depth of colonization by SAV, based on dozens of underwater surveys in wetlands in this region (Midwood pers. obs.). (iii) Upstream wetland (UP) habitat corresponded to all remaining wetlands that were hydrologically connected to the bay via surface water and that occur within a 2 km buffer of the shoreline. In vast majority of cases, beaver activity created conditions that separated the upper limit of the HM from the lower limit of the UP wetland habitat. Since these beaver ponds can be seasonally connected with downstream habitat, they can act as potential fish habitat for fish communities in the affected coastal wetlands and were therefore included in the inventory. Due to the large number of UP along the eastern and northern shore of Georgian Bay, there was not enough time to delineate UP habitat in all 81 images. UP was delineated for 21 IKONOS images covering the region from Parry Sound south to Severn Sound (Figure 3). These delineations were compared to existing wetland delineations in the OMNRs [28] Ontario Base Map (OBM) that corresponds to our UP habitat. We found the OBM wetland delineations provided comparable coverage (data not shown) and therefore the OBM wetland layer was directly incorporated into the MCWI for the remaining areas (60 images). We should note, however, that some of the wetlands incorporated into the MCWI from the OBM survey may not be directly connected by surface water and the extent of this error has not yet been determined.

To standardize variations among technicians working on this project, all technicians were first trained on the same five images. Only when they achieved an acceptable level of precision (greater than 85% similarity) were they allowed to contribute to the project. To further reduce technician bias associated with discerning the lower extent of LM, only one technician was assigned to digitize this habitat category for majority of the satellite images. There was less subjectivity associated with delineations of HM and UP, and, hence, more than one technician was assigned to these habitats. Once all polygons were digitized, a single technician went through the entire data set to ensure that edges between habitat categories did not overlap and that no wetlands had been missed.

We assigned wetlands to the quaternary watersheds (acquired as OBM from [28]; Figure 4) that surrounded them using a spatial join tool in ArcMap 9.2 (ESRI Inc., Redlands, CA, USA, 2006). When a wetland occurred on the boundaries of two or more watersheds, it was assigned to the watershed that held the majority of the wetland. If it was unclear to which particular watershed a wetland should be assigned, the wetland was assigned to both. This is the reason why there is a slight discrepancy (occurred in <5% of the wetlands) between total area when all wetlands are summed without regard to watershed origin and when they are pooled after they have been sorted by quaternary watershed.

2.4. Comparison of Differences in Inventories

The shapefile for the GLCWCI was available online from the Great Lakes Commission’s website [29] and was used to conduct a comparison of the MCWI and GLCWCI. The wetland layers created in the MCWI were used to clip out portions of the GLCWCI corresponding to the Coastal and Upstream regions. Our “Coastal” zone is the same as the “Lacustrine” class defined by the GLCWCI and can be used interchangeably. We note, however, that there is no category that matches our “Upstream” portion, but in the GLCWCI shapefile, there are regions that would be defined as “Upstream,” even though they were not actually classified as such. To assess differences between inventories, we used a GIS to calculate total area of wetlands for the respective inventories. The comparison included wetlands in the MCWI and the GLCWCI along the northern and eastern shore of Georgian Bay. We excluded all wetlands <2 ha in the MCWI to make this criterion consistent with that of the GLCWCI.

3. Results

3.1. McMaster Coastal Wetland Inventory

In total, 73 IKONOS and 8 QuickBird images were digitized, and each image was delineated separately for both the high-marsh and low-marsh habitats. For upstream wetland habitat, however, only 21 IKONOS images were delineated (south of Parry Sound) with the remaining data being filled in from the OBM wetland layer (Figure 3).

Despite our best efforts, we were unable to acquire appropriate imagery to fill one small gap in northern Georgian Bay (Figure 2). We know that wetlands exist in these gaps because we have conducted field sampling there (Chow-Fraser, unpub. data). Therefore, the estimate in this document should be considered a slight underestimate of the actual amount of coastal wetland habitat in northern Georgian Bay, and future efforts should be made to fill this gap with some other satellite media of the same vintage. Relative to the remainder of the shoreline in eastern and northern Georgian Bay, this gap in imagery amounts to only a small fraction of the shoreline and should be relatively easy to update as soon as appropriate imagery has been acquired.

In total, 3771 units (414 units > 2 ha) of low marsh, 6355 units (289 units > 2 ha) of high marsh, and 2603 units (883 units > 2 ha) of upland wetland are included in the MCWI. Size of wetlands in LM habitat varied a great deal, with a mean of 1.4 (± 12.0) ha and a median of 0.3 ha. By comparison, those in HM habitat were more uniform in size, with a mean of 0.5 (± 2.2) ha and median of 0.1 ha, and those in the UP habitat were larger, with a mean of 3.3 (± 9.0) ha and a median of 1.1 ha. UP habitat covered the largest area (8676 ha), followed by LM (5376 ha) and HM (3298 ha) (Table 2; Figure 4).

Along the eastern and northern shores of Georgian Bay, there are a total of 37 quaternary watersheds ranging in size from 564 ha (Giants Tomb) to 126103 ha (French River) (Figure 5; Table 1). The largest amount of wetland habitat (2394 ha) was found in the Moon-Musquash watershed (Table 1). When sorted by different type of habitat, however, we found that the Coldwater watershed was associated with the greatest amount of LM habitat (49 units with a total area of 797 ha) (Table 3). It was surprising that this LM habitat only accounted for 3.7% of the total Coldwater watershed area, when the LM habitat in Beausoleil-Severn Island accounted for 25.6% of the total watershed area; Islands in Beausoleil-Severn were also associated with the highest percentage of HM habitat (8.7%). The Eastern Coast Islands watershed contained 1035 units of HM, for a total area of 404 ha (Table 3). The Moon-Musquash River watershed contained the greatest amount of UP habitat, with 159 units and a total area of 1708 ha (Table 3), while the MacGregor-Sampson Islands watershed had the highest percentage of UP habitat (8.6%) of all 37 quaternary watersheds.

3.2. MCWI/GLCWCI Comparison

We compared the total amount of wetlands in both the MCWI and the GLCWCI with respect to wetlands >2 ha (Table 2). There were 1586 wetland units in the MCWI, covering an area of 13267 ha; by comparison, the GLCWCI only included 696 wetland units, covering a total area of 3660. Within the Coastal zone, the MCWI included more than twice as many LM (414 versus 170) and UP (883 versus 379) units and a greater number of HM units (289 versus 234). In terms of area, however, there was almost 14 times the area of LM (4044 versus 298 ha), greater than three times the amount of HM (1842 versus 587), and more than six times the amount of UP (7381 versus 1762 ha) habitat. The greater number of wetland units and area included in the MCWI is despite the inclusion of 1014 ha of wetland area that was unique to the GLCWCI (see Figure 6).

4. Discussion

The McMaster Coastal Wetland Inventory is currently the most comprehensive inventory of coastal wetlands for eastern and northern Georgian Bay. With the completion of this inventory, we can now update the total coastal wetland habitat area for the Canadian side of the Great Lakes. With the inclusion of the nonoverlapping areas of the MCWI, the Great Lakes contain 78405 ha of coastal wetland habitat in Canada, increasing the total for Lake Huron by 47.6% to 30882 ha. We have identified all coastal marshes in a region where complete data has been lacking and that have filled an important void in the distribution of coastal wetland habitat in the Great Lakes basin. This is an important advancement in the tools available for wetland managers where basin-wide decision making is essential for the future persistence of these habitats. By accessing a large collection of satellite imagery, we have been able to fill in major gaps in the Great Lakes Coastal Wetland Consortium Inventory along eastern and northern Georgian Bay that had been noted by the authors [10]. Coastal wetlands are dynamic systems that vary in wetland size and dominant vegetation type when water levels fluctuate seasonally and annually [21, 30, 31]. In order to create an inventory that provides consistent wetland coverage, it is essential that all delineations utilize imagery acquired during similar water-level conditions. If there is any temporal discrepancy in image acquisition, wetlands in images that were acquired during low water levels may have more HM habitat [31] than the same wetlands delineated using imagery acquired during higher water levels. All wetlands in the MCWI were digitized from high-resolution IKONOS and QuickBird satellite imagery acquired during a period of low water levels between 2002 and 2008 (mean of 176.12 m ± 0.13). This means that all wetland areas in our inventory are standardized and are directly comparable. Conversely, wetlands in the GLCWCI were digitized from aerial photos or satellite images acquired at different years (from 1983 to 1999; mean of 176.78 ± 0.29) and were not standardized to a consistent water level. Although wetlands digitized in the GLCWCI may provide information on wetland location and some UP boundaries, regional comparisons of LM and HM wetland area may not be possible due to the potential influence of interannual water level variation and resulting vegetation changes.

The MCWI and the GLCWCI also differ with respect to the level of detail provided by each inventory for eastern and northern Georgian Bay. The GLCWCI only identified 21% of wetlands (by area) available in the MCWI. The MCWI raises some important issues concerning the minimum mapping unit of inventories. The GLCWCI used a minimum mapping unit of 2 ha to identify coastal wetland habitat. In our objectives for the MCWI, we decided to identify all coastal marsh habitat possible with the resolution of our imagery. This proved to identify a unique characteristic of this region in that a majority of the coastal wetlands are <2 ha. A minimum mapping unit of <2 ha would exclude a large portion of the data in this project which was meant to identify critical fish habitat. Accordingly, we suggest that the GLCWCI should not be used to estimate the amount of coastal fish habitat in Georgian Bay because the amount of LM habitat in the MCWI was 13 times higher for wetlands >2 ha and 18 times higher for all wetlands regardless of size. As outlined previously, the major differences between inventories reflect how they were created. Authors of the GLCWCI indicated that they relied heavily on OWES-identified wetlands (minimum size of 2 ha) and availability of aerial photography to fill the considerable gaps in eastern and northern Georgian Bay. Unfortunately, there were more data available for upstream wetlands than for those in the coastal zone because the initial acquisitions of aerial photographs were for forest survey purposes. As a result, the authors also recognized that there were major gaps in coverage for the GLCWCI and that the missed areas likely contained a considerable number of coastal wetlands [10]. In this respect, the MCWI should have the most complete coverage given that the satellite images we acquired provided a seamless coverage of the entire 2 km coastline of eastern and most of northern Georgian Bay.

Herdendorf [4] identified approximately 1500 wetlands in the Great Lakes (total area of 1730 km2) that were sufficiently large to have local ecological importance. Among these, only one, Matchedash Bay, was found in eastern Georgian Bay. We feel that this is a severe underrepresentation of large wetlands of ecological importance in Georgian Bay and that Herdendorf’s list should be updated with information from this study. We speculate that size alone is not a sufficient criterion for determining ecological significance. The Ontario Wetland Evaluation System [12] indicates that small coastal wetlands (i.e., <2 ha) can be grouped together to form complexes if there is a biological or hydrological rationale for doing so. For Georgian Bay, many of the smaller wetlands could be grouped into complexes since they are often found close together (within 750 m). With completion of the MCWI, we are now in the position to create the complexes, once we have a better understanding of the role that small wetlands play in terms of ecosystem functions, such as providing suitable nursery and spawning habitat for the Lake Huron fishery.

Water levels in the Great Lakes, specifically in Lake Michigan-Huron, are expected to decline as a result of climate change [19, 20]. These changes will alter the distribution, areal coverage, and vegetation structure in the coastal wetlands of eastern Georgian Bay. Wetland habitat in the MCWI has been classified according to three unique habitat zones: low marsh, which is critical habitat for fish [6], high marsh which is critical habitat for marsh birds and turtles [7, 32], and upstream habitat, that plays a critical role in controlling water quality in downstream habitats (UP) [33]. Both Wei and Chow-Fraser [34] and Midwood and Chow-Fraser [35] utilized IKONOS satellite imagery to map different types of fish habitat in the LM portion of coastal wetlands. In the HM zone, Rokitnicki-Wojcik et al. [36] developed a method that can provide detailed maps for HM vegetation. These methods can now be applied to map all LM and HM habitat in the MCWI that was delineated with IKONOS imagery. This mapping should produce consistent, baseline maps of fish habitat as well as meadow habitat for majority of coastal wetlands in eastern and northern Georgian Bay. With the acquisition of new satellite imagery, changes in vegetation coverage can be monitored and linked to the observed changes in water level in Lake Michigan-Huron.

DeCatanzaro et al. [16] were able to use road density as a surrogate for human development. They found poorer water quality (i.e., increased nutrients, conductivity, and suspended solids) in coastal wetlands adjacent to quaternary watersheds that were associated with high road density. Majority of coastal wetlands in Georgian Bay occur within watersheds that have low road density except for two (Sturgeon River and Coldwater River) [16]. The highest density occurs in the Coldwater River watershed (16.1 m/ha), and this watershed also contains the largest single area of LM habitat, and this may mean that one of the largest fish spawning areas in Georgian Bay is currently threatened by human development. Eleven of the 32 watersheds that are not currently being impacted by human development have large chains of islands. These island watersheds represent ideal conservation sites since they have limited human access, except for some cottage development. In addition, despite the fact that these islands cover less than 5% of the total watershed area within the basin, they account for nearly a quarter of all coastal habitats in eastern and northern Georgian Bay. A first step towards conserving critical wetland habitat in Georgian Bay should be to protect these islands, and, currently, the Georgian Bay Land Trust (GBLT) has managed to acquire and protect islands covering over 250 ha (http://www.GBLT.org/).

In creating the MCWI, we have provided a consistent and accurate inventory of coastal wetlands in eastern and northern Georgian Bay under low water level conditions. This project took three years to complete with help of many GIS technicians. At all times, we tried to ensure that each technician was delineating at a consistent level of accuracy. We do acknowledge that small differences in the date of image acquisition, which created slight differences in ground feature colour, may result in some discrepancy in wetland delineation. This type of error in image interpretation is unavoidable, and we believe the resulting error does not significantly alter the accuracy of the MCWI. The incorporation of the OBM wetland layer into the MCWI was necessary in order to complete the project in a timely fashion. While we found that this layer provided a sufficient level of coverage for the UP portion of our inventory, manual delineation designed specifically to identify upland habitat which was connected via surface water may provide a more accurate inventory.

While we believe that the MCWI in its current form provides a useful and comprehensive tool that should be adopted and utilized by conservation managers, we know that it can be improved with further enrichment. First, the image gap in a small portion of the French River Delta needs to be filled and all wetlands in this area need to be delineated. Secondly, wetlands identified by the inventory need to be grouped into ecologically relevant complexes in accordance with the complexing rules outlined in the OWES [12] or with suitable modifications. Finally, we recommend that satellite imagery be acquired every five years for a statistically valid subset of the MCWI. This will allow researchers and managers to track general trends in areal wetland coverage change as water levels fluctuate.

Here, we show how to create a complex habitat-based inventory of coastal wetlands for a large expanse of eastern and northern Georgian Bay, Ontario, Canada. The application of this project beyond the scope of the Laurentian Great Lakes is widespread as managers are continually in need of cost-effective methods to produce high-quality and ecologically relevant geospatial data. Although aerial photography is considered the gold standard for habitat identification, in the context of the Great Lakes, the MCWI is able to provide a static view of habitat conditions across a significant portion of the entire basin, which has proven to be too costly in past projects (GLCWCI). We recommend that managers undertaking mapping projects at similar spatial scales as the MCWI and GLCWCI consider the benefits of having contiguous data coverage within a time scale where geographic comparisons are valid and regional differences due to succession are minimized.

Acknowledgments

The authors thank Georgian Bay Forever for providing a research grant to P. Chow-Fraser to undertake this three-year project and for giving them access to the IKONOS imagery. Many GIS technicians contributed to the final inventory, but Kristina Cimaroli undoubtedly made the most important contributions. They acknowledge additional funding in the form of an NSERC scholarship to D. Rokitnicki and an OGS scholarship to J. Midwood.