Shifting elasmobranch community assemblage at Cocos Island—an isolated marine protected area

Fishing pressure has increased the extinction risk of many elasmobranch (shark and ray) species. Although many countries have established no‐take marine reserves, a paucity of monitoring data means it is still unclear if reserves are effectively protecting these species. We examined data collected by a small group of divers over the past 21 years at one of the world's oldest marine protected areas (MPAs), Cocos Island National Park, Costa Rica. We used mixed effects models to determine trends in relative abundance, or probability of occurrence, of 12 monitored elasmobranch species while accounting for variation among observers and from abiotic factors. Eight of 12 species declined significantly over the past 2 decades. We documented decreases in relative abundance for 6 species, including the iconic scalloped hammerhead shark (Sphyrna lewini) (?45%), whitetip reef shark (Triaenodon obesus) (?77%), mobula ray (Mobula spp.) (?78%), and manta ray (Manta birostris) (?89%), and decreases in the probability of occurrence for 2 other species. Several of these species have small home ranges and should be better protected by an MPA, which underscores the notion that declines of marine megafauna will continue unabated in MPAs unless there is adequate enforcement effort to control fishing. In addition, probability of occurrence at Cocos Island of tiger (Galeocerdo cuvier), Galapagos (Carcharhinus galapagensis), blacktip (Carcharhinus limbatus), and whale (Rhincodon typus) sharks increased significantly. The effectiveness of MPAs cannot be evaluated by examining single species because population responses can vary depending on life history traits and vulnerability to fishing pressure.     

Total mercury content found in edible tissues of top predator fish from the Gulf of California,Mexico

Fish at the top of the food chain bioaccumulate and biomagnify toxic metals including mercury (Hg), at a higher rate than nonpredatory fish. However, although some top predator fish species are important in the human diet, the risks for consumers in Mexico are difficult to evaluate due to the scarce baseline information available. In the present article, data on the total mercury (THg) concentrations found in edible tissues of different species of sharks, rays, large pelagic fishes and groupers from the Gulf of California are presented and compared with national and international health guidelines. During 2003 and 2004, 73 shark dorsal tissue samples, 52 ray samples, 66 large pelagic fish samples and 16 grouper samples were collected at different sites along the Gulf of California. Samples were digested in a microwave system and analyzed by cold vapor atomic absorption spectrophotometry. Thirty samples (14%) showed concentrations of Hg above the national statutory limit of 1?mg?kg^?1, from which 25 samples were sharks and 3 were large pelagic fishes. A positive correlation was found between the natural log of THg [ln(THg)] and total length (TL cm) in species of sharks, rays and groupers. The Teacapán estuary in Sinaloa was detected as a possible source of Hg into the Gulf of California due to historic mining in the watershed. Mean Hg concentrations exceeded the regulatory limit of 1?mg?kg^?1 in the following highly commercial shark species: smooth hammerhead, pelagic thrasher, pacific sharpnose shark, dusky shark, scalloped hammerhead, and whitenose shark, in that order. It is recommended that preventive food advisories be issued with respect to these species until further investigations are conducted.     

Cryptic hammerhead shark lineage occurrence in the western South Atlantic revealed by DNA analysis

A cryptic lineage of hammerhead shark closely related to but evolutionarily distinct from the scalloped hammerhead (Sphyrna lewini) was recently documented in the western North Atlantic Ocean. Here, we demonstrate using nuclear and mitochondrial DNA sequences that this cryptic lineage also occurs in the western South Atlantic Ocean, extending its distribution >7,000 km from its only previously reported location. Our results also further validate the existence of this evolutionarily distinct hammerhead shark lineage. The southern hemisphere cryptic individuals were 1.6 and 5.8% divergent from S. lewini (sensu stricto) for the nuclear internal transcribed spacer 2 (ITS2) and mitochondrial control region loci, respectively, and formed a strongly supported, reciprocally monophyletic sister group to sympatric S. lewini. Coalescent analysis (ITS2 locus) yielded a divergence estimate of ~4.5 million years between S. lewini and the cryptic lineage. Given expanding concerns about overfishing of the large-bodied hammerhead sharks, this cryptic lineage needs to be formally recognized and incorporated into shark management and conservation planning to avoid the inadvertent, potential extirpation of a unique hammerhead lineage.     

Genetic evidence of cryptic speciation within hammerhead sharks (Genus Sphyrna)

Surveys of genetic variation within cosmopolitan marine species often uncover deep divergences, indicating historical separation and potentially cryptic speciation. Based on broad geographic (coastal eastern North America, Gulf of Mexico, western Africa, Australia, and Hawaii) and temporal sampling (1991–2003), mitochondrial (control region [CR] and cytochrome oxidase I [COI]) and nuclear gene (lactate dehydrogenase A intron 6 [LDHA6]) variation among 76 individuals was used to test for cryptic speciation in the scalloped hammerhead, Sphyrna lewini (Griffith and Smith). CR and COI gene trees confirmed previous evidence of divergence between Atlantic and Indo-Pacific scalloped hammerhead populations; populations were reciprocally monophyletic. However, the between-basin divergence recorded in the mtDNA genome was not reflected in nuclear gene phylogenies; alleles for LDHA6 were shared between ocean basins, and Atlantic and Indo-Pacific populations were not reciprocally monophyletic. Unexpectedly, CR, COI, and LDHA6 gene trees recovered a deep phylogenetic partition within the Atlantic samples. For mtDNA haplotypes, which segregated by basin, average genetic distances were higher among Atlantic haplotypes (CR: D _HKY=0.036, COI: D _GTR=0.016) than among Indo-Pacific haplotypes (CR: D _HKY=0.010, COI: D _GTR=0.006) and approximated divergences between basins for CR (D _HKY=0.036 within Atlantic; D _HKY=0.042 between basins). Vertebral counts for eight specimens representing divergent lineages from the western north Atlantic were consistent with the genetic data. Coexistence of discrete lineages in the Atlantic, complete disequilibrium between nuclear and mitochondrial alleles within lineages and concordant partitions in genetic and morphological characters indicates reproductive isolation and thus the occurrence of a cryptic species of scalloped hammerhead in the western north Atlantic. Effective management of large coastal shark species should incorporate this important discovery and the inference from sampling that the cryptic scalloped hammerhead is less abundant than S. lewini, making it potentially more susceptible to fishery pressure.     

Tiger shark (Galeocerdo cuvier) abundance and growth in a subtropical embayment: evidence from 7 years of standardized fishing effort

The tiger shark (Galeocerdo cuvier Peron and Lesueur 1822) is a widely distributed predator with a broad diet and the potential to affect marine community structure, yet information on local patterns of abundance for this species is lacking. Tiger shark catch data were gathered over 7 years of tag and release research fishing (1991–2000, 2002–2004) in Shark Bay, Western Australia (25°45′S, 113°44′E). Sharks were caught using drumlines deployed in six permanent zones (~3 km² in area). Fishing effort was standardized across days and months, and catch rates on hooks were expressed as the number of sharks caught h⁻¹. A total of 449 individual tiger sharks was captured; 29 were recaptured. Tiger shark catch rate showed seasonal periodicity, being higher during the warm season (Sep–May) than during the cold season (Jun–Aug), and was marked by inter-annual variability. The most striking feature of the catch data was a consistent pattern of slow, continuous variation within each year from a peak during the height of the warm season (February) to a trough in the cold season (July). Annual growth rates of recaptured individuals were generally consistent with estimates from other regions, but exceeded those for populations elsewhere for sharks >275 cm fork length (FL), perhaps because mature sharks in the study area rely heavily on large prey. The data suggest that (1) the threat of predation faced by animals consumed by tiger sharks fluctuates dramatically within and between years, and (2) efforts to monitor large shark abundance should be extensive enough to detect inter-annual variation and sufficiently intensive to account for intra-annual trends.     

Negligible evidence for regional genetic population structure for two shark species <Emphasis Type="Italic">Rhizoprionodon acutus</Emphasis> (Rüppell, 1837) and <Emphasis Type="Italic">Sphyrna lewini</Emphasis> (Griffith &; Smith, 1834) with contrasting biology

Biodiversity of sharks in the tropical Indo-Pacific is high, but species-specific information to assist sustainable resource exploitation is scarce. The null hypothesis of population genetic homogeneity was tested for scalloped hammerhead shark (Sphyrna lewini, n = 237) and the milk shark (Rhizoprionodon acutus, n = 207) from northern and eastern Australia, using nuclear (S. lewini, eight microsatellite loci; R. acutus, six loci) and mitochondrial gene markers (873 base pairs of NADH dehydrogenase subunit 4). We were unable to reject genetic homogeneity for S. lewini, which was as expected based on previous studies of this species. Less expected were similar results for R. acutus, which is more benthic and less vagile than S. lewini. These features are probably driving the genetic break found between Australian and central Indonesian R. acutus (F-statistics; mtDNA, 0.751–0.903, respectively; microsatellite loci, 0.038–0.047 respectively). Our results support the spatially homogeneous monitoring and management plan for shark species in Queensland, Australia.     

Diel movement patterns of the scalloped hammerhead shark (Sphyrna lewini) in relation to El Bajo Espiritu Santo: a refuging central-position social system

Summary Movement patterns of scalloped hammerhead sharks in the vicinity of El Bajo Espiritu Santo, a seamount in the Gulf of California, were determined by tracking by ultrasonic telemetry 13 sharks and marking 100 sharks. The 13 tracked sharks swam back and forth along the seamount ridge throughout the day. They did not swim in different directions to reduce swimming effort when currents changed from a parallel to a perpendicular orientation to the ridge. Sharks tracked up to 8 km away into the pelagic environment soon returned to the seamount. From such trackings and repeated observations of marked sharks over periods of several weeks, it is believed that most sharks disperse and return to the seamount in a rhythmical fashion. The separate departures of individual hammerheads in five paired trackings indicated that the sharks left the seamount either in small groups or singly. For these reasons, we argue that the social system of the scalloped hammerhead shark can be described as a refuging system.     

A robust-design formulation of the incidence function model of metapopulation dynamics applied to two species of rails

The incidence function model (IFM) uses area and connectivity to predict metapopulation dynamics. However, false absences and missing data can lead to underestimates of the number of sites contributing to connectivity, resulting in overestimates of dispersal ability and turnovers (extinctions plus colonizations). We extend estimation methods for the IFM by using a hierarchical Bayesian model to account both for false absences due to imperfect detection and for missing data due to sites not surveyed in some years. We compare parameter estimates, measures of metapopulation dynamics, and forecasts using stochastic patch occupancy models (SPOMs) among three IFM models: (1) a Bayesian formulation assuming no false absences and omitting site-year combinations with missing data; (2) a hierarchical Bayesian formulation assuming no false absences but incorporating missing data; and (3) a hierarchical Bayesian formulation allowing for imperfect detection and incorporating missing data. We fit the models to multiyear data sets of occupancy for two bird species that differ in body size and presumed dispersal ability but inhabit the same network of sites: the small Black Rail (Laterallus jamaicensis) and the medium-sized Virginia Rail (Rallus limicola). Incorporating missing data affected colonization parameters and led to lower estimates of dispersal ability for the Black Rail. Detection rates were high for the Black Rail in most years but moderate for the Virginia Rail. Incorporating imperfect detection resulted in higher occupancy and lower turnover rates for both species, with largest effects for the Virginia Rail. Forecasts using SPOMs were sensitive to both missing data and false absences; persistence in models assuming no false absences was more optimistic than from robust models. Our results suggest that incorporating false absences and missing data into the IFM can improve (1) estimates of dispersal ability and the effect of connectivity on colonization, (2) the scaling of extinction risk with patch area, and (3) forecasts of occupancy and turnover rates.     

Bayesian Small Area Models for Assessing Wildlife Conservation Risk in Patchy Populations

Abstract:  Species conservation risk assessments require accurate, probabilistic, and biologically meaningful maps of population distribution. In patchy populations, the reasons for discontinuities are not often well understood. We tested a novel approach to habitat modeling in which methods of small area estimation were used within a hierarchical Bayesian framework. Amphibian occurrence was modeled with logistic regression that included third-order drainages as hierarchical effects to account for patchy populations. Models including the random drainage effects adequately represented species occurrences in patchy populations of 4 amphibian species in the Oregon Coast Range (U.S.A.). Amphibian surveys from other locations within the same drainage were used to calibrate local drainage-scale effects. Cross-validation showed that prediction errors for calibrated models were 77% to 86% lower than comparable regionally constructed models, depending on species. When calibration data were unavailable, small area and regional models performed similarly, although poorly. Small area estimation models complement wildlife ecology and habitat studies, and can help managers develop a regional picture of the conservation status for relatively rare species.     

Flexible modeling for stock-recruitment relationships using Bayesian nonparametric mixtures

The stock and recruitment relationship is fundamental to the management of fishery natural resources. However, inferring stock-recruitment relationships is a challenging problem because of the limited available data, the collection of plausible models, and the biological characteristics that should be reflected in the model. Motivated by limitations of traditional parametric stock-recruitment models, we propose a Bayesian nonparametric approach based on a mixture model for the joint distribution of log-reproductive success and stock biomass. Flexible mixture modeling for this bivariate distribution yields rich inference for the stock-recruitment relationship through the implied conditional distribution of log-reproductive success given stock biomass. The method is illustrated with cod data from six regions of the North Atlantic, including comparison with simpler Bayesian parametric and semiparametric models.     

Lifting a veil on diversity: a Bayesian approach to fitting relative-abundance models.

Bayesian methods incorporate prior knowledge into a statistical analysis. This prior knowledge is usually restricted to assumptions regarding the form of probability distributions of the parameters of interest, leaving their values to be determined mainly through the data. Here we show how a Bayesian approach can be applied to the problem of drawing inference regarding species abundance distributions and comparing diversity indices between sites. The classic log series and the lognormal models of relative- abundance distribution are apparently quite different in form. The first is a sampling distribution while the other is a model of abundance of the underlying population. Bayesian methods help unite these two models in a common framework. Markov chain Monte Carlo simulation can be used to fit both distributions as small hierarchical models with shared common assumptions. Sampling error can be assumed to follow a Poisson distribution. Species not found in a sample, but suspected to be present in the region or community of interest, can be given zero abundance. This not only simplifies the process of model fitting, but also provides a convenient way of calculating confidence intervals for diversity indices. The method is especially useful when a comparison of species diversity between sites with different sample sizes is the key motivation behind the research. We illustrate the potential of the approach using data on fruit-feeding butterflies in southern Mexico. We conclude that, once all assumptions have been made transparent, a single data set may provide support for the belief that diversity is negatively affected by anthropogenic forest disturbance. Bayesian methods help to apply theory regarding the distribution of abundance in ecological communities to applied conservation.     

Assessing conservation status with extensive but low-resolution data: Application of frequentist and Bayesian models to endangered Athabasca River rainbow trout

Use of extensive but low-resolution abundance data is common in the assessment of species at-risk status based on quantitative decline criteria under International Union for Conservation of Nature (IUCN) and national endangered species legislation. Such data can be problematic for 3 reasons. First, statistical power to reject the null hypothesis of no change is often low because of small sample size and high sampling uncertainty leading to a high frequency of type II errors. Second, range-wide assessments composed of multiple site-specific observations do not effectively weight site-specific trends into global trends. Third, uncertainty in site-specific temporal trends and relative abundance are not propagated at the appropriate spatial scale. A common result is the propensity to underestimate the magnitude of declines and therefore fail to identify the appropriate at-risk status for a species. We used 3 statistical approaches, from simple to more complex, to estimate temporal decline rates for a designatable unit (DU) of rainbow trout in the Athabasca River watershed in western Canada. This DU is considered a native species for purposes of listing because of its genetic composition characterized as >0.95 indigenous origin in the face of continuing introgressive hybridization with introduced populations in the watershed. Analysis of abundance trends from 57 time series with a fixed-effects model identified 33 sites with negative trends, but only 2 were statistically significant. By contrast, a hierarchical linear mixed model weighted by site-specific abundance provided a DU-wide decline estimate of 16.4% per year and a 3-generation decline of 93.2%. A hierarchical Bayesian mixed model yielded a similar 3-generation decline trend of 91.3% and the posterior distribution showed that the estimate had a >99% probability of exceeding thresholds for an endangered listing. We conclude that the Bayesian approach was the most useful because it provided a probabilistic statement of threshold exceedance in support of an at-risk status recommendation.     

Predicting [corrected] extinction risk in spite of predator-prey oscillations.

Most population viability analyses (PVA) assume that the effects of species interactions are subsumed by population-level parameters. We examine how robust five commonly used PVA models are to violations of this assumption. We develop a stochastic, stage-structured predator-prey model and simulate prey population vital rates and abundance. We then use simulated data to parameterize and estimate risk for three demographic models (static projection matrix, stochastic projection matrix, stochastic vital rate matrix) and two time series models (diffusion approximation [DA], corrupted diffusion approximation [CDA]). Model bias is measured as the absolute deviation between estimated and observed quasi-extinction risk. Our results highlight three generalities about the application of single-species models to multi-species conservation problems. First, our collective model results suggest that most single-species PVA models overestimate extinction risk when species interactions cause periodic variation in abundance. Second, the DA model produces the most (conservatively) biased risk forecasts. Finally, the CDA model is the most robust PVA to population cycles caused by species interactions. CDA models produce virtually unbiased and relatively precise risk estimates even when populations cycle strongly. High performance of simple time series models like the CDA owes to their ability to effectively partition stochastic and deterministic sources of variation in population abundance.     

A Bayesian hierarchical framework for calibrating aquatic biogeochemical models

Model practitioners increasingly place emphasis on rigorous quantitative error analysis in aquatic biogeochemical models and the existing initiatives range from the development of alternative metrics for goodness of fit, to data assimilation into operational models, to parameter estimation techniques. However, the treatment of error in many of these efforts is arguably selective and/or ad hoc. A Bayesian hierarchical framework enables the development of robust probabilistic analysis of error and uncertainty in model predictions by explicitly accommodating measurement error, parameter uncertainty, and model structure imperfection. This paper presents a Bayesian hierarchical formulation for simultaneously calibrating aquatic biogeochemical models at multiple systems (or sites of the same system) with differences in their trophic conditions, prior precisions of model parameters, available information, measurement error or inter-annual variability. Our statistical formulation also explicitly considers the uncertainty in model inputs (model parameters, initial conditions), the analytical/sampling error associated with the field data, and the discrepancy between model structure and the natural system dynamics (e.g., missing key ecological processes, erroneous formulations, misspecified forcing functions). The comparison between observations and posterior predictive monthly distributions indicates that the plankton models calibrated under the Bayesian hierarchical scheme provided accurate system representations for all the scenarios examined. Our results also suggest that the Bayesian hierarchical approach allows overcoming problems of insufficient local data by “borrowing strength” from well-studied sites and this feature will be highly relevant to conservation practices of regions with a high number of freshwater resources for which complete data could never be practically collected. Finally, we discuss the prospect of extending this framework to spatially explicit biogeochemical models (e.g., more effectively connect inshore with offshore areas) along with the benefits for environmental management, such as the optimization of the sampling design of monitoring programs and the alignment with the policy practice of adaptive management.     

Catch equations: restoring the missing terms in the nominally generalized Baranov catch equation

Observational models for the catch of fish at age a (or size) at time t are fundamental equations in fisheries science, linking a population model with data. The well known Baranov catch equation (which assumes that fishing and natural mortalities are constant over both age and time) is a nominal basis of those most commonly used in fish stock assessment and fish population models (which assume that fishing and natural mortalities vary with both age and time). But, what should a catch equation look like, if the instantaneous rates of fishing and natural mortalities of fish of age a at time t vary with age a and time t? Without answering this question, use of those catch equations in fish stock assessment and population models renders their results uncertain. In this paper, I derive a general catch in number or in biomass equation as observational models of an age- and time-dependent model for a fish population by Taylor series expansion of, and by directly manipulating, a general catch integral, reduce it to commonly used catch equations, and compare the performance of 11 of them using data on the western king prawn Penaeus latisulcatus. I show that the nominal generalization of the Baranov catch equation misses several terms. In so doing, I derive the catch equations more accurately and restore these missing terms. Although almost all approximations overestimate the catch per recruit for older prawns, all commonly used catch equations and their extensions perform worse than theoretically sound representations of the general catch equation and their approximations. The age-specific bias of all models is <2.5, <18 and <90% for a time interval of sampling of 1, 7 and 30 days, respectively. Such large biases even for moderate values of the length of the time interval of sampling highlight a need for assessing the utility of commonly used catch equations for individual species.     

Hotspots within hotspots? Hammerhead shark movements around Wolf Island, Galapagos Marine Reserve

Are pelagic species such as sharks and tuna distributed homogenously or heterogeneously in the oceans? Large assemblages of these species have been observed at seamounts and offshore islands in the eastern tropical Pacific, which are considered hotspots of pelagic biodiversity. Is the species distribution uniform at these hotspots or do species aggregate at a finer spatial scale at these sites? We employed three techniques to demonstrate that the aggregations of scalloped hammerhead sharks, Sphyrna lewini, and other pelagic species were confined to the southeastern corner of Wolf Island in the Galapagos Marine Reserve. Coded ultrasonic transmitters were placed on individuals at this site and at another aggregation site at Darwin Island, separated from Wolf by 40 km, and they were detected by monitors moored at the southeastern corner of Wolf Island and rarely by monitors deployed at other sites around the island. Hammerhead sharks, carrying depth-sensing continual transmitters, were tracked for two-day periods in a vessel and shown to reside a disproportionately large fraction of their time at the southeastern corner. Visual censuses were carried out seasonally at the eight monitor sites at Wolf Island, recording the abundance of one species of tuna, four species of jacks, and a number of other species. The highest diversity and abundance of these species occurred in the southeastern corner of the island. Our results support the use of hammerhead sharks as indicator and umbrella species for pelagic hotspots on a fine scale.     

Estimation of spatiotemporal trends in bat abundance from mortality data collected at wind turbines

Renewable energy sources, such as wind energy, are essential tools for reducing the causes of climate change, but wind turbines can pose a collision risk for bats. To date, the population-level effects of wind-related mortality have been estimated for only 1 bat species. To estimate temporal trends in bat abundance, we considered wind turbines as opportunistic sampling tools for flying bats (analogous to fishing nets), where catch per unit effort (carcass abundance per monitored turbine) is a proxy for aerial abundance of bats, after accounting for seasonal variation in activity. We used a large, standardized data set of records of bat carcasses from 594 turbines in southern Ontario, Canada, and corrected these data to account for surveyor efficiency and scavenger removal. We used Bayesian hierarchical models to estimate temporal trends in aerial abundance of bats and to explore the effect of spatial factors, including landscape features associated with bat habitat (e.g., wetlands, croplands, and forested lands), on the number of mortalities for each species. The models showed a rapid decline in the abundance of 4 species in our study area; declines in capture of carcasses over 7 years ranged from 65% (big brown bat [Eptesicus fuscus]) to 91% (silver-haired bat [Lasionycteris noctivagans]). Estimated declines were independent of the effects of mitigation (increasing wind speed at which turbines begin to generate electricity from 3.5 to 5.5 m/s), which significantly reduced but did not eliminate bat mortality. Late-summer mortality of hoary (Lasiurus cinereus), eastern red (Lasiurus borealis), and silver-haired bats was predicted by woodlot cover, and mortality of big brown bats decreased with increasing elevation. These landscape predictors of bat mortality can inform the siting of future wind energy operations. Our most important result is the apparent decline in abundance of four common species of bat in the airspace, which requires further investigation.     

Highly directional swimming by scalloped hammerhead sharks,Sphyrna lewini,and subsurface irradiance,temperature, bathymetry,and geomagnetic field

The homing behavior of scalloped hammerhead sharks (Sphyrna lewini) to and fro between Espiritu Santo Seamount and Las Animas Island and the surrounding pelagic environment was studied to reveal their mechanisms of navigation in the oceanic environment. Four sharks were tagged with ultrasonic transmitters and tracked at the former location and one shark at the latter site during July, August, or September between 1981 and 1988. Hammerhead swimming movements were highly oriented: the mean coefficient of concentration (r) for sets of ten consecutive swimming directions recorded during eight homing movements by three hammerhead sharks ranged from 0.885 to 0.996. Drift within a current could not explain this directionality, since highly variable directions were recorded from a transmitter floating at the sea surface after becoming detached from a shark. Forward swimming momentum was an unlikely explanation, since highly directional swimming was maintained for a period of 32 min with only a gradual change in course. To maintain directionality over this period, an environmental property should be necessary for guidance. The hammerheads swam at night, with repeated vertical excursions ranging from 100 to 450 m deep, out of view of either the sea surface or the sea floor. The sharks' vertical diving movements were compared to distributions of spectral irradiance (relative to elasmobranch scotopic and photopic visual sensitivities), temperature, and current-flow directions in the water column. No relationships were evident between these properties and the sharks' oriented swimming movements. Movements of scalloped hammerhead sharks to and from a seamount were compared to topographic features in bathymetry and geomagnetic field leading away from the seamount. Sharks swam repeatedly over fixed geographic paths, and these paths occurred less often along submarine ridges and valleys than maxima and minima in the geomagnetic field. No significant difference existed between the degree of association of points from the sharks' tracks and points from track simulations and 20° changes in the slope of the depth record. On the other hand, significantly more points from the sharks' tracks were associated with slope changes in the magnetic intensity record than points from track simulations. A magnetic intensity gradient of 0.037 nanoteslas/m (nT/m) existed at 175 m depth, where a shark swam directionally, and this gradient was three times steeper than that measured at the sea surface and exceeded that recorded at a depth of 200 m. The hammerheads are hypothesized to find the seamount using geomagnetic topotaxis. The shark could be attracted to and move back and forth along ridges and valleys, features in the relief of magnetic field intensities occurring over a geographical area.     

Trends in the exploitation of South Atlantic shark populations

Approximately 25% of globally reported shark catches occur in Atlantic pelagic longline fisheries. Strong declines in shark populations have been detected in the North Atlantic, whereas in the South Atlantic the situation is less clear, although fishing effort has been increasing in this region since the late 1970s. We synthesized information on shark catch rates (based on 871,177 sharks caught on 86,492 longline sets) for the major species caught by multiple fleets in the South Atlantic between 1979 and 2011. We complied records from fishing logbooks of fishing companies, fishers, and onboard observers that were supplied to Brazilian institutions. By using exploratory data analysis and literature sources, we identified 3 phases of exploitation in these data (Supporting Information). From 1979 to 1997 (phase A), 5 fleets (40 vessels) fished mainly for tunas. From 1998 to 2008 (phase B), 20 fleets (100 vessels) fished for tunas, swordfishes, and sharks. From 2008 to 2011 (phase C), 3 fleets (30 vessels) fished for multiple species, but restrictive measures were implemented. We used generalized linear models to standardize catch rates and identify trends in each of these phases. Shark catch rates increased from 1979 to 1997, when fishing effort was low, decreased from 1998 to 2008, when fishing effort increased substantially, and remained stable or increased from 2008 to 2011, when fishing effort was again low. Our results indicate that most shark populations affected by longlines in the South Atlantic are currently depleted, but these populations may recover if fishing effort is reduced accordingly. In this context, it is problematic that comprehensive data collection, monitoring, and management of these fisheries ceased after 2012. Concurrently with the fact that Brazil is newly identified by FAO among the largest (and in fastest expansion) shark sub‐products consumer market worldwide.     

Defining trade-offs among conservation, profitability, and food security in the California current bottom-trawl fishery

Although it is recognized that marine wild-capture fisheries are an important source of food for much of the world, the cost of sustainable capture fisheries to species diversity is uncertain, and it is often questioned whether industrial fisheries can be managed sustainably. We evaluated the trade-off among sustainable food production, profitability, and conservation objectives in the groundfish bottom-trawl fishery off the U.S. West Coast, where depletion (i.e., reduction in abundance) of six rockfish species (Sebastes) is of particular concern. Trade-offs are inherent in this multispecies fishery because there is limited capacity to target species individually. From population models and catch of 34 stocks of bottom fish, we calculated the relation between harvest rate, long-term yield (i.e., total weight of fish caught), profit, and depletion of each species. In our models, annual ecosystem-wide yield from all 34 stocks was maximized with an overall 5.4% harvest rate, but profit was maximized at a 2.8% harvest rate. When we reduced harvest rates to the level (2.2% harvest rate) at which no stocks collapsed (<10% of unfished levels), biomass harvested was 76% of the maximum sustainable yield and profit 89% of maximum. A harvest rate under which no stocks fell below the biomass that produced maximum sustainable yield (1% harvest rate), resulted in 45% of potential yield and 67% of potential profit. Major reductions in catch in the late 1990s led to increase in the biomass of the most depleted stocks, but this rebuilding resulted in the loss of >30% of total sustainable yield, whereas yield lost from stock depletion was 3% of total sustainable yield. There are clear conservation benefits to lower harvest rates, but avoiding overfishing of all stocks in a multispecies fishery carries a substantial cost in terms of lost yield and profit.