Crop, Native Vegetation, and Biofuels: Response of White- Tailed Deer to Changing Management Priorities

The expansion of the cellulosic biofuels industry throughout the United States has broad-scale implications for wildlife management on public and private lands. Knowledge is limited on the effects of reverting agriculture to native grass, and vice versa, on size of home range and habitat use of white-tailed deer (Odocoileus virginianus). We followed 68 radio-collared female deer from 1991 through 2004 that were residents of DeSoto National Wildlife Refuge (DNWR) in eastern Nebraska, USA. The refuge was undergoing conversion of vegetation out of row-crop agriculture and into native grass, forest, and emergent aquatic vegetation. Habitat in DNWR consisted of 30% crop in 1991 but removing crops to establish native grass and wetland habitat at DNWR resulted in a 44% reduction in crops by 2004. A decrease in the amount of crops on DNWR contributed to a decline in mean size of annual home range from 400 ha in 1991 to 200 ha in 2005 but percentage of crops in home ranges increased from 21% to 29%. Mean overlap for individuals was 77% between consecutive annual home ranges across 8 years, regardless of crop availability. Conversion of crop to native habitat will not likely result in home range abandonment but may impact disease transmission by increasing rates of contact between deer social groups that occupy adjacent areas. Future research on condition indices or changes in population parameters (e.g., recruitment) could be incorporated into the study design to assess impacts of habitat conversion for biofuel production.

A New Birth-Interval Approach to Estimating Demographic Parameters of Humpback Whales

A demographic model is developed based on interbirth intervals and is applied to estimate the population growth rate of humpback whales (Megaptera novaeangliae) in the Gulf of Maine. Fecundity rates in this model are based on the probabilities of giving birth at time t after a previous birth and on the probabilities of giving birth first at age x. Maximum likelihood methods are used to estimate these probabilities using sighting data collected for individually identified whales. Female survival rates are estimated from these same sighting data using a modified Jolly–Seber method. The youngest age at first parturition is 5 yr, the estimated mean birth interval is 2.38 yr (SE = 0.10 yr), the estimated noncalf survival rate is 0.960 (SE = 0.008), and the estimated calf survival rate is 0.875 (SE = 0.047). The population growth rate (l) is estimated to be 1.065; its standard error is estimated as 0.012 using a Monte Carlo approach, which simulated sampling from a hypothetical population of whales. The simulation is also used to investigate the bias in estimating birth intervals by previous methods. The approach developed here is applicable to studies of other populations for which individual interbirth intervals can be measured.

Demography and population viability of polar bears in the Gulf of Boothia, Nunavut

We estimated demographic parameters and harvest risks for polar bears (Ursus maritimus) inhabiting the Gulf of Boothia, Nunavut, from 1976 to 2000. We computed survival and abundance from capture–recapture and recovery data (630 marks) using a Burnham joint live–dead model implemented in program MARK. Annual mean total survival (including harvest) was 0.889 ± 0.179 ( x ± 1 SE) for cubs, 0.883 ± 0.087 for subadults (ages 1–4), 0.919 ± 0.044 for adult females, and 0.917 ± 0.041 for adult males. Abundance in the last 3 yr of study was 1,592 ± 361 bears. Mean size of newborn litters was 1.648 ± 0.098 cubs. By age 7, 0.97 ± 0.30 of available females were producing litters. Harvest averaged 38.4 ± 4.2 bears/year in the last 5 yr of study; however, the 2002–2007 kill averaged 56.4 bears/yr. We used a harvested Population Viability Analysis (PVA) to examine impacts of increasing rates of harvest. We estimated the current population growth rate, λH, to be 1.025 ± 0.032. Although this suggests the population is growing, progressive environmental changes may require more frequent population inventory studies to maintain the same levels of harvest risk.

Improved Regression Estimation of a Multivariate Relationship with Population Data on the Bivariate Relationship

Regression coefficients specify the partial effect of a regressor on the dependent variable. Sometimes the bivariate or limited multivariate relationship of that regressor variable with the dependent variable is known from population-level data. We show here that such population- level data can be used to reduce variance and bias about estimates of those regression coefficients from sample survey data. The method of constrained MLE is used to achieve these improvements. Its statistical properties are first described. The method constrains the weighted sum of all the covariate-specific associations (partial effects) of the regressors on the dependent variable to equal the overall association of one or more regressors, where the latter is known exactly from the population data. We refer to those regressors whose bivariate or limited multivariate relationships with the dependent variable are constrained by population data as being ‘‘directly constrained.’’ Our study investigates the improvements in the estimation of directly constrained variables as well as the improvements in the estimation of other regressor variables that may be correlated with the directly constrained variables, and thus ‘‘indirectly constrained’’ by the population data. The example application is to the marital fertility of black versus white women. The difference between white and black women’s rates of marital fertility, available from population-level data, gives the overall association of race with fertility. We show that the constrained MLE technique both provides a far more powerful statistical test of the partial effect of being black and purges the test of a bias that would otherwise distort the estimated magnitude of this effect. We find only trivial reductions, however, in the standard errors of the parameters for indirectly constrained regressors.

Contributions of demography and dispersal parameters to the spatial spread of a stage-structured insect invasion

Stage-structured models that integrate demography and dispersal can be used to identify points in the life cycle with large effects on rates of population spatial spread, information that is vital in the development of containment strategies for invasive species. Current challenges in the application of these tools include: (1) accounting for large uncertainty in model parameters, which may violate assumptions of ‘‘local’’ perturbation metrics such as sensitivities and elasticities, and (2) forecasting not only asymptotic rates of spatial spread, as is usually done, but also transient spatial dynamics in the early stages of invasion. We developed an invasion model for the Diaprepes root weevil (DRW; Diaprepes abbreviatus [Coleoptera: Curculionidae]), a generalist herbivore that has invaded citrus-growing regions of the United States. We synthesized data on DRW demography and dispersal and generated predictions for asymptotic and transient peak invasion speeds, accounting for parameter uncertainty. We quantified the contributions of each parameter toward invasion speed using a ‘‘global’’ perturbation analysis, and we contrasted parameter contributions during the transient and asymptotic phases. We found that the asymptotic invasion speed was 0.02–0.028 km/week, although the transient peak invasion speed (0.03– 0.045 km/week) was significantly greater. Both asymptotic and transient invasions speeds were most responsive to weevil dispersal distances. However, demographic parameters that had large effects on asymptotic speed (e.g., survival of early-instar larvae) had little effect on transient speed. Comparison of the global analysis with lower-level elasticities indicated that local perturbation analysis would have generated unreliable predictions for the responsiveness of invasion speed to underlying parameters. Observed range expansion in southern Florida (1992–2006) was significantly lower than the invasion speed predicted by the model. Possible causes of this mismatch include overestimation of dispersal distances, demographic rates, and spatiotemporal variation in parameter values. This study demonstrates that, when parameter uncertainty is large, as is often the case, global perturbation analyses are needed to identify which points in the life cycle should be targets of management. Our results also suggest that effective strategies for reducing spread during the asymptotic phase may have little effect during the transient phase. Includes Appendix.