Modeling Dynamics of Habitat-Interior Bird Populations in Fragmented Landscapes

Abstract: A stochastic computer model was used to examine the effects of varying degrees of habitat fragmentation on the dynamics of a hypothetical population of forest-interior bid. The primary demographic parameter that influenced the population's dynamics was fecundity, which varied as a function of how far a birds territory was from an ecological edge. As our model landscape became more fragmented the proportion of forest habitat that was near edges increased geometrically, and the population's overall fecundity dropped as a result. The model demonstrates that impaired reproduction in a fragmented landscape is, by itself a sufficient disruption of the population's dynamics to generate population declines and shifts in distribution similar to those observed in the fragmented forests of southern Wisconsin. Without immigration of recruits from other regions where reproduction is better, habitat-interior populations in a severely fragmented landscape can become locally extinct.     

Interactive effects of prey and p,p'-DDE on burrowing owl population dynamics.

We used population models to explore the effects of the organochlorine contaminant p,p'-DDE and fluctuations in vole availability on the population dynamics of Burrowing Owls (Athene cunicularia). Previous work indicated an interaction between low biomass of voles in the diet and moderate levels of p,p'-DDE in Burrowing Owl eggs that led to reproductive impairment. We constructed periodic and stochastic matrix models that incorporated three vole population states observed in the field: average, peak, and crash years. We modeled varying frequencies of vole crash years and a range of impairment of owl demographic rates in vole crash years. Vole availability had a greater impact on owl population growth rate than did reproductive impairment if vole populations peaked and crashed frequently. However, this difference disappeared as the frequency of vole crash years declined to once per decade. Fecundity, the demographic rate most affected by p,p'-DDE, had less impact on population growth rate than adult or juvenile survival. A life table response experiment of time-invariant matrices for average, peak, and crash vole conditions showed that low population growth under vole crash conditions was due to low adult and juvenile survival rates, whereas the extremely high population growth under vole peak conditions was due to increased fecundity. Our results suggest that even simple models can provide useful insights into complex ecological interactions. This is particularly valuable when temporal or spatial scales preclude manipulative experimental work in the field or laboratory.     

Forest age structure as indicator of boreal forest sustainability under alternative management and fire regimes: A landscape level sensitivity analysis

Effective forest ecosystem-based management requires a thorough understanding of the interactions between anthropogenic and natural disturbance processes over larger spatial and temporal scales than stands and rotation ages. Because harvesting does not preclude fire, it is important to evaluate the combined effects of harvesting and fire on forest age structure, a coarse indicator of forest ecosystem state. We performed a sensitivity analysis of landscape scale effects of forest management (strategy, harvest rate and access cost) and fire regime (fire return interval and extent) in terms of combined impacts on forest stand age-class structure on a study area of 3.5 million hectares of boreal forest of Québec. A series of scenarios were simulated over 500 years and replicated 30 times using a previously reported spatially explicit landscape model. Within the parameter space of our sensitivity analysis, we found that harvest rate, fire return interval and management strategy were the most significant parameters affecting stand age-class distribution across the landscape. The former are not so surprising, given that they combine to produce an overall disturbance rate, but the latter shows that the resulting impact on age-class structure can be influenced to some degree through management objectives. A harvesting strategy of clearcutting for sustained timber supply, using a harvest rotation based on minimum merchantable age (approximately 100 years in this analysis), creates a trend for the stand age-class distribution away from the expected range of natural variation for the study area. Within the scope of our simulations, alternative management strategies with extended harvest rotation age proved the most robust forest management practice to absorb variations in fire regime.     

Use of Integrated Modeling to Enhance Estimates of Population Dynamics Obtained from Limited Data

Abstract:  Demographic data of rare and endangered species are often too sparse to estimate vital rates and population size with sufficient precision for understanding population growth and decline. Yet, the combination of different sources of demographic data into one statistical model holds promise. We applied Bayesian integrated population modeling to demographic data from a colony of the endangered greater horseshoe bats (Rhinolophus ferrumequinum) . Available data were the number of subadults and adults emerging from the colony roost at dusk, the number of newborns from 1991 to 2005, and recapture data of subadults and adults from 2004 and 2005. Survival rates did not differ between sexes, and demographic rates remained constant across time. The greater horseshoe bat is a long-lived species with high survival rates (first year: 0.49 [SD 0.06]; adults: 0.91 [SD 0.02]) and low fecundity (0.74 [SD 0.12]). The yearly average population growth was 4.4% (SD 0.1%) and there were 92 (SD 10) adults in the colony in year 2005. Had we analyzed each data set separately, we would not have been able to estimate fecundity, the estimates of survival would have been less precise, and the estimate of population growth biased. Our results demonstrate that integrated models are suitable for obtaining crucial demographic information from limited data.     

Population growth in snow geese: a modeling approach integrating demographic and survey information

There are few analytic tools available to formally integrate information coming from population surveys and demographic studies. The Kalman filter is a procedure that facilitates such integration. Based on a state-space model, we can obtain a likelihood function for the survey data using a Kalman filter, which we may then combine with a likelihood for the demographic data. In this paper, we used this combined approach to analyze the population dynamics of a hunted species, the Greater Snow Goose (Chen caerulescens atlantica), and to examine the extent to which it can improve previous demographic population models. The state equation of the state-space model was a matrix population model with fecundity and regression parameters relating adult survival and harvest rate estimated in a previous capture-recapture study. The observation equation combined the output from this model with estimates from an annual spring photographic survey of the population. The maximum likelihood estimates of the regression parameters from the combined analysis differed little from the values of the original capture-recapture analysis, though their precision improved. The model output was found to be insensitive to a wide range of coefficient of variation (CV) in fecundity parameters. We found a close match between the surveyed and smoothed population size estimates generated by the Kalman filter over an 18-year period, and the estimated CV of the survey (0.078-0.150) was quite compatible with its assumed value (approximately 0.10). When we used the updated parameter values to predict future population size, the model underestimated the surveyed population size by 18% over a three-year period. However, this could be explained by a concurrent change in the survey method. We conclude that the Kalman filter is a promising approach to forecast population change because it incorporates survey information in a formal way compared with ad hoc approaches that either neglect this information or require some parameter or model tuning.     

Defining Forest Fragmentation by Corridor Width: The Influence of Narrow Forest-Dividing Corridors on Forest-Nesting Birds in Southern New Jersey

In studies of forest fragmentation, a fundamental inconsistency exists in the distance criterion used to define the discreteness of forest fragments. We examined three types of ubiquitous, narrow, forest-dividing corridors for effects that influence the relative abundance and community composition of forest-nesting birds. Fixed-radius (100-meter) point counts were conducted on 54 transects established along three width classes of corridors: unpaved roads (8 meters wide), paved roads (16 meters wide), and powerlines (23 meters wide). Transect locations were distributed equally among corridor edge, forest margin 100 meters from corridor edge, and forest interior 300 meters from corridor edge. Forest-interior species of Neotropical migrants had significantly reduced relative abundances on edge transects along 16- and 23-meter corridors, compared with 8-meter corridors and with forest interior points along all three corridor-width classes. At a landscape scale, the consequences of apparently small reductions in forest area by the presence of narrow forest-dividing corridors may be cumulatively significant for abundances of forest-interior species. Brown-headed Cowbirds were more abundant than 20 of 21 forest-interior Neotropical migrants. We found surprisingly high abundances of cowbirds associated with narrow forest-dividing corridors, especially those with mowed grass. Corridor widths as narrow as 8 meters produce forest fragmentation effects in part by attracting cowbirds and nest predators to corridors and adjacent forest interiors. The most serious implication of this study is that narrow forest-dividing corridors may function as ecological traps for forest-interior Neotropical migrants. We suggest that these widespread corridors may be inconspicuous but important contributors to declines of forest-interior nesting species in eastern North America.     

Effects of Multiple Levels of Social Organization on Survival and Abundance

Abstract: Identifying how social organization shapes individual behavior, survival, and fecundity of animals that live in groups can inform conservation efforts and improve forecasts of population abundance, even when the mechanism responsible for group‐level differences is unknown. We constructed a hierarchical Bayesian model to quantify the relative variability in survival rates among different levels of social organization (matrilines and pods) of an endangered population of killer whales (Orcinus orca). Individual killer whales often participate in group activities such as prey sharing and cooperative hunting. The estimated age‐specific survival probabilities and survivorship curves differed considerably among pods and to a lesser extent among matrilines (within pods). Across all pods, males had lower life expectancy than females. Differences in survival between pods may be caused by a combination of factors that vary across the population's range, including reduced prey availability, contaminants in prey, and human activity. Our modeling approach could be applied to demographic rates for other species and for parameters other than survival, including reproduction, prey selection, movement, and detection probabilities.     

Variation in responses of late-seral herbs to disturbance and environmental stress

Clonal herbs that attain maximum development in late-seral forest are often assumed to have similar responses to disturbance and to be functionally equivalent. However, little is known about the demographic or physiological responses of these plants to disturbance or to the altered conditions of the post-disturbance environment. Following harvest of a mature coniferous forest, we compared abundance, demographic changes, and physiological acclimation of three clonal herbs (Asarum caudatum, Clintonia uniflora, and Pyrola picta) that differ in belowground morphology and leaf longevity. We measured ramet density, leaf area, and demographic variables (survival, clonal growth, flowering, and seedling establishment) before and for two years after harvest, and in adjacent undisturbed forest. Acclimation to increased solar radiation was assessed two years after harvest by measuring leaf mass per unit area (LMA) and chlorophyll a:b ratios of leaves produced in the current year. Although initial declines in abundance were similar, demographic responses indicate that patterns of recovery varied greatly among species. Two years after logging, ramet survival and clonal growth (production of new ramets) of Clintonia were greater in the harvest area than in the forest. Asarum had lower survival in the harvest area, but greater clonal growth, and Pyrola showed no difference in either survival or growth between environments. Only Asarum produced seedlings, although their survival was low in the harvest area. All species had higher LMA in the harvest area, but only Clintonia (with annual leaves) had a higher chlorophyll a:b ratio, suggesting the greatest potential for acclimation to increased light. Our results demonstrate that forest herbs with greater rhizome plasticity and shorter leaf duration have greater potential to acclimate after disturbance than those with rigid architectures and persistent leaves. Thus, species with comparable successional roles can vary substantially in their demographic and physiological responses to disturbance, with potential consequences for long-term recovery.     

An Age-Structured Demographic Model for the Endangered Stephens' Kangaroo Rat

Effective conservation of endangered species often is hampered by inadequate knowledge of demography. We extracted information on survival and fecundity from an 18-month, live-trapping study of Dipodomys stephensi , and from this we developed an age-structured demographic model to assess population viability. Adult Stephens' kangaroo rats persisted longer than juveniles, and adult females persisted longer than adult males. Disappearance rates were high in the first months after initial capture. Thereafter, the fraction of animals persisting decreased slowly and in an approximately linear fashion on a semilogarithmic scale, suggesting age-independent mortality factors such as predation. Juvenile persistence did not differ substantially between two years of strikingly different rainfall. Onset of breeding followed the start of winter rains. Length of the breeding season, average number of litters per female, and the fraction of first-year females breeding were much greater in the year of higher rainfall. We propose a birth-pulse demographic model for D. stephensi that distinguishes juvenile and adult age classes. Temporal environmental variation can be modeled adequately with a constant survivorship schedule and variable fecundity determined by yearly precipitation. Several issues should be resolved, however, before conservation decisions are based on the model. Better estimates of juvenile survivorship are critical, the quantitative relationship between precipitation and fecundity must be determined, and the potential for density dependence and source-sink population dynamics must be evaluated.     

Invasion dynamics of the varnish clam (Nuttallia obscurata): a matrix demographic modeling approach

Marine invaders have become a significant threat to native biodiversity and ecosystem function. In this study, the invasion of the varnish clam (Nuttallia obscurata) in British Columbia, Canada, is investigated using a matrix modeling approach to identify the life history characteristics most crucial for population growth and to investigate population differences. Mark-recapture analyses and field collections from 2003 to 2004 were used to determine individual growth, survival rates, and fecundity for two sites. A multi-state matrix model was used to determine population growth rates and to conduct sensitivity and elasticity analyses. A life table response experiment was also used to determine what life history stage contributed most to observed differences in population growth rates. Population survey data were used in conjunction with the matrix model to determine plausible recruitment levels and to investigate recruitment scenarios. Both populations are currently declining but are likely sustainable because of the pulsed nature of large recruitment events. Survival of larger clams (>40 mm) is the most important for population growth based on elasticity and sensitivity analyses. Adult survival also had the largest influence on observed differences between site-specific population growth rates. The two populations studied differed in recruitment dynamics; one experiencing annual recruitment with higher post-settlement mortality and the other, episodic recruitment and lower post-settlement mortality. The most influential factor for the successful invasion of the varnish clam appears to be survival of the larger size classes. Therefore, any process that decreases adult survival (e.g., predation, commercial harvest) will have the greatest impact on population growth.     

Effects of Harvest of Nontimber Forest Products and Ecological Differences between Sites on the Demography of African Mahogany

Abstract: The demographic impacts of harvesting nontimber forest products (NTFP) have been increasingly studied because of reports of potentially unsustainable harvest. Nevertheless, our understanding of how plant demographic response to harvest is altered by variation in ecological conditions, which is critical for developing realistic sustainable‐use plans, is limited. We built matrix population models to test whether and how variation in ecological conditions affects population responses to harvest. In particular, we examined the effect of bark and foliage harvest on the demography of populations of African mahogany (Khaya senegalensis) in two contrasting ecological regions of Benin, West Africa. K. senegalensis bark and foliage harvest significantly reduced its stochastic population growth rates, but ecological differences between regions had a greater effect on population growth rates than did harvest. The effect of harvest on population growth rates (Δλ) was slightly stronger in the moist than in the drier region. Life‐table response experiments revealed that the mechanism by which harvesting reduced λ differed between ecological regions. Lowered stasis (persistence) of larger life stages lead to a reduction in λ in the drier region, whereas lowered growth of all life stages lowered λ in moist region. Potential strategies to increase population growth rates should include decreasing the proportion of individuals harvested, promoting harvester‐owned plantations of African mahogany, and increasing survival and growth by promoting no‐fire zones in gallery forests. Our results show how population responses to harvest of NTFP may be altered by ecological differences across sites and emphasize the importance of monitoring populations over the climatic range in which they occur to develop more realistic recommendations for conservation.     

Effects of Harvesting Flowers from Shrubs on the Persistence and Abundance of Wild Shrub Populations at Multiple Spatial Extents

Abstract: Wildflower harvesting is an economically important activity of which the ecological effects are poorly understood. We assessed how harvesting of flowers affects shrub persistence and abundance at multiple spatial extents. To this end, we built a process‐based model to examine the mean persistence and abundance of wild shrubs whose flowers are subject to harvest (serotinous Proteaceae in the South African Cape Floristic Region). First, we conducted a general sensitivity analysis of how harvesting affects persistence and abundance at nested spatial extents. For most spatial extents and combinations of demographic parameters, persistence and abundance of flowering shrubs decreased abruptly once harvesting rate exceeded a certain threshold. At larger extents, metapopulations supported higher harvesting rates before their persistence and abundance decreased, but persistence and abundance also decreased more abruptly due to harvesting than at smaller extents. This threshold rate of harvest varied with species’ dispersal ability, maximum reproductive rate, adult mortality, probability of extirpation or local extinction, strength of Allee effects, and carrying capacity. Moreover, spatial extent interacted with Allee effects and probability of extirpation because both these demographic properties affected the response of local populations to harvesting more strongly than they affected the response of metapopulations. Subsequently, we simulated the effects of harvesting on three Cape Floristic Region Proteaceae species and found that these species reacted differently to harvesting, but their persistence and abundance decreased at low rates of harvest. Our estimates of harvesting rates at maximum sustainable yield differed from those of previous investigations, perhaps because researchers used different estimates of demographic parameters, models of population dynamics, and spatial extent than we did. Good demographic knowledge and careful identification of the spatial extent of interest increases confidence in assessments and monitoring of the effects of harvesting. Our general sensitivity analysis improved understanding of harvesting effects on metapopulation dynamics and allowed qualitative assessment of the probability of extirpation of poorly studied species.     

Robust population management under uncertainty for structured population models.

Structured population models are increasingly used in decision making, but typically have many entries that are unknown or highly uncertain. We present an approach for the systematic analysis of the effect of uncertainties on long-term population growth or decay. Many decisions for threatened and endangered species are made with poor or no information. We can still make decisions under these circumstances in a manner that is highly defensible, even without making assumptions about the distribution of uncertainty, or limiting ourselves to discussions of single, infinitesimally small changes in the parameters. Suppose that the model (determined by the data) for the population in question predicts long-term growth. Our goal is to determine how uncertain the data can be before the model loses this property. Some uncertainties will maintain long-term growth, and some will lead to long-term decay. The uncertainties are typically structured, and can be described by several parameters. We show how to determine which parameters maintain long-term growth. We illustrate the advantages of the method by applying it to a Peregrine Falcon population. The U.S. Fish and Wildlife Service recently decided to allow minimal harvesting of Peregrine Falcons after their recent removal from the Endangered Species List. Based on published demographic rates, we find that an asymptotic growth rate lambda > 1 is guaranteed with 5% harvest rate up to 3% error in adult survival if no two-year-olds breed, and up to 11% error if all two-year-olds breed. If a population growth rate of 3% or greater is desired, the acceptable error in adult survival decreases to between 1% and 6% depending of the proportion of two-year-olds that breed. These results clearly show the interactions between uncertainties in different parameters, and suggest that a harvest decision at this stage may be premature without solid data on adult survival and the frequency of breeding by young adults.     

Demographic Responses by Birds to Forest Fragmentation

Abstract:  Despite intensive recent research on the effects of habitat loss and fragmentation on bird populations, our understanding of underlying demographic causes of population declines is limited. We reviewed avian demography in relation to habitat fragmentation. Then, through a meta-analysis, we compared specific demographic responses by forest birds to habitat fragmentation, providing a general perspective of factors that make some species and populations more vulnerable to fragmentation than others. We obtained data from the scientific literature on dispersal, survival, fecundity, and nesting success of birds. Birds were divided into subgroups on the basis of region, nest site, biogeographical history, and migration strategy. Species most sensitive to fragmentation were ground- or open-nesters nesting in shrubs or trees. Residents were equally sensitive to fragmentation in the Nearctic and Palearctic regions, but Nearctic migrants were more sensitive than Palearctic migrants. Old World species were less sensitive than New World species, which was predicted based on the history of forest fragmentation on these two continents. Pairing success was the variable most associated with fragmentation, suggesting an important role of dispersal. Fledgling number or condition, timing of nesting, and clutch size were not associated with sensitivity to fragmentation, suggesting that negative fragmentation effects on birds do not generally result from diminished food resources with increasing level of fragmentation. Future studies on demographic responses of birds to habitat fragmentation would be more effective if based on a combination of measures that can distinguish among the demographic mechanisms underlying population changes related to habitat fragmentation.     

Bet hedging in a guild of desert annuals

Evolutionary bet hedging encapsulates the counterintuitive idea that organisms evolve traits that reduce short-term reproductive success in favor of longer-term risk reduction. It has been widely investigated theoretically, and many putative examples have been cited including practical ones such as the dormancy involved in microbe and weed persistence. However, long-term data on demographic variation from the actual evolutionarily relevant environments have been unavailable to test for its mechanistic relationship to alleged bet hedging traits. I report an association between delayed germination (a bet hedging trait) and risk using a 22-year data set on demographic variation for 10 species of desert annual plants. Species with greater variation in reproductive success (per capita survival from germination to reproduction x per capita fecundity of survivors) were found to have lower average germination fractions. This provides a definitive test using realistic data on demographic variance that confirms the life history prediction for bet hedging. I also showed that the species with greater long-term demographic variation tended to be the ones with greater sensitivity of reproductive success to variation among years in growing-season precipitation.     

The Influence of Model Structure on Conclusions about the Viability and Harvesting of Serengeti Wildebeest

We investigate how the viability and harvestability predicted by population models are affected by details of model construction. Based on this analysis we discuss some of the pitfalls associated with the use of classical statistical techniques for resolving the uncertainties associated with modeling population dynamics. The management of the Serengeti wildebeest (Connochaetes taurinus) is used as a case study. We fitted a collection of age-structured and unstructured models to a common set of available data and compared model predictions in terms of wildebeest viability and harvest. Models that depicted demographic processes in strikingly different ways fitted the data equally well. However, upon further analysis it became clear that models that fit the data equally well could nonetheless have very different management implications. In general, model structure had a much larger effect on viability analysis (e.g., time to collapse) than on optimal harvest analysis (e.g., harvest rate that maximizes harvest). Some modeling decisions, such as including age-dependent fertility rates, did not affect management predictions, but others had a strong effect (e.g., choice of model structure). Because several suitable models of comparable complexity fitted the data equally well, traditional model selection methods based on the parsimony principle were not practical for judging the value of alternative models. Our results stress the need to implement analytical frameworks for population management that explicitly consider the uncertainty about the behavior of natural systems.     

Relating chronic toxicity responses to population-level effects: A comparison of population-level parameters for three salmon species as a function of low-level toxicity

Standard laboratory toxicity tests assess the physiological responses of individual organisms to exposure to toxic substances under controlled conditions. Time and space restrictions often prevent the assessment of population-level responses to a toxic substance. Contaminants can affect various biological functions (e.g. growth, fecundity or behavior), which may alter different demographic traits, leading to population-level impacts. In this study, immune suppression, reproductive dysfunction and somatic growth impairment were examined using life history matrix models for coho salmon (Oncorhynchus kisutch), sockeye salmon (Oncorhynchus nerka) and chinook salmon (Oncorhynchus tshawytscha). Our intent was to gauge the relative magnitude of response to toxic effects among species and between life history stages, not provide a specific estimate of population growth rate or abundance. Effects due to immune suppression were modeled as reductions in age-specific survival. Toxic impacts on reproductive function were modeled as a 10% reduction in reproductive contribution for all reproductively mature age groups. Model runs that examined the effect of somatic growth reduction on population parameters incorporated both survival and reproductive impacts. All impacts were modeled as 10% reductions in the affected population demographic parameters. First-year survival and reproductive impacts produced similar population growth rates (λ), but resulted in different sensitivity and stable age distributions. Modeled somatic growth reduction produced additive effects on survival and reproduction. Toxic stressors producing similar changes in λ did not necessarily produce similar changes in the age distributions. Sensitivity and elasticity analyses demonstrated that changes to the first-year survival rate produced the greatest per-unit effect on λ for each species. Alteration in abundance of mature females also strongly influenced λ. Differences observed between species showed that the number of reproductive ages and time to reproductive maturity were important components for population-level responses. These results emphasize the importance of linking toxicity responses at low concentrations to the demographic traits they affect, and help to highlight the toxicity tests that are more suitable for assessing impacts on the focal species. Additionally, life history modeling is a useful tool for developing testable hypotheses regarding impacts on specific populations as well as for conducting comparisons between populations.     

Cost‐Effectiveness of Translocation Options for a Threatened Waterbird

Abstract: Reintroduction of captive‐reared animals has become increasingly popular in recent decades as a conservation technique, but little is known of how demographic factors affect the success of reintroductions. We believe whether the increase in population persistence associated with reintroduction is sufficient to warrant the cost of rearing and relocating individuals should be considered as well. We examined the trade‐off between population persistence and financial cost of a reintroduction program for Crested Coots (Fulica cristata). This species was nearly extirpated from southern Europe due to unsustainable levels of hunting and reduction in amount and quality of habitat. We used a stochastic, stage‐based, single‐sex, metapopulation model with site‐specific parameters to examine the demographic effects of releasing juveniles or adults in each population for a range of durations. We parameterized the model with data from an unsuccessful reintroduction program in which juvenile captive‐bred Crested Coots were released between 2000 and 2009. Using economic data from the captive‐breeding program, we also determined whether the strategy that maximized abundance coincided with the least expensive strategy. Releasing adults resulted in slightly larger final abundance than the release of nonreproductive juveniles. Both strategies were equally poor in achieving a viable metapopulation, but releasing adults was 2–4 times more expensive than releasing juveniles. To obtain a metapopulation that would be viable for 30 years, fecundity in the wild would need to increase to the values observed in captivity and juvenile survival would need to increase to almost unity. We suggest that the most likely way to increase these vital rates is by increasing habitat quality at release sites.     

Effect of stage-specific vital rates on population growth rates and effective population sizes in an endangered iteroparous plant

Effective population size (N(e)) determines the strength of genetic drift and can influence the level of genetic diversity a population can maintain. Assessing how changes in demographic rates associated with environmental variables and management actions affect N(e) thus can be crucial to the conservation of endangered species. Calculation of N(e) through demographic models makes it possible to use elasticity analyses to study this issue. The elasticity of N(e) to a given vital rate is the proportional change in N(e) associated with a proportional increase in that vital rate. In addition, demographic models can be used to study N(e) and population growth rate (λ) simultaneously. Simultaneous examination is important because some vital rates differ diametrically in their associations with λ and N(e). For example, in some cases increasing these vital rates increases λ and decreases N(e). We used elasticity analysis to study the effect of stage-specific survival and flowering rates on N(e), annual effective population size (N(a)), and λ in seven populations of the endangered plant Austrian dragonhead (Dracocephalum austriacum). In populations with λ ≥ 1, the elasticities of N(e) and N(a) were similar to those of λ. Survival rates of adults were associated with greater elasticities than survival rates of juveniles, flowering rates, or fecundity. In populations with λ < 1, N(e) and N(a) exhibited greater elasticities to juvenile than to adult vital rates. These patterns are similar to those observed in other species with similar life histories. We did not observe contrasting effects of any vital rate on λ and N(e); thus, management actions that increase the λ of populations of Austrian dragonhead will not increase genetic drift. Our results show that elasticity analyses of N(e) and N(a) can complement elasticity analysis of λ. Moreover, such analyses do not require more data than standard matrix models of population dynamics.     

Causasilty among wildfire,ENSO, timber harvest,and urban sprawl: The vector autoregression approach

This article demonstrates the applicability of vector autoregression (VAR) modeling in probing the causality relationships among wildfire, El Niño/Southern Oscillation (ENSO), timber harvest, and urban sprawl in the U.S. The VAR approach allows for the multi-directional, multi-faceted interactions among the variables concerned and enables us to portray the temporal impacts of ENSO, the volume of timber harvested, and urban sprawl on wildfire. The empirical analysis, though intended mainly for illustration, reveals that an individual factor may not affect wildfire activity (number of fires and area burned) when acting alone, but can significantly influence fire activity when coupled with other factors, and that wildfire activity has feedback effects on other variables. The impact of a change in ENSO, the volume of timber harvested, and urban population density on wildfire activity could last two decades with the most noticeable impact occurring in the initial 5–10 years. Though ENSO, timber harvest, and urban sprawl all Granger-cause wildfire activity, the impulse response functions show that wildfire activity is more responsive to urban population density than to the volume of timber harvested or ENSO. Thus, controlling urban sprawl represents another option for wildfire mitigation; and integrative wildfire management is essential.