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.

Eastern temperate North Pacific offshore killer whales (Orcinus orca): Occurrence, movements, and insights into feeding ecology

Beginning in the late 1980s, large groups of previously unidentified killer whales (Orcinus orca) were sighted off the west coast of Vancouver Island and in the Queen Charlotte Islands, British Columbia. Scientists working in this region produced two killer whale photo-identification catalogues that included both transient (mammal-eating) whales and 65 individual whales that investigators believed represented a distinct killer whale community (Ford et al. 1992, Heise et al. 1993). It was thought that these killer whales maintained a generally offshore distribution and were provisionally termed “offshores”; a term that has since been used as a population identifier for the eastern temperate North Pacific offshore killer whale population. Then in September 1992, 75 unidentified whales entered the Strait of Juan de Fuca just south and east of Victoria, British Columbia (Walters et al. 1992). Although most of these whales had not been seen before, two were matched to killer whales in the Queen Charlotte photo-identification catalogue (Ford et al. 1992, Heise et al. 1993) and were thus listed as “offshore” killer whales. During a similar time period, other large groups of killer whales, previously unidentified, were also being sighted off Alaska and California (Dahlheim et al. 1997; Nancy Black and Alisa Schulman- Janiger, unpublished data, respectively).

Killer Whale Predation on Belugas in Cook Inlet, Alaska: Implications for a Depleted Population

Killer whale predation on belugas in Cook Inlet, Alaska, has become a concern since the decline of these belugas was documented during the 1990s. Accordingly, killer whale sightings were compiled from systematic surveys, observer databases, and anecdotal accounts. Killer whales have been relatively common in lower Cook Inlet (at least 100 sightings from 1975 to 2002), but in the upper Inlet, north of Kalgin Island, sightings were infrequent (18 in 27 yr), especially prior to the 1990s. Beach cast beluga carcasses with teeth marks and missing flesh also provided evidence of killer whale predation. Most observed killer whale/beluga interactions were in the upper Inlet. During 11 of 15 observed interactions, belugas were obviously injured or killed, either through direct attacks or indirectly as a result of stranding. Assuming at least one beluga mortality occurred during the other four encounters, we can account for 21 belugas killed between 1985 and 2002. This would suggest a minimum estimate of roughly l/yr and does not include at least three instances where beluga calves accompanied an adult that was attacked.

Harbor Porpoise (Phocoena phocoena) Abundance in Alaska: Bristol Bay to Southeast Alaska, 1991-1993

Between 1991 and 1993, Alaska harbor porpoise (Phocoena phocoena) abundance was investigated during aerial surveys throughout much of the coastal and offshore waters from Bristol Bay in the eastern Bering Sea to Dixon Entrance in Southeast Alaska. Line-transect methodology was used, and only those observations made during optimal conditions were analyzed. Survey data indicated densities of 4.48 groups/100 km2, or approximately 3,531 harbor porpoises (95% C.I. 2,206-5,651) in Bristol Bay and 0.54 groups/100 km2, or 136 harbor porpoises (95% C.I. 11-1,645) for Cook Inlet. Efforts off Kodiak Island resulted in densities of 1.85 groups/100 km2, or an abundance estimate of 740 (95% C.I. 259-2,115). Surveys off the south side of the Alaska Peninsula found densities of 2.03 groups/100 km2 and an abundance estimate of 551 (95% C.I. 423-719). Surveys of offshore waters from Prince William Sound to Dixon Entrance yielded densities of 4.02 groups/100 km’ and an abundance estimate of 3,982 (95% C.I. 2,567-6,177). Combining all years and areas yielded an uncorrected density estimate of 3.82 porpoises per 100 km2, resulting in an abundance estimate of 8,940 porpoises (CV = 13.8%) with a 95% confidence interval of 6,746-11,848. Using correction factors from other studies to adjust for animals missed by observers, the total number of Alaska harbor porpoises is probably three times this number.

Antler Size of Alaskan Moose Alces alces gigas: Effects of Population Density, Hunter Harvest and Use of Guides

Moose Alces alces gigas in Alaska, USA, exhibit extreme sexual dimorphism, with adult males possessing large, elaborate antlers. Antler size and conformation are influenced by age, nutrition and genetics, and these bony structures serve to establish social rank and affect mating success. Population density, combined with anthropogenic effects such as harvest, is thought to influence antler size. Antler size increased as densities of moose decreased, ostensibly a density-dependent response related to enhanced nutrition at low densities. The vegetation type where moose were harvested also affected antler size, with the largest-antlered males occupying more open habitats. Hunts with guides occurred in areas with low moose density, minimized hunter interference and increased rates of success. Such hunts harvested moose with larger antler spreads than did non-guided hunts. Knowledge and abilities allowed guides to satisfy demands of trophy hunters, who are an integral part of the Alaskan economy. Heavy harvest by humans was also associated with decreased antler size of moose, probably via a downward shift in the age structure of the population resulting in younger males with smaller antlers. Nevertheless, density-dependence was more influential than effects of harvest on age structure in determining antler size of male moose. Indeed, antlers are likely under strong sexual selection, but we demonstrate that resource availability influenced the distribution of these sexually selected characters across the landscape. We argue that understanding population density in relation to carrying capacity (K) and the age structure of males is necessary to interpret potential consequences of harvest on the genetics of moose and other large herbivores. Our results provide researchers and managers with a better understanding of variables that affect the physical condition, antler size, and perhaps the genetic composition of populations, which may be useful in managing and modeling moose populations.

Space–time zero-inflated count models of Harbor seals

Environmental data are spatial, temporal, and often come with many zeros. In this paper, we included space–time random effects in zero-inflated Poisson (ZIP) and ‘hurdle’ models to investigate haulout patterns of harbor seals on glacial ice. The data consisted of counts, for 18 dates on a lattice grid of samples, of harbor seals hauled out on glacial ice in Disenchantment Bay, near Yakutat, Alaska. A hurdle model is similar to a ZIP model except it does not mix zeros from the binary and count processes. Both models can be used for zero-inflated data, and we compared space–time ZIP and hurdle models in a Bayesian hierarchical model. Space–time ZIP and hurdle models were constructed by using spatial conditional autoregressive (CAR) models and temporal first-order autoregressive (AR(1)) models as random effects in ZIP and hurdle regression models. We created maps of smoothed predictions for harbor seal counts based on ice density, other covariates, and spatio-temporal random effects. For both models predictions around the edges appeared to be positively biased. The linex loss function is an asymmetric loss function that penalizes overprediction more than underprediction, and we used it to correct for prediction bias to get the best map for space–time ZIP and hurdle models.

Spatial statistical models that use flow and stream distance

We develop spatial statistical models for stream networks that can estimate relationships between a response variable and other covariates, make predictions at unsampled locations, and predict an average or total for a stream or a stream segment. There have been very few attempts to develop valid spatial covariance models that incorporate flow, stream distance, or both. The application of typical spatial autocovariance functions based on Euclidean distance, such as the spherical covariance model, are not valid when using stream distance. In this paper we develop a large class of valid models that incorporate flow and stream distance by using spatial moving averages. These methods integrate a moving average function, or kernel, against a white noise process. By running the moving average function upstream from a location, we develop models that use flow, and by construction they are valid models based on stream distance. We show that with proper weighting, many of the usual spatial models based on Euclidean distance have a counterpart for stream networks. Using sulfate concentrations from an example data set, the Maryland Biological Stream Survey (MBSS), we show that models using flow may be more appropriate than models that only use stream distance. For the MBSS data set, we use restricted maximum likelihood to fit a valid covariance matrix that uses flow and stream distance, and then we use this covariance matrix to estimate fixed effects and make kriging and block kriging predictions.

Long-term recovery patterns of arctic tundra after winter seismic exploration

In response to the increasing global demand for energy, oil exploration and development are expanding into frontier areas of the Arctic, where slow-growing tundra vegetation and the underlying permafrost soils are very sensitive to disturbance. The creation of vehicle trails on the tundra from seismic exploration for oil has accelerated in the past decade, and the cumulative impact represents a geographic footprint that covers a greater extent of Alaska’s North Slope tundra than all other direct human impacts combined. Seismic exploration for oil and gas was conducted on the coastal plain of the Arctic National Wildlife Refuge, Alaska, USA, in the winters of 1984 and 1985. This study documents recovery of vegetation and permafrost soils over a two-decade period after vehicle traffic on snow-covered tundra. Paired permanent vegetation plots (disturbed vs. reference) were monitored six times from 1984 to 2002. Data were collected on percent vegetative cover by plant species and on soil and ground ice characteristics. We developed Bayesian hierarchical models, with temporally and spatially autocorrelated errors, to analyze the effects of vegetation type and initial disturbance levels on recovery patterns of the different plant growth forms as well as soil thaw depth. Plant community composition was altered on the trails by species-specific responses to initial disturbance and subsequent changes in substrate. Long-term changes included increased cover of graminoids and decreased cover of evergreen shrubs and mosses. Trails with low levels of initial disturbance usually improved well over time, whereas those with medium to high levels of initial disturbance recovered slowly. Trails on ice-poor, gravel substrates of riparian areas recovered better than those on ice-rich loamy soils of the uplands, even after severe initial damage. Recovery to pre-disturbance communities was not possible where trail subsidence occurred due to thawing of ground ice. Previous studies of disturbance from winter seismic vehicles in the Arctic predicted short-term and mostly aesthetic impacts, but we found that severe impacts to tundra vegetation persisted for two decades after disturbance under some conditions. We recommend management approaches that should be used to prevent persistent tundra damage.

Implantation of Subcutaneous Radio Transmitters in the Harbor Seal (Phoca vitulina)

Radio telemetry has become a standard tool for studying the behavior, physiology, life history traits, and population dynamics of marine mammals. Radio transmitters typically are attached to the hind flippers of pinnipeds or glued to the fur using marine epoxy or other cyanocrylare adhesives (Fedak et al. 1983, Bengtson 1993, Jeffries et al. 1993). Longterm data acquisition is difficult, however, because radio-flipper transmitters commonly tear from the webbing of the flipper and instruments that are glued to the fur are shed during the seasonal molt.

Distance software: design and analysis of distance sampling surveys for estimating population size

1. Distance sampling is a widely used technique for estimating the size or density of biological populations. Many distance sampling designs and most analyses use the software Distance. 2. We briefly review distance sampling and its assumptions, outline the history, structure and capabilities of Distance, and provide hints on its use. 3. Good survey design is a crucial prerequisite for obtaining reliable results. Distance has a survey design engine, with a built-in geographic information system, that allows properties of different proposed designs to be examined via simulation, and survey plans to be generated. 4. A first step in analysis of distance sampling data is modeling the probability of detection. Distance contains three increasingly sophisticated analysis engines for this: conventional distance sampling, which models detection probability as a function of distance from the transect and assumes all objects at zero distance are detected; multiple-covariate distance sampling, which allows covariates in addition to distance; and mark–recapture distance sampling, which relaxes the assumption of certain detection at zero distance. 5. All three engines allow estimation of density or abundance, stratified if required, with associated measures of precision calculated either analytically or via the bootstrap. 6. Advanced analysis topics covered include the use of multipliers to allow analysis of indirect surveys (such as dung or nest surveys), the density surface modeling analysis engine for spatial and habitat-modeling, and information about accessing the analysis engines directly from other software. 7. Synthesis and applications. Distance sampling is a key method for producing abundance and density estimates in challenging field conditions. The theory underlying the methods continues to expand to cope with realistic estimation situations. In step with theoretical developments, state-of- the-art software that implements these methods is described that makes the methods accessible to practicing ecologists.

Double-Observer Line Transect Methods: Levels of Independence

Double-observer line transect methods are becoming increasingly widespread, especially for the estimation of marine mammal abundance from aerial and shipboard surveys when detection of animals on the line is uncertain. The resulting data supplement conventional distance sampling data with two-sample mark–recapture data. Like conventional mark–recapture data, these have inherent problems for estimating abundance in the presence of heterogeneity. Unlike conventional mark–recapture methods, line transect methods use knowledge of the distribution of a covariate, which affects detection probability (namely, distance from the transect line) in inference. This knowledge can be used to diagnose unmodeled heterogeneity in the mark–recapture component of the data. By modeling the covariance in detection probabilities with distance, we show how the estimation problem can be formulated in terms of different levels of independence. At one extreme, full independence is assumed, as in the Petersen estimator (which does not use distance data); at the other extreme, independence only occurs in the limit as detection probability tends to one. Between the two extremes, there is a range of models, including those currently in common use, which have intermediate levels of independence. We show how this framework can be used to provide more reliable analysis of double-observer line transect data. We test the methods by simulation, and by analysis of a dataset for which true abundance is known. We illustrate the approach through analysis of minke whale sightings data from the North Sea and adjacent waters.

Health Assessment of Weddell Seals, Leptonychotes weddellii, in McMurdo Sound, Antarctica

The demography of Weddell seals in eastern McMurdo Sound, Antarctica, has been well studied during the past three decades (e.g. Stirling 1971; Siniff et al. 1977; Testa and Siniff 1987; Hastings and Testa 1998; Gelatt et al. 2001). Detailed life-history data are available on thousands of seals tagged as pups in McMurdo Sound, making this population a rich resource for wildlife health studies because health parameters can be evaluated in the light of reproductive histories and genetic relationships of several generations of tagged seals. Recently, evidence of exposure to diseases generally associated with domestic animals and feral wildlife has been detected in Antarctic wildlife (Austin and Webster 1993; Olsen et al. 1996; Gardner et al. 1997; Retamal et al. 2000; Foster et al. 2002) and this has generated concern and debate regarding the risks of disease introduction to Antarctic wildlife. Antibodies to viruses that have caused large die-offs in phocids in other areas of the world have been detected in Weddell seals (Bengtson et al. 1991), and there is a historical report of a mass die-off of crabeater seals that may have had a viral etiology (Laws and Taylor 1957).

Development and application of DNA techniques for validating and improving pinniped diet estimates

Polymerase chain reaction techniques were developed and applied to identify DNA from .40 species of prey contained in fecal (scat) soft-part matrix collected at terrestrial sites used by Steller sea lions (Eumetopias jubatus) in British Columbia and the eastern Aleutian Islands, Alaska. Sixty percent more fish and cephalopod prey were identified by morphological analyses of hard parts compared with DNA analysis of soft parts (hard parts identified higher relative proportions of Ammodytes sp., Cottidae, and certain Gadidae). DNA identified 213 prey occurrences, of which 75 (35%) were undetected by hard parts (mainly Salmonidae, Pleuronectidae, Elasmobranchii, and Cephalopoda), and thereby increased species occurrences by 22% overall and species richness in 44% of cases (when comparing 110 scats that amplified prey DNA). Prey composition was identical within only 20% of scats. Overall, diet composition derived from both identification techniques combined did not differ significantly from hard-part identification alone, suggesting that past scat-based diet studies have not missed major dietary components. However, significant differences in relative diet contributions across scats (as identified using the two techniques separately) reflect passage rate differences between hard and soft digesta material and highlight certain hypothesized limitations in conventional morphological-based methods (e.g., differences in resistance to digestion, hard part regurgitation, partial and secondary prey consumption), as well as potential technical issues (e.g., resolution of primer efficiency and sensitivity and scat subsampling protocols). DNA analysis of salmon occurrence (from scat soft-part matrix and 238 archived salmon hard parts) provided species-level taxonomic resolution that could not be obtained by morphological identification and showed that Steller sea lions were primarily consuming pink (Oncorhynchus gorbuscha) and chum (Oncorhynchus keta) salmon. Notably, DNA from Atlantic salmon (Salmo salar) that likely originated from a distant fish farm was also detected in two scats from one site in the eastern Aleutian Islands. Overall, molecular techniques are valuable for identifying prey in the fecal remains of marine predators. Combining DNA and hard-part identification will effectively alleviate certain predicted biases and will ultimately enhance measures of diet richness, fisheries interactions (especially salmon-related ones), and the ecological role of pinnipeds and other marine predators, to the benefit of marine wildlife conservationists and fisheries managers.

Biology and Impacts of Pacific Island Invasive Species. 5. Eleutherodactylus coqui, the Coqui Frog (Anura: Leptodactylidae)

The nocturnal, terrestrial frog Eleutherodactylus coqui, known as the Coqui, is endemic to Puerto Rico and was accidentally introduced to Hawai‘i via nursery plants in the late 1980s. Over the past two decades E. coqui has spread to the four main Hawaiian Islands, and a major campaign was launched to eliminate and control it. One of the primary reasons this frog has received attention is its loud mating call (85–90 dB at 0.5 m). Many homeowners do not want the frogs on their property, and their presence has influenced housing prices. In addition, E. coqui has indirectly impacted the floriculture industry because customers are reticent to purchase products potentially infested with frogs. Eleutherodactylus coqui attains extremely high densities in Hawai‘i, up to 91,000 frogs ha-1, and can reproduce year-round, once every 1–2 months, and become reproductive around 8–9 months. Although the Coqui has been hypothesized to potentially compete with native insectivores, the most obvious potential ecological impact of the invasion is predation on invertebrate populations and disruption of associated ecosystem processes. Multiple forms of control have been attempted in Hawai‘i with varying success. The most successful control available at this time is citric acid. Currently, the frog is established throughout the island of Hawai‘i but may soon be eliminated on the other Hawaiian Islands via control efforts. Eradication is deemed no longer possible on the island of Hawai‘i.

Comparison of Endocrine Response to Stress Between Captive-Raised and Wild-Caught Bighorn Sheep

Stress hormones in Rocky Mountain bighorn sheep (Ovis canadensis canadensis), produced in response to environmental changes, road development, or high population density, may impact their immune systems to a threshold level that predisposes them to periodic, large-scale mortality. We compared the stress response to a novel environmental situation and repeated handling between bighorn sheep born and raised in captivity (CR) and bighorn sheep born in the wild (WC) and brought into captivity. We measured plasma epinephrine, norepinephrine, cortisol, and fecal glucocorticoid metabolites (FGM). Three weeks after each group’s arrival we used a one-time drop-net event to elicit an acute stress response, and we collected blood samples from each sheep over 35 minutes, as well as one fecal sample. We collected blood and fecal samples from both groups on 7 other occasions over the subsequent 6 months. We also collected fecal samples from the pen at approximately 24-hour intervals for 3 days following every handling event to monitor the stress response to handling. We found that CR sheep had a stronger autonomic nervous system response than WC sheep, as measured by epinephrine and norepinephrine levels, but we found a very similar hypothalamic–pituitary–adrenal axis (HPA) response, measured by cortisol levels, to the acute stress event of a drop-net restraint. We also found that once the WC sheep had acclimated, as indicated by the return to the initial baseline FGM levels within 12 weeks, the CR and WC groups’ HPA responses to sampling events were not significantly different from one another. Fecal samples can provide a noninvasive mechanism for managers to monitor baseline FGM for a given herd. Using long-term monitoring of FGM rather than values from a single point in time may allow managers to correlate these levels to outside influences on the herd and better understand the impacts of management changes, population density, or increased human developments on the health of the sheep population.