Abstract

We used aerial photography, field measurements, and bird surveys to evaluate 7 Ohio mitigation wetlands for their capacity to support avian guilds at both local and landscape scales. At the local scale, we assessed each wetland with habitat suitability indices (HSI) for eight wetland-dependent bird species as indicators for four guilds: wading, diving, dabbling, and emergent dependent. We characterized landscapes within 2.5 km of each wetland by measuring the buffer width, road density, connectedness, and anthropogenic land development. The changes in landscape variables over time were determined by comparison of aerial photos taken near the time of wetland construction and near the time of this study. Bird abundance data were poorly correlated with HSI scores but were well described with logistic models of buffer width, wetland area, and road density. Our results suggest that landscape variables are better predictors of bird abundance than HSI scores for these guilds in these wetlands.

1. Introduction

Habitat for wetland-dependent species in the United States has declined in quality and quantity over a long history of wetland drainage, filling, and impairment [1, 2]. Recently, national policies of mandatory mitigation have slowed the rate of wetland loss [2], and for the period of 1998–2004 Dahl [3] estimated a net gain in total wetland area for the first time. Through wetland creation and restoration the United States has thus achieved and exceeded the goal of no-net-loss, but there continue to be questions about the ecological effectiveness of mitigation wetlands [4]. In particular, the capacity of created and restored wetlands to support diverse biota is inconsistent. Many variables can influence the establishment of a diverse community or certain desirable species, and it is now clear that biotic restoration is dependent not only on characteristics of the constructed ecosystem but also on its placement within the landscape [5, 6]. The relative importance of local and landscape factors remains a central question in restoration ecology.

Created and restored wetlands are typically designed for the establishment of appropriate biotic communities and abiotic conditions within the zone of inundation. These local attributes are planned according to prescribed performance criteria (e.g., diversity of the native plant community, basin depth, or hydrologic regime) and for the ecosystem services they provide (e.g., flood retention capacity, habitat provision). For example, local attributes are important in the restoration or creation of wetland-dependent bird habitat. Fairbairn and Dinsmore [7] found coverage of emergent vegetation, total wetland area, and perimeter-to-area ratio to be predictive of bird community diversity in prairie pothole wetlands. Other studies [810] have shown correlation of local wetland characteristics with the occurrence of particular bird species, which are used as indicators of habitat quality. Habitat suitability models have been used in species-based ecosystem evaluation and management at the local scale [11, 12].

Of course, local attributes occur within a landscape matrix, and the landscape is a critical factor in wetland restoration or creation for avian habitat. Numerous studies provide evidence that landscape attributes like road density [13], proximity to other ecosystems [7, 14], and anthropogenic development [15] affect local bird diversity. In general wetland bird diversity has been shown to decrease in response to increased road density, habitat fragmentation and isolation, and anthropogenic land use in the landscape [13, 16, 17]. And, as with local attributes, landscape factors have been associated with the presence of particular bird species [18, 19].

The establishment of a suitable habitat for indicator bird species or for general avian diversity are noble goals of wetland creation and restoration. However, analysis of habitat suitability for bird guilds may provide a better assessment of the wetland’s capacity to support a functional community [20, 21]. Guilds offer certain advantages over the use of single species or general diversity in ecosystem assessment, including the ability to analyze the viability of functional niches [21]. A variety of bird guilds have been used in the analysis of wetland habitat quality including those based on specialist-to-generalist gradient [22], response to landscape disturbance [23], migratory status [24], habitat preference [25], and feeding preference [26]. Verner [27] proposed the use of bird guilds defined according to zones of habitat as an indicator of an ecosystem’s capacity to support wildlife. This approach may be particularly useful in the evaluation of mitigation wetlands, as it would allow for a direct assessment of the efficacy of wetland design.

In this study, we use physical aspects of restored and created wetlands as predictors of structural suitability for eight bird species in four habitat-zone guilds: wading, diving, dabbling, and emergent dependent. We assess the extent to which the construction of suitable local habitat characteristics corresponds with the actual abundance of birds from the respective guilds. We also analyze the relative effectiveness of structural wetland attributes and landscape factors in the prediction of avian guild presence in restored and constructed wetlands.

2. Methods

We studied local attributes, landscape attributes, and bird abundance in seven active mitigation bank wetlands in Ohio, USA (Figure 1). In 2006 these wetlands ranged from one to thirteen years since construction ( 7 . 8 ± 1 . 5 years, m e a n ± S E ), with the oldest constructed in 1993 and the youngest in 2005. The wetland mitigation bank objectives and site plans were obtained from the Army Corps Districts of Huntington and Buffalo. To facilitate spatial analyses, we obtained aerial photographs of each wetland and its surrounding landscape. Each wetland was visited once between early June and late September of 2006 for verification of aerial photograph analyses and for measurement of local attributes. Wetlands were visited on two additional occasions in 2006 and 2007 for avian survey.

Analysis of habitat at the local scale included a bird habitat field survey. Eight waterbird species were chosen as proxies to represent four guilds: two wading species, the great blue heron (Ardea herodias), and great egret (Ardea alba); two diving species, the lesser scaup (Aythya affinis) and American coot (Fulica americana); two dabbling ducks, the blue-winged teal (Anas discors) and wood duck (Aix sponsa); and two emergent-dependent bird species, the marsh wren (Cistothorus palustris) and red-winged blackbird (Agelaius phoeniceus). These species were selected according to their primary feeding and nesting guild, their geographic range, and the availability of a US Fish and Wildlife Service Habitat Suitability Index (HSI) model for each (corresponding references given Table 6). The HSI scores are based on a zero-to-one scale in which a score of one indicates ideal habitat and zero indicates an absence of suitable habitat. In order to complete the Habitat Suitability Indices for each guild, several field measurements were required including the cover, growth form, and height of dominant vegetation, average water depth, degree of vegetation interspersion, average water depth under emergent herbaceous and woody vegetation, emergent plant identification, as well as buffer width and the capacity for foraging and nesting in the regional proximity of the wetland (Table 6). HSI scores were averaged for each guild as an estimate of guild suitability in each wetland.

A United States Geological Survey (USGS) Digital Orthophoto Quadrangle (DOQ) from a time prior to or shortly after the construction of each mitigation bank (1988–2001) was obtained for each site through the online database TerraServer. We digitized these “early” landscape photographs to a radius 2.5 km from the center of each wetland complex using Global Mapper v. 7 (Global Mapper Software, LLC). Land uses were organized into four categories: agriculture, residential/urban, natural area (wooded, grassland, wetland), and open water. “Recent” (January 2004 to July 2006) aerial photographs were obtained for each site through TerraServer or from MetroFlyer Aerial and Architectural Photography (Table 1) and georectified with Global Mapper v.7. For each wetland, we digitized the surrounding landscape within a radius of 2.5 km on a recent photograph using the same categories as the DOQ. A second recent photograph of higher resolution was used to digitize basin and buffer features. The early and recent aerial photographs used for landscape analyses differed in age by 9 . 9 ± 0 . 3  y ( m e a n ± S E ). As adjusted for wetland area, these photographs were compared to determine the rate of residential development and road density in the landscape surrounding the wetland, and the connectedness of the wetland with other natural areas.

From the aerial photograph analyses we computed five landscape variables: average buffer width, connectedness, road density, rate of residential development, and landscape development index. Buffer width was defined as the average distance from the edge of the wetland to the nearest anthropogenic land use (e.g., farm, road, or residential area) in each of the four cardinal directions. Following Findlay and Houlahan [17] we determined road density by dividing the total length of road within the 2.5 km radius by the area within the radius, including the wetland itself. We did not differentiate roads by type or traffic load. We measured connectedness as the ratio of area in a conterminous polygon of natural area (e.g., forest, wetland, grassland, and open water) around and including the wetland to total area in a 2.5 km radius from the center of the wetland after Whited et al. [13]. Using the early and recent aerial photographs the change in connectedness and road density per year were determined by dividing the early to recent difference by the number of years elapsed. We used the Landscape Development index (LDI) of Brown and Vivas [37] as adapted to Ohio by Mack [28] to evaluate the landscape in a 1 km radius from each wetland. The LDI is a weighted average of coefficients that rank various land use classes from one (natural system) to ten (highly modified urban environment) to quantify the anthropogenic stress in the wetland landscape.

Each wetland was surveyed for birds on two occasions, once in summer 2006 and once in summer 2007. Birds were surveyed with the round count method [38], in which the wetland is circumnavigated by the surveyor. Each survey took place on a single day, beginning one half hour before sunrise and continuing until circumnavigation was complete; survey effort was thus related to wetland area. We categorized bird abundance according to habitat preferences for nesting and foraging according to Ehrlich et al. [39]. Wetland-dependent birds were classified by dependence on mud flats and shallow emergent zones (wading birds), deep open water (diving birds), shallow open water (dabbling birds), and emergent vegetation (emergent-dependent birds). Bird species without a clear wetland association were not considered in this study.

We tested mean HSI statistics for similarity among the eight waterbird species with one-way analysis of variance (ANOVA; SPSS 15.0). Significant differences among mean HSI scores were determined with Tukey’s pairwise comparisons. The significance of difference between the mean area of natural, agricultural, open water, and residential/urban land in early and recent landscape surrounding each wetland was determined with paired t-tests. We also compared the early and recent buffer width, connectedness, road density, and LDI with paired t-tests. Pearson’s r was used to describe correlations among bird survey data, HSI variables, landscape variables, and wetland age. We used curve fitting analyses to determine appropriate data transformation for linear regression. Significant predictors were transformed according to the logistic equation ( 𝐾 1 + 𝛽 0 𝛽 1 𝑥 ) 1 , where 𝐾 is the curve’s upper boundary value, 𝛽 0 and 𝛽 1 are regression coefficients, and x is the predictor value. Significance for all tests was evaluated at 𝛼 = 0 . 0 5 ; marginal significance is noted for correlation and regression models at 𝛼 < 0 . 1 0 .

3. Results

The wetlands considered in this study encompass 536 ha, with a mean (±SE) area of 7 6 . 6 ± 2 5 . 4  ha. Upland buffer between the wetland and the nearest anthropogenic land use averaged 2 5 5 ± 9 2  m. Based on our analyses of early aerial photographs, the wetlands were constructed in landscapes (2.5 km radius) that are primarily agriculture ( 6 0 ± 4 % ) and natural ( 2 4 ± 3 % ; predominantly wooded) land use, with 6 ± 1 % of the surrounding landscape in open water and 8 ± 2 % in residential or urban development. The land area surrounding the mitigation bank wetlands in natural use ( 𝑃 = 0 . 4 6 ), agriculture ( 𝑃 = 0 . 4 4 ), and open water ( 𝑃 = 0 . 2 7 ) were statistically unchanged from early to recent photograph while land area in residential or urban development increased by an average of 2 1 ± 8  ha, from 7 . 5 ± 2 % in the early photograph to 8 . 6 ± 2 % in the recent photograph ( 𝑃 = 0 . 0 7 ).

Mean Habitat Suitability Index values differ significantly among waterbird species ( 𝐹 7 , 5 5 = 4 . 3 ; 𝑃 < 0 . 0 0 1 ) and among wetlands ( 𝐹 6 , 5 5 = 2 . 9 ; 𝑃 = 0 . 0 2 ; Table 2). Across all wetlands, habitat suitability was higher for wading ( 0 . 4 0 ± 0 . 1 2 ) and diving ( 0 . 4 4 ± 0 . 0 6 ) guilds than for dabbling ( 0 . 2 3 ± 0 . 0 3 ) and emergent dependent ( 0 . 2 9 ± 0 . 0 8 ), though the difference was not statistically significant ( 𝐹 3 , 5 5 = 1 . 5 ; 𝑃 = 0 . 2 3 ). Connectedness to natural areas in the landscapes surrounding the wetlands in this study ranged from moderate connectedness to relative isolation (Table 3). The average connectedness at the time of this study was 1 2 . 5 ± 3 . 9 % , an increase from the early photographs ( 1 1 . 4 ± 3 . 8 % ; 𝑃 = 0 . 0 3 ). Road density at the time of this study ( 1 1 ± 0 . 8  m/ha) had not changed significantly since wetland construction ( 𝑃 = 0 . 3 ). The slight increase in residential development from early to recent photographs corresponds with an average annual increase of 2 . 0 ± 1 . 1  ha/yr in residential or urban area in the landscape within a 2.5 km radius of these wetlands. The Land Use Development Index (determined for a 1 km radius from each wetland) ranged from 1.9 to 3.8 for these landscapes, with a mean of 3 . 1 ± 0 . 3 .

Thirty-five wetland-dependent bird species were represented by 1437 individuals in the seven wetlands surveyed (Table 7). Emergent-dependent birds were the most abundant (618), followed by dabbling (589), diving (157), and wading birds (73). Bird abundance by wetland ranged widely, from 444 at Slate Run to 17 at Castalia Quarry. Abundances of the eight representative species in this study were poorly correlated with the corresponding HSI at these seven wetlands. With the exception of blue-winged teal ( 𝑟 = 0 . 7 ), no species HSI correlation value exceeded 0.35 and no guild suitability score exceeded 0.44. Landscape variables were more strongly correlated with bird abundances than HSI values (Table 4). For all four guilds, buffer width was positively correlated with abundance (Table 4). For dabbling and wading birds, total wetland area was positively correlated with abundance while road density and LDI were negatively correlated (Table 4). Wetland age was not significantly correlated with the abundance of birds in any category. The relationships of predictor variables with guild abundance were nonlinear and best described with logistic models. Regression model parameters for significant predictors are given in Table 5. Buffer width, wetland area, and road density were the only significant predictors of bird abundance for any guild.

4. Discussion

In this study we attempt to discern attributes of mitigation wetlands that are useful predictors of bird abundance. At the local scale, characteristics of created and restored wetlands such as water depth, hydrologic regime, substrate, and vegetation form, cover, and interspersion are ostensibly all design criteria. These characteristics are also components of the Habitat Suitability Indices for bird species in this study (Table 6). We test the use of HSI scores as predictors of bird abundance and find them to be ineffective. Though all eight species occurred in at least one wetland in the study, only the abundance of the blue-winged teal was predicted reasonably well by HSI scores ( 𝑟 = 0 . 7 , 𝑃 = 0 . 0 4 ). The blue-winged teal HSI, more than any other HSI used in this study, is strongly predicated on landscape attributes such as regional wetland density, proximity, and area. The generally poor correlation of the other HSI scores with bird abundance suggests that the design of created and restored wetlands to local HSI specifications is alone insufficient to attract and support particular bird species. This may be a result of imperfect indices [4043] or confounding factors of landscape, stress, and disturbance. Similarly, our proxy for guild suitability was an unsatisfactory predictor of actual guild abundance. Other researchers [23, 44] have provided evidence of general response in bird guilds to changing environmental conditions. It may be that habitat-zone guild suitability is an application for which HSI is not well suited, or it may be that the suitability of local conditions for particular guilds were overwhelmed by landscape factors in these wetlands.

While the predictive power of our study is limited by sample size, the data suggest that landscape-level variables are better predictors of avian abundance than HSI scores. The best predictors of bird abundance in this study were wetland area and buffer width. Larger wetlands support greater numbers of dabbling and wading birds—an effect that might be expected according to the theory of island biogeography [45]. Other authors have found the total wetland area to be a significant predictor of bird species richness and density [7, 17]. Our regression results indicate that some Ohio mitigation bank wetlands are large enough to support substantial populations of dabbling and wading birds, though area appears to be a limiting factor for wetlands smaller than 25 ha. The influence of buffer width has routinely been cited for its role in nutrient dynamics in wetlands [46] and its importance for herptile populations [47, 48], though its effects on bird populations are not as well studied. According to Semlitsch and Jensen [49], numerous studies have recommended a terrestrial buffer width of 30 to 60 m as a means of protecting the wetland from landscape stressors. By this standard, the wetlands in our study are well buffered, with a mean width of over 250 m. Regression models in this study show a potential beneficial effect of buffer width up to 1000 m for diving and dabbling birds and up to 2000 m for emergent-dependent birds.

In the wetlands we surveyed, road density and LDI had a weak negative relationship with bird abundance. Findlay and Houlahan [17] indicated that road density is negatively correlated with wetland bird species richness in Ontario. While they find the greatest road density effect at a radius of 500 m, the effect is significant to 2000 m and corresponds with a 14% decline in wetland bird species richness for each 2 m/ha increase in paved road surface. Our regression model for wading birds predicts a 40% loss of abundance for each 2 m/ha increase in road density. The average road density for wetlands in our study, 1 3 ± 1  m/ha, is well above the level at which Findlay and Houlahan [17] noted this effect, suggesting that road density is likely a stress on these wetlands. Anthropogenic stressors are embedded in Mack’s [28] LDI for Ohio. Mack [28] found a strong positive correlation of the LDI with the Vegetation Index of Biotic Integrity (VIBI; [50]) for Ohio wetlands, indicating that landscape-level anthropogenic stress is associated with the wetland plant community. For our indicators of the avian community, LDI was a significant predictor for the abundance of only the wading guild.

Our results do not show a strong relationship of bird abundance with the degree of connectedness with adjacent natural areas. The findings of Whited et al. [13] imply that the positive correlation between natural area connectedness and waterbird species richness is strongest at the 2.5 km interval. At this spatial scale, increases in connectedness were positively associated with the number of waterbird species present in the wetland. In theory, connectedness of natural areas decreases the impact of surrounding land use stresses, such as non-point source pollution or invasive species, by minimizing the length of the perimeter in proportion to the core habitat and reducing the edge effect [13]. Connectedness also increases habitat patch diversity and allows for dispersion among local, native populations within the connected natural area [13]. It may be that the species in our study are not particularly sensitive to connectedness, or the influence of connectedness on bird abundance may be masked by wetland area, buffer width, and landscape stressors.

Given these results on habitat characteristics, landscape features, and bird abundance, we can comment on the capacity of the wetlands in our study to support wetland-dependent bird species. Bird abundance is associated with wetland area, buffer width, and landscape stress to a greater degree than it is associated with local attributes. Some restored or created wetlands, including some in this study, are large, well-buffered ecosystems that have been placed in relatively low-stress environments. Others, however, provide avian habitat that may be limited by the wetland’s configuration or placement within the landscape. Local wetland characteristics that are critical for breeding and survivorship—such as vegetation cover, interspersion, foraging habitat, and nesting sites—are likely to mature and improve as the restored and created wetlands age, but anthropogenic encroachment in the landscape, even at the relatively slow rate of development we note in this study, is likely to become more stressful over time. For this reason, construction of large, well-buffered wetlands in carefully selected locations is a critical consideration in the effort to maximize the functional replacement of avian habitat in restored and created wetlands.

Acknowledgments

We thank Melissa Kennedy, Jenna Lawson, Alexander Mora, Chelsea Mora, and Jessica Piispanen for their field assistance. This research was funded by the Anderson Scholarship Endowment of Denison University, the Denison University Research Foundation, and the Ohio Biological Survey.