Abstract

Some avian species in the southeastern United States are declining, and population decreases may arise from changes in vegetation type area or structural condition. Our objective was to compare abundance of conservation priority bird species with landscape variables. We found, even in the highly forested Coastal Plain of Georgia, that areal extent and core area of cover types were related to abundance for certain bird species. Acadian flycatcher and field sparrow had models that incorporated positive area variables. Downy woodpecker, northern parula, orchard oriole, prairie warbler, and summer tanager had models that included positive area and edge associations with varying scales and vegetation types. Edge appeared in models for red-bellied woodpecker, blue jay, and brown-headed cowbird. More than half of all species did not have models that met prediction thresholds. Systematic assessment of area requirements for declining species provides information for management, conservation, and research.

1. Introduction

Populations of certain bird species in the southeastern United States are declining, particularly disturbance-dependent species associated with grasslands, shrublands, and open forests [1]. Population trends may arise from land use changes in vegetation type area or structural condition. Stand elements, including vegetation composition and structure, can affect a variety of bird species, such as cavity nesters that require older trees or early successional species that need an open canopy and midstory. Thus, availability of stand-level structural elements may also be related, cumulatively at the landscape scale, to abundance of some declining bird species.

Most bird-landscape studies have taken place outside of the Southeast. Compared to regions with forests fragmented by agriculture and urbanization, Coastal Plain landscape research on breeding birds has been equivocal, perhaps because many patches of one forest type or stage are enclosed within forest of another type or stage [2, 3]. Such studies include Krementz and Christie [4], who detected no effect of clearcut size on species richness or juvenile to adult ratios in birds captured in mist nets. Aquilani and Brewer [5] determined that wood thrush (see Table 2 for scientific names) nest success was greatest in large fragments of vegetation types and least near clearcut edges, primarily due to varying predation levels. Edge increased nest predation and negatively affected indigo bunting nesting success [6], but edge did not depress Acadian flycatcher nest survival [7]. In investigations of bottomland hardwood widths, Hodges and Krementz [8] and Kilgo et al. [9] found that species richness increased with riparian width.

Given that there is incomplete regional knowledge about avian habitat requirements, exploratory models can contribute valuable information about landscape metrics associated with avian presence for conservation management and research. Habitat selection involves multiple scales, or at least changes depending on observation scale, and may vary by region [10, 11]. Regional habitat modeling at different scales for birds that are declining may help establish area sensitivity classifications. Our objective was to identify land class variables at varying scales that potentially predict abundance of priority avian species in Coastal Plain Georgia.

2. Methods

2.1. Study Area

Southeastern Georgia is in the Coastal Plain, a low, flat physiographic region of the southeastern United States. Coastal Plain vegetation consists of upland forests interspersed with wetlands and poorly drained flatwoods. Land use consists of row crop agriculture and intensively managed pine forest, with urbanization along the coast [12]. In 1999, Georgia had 9.5 million ha of timberlands, 60–65% of the State’s total area, including 2.5 million ha in planted pine and 1.85 million ha in natural pine [13]. Forested areas are young; 3.3 million ha are in the seedling sapling stage (less than 12.7 cm dbh; [13]). Disturbance agents include fire, hurricanes, tornadoes, floods, and ice storms.

2.2. Data Sets

We combined the North American Breeding Bird Survey (BBS; [14]), coordinated by USGS Patuxent Wildlife Research Center, with the Georgia Gap Analysis land cover grid (GA-GAP; Kramer et al. 2003 [12]), to correlate bird species abundance with spatial metrics. Breeding Bird Survey routes are approximately 40 km long and consist of 50 points that are 0.8 km apart. Volunteers record birds within a 400 m radius during 3 minutes at each point. One bias of BBS is that surveys occur alongside roads. However, there are roads throughout Georgia, where roadless areas may be limited [15]. The GA-GAP classified 30-meter resolution Landsat Thematic Mapper satellite imagery, using 1996–1998 imagery. There are 44 land cover classes with an overall accuracy of 75.5%. The satellite imagery occurred before the accuracy assessment, which contributed to error rate.

From the BBS database, we selected all routes in Georgia’s Coastal Plain and Flatwoods with 3 survey years under approved conditions during 1995 to 1999 ( 𝑛 = 2 7 routes). We divided the routes into 5, 10-stop partial routes about 8 km in length and selected the straightest (i.e., covering greatest area) partial route from each end (i.e., not the middle partial route). We placed buffered extents of 0.5, 1, 2, and 4 km around each partial route. We eliminated 2 partial routes, due to overlap at the 4 km buffer extent, leaving 52 route segments. For all these operations, we used ArcGIS (ESRI, version 9.0, Redlands, CA, USA).

We retained, with some reclassification, 8 GA-GAP land cover classes for analysis: (1) hardwood forests (2) hydric hardwoods (bottomland hardwood, cypress-gum swamp, evergreen forested wetland), (3) clearcut (recent clearcuts, sparse vegetation, and other early successional areas, (4) pasture/hay, (5) mixed forest, (6) managed pine (loblolly shortleaf, loblolly slash), (7) longleaf pine, and (8) shrub (sandhill, shrub wetland). We clipped the reclassified grid using the buffered partial route shapes, creating grids of each partial route at 4 buffer distances.

We used FRAGSTATS [16] to compute 7 spatial metrics for each class type. Metrics for modeling included area (AREA; mean patch area, depends on patch size and number), core area (CORE; mean core area of patch, excludes 90 m buffer from edge), cohesion (COHESION; connectivity of class type), edge density (ED; edge length of patch standardized by area), and the interspersion and juxtaposition index (IJI; class type intermixing). Additionally, for landscape descriptive statistics (Table 3), we calculated percentage of landscape (PLAND; proportional abundance of class type, standardized by area) and core percentage landscape (CPLAND; proportional abundance of class type core area, excludes 90 m buffer from edge) for each cover class.

We chose bird species if they scored as regionally important by Partners in Flight for the southeastern Coastal Plain region [17] (Table 2). We also included brown-headed cowbird, a nest parasite, and blue jay, a nest predator, because of their possible impact on declining species. Then, we averaged BBS counts for each species by year (i.e., mean of 3 years) and partial route, to calculate a species mean; we did not keep species with means below 0.20. Routes were categorized as low abundance for less than the mean and higher abundance for greater than or equal to the mean.

2.3. Statistical Analyses

We randomly used 37 partial routes for modeling, while reserving 15 partial routes for validation. Although there was little correlation, we removed one variable for each pair that was at least 70% correlated (PROC CORR; SAS software, version 9.1, Cary, North Carolina, USA) based on the following order to retain: AREA, CORE, ED, COHESION, and IJI due to our interest in detecting species that may be vulnerable to areal loss. Then, for each species and extent, we selected the 5 best fitting, one to 4 variable models, using logistic regression with score selection (PROC LOGISTIC). We evaluated these candidate models with Akaike’s Information Criteria corrected for small sample size (AICc).

To assess model accuracy, we used all models within 2 AICc units of the least AICc value to predict lesser or greater abundance for 15 model validation routes (PROC LOGISTIC). We classified model fit as correct for a route if predicted probability was greater than or equal to 50% and bird abundance mean fell within the greater abundance category, or alternatively if probability was less than 50% and bird abundance mean was within the lesser abundance category. Final best model selection incorporated models with the greatest prediction rate from models that correctly predicted at least 11 out of 15 routes at each buffer extent. We removed models with more variables if there was a nested smaller model that predicted equally well at the same extent. Lastly, we used the SAS 𝑐 statistic, which measures model accuracy by the area under the curve (AUC) for a Receiver Operating Characteristic (ROC) curve at the 0.5 threshold, to determine how well the models fit all 52 routes, and we eliminated any models with a c statistic below 0.80.

To determine if model fit would increase by adjusting for spatial variability, we compared the final models for each species using GLIMMIX with and without a spatial covariance structure. For GLIMMIX covariance parameter estimates, we used estimates from a variogram for each model, using the residuals from PROC LOGISTIC.

3. Results

Five species had models that contained a mixture of positively associated area and edge variables (Table 1). Downy woodpecker greater abundance was associated positively with clearcut area and pasture edge density. Northern parula models combined positive variables of hydric hardwood edge density with mixed forest area at 0.5, 1, and 2 km extents along with clearcut area at 1 and 2 km. Orchard oriole models consisted of positive variables that included pasture edge density at 1 and 4 km, hydric hardwood edge density and area, and core area of longleaf pine. Prairie warbler models encompassed positive variables of hydric hardwood area and shrub edge. Summer tanager was linked to positive model variables of pasture edge density at 1, 2, and 4 km extents, longleaf pine area and edge density, and managed pine area.

Five species had models that did not contain a mixture of positive area and edge variables and involved only one model, and thus only one buffer extent (Table 1). Acadian flycatcher abundance was associated positively with mixed forest area. Field sparrow models included positive model variables of shrub core at multiple extents. Three species had models that contained edge density but not area variables. Blue jay was related positively to pasture edge density. Brown-headed cowbirds were associated positively with longleaf pine edge density. Red-bellied woodpecker model variables were related positively to longleaf pine edge density. Models for pine warbler did not have any positive area or edge variables.

Brown-headed nuthatch, brown thrasher, Carolina chickadee, Carolina wren, eastern kingbird, eastern towhee, eastern wood-pewee) indigo bunting, northern bobwhite, red-headed woodpecker, pileated woodpecker, white-eyed vireo, and wood thrush did not have any models that correctly predicted the minimum 11 out of 15 on the validation routes or did not meet the threshold 0.80 for 𝑐 statistic value.

For seven species, a spatial covariance structure did not produce better model fit. Prairie warbler and northern parula models improved with a spatial term, and for summer tanager, the models at 2 and 4 km improved with a spatial term. The field sparrow second and fourth models had greater fit with a spatial term.

4. Discussion

Even in the highly forested Coastal Plain core area, areal extent, edge density, and other landscape characteristics were important predictors of abundance for certain species, and the Coastal Plain may not provide enough continuous vegetation for area-sensitive species (Table 3). Despite the limited area sizes, area was still prevalent in models for most bird species. Area was a model variable for Acadian flycatcher, a species which may be area sensitive [18]. Field sparrow also appeared to respond to area, and they may avoid edges when interior dense vegetation is available [19], prefer to nest away from edge [20], or have reduced brood parasitism rates away from edge [21, 22]. Downy woodpecker, northern parula, orchard oriole, and summer tanager exhibited area as a model variable, which helps support positive area findings in previous studies [23, 24].

A mixture of area and edge variables developed in most models. One forest type adjacent to another forest type may expand the functional area of each forest, and both edge and area would become important. Species associated with dense shrubby vegetation, which is characteristic of edges, may be more likely to respond to both edge and area. Alternatively, edge and area model combinations may reflect that scale and vegetation type influence area relationships. For example, a species could have, as a model variable, edge density linked with a vegetation type that the species avoids, if the edge represented lower abundance of the vegetation type. Occurrence of area and edge of the same vegetation type may show more clearly that the species is responding directly to the vegetation type, rather than the particular metric. Orchard oriole had models that included hydric hardwood area at 1 km and hydric hardwood edge density at 4 km, whereas summer tanager models incorporated longleaf pine area at 1 and 2 km and longleaf pine edge density at 4 km. However, these are separate models and thus could represent distinctive selection at different scales.

Brown-headed cowbirds were associated with their breeding habitat of an ecotone [25], specifically a longleaf pine edge, and away from larger grassy or woody extents. Perhaps longleaf pine edges provide a balance between breeding opportunities and access to feeding areas. Likewise, blue jay was linked to pasture edge density and thus may be a problem for birds nesting in ecotones.

The models did produce some interesting bird-vegetation type associations, compared to documented associations [26], some of which may provide new insight, and most may result from vegetation type classification error in the Gap land cover grid, bias in the BBS surveys, or spurious results from modeling. Clearcut areas, although a critical vegetation type, had an accuracy rate of only about 50% for the Gap land cover. The prairie warbler model did not contain clearcuts, rather hydric hardwoods and shrub, and the orchard oriole model also had hydric hardwoods, rather than open areas; this may indicate a Gap misclassification of clearcuts, and replacement of the hydric hardwoods vegetation type by clearcuts may correct these models. Northern parula had models containing clearcut area and forest edge. If the clearcut vegetation represents hardwoods, then the model is uncomplicated. Alternatively, northern parula use gaps, particularly after breeding [27], which might account for forest edge in the model, and to some extent, clearcuts if they represent small gaps. Field sparrow is associated with shrub/scrub [28], and although field sparrow models involved shrub area, the models also contain a negative association to pasture edge density, suggesting that field sparrows are not using edge as much as they are using either core area or else an abundance of pasture.

Scale affects which variables will be represented, and therefore scale choice is important. Each buffer extent included models for 3–5 species, and 8 species only had models at one extent. For species with models at multiple extents, there was variable overlap among extents, which may reflect continuity in habitat selection. Also, in part scale differences may contribute to conflicting results in the literature about area sensitivity of some species.

Species that did not have models that met our criteria could be selecting sites based on different factors than modeled vegetation types or landscape metrics. Carolina chickadee, eastern towhee, northern bobwhite, and pine warbler may have area or edge requirements [23, 29, 30], but modeling did not reflect these metrics. Stand scale elements, such as vegetation structure and composition including residual trees, may be particularly important for these species. Species that should be associated with a specific but rare vegetation type, such as brown-headed nuthatch and longleaf pine, may be using another vegetation type with appropriate structure [31]. For example, high density pine forests often lack a developed understory and pine forests without disturbance may contain a hardwood midstory rather than a herbaceous ground layer and shrubby understory, potentially limiting abundance of some species typically associated with open-canopy pine forests [32].

These models provide basic associations of bird species with vegetation type metrics in a region where relatively little is known about landscape relationships. Edge variables proved important for two species, blue jay and brown-headed cowbird, known to affect populations of some avian species through nest predation and parasitism, respectively. Area or a combination of area and edge is important for seven species of conservation concern, despite small areal extents and little core area in the Coastal Plain of Georgia. For other declining species, stand-scale metrics may be more critical, particularly when vegetation types are fragmented. Further research by modelers, field scientists and managers can contribute to understanding habitat that sustains avian populations.

Across a broader landscape, management that increases patch sizes of vegetation types may enhance habitat suitability for a large suite of species. Conversion to agriculture and large-scale timber harvesting in the nineteenth and early twentieth century affected much of the southern landscape, including Georgia [33]. Approximately 9.5% of southeastern timberlands are publicly owned, in a patchwork of sizes and managed by various agencies [34]. Public land managers and private land owners, in collaboration, are in position to develop an appropriate mix of vegetation types and components across the landscape that optimizes biodiversity and addresses regional conservation concerns. Such a mix should include mature pine and hardwood forests, as well as early successional grasslands and shrublands. Present day forest management, with careful planning, can help coordinate development of continuous expanses of similar vegetation types while conserving low-contrast borders, to limit interior exposure to nest predators and parasites. Because there now are harvest size regulations, clearcut areas averaged 92 ha (20–600 ha) on commercial lands and 36 ha in public forests [35] in the past and more recent cut sizes ranged from about 16 to 32 ha [36]. The Sustainable Forestry Initiative requires the average size of cuts to not exceed 49 ha [37]. Nevertheless, large areas of similar stages can remain close through interspersion of harvests that are slightly offset in time; larger harvests and clustered cuts will retain more extensive areas of similar vegetation age, from regeneration to harvest stages, than multiple dispersed clearings.

Acknowledgments

This paper was funded by The National Council for Air and Stream Improvement, Inc., Federal Aid in Wildlife Restoration (Project W-48, Study 57), Mississippi Department of Wildlife, Fisheries and Parks, Weyerhaeuser Company, International Paper Company, MeadWestvaco Corporation, and Boise Cascade Corporation. Mississippi State University Department of Wildlife and Fisheries provided facilities. This is publication XXX of the Forest and Wildlife Research Center, Mississippi State University.