Utilizing Occupancy Models and Platforms-of-Opportunity to Assess Area Use of Mother-Calf Humpback Whales

The Hawaiian Islands, and particularly the Maui 4-island region, are a critical breeding and calving habitat for humpback whales (Megaptera novaeangliae) belonging to the Hawaii distinct population segment. Our aims were to test the use of platforms-of-opportunity to determine trends in mother-calf pod use of the region and to present opportunistic platforms as an alternative method of long-term, cross-seasonal monitoring. Data were collected from whale watching vessels over a 4-year period and analyzed using occupancy models to determine the probability of habitat use of pods with calves and pods without calves within the study area. Detection probability was influenced by survey effort and month for all pod types with detection of adult only pods further influenced by year. Pods with a calf showed a preference for shallow (<100 meters) low latitude waters (<20.7 ̊N), while pods without a calf preferred deeper waters (>75 meters). Results presented here align with previous work, both in Hawaii and in other breeding grounds, which show a distinct segregation of mothers with a calf from other age-classes of humpback whales. The need for long-term continuous monitoring of cetacean populations is crucial to ensure species conservation. Data collected aboard platforms-of-opportunity, as presented here, provide important insight on humpback whale spatial and temporal distribution, which are essential for species protection and management.


Introduction
For the humpback whale (Megaptera novaeangliae), mother-calf associations are considered the strongest bond within an otherwise solitary species, and maintenance of this association is crucial to calf survival and eventual recruitment to the population [1] [2].As such, understanding how and when mother-calf humpback whale pods use a specific area is important for the ecology of the species as well as its conservation and management, particularly for the Hawaiian population segment which was removed from the Endangered Species Act in 2016 [3].Winter aggregations of humpback whales in Hawaiian waters belong to the Hawaii distinct population segment that utilize the main Hawaiian Islands as breeding and calving grounds and constitutes a significant portion: ~10,000 [4] of the global population: ~90,000 [5].The highest densities of Hawaiian humpbacks occur within the Maui 4-island region [6] [7], consisting of Maui, Lana'i, Moloka'i, and Kaho'olawe.Spatial segregation of mother-calf pods within the Hawaiian breeding grounds has been observed (e.g., [8] [9] [10] with the Maui 4-island region being the preferred habitat for mother-calf pods in this population [11] [12]).
Research on humpback whales in Brazil, Madagascar, Ecuador, and Costa Rica found that mother-calf pairs used shallower waters than groups without calves [10] [13] [14] [15].Mother-calf pods may prefer shallower waters to avoid energetically costly situations such as aggressive attention from male suitors, rough waters that complicate nursing, or competition groups which may put the calf in danger [9] [10] [13] [14] [15] [16].Observations of mother-calf pairs in Hawaii have reported similar findings across all islands.Pods containing mothers and calves near the island of Hawaii were observed in waters significantly shallower than groups of one or two adults [9].Similarly [17] observed pods with a calf in shallower waters than other pod types observed on the north shore of Kaua'i.
Based on observations of pods with calves in the deeper waters of the Au'au channel between the islands of Maui and Lana'i, It has been suggested females with calves avoid the shallowest areas within 2 kilometers of the coastline [18].Although this pattern contrasts other studies of mother-calf habitat preference for shallow waters, limitations were expressed by authors including a short survey season and timing of the study [18], potentially contributing to the contradicting results from other studies in the Hawaiian Islands and around the world [9] [10] [13] [14] [15] [16].Observed variation in mother-calf habitat preference among breeding areas could be caused by differences in behavioral context, anthropogenic stress, or sampling methodologies [18].Knowledge of mother-calf distribution throughout the Maui 4-island region could inform assessment of restrictions, including but not limited to vessel speed and area use, that may be placed to reduce animal disturbance and harm.This study compares the spatial distribution of humpback whale pods with and without calves over a 4-year period throughout the winter breeding season (December to April).We present on the analysis of humpback whale sighting data collected from platforms-of-opportunity to relate humpback whale pres-Open Journal of Marine Science ence/absence to environmental characteristics and time of year in driving the observed distribution within the Maui 4-island region, utilizing occupancy modelling techniques [19] [20] [21].This study, besides furthering the understanding of humpback whale distributions, presents opportunistic surveys as an effective method of continuous population monitoring, aiding in conservation efforts and species management.

Study Area and Sampling Units
The Maui 4-island region, which extends from 21˚0'N to 21˚35'N and 156˚50'W to 156˚25'W, is semi-enclosed by the islands of Maui, Moloka'i, Lana'i, and Kaho'olawe and located within the Hawaiian Islands Humpback Whale National Marine Sanctuary.The water depths within this area range from 1 -270 m, and the bottom topography consists of drowned reef features and sandy, flat concentric basins [22].
The study area was divided into 1km x 1km grid cells that covered 654 km 2   (Figure 1).The majority (56%) of water depth in the area ranged between 50 and 120 m, with sea surface temperature remaining fairly constant throughout the season at 25˚C -26˚C.

Data Collection
Humpback whale sighting data were collected from multiple whalewatching trips departing from both Ma'alaea and Lahaina Harbors daily between December 1, 2013 and April 30, 2017.Each whalewatching trip lasted 2 hours and followed a non-systematic survey design.Vessel speeds did not exceed 15 knots and whalewatches were conducted in Beaufort sea states of 6 or less.
Scanning for whales was completed by two trained observers and began once boats departed the harbor.Scanning was conducted by naked eye and with binoculars looking for visual cues of humpback whale presence.Survey effort (GPS track) and sighting details were recorded using the Whale and Dolphin Tracker (WDT) application [23].When a whale was sighted, the vessel approached the pod to a distance ≥ 91 meters (100 yards), in accordance with federal regulations.Upon approach and subsequent observation, information on group size, composition, and behavior were recorded.To ensure accurate group composition and size estimates, only sightings where group composition could be confirmed were included in subsequent analysis.Encounter location (latitude and longitude) was recorded using WDT when the vessel was ≤150 meters from the focal pod, as this distance provided sufficient geospatial resolution for the subsequent analysis and remained outside the 100 yard approach limit.

Data Analysis
Whalewatches are non-systematic surveys that do not incorporate transect lines, equal coverage of the sampling area, or distance sampling techniques.As such, a full distance sampling analysis generally utilized for cetacean surveys is not appropriate for this type of dataset.Instead, the authors utilize occupancy modelling which allows for the correction of biases inherent in opportunistic data collection, without requiring distance sampling techniques.Occupancy modelling relies on two components: site occupancy ( i ψ ), which is the probability that a species occupies the sampling site i, and is modelled as a function of site-specific covariates; and detectability ( ij p ), which is the probability that a species is de- tected during survey j at site i, and is modelled as a function of either site-specific or survey-specific covariates.Each are estimated using the logit-link function [24] [25]: ( ) Occupancy logit ln 1 Detection Probability logit ln 1 , where ln is the natural logarithm, x represents site-specific covariates, y represents survey-specific covariates, α and β are the estimated regression coefficients.
Although rarely used for cetacean surveys, occupancy models, with a slight Open Journal of Marine Science adjustment in interpretation, can be used as a viable option for analyzing habitat preferences of humpback whales using sighting data collected on WDT.The adjustments to parameter interpretation and details on how assumptions for occupancy modelling were met in the context of this dataset are detailed below.
Similar to the approach presented in previous studies [25], this study examined humpback whale use as opposed to occupancy, and the assumption of geographic closure is not met.To account for this we interpret model output as the probability of humpback whale use, as opposed to occupancy, which in conjunction with the whales' ability to move in and out of the survey sites allows us to relax the assumption of geographic closure [24].
The likelihood of false identifications of humpback whales is minimal given the size and nature of humpback whales.Additionally, there are no other mysticete whales that are commonly found in the Maui 4-island region.To ensure that the survey adhered to the assumption of no false identifications, only sightings where the whalewatching vessel approached the pod and positively confirmed the species as humpback whales were used.
To ensure detection of a humpback whale at one site was independent of detection at another, the final data were randomly subset to include only a single whalewatch trip from each harbor for each survey day.This subset was then utilized in subsequent analysis and the distance between the harbors ensured a humpback whale was not detected at multiple sites during a survey day.
To ensure there was no unmodelled heterogeneity, site and survey covariates known to impact humpback whale detection and occupancy were recorded in the field or calculated post survey and detailed below:

Site Covariates
Each site was characterized by the following variables: water depth, slope, distance to shore, latitude, longitude, and relative bottom terrain ruggedness.Mean site water depth was expressed in meters taken from Hawaii Mapping Research Group (HMRG) Bathymetry Synthesis data set (Resolution: 50 m) [26].The mean slope at each site was expressed in degrees (0˚ -90˚) and determined using the slope tool in ArcMap [27].Distance from shore was expressed as positive value in kilometers and determined using the near tool in ArcMap [27] to measure the distance from the center of each site to the nearest shoreline.Relative bottom terrain ruggedness [28] was calculated from bathymetry dataset [26] in ArcMap [27].Site covariates were tested for pairwise correlations using the stats package in R [29].To account for non-normality in site covariates, the Spearman correlation coefficient (r s ) was used to assess correlations.If site covariates were highly correlated (r s ≥ 0.7) [30], the most biologically relevant variable was retained.

Survey Covariates
Potential effects on detection included month within the survey season (Janu-Open Journal of Marine Science ary-April), survey effort/grid cell expressed as kilometers, and year.In this study, the effects of pod size and whale behavior were not accounted for since these variables could not readily be incorporated into the model.

Data Preparation
The data set was divided into two categories and detection histories created for each: one for pods without a calf present and one for pods with a calf present for each month during whale season (December, January, February, March, and April) from 2013-2017.For adult pods, a grid cell was given a "1" for detection when the cell was only occupied by pods that did not contain calves throughout the survey period.Sites that were sampled, but did not detect any pods, or detected pods with calves were given a "0".For calf pods, a grid cell was given a "1" for detection when the cell was occupied at least once by pod(s) that contained a calf throughout the survey period.A "0" was assigned to sites that were sampled, but did not detect any pods, or detected pods that did not contain a calf.
Data were analyzed using the unmarked package [31] in R to fit dynamic occupancy models [32] to humpback whale presence absence data.Two models were developed, one for pods without a calf present and one for pods with calves present (here after referred to as "adult model" and "calf model" respectively).In these analyses, data on "pod" presence and absence within a particular site are used to draw conclusions on humpback whale habitat preference for pods with and without a calf.The final adult and calf model predictions should be interpreted as the probability of humpback whale adult pod and humpback whale calf pod use within the Maui 4-island region.
Akaike Information Criterion (AIC) values [33] of fitted models were used to assess candidate models and obtained according to the following formula: where L is the maximum log-likelihood and K represents the number of parameters in the fitted model.The change in AIC (ΔAIC) was used to rank the candidate models and calculated using the following formula: To test model fit, the MacKenzie and Bailey Goodness-of-fit Test in R (mb.gof.test) using 10,000 bootstrap samples was run and c-hat calculated to Open Journal of Marine Science check for over dispersion [32] [34].

Results
From December 2013 to April 2017 a total of 1089 whalewatching trips were completed accounting 13,878 nautical miles of survey effort.During this time period humpback whale pods without calves were detected 786 times in 309 of the 654 grid cells sampled, while pods with calves were detected 682 times in 256 of the 654 grid cells sampled.Over the 4 months of observation, detections were highest in January for pods without calves (n = 250) and March for pods with calves (n = 201) (Table 1).
The most supported model of detection probability for the calf model and the adult model both included effects of survey effort and month (Table 2 and Table 3).Including a yearly effect in detection probability significantly improved model fit for the adult model only (Table 2 and Table 3).
The probability of detection for both the adult model and calf model increased with effort (Table 4).Progression from early whale season (December) to late whale season (April) saw an increase in probability of detection for the calf model and a decrease in probability of detection for the adult model (Table 4; Figure 2):   Correlation analysis found depth and distance from shore as well as latitude and longitude to be correlated, with values of 0.76 and 0.73 respectively.Between depth and distance from shore, depth was retained as the most biologically relevant parameter given that there is no clear on-shore off-shore gradient, and the surrounding 4 islands can confound this measurement with no easily explained proxy for distance from shore relating to humpback whale use within this particular study area.Furthermore, the bathymetry within the Maui 4-island region is complex with areas or large relatively uniform depth as well as steep gradients and depth contours which may impact how adult and calf humpback whale pods utilize the region.Latitude was favored over longitude, given the larger north-south gradient within the study area.Once correlated variables were removed, the following site covariates remained: depth, slope, roughness, and latitude.Depth ranked at the top for both adult and calf models (Table 5 and Table  6).Including only a depth covariate resulted in the best fitting model for pods without calves (Table 5).
The top ranked model for pods with calves included site covariates for depth,  slope, and latitude (Table 6).However, including slope only marginally reduced the AIC value (ΔAIC = 0.43).To ensure the most parsimonious model was used, the final model selected included covariates for depth and latitude only.
There was no evidence of over dispersion for the best adult and calf models ( ĉ < 1).The Mackenzie-Bailey (2004) chi-square test found both models predicted the data well (calf model p-value = 0.58; adult model p-value = 0.48).
Depth was the most important variable for determining humpback whale use, appearing in all of the top ranking models and being absent from the lowest ranking models (Table 5 and Table 6).Depth showed a significant positive relationship with humpback whale use for the adult model and a significant negative relationship with humpback whale use for the calf model (Table 7).The probability of use increased with increasing depth for adult pods, with the highest probability occurring in waters with a depth of 75 meters or more (Figure 3): The reverse was observed for calf pods, with the probability of use remaining high at depths of <100 meters, after which the probability of use drops quickly (Figure 4).In addition to depth, calf pods were also influenced by latitude, with the probability of use remaining high at lower latitude (<20.7˚N) for calf pods, after which probability of use decreases (Figure 5):

Discussion
This study compares the spatial distribution of humpback whale pods with and without calves over a 4-year period throughout the winter breeding season (December to April).The Maui 4-island region is a critical breeding habitat for humpback whales belonging to the Central North Pacific stock.This study, besides furthering the understanding of humpback whale distribution and habitat use in this region, presents opportunistic surveys as an effective method of continuous population monitoring aiding in conservation efforts and species management.

Depth and Latitude Effects
The results from this study provide strong evidence that the typical patterns of distribution according to age-class and group composition that are seen in other humpback whale breeding areas are also evident in the Maui 4-island region.This study agrees with previous work, both in Hawaii and in other breeding grounds, that mother-calf humpback whale pods demonstrate a habitat preference for shallow, nearshore waters.Our results show that there is an inverse relationship in the probability of area use between pods with a calf and pods without a calf present.Pods containing a calf were likely to be found in shallower waters, while pods containing no calves were likely to be found in deeper waters.The data collected had high geographic Open Journal of Marine Science and temporal coverage, which lead to strong confidence in the predicted results.
The effect of latitude suggests that pods with a calf present were most likely to use the southern portion of the survey area, specifically areas to the south of Ma'alaea Bay.This southern area has similar depth characteristics as other habitat areas along the northwest coast of Maui but may represent a region which is more protected from rough weather and anthropogenic disturbances.Rough weather could include wind or swells that might complicate nursing or increase energetic demands of swimming [9] [35] [36].Northeast trade winds dominate weather patterns in Maui and throughout the Main Hawaiian Islands [37], and a north-northwest swell is predominant during the winter months while the humpback whales utilize the Hawaiian Islands [38].Within the Maui 4-island region, this southern area may represent one of the most protected habitats in terms of wind and swell due to its southeast location in the lee of Maui and Haleakala volcano.This observed trend may also represent an avoidance strategy of commercial shipping traffic, which generally use the larger channels (Au'au, Kealaikahiki, Alalakeiki) and bypass south Maui [39].
There are no long term studies examining the habitat preference of mother-calf pairs in Maui leeward waters.Previous studies have reported mother-calf pairs preferring nearshore waters [40], locations < 0.5 km of shoreline along west Maui [41], and shallower waters [16].Previous research reported contradicting trends to previous studies [18]; however, that study was limited both geographically and temporally (examining only a portion of the Au'au Channel and only in March).Our model is strengthened by our large data set: four months (December to April) of daily whale sightings spread out over a large survey area.The model predicted pods containing a calf to occur just inside of the bathymetry gradient occurring in the Au'au Channel.This aligns with the finding that depth is a significant predictor for use of a grid cell by both pods with a calf and without a calf.

Seasonal Effects
Humpback whales tend to migrate based on age and sex classes [42], with sub-adults being first, then mature adults and mother-calf pairs being the final group to migrate [40] [43].Similarly, the composition of pods sighted in this study varied by month.The predominant group type was non-calf pods in December and January, while those with calves were predominantly observed in February, March, and April.The continued long-term, cross-seasonal monitoring presented here will allow for the determination of trends and for monitoring of potential changes and seasonal shifts.

Using Opportunistic Data
The use of opportunistic data is becoming increasingly important as an inexpensive, accessible method of cetacean monitoring.WDT is no exception, having recorded over 35,000 georeferenced sightings of 11 different species at no addi-Open Journal of Marine Science tional cost for survey effort.The ability to track large scale changes in marine environments can be aided by opportunistic data, such as those collected using WDT.However, it is important to note the use of opportunistic platforms can potentially introduce sources of variation and bias which need to be accounted for.Our occupancy models to predict humpback whale use from WDT data aligns with findings from previous research showing distinct area use of mother-calf pods within the Maui 4-island region [18] [40] [41].The standardization of data collection methodologies implemented with WDT and the use of occupancy models to account for biases in detection helps reduce the noise inherent in opportunistic data collection.The results presented here represent the first use of occupancy models on opportunistic data and showcase the potential use of such data to make ecological and management inferences.

Implications
Understanding patterns of habitat use and distribution among differing social groups of humpback whales identifies areas where potential anthropogenic effects could impact these animals.The preference of mothers with calves for shallower waters and nearshore areas overlaps with areas that are frequently used for fishing, whalewatching, and other recreational activities.Furthering our understanding of habitat preference and use, in conjunction with continued long-term monitoring as presented here, is of critical importance to developing biologically relevant conservation initiatives and management plans.Minimizing disturbance to mother-calf pods will likely contribute to the continued recovery of the Hawaiian population segment following its recent delisting from the United States Endangered Species Act.This study region largely comprises the Hawaiian Islands Humpback Whale National Marine Sanctuary, which has some restrictions in place to manage this population of whales: a minimum approach distance by water and air and restrictions on certain types of high-speed vessels (parasailing and jet-skis) between December and April.A greater understanding of habitat preference by mother-calf pairs in this important breeding ground could lead to informed management practices, including but not limited to, timed area closures or speed restrictions in the areas and times when calves are more likely to be present.

Figure 1 .
Figure 1.Map of study area showing bathymetry and 1 km × 1 km grid cells sampled from December 1, 2013 to April 30, 2017.
The most parsimonious model, having the lowest AIC value, was selected as the best fit model.Initial models included all sampling and site covariates, which were further subset and in different combinations based on model AIC values and weights (AIC wt ).AIC wt was calculated for model m, using the following formula:

Figure 2 .
Figure 2. Relationship of predicated detection probability with month from the adult (blue) and calf (red) models for humpback whales detected in the Maui 4-island.The error bars represent the 95% confidence interval.

Figure 3 .
Figure 3. Relationship of probability of whale use with varying depth for humpback whale pods without calves predicted from the adult model.The blue band represents the 95% confidence interval.

Figure 4 .
Figure 4. Relationship of probability of whale use with varying depth for humpback whale pods with calves predicted from the calf model.The blue band represents the 95% confidence interval.

Figure 5 .
Figure 5. Relationship of probability of whale use with varying latitude for humpback whale pods with calves predicted from the calf model.The blue band represents the 95% confidence interval.

Figure 6 .
Figure 6.Predicted probability of whale use for pods without calves (a) and pods with calves (b).Note: Points represent centroid of predicted site (1 km × 1 km) and were utilized to allow for visualization of varying use with bathymetry.

Table 1 .
Detection histories of humpback whale pods with and without calves sighted within a 1 km × 1 km site in the Maui 4-island region, Hawaii from December 2013 to January 2017.

Table 2 .
Model selection results of detection probability for humpback whale pods without calves using the Maui 4-island region from December 2013 to January 2017.
Note: Table values represent model AIC, the change in AIC for the best model (ΔAIC), the Akaike weights (AIC wt ), and number of parameters (K).

Table 3 .
Model selection results of detection probability for humpback whale pods with calves using the Maui 4-island region from December 2013 to January 2017.

Table 4 .
Untransformed estimates and standard errors of detection covariates for constant habitat model for humpback whale pods with and without calves in the Maui 4-island region, Hawaii.

Table 5
. AIC values for multi-season occupancy models describing the influence of environmental factors on the occurrence of humpback whale pods without calves in the Maui 4-island region from December 2013 to January 2017.

Table 6 .
AIC values for multi-season occupancy models describing the influence of environmental factors on the occurrence of humpback whale pods with calves in the Maui 4-island region from December 2013 to January 2017.

Table 7 .
Untransformed estimates and standard errors of site covariates for the probability of humpback whale use, after accounting for imperfect detections, for pods with and without calves in the Maui 4-island region, Hawaii.