Satellite tracking reveals a dichotomy in migration strategies among juvenile loggerhead turtles in the Northwest Atlantic

Few data are available on the movements and behavior of immature Atlantic loggerhead sea turtles (Caretta caretta) from their seasonal neritic foraging grounds within the western north Atlantic. These waters provide developmental habitat for loggerheads originating from several western Atlantic nesting stocks. We examined the long-term movements of 23 immature loggerheads (16 wild-caught and seven headstart turtles) characterizing their seasonal distribution, habitat use, site fidelity, and the oceanographic conditions encountered during their migrations. We identified two movement strategies: (1) a seasonal shelf-constrained north–south migratory pattern; and (2) a year-round oceanic dispersal strategy where turtles travel in the Gulf Stream to the North Atlantic and their northern dispersal is limited by the 10–15°C isotherm. When sea surface temperatures dropped below 20°C, neritic turtles began a migration south of Cape Hatteras, North Carolina (USA) where they established fidelity to the waters between North Carolina’s Outer Banks and the western edge of the Gulf Stream along outer continental shelf. Two turtles traveled as far south as Florida. Several turtles returned to their seasonal foraging grounds during subsequent summers. Northern movements were associated with both increased sea surface temperature (>21°C) and increased primary productivity. Our results indicate strong seasonal and interannual philopatry to the waters of Virginia (summer foraging habitat) and North Carolina (winter habitat). We suggest that the waters of Virginia and North Carolina provide important seasonal habitat and serve as a seasonal migratory pathway for immature loggerhead sea turtles. North Carolina’s Cape Hatteras acts as a seasonal “migratory bottleneck” for this species; special management consideration should be given to this region. Six turtles spent time farther from the continental shelf. Three entered the Gulf Stream near Cape Hatteras, traveling in the current to the northwest Atlantic. Two of these turtles remained within an oceanic habitat from 1 to 3 years and were associated with mesoscale features and frontal systems. The ability of large benthic subadults to resume an oceanic lifestyle for extended periods indicates plasticity in habitat use and migratory strategies. Therefore, traditional life history models for loggerhead sea turtles should be reevaluated.     

Multiple trans-Arctic passages in the red alga Phycodrys rubens: evidence from nuclear rDNA ITS sequences

In order to investigate how episodes of geological and climatic change have influenced the distribution and evolutionary diversification of Arctic to cold temperate-North Atlantic seaweed species, intraspecific genetic variation was analyzed among isolates of the sublittoral, benthic red alga Phycodrys rubens (collected between June 1992 and January 1994). Rooted phylogenetic analyses of nuclear ribosomal DNA internal transcribed spacer (ITS) sequences and the plastid encoded Rubisco spacer sequences suggest that P. rubens invaded the North Atlantic from the Pacific shortly after the opening of the Bering Strait (3 to 3.5 million years ago), colonizing both the western and eastern Atlantic coasts. Based on these data we further hypothesize that P. rubens survived along the European coasts during the more recent Pleistocene glaciations, while becoming locally extinct along the North American Atlantic coasts. Following retraction of the last ice sheet, the western Atlantic coast was colonized a second time from the Pacific. The presence of two distinct genetic types (based on ITS and Rubisco sequences) along the European coasts is postulated to be a result of isolation and subsequent differentiation. This is likely because ice-free areas are known to have existed in northern Scotland and Norway during the last glaciation. The presence of an East Atlantic genetic type along the West Atlantic coast is believed to be a recent introduction (caused by human activity) of P. rubens to Newfoundland.     

Upper-slope distributions of oceanic lanternfishes (family: Myctophidae)

In order to determine the effects of bottom depth and water column structure on slope distributions, lanternfishes in the Cape Canyon and Cape Point Valley (eastern South Atlantic) and from the Feni Ridge and Hebridean Terrace (eastern North Atlantic) were sampled with commercial midwater nets fished 10 m above the bottom at preselected depths. Fifty-one species were taken in the eastern South Atlantic and 17 species in the eastern North Atlantic. Analyses of the samples included clustering and multi-dimensional scaling ordination of the root-root transformed densities using the Bray-Curtis similarity index. In the eastern South Atlantic (March 1988), where the horizontal and vertical temperature structuring of the water column was marked, oceanic myctophids were not taken shallower than the 300 m-depth horizon and less than one-third of the total number of species occurred inside the 500 m-depth horizon. Temperatures and bottom depths correlated with the downslope distributions of the species. Down-slope zonation was apparent. In the eastern North Atlantic (May 1983), where the temperature structure of the water column was similar at all stations, no oceanic species were taken at the 200 m-depth horizon and more than one-half of the total number of species was taken at the 300 m-depth horizon. Species distributions correlated with depth, and downslope zonation was not evident. The data sets, as well as information from the West Florida slope, suggest that both the depth and the regional temperature structure of the water column effect the upper slope distributions of oceanic myctophids. Long-slope, ribbon distribution patterns (generalized down-slope zonation) in myctophids may therefore be localized artifacts produced by the interaction of these factors.     

Amphi-Atlantic phylogeography of direct-developing lineages of Lasaea, a genus of brooding bivalves

Direct-developing intertidal Lasaea spp. occur in the North Atlantic as both continental margin and oceanic island populations. We conducted a molecular phylogenetic analysis of representative populations in order to test colonization hypotheses for North Atlantic oceanic islands. Thirty individuals each were collected in 1995 and 1996 from two continental putative source populations (Florida, Iberia) and two oceanic island populations (Bermuda, Azores). They were sequenced for a 462 nucleotide portion of the mitochondrial large ribosomal subunit (16S) gene. No amphi-Atlantic genotypes were detected: Bermudan lineages co-clustered exclusively with Floridian congeners, and Azorean samples formed an exclusive clade with Iberian haplotypes. Our data indicate that geographical proximity to continental source populations is a better predictor of phylogenetic relationships in North Atlantic Lasaea spp. than present-day oceanic surface circulation patterns. The phylogenetic trees generated are not consistent with colonization of oceanic islands by indirect-developing ancestral lineages or by truly trans-oceanic rafting events. However, they are consistent with predicted topologies resulting from limited (≤ 2000 km), long-distance colonization by rafting (against present-day circulation patterns in the case of the Azores) and from anthropogenic introductions. Received: 17 December 1998 / Accepted: 7 June 1999     

Global population divergence of the sea star Hippasteria phrygiana corresponds to the onset of the last glacial period of the Pleistocene

Genetic structure and connectivity of populations of the globally distributed and eurybathic sea star Hippasteria phrygiana (Parelius 1768) were studied in 165 individuals sampled from three oceanic regions: the North Pacific Ocean, the South Pacific Ocean (considered to include the adjacent regions of the Southern Ocean and the southern Indian Ocean) and the North Atlantic Ocean. A nuclear gene region (ATP synthase subunit α intron #5, ATPSα) and a mitochondrial gene region (cytochrome oxidase subunit I, COI) were amplified and sequenced. Significant heterogeneity was detected in an AMOVA analysis among the three sampled oceanic regions for COI, but not for ATPSα. Neither gene showed significant genetic heterogeneity within the North Atlantic, as assessed by ΦST values. Significant heterogeneity was detected for COI (but not ATPSα) in the North Pacific, but the converse was true in the South Pacific. Coalescent simulations suggested that the three regions have been diverging with little or no gene flow for the past 50–75,000 years, a time frame that corresponds to the onset of the last glacial period of the Pleistocene. A possible genetic signature of recent population expansion (or non-neutrality) was detected for each gene in the North Pacific, but not in the other two oceanic regions.     

Oceanographic characteristics of areas used by Cory’s shearwaters during short and long foraging trips in the North Atlantic

Many breeding seabird species are central-place foragers and constrained to find productive prey patches within their foraging ranges. We assessed how different populations of a pelagic seabird species, the Cory’s shearwater Calonectris diomedea, breeding in oceanic and neritic conditions, cope with these constraints in the North Atlantic, during both incubation and chick-rearing periods. We analysed 237 foraging trips to study the movements and oceanographic characteristics of foraging habitats of seven different populations of Cory’s shearwaters. Generally, oceanic populations exhibited higher foraging effort, by travelling more time and to more distant areas, and larger home ranges and feeding areas, than the neritic population (i.e. breeding on an island within the Portuguese continental platform). On their short trips (i.e. ≤4 days), birds from the different populations fed mostly in shallower waters around the colony. During long trips (i.e. ≥5 days), feeding areas of both oceanic and neritic populations were characterized by high concentration values of chlorophyll-a, low sea-surface temperature and shallower habitats, with oceanic populations of the Azores exploiting areas north of the islands over known seamounts and frontal regions. Birds from other oceanic population (Selvagens) also exploited the African continental shelf system on their long trips. The home ranges of the different populations overlapped widely, but there was a general spatial segregation in terms of the core feeding areas at the population level. Core feeding areas and areas of foraging overlap between different populations should be important to inform conservation management measures, such as the definition of Marine Important Bird Areas for seabirds over the North Atlantic.     

Identifying oceanic foraging grounds of sea turtles in the Atlantic using lead isotopes

Many species of marine organisms go through ontogenetic shifts that occur in unknown or inaccessible locations. Finding these areas is crucial to understand connectivity and resilience of populations, both of which have conservation implications. When extrinsic markers are not suitable to track organisms, intrinsic markers can be useful to infer the location of inaccessible areas where these cryptic stages occur. Our study focuses on the location of oceanic foraging areas of the cryptic early juvenile stage of green turtles, Chelonia mydas, in the Atlantic. Due to the small size of hatchlings, the use of telemetry is limited to short periods of time and small spatial ranges, which do not allow determining the location of oceanic foraging areas. We used lead (Pb) stable isotopes to determine the possible location of oceanic foraging areas of small green turtles in the Atlantic Ocean. Pb stable isotope ratios in the scute tissue deposited when turtles were in the oceanic habitat were compared to ratios of major sources of lead in the Atlantic and oceanic areas in the Atlantic to determine the location of oceanic foraging grounds. The Pb isotope ratios in the scute of oceanic-stage green turtles indicated that turtles use different regions in the Atlantic and that they are capable of transatlantic migrations. We compare the oceanic locations identified by this study with those suggested by two previous studies.     

Concentrations of V,Cr, Mn,Fe, Ni,Co, Cu,Zn, As and Cd in mesopelagic crustaceans from the North East Atlantic Ocean

Specimens of oceanic decapods, mysids and euphausiids, collected from the North East Atlantic Ocean during July 1985, were analysed for V, Cr, Mn, Fe, Ni, Co, Cu, Zn, As and Cd. Measurement of the metals was carried out using two techniques: inductively-coupled plasma-source mass spectrometry (ICP-MS) and atomic absorption (AAS). With the exception of the iron data, good agreement was observed between the data from the two techniques. Baseline-metal data are presented for the species measured, along with any effects of animal size on metal concentration. The data presented for oceanic animals from different trophic levels are important in studies of metal fluxes and for the assessment of markers within the food web in the oceans.     

Song of the humpback whale — Population comparisons

Summary Humpback whale songs recorded on tropical calving grounds exhibit different dialects depending on the oceanic basin. Songs sampled simultaneously from two populations in the North Pacific (Hawaii and Mexico) were essentially identical. These North Pacific songs were clearly different from the song type shared by two populations in the North Atlantic (Cape Verde Islands and West Indies). Songs from the Southern Hemisphere (Tonga) represent a third distinct dialect. Our evidence shows that, despite annual change in song organization, significant differences in humpback song occur between isolated ocean basins, while only subtle differences exist within an oceanic population (Hawaii and Pacific Mexico).     

Hierarchical modeling of count data with application to nuclear fall-out

We propose a method for a Bayesian hierarchical analysis of count data that are observed at irregular locations in a bounded domain of R². We model the data as having been observed on a fine regular lattice, where we do not have observations at all the sites. The counts are assumed to be independent Poisson random variables whose means are given by a log Gaussian process. In this article, the Gaussian process is assumed to be either a Markov random field (MRF) or a geostatistical model, and we compare the two models on an environmental data set. To make the comparison, we calibrate priors for the parameters in the geostatistical model to priors for the parameters in the MRF. The calibration is obtained empirically. The main goal is to predict the hidden Poisson-mean process at all sites on the lattice, given the spatially irregular count data; to do this we use an efficient MCMC. The spatial Bayesian methods are illustrated on radioactivity counts analyzed by Diggle et al. (1998).     

Magnetic orientation by hatchling loggerhead sea turtles (Caretta caretta) from the Gulf of Mexico

Loggerhead sea turtle hatchlings emerge from nests on either the east or west coast of the South Florida peninsula and then migrate offshore in opposite directions. Under laboratory conditions, magnetic cues induce east coast hatchlings to swim in directions that promote their transport by oceanic surface currents, such as the North Atlantic gyre. However, the surface currents used by west coast hatchlings are unknown. We examined the responses of west (Sarasota) hatchlings to magnetic cues in the Gulf of Mexico, the Florida Straits, and the Gulf Stream to determine their (1) likely migratory routes (2) orientation where currents lead into the Atlantic Ocean, and (3) orientation adjacent to Florida’s east coast. The results suggest that migration inside Gulf waters may be circuitous, that the turtles respond appropriately to enter Atlantic waters, and that orientation along Florida’s east coast probably promotes transport by the Gulf Stream into the North Atlantic gyre.     

Afro-Eurasia and the Americas present barriers to gene flow for the cosmopolitan neustonic nudibranch Glaucus atlanticus

Pelagic species have been traditionally thought to occupy vast, genetically interconnected, geographic ranges in an essentially homogeneous environment. Although this view has been challenged recently for some mesopelagic planktonic taxa, the population structure of hyponeustonic (surface-drifting) species remains unknown. Here, we test the hypothesis of panmixis in Glaucus atlanticus, a cosmopolitan neustonic nudibranch, by assessing the genetic differentiation of multiple representatives from a global neustonic sampling effort. Specimens were collected from all subtropical oceanic gyre systems (North Atlantic, South Atlantic, North Pacific, South Pacific, and Indian Ocean). We sequenced a fragment of the mitochondrial cytochrome oxidase I gene for 98 individuals and performed population structure, differentiation (analysis of molecular variance, spatial analysis of molecular variance, F _ST, Jost’s D), and molecular clock analyses. Our results indicate that G. atlanticus is not globally panmictic, but that populations appear to be panmictic within ocean basins. We detected several topologically ectopic haplotypes in the Atlantic Ocean, but the molecular clock analysis indicates that these have diverged from closely related Indo-Pacific haplotypes over 1.2 MYA, coinciding with cooling in waters around in the southern tip of Africa and resulting oceanographic changes. These data and the fact that G. atlanticus is not known from polar latitudes suggest that gene flow between ocean basins is hindered by physical barriers (supercontinents) and water temperatures in the Arctic and Southern Oceans.     

Changes in size and trends of North American sea duck populations associated with North Pacific oceanic regime shifts

Broad-scale multi-species declines in populations of North American sea ducks for unknown reasons is cause for management concern. Oceanic regime shifts have been associated with rapid changes in ecosystem structure of the North Pacific and Bering Sea. However, relatively little is known about potential effects of these changes in oceanic conditions on marine bird populations at broad scales. I examined changes in North American breeding populations of sea ducks from 1957 to 2011 in relation to potential oceanic regime shifts in the North Pacific in 1977, 1989, and 1998. There was strong support for population-level effects of regime shifts in 1977 and 1989, but little support for an effect of the 1998 shift. The continental-level effects of these regime shifts differed across species groups and time. Based on patterns of sea duck population dynamics associated with regime shifts, it is unclear if the mechanism of change relates to survival or reproduction. Results of this analysis support the hypothesis that population size and trends of North American sea ducks are strongly influenced by oceanic conditions. The perceived population declines appear to have halted >20 years ago, and populations have been relatively stable or increasing since that time. Given these results, we should reasonably expect dramatic changes in sea duck population status and trends with future oceanic regime shifts.     

Mitochondrial evolution and phylogeography in the hydrozoan<Emphasis Type="Italic"> Obelia geniculata</Emphasis> (Cnidaria)

The distribution and genetic structure of many marine invertebrates in the North Atlantic have been influenced by the Pleistocene glaciation, which caused local extinctions followed by recolonization in warmer periods. Mitochondrial DNA markers are typically used to reconstruct species histories. Here, two mitochondrial markers [16S rDNA and cytochrome c oxidase I (COI)] were used to study the evolution of the widely distributed hydrozoan Obelia geniculata (Linnaeus, 1758) from the North Atlantic and the Pacific and, more specifically, in the context of North Atlantic phylogeography. Samples were collected from six geographic localities between 1998 and 2002. Hydroids from the North Atlantic, North Pacific (Japan), and South Pacific (New Zealand) are reciprocally monophyletic and may represent cryptic species. Using portions of the 16S rDNA and COI genes and the date of the last trans-Arctic interchange (3.1–4.1 million years ago), the first calibrated rate of nucleotide substitutions in hydrozoans is presented. Whereas extremely low substitution rates have been reported in other cnidarians, mainly based on anthozoans, substitution rates in O. geniculata are comparable to other invertebrates. Despite a life history that ostensibly permits substantial dispersal, there is apparently considerable genetic differentiation in O. geniculata. Divergence estimates and the presence of unique haplotypes provide evidence for glacial refugia in Iceland and New Brunswick, Canada. A population in Massachusetts, USA, appears to represent a relatively recent colonization event.Communicated by J.P. Grassle, New Brunswick     

Gridded meteorological data as a resource for mechanistic macroecology in coastal environments

Gridded weather data were evaluated as sources of forcing variables for biophysical models of intertidal animal body temperature with model results obtained using local weather station data serving as the baseline of comparison. The objective of the study was to determine which gridded data are sufficient to capture observed patterns of thermal stress. Three coastal sites in western North America were included in this analysis: Boiler Bay, Oregon; Bodega Bay, California; and Pacific Grove, California. The gridded data with the highest spatial resolution, the 32-km North American Regional Reanalysis (NARR) and the 38-km Climate Forecasting System Reanalysis (CFSR), predicted daily maximum intertidal animal temperature most similarly to the local weather Station data. Time step size was important for variables that change rapidly throughout the day, such as solar radiation. There were site-based differences in the ability of the model to predict daily maximum intertidal animal temperature, with the gridded data predictions being the closest to local weather station predictions in Boiler Bay, Oregon. In a review of gridded data used as part of ecological studies, there was broad use of the data across subject areas and ecosystems so the recent improvements in the spatial (from 2 degrees to 32 km) and temporal scales (from 6 hours to 1 hour) of gridded data will further add to the applicability within the ecological community particularly for mechanistic studies.     

Closed population estimation models and their extensions in Program MARK

Program MARK provides > 65 data types in a common configuration for the estimation of population parameters from mark-encounter data. Encounter information from live captures, live resightings, and dead recoveries can be incorporated to estimate demographic parameters. Available estimates include survival (S or ϕ), rate of population change (λ), transition rates between strata (Ψ), emigration and immigration rates, and population size (N). Although N is the parameter most often desired by biologists, N is one of the most difficult parameters to estimate precisely without bias for a geographically and demographically closed population. The set of closed population estimation models available in Program MARK incorporate time (t) and behavioral (b) variation, and individual heterogeneity (h) in the estimation of capture and recapture probabilities in a likelihood framework. The full range of models from M ₀ (null model with all capture and recapture probabilities equal) to M _tbh are possible, including the ability to include temporal, group, and individual covariates to model capture and recapture probabilities. Both the full likelihood formulation of Otis et al. (1978) and the conditional model formulation of Huggins (1989, 1991) and Alho (1990) are provided in Program MARK, and all of these models are incorporated into the robust design (Kendall et al. 1995, 1997; Kendall and Nichols 1995) and robust-design multistrata (Hestbeck et al. 1991, Brownie et al. 1993) data types. Model selection is performed with AICc (Burnham and Anderson 2002) and model averaging (Burnham and Anderson 2002) is available in Program MARK to provide estimates of N with standard error that reflect model selection uncertainty.     

Particulate protein-nitrogen in North Atlantic surface waters

Depth profiles of particulate protein-nitrogen at 4 oceanic and 2 upwelling stations in the North Atlantic Ocean were measured by a new fluorometric method. The protein-nitrogen in the upper 20 m ranged from 0.19 to 1.61 μg-at N/1 at the oceanic stations and from 0.43 to 3.54 μg-at/1 at the upwelling stations. The mean values in the euphotic zone were 0.54 μg-at N/1 for the oceanic stations and 1.70 μg-at N/1 for the upwelling stations. The ratio of protein-nitrogen to chlorophyll at the two sets of stations was 2.83 and 0.54 μg-at N/μg chlorophyll, respectively. Regression analysis of the pooled data yielded a detritus and zooplankton-free ratio of 0.38 μg-at N:μg chlorophyll. Calculations of the phytoplankton protein-nitrogen, based on this ratio, suggest that in the oceanic water only 20% of the sestonic protein-nitrogen is associated with the phytoplankton. In the upwelling waters, the phytoplankton may account for 65% of the sestonic proteinnitrogen.     

Kernel-based home range method for data with irregular sampling intervals

Studies of habitat selection and movements often use radio-tracking data for defining animal home ranges. Home ranges (HR) can be approximated by a utilization density distribution (UD) that instead of assuming uniform use of areas within HR boundary provides a probabilistic measure of animal space use. In reality, radio-tracking data contain periods of frequent autocorrelated observations interspersed with temporally more independent observations. Using such temporally irregular data directly may result in biased UD estimates, because areas that have been sampled intensively receive too much weight. The problem of autocorrelation has been tackled by resampling data with an appropriate time interval. However, resampling may cause a large reduction in the data set size along with a loss of information. Evidently, biased UD estimates or reduction in data may prejudice the results on animal habitat selection and movement. We introduce a new method for estimating UDs with temporally irregular data. The proposed method, called the time kernel, accounts for temporal aggregation of observations and gives less weight to temporally autocorrelated observations. A further extension of the method accounts also for spatially aggregated observations with relatively low weights given to observations that are both temporally and spatially aggregated. We test the behaviour of the time kernel method and its spatiotemporal version using simulated data. In addition, the method is applied to a data set of brown bear locations.     

A spatially explicit estimate of the prewhaling abundance of the endangered North Atlantic right whale

The North Atlantic right whale (NARW) (Eubalaena glacialis) is one of the world's most threatened whales. It came close to extinction after nearly a millennium of exploitation and currently persists as a population of only approximately 500 individuals. Setting appropriate conservation targets for this species requires an understanding of its historical population size, as a baseline for measuring levels of depletion and progress toward recovery. This is made difficult by the scarcity of records over this species’ long whaling history. We sought to estimate the preexploitation population size of the North Atlantic right whale and understand how this species was distributed across its range. We used a spatially explicit data set on historical catches of North Pacific right whales (NPRWs) (Eubalaena japonica) to model the relationship between right whale relative density and the environment during the summer feeding season. Assuming the 2 right whale species select similar environments, we projected this model to the North Atlantic to predict how the relative abundance of NARWs varied across their range. We calibrated these relative abundances with estimates of the NPRW total prewhaling population size to obtain high and low estimates for the overall NARW population size prior to exploitation. The model predicted 9,075–21,328 right whales in the North Atlantic. The current NARW population is thus <6% of the historical North Atlantic carrying capacity and has enormous potential for recovery. According to the model, in June–September NARWs concentrated in 2 main feeding areas: east of the Grand Banks of Newfoundland and in the Norwegian Sea. These 2 areas may become important in the future as feeding grounds and may already be used more regularly by this endangered species than is thought.     

Upper temperature tolerances of some European molluscs. I Tellina fabula and T. tenuis

The upper temperature tolerances of two European species of Tellina from populations near the northern and southern limits of their geographical range are compared by means of median lethal temperature (LT₅₀), and median burial temperature (BT₅₀) determination for periods of exposure up to 96 h. The influence of previous acclimation temperature, and salinity within a limited range, on these determinations is also considered. T. fabula Gmelin, which is found in deeper water throughout its range, shows a lower thermal tolerance than T. tenuis da Costa. Both species display a significant effect of previous acclimation temperatures, the effect on the burrowing responses being greater than that shown by LT₅₀, and in both species the thermal tolerance of individuals from Mediterranean populations is greater than that of individuals from North Atlantic populations.