Seasonal Variation in Reception of Fin Whale Calls at Five Geographic Areas in the North Pacific

In late August 1991 scientists at the National Oceanic and Atmospheric Administration’s (NOAA) National Marine Mammal Laboratory (NMML) and Pacific Marine Environmental Laboratory (PMEL) began a pilot study to investigate the capability of hydrophones from the US. Navy’s fixed array system to detect large whales in the North Pacific by passive reception of their calls. PMEL had previously established a direct data link from five bottom-mounted arrays of the Navy SOSUS (Sound Surveillance System), via the Naval Oceanographic Processing Facility (NOPF) at Whidbey Island, Washington, to study low-level seafloor seismicity (Fox et al. 1994). PMEL subsequently provided NMML tapes of SOSUS hydrophone data from which whale calls were analyzed. As in an analogous study conducted in the North Atlantic (Nishimura and Conlon 1994, Clark 1995, Mellinger and Clark 1995), calls attributable to whales were received at each SOSUS site at rates that varied seasonally (Anonymous 1996).

Use of chemical tracers in assessing the diet and foraging regions of eastern North Pacific killer whales

Top predators in the marine environment integrate chemical signals acquired from their prey that reflect both the species consumed and the regions from which the prey were taken. These chemical tracers—stable isotope ratios of carbon and nitrogen; persistent organic pollutant (POP) concentrations, patterns and ratios; and fatty acid profiles—were measured in blubber biopsy samples from North Pacific killer whales (Orcinus orca) (n = 84) and were used to provide further insight into their diet, particularly for the offshore group, about which little dietary information is available. The offshore killer whales were shown to consume prey species that were distinctly different from those of sympatric resident and transient killer whales. In addition, it was confirmed that the offshores forage as far south as California. Thus, these results provide evidence that the offshores belong to a third killer whale ecotype. Resident killer whale populations showed a gradient in stable isotope profiles from west (central Aleutians) to east (Gulf of Alaska) that, in part, can be attributed to a shift from off-shelf to continental shelf-based prey. Finally, stable isotope ratio results, supported by field observations, showed that the diet in spring and summer of eastern Aleutian Island transient killer whales is apparently not composed exclusively of Steller sea lions.

Use of chemical tracers in assessing the diet and foraging regions of eastern North Pacific killer whales

Top predators in the marine environment integrate chemical signals acquired from their prey that reflect both the species consumed and the regions from which the prey were taken. These chemical tracers—stable isotope ratios of carbon and nitrogen; persistent organic pollutant (POP) concentrations, patterns and ratios; and fatty acid profiles—were measured in blubber biopsy samples from North Pacific killer whales (Orcinus orca) (n = 84) and were used to provide further insight into their diet, particularly for the offshore group, about which little dietary information is available. The offshore killer whales were shown to consume prey species that were distinctly different from those of sympatric resident and transient killer whales. In addition, it was confirmed that the offshores forage as far south as California. Thus, these results provide evidence that the offshores belong to a third killer whale ecotype. Resident killer whale populations showed a gradient in stable isotope profiles from west (central Aleutians) to east (Gulf of Alaska) that, in part, can be attributed to a shift from off-shelf to continental shelf-based prey. Finally, stable isotope ratio results, supported by field observations, showed that the diet in spring and summer of eastern Aleutian Island transient killer whales is apparently not composed exclusively of Steller sea lions.

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 methods for plot-based sampling of wildlife populations

Classical sampling methods can be used to estimate the mean of a finite or infinite population. Block kriging also estimates the mean, but of an infinite population in a continuous spatial domain. In this paper, I consider a finite population version of block kriging (FPBK) for plot-based sampling. The data are assumed to come from a spatial stochastic process. Minimizing mean-squared-prediction errors yields best linear unbiased predictions that are a finite population version of block kriging. FPBK has versions comparable to simple random sampling and stratified sampling, and includes the general linear model. This method has been tested for several years for moose surveys in Alaska, and an example is given where results are compared to stratified random sampling. In general, assuming a spatial model gives three main advantages over classical sampling: (1) FPBK is usually more precise than simple or stratified random sampling, (2) FPBK allows small area estimation, and (3) FPBK allows nonrandom sampling designs.

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.

Review of Applied Spatial Data Analysis with R by R. S. Bivand, E. J. Pebesma, and V. Gomez-Rubio

Most authors struggle to pick a title that adequately conveys all of the material covered in a book. When I first saw Applied Spatial Data Analysis with R, I expected a review of spatial statistical models and their applications in packages (libraries) from the CRAN site of R. The authors’ title is not misleading, but I was very pleasantly surprised by how deep the word “applied” is here. The first half of the book essentially covers how R handles spatial data. To some statisticians this may be boring. Do you want, or need, to know the difference between S3 and S4 classes, how spatial objects in R are organized, and how various methods work on the spatial objects? A few years ago I would have said “no,” especially to the “want” part. Just let me slap my EXCEL spreadsheet into R and run some spatial functions on it. Unfortunately, the world is not so simple, and ultimately we want to minimize effort to get all of our spatial analyses accomplished. The first half of this book certainly convinced me that some extra effort in organizing my data into certain spatial class structures makes the analysis easier and less subject to mistakes. I also admit that I found it very interesting and I learned a lot.

Accounting for uncertainty in ecological analysis: the strengths and limitations of hierarchical statistical modeling

Analyses of ecological data should account for the uncertainty in the process(es) that generated the data. However, accounting for these uncertainties is a difficult task, since ecology is known for its complexity. Measurement and/or process errors are often the only sources of uncertainty modeled when addressing complex ecological problems, yet analyses should also account for uncertainty in sampling design, in model specification, in parameters governing the specified model, and in initial and boundary conditions. Only then can we be confident in the scientific inferences and forecasts made from an analysis. Probability and statistics provide a framework that accounts for multiple sources of uncertainty. Given the complexities of ecological studies, the hierarchical statistical model is an invaluable tool. This approach is not new in ecology, and there are many examples (both Bayesian and non-Bayesian) in the literature illustrating the benefits of this approach. In this article, we provide a baseline for concepts, notation, and methods, from which discussion on hierarchical statistical modeling in ecology can proceed. We have also planted some seeds for discussion and tried to show where the practical difficulties lie. Our thesis is that hierarchical statistical modeling is a powerful way of approaching ecological analysis in the presence of inevitable but quantifiable uncertainties, even if practical issues sometimes require pragmatic compromises.

A General Framework for the Analysis of Animal Resource Selection from Telemetry Data

We propose a general framework for the analysis of animal telemetry data through the use of weighted distributions. It is shown that several interpretations of resource selection functions arise when constructed from the ratio of a use and availability distribution. Through the proposed general framework, several popular resource selection models are shown to be special cases of the general model by making assumptions about animal movement and behavior. The weighted distribution framework is shown to be easily extended to readily account for telemetry data that are highly auto-correlated; as is typical with use of new technology such as global positioning systems animal relocations. An analysis of simulated data using several models constructed within the proposed framework is also presented to illustrate the possible gains from the flexible modeling framework. The proposed model is applied to a brown bear data set from southeast Alaska.

Results from the April 2009 Gulf of Alaska Line Transect Survey (GOALS) in the Navy Training Exercise Area

Little is known about the present-day occurrence of cetaceans found in offshore waters in the Gulf of Alaska; however, whaling records and a few recent surveys have shown this area to be important habitat. The U.S. Navy maintains a maritime training area in the central Gulf of Alaska, east of Kodiak Island, and has requested additional information on marine mammal presence and use of this area. To describe the occurrence and distribution of marine mammals in and around the U.S. Navy training area, a line transect visual and acoustic survey was conducted 10-20 April 2009 from the NOAA ship Oscar Dyson. The primary survey area encompassed nearshore and offshore pelagic waters of the central Gulf of Alaska. Survey lines were designed to provide equal coverage of the nearshore and offshore habitat.

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.

Regional differences in the spatial and temporal heterogeneity of oceanographic habitat used by Steller sea lions

Over the past three decades, the decline and altered spatial distribution of the western stock of Steller sea lions (Eumetopias jubatus) in Alaska have been attributed to changes in the distribution or abundance of their prey due to the cumulative effects of fisheries and environmental perturbations. During this period, dietary prey occurrence and diet diversity were related to population decline within metapopulation regions of the western stock of Steller sea lions, suggesting that environmental conditions may be variable among regions. The objective of this study, therefore, was to examine regional differences in the spatial and temporal heterogeneity of oceanographic habitat used by Steller sea lions within the context of recent measures of diet diversity and population trajectories. Habitat use was assessed by deploying satellite-depth recorders and satellite relay data loggers on juvenile Steller sea lions (n = 45) over a five-year period (2000–2004) within four regions of the western stock, including the western, central, and eastern Aleutian Islands, and central Gulf of Alaska. Areas used by sea lions during summer months (June, July, and August) were demarcated using satellite telemetry data and characterized by environmental variables (sea surface temperature [SST] and chlorophyll a [chl a]), which possibly serve as proxies for environmental processes or prey. Spatial patterns of SST diversity and Steller sea lion population trends among regions were fairly consistent with trends reported for diet studies, possibly indicating a link between environmental diversity, prey diversity, and distribution or abundance of Steller sea lions. Overall, maximum spatial heterogeneity coupled with minimal temporal variability of SST appeared to be beneficial for Steller sea lions. In contrast, these patterns were not consistent for chl a, and there appeared to be an ecological threshold. Understanding how Steller sea lions respond to measures of environmental heterogeneity will ultimately be useful for implementing ecosystem management approaches and developing additional conservation strategies.

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.

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.