First Direct Evidence for Natal Wintering Ground Fidelity and Estimate of Juvenile Survival in the New Zealand Southern Right Whale Eubalaena australis

Abstract

Juvenile survival and recruitment can be more sensitive to environmental, ecological and anthropogenic factors than adult survival, influencing population-level processes like recruitment and growth rate in long-lived, iteroparous species such as southern right whales. Conventionally, Southern right whales are individually identified using callosity patterns, which do not stabilise until 6–12 months, by which time the whale has left its natal wintering grounds. Here we use DNA profiling of skin biopsy samples to identify individual Southern right whales from year of birth and document their return to the species’ primary wintering ground in New Zealand waters, the Subantarctic Auckland Islands. We find evidence of natal fidelity to the New Zealand wintering ground by the recapture of 15 of 57 whales, first sampled in year of birth and available for subsequent recapture, during winter surveys to the Auckland Islands in 1995–1998 and 2006–2009. Four individuals were recaptured at the ages of 9 to 11, including two females first sampled as calves in 1998 and subsequently resampled as cows with calves in 2007. Using these capture-recapture records of known-age individuals, we estimate changes in survival with age using Cormack-Jolly-Seber models. Survival is modelled using discrete age classes and as a continuous function of age. Using a bootstrap method to account for uncertainty in model selection and fitting, we provide the first direct estimate of juvenile survival for this population. Our analyses indicate a high annual apparent survival for juveniles at between 0.87 (standard error (SE) 0.17, to age 1) and 0.95 (SE 0.05: ages 2–8). Individual identification by DNA profiling is an effective method for long-term demographic and genetic monitoring, particularly in animals that change identifiable features as they develop or experience tag loss over time.

About the Authors:

E. L. Carroll

* E-mail: [email protected]

Affiliations School of Biological Sciences, University of Auckland, Auckland, New Zealand, Scottish Oceans Institute, University of St Andrews, East Sands, St Andrews, Scotland

R. M. Fewster

Affiliation: Department of Statistics, University of Auckland, Auckland, New Zealand

S. J. Childerhouse

Affiliation: Blue Planet Marine, Kingston, Hobart, Tasmania, Australia

N. J. Patenaude

Affiliation: Collégial International Sainte-Anne, Montréal, Québec, Canada

L. Boren

Affiliation: New Zealand Department of Conservation, Wellington, New Zealand

C. S. Baker

Affiliations School of Biological Sciences, University of Auckland, Auckland, New Zealand, Marine Mammal Institute and Department of Fisheries and Wildlife, Hatfield Marine Science Center, Oregon State University, Newport, United States of America

Introduction

Effective management of populations and species requires precise estimates of key demographic parameters, such as survival and rate of population change. While critical in the context of endangered and threatened species [1], such demographic parameters typically require long-term studies in which individuals are tracked through different life history stages. Such studies rely on individual identification using features that remain intact as the individual ages. Commonly used methods for identifying individuals include implanted or attached tags (e.g. PIT tagging [2]), genetic markers [3,4], assessment of phenotypes that differ between individuals such as acoustic calls [5,6] and photo-identification (photo-ID) of natural markings. Such natural markings can include part of the animal’s body, such as primate faces [7] and lion (Panthera leo) whisker patterns [8], or a composite of distinctive features, such as those used to identify Asian elephants (Elephas maximus [9]), tigers (Panthera tigris [10]), and marine species including various sharks and rays [11]. Computerised systems are increasingly used to recognise and match the phenotypes of individuals to create capture histories for use in demographic and ecological studies [12].

Many typical methods for identifying individuals may not be applicable to juveniles, particularly in large, long-lived species where the identifiable features of individuals can change with age. For example, tiger stripe patterns [10] and fluke patterns of humpback whale (Megaptera novaeangliae [13]) are not suitable for identifying young of the year and/or juveniles. In addition, cumulative marks may not be suitable for identification of juveniles. For example, bottlenose dolphins (Tursiops spp.) can be identified by scars, nicks or scratches on the dorsal fin [14] that are acquired by dolphins with age, therefore, adults are more likely to have such markings that make them individually identifiable.

There are exceptions where animals are identifiable from birth through natural markings (e.g., badgers Meles meles; [15], gray whales Eschrichtius robustus; [16]). However, it is more common that fewer juveniles than adults are individually identifiable, and those juveniles that are recognisable have markings that are distinctive or atypical [17], or have anomalies like albinism [18]. Furthermore, there can be heterogeneity between juveniles and adults in the proportion of individuals with distinctive markings, leading to juveniles being excluded from capture-recapture studies [17,19].

Monitoring juvenile movement and survival can be important to the management and conservation of species. Juvenile survival can show a greater sensitivity to environmental conditions than adult survival [20,21]. Furthermore, juveniles can make up a considerable proportion of a population [10], particularly in those populations that are growing [22]. Juvenile survival could therefore provide an early indicator of potential impacts on population growth rates and viability [23].

Demographic parameters for juveniles have been estimated indirectly using models such as life history or mortality tables [24], state space models [25] and matrix models [26]. Indirect methods do not track individual juveniles but instead deduce survival rates from the estimated numbers of animals in each age class and mean number of offspring per adult, relying on each individual’s progression through the demographic classes being largely constant and predictable. In cetaceans, other indirect methods for assessing calf and juvenile survival involve monitoring the association of offspring with photo-identifiable adult females, presumed to be their mother (e.g., [27]) and by evaluating the inter-calf interval in baleen whales. In the latter, a shorter than expected time interval between calving events implies that a calf has been lost and can be used as an indicator of calf mortality [28,29].

Direct methods, where individuals are marked and tracked, have the advantage of enabling fine-scale estimation such as differential survival between male and female juveniles, and are better able to accommodate intermittent field effort than indirect methods. Direct methods also allow juvenile survival to be estimated in species where juveniles take a long and variable time to reappear as adults in the sampled population, as is commonly the case for migratory whales and long-lived seabirds sampled at breeding grounds.

Here we focus on the use of genetic markers to monitor populations, which is becoming an increasingly important practice, particularly when there are conservation and management considerations [30,31]. Identification of individuals through genetic markers has several advantages over other methods, particularly for long-lived species. Unlike some forms of natural markings, genetic tags are stable and do not change over time. This makes them ideal to use for identifying individuals from birth.

Our study population is the New Zealand Southern right whale (Eubalaena australis), which was the subject of extensive commercial whaling that killed over 30,000 whales between 1827 and 1970 [32]. The contemporary New Zealand population is genetically distinct from Australian populations, and inhabits two wintering grounds where calving and breeding occur: the primary habitat in the New Zealand Subantarctic (centred at Port Ross, Auckland Islands) and the secondary habitat around mainland New Zealand (North and South Islands) [33].

Traditionally, photo-identification has been used to successfully identify non-calf right whales since the 1970s and focuses on callosity patterns found on the lip and rostrum, crenulations along the lower lip, and scars and unusual pigmentation on the head or body [34,35]. Such patterns do not typically stabilise until after 6 months of age, meaning dependent calves cannot normally be reliably identified from natural markings on their natal wintering ground [36]. Here we attempt to bridge this knowledge gap by identifying whales in their first 6 months of life using DNA profiles, comprising mitochondrial DNA haplotype, genetically identified sex and multilocus microsatellite genotype [37]. These data are used to investigate natal wintering ground fidelity for juvenile age classes and to provide the first estimate of juvenile survival in the New Zealand right whale population using a Cormack-Jolly-Seber (CJS) model [38].

Methods

Field and laboratory methods

Previous aerial and boat-based survey work indicates that the primary New Zealand wintering grounds of Southern right whales covers an area of approximately 20 km², limited to the waters of Port Ross, in the Subantarctic Auckland Islands (50°32' S, 166°15' E) and nearby coastal waters (Fig 1) [39,40]. All age and demographic classes are present in this wintering ground, meaning samples taken here are likely to be representative of the overall population [41,42]. Surveys were designed to coincide with peak abundance of Southern right whales in mid-July to early August [41,43] and were conducted during a total of 8 austral winters from 1995–1998 and 2006–2009.

[Figure omitted. See PDF.]

Fig 1. Location of primary Southern right whale breeding grounds and inset of study area at the Auckland Islands, New Zealand Subantarctic.

Map of four of the primary Southern right whale wintering grounds including Península Valdés, Argentina (ARG), South Africa (SAF), Australia (AUS) and the New Zealand Subantarctic (NZSA). Inset shows a map of the Auckland Islands, with a star showing the location of Port Ross, the primary survey area.

https://doi.org/10.1371/journal.pone.0146590.g001

As previously described [39], surveys were conducted from small vessels (4.6–5.2m). To avoid biases in sampling, search effort was distributed approximately evenly across the study area and attempts were made to approach every group sighted as described by Carroll et al. [44]. Skin biopsy samples were collected using small, stainless steel biopsy darts deployed from a crossbow in 1995–1998 [45] or a modified veterinary capture rifle in 2006–2009 [46] as previously described [33,47]. Biopsy samples were collected during the 2006–2009 field seasons under Department of Conservation (DOC) Marine Mammal Research permit and University of Auckland Animal Ethics Committee approved protocol to C. S. Baker. Biopsy samples were collected during the 1995–1998 field seasons under permit from the New Zealand DOC to C. S. Baker and N. Gales and University of Auckland Animal Ethics Committee approved protocol to C. S. Baker. Skin biopsy samples were preserved in 70% ethanol on location and transferred to the University of Auckland for storage at -20°C.

In the field, sampled whales were classified into two age groups based on their body length: non-calf or calf. The latter was defined as a whale whose portion of body visible at the surface was less than half of the length of an accompanying adult, and such an animal was assumed to be in its first year of life (age 0). A whale in close association with a calf was assumed to be its mother and was noted as a cow. In the subsequent modeling exercises, we define juveniles to be whales aged between 1 and 8 years. Age assignments were based solely on the progression of known-aged whales first captured as calves: in the field juveniles and adults are not consistently distinguishable. As reproductive maturity in the Southern right whale occurs by 9 years of age in females, we considered 8 years as the conservative maximum estimate of juvenile age [48–50].

DNA was extracted using standard proteinase K digestion and phenol/chloroform methods [51], as modified for small samples by Baker et al. [52]. DNA profiles, comprising genetically identified sex, mitochondrial control region (mtDNA) haplotypes and multi-locus genotypes, were constructed following protocols previously described [44]. Briefly, the sex of the sampled whale was identified by amplification of the male-specific SRY gene, multiplexed with an amplification of the ZFY/ZRX region as positive control [53,54]. Sequencing of the mitochondrial control region (500 bp) was conducted using the primers dlp1.5 [55] and tphe [33], the former modified with a 5’-M13 primer extension to facilitate subsequent sequencing reactions. The subsequent PCR product (~950 bp) was purified for sequencing with ExoSAP-IT (USB, Cleveland, Ohio, USA) and sequenced using BigDye Dye Terminator Chemistry (Applied Biosystems, Foster City, California, USA) on an ABI 3730 or ABI 3130 (Applied Biosystems). Sequences were aligned and edited in Sequencher v4.2 (Gene Codes Corporation, Ann Arbor, Michigan, USA) or Geneious [56], and haplotypes were identified from a 500 bp consensus region using haplotype codes established by Carroll et al. [33].

We undertook microsatellite genotyping using thirteen microsatellite loci (EV1, EV37, and EV14 [57]; GATA28 and GATA98 [58]; RW18, RW31, RW410, and RW48 [59]; GT23 [60]; TR3G1, TR3G2 and TR3F4 [61]), each amplified in individual 10 μL PCR reactions as previously described by Carroll et al. [33]. Each 96-well tray included a set of seven standard samples as an internal control to ensure consistent allele sizing and a negative control to detect contamination. Amplicons from 4–6 loci were co-loaded for capillary electrophoresis with an ABI 3730 or an ABI 3130. Alleles were sized with Genemapper v 4.0 (Applied Biosystems) and all automated calling was confirmed by visual inspection.

Identification of matching genotypes, i.e., recaptures, was conducted as described in Carroll et al. [33]. Briefly, a quality control measure of amplifying at a minimum of 9 loci was employed and matching genotypes were identified using CERVUS v3.0 [62]. As a precaution against false exclusion due to allelic dropout and other genotyping errors [63,64], the initial comparison allowed for mismatches at up to three loci.

We make the assumption that any calf captured at the Auckland Islands was born there, given the remote location and comparatively poor swimming capacity of right whale calves in their first few weeks of life [65]. We consider the calf to be recaptured in the Auckland Islands if it was sampled in a subsequent year.

Evaluating model assumptions

To estimate survival, we used an open capture-recapture model, the Cormack-Jolly-Seber (CJS) model. This capture-recapture model has several underlying assumptions (for a review see [19]), including (1) marks are permanent and can be correctly identified by researchers; (2) there is no behavioral response to capture; (3) there are equal survival and capture probabilities within demographic classes and (4) the population has geographic closure [66].

Individuals were identified using DNA profiles, which included a multilocus microsatellite genotype comprising an average of 12 loci per individual [37]. As previously reported, the least variable 11 loci provided a probability of identity (P_ID), or the chance that two individuals will have the same genotype by chance, of 7.9 x 10⁻¹⁴ [67], and a probability of identity for siblings, or the probability that two closely related individuals will have the same genotype by chance, of 1.73 x 10⁻⁵ [68]. Given these low probabilities, we considered that the suite of microsatellite loci could be used to confidently differentiate between even closely related individuals in a population estimated to have a superpopulation size of 2169 whales over the period 1995–2009 [37]. The low per allele error rate of 0.61% [37] and relaxed matching process means that it is unlikely that genotyping error would lead to a missed recapture event.

To evaluate assumption (2), behavioral response to capture, and to test for heterogeneity in survival probability between individuals, most commonly attributed to transiency, we used the program U-CARE [69]. To evaluate assumption (3), we tested whether there were differences in the recapture rates between males and females first sampled as calves and between the two sampling periods with Fisher’s exact test. Regarding assumption (4), the same wintering area, the Port Ross region of the Auckland Islands, was surveyed in each of the 8 expeditions [37].

The information from the recapture of adult individuals over time has been used to estimate the abundance of the population in both 1998 and 2009 [37,44], and to track movements of individuals between mainland New Zealand and the New Zealand Subantarctic [70]. However, existing analyses of recapture data, while being sex-specific, have not been linked to age classes. Furthermore, as calf survival is thought to be lower than that of adults in the Southern right whale [71], this age class was excluded from all previous analyses. In this study, all modeling was done using known-age whales that were first sampled in year of birth, i.e. first sampled as calves.

Mark recapture modeling

We used the CJS model, coded and implemented in R [72], to estimate the annual survival probability for calves, covering the year from age 0 to age 1, and for subsequent juvenile and adult age classes. The annual probability of apparent survival (Φ) was modelled as age-invariant Φ(.), or allowed to differ in the first year of life but remain constant thereafter (Φ_A2), following previous studies in other right whale populations which suggested that survival in the first year of life might be lower than that in subsequent years [71].

In addition, we considered a model denoted Φ_(curve) for continuous change in survival probability with age, in contrast with discrete age categories. We developed this model because Southern right whale calves grow rapidly in the first year of life [73] and the closely related North Atlantic right whale is estimated to reach 75% of its adult size by age one due to high maternal investment [74]. We postulate that survival may increase through the juvenile years due to the increased protection that size brings from predators such as killer whales and large sharks, eventually reaching a plateau at the adult survival level [48–50]. We used the following family of curves to model the relationship between age and apparent survival:where κ and σ are parameters to be estimated, and Φ(a) denotes the probability of apparent survival from age a − 1 years to age a. Parameter κ (0 < κ < 1) specifies the final adult survival rate at the plateau, and parameter σ (σ > 0) determines how quickly the survival curve reaches the adult level, with small values of σ indicating that juveniles quickly attain adult survival probabilities.

For each of the standard CJS and Φ_(curve) models, we explored models with constant capture probability (p(.)), and capture probability that varied with year or capture occasion (p(t)), and according to survey period (1995–1998 and 2006–2009), denoted p(90s,00s). In total we fitted nine models, combining each of three choices for survival (Φ(.), Φ_A2, and Φ_(curve)) with each of three choices for capture probability (p(.), p(t), and p(90s,00s)). Model-fitting was conducted by maximum likelihood, and we checked that our R results were identical to those from program MARK [75] for the six models with discrete age classes that could be constructed in program MARK.

In view of the small sample size of 66 capture histories, we used a bootstrap analysis [76] to construct confidence intervals. We randomly sampled the 66 capture histories with replacement to create 1000 bootstrapped datasets, each containing 66 capture histories and each representing animals first captured as calves. The model of interest was fitted to each bootstrapped dataset, generating 1000 estimates of each parameter including the survival probabilities Φ(a) for ages a = 1, 2, …14, with age 14 being the oldest possible age of animals in the study from 1995–2009. Bootstrapped 95% confidence intervals were taken from the 2.5% to the 97.5% quantiles of the 1000 estimates. We checked the bootstrap confidence intervals against analytic confidence interval computations, and found that in most cases the bootstrap confidence intervals were slightly wider. We used the bootstrap approach because it does not rely upon the large sample sizes that are assumed for analytic calculations, and because it is able to generate meaningful estimates of uncertainty in cases where analytic calculations fail.

The Akaike Information Criterion (AIC [77]), corrected for small sample sizes (AICc [78]), was used to assess support for each model [79]. As several models were closely ranked by AICc, we conducted an overall bootstrap analysis to give a single model-averaged estimate and confidence intervals for each survival probability, incorporating model-selection uncertainty. This approach is recommended by Buckland et al [80,81]. We generated 1000 bootstrapped datasets as before. To each bootstrapped dataset, we fitted all nine models, selected the best model based on AICc, and used the selected model to derive the age-based estimates of survival corresponding to that bootstrap iteration, namely Φ(a) for ages a = 1, 2, …14. The model-averaged estimate and standard error (SE) of each Φ(a) is given by the mean and standard deviation of the 1000 bootstrapped estimates, and the model-averaged 95% confidence interval for Φ(a) is taken from the 2.5% to the 97.5% quantiles of the 1000 estimates.

Results

Sample sizes and between-year recaptures

Over the 8 years of field work (1995–1998 and 2006–2009), 66 samples from calves produced DNA profiles that passed the quality control criterion of amplifying at a minimum of 9 of 13 loci. Of these, 34 were male and 32 were female; sample size and recaptures by years are shown in Table 1.

[Figure omitted. See PDF.]

Table 1. The number of between-year recaptures of Southern right whales first sampled as dependent calves during austral winter field surveys at the Auckland Islands from 1995–1998 and 2006–2009.

https://doi.org/10.1371/journal.pone.0146590.t001

During the 1995–1998 field survey period, 12 calves were sampled: 6 females and 6 males. Of these, 3 females and 1 male were recaptured as presumed adults in the 2006–2009 field surveys. Two female calves that were sampled in 1998 were subsequently recaptured as cows with calves in 2007. As female maturity in right whales is estimated at between 7 and 9 years [50,71,82], it is likely these records represent the first parturition for these females. Another female, first sampled as a calf in 1995, was subsequently recaptured in 1996, 1998 and 2006: each time without a calf. The male calf, first captured in 1996, was recaptured as a presumed adult, 11 years later, in 2007.

There were 54 dependent calves sampled during the 2006–2009 survey period: 26 females and 28 were males. Of these calves, 9 were sampled in 2009 and not available for subsequent recapture. Of the 45 available for recapture in at least one year, 9 were sampled in two years (20%; 4 males and 5 females), and one male was sampled in three years (2%). The P_ID of these matches was between 5.76 x 10⁻²³ and 4.70 x 10⁻¹⁴ (average 3.96 x 10⁻¹⁵), based on between 8 and 13 loci (average of 11 loci, S1 Table shows matching genotypes).

Estimates of calf survival

The U-CARE tests indicated no evidence of violations in the assumptions of mark-recapture models: both Test 3.SR (two sided test, p = 0.73), which tests for equal survival probability, and Test 2.CT (two sided test, p = 0.64), which tests for behavioural response to capture, were non-significant. There was no evidence of heterogeneity in capture probability linked to sex in the juvenile dataset. The recapture rate between females (5/28) and males (8/29) was not significantly different (2009 samples excluded as they could not be recaptured; Fisher’s exact test, p = 0.53). Similarly, the recapture rate of individuals sampled in the differing sampling periods was not significantly different (1995–1998: 4/12 and 2006–2009: 10/45, Fisher’s exact test p = 0.46). Therefore data were pooled for subsequent analyses.

Based on AICc, the Φ(.)p(.) CJS model was most highly ranked and produced an estimate of apparent annual survival of 0.96 for all ages (95% confidence intervals 0.85, 1.00: Table 2). However, several models were closely ranked (ΔAICc 1.75–2.16), including models that permit survival to vary with age. The fit of these models indicated that only the first two years of life has a decrease in survival probability relative to adults or other juvenile classes: the Φ_(curve) model had the flexibility to allow survival to take several years to plateau at the adult level, but it estimated the plateau to be reached by age three (see reference curves in Fig 2). This was also borne out by the bootstrap analysis, which estimated apparent annual survival to age one to be 0.87 (SE 0.17), increasing to 0.95 (SE 0.05) in subsequent years (Table 3). However, most of the models (80%) selected by the 1000 bootstrap replicates had age-invariant survival, indicating little support for a lower survival in the first year (Table 2).

[Figure omitted. See PDF.]

Fig 2. Family of curves used to model annual apparent survival as a function of age in the Φ(curve) models.

The models follow the equations, , where κ represents adult survival (κ = 0.98 in curves shown here) and σ controls the speed at which adult survival is attained. Curves where σ<1.0 fit the data best, indicating survival plateaued at a young age, compared with σ = 3.0 where survival plateaus at around age 8.

https://doi.org/10.1371/journal.pone.0146590.g002

[Figure omitted. See PDF.]

Table 2. Cormack-Jolly-Seber model estimates of apparent annual survival and 95% confidence intervals (95% CI) for juvenile New Zealand Southern right whales.

https://doi.org/10.1371/journal.pone.0146590.t002

[Figure omitted. See PDF.]

Table 3. Model averaged estimates of apparent annual survival of Southern right whales from New Zealand.

https://doi.org/10.1371/journal.pone.0146590.t003

It is notable that there is a hierarchical nature to the model ranking, with the models with constant capture probability ranking highest, followed by those models with survey-dependent capture probability and lastly models with time-dependent capture probabilities (Table 2). Within each capture probability configuration, the time-invariant survival model ranks first, the first-year survival model ranks second and the curve-based model ranks third.

Discussion

Here we have successfully used genetic tags in the form of DNA profiles to identify recaptures of juvenile Southern right whales first caught as calves, for the purpose of investigating juvenile survival and fidelity. In doing so, we highlight the use of genetic monitoring and maintaining a catalogue of DNA profiles as a method to address a gap in methodology: a method to reliably identify juveniles in species where natural markings are not sufficient to identify young individuals or where animals change as they develop.

Direct evidence of fidelity to natal wintering grounds in New Zealand Southern right whales

The recapture of male and female calves first sampled in the 1990s returning to the New Zealand Subantarctic up to 11 years later provides first direct evidence of long-term, maternally-directed fidelity to this wintering ground. This is consistent with long-term studies in conspecific populations that are based on the resighting of the small proportion of calves that are photo-identifiable through white dorsal pigmentation anomalies. For example, Best et al. [83] estimated that at least 93.4% of female calves born in the South African wintering ground return there to have their first calf. In Península Valdés, Argentina, 54 of 92 (59%) whales first identified as calves from 1971–1989 were resighted at the wintering ground by 1990 [84]. This study also suggested whales were more likely to be resighted when they were juveniles aged 1–4 (34% of possible resightings were observed) than as adults aged 4+ (23% of possible resightings observed). A comparison between the New Zealand and Argentinean wintering grounds is not possible as the eight-year gap between the two sets of consecutive surveys limited the opportunity to recapture calves for three years after birth. The only years where this was possible was in 1995 and 2006; of the 13 calves captured in these years, 5 (38%) were recaptured in the subsequent three years.

Estimating juvenile survival in New Zealand Southern right whales

Here we provide the first estimate of juvenile survival for the New Zealand population of Southern right whales. Using known-age whales and a flexible modelling framework, we were able to investigate changes in survival with age using discrete and continuous age classes. Overall, our bootstrapped results suggest apparent annual juvenile survival (ages 2 to 8) is high at 0.95 (SE 0.05), with only weak evidence for a decrease in apparent survival in the first year of life to 0.87 (SE 0.17).

The sample size used in this study was small (n = 66, 57 of which were available for recapture), a problem that is difficult to avoid in a population that is recovering from exploitation. This small sample size has implications for our ability to test for violations in model assumptions and impacts the precision of our estimates. For example, we have found low support for a decrease in survival in the first year of life in our model-average exercise: only 20% of bootstrap replicates selected a model that supported such a decrease. The small sample size means that we might not have had the power to detect such a pattern, which is a common mammalian trait [85]. Our ability to detect a difference in the recapture rate between males and females is also limited by sample size. However, while reproductively mature females show a clear pattern in recapture linked to reproductive cycle [37], there is no reason to expect such a pronounced difference in the pattern or rate of recaptures in non-reproductive females and males.

Finally, we did not investigate or account for individual-level heterogeneity. This could be investigated using a random effects model (e.g., [86]), but such models require larger sample sizes than we have available. Survival estimates from the CJS model are relatively robust to individual-level heterogeneity in capture probabilities [19,87]. Another potential problem is that we are estimating apparent survival, a composite of true survival and fidelity. Although we have presented evidence of long-term fidelity in the New Zealand population, there can be plasticity in this philopatric behaviour [88]. Therefore, we could be under-estimating the survival rate if fidelity is lower than expected and there is permanent emigration. Given the estimates of survival are high, this does not appear to be a major concern.

Despite our limitations, our estimates and precision are comparable to those from other studies on baleen whales, including Southern right whales. For example, first year female survival of Southern right whale was estimated at between 0.713 (95% confidence intervals 0.529, 0.896) and 0.956 (SE 0.051), using a multi-state matrix model that integrates capture histories of reproductively mature females captured at the South African wintering ground [48,71]. Female survival for ages 2+ was estimated to be 0.988 (SE 0.001) using the same model [48]. Estimates for the Península Valdés population using a similar modelling procedure suggested survival rate to reproductive maturity was 0.92 (SE 0.11) for females, implying an average annual mortality rate of calves and juveniles of about 0.01 (SE 0.01) [82].

Although estimates of juvenile survival in cetacean species are rare, our estimate of apparent survival in New Zealand right whales is higher than those available for gray whales: a study based on the capture-recapture of 129 photo-identified western gray whales (Eschrichtius robustus), which used Pollock's robust design model, estimated first year calf survival to be 0.701 (SE 0.0994, 95% confidence intervals 0.492, 0.850) and non-calf survival to be 0.951 (SE 0.0135, 95% confidence intervals 0.917, 0.972) [16]. Gabriele et al. [89] used an indirect method to calculate calf mortality in humpback whales by estimating the proportion of mothers identified with a calf in Hawai'i but without a calf in Alaska. This gave a point estimate of calf mortality of between 0.15 (95% confidence intervals 0.032, 0.378) and 0.241 (95% confidence intervals 0.103, 0.434), depending on whether all calf absences were attributed to mortality or a portion were attributed to weaning. Using an indirect method similar to that of Brandão et al. [48] that incorporates the calving intervals of humpback whales from the Gulf of Maine into a maximum likelihood framework, Barlow & Clapham [90] estimated humpback whale calf survival to be 0.875 (SE = 0.047). In an example from the small cetacean literature concerning bottlenose dolphins in Shark Bay, Australia, Mann et al. [27] estimated calf mortality by monitoring the association of offspring with photo-identifiable adult females, presumed to be their mother. In this study, 44% of bottlenose dolphin calves died by age 3, and a study of bottlenose dolphins at Doubtful Sounds, New Zealand, using similar methodology found calf survival to age 3 was 40% [91].

The high survival rate presented here and in the literature generally for Southern right whales may be due to analytical or biological effects. There is evidence to suggest that in the closely related species, the North Atlantic right whale, calves have a higher level of maternal investment than other baleen whale species, such as prolonged lactation period during which the calf can reach 75% of its expected adult length [74]. This optimised size at weaning may facilitate a higher juvenile survival in North Atlantic right whales, and the related Southern right whale, compared with other cetacean species.

Genotyping error and individual identification

Skin biopsy samples generally produce DNA of quality superior to that of typical ‘non-invasive’ sources of DNA. This is shown by the low error rate of the overall DNA profile dataset used in this study (0.61% per allele). This is similar to error rates estimated in studies that used tissue samples (0.8% for tissue [92],1.3% per locus and 0.8% per allele [93]). Studies which use non-invasively collected samples often have higher estimates of error, for example, error rates of 2.0%–11.3% and 2.55%–18.7% have been reported for studies using DNA extracted from faeces and hair, respectively [92,94,95]. We also attempted to minimise genotyping error from the start of the study by following published recommendations e.g., assessing DNA quantity and quality and ensuring samples are tracked accurately through the genotyping process [92,96].

However, even low error rates across large datasets, particularly those based on samples of varying quality, mean that genotyping error needs to be assessed and accounted for. Here we did this by estimating the genotype error rate and by matching samples in a relaxed framework. By identifying soft matches–those that match at multiple loci but that mismatch at up to three loci–we identify recaptures while allowing for genotyping error in the dataset. A study of misidentification error based on the adult captures from these field seasons, which were genotyped and matched by the same protocols, failed to find any evidence of matching errors [97].

Linking genetic samples to independent datasets of other tag types, such as natural markings or attached tags, has the potential to be an even more powerful method. Indeed, combining different tag types allows for animals to be tracked by different research groups over space and time, and is known to be statistically more informative than using a single tag [98,99].

Conclusion

The use of genetic tags allows individuals to be identified from birth, and is suitable for use with intermittent field effort, as was the case in this study. Genetic tagging is just one of many and varied methods to investigate fidelity and survival in calf and juvenile age classes. In the Southern right whale alone, there have been several indirect methods used, in addition to our direct method based on genetic tagging. This is an important topic as environmental conditions at high-latitude feeding grounds are hypothesised to affect population dynamics through reproductive success and juvenile survival in the Southern right whale [100,101]. The use of genetic monitoring means that there is no tag loss and all individuals can be identified from birth. We suggest integrating this method into the toolkit used to monitor the New Zealand Southern right whale and other baleen whale populations. In addition, increased sampling at both winter breeding and summer feeding grounds would enable more precise estimation of survival, and might enable application of covariate models to investigate impacts of environmental variables and specific feeding grounds.

Supporting Information

[Figure omitted. See PDF.]

S1 Appendix. Capture histories for southern right whales first captured as calves during field trips to the Auckland Islands, New Zealand Subantarctic in 1995–1998 and 2006–2009.

https://doi.org/10.1371/journal.pone.0146590.s001

(INP)

S1 Table. DNA profiles for southern right whales first captured as calves and recaptured in at least one subsequent year, during field trips to the Auckland Islands, New Zealand Subantarctic in 1995–1998 and 2006–2009.

https://doi.org/10.1371/journal.pone.0146590.s002

(DOCX)

Acknowledgments

Biopsy samples were collected during the 1995–1998 field seasons under permit from the New Zealand Department of Conservation (DOC) to C. S. Baker and N. Gales and University of Auckland Animal Ethics Committee approved protocol to C. S. Baker. Biopsy samples were collected during the 2006–2009 field seasons under DOC Marine Mammal Research permit and University of Auckland Animal Ethics Committee approved protocol to C. S. Baker. We thank G. Dunshea for assistance in planning and executing the 2006–2007 field seasons. Logistic support was provided by the Southland DOC Conservancy, the University of Auckland, and the Australian Antarctic Division. We thank Captain Steve Kafka and the crew of the Evohe for their help and support in the field. We thank B. Skerry, M. Handler, M. Vant, K. Munkres, S. Kraus, R. Rolland, L. Chilvers, C. Olavarría, P. McClelland, R. Cole, S. Cockburn, S. Smith, L. Douglas, W. Rayment, and T. Webster for help in the field. Thanks to R. Constantine for helpful comments on this manuscript.

Author Contributions

Conceived and designed the experiments: ELC RMF CSB. Performed the experiments: ELC SJC NJP LB CSB. Analyzed the data: ELC RMF. Contributed reagents/materials/analysis tools: ELC RMF SJC NJP LB CSB. Wrote the paper: ELC RMF SJC NJP LB CSB.

Citation: Carroll EL, Fewster RM, Childerhouse SJ, Patenaude NJ, Boren L, Baker CS (2016) First Direct Evidence for Natal Wintering Ground Fidelity and Estimate of Juvenile Survival in the New Zealand Southern Right Whale Eubalaena australis. PLoS ONE 11(1): e0146590. https://doi.org/10.1371/journal.pone.0146590

References

1. Norris K. Managing threatened species: the ecological toolbox, evolutionary theory and declining-population paradigm. J Appl Ecol. 2004;41: 413–426.

2. Gibbons JW, Andrews KM. PIT Tagging: simple technology at its best. Bioscience. 2004;54: 447.

3. Palsbøll PJ, Allen J, Bérubé M, Clapham PJ, Feddersen TP, Hammond PS, et al. Genetic tagging of humpback whales. Nature. 1997;388: 767–769. pmid:9285587

4. Taberlet P, Luikart G. Non-invasive genetic sampling and individual identification. Biol J Linn Soc. 1999;68: 41–55.

5. Charrier I, Pitcher BJ, Harcourt RG. Vocal recognition of mothers by Australian sea lion pups: individual signature and environmental constraints. Anim Behav. 2009;78: 1127–1134.

6. Ceugniet M, Izumi A. Vocal individual discrimination in Japanese monkeys. J Primatol. 2004;45: 119–28.

7. Loos A, Ernst A. An automated chimpanzee identification system using face detection and recognition. EURASIP J Image Video Process. 2013;1: 49.

8. Ogutu JO, Piepho H-P, Dublin HT, Reid RS, Bhola N. Application of mark-recapture methods to lions: satisfying assumptions by using covariates to explain heterogeneity. J Zool. 2006;269: 161–174.

9. Goswami VR, Madhusudan MD, Karanth KU. Application of photographic capture-recapture modelling to estimate demographic parameters for male Asian elephants. Anim Conserv. 2007;10: 391–399.

10. Karanth KU, Chundawat RS, Nichols JD, Kumar NS. Estimation of tiger densities in the tropical dry forests of Panna, Central India, using photographic capture–recapture sampling. Anim Conserv. 2004;7: 285–290.

11. Marshall A, Pierce S. The use and abuse of photographic identification of sharks and rays. J Fish Biol. 2012;80: 1361–1379. pmid:22497388

12. Kühl HS, Burghardt T. Animal biometrics: quantifying and detecting phenotypic appearance. Trends Ecol Evol. 2013;28: 432–41. pmid:23537688

13. Hammond PS. Capturing whales on film—estimating cetacean population parameters from individual recognition data. Mamm Rev. 1990;20: 17–22.

14. Wursig B, Jefferson T. Methods of photo-identification for small cetaceans. Rep Int Whal Comm Spec Issue. 1990;12: 43–53.

15. Dixon DR. A non-invasive technique for identifying individual badgers Meles meles. Mamm Rev. 2003;33: 92–94.

16. Bradford A, Wade P, Weller D, Burdin A, Ivashchenko Y, Tsidulko G, et al. Survival estimates of western gray whales Eschrichtius robustus incorporating individual heterogeneity and temporary emigration. Mar Ecol Prog Ser. 2006;315: 293–307.

17. Forcada J, Aguilar A. Use of photographic identification in capture-recapture studies on Mediterranean monk seals. Mar Mammal Sci. 2000;16: 767–793.

18. Burns D, Brooks L, Clapham P, Harrisson P. Between-year synchrony in migratory timing of individual humpback whales, Megaptera novaeangliae. Mar Mammal Sci. 2012;29: 228–235.

19. Pollock KH, Nichols JD, Brownie C, Hines JE. Statistical inference for capture-recapture experiments. Wildl Monogr. 1990;107: 1–97.

20. Tarwater C, RIcklefs R, Maddox J, Brawn J. Pre-reproductive survival in a tropical bird and its implications for avian life histories. Ecology. 2011;92: 1271–1281. pmid:21797155

21. Galliard J-M, Festa-Bianchet M, Yoccoz N, Loison A, Toigo C. Temporal variation in fitness components and population dynamics of large herbivores. Annu Rev Ecol Syst. 2000;31: 367–393.

22. Taylor BL, Chivers S, Larese J, Perrin N. Generation length and percent mature estimates for IUCN assessments of cetaceans. Administrative report LJ-07-01, National Marine Fisheries Service, Southwest Fisheries Science Center. 2007.

23. Dybala KE, Gardali T, Eadie JM. Dependent vs. independent juvenile survival: contrasting drivers of variation and the buffering effect of parental care. Ecology. 2013;94: 1584–93. pmid:23951718

24. Ebert T. Plant and Animal Populations. New York: Academic Press; 1998.

25. Buckland ST, Newman KB, Thomas L, Koesters NB. State-space models for the dynamics of wild animal populations. Ecol Modell. 2004;171: 157–175.

26. Neubert MG, Caswell H. Density dependent vital rates and their population dynamic consequences. J Math Biol. 2000;41: 103–141. pmid:11039693

27. Mann J, Connor RC, Barre LM, Heithaus MR. Female reproductive success in bottlenose dolphins (Tursiops spp.): life history, habitat, provisioning, and group-size effects. Behav Ecol. 1998;11: 210–219.

28. Weinrich M, Corbelli C. Does whale watching in Southern New England impact humpback whale (Megaptera novaeangliae) calf production or calf survival? Biol Conserv. 2009;142: 2931–2940.

29. Best P, Reeb D, Rew M, Palsbøll P, Schaeff C, Brandão A. Biopsying southern right whales: their reactions and effects on reproduction. J Wildl Manage. 2005;69: 1171–1180.

30. Schwartz MK, Luikart G, Waples R. Genetic monitoring as a promising tool for conservation and management. Trends Ecol Evol. 2007;22: 25–34. pmid:16962204

31. Baker CS. A truer measure of the market: the molecular ecology of fisheries and wildlife trade. Mol Ecol. 2008;17: 3985–3998. pmid:18643915

32. Carroll EL, Jackson JA, Paton D, Smith TD. Two intense decades of 19th century whaling precipitated rapid decline of right whales around New Zealand and east Australia. PLoS One. 2014;9: e93789. pmid:24690918

33. Carroll EL, Patenaude NJ, Alexander AM, Steel D, Harcourt R, Childerhouse S, et al. Population structure and individual movement of southern right whales around New Zealand and Australia. Mar Ecol Prog Ser. 2011;432: 257–268.

34. Best P. Natural marking and their use in determining calving intervals in right whales off South Africa. South African J Zool. 1990;25: 114–123.

35. Payne R, Brazier O, Dorsey E, Perkins J, Rowntree V, Titus A. External features in southern right whales (Eubalaena australis) and their use in identifying individuals. In: Payne R, editor. Communication and behaviour of whales. Boulder, Colardo: Westview Press; 1983. pp. 371–445.

36. Kraus SD, Moore MJ, Price CA, Crone M, Watkins W, Winn H, et al. The use of photographs to identify individual North Atlantic right whales (Eubalaena glacialis). Rep Int Whal Comm Spec Issue. 1986;10: 145–151.

37. Carroll EL, Childerhouse S, Fewster R, Patenaude NJ, Steel D, Dunshea G, et al. Accounting for female reproductive cycles in a superpopulation capture recapture framework. Ecol Appl. 2013;23: 1677–1690. pmid:24261048

38. Cormack R. Estimates of survival from the sighting of marked animals. Biometrika. 1964;51: 429–438.

39. Patenaude NJ, Baker CS. Population status and habitat use of southern right whales in the Subantarctic Auckland Islands of New Zealand. J Cetacean Res Manag Spec Issue. 2001;10: 111–116.

40. Rayment W, Davidson A, Dawson S, Slooten E, Webster T. Distribution of southern right whales on the Auckland Islands calving grounds. New Zeal J Mar Freshw Res. 2012;46: 431–436.

41. Patenaude NJ. Demographic and genetic status of right whales at the Auckland Islands, New Zealand. PhD thesis. University of Auckland, Auckland, New Zealand. 2002.

42. Carroll EL. Return of the Right Whale: Assessment of abundance, population structure and geneflow in the New Zealand southern right whale. PhD Thesis, University of Auckland, New Zealand. 2011.

43. Cawthorn M. New Zealand progress report on cetacean research, April 1991 to April 1992. Rep Int Whal Comm. 1993;43: 286–288.

44. Carroll EL, Patenaude NJ, Childerhouse SJ, Kraus SD, Fewster RM, Baker CS. Abundance of the New Zealand subantarctic southern right whale population estimated from photo-identification and genotype mark-recapture. Mar Biol. 2011;158: 2565–2575.

45. Lambertsen R. A biopsy system for large whales and its use for cytogenetics. J Mammal. 1987;68: 443–445.

46. Krützen M, Barré L, Möller L, Heithaus M, Simmer C, Sherwin WB. A biopsy system for small cetaceans; darting success and wound healing in Tursiops spp. Mar Mammal Sci. 2002;18: 863–878.

47. Patenaude NJ, Baker CS, Gales N. Observations of southern right whales on New Zealand’s subantarctic wintering grounds. Mar Mammal Sci. 1998;14: 350–355.

48. Brandão A, Butterworth D, Muller A, Best PB. Application of a photo-identification based assessment model to southern right whales in South African waters. Report SC/64/BRG24 presented to the Scientific Committee of the International Whaling Commission, Cambridge, UK. 2013. Available from iwc.int/home

49. Cooke J, Rowntree V, Payne R. Analysis of inter-annual variation in reproductive success of South Atlantic right whales (Eubalaena australis) from photo-identification of calving females observed off Península Valdés. Report SC/55/O23 presented to the Scientific Committee of the International Whaling Commission, Cambridge, UK. 2003. Available from iwc.int/home

50. Burnell SR. Estimates of demographic parameters of southern right whales off Australia. Report SC/60/BRG12 presented to the Scientific Committee of the International Whaling Commission, Cambridge, UK. 2008. Available from iwc.int/home

51. Sambrook J, Fritsch EF, Maniatis T. Molecular Cloning: A Laboratory Manual 2nd ed. Cold Spring Harbor, New York: Cold Spring Harbor Laboratory Press; 1989.

52. Baker CS, Slade RW, Bannister JL, Abernethy R, Weinrich M, Lien J, et al. Hierarchical structure of mitochondrial DNA gene flow among humpback whales, Megaptera novaeangliae, world-wide. Mol Ecol. 1994;3: 313–327. pmid:7921358

53. Aasen E, Medrano J. Amplification of the ZFX and ZFY genes for sex identification in humans, cattle, sheep and goats. Nat Biotechnol. 1990;8: 1279–1281.

54. Gilson A, Syvanen M, Levine K, Banks J. Deer gender determination by polymerase chain reaction: validation study and application to tissues, bloodstains, and hair forensic samples from California. Calif Fish Game. 1998;84: 159–169.

55. Baker CS, Medrano-Gonzalez L, Calambokidis J, Perry A, Pichler F, Rosenbaum H, et al. Population structure of nuclear and mitochondrial DNA variation among humpback whales in the North Pacific. Mol Ecol. 1998;7: 695–707. pmid:9640650

56. Biomatters. Geneious v6.1. Available: http://www.geneious.com. 2012.

57. Valsecchi E, Amos W. Microsatellite markers for the study of cetacean populations. Mol Ecol. 1996;5: 151–156. pmid:9147690

58. Palsbøll P, Bérubé M, Larsen AH, Jørgensen H. Primers for the amplification of tri- and tetramer microsatellite loci in baleen whales. Mol Ecol. 1997;6: 893–895. pmid:9301078

59. Waldick RC, Brown MW, White BN. Characterization and isolation of microsatellite loci from the endangered North Atlantic right whale. Mol Ecol. 1999;8: 1763–1765. pmid:10583845

60. Bérubé M, Jørgensen H, McEwing R, Palsbøll P. Polymorphic di-nucleotide microsatellite loci isolated from the humpback whale Megaptera novaeangliae. Mol Ecol. 2000;9: 2181–2183. pmid:11123643

61. Frasier TR, Rastogi T, Brown MW, Hamilton P, Kraus SD, White BN. Characterization of tetranucleotide microsatellite loci and development and validation of multiplex reactions for the study of right whale species (genus Eubalaena). Mol Ecol Notes. 2006;6: 1025–1029.

62. Kalinowski S, Taper M, Marshall T. Revising how the computer program CERVUS accommodates genotyping error increases success in paternity assignment. Mol Ecol. 2007;9: 801–888.

63. Waits JP, Leberg PL. Biases associated with population estimation using molecular tagging. Anim Conserv. 2000;3: 191–199.

64. Waits LP, Luikart G, Taberlet P. Estimating the probability of identity among genotypes in natural populations: cautions and guidelines. Mol Ecol. 2001;10: 249–256. pmid:11251803

65. Thomas PO, Taber SM. Mother-infant interaction and behavioral development in southern right whales, Eubalaena australis. Behaviour. 1984;88: 42–60.

66. White GC, Anderson DR, Burnham KP, Otis DL. Capture-recapture and removal methods for sampling closed populations. Los Alamos National Laboratory LA-8787-NERP; 1982.

67. Paetkau D, Strobeck C. Microsatellite analysis of genetic variation in black bear populations. Mol Ecol. 1994;3: 489–495. pmid:7952329

68. Evett I, Weir B. Interpreting DNA evidence: Statistical genetics for forensic scientists. Sinauer, Sunderland; 1998.

69. Choquet R, Lebreton J-D, Gimenez O, Reboulet A-M, Pradel R. U-CARE: Utilities for performing goodness of fit tests and manipulating CApture-REcapture data. Ecography. 2009;32: 1071–1074.

70. Carroll EL, Rayment W, Alexander AE, Baker CS, Patenaude NJ, Steel D, et al. Reestablishment of former wintering grounds by the New Zealand southern right whales. Mar Mammal Sci. 2014;30: 206–220.

71. Brandão A, Best P, Butterworth D. Estimates of demographic parameters for southern right whales off South Africa from survey data 1979 to 2006. Report SC/62/BRG30 presented to the Scientific Committee of the International Whaling Commission, Cambridge, UK. 2010. Available: https://iwc.int/home

72. R Development Core Team. R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing. Available: http://www.R-project.org; 2013.

73. Best P, Ruther H. Aerial photogrammetry of southern right whales Eubalaena australis. J Zool. 1992;228: 595–614.

74. Fortune SME, Trites AW, Perryman WL, Moore MJ, Pettis HM, Lynn MS. Growth and rapid early development of North Atlantic right whales (Eubalaena australis). J Mammal. 2012;93: 1342–1354.

75. White GC, Burnham KP. Program MARK: Survival estimation from populations of marked animals. Bird Study. 1999;46 S1: 120–138.

76. Efron B, Tibshirani RJ. An introduction to the bootstrap. New York: Chapman and Hall; 1993.

77. Akaike H. Information theory as an extension of the maximum likelihood principle. In: Petrov BN, Csaki F, editors. Second international symposium on information theory. Budapest, Hungary: Akademiai Kiado; 1973.

78. Hurvich C, Tsai C- L. Regression and time series model selection in small samples. Biometrika. 1989;76: 297–307.

79. Burnham KP, Anderson DR. Multimodel inference: Understanding AIC and BIC in model selection. Sociol Methods Res. 2004;33: 261–304.

80. Buckland ST, Burnham KP, Augustin NH. Model Selection: An Integral Part of Inference. Biometrics. 1997;53: 603–618.

81. Buckland S, Anderson DR, Burnham KP, Laake JL, Borchers DL, Thomas L. Introduction to distance sampling: Estimating abundance of biological populations. Oxford University Press; 2001.

82. Cooke J, Rowntree V, Payne R. Estimates of demographic parameters for southern right whales (Eubalaena australis) observed off Península Valdés, Argentina. J Cetacean Res Manag Spec Issue. 2001;2: 125–132.

83. Best PB. Coastal distribution, movements and site fidelity of right whales Eubalaena australis off South Africa, 1969–1998. South African J Mar Sci. 2000;22: 43–55.

84. Rowntree V, Payne R, Schell D. Changing patterns of habitat use by southern right whales (Eubalaena australis) on their nursery ground at Península Valdés, Argentina, and in their long-range movements. J Cetacean Res Manag Spec Issue. 2001;12: 133–143.

85. Caughley G. Mortality Patterns in Mammals. Ecology. 1966;47: 906–918.

86. Ford JH, Bravington M V., Robbins J. Incorporating individual variability into mark-recapture models. Methods Ecol Evol. 2012;3: 1047–1054.

87. Amstrup S, McDonald T, Manly B. Handbook of capture-recapture analysis. Princeton, New Jersey, USA: Princeton University Press; 2005.

88. Best P, Payne R, Rowntree V, Palazzo J, Do Carmo Both M. Long-range movements of South Atlantic right whales, Eubalaena australis. Mar Mammal Sci. 1993;9: 227–234.

89. Gabriele CM, Straley JM, Mizroch S, Baker CS, Craig AS, Herman LM, et al. Estimating the mortality rate of humpback whale calves in the central North Pacific Ocean. Can J Zool. 2001;79: 589–600.

90. Barlow J, Clapham P. A new birth interval approach to estimating demographic parameters of humpback whales. Ecology. 1997;78: 535–546.

91. Henderson SD, Dawson SM, Currey RJC, Lusseau D, Schneider K. Reproduction, birth seasonality, and calf survival of bottlenose dolphins in Doubtful Sound, New Zealand. Mar Mammal Sci. 2014;30: 1067–1080.

92. Bonin A, Bellemain E, Eidesen PB, Pompanon F, Brochmann C, Taberlet P. How to track and assess genotyping errors in population genetic studies. Mol Ecol. 2004;13: 3261–3273. pmid:15487987

93. Haaland O, Glover K, Sieliussen B, Skaug HJ. Genotyping errors in a calibrated DNA register: implications for identification of individuals. BMC Genet. 2011;12: 36. pmid:21507252

94. Broquet T, Menard N, Petit E. Noninvasive population genetics: a review of sample source, diet, fragment length and microsatellite motif effects on amplifcation success and genotyping error rates. Conserv Genet. 2007;8: 249–260.

95. Perez T, Vazquez F, Naves J, Fernandez A, Corao A, Abornoz J, et al. Non-invasive genetic study of the endangered Cantabrian brown bear (Ursus arctos). Conserv Genet. 2009;10: 291–301.

96. Morin P, Martien K, Archer F, Cipriano F, Steel D, Jackson JA, et al. Applied conservation genetics and the need for quality control and reporting of genetic data used in fisheries and wildlife management. J Hered. 2010;101: 1–10. pmid:19959596

97. Vale RTR, Fewster RM, Carroll EL, Patenaude NJ. Maximum likelihood estimation for model M t,α for capture-recapture data with misidentification. Biometrics. 2014;70: 962–971. pmid:24942186

98. Witteveen B, Straley J, Chenoweth E, Baker CS, Barlow J, Matkin C, et al. Using movements, genetics and trophic ecology to differentiate inshore from offshore aggregations of humpback whales in the Gulf of Alaksa. Endanger Species Res. 2011;14: 217–225.

99. Fewster RM, Jupp PE. Information on parameters of interest decreases under transformations. J Multivar Anal. 2013;120: 34–39.

100. Leaper R, Cooke J, Trathan P, Reid K, Rowntree V, Payne R. Global climate drives southern right whale (Eubalaena australis) population dynamics. Biol Lett. 2006;2: 289–292. pmid:17148385

101. Rowntree V, Uhart M, Sironi M, Chirife A, Di Martino M, La Sala L, et al. Unexplained recurring high mortality of southern right whale Eubalaena australis calves at Península Valdés, Argentina. Mar Ecol Prog Ser. 2013;493: 275–289