Energetic Effects of Whale-Watching Boats on Humpback Whales on a Breeding Ground

Interactions between whale-watching boats and cetaceans can lead to changes in their behavior. From a management perspective, it is important to understand how this type of disturbance can be translated into physiological effects, such as changes in their energetic metabolism. Humpback whales (Megaptera novaeangliae) typically do not feed while in breeding grounds, thus they depend on finite energy reserves. The effect of whale-watching boats on the energetic metabolism of humpback whales, in the breeding ground of northern Peru (4°10′35″S, 81°08′03″W) was evaluated. Groups of humpback whales were tracked from land, under the following scenarios: with, without, and before-during-after the presence of whale-watching boats. Mass-specific cost of transport (COT) was used as a proxy of energetic efficiency and calculated from swimming speed and breath frequency estimations. No differences were detected in breath frequency, swimming speed, and COT when comparing whales with and without boats. However, in the presence of boats, swim speed increased, and COT decreased as the number of boats increased. Exponential increment in breathing frequency at higher swimming speed was not detected. The absence of swimming speeds beyond the assumed optimal range suggested no shifts into metabolic inefficiency. Our results suggest optimal swimming speed between 2 and 4.05 m/s, representing COT values between 0.020 and 0.041 J × (kg × m)–1. In light of our results, we encourage the implementation of regulations of the activity, particularly limiting the number of boats interacting with the same group of humpback whales.


INTRODUCTION
Whale watching, the observation of dolphins and whales in nature, is a growing economic activity in oceanic and coastal waters in many regions of the world. As whale watching grows, several studies have demonstrated the negative consequences of this activity on the behavior of cetacean species (reviewed in Senigaglia et al., 2016). Effects have been reported for small and large cetacean species including alterations in swimming speed, direction, breathing frequency, and overall behavior (e.g., Noren et al., 2009;Christiansen et al., 2010;Stamation et al., 2010;Senigaglia et al., 2016;García-Cegarra et al., 2019). Studies have linked those behavioral changes into effects on energy budget and metabolism of the species (e.g., Williams et al., 2006;Christiansen et al., 2013Christiansen et al., , 2014a. Physiological responses are essential for a better understanding of organismal and population consequences of the disturbance caused by whale-watching boats (Costa, 2012;New et al., 2015;Pirotta et al., 2018).
Mating, breeding, and migration are highly energydemanding activities for cetacean species. Southern right whales (Eubalaena australis) lose on average of 25% of their body volume during the breeding season . To minimize the rate of decline in body condition and optimizing calf growth during migration to their feeding ground, lactating humpback whale females reduce their metabolic rate to half that of adults in foraging grounds Nielsen et al., 2019). During breeding, fin whales (Balaenoptera physalus) consume between 19 and 26% of their energy reserves (Lockyer, 1981(Lockyer, , 1984. Significant reductions in net energy intake and/or increases in energy expenditure can lead to changes in body condition (Frid and Dill, 2002). Fetal growth (Christiansen et al., 2014b) and calf body condition  of mysticeti whales is significantly determined by the body condition of the maternal female. Energetic consequences of behavioral changes could lead to long-term reductions in body condition, reproductive success (fitness), leading to negative population consequences . Whale watching can disturb critical behaviors such as lactation thus, reducing the energy transfer from the mother to the calf, affecting growth rates, which can have consequences in migratory timing and heat loss. Overall, these effects may lead to negative consequences for the long-term individual survival, reproduction success, and recruitment into the population (Bejder, 2005;Lusseau, 2006;Nowacek et al., 2016).
Studying the physiology of large sized, free-ranging cetaceans is difficult due to the methodological and logistical constraints. However, based on indirect estimations several studies have showed how behavioral changes can translate into physiological variability (Christiansen and Lusseau, 2015;Pirotta et al., 2018). For example, a 23.2% increase in energy expenditure was detected for traveling minke whales (Balaenoptera acutorostrata) due to increasing breath frequencies during interactions with whalewatching boats (Christiansen et al., 2014a). Energy intake of minke whales and killer whales (Orcinus orca) was reduced by 42 and 18% because approaching whale-watching boats induced a reduction of feeding times (Williams et al., 2006;Christiansen et al., 2013). So far, no changes have been reported in the energetic budget of humpback whales (Megaptera novaeangliae) in feeding grounds during whale watching (Di Clemente et al., 2018). Understanding behavior variability and underlying changes in physiology can provide knowledge to establish links between short and long-term consequences of disturbance on species of cetaceans (New et al., 2015;Pirotta et al., 2018).
The humpback whale is one of the most popular species for whale watching (O'Connor et al., 2009). A suite of negative effects in response to whale-watching boats, such as the alteration of surfacing and diving behavior, aerial activity, acoustic behavior, and swimming speed (e.g., Corkeron, 1995;Scheidat et al., 2004;Sousa-Lima and Clark, 2008;Stamation et al., 2010;García-Cegarra et al., 2019) have been documented on the species. Humpback whales migrate from polar and temperate feeding grounds to tropical and subtropical breeding grounds (Dawbin, 1966). In the latter, humpbacks whales typically do not feed (Chittleborough, 1965), and depend on the energetic reserves stored in the blubber layer, muscles, visceral organs, and bones (Nordøy et al., 1995;Gunnlaugsson et al., 2020) acquired during the feeding season. Although feeding events have been observed in breeding and migratory routes (Stamation et al., 2007;Frisch-Jordán et al., 2019;, energy intake is limited or absent during the breeding season. Here, we investigated the effect of whale-watching boats on the energy consumption and efficiency of humpback whales in the coast of northern Peru, Southeast Pacific. This population is also known as stock G (IWC, 1998), and its breeding ground ranges from northern Peru to Costa Rica (Scheidat et al., 2000;Acevedo et al., 2017;Valdivia et al., 2017) and possibly up to Nicaragua . Respiratory rates have been used to estimate the mass-specific cost of transport (COT), for large animals under undisturbed and disturbed conditions (Williams and Noren, 2009;Langman et al., 2012;Christiansen et al., 2014a;Maresh et al., 2014). This metric includes the mass and transport as a relative measure of the metabolic rate. COT can provide insights about the overall energy expenditure during movements, and how efficiently the energy is used relative to the distance traveled (Tucker, 1970). An optimal range of transport can then be estimated based on the range of swimming speeds where the COT is reduced to its minimum. Thus, variation in COT could be a useful proxy to understand the physiological effects of whale watching on cetaceans. We predicted that approaching whale-watching boats to whales would lead to an increase in energy consumption by increasing swimming speeds and breath frequencies.

Study Area
Daily land-based surveys were performed from a rocky cliff named "Cerro La Mesa" at 31 m above the sea level. From this position, we covered a panoramic view of 7 km of radius (153km 2 of area) of the coastal area between El Ñuro (4 • 13 01 S, 81 • 10 35 W) and Máncora (4 • 06 26 S, 81 • 02 50 W), where whale-watching activities are conducted (Figure 1). Observation of whales and boats were performed daily, between 07:00 and 10:00 from July 17th to August 29th and from September 5th to October 15th in 2016. During the study period, weather conditions were favorable during 95% of the survey days with visibility of 6 to10 km. Sea conditions allowed the tracking of whales throughout the study period with Beaufort states between 0 (34%) and 1 (66%). The area has been described as the Southern limit of the breeding ground of Stock G with the season ranging from mid-July to the end of October (Pacheco et al., 2009;Guidino et al., 2014).

Sampling Methods
Groups of humpback whales were recorded by continuous focal group follows (Altmann, 1974) using a digital theodolite (Nikon NPL-322, Nikon Trimble, Tokyo, Japan) and 15×50 Nikon binoculars. A group was defined as individuals of whales within 100 m distance of each other, moving in the same direction and displaying almost synchronized diving and movement patterns (Whitehead, 1983;Mobley and Herman, 1985). Calves were identified as whales with sizes ranging from one to two-thirds of the size to their accompanying adult, assumed to be the mother (Herman and Antinoja, 1977;Mobley and Herman, 1985). Groups of humpback whales were classified into two main categories: groups with calves and groups without calves. Mothercalf and mother-calf with one or more escorts were considered as groups with calves, while groups including only adults or subadults were considered as groups without calves (Supplementary Table 1 provides a description of all group categories). Tracking of focal groups started when both observers (DV and AGC) spotted humpback whales with binoculars. Once the group was spotted and followed unequivocally for more than 10 min, the group was chosen as a focal group for tracking. The type of group and number of individuals was defined at the initial sighting and confirmed during the tracking. If the focal group splitted, one of the groups was chosen to continue tracking. When groups merged, the tracking continued for the focal group. However, this type of variability in the dynamics of the groups was not included in the analysis of this study.
Horizontal and vertical angles were recorded, based on a georeferenced reference point, for each emergence of the focal group using the digital theodolite. Using the sinus theorem, elevation, and coordinates of the reference point, the angles were converted into geographic coordinates and geographic tracks for each focal group. The specific details of the methodology are provided in García-Cegarra et al. (2019) including the geographic position of the tracks. A preliminary analysis of track accuracy was carried out by tracking a moving boat and recording its position every 5 min. A measurement error of 35 m for distances >4.5 km was estimated (Romero, 2015) and used to correct all geographic calculations. For groups with more than one adult, the first adult emerging to the surface was considered as reference for the overall group movement. The total distance traveled by each group was calculated as the sum of distances between each recorded surfacing location. Simultaneously, the number of breaths of each whale, the number of boats (when present), and the time surfacing of the whale group were recorded.
The distance between focal groups and whale-watching boats was measured using the digital theodolite, following the same method as for the focal groups. We considered an interaction threshold of 400 m distance between whales and boats following Baker and Herman (1989). However, behavioral changes can occur before boats enter this threshold distance (e.g., Watkins, 1986;Baker and Herman, 1989;Sprogis et al., 2020a). Sprogis et al. (2020b) demonstrated that vessel noise drives the behavioral response of humpback whales to boats. Although, we recorded the number of boats and their timing with the whales, the limitations imposed by the distance of the land-based survey, precluded us to unravel the role of noise of the boats on the behavior of the whales.
To measure the swimming speed and breath frequency of whales in the presence and absence of whale-watching boats, a minimum of 10 min of interaction was determined as lower time threshold. This value was determined to increase the precision of the measurements and to reduce the amount of interaction events inadequate for estimations. For example, for whales emerging every 8-9 min, the presence of a whale-watching boat for less than 10 min would only allow to register one emergence. This would only provide one record of breath thus traveled distance, and swimming speed cannot be estimated. Even though, whale responses to approaching boats may occur in a shorter period of time (<10 min), we assumed that these short burst events are unlikely to have an important impact on the overall energy efficiency. To avoid underestimation and bias in traveled distance, swimming speed and breath frequency estimations, only relative linear tracks were used in the analysis (Christiansen et al., 2014a). Tracks with whales logging at the surface were also removed from the analysis. Based on the presence and absence of whale-watching boats, focal groups were divided into three independent scenarios: before, during, and after the presence of whale-watching boats. Data from groups tracked beforeduring, during-after or before-during-after the presence of boats were divided into two and/or three dependent observations, respectively. For example, a before-during track was divided into two data observations, one going to the before scenario and the other to the during scenario. Tracks with the boat always present or absent provided one measurement each, while tracks before-during; during-after or before-during-after, provided 2 or 3 measurements, respectively.
The average swimming speed (S) of the focal groups was calculated by dividing the sum of the distances traveled between all emergence positions by the total time of the track. Breath frequencies (f ) of humpback whales was calculated by dividing the total number of breaths per individual by the total time of the track. For groups with two or more adult whales, the total number of breaths was recorded for the group and finally divided by the number of adults to obtain the individual breath frequency. For groups with calves, breath frequencies and swimming speeds were calculated based only on the adult's behavior. This study was carried out under the approval of the Comité de Ética de Investigación Científica de la Universidad de Antofagasta, Chile (CEIC REV N • 039/2017 and 7298/2015).

Mass-Specific Cost of Transport Estimation
The calculated swimming speed and breath frequency were used to calculate the mass-specific cost of transport (COT) following the methodology described in Christiansen et al. (2014a). Supplementary Table 2 shows a summary of the parameters, equations, theoretical values, and references used for the estimation of COT in humpback whales.
The respiratory volume per minute, V min (l × min −1 , the volume of air breathed per minute) was calculated from the measured breath frequency, f (breaths × min −1 ) and the tidal volume, V t (l × breath −1 , volume of air inhaled per respiratory cycle. The tidal volume for humpback whale was assumed to be 60% of the volume or lung capacity, V c (Wahrenbrock et al., 1974;Blix and Folkow, 1995).
The lung capacity (V c ) was estimated from the body mass of the humpback whale (Dolphin, 1987).
The body mass, W (kg), was estimated from the average length (Lockyer, 1976). L refers to the length of the whale in meters. The average size (12.2 m) for both sexes based on humpback whale catch data of Peru, between 1961 and 1966, (Ramírez, 1988) was used for this purpose, resulting in average weight for individuals of both sexes of 25 317 kg.
The oxygen consumption rate, V O2 (l O2 × min −1 , volume of oxygen consumed per minute) was estimated from the respiratory volume per minute (V t ), the oxygen concentration in the air, P O2 (0.21, ratio of O 2 in the inspired air) and the oxygen extraction rate, E O2 . The oxygen extraction rate from the air in the lungs during a respiratory cycle was established at a value of 45% (Wahrenbrock et al., 1974;Blix and Folkow, 1995).
The metabolic rate, MR (J × min −1 ), energy consumed per minute, was estimated by transforming the oxygen consumption rate V O2 , into energy units using the conversion factor 20.1 kJ × l O2 −1 .
Only adult whales were used in the analyses, to avoid confounding effects of growth on the metabolic rates. COT for calf groups was estimated from the breath frequency and Frontiers in Marine Science | www.frontiersin.org -Nielsen, 1972) showing the relationship between speed and cost of transport based on energetic efficiency. Red dashed lines represent the theoretic division of the curve based on metabolic efficiency. Increasing swimming speed, increase the energetic efficiency by reducing the COT. This reduction continues until the optimal range of transport, where the maximum efficiency and minimum COT are reached. Finally, when speeds continue to increase, COT increases rapidly leading to an increasing energetic inefficiency.
swimming speed of the adult individual in each group (e.g., the mother). Even though this data would not reflect newborn metabolism, we assumed that the changes in behavior and metabolism of the adult/s in a group with calves could be an indirect proxy of the stress that the calf faces during whale watching.

Effect of Whale-Watching Boats
We hypothesized that breath frequency and swimming speed will increase while COT decreases between before and during scenarios and return to normal conditions after the boats have left. Higher number of boats and longer times spent with the whale groups will increase the mentioned effect in the three variables during the activity. To test the differences in the breath frequency, swimming speed and mass-specific cost of transport, three regression models were developed using an estimation procedure based on generalizing estimating equations for the generalized linear models (glm) and for the normal model. Normality of the variables was obtained after a log transformation for breath frequency and swimming speed and a logit transformation for COT. Observations of the same groups of whales (e.g., Before-During) were not independent, thus an exchangeable correlation structure was used for these observations. Each model tested the effects of the scenario (before, during, after) and calf presence on the breath frequency, swimming speed and COT. We analyzed the effect of the presence of calf on each of the models as several studies have already reported significant differences in behavior between groups with and without calves (Corkeron, 1995;Stamation et al., 2010;Craig et al., 2014;Sprogis et al., 2020a). The effect of the number of boats and time spent with the group of whales in the three variables was tested for the during and after scenarios. The estimations were performed using the geepack (Højsgaard et al., 2006) package in R (R Core Team, 2020) version 4.0.2.

Energy Efficiency
Non-linear regressions were performed between COT and swimming speeds values to determine whether these variables follow a U-shaped curvilinear relationship (Schmidt- Nielsen, 1972) as it has been documented for other marine mammals species (Figure 2, e.g., Williams et al., 1993;Otani et al., 2001;Rosen and Trites, 2002;Williams and Noren, 2009). The curve indicates that optimal speed ranges of transport can be detected when COT reaches its minimum values due to higher efficiency in the consumption of energy for displacement. However, when swimming speed increase beyond this range, the breath frequency increases disproportionally due to an increasing oxygen demand. This drastically reduces the energetic metabolism efficiency during transport and can be observed as an inflection point followed by an exponential increase of the COT values. Similarly, when analyzing the link between swimming speed and breath frequency, the onset of energetic inefficiency can be observed by an exponential increase of the breath frequency when speed exceed its optimal range (Williams et al., 1993;Yazdi et al., 1999). Linear regressions between swimming speed and breath frequency were performed to determine the relationship between the two variables. All tests and statistical analyses were performed using R version 4.0.2.

RESULTS
A total of 167.4 h of survey were completed. 132 h were used to visually follow focal groups; 39% (51.5 h) without and 61% (80.5 h) with whale-watching boats. A total of 412 humpback whales were followed in 173 focal groups: 91 groups without calves and 82 groups with calves ( Table 1). From the 173 focal groups, 67 and 48 whale groups were tracked with boats always absent and absent, respectively. Also, 29, 20, and 9 groups were tracked before and during; during and after; and before, during, and after the presence of boats, respectively. Dividing groups into their respective encounter scenarios and adding the tracks where boats were always present (during) or absent (before), a total of 105, 106, and 29 breath frequency, swimming speed and COT measurements were, respectively, obtained for the encounter scenarios before, during, and after ( Table 1). A mean of 3.13 (SD = 1.52; range: 1-9) and a maximum of 9 whale-watching boats were observed with the focal groups. On average boats followed groups of whales during 48.32 min (SD = 25.15 min, range: 13-125 min) and keeping an average minimum distance of 39.26 m (SD = 47.21 m, range: 6-125 m). Significant differences between groups with calves and without calves were found for breath frequency and swimming speed. Based on significant coefficients, groups without calves registered breath frequencies and swimming speeds 18.43 and 25.78% higher, respectively, than groups with calves (Tables 2, 3).

Effects of Whale-Watching Boats
Breath frequency, swimming speed, and COT, did not change significantly between before and during the presence of whalewatching boats. Likewise, no significant differences in breath   frequency, swimming speed and COT were found between before and after, or between during and after ( Table 3). The number of boats present, and the duration of interactions had a significant effect on breath frequency, swim speed and COT. Each additional boat with the whale group (e.g., from 2 to 3 boats) led to a 7% increase of the median swimming speed and a 10% reduction of the median COT. Furthermore, boats spending more time with a whale group led to a reduction of the speed of the groups, as each additional minute spent with the whale groups represented a significant 0.5% decrease in the median of swimming speed. This effect was present even after the boat left, as each additional minute spent with the whale group lead to a significant 0.4% in the median of swimming speed (Tables 2, 3).

Energetic Efficiency
Groups with and without calves showed the same pattern; high COT values at low swim speeds (Figure 3). Increasing swimming speeds led to a significant decrease in the COT, following a power function in groups with calves ( Figure 3A, COT = 0.069 × S −0.95 , pseudoR 2 = 0.63) and without calves ( Figure 3B, COT = 0.085 × S −1.05 , pseudoR 2 = 0.53). Changes in swimming speed explained between 53 and 63% of the variation of the COT values. Minimum COT was registered between 2 and 3.2 m × s −1 for groups with calves and between 2 and 4.05 m × s −1 for groups without calves, determining massspecific costs of transport between 0.023-0.036 and 0.020-0.041 J × [kg × m] −1 , respectively. However, the true minimum COT values remain unknown as no inflection point, where COT is excepted to increase, was detected during the movement of groups with and without calves.
While the presence of whale-watching boats did not lead to significant increases in the breath frequency of humpback whales, breath frequency increased linearly with swimming speed (Figure 4 f= 0.57+0.07 × S; R 2 = 0.018; F-stat = 2.17; p-value = 0.03), at a rate of 0.07 breaths × min −1 for every m × s −1 increase in swim speed. However, breath frequency in humpback whale groups might be influenced by other factors, as only 1.8% (R 2 : 0.018) of its variability was explained by changes in swimming speed.

DISCUSSION
Behavioral responses to whale-watching boats can potentially lead to an increase in the metabolic rate in cetaceans due to an increase in breath frequency and swimming speed (e.g., Christiansen et al., 2014a). Our results suggest that the mere presence of whale-watching boats does not lead to significant behavioral changes, but as the number of boats increased, so did the breath frequency and swim speed of the whales. Energy expenditure during transport remains near optimal values. No increments in mass-specific cost of transport were recorded and breath frequencies continued to increase linearly within the range of measured swimming speed.
groups with calves (Carvalho-Gonçalves et al., 2018;Bejder et al., 2019). Mother and calf groups face physical and physiological challenges, because calves are learning social skills, have less muscle strength, and lung size compared to adults. This reduces their escape response capacity (e.g., fast swimming) when facing natural predators such as killer whales, harassment of competing adult males or whale-watching disturbance. In breeding grounds, calf groups inhabit shallow and calm waters to reduce disturbance by competitive adults (Smultea, 1994) and optimizes the energy transfer from the mother to the newborn during nursing Sullivan and Cartwright, 2009;Videsen et al., 2017). Our minimal COT estimations for a 25,317 kg adult humpback whale, ranges around 0.023 J × [kg × m] −1 for groups with calves and 0.02 J × [kg × m] −1 for groups without calves. We developed a log-linear regression between body mass and COT values, using published values for other marine mammals including our results. Our estimates fitted into this regression (Figure 5, COT = 4.14-0.420 × log(W), R 2 = 0.75, F-stat = 40.35, p = 3.653e −5 ), confirming that they occur within the expected range for an animal of 25,000 kg.

The Effects of Whale-Watching Boats
In contrast with results elsewhere (e.g., Corkeron, 1995;Scheidat et al., 2004;Schaffar et al., 2010;Stamation et al., 2010), the sole presence of whale-watching boats did not trigger changes in the behavior of humpback whales. However, each additional boat led to a significant 7% increase of the swimming speed, and consequently a 10% COT reduction. Similarly, a high number of whale-watching boats at less than 400 m induced killer whales to follow a more sinuous swimming path (Williams et al., 2002;Williams and Ashe, 2007) together with an increase in their swimming speed (Williams et al., 2002). Arguably a higher number of boats could result in a higher noise level. Higher noise levels can induce strong behavioral changes in humpback whales (Sprogis et al., 2020b). This can explain the differences in response of humpback in function of the numbers of boats observed in our study. Also, a high number of boats may lead to a lower degree of compliance with the voluntary guidelines for whale watching proposed for this region (Pacheco et al., 2011). Boats positioning themselves closer to the whales presumably to ensure the satisfaction of the tourist (García-Cegarra and . Although this study did not gather information on specific features of the boats such as the type of engines, we recognize that different boat-engine configurations can lead to different noise levels, which may finally translate into different levels of disturbance. The time that boats spent with whales led to significant decreases in the swimming speed of groups of whales during and after the presence of the boat. The latter effect may be related to whale-watching boats performance. Groups of whales swimming at high speed will usually be sighted by boats for less time because they are more difficult to follow from a tour boat. The data presented here fitted the first half of the typical U-shaped relationship between COT and swimming speeds. Increasing swimming speed, implies a reduction of breaths per kilometer traveled, hence reducing COT values and reaching a minimum when reaching the optimal range of transport. The second half of the U-shaped curve was not observed, as the recorded maximum swimming speeds (4.05 m × s −1 ) did not exceed the optimal range. Similarly, studies on harbor porpoises (Phocoena phocoena) (Otani et al., 2001), killer whales (Williams and Noren, 2009), and minke whales (Christiansen et al., 2014a) did not register the second half of the U-shaped curve. Conversely, Williams et al. (1993) and Yazdi et al. (1999) described the complete curve for bottlenose dolphins (Tursiops truncatus), showing that COT increases when reaching speeds beyond the optimal range. Both studies were performed by dolphins trained to reach specific speeds, in confined (e.g., dolphinarium) and open water (e.g., following boats). Our results suggest that such speeds may not be reached in natural conditions and individuals will tend to maintain energy efficiency during travel even when facing disturbance.
The optimal swimming speed estimated for humpback whales ranged between ca. 2 and the maximum speed recorded in our study, 4.05 m × s −1 , generating minimum COT values between 0.020 and 0.041 J × [kg × m] −1 . As no inflection point and subsequent increase of the COT were recorded, energetic efficiency was maintained during tracking. However, since the second half of the U-shaped curve was not detected, the upper limit of the optimal swimming speed range cannot be determined with accuracy. Records of humpback whales swimming at more than 4.05 m × s −1 would be needed to determine if the optimal range continues beyond our estimated range. Additionally, breath frequencies increased linearly with increasing swimming speed, confirming persistence of energetic efficiency (Williams et al., 1993) over the range of swimming speeds recorded. Similar results have been described for killer (Williams and Noren, 2009) and minke whales (Christiansen et al., 2014a). Performing a quadratic regression, a possible inflection and increase of COT was explored (Figure 6). An inflection point can be noted between 2.3 and 2.4 m × s −1 and an optimal swimming speed range between 2 and 2.6 m/s. However, the low number of groups of whales swimming at speeds exceeding this range (S>2.6 m/s), increases notably the standard error reducing prediction accuracy. The optimal range could extend further, or on the contrary higher swimming speed could directly lead to energetic inefficiency. Experiments carried in controlled conditions and with small cetaceans (Williams et al., 1993;Yazdi et al., 1999), allowed to exceed their optimal range. This performance is unlikely to occur in large cetaceans like the humpback whale in nature.
Cetaceans may respond to human disturbance as they do against natural predators (e.g., killer whales) (Christiansen and Lusseau, 2012). Some species of baleen whales maintain high and sustained speeds to avoid killer whale attacks (Ford et al., 2005;Ford and Reeves, 2008). When chased by orcas, minke whales can keep high velocities for several hours (ca. 8.5 h) over large distances (ca. 18 km) (Ford et al., 2005). Humpback whales may physical defense themselves when confronting predators. Mothers would defend their calf when facing attacks from killer whales (Pitman et al., 2017). Species of cetacean that fight predators, tend to be less hydrodynamic but with a better ability to maneuver and with robust bodies with callosities that could be used as weapons or amour (Ford and Reeves, 2008). However, they can also sustain high speed when fleeing from predators. In humpback whales, the presence of a single boat may not trigger an escape response, however, several boats may elicit a fast response.
Humpback whales of Breeding Stock G face other anthropogenic stressors such as entanglement with fishing gears, shipping noise, and vessel collision, throughout their breeding and feeding grounds (O'Connor et al., 2009). The effect of whale-watching interaction cannot be considered only as a punctual and occasional event, because such repeated anthropogenic stressor events may occur for the same individual or group in addition to natural events (e.g., escape from predation, intraspecific competition). Vulnerable groups, such as mother and calf groups (Stamation et al., 2010;García-Cegarra et al., 2019) move slowly and are usually found closer to the coast, being easily approached by several whale-watching boats (García-Cegarra et al., 2019). Whale watching is a growing industry in Peru (Guidino et al., 2020) but is not regulated yet (Pacheco et al., 2011). We urge the establishment of regulations particularly measurements that controls the time and number of boats per group of whales especially in mother and calf groups (García-Cegarra et al., 2019).

DATA AVAILABILITY STATEMENT
The minimum data to reproduce the analysis supporting the conclusions of this article is available in the Supplementary  Table 3. Further inquiries can be directed to the corresponding author/s.

ETHICS STATEMENT
The animal study was reviewed and approved by the Comité de Ética de Investigación Científica de la Universidad de Antofagasta, Chile (CEIC REV N • 039/2017 and 7298/2015).

AUTHOR CONTRIBUTIONS
DV and AP conceived the manuscript. AG-C provided most of the funds during the field works. DV and AG-C collected the data with the advice and guidance of AP. DV analyzed and interpreted the data and results of the study. DV led the preparation of the manuscript with the guidance of AP and AG-C. All authors edited the manuscript.

FUNDING
This study was funded by the Rufford Foundation via Rufford Small Grants for Nature Conservation (RSG: 15903-1). AG-C was also supported by a Ph.D. Scholarship from the Chilean National Commission for Scientific and Technological Research (CONICYT/63140172-2014).
Alexander Alburquerque for taking us every day to the landbased observation point. We very much appreciate the comments and corrections to this manuscript made by Nicola L. Ransome, Sebastian Uhlmann and two reviewers that help us to improve an early version of this manuscript.