Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Refining estimates of availability bias to improve assessments of the conservation status of an endangered dolphin

  • Federico Sucunza ,

    Roles Conceptualization, Data curation, Formal analysis, Methodology, Project administration, Writing – original draft

    fsucunza@gmail.com

    Affiliations Laboratório de Ecologia Comportamental e Bioacústica, Programa de Pós-graduação em Ecologia, Universidade Federal de Juiz de Fora, Juiz de Fora, Minas Gerais, Brazil, Instituto Aqualie, Juiz de Fora, Minas Gerais, Brazil, Grupo de Estudos de Mamíferos Aquáticos do Rio Grande do Sul, Torres, Rio Grande do Sul, Brazil

  • Daniel Danilewicz,

    Roles Conceptualization, Funding acquisition, Methodology, Writing – review & editing

    Affiliations Instituto Aqualie, Juiz de Fora, Minas Gerais, Brazil, Grupo de Estudos de Mamíferos Aquáticos do Rio Grande do Sul, Torres, Rio Grande do Sul, Brazil, Universidade Estadual de Santa Cruz, Ilheus, Bahia, Brazil

  • Marta Cremer,

    Roles Methodology, Resources, Writing – review & editing

    Affiliation Laboratório de Nectologia, Departamento de Ciências Biológicas, Universidade da Região de Joinville, São Francisco do Sul, Santa Catarina, Brazil

  • Artur Andriolo,

    Roles Conceptualization, Funding acquisition, Methodology, Writing – review & editing

    Affiliations Laboratório de Ecologia Comportamental e Bioacústica, Programa de Pós-graduação em Ecologia, Universidade Federal de Juiz de Fora, Juiz de Fora, Minas Gerais, Brazil, Instituto Aqualie, Juiz de Fora, Minas Gerais, Brazil

  • Alexandre N. Zerbini

    Roles Conceptualization, Funding acquisition, Methodology, Supervision, Writing – review & editing

    Affiliations Laboratório de Ecologia Comportamental e Bioacústica, Programa de Pós-graduação em Ecologia, Universidade Federal de Juiz de Fora, Juiz de Fora, Minas Gerais, Brazil, Instituto Aqualie, Juiz de Fora, Minas Gerais, Brazil, Marine Mammal Laboratory, Alaska Fisheries Science Center, National Oceanic and Atmospheric Administration, Seattle, Washington, United States of America, Cascadia Research Collective, Olympia, Washington, United States of America, Marine Ecology and Telemetry Research, Seabeck, Washington, United States of America

Abstract

Estimation of visibility bias is critical to accurately compute abundance of wild populations. The franciscana, Pontoporia blainvillei, is considered the most threatened small cetacean in the southwestern Atlantic Ocean. Aerial surveys are considered the most effective method to estimate abundance of this species, but many existing estimates have been considered unreliable because they lack proper estimation of correction factors for visibility bias. In this study, helicopter surveys were conducted to determine surfacing-diving intervals of franciscanas and to estimate availability for aerial platforms. Fifteen hours were flown and 101 groups of 1 to 7 franciscanas were monitored, resulting in a sample of 248 surface-dive cycles. The mean surfacing interval and diving interval times were 16.10 seconds (SE = 9.74) and 39.77 seconds (SE = 29.06), respectively. Availability was estimated at 0.39 (SE = 0.01), a value 16–46% greater than estimates computed from diving parameters obtained from boats or from land. Generalized mixed-effects models were used to investigate the influence of biological and environmental predictors on the proportion of time franciscana groups are visually available to be seen from an aerial platform. These models revealed that group size was the main factor influencing the proportion at surface. The use of negatively biased estimates of availability results in overestimation of abundance, leads to overly optimistic assessments of extinction probabilities and to potentially ineffective management actions. This study demonstrates that estimates of availability must be computed from suitable platforms to ensure proper conservation decisions are implemented to protect threatened species such as the franciscana.

Introduction

Estimates of abundance are fundamental to properly assess the conservation status and to effectively design and implement management actions of threatened species or populations [1]. However, field-based research is commonly plagued by imperfect detection of the subjects of interest, often resulting in biased estimates of population parameters [2, 3]. For many marine species (e.g., cetaceans and chelonians) line transect aerial surveys constitute one of the most effective methods to estimate abundance [4]. In order to compute unbiased estimates, line transect methods assumes that detection on the survey trackline is certain (i.e., g(0) = 1) [2], but in aerial surveys visibility bias (i.e., g(0) < 1) often leads to underestimation of animals abundance [5, 6]. Marsh and Sinclair [7] coined the terms perception and availability bias to differentiate two forms of visibility bias during aerial surveys. Perception bias occurs when animals are available to be seen but are missed by the observers while availability bias correspond to animals that are missed because they are submerged and unavailable.

Availability bias can be a substantial source of bias for aerial surveys of marine mammals [5, 8, 9]. The high speed of the aircrafts reduces the amount of time a given area of the ocean is within the observer’s detectable view, resulting in a high proportion of undetected animals in this type of surveys. For example, Laake et al. [5] estimated that as many as 70% of harbor porpoises (Phocoena phocoena) would not be available to observers in an fixed-wing aircraft during aerial surveys in coastal waters of Washington State (USA). Accounting for this source of bias has, therefore, great relevance for population estimates, particularly of threatened species.

The franciscana, Pontoporia blainvillei, is a small dolphin endemic to coastal waters in the western South Atlantic Ocean. It occurs in waters typically shallower than 30m [10] from Itaúnas, Brazil (18º 25’S) to Golfo San Matías, Argentina (41º 10’S) [11, 12]. The species is considered the most threatened small cetacean in South America due to high, and possibly unsustainable, bycatch levels as well as increasing habitat degradation [13]. Incidental catches in fishing gear have been reported along most of the species’ range for almost seventy years [1315]. In order to guide conservation and management actions at a regional basis, the franciscana range was divided into four zones known as 'Franciscana Management Areas'—FMAs: two in southeastern and southern Brazil (FMA I and II), one in southern Brazil and Uruguay (FMA III) and one in Argentina (FMA IV) [16, 17]. The franciscana is currently listed as Vulnerable by the IUCN Red List of Threatened Species [18] and Critically Endangered by the Brazilian Government [19].

Because the franciscana is typically difficult to see in the wild [20], computing abundance estimates of this species has long been regarded as a challenge for conservation biologists [21, 22]. Line transect aerial surveys have been considered the most efficient method to estimate abundance of this species [23] and have been applied across almost their entire range [22, 2426]. Early estimates of franciscana abundance from aerial surveys [22, 24, 25] have been considered unreliable because, among other issues, they did not properly accounted for visibility bias [25, 27]. Dual platform experiments showed that franciscana density estimates from boat surveys are nearly 4–5 times greater than those obtained with aerial surveys, demonstrating that visibility bias (as well as underestimation of group sizes) represents a significant problem in aerial surveys for the species [28]. However, because this experiment was unable to separate the contribution of perception and availability bias to total visibility bias, further research was needed to better understand how these two sources of bias influence estimates of abundance.

To date, estimates of availability bias used to correct franciscana aerial surveys were computed from dive parameters collected from boats or land-based sites [24, 29]; and therefore may not properly reflect the amount of time franciscana groups are available to observers conducting aircraft-based surveys [25]. The main objective of this study is to estimate correction factors for availability bias using aerial platforms. In particular, helicopter surveys were used to estimate surface and dive intervals of franciscana groups as perceived from the air in order to (1) evaluate potential biological and environmental variables influencing dive parameters and (2) assess the magnitude of the differences in availability computed using data from surface-based and aerial platforms.

Methods

Study area

Aerial surveys were conducted in Babitonga Bay (26º 16’S, 048º 42’W), state of Santa Catarina, southern Brazil from 23 to 31 January 2014 (Fig 1). There is a number of advantages in conducting this study in this region. First, in Babitonga Bay, franciscanas predictably occur in relatively large densities throughout the year and show limited or no avoidance to small boats [30]. Second, group sizes seen in that region (mean = 2.90, range = 1 to 7 individuals, [28]) are representative of those seen throughout the range of the species (e.g., [11, 29]). Finally, the bay is relatively protected and therefore provides good weather conditions (e.g., relatively calm waters) for sighting surveys. In addition, previous studies to understand visibility bias have been conducted in this area [28], favoring a direct comparison of new information with existing data.

thumbnail
Fig 1. Study area.

Map of Babitonga Bay, State of Santa Catarina (red), southern Brazil. The inset shows the realized trackline effort and franciscana (Pontoporia blainvillei) sightings from helicopter surveys.

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

Data collection

Visual surveys were conducted from a Robinson R44 four-seat helicopter at 500 ft (150 m), an altitude consistent with that flown during aerial surveys to estimate abundance of this species (e.g., [25, 28]). Surveys had an average duration of 4 h and were carried out during the morning, when wind conditions were calm, optimal for this type of survey. The doors of the helicopter were detached during the study to maximize visibility for the observers. Two scientists with substantial experience in aerial surveys and familiar with the identification of franciscanas searched for groups of dolphins on the left side of the helicopter. Once a group was detected, the pilot hovered over it and each observer recorded surfacing and dive times independently. A group was defined as an aggregation of dolphins in close proximity of each other, moving in the same direction and in apparent association [31]. Visual and acoustic communication between the two observers were used to minimize simultaneous recording of the same group.

A record of sequences of surface-dive cycles started at the moment a group that had been detected became unavailable (dove), immediately after first being sighted, and ended when the observer lost the group. The first post-sighting interval at the surface was discarded from the analysis because it was not possible to ensure that the group was detected immediately after it became available. Inclusion of the first post-sighting interval could lead to underestimation of the surfacing interval and, consequently, to biased estimates of the length of the surface-dive cycle. Because the perspective from the helicopter makes it difficult distinguish individual franciscanas, observers recorded surfacings of any animal in the group. For each sampling period, the following information was registered: i) number of individuals visually available at, or near, the surface, ii) group size, iii) presence of calves and iv) a subjective assessment (categorized as: high, moderate, and low) of the level of confidence the observers had that they tracked the same group while collecting surfacing and dive times. In addition, the observers recorded ad libitum behavioral observations [32]. Based on the experience of the observers, "unusual" behaviors like an abrupt dive, a sudden turn away from the helicopter or an increase in swimming speed, were subjectively classified as a disturbance response to the presence of the helicopter.

Environmental conditions, including sea state (Beaufort scale), intensity and angle of glare, water color, and an assessment of the overall visibility conditions were recorded at the beginning and end of each sample or when the conditions changed. Data were recorded on audio digital recorders and time was synchronized to a GPS to allow for geo-referencing of each observation. Depth and water transparency (measured with a Secchi disc) at the location of each sighting were also recorded from boats operating in the same area and in radio communication with the helicopter. Tide level was obtained from the tide chart of the Port of São Francisco do Sul (26º 14’S, 048º 37’W).

Data analysis

Definition of dive parameters.

A surfacing interval was defined as the period of time in which at least one franciscana in the group was visually available to the observers in the helicopter at or near the surface. A diving interval was defined as the period of time in which all individuals of the group were not visible to the observers in the helicopter. A surface-dive cycle was defined as the period of time from the beginning of one surfacing to the next [6]. The proportion of time franciscana groups were visually available for detection from an aerial platform, hereafter referred to as “proportion at surface”, was computed as si/si+di, where si and di correspond, respectively, to the surfacing interval and diving interval times of the surface-dive cycle i (i = 1, 2,…, n surface-dive cycles).

Mixed-effects models.

Generalized linear mixed-effects models (GLMMs) [33,34] were used to evaluate the effects of biological and environmental predictors on the proportion at surface (the response variable) using package nlme [35] in the open source software R [36]. Because the same group may have been resampled during records of surface-dive cycles, individual franciscana groups were fitted as a random effect. The full set of models included all possible permutations of five predictor variables, which potentially affect the proportion at surface. This potential effect is interpreted here as the influence of the predictor variable on the behavior of the dolphins and/or the observers' ability to detect an individual. "Group size" describe the total number of individuals in the group. Because sample sizes were insufficient to fit models for all group size values, this variable was grouped into two categories small (groups with 4 or less individuals) and large (groups with more than 4 individuals). "Presence of calves" indicate the presence or absence of calves in the group. "Depth" and “Water Transparency” describe, respectively, the depth of the water column and the transparency at the location where the group was detected. "Tidal level" represents the level of the tide when surface-dive cycles were recorded. The proportion at surface was log-transformed to reduce the heterogeneity of the variance and to approximate the response variable to a normal distribution. Because multiple surface-dive cycles recorded for each franciscana group were temporally and spatial correlated, consideration was given to include an autocorrelation structure in the mixed effects models. Global models with (g1) and without (g2) autocorrelation structures were compared using the Akaike's Information Criterion AIC. Correlation structures of order 1 (AR-1) or an autocorrelation-moving average of order (p,q) were considered. The tested number of autoregressive parameters (p) and moving average parameters (q) varied between 1–2 and 0–1, respectively. To test whether the most supported autocorrelation structure, i.e., autocorrelation structure of order 1 (AR-1), was needed, the global model with (g1) and without (g2) the AR-1 structure were compared using AIC. Because the inclusion of the AR-1 structure resulted in a considerable improvement in AIC (g1 AIC = 515.792, g2 AIC = 520.344), all proposed models included the AR-1 structure.

The GLMM model selection was based on AIC, and the model with the smallest AIC value was considered the most parsimonious approximating model within the full set of models [37]. Akaike weights wi were calculated for each model as a representation of the probability of the model be the actual “best model” within the full set of models [37]. Inference about the relative importance (RI) of each predictor variable in determining the proportion at surface was based on the sum of Akaike weights of each variable across all candidate models containing the variable, and ranged from 1 (most important) to 0 (least important). Model averaging was conducted across the full set of models, and model-averaged parameters were estimated for each predictor variable, with unconditional standard errors incorporating model uncertainty [37]. Model averaging was performed using the package MuMIn [38].

Estimation of availability.

To estimate the probability of one franciscana group be visually available within the visual range of a passing observer in a fixed-wing aircraft, or availability (Pr) of franciscana groups, the method proposed by Laake et al. [5] was used: (1) where E(s), w(x), and E(d), correspond, respectively, to the mean time of each individual surfacing interval, the window of time during which a franciscana group is in the observer’s view at a distance x and the mean time of each individual diving interval. Availability was estimated for all groups together and for small and large groups.

w(x) was calculated following the McLaren’s [39] model: (2) where v is the mean aircraft velocity used in franciscana aerial surveys (v = 50 m/s), r is the radius of the observer’s view area (r = 300 m) and x is the distance from the survey line, assumed to correspond to zero for the computation of availability.

Standard errors and confidence intervals were estimated using a nonparametric bootstrap procedure [40]. For each bootstrap replicate, franciscana groups were sampled with replacement to match the original number of franciscana groups recorded (n = 101, see below) during the surveys and an estimate of Pr was obtained as in Eq 1. This process was repeated 1,000 times and the mean Pr and its standard error were computed from these replicates using standard formulae. Ninety-five percent confidence intervals were estimated by taking the 2.5 and 97.5 percentiles of the bootstrap replicates. Availability bias is equal to the complement of Pr.

Effects of platform type on the estimation of dive parameters.

To assess the effect of using dive parameters recorded from surface-based and aerial platforms to compute availability, the ratios of dive parameters (i.e., E(s) and E(d)) recorded in this study relative to those recorded from surface-based platforms [24, 29] was computed. Surface-based observations of franciscanas were obtained from shore-based platforms and boats in Anegada Bay, Argentina (see [29] for methods).

In addition, estimates of availability obtained with dive parameters from surface-based platforms (Pr[surface]) and applied to correct for availability bias in previously estimates of franciscana abundance from aerial surveys [22, 24, 26], were recomputed with data provided in this study (Pr[aerial]). For comparison, the same values and analytical methods as presented in the original studies were used. Secchi et al. [22] used dive parameters from Bordino et al. [29] and studies by Crespo et al. [24] and Zerbini et al. [26] used dive parameters from Crespo et al. [24]. Availability in all these studies was calculated with the method proposed by Barlow et al. [41].

Results

A total of 15 hours were flown during this experiment. Dolphins were well defined as visually available at some depth beneath the surface (Fig 2). A total of 337 complete surface-dive cycles was recorded for 120 franciscana groups in the original dataset. Franciscana groups appeared not to react to the presence of the helicopter very often. In only 2% (n = 6) of all recorded surface-dive cycles a response (e.g., an abrupt dive) was observed. After the removal of cycles were a reaction was documented as well as those classified as moderate/low confidence, the dataset for analysis was reduced to 248 complete surface-dive cycles from a total of 101 groups (mean = 2.45 cycles/group, SE = 1.99, range = 1–10). Biological and environmental variables recorded during this study are summarized in Table 1.

thumbnail
Fig 2. Franciscana (Pontoporia blainvillei) group.

Group classified as available to detection from the air.

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

thumbnail
Table 1. Summary of biological and environmental variables recorded in this study and tested in the generalized mixed-effects models.

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

Mixed-effects models

The most parsimonious approximating model to assess the influence of biological or environmental variables to the proportion at surface only included group size as the predictor variable (Table 2). This model suggested that time spent at the surface was positively correlated with the size of the group (p = 0.004). Group size was the most important predictor variable (RI = 0.96), followed by tidal level (RI = 0.36), water transparency (RI = 0.31), depth (RI = 0.30) and presence of calves (RI = 0.27) (Table 3). Except for group size, the other predictor variables had little influence in the estimated proportion at surface. The response variable slightly increased with increasing tide and depth and decreased with increasing water transparency and with the presence of calves in the group, but the effect of these variables was not statistically significant (Table 3).

thumbnail
Table 2. Most supported mixed-effects models (ΔAIC≤ 2) used to assess the influence of biological and environmental variables in the proportion at surface of franciscana groups.

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

thumbnail
Table 3. Model-averaged predictor coefficients and relative importance (RI).

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

Dive parameters and estimation of availability

Surfacing intervals varied from 1.03 seconds to 51.74 seconds (median = 16, mean = 16.10, SE = 9.74) (Fig 3A) and diving intervals varied from 0.46 seconds to 114.89 seconds (median = 36.86, mean = 39.77, SE = 29.06) (Fig 3B). The average proportion of time at surface was estimated at 0.36 (median = 0.29, SE = 0.22). The estimated window of time (w(x = 0)) was 6 seconds, which resulted in an estimation of availability of 0.39 (SE = 0.01) for all group sizes combined. Availability of large groups was significantly greater than that of small groups (Table 4).

thumbnail
Fig 3.

Histograms of surfacing (A, solid line represents the mean and the median) and dive (B, solid line represents the mean, and dashed line the median) times of franciscana (Pontoporia blainvillei) groups recorded from helicopter surveys in Babitonga Bay, southern Brazil.

https://doi.org/10.1371/journal.pone.0194213.g003

thumbnail
Table 4. Availability of franciscana (Pontoporia blainvillei) groups estimated for all groups together and for small and large groups separately assuming a window of time (w[x]) equal to 6 seconds.

https://doi.org/10.1371/journal.pone.0194213.t004

Effects of platform type on the estimation of dive parameters

The mean surface (E(s)) and dive (E(d)) times computed from the helicopter were greater than E(s) and E(d) computed from surface platforms at Anegada Bay, Argentina (Table 5). As a consequence, the use of dive parameters obtained at the surface level led to negatively biased estimates of availability of franciscanas for observers surveying from an aerial platform. The magnitude of the bias can be quite substantial, reaching as much as 46% depending on which dive parameters were used in the availability estimates (Table 6).

thumbnail
Table 5. Mean (SE in parenthesis) surfacing (E(s)) and diving (E(d)) times of franciscana groups recorded from different observation platforms (S = surface [land/boat], A = aerial [helicopter]).

https://doi.org/10.1371/journal.pone.0194213.t005

thumbnail
Table 6. Comparison between estimates of availability obtained with dive parameters from land/boat based surveys (Pr[surface], 20, 22, 24) and those computed with data provided in this study (Pr[aerial]) for aerial surveys conducted in three different FMAs.

https://doi.org/10.1371/journal.pone.0194213.t006

Discussion

Helicopter surveys in Babitonga Bay allowed for successful estimation of availability of franciscana groups. Dolphins were clearly defined as visually available or unavailable, facilitating the unambiguous estimation of the availability process [6]. However, the results presented in this study make the assumption that the helicopter did not influence the diving behavior of the franciscanas apart from those groups for which obvious reactions were observed. Because these responses were so infrequent (just 2% of our sample) and because behavior patterns were similar to those observed from boats in the same area [30, 43] and from the plane during aerial surveys, we conclude that the helicopter had a negligible influence on the behavior of the franciscanas for the majority of the dive cycles measured here. Experiments similar to the one conducted in this study with other species of dolphins of similar total body length of the franciscana (i.e., Phocoena phocoena [41] and Cephalorhyncus hectori [8]) indicated no evidence that the dolphins were disturbed by the presence of the helicopter supporting our conclusion. If our assumption is violated, however, estimates of dive parameters and availability presented here are biased. Studies (e.g., from surface platforms) comparing behavioral parameters such as time at surface and dive time in the presence and absence of the helicopter could be used to evaluate this assumption further.

Results from the GLMMs indicated that group size was the most important factor influencing the proportion at surface of a franciscana group. Because a surfacing interval was defined as the period of time in which at least one individual in a group was visually available to the observers at, or near, the surface, it is expected that an increase in the group size will increase the likelihood that the group is visually available to be detected. The positive effect of the group size on the availability observed in this study is not surprising and had been previously reported for other species of dolphins (e.g., [8]) and whales (e.g., [44]).

Data from time-depth recorders (TDRs) have been used to correct for availability bias in previous abundance estimates from aerial surveys [9, 45], but the use of these type of data has to be carefully considered depending on the characteristics of the species being surveyed. TDR data may provide information on the probability that the tagged individual is at, or near, the surface and available for detection, but they do not account for the behavior of all individuals in a group (in case animals are seen in groups). Lack of synchronous behavior may be a major concern about the use of data from TDRs to estimate availability because even when the tagged individual is underwater (and potentially unavailable), other individuals in the group may be, making that group detectable. If the size of the group has an effect on the probability of that group being available, such as observed in this study, estimates of availability computed from TDR dive data could misrepresent the true availability process of animals that do not swim in perfect coordination when in a group. Here, the average availability of small groups was significantly lower than that of large groups. However, the sample of large groups was small (n = 9, only about 10% of the total groups recorded), indicating that a more robust evaluation of the group size effect would be appropriate as more experiments are done and larger samples become available.

Use of the TDR data to assess availability of cetaceans can also be compromised by other factors such as water transparency and the characteristics of the body of the animal. In some circumstances, animals with conspicuous features (e.g., large body size and coloration that contrasts well with the water) could be detected from an aerial platform well below they break the surface or after they submerge following a surfacing. In these situations, the TDR data may indicate animals are unavailable when they are still visible. Experiments using TDRs attached to physical models, as the ones proposed by Pollock el at. [9], are important to evaluate these effects.

Another important aspect to consider with respect to availability of groups of different sizes is the need to accurately estimate the number of individuals in the group. Generally, determining the exact number of individuals is difficult and thus group size needs to be estimated. Such estimates are prone to biases from a number of factors including distance from the observer, behavior, body features (e.g., coloration and size), weather or sea conditions, and the size of the group itself (e.g., it is easier to estimate the number of animals in a group of 10 versus a group of 500 individuals). Therefore, before applying correction factors that are specific to group sizes or group size categories, the investigator shall ensure that the number of individuals in a group is accurately estimated during the survey. In this study, we assumed group size were estimated accurately because the helicopter provides an exceptional view to count individuals within the range of group size that franciscanas are observed in Babitonga Bay. Our assumption is further supported by group size estimates independently obtained during boat surveys [28]. The mean size and range of groups (mean = 2.90, SE = 1.24, range = 1–7) estimated by these independent surveys is identical to the ones computed from the helicopter.

Some of the environmental predictors examined in this study have been demonstrated to affect availability of marine mammals (e.g., water transparency, depth—[9, 46]). However, in the present experiment, they were neither included in the most parsimonious approximating model nor they seemed to influence the estimation of availability to a large extent. This can be explained, at least partially, by the relatively narrow range of the values recorded for some of these covariates in this study (e.g., water depth only varied between 4.5 and 12m). These ranges were limited because Babitonga Bay is a relatively confined environment and/or because the survey was conducted in a relatively short period of time (e.g, in the case of water transparency as it is expected to change during dry and rainy seasons). While many of the covariates analyzed in this study are representative of much of the habitat throughout the range of the franciscana, they may not have been sufficient to influence the behavior of the dolphins and/or the observers' capacity to detect an individual.

The effect of some environmental variables (e.g., Beaufort sea state) on estimates of the availability of marine mammals can be minimized by restricting the surveys to optimal survey conditions, however some variables such as water transparency and depth can vary across the species habitat and can change within few minutes of aerial survey time [9] (e.g., as the sampling platform moves between inshore and offshore habitats). This is the case for the franciscana, a species inhabiting habitats with very turbid waters near river mouths as well as clearer waters in the open ocean [24]. While the estimates of availability provided here represent a significant progress towards computing more reliable correction factors for estimates of abundance based on aerial surveys, it is important to recognize that their use could lead to bias in abundance if availability varies significantly across these habitats [47]. In this sense, experiments similar to the one conducted here are recommended for other habitats along the range of the franciscana where environmental characteristics differ from those at Babitonga Bay. These experiments would contribute to better understand the influence of spatially and temporally-varying environmental predictors on the availability of franciscanas to aerial platforms and lead to more robust estimates of availability.

Account for biological and/or environmental-specific availability correction factors in abundance estimates is challenging and highly dependent on proper field methods and robust survey designs. Dynamic environmental variables that can change dramatically while flying (e.g., water transparency) as well as factors that require estimation (e.g., group size), are of great concern when availability correction for different levels of these covariates is expected and when fine scale data on the variability of these covariates are not available. However, fine-scale correction factors may not be realistic on many surveys, and a lack of fine scale measures can be accommodated depending principally on the spatial extend of the survey [47]. If the survey area can be stratified in a way that availability vary among strata, but remain relatively uniform within it, broad categories of stratification could result in substantial improvement on estimates [47, 48]. Abundance estimates corrected for biological and/or environmental-specific availability estimates can be computed using the methods proposed by Pollock et al. [9].

A comparison of the dive parameters sampled from a helicopter in Babitonga Bay [this study] with those obtained from surface platforms in Anegada Bay, Argentina [24, 29] highlights the importance of collecting dive data from aerial platforms in order to develop appropriate availability correction factors for estimates of abundance computed from aerial surveys. Although regional differences in the behavior of franciscanas and in the environmental characteristics between the two areas could influence these results, it is clear that animals not only become available to observers in the air for some time before they break the surface, but also remain visible between short dives [8, 41, this study]. Therefore, the surfacing interval estimated by platforms at the water level (e.g., land or a boat) is expected to be much shorter. In fact, results of this study showed that the surfacing interval recorded from the helicopter were, on average, 13 times greater than the surfacing time recorded from a surface-based platform. The use of surfacing-diving intervals computed from a surface platform may result in (often substantially) negatively biased estimates of availability for aerial platforms. More important, correction factors developed from such estimates of availability would lead to an overestimation of abundance in aerial surveys.

Estimates of franciscana abundance to data have been computed for all FMAs [22, 2426, 49]. Those for FMA II [26], FMA III [22, 25] and FMA IV [24] have been corrected for availability bias using correction factors computed from the surface-dive intervals obtained through surface platforms. Because the results of this study show that these correction factors were negative biased (by as much as 46%), re-computed abundance estimates for these regions is appropriate, taking the availability estimates provided here into account.

Recalculations are not attempted in this study for the estimates from the FMA III [22, 25] or the FMA IV [24] because these did not account for perception bias and are therefore biased irrespective of how correction for availability is calculated. On the other hand, a re-evaluation of the estimate for the FMA II [26] is possible because an independent estimate of perception bias is available. The estimate originally provided for this region in 2008 (N = 8,525, CV = 0.34) would be reduced by 38% (N = 6,146, CV = 0.35) when availability computed with surface-dive cycles sampled from the aerial platform is used. The estimation of bycatch removals illustrates the management implications of this bias for the franciscanas in FMA II. The estimated bycatch in this area ranged between 300 and 500 individuals per year in the 2000s [14, 27, 50]. This represented an annual removal of 3.5–5.8% of the population as originally estimated by Zerbini et al. [26]. However, removal was substantially larger (4.8–8.6%/year) when calculated using the estimate corrected for availability computed with surface-dive data from the aerial platform, clearly indicating a significantly greater conservation problem.

Conclusions

The results presented here show that availability bias can account for a large proportion of missing animals during aerial surveys (nearly 60%). In addition, they demonstrate that estimation of availability correction factors for aerial surveys should be computed with data collected from aerial platforms, not from land or boats. Although the results provided here can be refined, they represent a substantial improvement on availability correction factors previously used for franciscana aerial surveys [22, 2426]. The conservation of the franciscana is a matter of great concern. The species is listed as vulnerable by the IUCN [18] and mortality resulting from incidental catches is thought to be unsustainable through its range [13, 14]. Therefore, reliable data on abundance is needed to assess the status of all franciscana population and to establish proper conservation actions. In this sense, the use of the availability correction factors presented here should be encouraged to correct for availability bias of fransicana in aerial surveys and, consequently, obtain more realistic estimates of the impact of bycatch and other human activities to this endangered species.

Acknowledgments

The authors are grateful to all colleagues who assisted in field activities: Ana K. Moraes, Annelise C. Holz, Beatriz Schulze, Camila M. Sartori, Natacha Zimmermann, Pamela S. Actis, Renan Paitach. We thank Dr. Jeff Laake for assistance with statistical analysis. The pilot of the helicopter Paulo Camejo was key to the success of this study.

References

  1. 1. Williams BK, Nichols JD, Conroy MJ. Analysis and management of animals populations: Modeling, estimation, and decision making. San Diego: Academic Press; 2001.
  2. 2. Buckland ST, Anderson DR, Burnham KP, Laake JL, Borchers DL, Thomas L. Introduction to distance sampling: Estimating abundance of biological populations. Oxford: Oxford University Press; 2001.
  3. 3. Schlossberg S, Chase MJ, Griffin CR. Testing the accuracy of aerial surveys for large mammals: An experiment with african savanna elephants (Loxodonta africana). PloS one. 2016; 11: e0164904. pmid:27755570
  4. 4. Thomas L, Buckland ST, Rexstad E, Laake JL, Strindberg S, Hedley SL, et al. Distance software: design and analysis of distance sampling surveys for estimating population size. J App Ecol. 2010; 47: 5–14. pmid:20383262
  5. 5. Laake JL, Calambokidis J, Osmek SD, Rugh DJ. Probability of detecting harbor porpoise from aerial surveys: estimating g(0). J Wildl Manage. 1997; 61: 63–75.
  6. 6. Laake JL, Borchers DL. Methods for incomplete detection at distance zero. In: Buckland ST, Anderson DR, Burnham KP, Laake JL, Borchers DL, Thomas L, editors. Advanced distance sampling: Estimating abundance of biological populations. Oxford: Oxford University Press; 2004. pp. 108–189.
  7. 7. Marsh H, Sinclair DF. Correcting for visibility bias in strip transect aerial surveys of aquatic fauna. J Wildl Manage. 1989; 53: 1017–1024.
  8. 8. Slooten E, Dawson SM, Rayment WJ. Aerial surveys for coastal dolphins: abundance of hector's dolphins off the South Island West Coast, New Zealand. Mar Mammal Sci. 2004; 20: 477–490.
  9. 9. Pollock KH, Marsh H, Lawler IR, Alldredge MW. Estimating animal abundance in heterogeneous environments: an aplication to aerial surveys for dugongs. J Wildl Manage. 2006; 70: 255–262.
  10. 10. Danilewicz D, Secchi ER, Ott PH, Moreno IB, Bassoi M, Borges-Martins M. Habitat use patterns of franciscana dolphins (Pontoporia blainvillei) off southern Brazil in relation to waters depth. J Mar Biol Assoc UK. 2009; 89: 943–949. https://doi.org/10.1017/S002531540900054X
  11. 11. Crespo EA, Harris G, Gonzáles R. Group size and distributional range of the franciscana, Pontoporia blainvillei. Mar Mammal Sci. 1998; 14: 845–849.
  12. 12. Siciliano S, Di Beneditto APM, Ramos RMA. A toninha, Pontoporia blainvillei, (Gervais & d'Orbigny, 1844) (Mammalia, Cetacea, Pontoporiidae), nos estados do Rio de Janeiro e Espírito Santo, costa sudeste do Brasil: caracterização dos hábitats e fatores de isolamento das populações. Boletim do Museu Nacional, Zoologia. 2002; 476: 1–5.
  13. 13. Secchi ER, Ott PH, Danilewicz D. Effects of fishing bycatch and the conservation status of the franciscana dolphin, Pontoporia blainvillei. In: Gales N, Hindell M, Kirkwood R, editors. Marine mammals: Fisheries, tourism and management issues. Collingwood: CSIRO Publishing; 2003. pp. 174–191.
  14. 14. Ott PH, Secchi ER, Moreno IB, Danilewicz D, Crespo EA, Bordino P, et al. Report of the working group on fishery interactions. Lat Am J Aquat Mamm. 2002; 1: 55–64. http://dx.doi.org/10.5597/lajam00008
  15. 15. Van Erp I. In quest of the La Plata dolphin. Pacific Discovery. 1969; 22: 18–24.
  16. 16. Secchi ER, Danilewicz D, Ott PH. Applying the phylogeographic concept to identify franciscana dolphin stocks: Implications to meet management objectives. J Cetacean Res Manag. 2003; 5: 61–68.
  17. 17. Anonymous, 2015. Report of the VIII Workshop for the Research and Conservation of the 477 Franciscana Dolphin (Pontoporia blainvillei) in the Soutwestern Atlantic. UNEP/CMS, 478 São Francisco do Sul, Brazil.
  18. 18. Zerbini, A.N., Secchi, E., Crespo, E., Danilewicz, D. & Reeves, R. 2017. Pontoporia blainvillei. In: IUCN 2017. The IUCN Red List of Threatened Species. 2017. Available: www.iucnredlist.org. Acessed 30 January 2018.
  19. 19. MMA. Lista nacional oficial de espécies da fauna ameaçada de extinção—Mamíferos, Aves, Répteis, Anfíbios e Invertebrados Terrestres. In: Portaria MMA no. 444, de 17 de dezembro de 2014, Brasil. Available from http://www.icmbio.gov.br/cepsul/images/stories/legislacao/Portaria/2014/p_mma_444_2014_lista_esp%C3%A9cies_ame%C3%A7adas_extin%C3%A7%C3%A3o.pdf. Accessed 27 July 2017.
  20. 20. Crespo EA. Fransicana dolphin Pontoporia blainvillei. In: Perrin WF, Würsig B, Thewissen JGM, editors. Encyclopedia of marine mammals. 2nd ed. San Diego: Academic Press; 2009. pp. 466–469.
  21. 21. Praderi R, Pinedo MC, Crespo EA. Conservation and management of Pontoporia blainvillei in Uruguay, Brazil and Argentina. In: Perrin WF, Brownell RL Jr., Zhou K, Liu J, editors. Biology and Conservation of the River Dolphins. Hong Kong: IUCN; 1989. pp. 52–56.
  22. 22. Secchi ER, Ott PH, Crespo E, Kinas PG, Pedraza S, Bordino P. A first estimate of franciscana (Pontoporia blainvillei) abundance off southern Brazil. J Cetacean Res Manag. 2001; 3: 95–100.
  23. 23. Crespo EA, Secchi ER, Rosa LD, Kinas PG, Danilewicz D, Bordino P. Report of the working group on abundance estimates. Lat Am J Aquat Mamm. 2002; 1: 65–66. http://dx.doi.org/10.5597/lajam00009
  24. 24. Crespo EA, Pedraza SN, Grandi MF, Dans SL, Garaffo GV. Abundance and distribution of endangered Franciscana dolphins in Argentina waters and conservation implications. Mar Mammal Sci. 2010.; 26: 17–35.
  25. 25. Danilewicz D, Moreno IB, Ott PH, Tavares M, Azevedo AF, Secchi ER, et al. Abundance estimate for a threatened population of franciscana dolphins in southern coastal Brazil: uncertainties and management implications. J Mar Biol Assoc UK. 2010; 90: 1659–1666. https://doi.org/10.1017/S0025315409991482
  26. 26. Zerbini AN, Secchi ER, Danilewicz D, Andriolo A, Laake JL, Azevedo AF. Abundance and distribution of the franciscana (Pontoporia blainvillei) in the Franciscana Management Area II (southeastern and southern Brazil). International Whaling Commission, Scientific Committee Paper SC/62/SM7. 2010.
  27. 27. International Whaling Commission. Report of the sub-committe on samll cetacean. Annex L. J Cetacean Res Manag. 2005; Supplement: 307–326.
  28. 28. Zerbini AN, Danilewicz D, Secchi ER, Andriolo A, Flores PAC, Cremer MJ, et al. Assessing bias in abundance estimates from aerial surveys to improve conservation of threatened franciscana dolphins: preliminary results from a survey conducted in southern Brazil. International Whaling Commission, Scientific Committee Paper SC/63/SM9. 2011.
  29. 29. Bordino P, Thompson G, Iñíguez M. Ecology and behaviour of the franciscana (Pontoporia blainvillei) in Bahía Anegada, Argentina. J Cetacean Res Manag. 1999; 1: 213–222.
  30. 30. Cremer MJ, Simões-Lopes PC. Distribution, abundance and density estimates of franciscanas, Pontoporia blainvillei (Cetacea:Pontoporiidae), in Babitonga bay, southern Brazil. Rev Bras Zool. 2008; 25: 397–402. http://dx.doi.org/10.1590/S0101-81752008000300003
  31. 31. Shane SH. Behavior and ecology of the bottlenose dolphin at Sanibel Island, Florida. In: Leatherwood S, Reeves RR, editors. The bottlenose dolphin. San Diego: Academic Press, 1990. pp. 245–265.
  32. 32. Altmann J. Observational study of behavior: Sampling methods. Behaviour. 1974; 49:227–267. http://www.jstor.org/stable/4533591 pmid:4597405
  33. 33. Pinheiro JC, Bates DM. Mixed-Effects Models in S and S-Plus. New York: Springer; 2000.
  34. 34. Zuur AF, Ieno EN, Walker NJ, Saveliev AA, Smith GM. Mixed Effects Models and Extensions in Ecology with R. New York: Springer; 2009.
  35. 35. Pinheiro J, Bates D, DebRoy S, Sarkar D, R Core Team. nlme: Linear and Nonlinear Mixed Effects Models. R package version 3.1–131. 2017.
  36. 36. R Development Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. 2017. Available: http://www.Rproject.org/
  37. 37. Burnham KP, Anderson DR. Model selection and multimodel inference: a practical information-theoretic approach. 2nd ed. New York: Springer; 2002.
  38. 38. Barton K. MuMIn: Multi-Model Inference. R package version 1.40.0. 2017.
  39. 39. McLaren IA. Methods of determining the numbers and availability of ringed seals in the eastern Canadian Arctic. Arctic. 1961; 14: 162–175.
  40. 40. Manly BFJ. Randomization, Bootstrap and Monte Carlo Methods in Biology: 3rd ed. Boca Raton: Chapman & Hall/CRC; 2006.
  41. 41. Barlow J, Oliver CW, Jackson TD, Taylor BL. Harbor porpoise, Phocoena phocoena, abundance estimation for California, Oregon, and Washington: II. aerial surveys. Fishery Bulletin. 1988; 86: 433–444.
  42. 42. Bordino P. Diving behavior. In: Crespo EA. Key aspects of the biology and ecology of the Franciscana dolphin (Pontoporia blainvillei). Final report to CMS; 2004. pp. 18–23.
  43. 43. Cremer MJ, Simões-Lopes PC. The occurrence of Pontoporia blainvillei (Gervais & d'Orbigny) (Cetacea, Pontoporiidae) in an estuarine area in southern Brazil. Rev Bras Zool. 2005; 22: 717–723. http://dx.doi.org/10.1590/S0101-81752005000300032
  44. 44. Hain JHW, Ellis SL, Kenney RD, Slay CK. Sightability of right whales in coastal waters of the southeastern United States with implications for the aerial monitoring program. In: Garner GW, Amstrup SC, Laake JL, Manly BFJ, McDonald LL, Robertson DG, editors. Marine Mammal Survey and Assessment Methods. Rotterdam: A. A. Balkema; 1999. pp. 191–207.
  45. 45. Heide-Jørgensen MP, Laidre K, Borchers D, Samarra F, Stern H. Increasing abundance of bowhead whales in West Greenland. Biol Lett. 2007; 3: 577–580. pmid:17686750
  46. 46. Hagihara R, Jones RE, Grech A, Lanyon JM, Sheppard JK, Marsh H. Improving population estimates by quantifying diving and surfacing patterns: A dugong example. Mar Mammal Sci. 2014; 31: 348–366.
  47. 47. Thomson JA, Cooper AB, Burkolder DA, Heithaus MR, Dill LM. Heterogeneous pattern of availability for detection during visual surveys: Spatiotemporal variation in sea turtle dive-surfacing behaviour on a feeding ground. Methods Ecol Evol. 2012; 3: 378–387.
  48. 48. Thomson JA, Cooper AB, Burkholder DA, Heithaus MR, Dill LM. Correcting for heterogeneous availability bias in surveys of long-diving marine turtles. Biol Conserv. 2013; 165: 154–161.
  49. 49. Danilewicz D, Zerbini AN, Andriolo A, Secchi ER, Sucunza F, Ferreira E, et al. Abundance and distribution of an isolated population of franciscana dolphins (Pontoporia blainvillei) in southeastern Brazil: red alert for FMA I? International Whaling Commission, Scientific Committee Paper SC/64/SM17. 2012.
  50. 50. Bertozzi CP, Zerbini AN. Incidental mortality of franciscana (Pontoporia blainvillei) in the artisanal fishery of Praia Grande, São Paulo State, Brazil. Lat Am J Aquat Mamm. 2002; 1: 153–160.