Estimating Population Size of Elusive Animals with DNA from Hunter-Collected Feces: Four Methods for Brown Bears

Abstract:  Noninvasive genetic methods can be used to estimate animal abundances and offer several advantages over conventional methods. Few attempts have been made, however, to evaluate the accuracy and precision of the estimates. We compared four methods of estimating population size based on fecal sampling. Two methods used rarefaction indices and two were based on capture-mark-recapture (CMR) estimators, one combining genetic and field data. Volunteer hunters and others collected 1904 fecal samples over 2 consecutive years in a large area containing a well-studied population of brown bears ( Ursus arctos ). On our 49,000-km² study area in south-central Sweden, population size estimates ranged from 378 to 572 bears in 2001 and 273 to 433 bears in 2002, depending on the method of estimation used. The estimates from the best model in the program MARK appeared to be the most accurate, based on the minimum population size estimate from radio-marked bears in a subsection of our sampling area. In addition, MARK models included heterogeneity and temporal variation in detection probabilities, which appeared to be present in our samples. All methods, though, incorrectly suggested a biased sex ratio, probably because of sex differences in detection probabilities and low overall detection probabilities. The population size of elusive animals can be estimated reliably over large areas with noninvasive genetic methods, but we stress the importance of an adequate and well-distributed sampling effort. In cases of biased sampling, calibration with independent estimates may be necessary. We recommend that this noninvasive genetic approach, using the MARK models, be used in the future in areas where sufficient numbers of volunteers can be mobilized.     

Effectiveness of Corridors Relative to Enlargement of Habitat Patches

Abstract:  The establishment of biological corridors between two otherwise isolated habitat patches is a common yet contentious strategy for conserving populations in fragmented landscapes. We compared the effectiveness of corridors with the effectiveness of an alternate conservation strategy, the enlargement of existing habitat patches. We used a spatially explicit population model that simulated population size in two kinds of patches. One patch had a corridor that connected it to a larger "source" patch and the other patch was unconnected and enlarged at the periphery by an area the same size as the corridor. Patch isolation, corridor width, patch size, and the probability that individuals would cross the border from habitat to matrix were varied independently. In general, population size was greater in enlarged patches than in connected patches when patches were relatively large and isolated. Corridor width and the probability of crossing the border from habitat to matrix did not affect the relative benefit of corridors versus patch enlargement. Although biological corridors may mitigate potential effects of inbreeding depression at long time scales, our results suggest that they are not always the best method of conserving fragmented populations.     

Evaluating the Effectiveness of Corridors: a Genetic Approach

Abstract: The effectiveness of corridors in maintaining dispersal in fragmented landscapes is a question of considerable conservation and ecological importance. We tested the efficacy of corridors as residual landscape structures in maintaining population structure in the red-backed vole ( Clethrionomys gapperi ), a closed-canopy specialist, and the deer mouse (   Peromyscus maniculatus ), a habitat generalist. In coniferous forests managed for timber production in northeastern Washington, we sampled pairs of populations in three landscape classes: (1) contiguous landscapes, in which sites were located completely within a matrix of closed-canopy forest; (2) corridor landscapes, in which sites were connected by a corridor of closed-canopy forest; and (3) isolated landscapes, in which sites were separated from one another by clearcut or young regeneration stands. For each species, we used four microsatellite loci to quantify genetic distance between population pairs. Nei's genetic distance (   D _s  ) increased from smallest to largest in the order of contiguous, corridor, and isolated landscapes for C. gapperi. For P. maniculatus, genetic distances across landscape configurations were not significantly different. The differences between the two species indicate that they respond differently to the presence of forest corridors. In managed forests, corridors between unlogged habitats appear to maintain higher population connectivity for C. gapperi than landscapes without corridors.     

Estimating Population Size with Noninvasive Capture-Mark-Recapture Data

Abstract:  Estimating population size of elusive and rare species is challenging. The difficulties in catching such species has triggered the use of samples collected noninvasively, such as feces or hair, from which genetic analysis yields data similar to capture-mark-recapture (CMR) data. There are, however, two differences between classical CMR and noninvasive CMR. First, capture and recapture data are gathered over multiple sampling sessions in classical CMR, whereas in noninvasive CMR they can be obtained from a single sampling session. Second, because of genotyping errors and unlike classical CMR, there is no simple relationship between (genetic) marks and individuals in noninvasive CMR. We evaluated, through simulations, the reliability of population size estimates based on noninvasive CMR. For equal sampling efforts, we compared estimates of population size N obtained from accumulation curves, a maximum likelihood, and a Bayesian estimator. For a closed population and without sampling heterogeneity, estimates obtained from noninvasive CMR were as reliable as estimates from classical CMR. The sampling structure (single or multiple session) did not alter the results, the Bayesian estimator in the case of a single sampling session presented the best compromise between low mean squared error and a 95% confidence interval encompassing the parametric value of N in most simulations. Finally, when suitable field and lab protocols were used, genotyping errors did not substantially bias population size estimates (bias < 3.5% in all simulations). The ability to reliably estimate population size from noninvasive samples taken during a single session offers a new and useful technique for the management and conservation of elusive and rare species.     

Linking noninvasive genetic sampling and traditional monitoring to aid management of a trans-border carnivore population

Noninvasive genetic sampling has been embraced by wildlife managers and ecologists, especially those charged with monitoring rare and elusive species over large areas. Challenges arise when desired population measures are not directly attainable from genetic data and when monitoring targets trans-border populations. Norwegian management authorities count individual brown bears (Ursus arctos) using noninvasive genetic sampling but express management goals in the annual number of bear reproductions (females that produce cubs), a measure that is not directly available from genetic data. We combine noninvasive genetic sampling data with information obtained from a long-term intensive monitoring study in neighboring Sweden to estimate the number of annual reproductions by females detected within Norway. Most female brown bears in Norway occur near the border with neighboring countries (Sweden, Finland, and Russia) and their potential reproduction can therefore only partially be credited to Norway. Our model includes a simulation-based method that corrects census data to account for this. We estimated that 4.3 and 5.7 reproductions can be credited to females detected with noninvasive genetic sampling in Norway in 2008 and 2009, respectively. These numbers fall substantially short of the national target (15 annual reproductions). Ignoring the potential for home ranges to extend beyond Norway's borders leads to an increase in the estimate of the number of reproductions by -30%. Our study shows that combining noninvasive genetic sampling with information obtained from traditional intensive/invasive monitoring can help answer contemporary management questions in the currency desired by managers and policy makers. Furthermore, combining methodologies and thereby accounting for space use increases the accuracy of the information on which decisions are based. It is important that the information derived from multiple approaches is applicable to the same focal population and that predictions are cross-validated. When monitoring and management are constrained to administrative units, census data should be adjusted by discounting portions of individual space utilization that extend beyond the focal jurisdiction. Our simulation-based approach for making such an adjustment may be useful in other situations where management authorities target portions of trans-border populations.     

Demographic Connectivity for Ursid Populations at Wildlife Crossing Structures in Banff National Park

Wildlife crossing structures are one solution to mitigating the fragmentation of wildlife populations caused by roads, but their effectiveness in providing connectivity has only been superficially evaluated. Hundreds of grizzly (Ursus arctos) and black bear (Ursus americanus) passages through under and overpasses have been recorded in Banff National Park, Alberta, Canada. However, the ability of crossing structures to allow individual and population‐level movements across road networks remains unknown. In April 2006, we initiated a 3‐year investigation into whether crossing structures provide demographic connectivity for grizzly and black bears in Banff National Park. We collected hair with multiple noninvasive methods to obtain genetic samples from grizzly and black bears around the Bow Valley. Our objectives were to determine the number of male and female grizzly and black bears that use crossing structures; examine spatial and temporal patterns of crossings; and estimate the proportions of grizzly and black bear populations in the Bow Valley that use crossing structures. Fifteen grizzly (7 female, 8 male) and 17 black bears (8 female, 9 male) used wildlife crossing structures. The number of individuals detected at wildlife crossing structures was highly correlated with the number of passages in space and time. Grizzly bears used open crossing structures (e.g., overpasses) more often than constricted crossings (e.g., culverts). Peak use of crossing structures for both bear species occurred in July, when high rates of foraging activity coincide with mating season. We compared the number of bears that used crossings with estimates of population abundance from a related study and determined that substantial percentages of grizzly (15.0% in 2006, 19.8% in 2008) and black bear (17.6% in 2006, 11.0% in 2008) populations used crossing structures. On the basis of our results, we concluded wildlife crossing structures provide demographic connectivity for bear populations in Banff National Park. Conectividad Demográfica para Poblaciones de Úrsidos en Estructuras para Cruce de Vida Silvestre en el Parque Nacional Banff     

Multiple data sources improve DNA-based mark-recapture population estimates of grizzly bears.

A fundamental challenge to estimating population size with mark-recapture methods is heterogeneous capture probabilities and subsequent bias of population estimates. Confronting this problem usually requires substantial sampling effort that can be difficult to achieve for some species, such as carnivores. We developed a methodology that uses two data sources to deal with heterogeneity and applied this to DNA mark-recapture data from grizzly bears (Ursus arctos). We improved population estimates by incorporating additional DNA "captures" of grizzly bears obtained by collecting hair from unbaited bear rub trees concurrently with baited, grid-based, hair snag sampling. We consider a Lincoln-Petersen estimator with hair snag captures as the initial session and rub tree captures as the recapture session and develop an estimator in program MARK that treats hair snag and rub tree samples as successive sessions. Using empirical data from a large-scale project in the greater Glacier National Park, Montana, USA, area and simulation modeling we evaluate these methods and compare the results to hair-snag-only estimates. Empirical results indicate that, compared with hair-snag-only data, the joint hair-snag-rub-tree methods produce similar but more precise estimates if capture and recapture rates are reasonably high for both methods. Simulation results suggest that estimators are potentially affected by correlation of capture probabilities between sample types in the presence of heterogeneity. Overall, closed population Huggins-Pledger estimators showed the highest precision and were most robust to sparse data, heterogeneity, and capture probability correlation among sampling types. Results also indicate that these estimators can be used when a segment of the population has zero capture probability for one of the methods. We propose that this general methodology may be useful for other species in which mark-recapture data are available from multiple sources.     

Demographic Analyses of a Hunted Black Bear Population with Access to a Refuge

Studying one of two bear species not experiencing widespread population decline, provides insight into the population responses of the six bear species that are in decline and into responses of other long-lived species for which data are difficult to collect. Black bear ( Ursus americanus ) sanctuaries were established in North Carolina (U.S.) in 1971 to protect core populations of bears and to provide dispersing bears for hunting. Population index values, derived from counts of bears visiting bait stations, were significantly greater inside the Pisgah Bear Sanctuary than outside and were greater along trails than along roads. Survivorship of bears outfitted with transmitter collars was greater for sanctuary bears alone than for sanctuary plus non-sanctuary bears. Monte Carlo analyses of Leslie matrices showed that the bear population in the sanctuary would be stable if cub survivorship, p ₀, was about O.7, and the population in the sanctuary plus the surrounding area would be stable if p ₀ was about 0.83. Estimates of litter survivorship in North Carolina indicate, however, that p ₀ can not exceed O.71. Overall, the matrix analyses indicated an ultimate population decline in the total bear population (sanctuary plus surrounding area). The population index of the bait station did not show a discernible decline. The Pisgah Bear Sanctuary provides dispersing bears for hunters and provides some protection for the resident bears. The sanctuary may not, however, provide resident bears with enough protection to maintain a viable breeding population within its boundaries. Reducing human access to bears and their habitat appears crucial, either by making large sanctuaries or by eliminating roads.     

Regional Consequences of Local Population Demography and Genetics in Relation to Habitat Management in Gentiana pneumonanthe

Abstract:  A joint demographic and population genetics stage-based model for a subdivided population was applied to Gentiana pneumonanthe , an early successional perennial herb, at a regional (metapopulation) scale. We used numerical simulations to determine the optimal frequency of habitat disturbance (sod cutting) and the intensity of gene flow among populations of G. pneumonanthe to manage both population viability and genetic diversity in this species. The simulations showed that even small populations that initially had near-equal allele frequencies could, if managed properly through sod cutting every 6 to 7 years, sustain their high genetic variation over the long run without gene flow. The more the allele frequencies in the small populations are skewed, however, the higher the probability that in the absence of gene flow, some alleles will be lost and within-population genetic variation will decrease even under proper management. This implies that although local population dynamics should be the major target for management, regional dynamics become important when habitat fragmentation and decreased population size lead to the loss of local genetic diversity. The recommended strategy to improve genetic composition of small populations is the introduction of seeds or seedlings of nonlocal origin.     

Temporal Analysis of Genetic Structure to Assess Population Dynamics of Reintroduced Swift Foxes

Reintroductions are increasingly used to reestablish species, but a paucity of long‐term postrelease monitoring has limited understanding of whether and when viable populations subsequently persist. We conducted temporal genetic analyses of reintroduced populations of swift foxes (Vulpes velox) in Canada (Alberta and Saskatchewan) and the United States (Montana). We used samples collected 4 years apart, 17 years from the initiation of the reintroduction, and 3 years after the conclusion of releases. To assess program success, we genotyped 304 hair samples, subsampled from the known range in 2000 and 2001, and 2005 and 2006, at 7 microsatellite loci. We compared diversity, effective population size, and genetic connectivity over time in each population. Diversity remained stable over time and there was evidence of increasing effective population size. We determined population structure in both periods after correcting for differences in sample sizes. The geographic distribution of these populations roughly corresponded with the original release locations, which suggests the release sites had residual effects on the population structure. However, given that both reintroduction sites had similar source populations, habitat fragmentation, due to cropland, may be associated with the population structure we found. Although our results indicate growing, stable populations, future connectivity analyses are warranted to ensure both populations are not subject to negative small‐population effects. Our results demonstrate the importance of multiple sampling years to fully capture population dynamics of reintroduced populations. Análisis Temporal de la Estructura Genética para Evaluar la Dinámica Poblacional de Zorros (Vulpes velox) Reintroducidos     

Evaluating multispecies landscape connectivity in a threatened tropical mammal community

Habitat corridors are important tools for maintaining connectivity in increasingly fragmented landscapes, but generally they have been considered in single‐species approaches. Corridors intended to facilitate the movement of multiple species could increase persistence of entire communities, but at the likely cost of being less efficient for any given species than a corridor intended specifically for that species. There have been few tests of the trade‐offs between single‐ and multispecies corridor approaches. We assessed single‐species and multispecies habitat corridors for 5 threatened mammal species in tropical forests of Borneo. We generated maps of the cost of movement across the landscape for each species based on the species’ local abundance as estimated through hierarchical modeling of camera‐trap data with biophysical and anthropogenic covariates. Elevation influenced local abundance of banded civets (Hemigalus derbyanus) and sun bears (Helarctos malayanus). Increased road density was associated with lower local abundance of Sunda clouded leopards (Neofelis diardi) and higher local abundance of sambar deer (Rusa unicolor). Pig‐tailed macaque (Macaca nemestrina) local abundance was lower in recently logged areas. An all‐species‐combined connectivity scenario with least‐cost paths and 1 km buffers generated total movement costs that were 27% and 23% higher for banded civets and clouded leopards, respectively, than the connectivity scenarios for those species individually. A carnivore multispecies connectivity scenario, however, increased movement cost by 2% for banded civets and clouded leopards. Likewise, an herbivore multispecies scenario provided more effective connectivity than the all‐species‐combined scenario for sambar and macaques. We suggest that multispecies habitat connectivity plans be tailored to groups of ecologically similar, disturbance‐sensitive species to maximize their effectiveness. Evaluación de la Conectividad de Terrenos Multiespecie en una Comunidad Tropical de Mamíferos     

Demographic Processes Underlying Population Growth and Decline in Salamandra salamandra

Abstract:  Human activity commonly has negative impacts on wildlife. Often, however, only a single element of the life cycle is affected, and it is unclear whether such effects translate into effects on population growth. This is particularly true for research into the causes of global amphibian declines, where experimental research focuses primarily on the aquatic larval stages but theory suggests these stages have only minor importance for population growth. We used data from long-term mark-recapture studies of two natural populations of the salamander Salamandra salamandra to confirm the predictions of population models. One population remained stable (i.e., stationary) throughout the 20 years of the study whereas the other declined to local extinction. We used mark-recapture models to break down population growth rate into its two main components, recruitment and adult survival. Survival of postmetamorphic salamanders was constant over time in the stable population, whereas the declining population was characterized by a decrease in survival and constant recruitment. Population growth was most sensitive to variation in adult survival. Current amphibian research focuses on preadult stages, and researchers assume recruitment is the most important determinant of population growth. This may not be the case. A better understanding of amphibian population dynamics is possible only through the integration of experiments, theory, and data from natural populations. Our results also suggest that amphibian conservation efforts should focus on all stages of the life cycle and their associated habitats.     

A Meta‐Analytic Review of Corridor Effectiveness

Abstract: Using corridors for conservation is increasing despite a lack of consensus on their efficacy. Specifically, whether corridors increase movement of plants and animals between habitat fragments has been addressed on a case‐by‐case basis with mixed results. Because of the growing number of well‐designed experiments that have addressed this question, we conducted a meta‐analysis to determine whether corridors increase movement; whether corridor effectiveness differs among taxa; how recent changes in experimental design have influenced findings; and whether corridor effectiveness differs between manipulative and natural experiments. To conduct our meta‐analysis, we analyzed 78 experiments from 35 studies using a conservative hierarchical Bayesian model that accounts for hierarchical and sampling dependence. We found a highly significant result that corridors increase movement between habitat patches by approximately 50% compared to patches that are not connected with corridors. We found that corridors were more important for the movement of invertebrates, nonavian vertebrates, and plants than they were for birds. Recent methodological advances in corridor experiments, such as controlling for the area added by corridors, did not influence whether corridors increased movement, whereas controlling for the distance between source and connected or unconnected recipient patches decreased movement through corridors. After controlling for taxa differences and whether studies controlled for distance in experimental design, we found that natural corridors (those existing in landscapes prior to the study) showed more movement than manipulated corridors (those created and maintained for the study). Our results suggest that existing corridors increase species movement in fragmented landscapes and that efforts spent on maintaining and creating corridors are worthwhile.     

Identifying habitat patches and potential ecological corridors for remnant Asiatic black bear (Ursus thibetanus japonicus) populations in Japan

The Japanese National Biodiversity Strategy 2010 calls for the creation of ecological networks as a biodiversity conservation policy. However, there is an obvious lack of information on the spatial distribution of many species and a lack of scientific methods for examining habitat requirements to establish the need for constructing these networks for target species. This study presents a quantitative method for assessing the need for ecological networks through modeling the potential geographic distributions of species based on a case study of local populations of Asiatic black bear (Ursus thibetanus japonicus) in Fuji and Tanzawa, Japan. A total of 1541 point records of occurrences of Asiatic black bears and 11 potential predictors were analyzed in a GIS environment. After a predictive distributional map was obtained using the Maximum Entropy (MaxEnt) algorithm, a gap analysis was carried out and population size was estimated. Approximately 24% of the bear's predicted habitat area fell within a wildlife protection area, 2% within a nature reserve, and 37% within natural parks. Conservation forest comprised 54% of the total area of predicted habitat; of this, national forest comprised 2%, and private and communal forest comprised 37%. The total estimated Asiatic black bear population in this region was 242, with 179 individuals in the Fuji local population, 26 in the Tanzawa local population, and 37 in the corridor patch between the two local populations. Our study also found a potential corridor connecting the Fuji and Tanzawa local populations, as well as potential habitat corridors in the Fuji region containing subpopulations on Mt. Fuji (119 individuals) and Mt. Kenashi (53 individuals). An additional subpopulation on Mt. Ashitaka (7 individuals) is isolated and not fully protected by a zoning plan. Mt. Furo's subpopulation is considered to be almost extinct, although black bears were observed here until 2002 based on the report by Mochizuki et al. (2005). The total black bear population of the Fuji-Tanzawa region is considered to be “endangered”; thus, an adequate population size might be difficult to maintain even if this region were to be internally connected by means of an ecological network.     

Relationship between Population Size and Fitness

Abstract:  Long-term effective population size, which determines rates of inbreeding, is correlated with population fitness. Fitness, in turn, influences population persistence. I synthesized data from the literature concerning the effects of population size on population fitness in natural populations of plants to determine how large populations must be to maintain levels of fitness that will provide adequate protection against environmental perturbations that can cause extinction. Integral to this comment on what has been done and what needs to be done, sThe evidence suggests that there is a linear relationship between log population size and population fitness over the range of population sizes examined. More importantly, populations will have to be maintained at sizes of >2000 individuals to maintain population fitness at levels compatible with the conservation goal of long-term persistence. This approach to estimating minimum viable population size provides estimates that are in general agreement with those from numerous other studies and strengthens the argument that conservation efforts should ultimately aim at maintaining populations of several thousand individuals to ensure long-term persistence.     

Testing the effectiveness of capture mark recapture population estimation techniques using a computer simulation with known population size

Estimation of small mammal population sizes is important for monitoring ecosystem condition and for conservation. Here, we test the accuracy of standard methods of population size estimation using Capture-Mark-Recapture (CMR) on a simulated population of agents. The use of a computer simulation allows complete control of population sizes and behaviors, thereby avoiding assumptions that may be violated in real populations. We find that the recommended protocol for CMR sampling, using uniformly distributed traps, consistently overestimates population sizes by as much as 100% when studies are conducted over only two trapping periods. More than 20 trapping periods are required before this method, or that of placing traps randomly, gives an accurate estimation of population size (i.e., within a 95% confidence limit of the actual value). Non-random sampling, by placing traps on runways used by small mammals, produces the most accurate, and least variable, estimates of population. However, we show that around 10 trapping periods are still required to produce an accurate population estimate using this method. Given that most real populations do not comply with the ‘ideal’ assumptions made by CMR, we suggest that population estimates based on CMR may be fundamentally flawed, and recommend that protocols for CMR population estimation methods may need revising.     

Intraspecific Chromosome Number Variation: a Neglected Threat to the Conservation of Rare Plants

Abstract: The effectiveness of rare plant conservation will increase when life history, demographic, and genetic data are considered simultaneously. Inbreeding depression is a widely recognized genetic concern in rare plant conservation, and the mixing of genetically diverse populations in restoration efforts is a common remedy. Nevertheless, if populations with unrecognized intraspecific chromosome variation are crossed, progeny fitness losses will range from partial to complete sterility, and reintroductions and population augmentation of rare plants may fail. To assess the current state of cytological knowledge of threatened and endangered plants in the continental United States, we searched available resources for chromosome counts. We also reviewed recovery plans to discern whether recovery criteria potentially place listed species at risk by requiring reintroductions or population augmentation in the absence of cytological information. Over half the plants lacked a chromosome count, and when a taxon did have a count it generally originated from a sampling intensity too limited to detect intraspecific chromosome variation. Despite limited past cytological sampling, we found 11 plants with documented intraspecific cytological variation, while 8 others were ambiguous for intraspecific chromosome variation. Nevertheless, only one recovery plan addressed the chromosome differences. Inadequate within‐species cytological characterization, incomplete sampling among listed taxa, and the prevalence of interspecific and intraspecific chromosome variation in listed genera, suggests that other rare plants are likely to have intraspecific chromosome variation. Nearly 90% of all recovery plans called for reintroductions or population augmentation as part of recovery criteria despite the dearth of cytological knowledge. We recommend screening rare plants for intraspecific chromosome variation before reintroductions or population augmentation projects are undertaken to safeguard against inadvertent mixtures of incompatible cytotypes.     

Characterizing Core and Corridor Use by Nubian Ibex in the Negev Desert, Israel

Abstract: Corridors have become a key element of conservation efforts. To identify and characterize corridors used by Nubian ibex ( Capra ibex nubiana ), we analyzed sighting data recorded for the past 20 years in the Israel Nature Reserves Authority data banks. We categorized each recorded sighting as belonging to a corridor or a core zone based on the total number of sightings in its vicinity. We identified three main core populations, a natural corridor connecting two of them, and a corridor that was not evident connecting the third population. Ibex inside and outside core zones were sighted on steeper terrain than expected by chance; this was more pronounced inside the core zones than outside them. We hypothesize that because ibex outside core zones must move rapidly and directionally across unfamiliar habitats, they must use more moderate terrain. Although ibex sightings in core zones were mostly in the vicinity of water, we found no relationship between the location of the sightings and proximity to water sources in the corridor. Hence, water does not appear to be an important factor in movement through corridors. It was more common to observe ibex out of the core zones during the summer. Males, which can be twice the size of females, were found traveling alone in corridors more often than in core zones. There was no difference between males and females in the steepness of terrain in which they were sighted outside core zones. Our data show that protecting ibex habitat in core zones and corridors is important to ibex conservation in Israel's arid zones. In addition, protecting this habitat may benefit other rock-dwelling species in the area.     

Testing the Environmental Kuznets Curve Hypothesis with Bird Populations as Habitat-Specific Environmental Indicators: Evidence from Canada

Abstract:  The traditional environmental Kuznets curve (EKC) hypothesis postulates that environmental degradation follows an inverted U-shaped relationship with gross domestic product (GDP) per capita. We tested the EKC hypothesis with bird populations in 5 different habitats as environmental quality indicators. Because birds are considered environmental goods, for them the EKC hypothesis would instead be associated with a U-shaped relationship between bird populations and GDP per capita. In keeping with the literature, we included other variables in the analysis—namely, human population density and time index variables (the latter variable captured the impact of persistent and exogenous climate and/or policy changes on bird populations over time). Using data from 9 Canadian provinces gathered over 37 years, we used a generalized least-squares regression for each bird habitat type, which accounted for the panel structure of the data, the cross-sectional dependence across provinces in the residuals, heteroskedasticity, and fixed- or random-effect specifications of the models. We found evidence that supports the EKC hypothesis for 3 of the 5 bird population habitat types. In addition, the relationship between human population density and the different bird populations varied, which emphasizes the complex nature of the impact that human populations have on the environment. The relationship between the time-index variable and the different bird populations also varied, which indicates there are other persistent and significant influences on bird populations over time. Overall our EKC results were consistent with those found for threatened bird species, indicating that economic prosperity does indeed act to benefit some bird populations.