Intra‐season variations in distribution and abundance of humpback whales in the West Antarctic Peninsula using cruise vessels as opportunistic platforms

Abstract Fine‐scale knowledge of spatiotemporal dynamics in cetacean distribution and abundance throughout the Western Antarctic Peninsula (WAP) is sparse yet essential for effective ecosystem‐based management (EBM). Cruise vessels were used as platforms of opportunity to collect data on the distribution and abundance of humpback whales (Megaptera novaeangliae) during the austral summer of 2019/2020 in a region that is also important for the Antarctic krill (Euphausia superba) fishery, to assess potential spatiotemporal interactions for future use in EBM. Data were analyzed using traditional design‐based line transect methodology and spatial density surface hurdle models fitted using a set of physical environmental covariates to estimate the abundance and distribution of whales in the area, and to describe their temporal dynamics. Our results indicate a rapid increase in humpback whale abundance in the Bransfield and Gerlache Straits through December, reaching a stable abundance by mid‐January. The distribution of humpback whales appeared to change from a patchier distribution in the northern Gerlache Strait to a significantly concentrated presence in the central Gerlache and southern Bransfield Straits, followed by a subsequent dispersion throughout the area. Abundance estimates agreed well with previous literature, increasing from approximately 7000 individuals in 2000 to a peak of 19,107 in 2020. Based on these estimates, we project a total krill consumption of between 1.4 and 3.7 million tons based on traditional and contemporary literature on per capita krill consumption of whales, respectively. When taken in the context of krill fishery catch data in the study area, we conclude that there is minimal spatiotemporal overlap between humpback whales and fishery activity during our study period of November–January. However, there is potential for significant interaction between the two later in the feeding season, but cetacean survey efforts need to be extended into late season in order to fully characterize this potential overlap.

Our results indicate a rapid increase in humpback whale abundance in the Bransfield and Gerlache Straits through December, reaching a stable abundance by mid-January.
The distribution of humpback whales appeared to change from a patchier distribution in the northern Gerlache Strait to a significantly concentrated presence in the central Gerlache and southern Bransfield Straits, followed by a subsequent dispersion throughout the area. Abundance estimates agreed well with previous literature, increasing from approximately 7000 individuals in 2000 to a peak of 19,107 in 2020.
Based on these estimates, we project a total krill consumption of between 1.4 and 3.7 million tons based on traditional and contemporary literature on per capita krill consumption of whales, respectively. When taken in the context of krill fishery catch data in the study area, we conclude that there is minimal spatiotemporal overlap between humpback whales and fishery activity during our study period of November-January. However, there is potential for significant interaction between the two later in the feeding season, but cetacean survey efforts need to be extended into late season in order to fully characterize this potential overlap.

| INTRODUC TI ON
As a result of overexploitation, fishing has shifted focus to areas beyond national jurisdictions (Kawaguchi & Nicol, 2007;Nicol et al., 2012). The fishery for Antarctic krill (Euphausia superba, henceforth krill) is the largest fishery in the Southern Ocean in terms of biomass (Nicol et al., 2012). Over the preceding three decades, its spatial footprint has contracted around the Western Antarctic Peninsula (WAP; Krüger, 2019;Santa Cruz et al., 2018), and in combination with a steady increase in catch levels that now represent the highest the fishery has ever seen, concerns about fisheries impacts on the broader marine ecosystem are being raised.
Currently, the fishery is managed by the Commission for the Conservation of Antarctic Marine Living Resources (CCAMLR) using a precautionary catch limit of 620,000 tons, distributed between subareas across the Drake Passage, Antarctic Peninsula, and Scotia Sea (Figure 1; Hill et al., 2016). Additional quotas are available in East Antarctic waters; however, the fishery has not operated there at commercially meaningful levels for several decades (CCAMLR Secretariat, 2021). Despite the fishery operating almost exclusively in CCAMLR management Area 48 (covering the WAP and Scotia Sea) since the turn of the millennia, two large-scale krill surveys in this area conducted two decades apart showed remarkably similar biomass estimates (ca. 60 million tons; Macaulay et al., 2019). Thus, given that the fishery currently takes less than 0.75% of the estimated biomass in area 48, there is increasing interest in expanding the quota beyond the current trigger level (Nicol et al., 2012).
In addition to pressures from fishing, rapid climate warming (Ducklow et al., 2013;Meredith & King, 2005) may potentially affect the spatial distribution, abundance, and recruitment of krill (Atkinson et al., 2004, which in turn impact krill-dependent predators in the WAP area (Krüger et al., 2020;Watters et al., 2020).
The interactive effects of climate-driven modification of the WAP on krill dynamics and various population trajectories of upper trophic predators are likely to be complex Trivelpiece et al., 2011). In light of this, CCAMLR is currently looking to implement ecosystem-based management (EBM) to adaptively manage catch limits and spatial allocations in a manner that avoids negatively impacting krill-dependent aspects of the ecosystem (SC-CAMLR, 2011). Implicit to this is the need for ecological information on krillconsuming predators.
For more than two centuries, pinnipeds and large cetaceans were subjected to uncontrolled harvesting (Trathan & Reid, 2009), with some species driven to near extinction (Tulloch et al., 2018). Surma et al. (2014) estimated that this reduction of large-bodied krill predators in the ecosystem may have resulted in the competitive release of krill, sensu the Krill Surplus Hypothesis (Laws et al., 1977).
Under this hypothesis, the standing stock of krill is estimated to have increased approximately 17%-25% in response to declines in krill predators such as seals and whales, allowing population increases of less competitive krill predators with shorter generation times K E Y W O R D S Antarctic Peninsula, density surface modeling, ecosystem interactions, humpback whale, platforms of opportunity, spatiotemporal variation

T A X O N O M Y C L A S S I F I C A T I O N
Spatial ecology; Trophic interactions F I G U R E 1 The study area for analysis, marked white, survey transects for the different cruises (C1 (25.11.19-12.12.19), C2 (16-27.12.29), and C3 (12-18.01.20)) and humpback whale observations (black dots) from all cruises. The Bransfield and Gerlache (indicated by arrow) Straits and Drake Passage are shown. Inlet map shows study area and tracklines in relation to the South American and Antarctic continents. Map produced using Quantarctica (Matsuoka et al., 2021) such as penguins (particularly gentoo Pygoscelis papua, adélie P. adeliae, and chinstrap penguins P. antarcticus) and pinnipeds (particularly crabeater seals Lobodon carcinophagus and Antarctic fur seals Arctocephalus gazella; Surma et al., 2014). Following the moratorium prohibiting the commercial catch of whales, set by the International Whaling Committee in 1986, many cetacean populations are showing signs of recovering toward their preharvest levels (Bettridge et al., 2015;Thomas et al., 2016;Tulloch et al., 2018;Zerbini et al., 2019). Consequently, krill consumption by cetaceans in the CCAMLR Area 48 is likely to have increased in the postwhaling era in line with rebounding populations. The krill consumption by cetaceans in area 48 during the austral summer season of 1999/2000 was estimated to 1.6-2.7 million tons (Reilly et al., 2004) but given potential and observed recovery rates of cetacean species the past four decades (Zerbini et al., 2019), this consumption estimate may be outdated.
Subsequently, recovering populations of these major krill consumers may exert significant top-down pressure on krill, with potential competitive consequences for other predators such as penguins, none of which are currently considered in a management context by CCAMLR.
The humpback whale (Megaptera novaeangliae) is a major krill consuming predator in the Southern Ocean. Most individuals migrate annually from low-latitude breeding grounds to high-latitude feeding grounds, where they are assumed to consume 83% (Lockyer et al., 1981) of their annual caloric intake during their summer foraging period. Contemporaneous estimates of up to 7000 humpback whales in the northern Antarctic peninsula during summer are now two decades old (Hedley et al., 2001), and were estimated to consume approximately 417,000-806,000 tons of krill biomass in the Antarctic Peninsula alone (Reilly et al., 2004). A recent study of humpback whale reproductive success indicated that numbers in the WAP may be increasing (Pallin et al., 2018), suggesting prolonged recovery from historical harvesting. This is in line with recent population assessments, which indicate that populations in both the western South Atlantic and the eastern South Pacific have recovered to ca. 90% of pre-exploitation levels (Johnston et al., 2020;Zerbini et al., 2019).
Cetacean abundance is traditionally estimated from line transect surveys using distance sampling as the general analytical framework (Buckland et al., 2005). During line transect surveys, whales are counted along systematic predetermined transects, arranged across the study region to minimize sampling bias in coverage. Distance sampling is then used to estimate a detection probability function that can be applied to estimate whale density (number km −2 or nm −2 ) and subsequently estimating the total abundance of whales in the area of interest by simply multiplying the density by the area.
The latter assumes that the survey covers representative areas of the entire region (Buckland et al., 2005) and is generally referred to as a "design-based" approach (Barry & Welsh, 2001). Density Surface Modelling (DSM) is a two-stage approach for estimating spatially varying density from distance-sampling data (Miller et al., 2013), and is an appropriate methodology to overcome the assumptions implicit to design-based approaches. DSM explores potential relationships between animal presence and environmental covariates (Friedlaender et al., 2006Herr et al., 2016) to provide potentially more ecologically reliable estimates of abundance throughout a wider area around the survey transects. Extending DSM further to incorporate presence-absence data with abundance responses to covariates has led to the development of Density Surface Hurdle Models (DSHMs; , with both DSM and DSHM commonly referred to as "model-based" approaches. Optimal line-transect surveys typically use Research Vessels fitted with double sighting platforms, following predetermined systematic transects arranged across the study region to ensure unbiased sampling coverage (Buckland et al., 2005). However, this requires dedicated vessel and personnel, a sizeable economic burden that makes such surveys costly and, as a result, they are only carried out relatively infrequently. Thus, design-based surveys of cetaceans tend not to provide information on the spatiotemporal aspects of cetacean abundance and distribution within a season. Given that most cetaceans species perform major seasonal migrations, from lower latitude breeding grounds to higher latitude feeding grounds, they are likely to be neither temporally or spatially stationary, and thus presumably will affect the predator-prey and predator-fishery interactions in both dimensions.

The International Association of Antarctic Tour Operators
(IAATO) is the organization responsible for coordinating the use of Antarctica as a tourist destination. In the past two decades, commercially driven tour operators have increased ten-fold in the Antarctic Peninsula (Bender et al., 2016), and as a consequence provides infrastructure suitable for supporting scientific research, for example as transport between monitoring stations or as opportunistic observation vessels. The use of tourist vessels as surveying platforms not only represents a fraction of the cost of research vessels, but also provides the considerable advantage of following highly repeatable transects multiple times throughout the season. These transects, however, are neither completely systematic nor completely random and lack the same systematic coverage obtained by dedicated research vessels. To address this, model-based methods Miller et al., 2013) have been developed and become the preferred tool (Williams et al., 2006), compared to traditional methods (Buckland et al., 2005) for assessing cetacean abundance from nonrandom survey data.
Developing EBM for the krill fishery requires understanding and mitigating the risks to nontarget species, both through direct and indirect competitive interactions. Given the poor understanding of the spatiotemporal dynamics in humpback whale abundance and distribution in the WAP area, we used Antarctic cruise vessels as platforms of opportunity to undertake repeated line transect surveys of humpback whales across the Bransfield Strait to develop time-evolving abundance and distribution estimates throughout the summer season. While our overarching goal was to assess the suitability of platforms of opportunity as cetacean monitoring platforms, our specific goals were to estimate intra-season variation in (1) humpback whale density and abundance, (2) distribution in relation to environmental drivers, (3) calculate temporally discrete krill consumption estimates, and (4) identify temporal overlap between key ecosystem components. Finally, we draw together our spatiotemporally resolved consumption estimates, with the assumption that whales migrate out of the area in the same way they arrived, into a conceptual model of the overlap in space and time between humpback whales, centrally foraging penguins, and the fishery.  (Table 1). The study area polygon was constructed to coincide with the main area of interest of the regional krill fishery. For analytical purposes, the five trips were grouped into three temporal windows. The first and second trips by MS Midnatsol were grouped to "C1," in order to ensure a sufficient number of sightings for analysis (Buckland et al., 2015).

| Data collection
The third MS Midnatsol trip and first MS Fram trip were grouped to "C2" due to a temporal overlap, while the final cruise of MS Fram represents "C3." On-effort periods lasted from approximately 05:00 to 22:00 local time and consisted of two or more observers, one as the dedicated recorder and the other(s) as observer(s). Two or more team members were typically stationed on effort, each dedicated as observer or recorder, but in areas of low density of animals, only one member on effort serving as both observer and recorder. During food breaks and necessary rests, effort was either maintained by one observer or ended when no observer was on duty. To the extent possible, these periods were timed to coincide with periods of low-density animal sightings to minimize the potential of missed sightings. Effort was halted in difficult sighting conditions (Beaufort Sea State >5, or visibility ≲300 m) until conditions improved. Dedicated observations were carried out from the navigation bridge of each vessel, covering either the port or starboard forward quarter (i.e., 0°-270° for port and 0°-90° for starboard, relative to the ship bearing) by scanning for cues using the naked eye, followed by binocular confirmation. Cetaceans were mostly spotted by their blows, but sometimes also through other cues such as breaching or fluking. Immediately following the sighted cue, the recorder would measure radial distance either using a distance stick (Todd et al., 2015) or binoculars (Opticron Marine PS II 7x50/C with integrated distance reticles). Sighting distance was similar for the two ships; bridge height was 13.5 m for MS Midnatsol and 13.2 m for MS Fram, yielding distances to the horizon of 8.25 and 8.17 nm for the two vessels, respectively, using an average eye height of 1.7 m.
The radial angle was measured using an angle board. Radial distance, angle (degrees), and eye height were later used to calculate the perpendicular distance of an animal from the transect line. Additional variables were recorded, such as species, group size, swim direction, and behavior. Photographs of animals were attempted if conditions allowed, but they were not prioritized. Species identifications were based on the shape of the dorsal fin, blow shape and frequency, and species-specific behavior based on relevant literature (Shirihai, 2006(Shirihai, , 2008. In the case of species uncertainty, observations were placed in general categories, such as "large baleen whale" or "like Antarctic minke."

| Data preparation and analysis
The density and abundance of humpback whales were estimated using model-based density surface hurdle modeling (DSHM; . Model-based estimates were achieved by TA B L E 1 List of the three cruises and five trips with observer effort; M1-3 represent the three trips on MS Midnatsol, and F1-2 represent the two trips onboard MS Fram; Cruise C1 (25.11.19-12.12.19), C2 (16-27.12.29) and C3 (12-18.01.20) is the identification of some combined trips for analysis; along with their respective transect length in nautical miles (nm); SA (%) is the area coverage as a percentage of the survey area; Obs(ind; n) is the number of observations and number of individuals in parenthesis; encounter rate is measured in observations per nautical mile fitting a detection function (the probability of animal detection as a function of perpendicular distance) followed by fitting generalized additive models (GAMs; Wood et al., 2016) to explore the relationship of animal presence and abundance to environmental covariates. This relationship was then used to predict animal density in the study area limited by the range of observed environmental values fitted in the GAM Kosicki, 2020;Mannocci et al., 2014). In addition, standard distance sampling estimates were also calculated, which can be found in the Appendix S1.
The detection function was fitted to the data using the combined observations from all cruises (C1, C2, and C3) in the study area for best representation of the overall detection probability. To find the best fitting detection function, multiple candidate models were tested.
Key functions include half-normal (HN), hazard rate (HR), and uniform (U; Buckland et al., 2015). Additionally, candidate models included cosine and Hermite polynomial series expansion to find the best fit.
Various covariates that may affect sightings probability were also included, such as group size, observer bias, and weather conditions. To reduce the effect of outliers when fitting the detection function, the distance data were truncated at 15% to improve model fit and remove observational outliers (Buckland et al., 2015). Model selection was based on Akaike Information Criterion (AIC; Sakamoto et al., 1986), goodness-of-fit, and Cramér-von Mises test (Darling, 1957).
Model-based results derived using DSHM protocols   covariates and validations statistics can be found in the Appendix S1.
Data were processed and handled using R Studio (RStudio Team, 2020). Detection function and design-based distance sampling abundance and density estimates were estimated using the "Distance" package (Miller et al., 2019), while model-based abundance and distribution estimates were obtained using the "dshm" package implementing the DSHM ).

| Prey consumption and temporal overlap
Daily per capita consumption estimates were derived using both traditional estimates (Kleiber, 1961 Clarke, 1980).
To estimate temporal trends, we fit our data points to a sigmoid curve representing animal arrival to the feeding ground and defined the beginning of the feeding season as when 50% of the population had arrived, based on our maximum predicted abundance. One hundred and twenty days (Lockyer et al., 1981) later, the sigmoid curve was mirrored, assuming that humpback whales leave the area at the same rate as they arrived. To calculate seasonal consumption, we used the mean daily abundance, and scaled it with daily individual consumption estimates from our traditional (Reilly et al., 2004) and contemporary (Acevedo & Urbán, 2021) sources and multiplying this to match the 120-day feeding season assumption.

| Hurdle models
The best submodels (both the binomial PA and zero-inflated Poisson AB), which included surface temperature (SST), temperature gradient (TG), and oceanographic coherent structures (FSLE) and the baseline detection function, explained 30.5 and 33.7% of the overall deviance for PA and AB, respectively (Figure 3). A complete list of all models run for the DSHM along with deviance explained and ΔAIC is shown in Appendix S1: Tables S1 and S2 for PA and AB submodels, respectively.
The presence of whales was low irrespective of sea surface temperature variations below 2°C but increased with increasing temperatures above 2°C (Figure 3). Humpback whales showed an increasing affinity to FSLE values of −0.10 to −0.25 but showed no distinct trends in values closer to zero. Finally, our data showed a nonlinear increase in the probability of presence toward higher temperature gradients. When humpbacks were present, their abundance was positively correlated to sea surface temperature, showing a rapid increase with increasing temperatures, with a peak at ~1.5°C.
The species also showed a nonlinear increasing affinity toward higher TG areas, while the response to changes in FSLE was weaker than for the other two covariates. Additional test statistics can be found in Appendix S1.

| Abundance and density estimates
Total abundance (i.e., the sum of predicted densities per grid cell) deriving from the DSHMs had a mean of 13,964 individuals, ranging from 4124 to 19,107 individuals throughout the study period (

| Consumption estimates
Mean daily krill consumption for humpback whales throughout the study period was estimated to be 12,205 kg by traditional sources (Reilly et al., 2004) and 31,601 kg by contemporary sources (Acevedo & Urbán, 2021; Table 3). Intuitively, seasonal variation of consumption is parallel to the seasonal trends in abundance mentioned above. However, overall seasonal consumption was estimated to range from 1.4 to 3.7 million tons, respective of literature.

| DISCUSS ION
We provide the first assessment of within-season spatial and temporal trends of humpback whale distribution in the WAP (Andrews- Goff et al., 2018;Weinstein & Friedlaender, 2017), and provide an update to Hedley et al. (2001)

| Spatiotemporal variation in humpback densities
Using oceanographic covariates as proxies or cues for Antarctic krill to predict distribution of baleen whales such as minke (Balaenoptera bonaerensis; Friedlaender et al., 2006, Ainley et al., 2012, fin (Balaenoptera physalus; Herr et al., 2016) and humpback whales (Friedlaender et al., 2006Herr et al., 2016;Williams et al., 2006) is common practice. Nutrient-rich Circumpolar Deep Water (Prézelin et al., 2000) heavily influences the Gerlache and Bransfield Straits (Ballerini et al., 2014), entering the study area from the west and subsequently mixing with shelf water to create Transitional Zonal Waters with Bellingshausen Sea influence (TBW; Garcıá et al., 2002). Our results indicate that humpback whales tend to be more abundant in this warmer, fresher water mass, compared to the colder, saltier water mass with Weddell Sea influence entering the study area from the northeastern part of the Bransfield Strait (Garcıá et al., 2002), indicating higher densities of krill in these waters. The presence and abundance of humpback whales declined with increasing FSLE until a lower threshold (~−0.5) suggesting that humpback whales prefer areas of higher particle retention. Such regions are known to attract predators, given that they tend to be areas with higher retention of primary productivity, such as frontal-and eddie systems, both of which are found in the Gerlache and Bransfield Straits (Anadón & Estrada, 2002;Sangrà et al., 2011). Combined, our study underlines the importance of western TBW influence and submesoscale hydrographic structures.
The hurdle models predicted an overall low density early in the season, with hotspots mostly centered on the northern Gerlache Strait, followed by high aggregations of humpback whales in the northern Gerlache Strait and southern-and central Bransfield Strait during the "mid-season," ensued by a dispersion of animals throughout the study area towards the end of January. Our findings emphasize the importance of the northern Gerlache and western Bransfield Straits as a hotspot for foraging humpback whales through the first half of the summer feeding season. These results align well with regional Antarctic krill spawning events taking place in the Gerlache Strait in late December (Huntley & Brinton, 1991), and humpback whales taking advantage of the local spawning hotspots and nursery areas for immature krill (Cleary et al., 2016;Perry et al., 2019), exhibiting size-dependent predation, preferring smaller size ranges of krill (≤34 mm; Friedlaender et al., 2008Friedlaender et al., , 2009Santora et al., 2010).
Future work should examine whether these results are consistent over time.
Inherently, cruise vessels will not have the same spatial coverage as dedicated research vessels. Model-based approaches such as two-part DSHMs work to overcome this by using habitat models to predict distribution throughout the area of interest (Gowan & Ortega-Ortiz, 2014;Scott-Hayward et al., 2015;Waggitt et al., 2020;Warwick-Evans et al., 2021). This methodology is very beneficial for spatial management analyses in remote areas using limited datasets but does include more uncertainty compared to more traditional methods (Buckland et al., 2005). Comparing our transects ( Figure 1), C1 has the least coverage, mostly sailing in the western parts of the study area, while the C2 and C3 transects also cover northern and eastern areas. Thus, C1 predictions may have a higher level of uncertainty, specifically in eastern areas. However, these predictions were made within the confines of observed environmental values . Further, the insight gained in the spatiotemporal aspect of the ecosystem agrees well with previous studies using similar analyses (Friedlaender et al., 2021;Herr et al., 2016;Williams et al., 2006), and will be useful in spatial management and conservation efforts.

| Abundance and recovery
The Southern Ocean populations of humpback whales are increasing significantly (Pallin et al., 2018;Tulloch et al., 2018), some of which have recovered to near pre-exploitation levels (Zerbini et al., 2019).
Comparable abundance estimates in the area stem from Hedley et al. (2001), concluding with an abundance of approximately 7000 individuals, which is ca. 12,000 lower than the highest abundance in this study (19,107 individuals). This corresponds to an annual population growth of 5.1%, which agrees with the estimates of between 3.1% and 11.8% found elsewhere in the literature (Stevick et al., 2003;Zerbini et al., 2010).
Marine ecosystems are highly dynamic, and this region is no exception. The WAP ecosystem displays significant seasonal dynamics in primary production, with a consequent high variability in the abundance and distribution of upper trophic predators such as migratory baleen whales. Some Antarctic feeding grounds are used by several subpopulations of humpback whales, such as the northern Antarctic Peninsula (Albertson et al., 2018;Dalla Rosa et al., 2008;Robbins et al., 2011), and these areas may be experiencing compounding predatory pressure on krill with a collective subpopulation growth rather than one singular. As a consequence, abundance estimates and population growth will be higher in the feeding grounds compared to the breeding grounds. This has implications for management; when managing a fishery in the feeding ground, data from breeding grounds may give an insight as to the status of upper trophic competition but will underestimate the abundance and growth rate of the population in the management area.

Relative abundance (N)
Traditional ( The dynamic nature of these animals underlines the need for both fine-and coarse-scale data collection on krill and its predators at regular intervals in future in order to appropriately determine the degree and significance of predatory pressure. Our approach, using cruise vessels as observation platform, supports earlier work by Williams et al. (2006) in identifying a cost-effective mechanism to monitor marine mammals and seabirds.

| Potential interactions between humpback whales and fisheries
With the contraction of the Antarctic krill fishery around the Antarctic Peninsula and Scotia Sea over recent decades, there is an increasing concern over interactions with and impact of the fishery and other upper trophic predators such as penguins and seals. In the context of direct competition, both spatial and temporal overlap between the components need to be in place (Hinke et al., 2017).
Comparing our results spatially with average fishing activity, the past decade (Appendix S1: Figure S4), fishing activity and humpback whales overlap along the frontal systems of the central Bransfield Strait and the northeast Gerlache Strait.

| Consumption estimates and temporal overlap
Traditional estimates used to quantify the energetic needs of animals have been based on the Kleiber Equation (Kleiber, 1961), addressing basal metabolic rate. However, recent studies have incorporated detailed physical and biological parameters by high-resolution using high-resolution tagging of humpback whales and in situ hydroacoustic prey-field density estimates to more completely assess the caloric requirements of cetaceans (Acevedo & Urbán, 2021;Smith et al., 2015). Additionally, Acevedo and Urbán (2021) Figure 5a). However, the latter half of their feeding season, presumably lasting from February to April, has yet to be properly quantified at a larger scale. The 120-day feeding season assumption has been useful but may be outdated.
Humpback whale telemetry studies have indicated that movement patterns and departure dates in the WAP may be dictated by sea-ice cover, with a mean departure date of May 12th and some individuals remaining in the northern WAP until July (Weinstein & Friedlaender, 2017 (Figure 5b).
This species, among others, has flourished in the absence of large baleen whales (Fraser et al., 1992), but have in recent years been experiencing significant declines in population sizes throughout their range (Dunn et al., 2019;Naveen et al., 2012). While the direct cause of this decline is currently unknown, some literature attributes this deterioration to a combination of climate change, pollution, fishing activity, and habitat loss (Krüger et al., 2020;Trathan et al., 2015). Surely some of these factors may be important, but by excluding important predatory interactions such as predation by humpback whales, there is little hope to unraveling the casual mechanisms behind the declining performance in penguins in the WAP area. Additionally, the Antarctic is at the forefront of climate change effects (Meredith & King, 2005), but the ability to correctly measure the biotic effects of climate change is severely limited without a complete understanding of the ecosystem dynamics. This begs the question of how the ecosystem will respond to an increasing abundance of mass krill consumers, such as the recovering cetacean populations.

| CON CLUS ION
Our study demonstrates the utility of platforms of opportunity such as tourist vessels as a viable means to collect data in support of EBM in the WAP area. Platforms of opportunity may potentially be the only logistically and financially feasible approach to gather crucial fine-scale data regarding cetacean abundance, at least until satellite-derived abundance estimates (Bamford et al., 2020) are reliable enough to be considered suitable from a management perspective. We have provided some insight of the spatial and temporal distribution of humpback whales in the WAP during the first half of their feeding season and found that the northern Gerlache Strait is an important feeding area for humpback whales, with the majority of animals arriving in late December. Our results suggest limited spatiotemporal overlap between the humpback whales and the Antarctic Krill fishery in the area. However, if the assumed seasonal distribution, in which whales occur in high numbers between mid-December to late-April, is correct, then there is potential for interaction between fisheries and humpback whales in April. Future efforts should therefore focus on subsequent months and with regular intervals, to fully quantify the extent of the humpback whale feeding season, as well as other regions with potential for overlap with fishery activity. Our study supports the contention that humpback whale populations are recovering from historical harvesting, that their feeding season directly coincides with critical reproductive dates of chinstrap penguins and that cetaceans are a component of the ecosystem that requires greater consideration in the context of fisheries management in the Antarctic Peninsula.

ACK N OWLED G M ENTS
We thank Hurtigruten Expeditions and their Chief Scientist Dr.
Verena Meraldi for providing field support and logistics. We thank the Captains Ole Johan Andreassen, Bjørn Ivar Pedersen, and Martin Iversen, along with the wonderful crew and expedition teams, for their hospitality on board MS Fram and MS Midnatsol.

CO N FLI C T O F I NTE R E S T
The authors have no conflict of interest to declare.

DATA AVA I L A B I L I T Y S TAT E M E N T
All observational data from the two survey vessels can be found at https://doi.org/10.5061/dryad.rfj6q 57c9.