# Modelling the influence of environmental factors on the acoustic presence of blue whale populations in the southern Indian Ocean

Mathilde Michel, Baptiste Alglave, Maxime Olmos, Maëlle Torterotot, Auriane Virgili, Salomé Martin-Marin, Jean-Yves Royer, Flore Samaran

<figure><img src="image_1.png" /></figure>

## To cite this version:

Mathilde Michel, Baptiste Alglave, Maxime Olmos, Maëlle Torterotot, Auriane Virgili, et al. Modelling the influence of environmental factors on the acoustic presence of blue whale populations in the southern Indian Ocean. Scientific Reports, 2025, 15 (1), 23218 (16 p.). 10.1038/s41598-025-02941-9. hal-05206173

HAL Id: hal-05206173

https://hal.science/hal-05206173v1

Submitted on 13 Aug 2025

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

Check for updates

OPEN

# Modelling the influence of environmental factors on the acoustic presence of blue whale populations in the southern Indian Ocean

Mathilde Michel<fnref n="1" />✉, Baptiste Alglave², Maxime Olmos³, Maëlle Torterotot<fnref n="1" />, Auriane Virgili⁴, Salomé Martin-Marin<fnref n="1" />, Jean-Yves Royer⁵ & Flore Samaran<fnref n="1" />

Blue whales in the Indian Ocean have been severely depleted by previous extensive commercial whaling. A good understanding of their spatio-temporal distribution is crucial for conservation. The songs of three blue whale acoustic populations - Antarctic blue whales (*Balaenoptera musculus* intermedia, ANT BW) and pygmy blue whales (*B. musculus brevicauda*) from the Southeast (SEIO PBW) and Southwest Indian Ocean (SWIO PBW) - were analyzed using 13 years of passive acoustic recordings from 10 sites in the southwest Indian Ocean. Although blue whale vocalizations comprise both songs and non-song calls (e.g., D-calls), the present study concentrates on the examination of songs. Generalized additive models (GAMs) were used to relate acoustic presence, measured by the number of positive minutes per day (averaged weekly), to environmental drivers such as sea surface temperature (SST), chlorophyll-a concentrations, and sea ice extent. These models allowed predictions of blue whale acoustic presence across the region. Empirical orthogonal functions (EOFs) were applied for dimensionality reduction to identify key habitats, including the Kerguelen Plateau and Madagascar Basin, which may serve as important feeding and resting zones based on acoustic presence and environmental data. Antarctic blue whales were predominantly detected in austral winter and spring, associated with lower SST and higher chlorophyll-a. In contrast, SEIO and SWIO pygmy blue whales were more frequent in summer and autumn, with some overlap suggesting ecological interactions. These findings lay the groundwork for targeted conservation efforts to protect critical blue whale habitats in a rapidly changing ocean.

**Keywords** Spatio-temporal modelisation, Blue whale, Acoustic recording

Blue whale (*Balaenoptera musculus*) populations in the southern hemisphere have been severely depleted by previous extensive commercial whaling. As a consequence, the International Whaling Commission (IWC) established the Indian Ocean Whale Sanctuary (IOWS) in 1979, prohibiting commercial whaling in the region<fnref n="1" />. The rare nature of blue whales, combined with the vast and remote expanses of their habitat, complicates efforts to observe them directly. Given these challenges, passive acoustic monitoring has emerged as the most effective method to study these cetaceans, allowing researchers to track their presence and behavior through sounds they produced²⁻⁷. This approach is particularly appropriate for blue whales, which are highly vocal animals characterized by long low-frequency sounds⁸.

Blue whale sounds are broadly categorized into two types: songs and non-song calls. Songs are complex long-duration calls and are made up of repeated units organized in series, separated by regular intervals (intercall interval, ICI)⁵. They are recorded in all seasons in both feeding and breeding areas and are characterized by a stereotype time-frequency structure, allowing differentiation between populations³,⁴. In contrast, non-song calls, such as D-calls, are shorter, frequency-modulated vocalizations typically associated with foraging

<footnote n="1">¹ENSTA, IPP, CNRS Lab-STICC, 29200 Brest, France. ²Université Bretagne Sud, Lab-STICC, 56000 Vannes, France. ³DECOD (Ecosystem Dynamics and Sustainability), L'Institut Agro, Ifremer, INRAE, 29280 Plouzané, France. ⁴Share the Ocean, 56870 Larmor-Baden, France. ⁵Geo-Ocean, CNRS, Université de Brest, Ifremer, 29280 Plouzané, France. ✉email: mathilde.michel@ensta.fr</footnote>
behavior⁹. D-calls exhibit highly variable frequencies and are produced by both males and females, without a clear geographic variation in their time-frequency shape¹⁰⁻¹³.

Temporally stable differences in songs provide characteristics that can be used to distinguish the two subspecies and to define several acoustic populations of pygmy blue whales in the Indian Ocean¹⁴. Although specific to each subspecies and acoustic population, blue whale songs share similar acoustic properties: the duration of song units exceeds 15 seconds, their amplitude is greater than 180 dB re 1 µPa at 1 m, and their pitch ranges from 20 to 100 Hz. Songs can be composed of multiple units, depending on the subspecies or acoustic population, and intercall intervals are periodic and do not exceed a few minutes. The year-round recording of blue whale vocal activity in various locations suggests that individuals produce acoustic emissions continuously throughout the year at the population level³,⁴,⁷,¹⁵. Furthermore, long-term analyzes have revealed a gradual decline in the frequency of blue whale calls in recent decades, although the underlying causes remain uncertain¹⁶. The differences in vocal repertoire between the blue whale subspecies and the acoustic populations of pygmy blue whales provide a unique means of tracking their migratory movements throughout an ocean basin⁷,¹⁵,¹⁷,¹⁸.

The southern region of the Indian Ocean hosts seasonally several acoustic populations of blue whales, including different subspecies and acoustic populations. Two subspecies are regularly present in this area: the Antarctic blue whale (B. m. intermedia) and the pygmy blue whale (B. m. brevicauda)¹⁹. There are notable differences in the way populations use their habitats, each of which exhibits distinct migration patterns. Antarctic blue whales (ANT BW) follow a seasonal migration from subantarctic regions in austral autumn to subtropical areas in austral winter, before returning to subantarctic waters in austral spring (September to November)²,³,⁵,²⁰.

The acoustic populations of pygmy blue whales in the Southeast Indian Ocean (SEIO PBW) and the Southwest Indian Ocean (SWIO PBW) show distinct distribution patterns within the Indian Ocean and are predominantly observed during the austral autumn and winter seasons³,⁴,⁷,¹⁵. The SEIO PBW acoustic population is recorded mainly along the southwest coast of Australia, predominantly between March and May²¹,²², while the SWIO PBW acoustic population is found in regions further west, within the southwestern Indian Ocean²³⁻²⁵. However, recent findings reveal that these acoustic populations have been identified beyond their previously recognized distributional limits. SEIO PBW have been observed farther west than expected, in proximity to Kerguelen Island³,⁵,⁷, while SWIO PBW songs have been detected further south and west, near the Crozet Archipelago and southwest of Amsterdam Island³,⁷.

Oceanographic parameters, including sea surface temperature, chlorophyll-a concentration, and sea ice extent, are critical in determining prey availability and subsequently influence whale migratory movements and foraging behavior. Previous research has concentrated primarily on modelling the associations between these environmental variables and the distribution of cetaceans on extensive spatial scales, predominantly using visual or aerial observation data²⁶⁻²⁸. Other studies have used passive acoustic data, but have been limited to smaller spatial and temporal scales, often focusing only on correlations between environmental variables and acoustic presence²⁹,³⁰. Building on this foundation, this study seeks to extend these analyses to a larger geographic scale and a longer time frame. By integrating passive acoustic monitoring with oceanographic data, we provide comprehensive predictions for the entire study region, thus enhancing our understanding of key blue whale habitats and supporting conservation efforts.

Specifically, this study aims to improve the characterization of blue whale habitat based on 13 years of acoustic recordings collected in the southern Indian Ocean. Using a generalized additive model (GAM) framework, we evaluated the influence of oceanographic covariates on the acoustic presence, defined as an acoustic intensity metric based on the number of positive detection minutes per day, averaged weekly. This modelling approach allows us to (i) identify key environmental drivers of whale distribution and assess their non-linear effects, and (ii) predict spatio-temporal variations in blue whale acoustic presence, including inter-annual and seasonal patterns. In a subsequent phase, we will delineate functional areas using Empirical Orthogonal Functions (EOFs), a spatio-temporal dimension reduction technique. The identified spatial and temporal patterns are then compared with the existing literature to ensure alignment with established knowledge. These results provide valuable information for the conservation of this iconic species.

# Methods

## Area Description

The study area in the southern Indian Ocean extends from 20°S to 60°S latitude and from 45°E to 90°E longitude. In this region, the SST varies significantly, ranging from about 28 °C at lower latitudes to nearly 0 °C near the South Pole. This gradient is influenced by the separation of the water masses at the subantarctic and polar fronts. The subantarctic front, at about 45°S, separates warm mid-latitude waters from cooler Antarctic waters, improving biological productivity³¹. Further south, between 50°S and 60°S, the polar front separates the subantarctic from the polar waters, promoting the formation of sea ice³². Antarctic sea ice expands seasonally, peaking in the austral winter and covering much of the southern Indian Ocean. These variations correspond to the natural cycles of sea ice formation and melting dictated by the regional climate. The chlorophyll-a distribution shows spatial variability, with high concentrations in nutrient-rich colder waters. On the Kerguelen Plateau, concentrations peak at 3.96 mg/m³ due to nutrient accumulation and upwelling³³. In contrast, lower latitudes have lower concentrations, around 0.02 mg/m³, due to subsurface features that affect currents and nutrient transport³⁴. Phytoplankton blooms on the polar front are highly variable in space and time, especially in December, and are associated with increased solar radiation and an increase in SST in the southern ocean³⁵.

## Acoustic data collection

The acoustic data used in this study were collected by a wide network of autonomous hydrophones called the Hydro-Acoustic Observatory of Seismicity and Biodiversity (OHA-SIS-BIO)³⁶. This network was deployed
in the southern Indian Ocean from December 2009 to February 2023 (see Fig. 1). Initially and until 2016, it included five sites: south of la Réunion Island (MAD), north of the Crozet Archipelago (NCRO), west of Kerguelen Island (WKER), southwest (SWAMS) and northeast (NEAMS) of the Saint-Paul and Amsterdam Islands. An additional site south of the Southeast Indian Ridge (SSEIR) was installed in 2014. In 2017, the geometry of the network was slightly modified to improve coverage in the northern regions and better pinpoint the location of southern blue whale feeding grounds. Five new sites were added: near the junction of the three Indian Ocean seafloor spreading ridges (RTJ), west (MAD-W) and east (MAD-E) of the Madagascar Basin, south of the southwest Indian ridge (SSWIR) and south of the Kerguelen Plateau (ELAN), while the MAD and NCRO sites were no longer deployed.

Each mooring consisted of an anchor, an acoustic release, an adjustable line, and a submerged buoy that houses the recording system. The instruments were deployed at depths ranging from 1000 to 1300 m on the axis of the Sound Fixing and Ranging (SOFAR) channel, which facilitates the transmission of low-frequency sounds over long distances³⁷. The recording was continuous at 240 Hz sampling with a 24-bit analog-to-digital conversion. The hydrophones used had a sensitivity of -163.5 dB re 1V/µPa, except for the moorings at SWAMS in 2015 and 2016, where the sensitivity was -168.6 dB re 1V/µPa. The data collected were stored locally on each device³⁸. The moorings were maintained during the annual scientific cruise of the R/V Marion Dufresne in the French Southern and Antarctic Territories in January/February each year. The database used in this study contains almost continuous acoustic recordings over 13 years, from 2010 to 2022.

## Acoustic data analysis

In this study, blue whale song detections from 2010 to 2018 were taken from Torterotot et al.⁷, and from 2019 to 2022 were processed using the same detection method, as follows.

A detection algorithm using dictionary learning and sparse representation was applied to automatically detect blue whale songs³⁹. This technique uses large dictionaries to capture song variability and uses linear combinations of dictionary entries to represent song complexity. Without relying on a fixed template, the algorithm builds dictionaries directly from temporal call signals extracted from the data and reconstructs

<figure><img src="image_4.png" /><figcaption>**Fig. 1.** Locations of the mooring sites of the OHA-SIS-BIO hydrophone network. The blue square represents the Kerguelen Plateau region. The map was created using R (version 4.1.2; https://www.r-project.org/)⁵¹.</figcaption></figure>
observed signals with sparse combinations of waveform entries. Higher-quality reconstructions correspond to higher similarity scores³⁹,⁴⁰. The signal used for recognition was the entire call for the ANT BW and a harmonic of a single unit of the call for both pygmy blue whales (see the white frames in Fig. 2). These harmonics were chosen for their high amplitude and frequency range that did not overlap with other blue whale songs.

The performance of this detection algorithm was evaluated on a substantial subset of manually annotated data spanning different seasons and locations within the OHA-SIS-BIO array. By setting the detection threshold to allow for one false positive per hour, we achieved a 90% recall for all calls with an SNR greater than 0 dB. A detailed review of all detections at the WKER location throughout the 2015 recordings confirmed that the number of false positives was very low, most of them caused by calls masked within the chorus⁴¹. The year-to-year frequency decrease in blue whale calls is very small (from -0.12 to .0.32 Hz/yr depending on the species)¹⁵. As a result, the dictionary from one year can still be applied effectively to the adjacent year. Although some calls may be missed due to a slight drop in frequency, most will still be detected successfully. A new dictionary was then built for each year of data, based on the calls detected by the adjacent year's detector. This approach effectively accounts for the frequency shift in the calls. A detailed evaluation of the algorithm performance in the dataset used is provided by Torterotot et al. (2019)⁴¹.

## Environmental variables

In this study, three dynamic environmental variables were analyzed to explain the acoustic presence of blue whales derived from acoustic detections. These variables were chosen because they effectively characterize the environment in the study area, as discussed in the “Area description” section.

* Sea Surface Temperature (SST) Changes in sea surface temperature can signal shifts in regional biological and chemical processes that affect productivity. The SST data in this study come from the Copernicus Marine global ocean reanalysis (DOI: 10.48670/moi-00021), which covers the period from 1992 to 2024, with a high spatial resolution of 1/12° (approximately 8 km). The reanalysis includes observations from satellites such as TOPEX POSEIDON and ERS-1 and models of 3D thermodynamic and dynamic variables. SST is measured in degrees Celsius (°C) with daily temporal resolution. Daily data were extracted and averaged over a 50 km area around each hydrophone, assuming the songs originated in that region. In addition, these averages were calculated weekly to align with the response variable under study. Figure 3 summarizes the SST data for the study area for each season.

* Chlorophyll-a concentration Chlorophyll-a concentration, a common indicator of primary productivity, is considered as a distal rather than proximal variable for prey such as zooplankton⁴². Although proximal variables typically better explain species distributions, the environmental variables chosen reflect prey distribution patterns⁴². Chlorophyll-a data were obtained from Mercator-Ocean, from January 1, 1993, to February 29, 2024, with a resolution of 1/4° and concentrations reported in mg/m³ (DOI: 10.48670/moi-00019).

<figure><img src="image_5.png" /><figcaption>**Fig. 2.** Spectrograms of (a) SEIO PBW songs (b) SWIO PBW songs, and (c) ANT BW songs from the southern Indian Ocean. The units selected for the detector are framed in white.</figcaption></figure>
were collected daily and averaged over a 50 km area per site to account for spatial variability. An overview of these data for the study area is presented in Fig. 3, which summarizes the chlorophyll-a data for each season.

* Sea Ice Concentration Sea ice significantly affects marine habitats, the availability of whale prey, and the migration of whales<sup>28,43</sup>. The sea ice concentration data come from the NOAA/NSIDC Climate Data Record, Version 4<sup>44</sup>, based on passive microwave measurements following NOAA criteria. It combines estimates from NASA Team and Bootstrap algorithms, using the highest precision<sup>45</sup>. Data are available at 25 km x 25 km resolution, daily and monthly. We analyze sea ice in a specific area (50°E to 85°E, -70°S to -10°S) by calculating pixels that exceed a concentration of 15% ice, a standard threshold to distinguish sea ice from open water<sup>32</sup>. An overview of these data for the study area is presented in Fig. 3, which summarizes the sea ice concentration data for each season.

* Geographic coordinates Geographic coordinates of the study sites were included in the models to capture spatial variability not captured by the other covariates and to account for large-scale distributions resulting from the geographical extent of the subspecies<sup>5</sup>.

* Temporal Parameters: Months and years Months and years were used as temporal predictors in the models to account for seasonal and annual variations in the acoustic presence of an acoustic population. These variations reflect the migratory patterns of the species and the inter-annual variability<sup>7</sup>.

In this study, the austral seasons are defined as follows: summer includes December, January and February; autumn includes March, April, and May; winter includes June, July, and August; and spring includes September, October, and November. To simplify the visual representation of the seasonal patterns in Figs. 3 and 5, a single representative month was selected for each season. Specifically, January was designated for summer, April for autumn, July for winter, and November for spring.

## Habitat modelling

The environmental variables (described above) may not have a linear relationship with the spatio-temporal distribution of the whales. Consequently, we modeled the acoustic presence of the whales with a spatio-temporal GAM<sup>46</sup>. This approach accounts for non-linear relationships with covariates, can deal with zero inflation through non-standard observation models (e.g. tweedie distribution), and models spatio-temporal correlations through spatial and temporal splines<sup>47,48</sup>.

To fit the model, we used observations from the average number of minutes with positive detections per day (DPM), averaged weekly. Observations were zero-inflated positive with heavy tails. Consequently, we tested two different likelihood distributions that can be considered for count data that exhibit overdispersion and excess zeros (Tweedie or Negative Binomial, Table 1) to model the observations. The averaged number of minutes with positive DPM $\mu_i$ is modeled through the equation:

$$g(\mu_i) = \mu_0 + \beta_{\text{Month}} \cdot \text{Month}_i + \beta_{\text{Year}} \cdot \text{Year}_i + f_1(\text{Chla}_i) + f_2(\text{SST}_i) + f_3(\text{MSI}_i) + f_4(\text{lon}_i, \text{lat}_i)$$

where g is a link function (here log), $\mu_0$ is the intercept, the linear terms related to months and years are $\beta_{\text{Month}}$ and $\beta_{\text{Year}}$, $f_j, j \in \{1, \dots, 4\}$ are the smoothed functions related to the covariates chlorophyll-a (Chla), sea surface temperature (SST), mean sea ice (MSI), and the spatial effect related to longitude (lon) and latitude (lat). Collinearity among environmental covariates (Chla, SST and MSI) was assessed using pairwise correlation coefficients. A strong negative correlation was found between Chla and SST (r = -0.81), but both variables were retained in the model due to their ecological relevance and because AIC-based model selection supported their inclusion.

From this model, we provide spatio-temporal maps of the acoustic presence of the whales by predicting $\hat{\mu}$ over a 50 km² spatial grid that covers the spatial domain at each time step. Different models were tested to investigate the spatio-temporal variations in the acoustic data (Table 1). In the M2 reference model, we modeled acoustic data as a function of months, years, environmental covariates, and spatial effects related to longitude and latitude. Two models of lower complexity were also considered (M0, M1, Table 1). We used Akaike's information criterion (AIC) for model selection as a measure of the goodness of fit penalized by the number of parameters to account for model parsimony<sup>49</sup>. To avoid overfitting, we split the data into 80% training and 20% validation sets, and cross-validation was performed using Spearman's rank correlation ($R_{\text{Spearman}}$) to assess model performance. The results showed improved performance for more complex models that incorporated temporal and spatial effects. Furthermore, the contribution of each covariate to the explained variance in the GAM was assessed by comparing the deviation of the full model with reduced models, where one covariate was removed at a time. This approach enabled calculation of the proportion of variance explained by each covariate as the relative change in deviation, normalized by the null deviation. The asterisks (*) indicate significant contributions (p < 0.001). These significance values are provided to indicate the importance of each covariate in explaining the variance in the model (Fig. Supplementary 4). Finally, the coefficients of variation of the model predictions were calculated as the standard deviation divided by the mean of the predictions (Fig. Supplementary 5). The model was implemented using the 'mgcv' package<sup>50</sup> in R<sup>51</sup>.

## EOFs to identify functional areas

To identify functional areas based on GAM spatio-temporal predictions, we relied on standard spatio-temporal dimension reduction methods known as Empirical Orthogonal Functions (EOFs)<sup>52,53</sup>.

EOFs are derived from climate science to extract the main spatial and temporal patterns in spatio-temporal data. They are principal component analysis (PCA) in a matrix where locations are individuals and time steps are variables. They have been widely used in marine ecology to analyze ecosystem variability, species distribution, and phenology<sup>53,54</sup>. Here, the spatio-temporal data are the model predictions $\mu$ at each location of the spatial
<figure><img src="image_6.png" /><figcaption>b)</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>c)</figcaption></figure>

<figure><img src="image_8.png" /></figure>

**Fig. 3.** Monthly averages of covariates over the entire study area in the Indian Ocean for the four seasons, represented by January (summer), April (autumn), July (winter), and November (spring), for the years 2014 and 2022. The variables represented are: (a) SST (°C), (b) log (Chl-a) (mg/m³), and (c) MSI (%). The red dots indicate the locations of the hydrophones used in this study.

grid for each time step (month of each year). They are stacked in a matrix $n \times m$ where n is the number of spatial locations (here X) and m is the number of time steps (here Y).

EOFs summarize information from a sequence of maps into a smaller set of *r* maps, denoted $p_j = (p_j(x_1), \dots, p_j(x_n))$ where *n* is the number of spatial locations. These terms (called spatial factors) best represent the spatio-temporal variability patterns of the spatio-temporal variable. The first spatial factor of the
<table><caption>**Table 1.** Summary of the hypotheses tested, the associated model configurations, and AIC values attributed to each model for each acoustic population: ANT BW, SEIO PBW, and SWIO PBW. Explained deviance (%) quantifies the proportion of variation explained by the model. ΔAIC is the difference in the AIC score between the best model (*) and the model being compared.</caption><thead><tr><th>Models</th><th>Spatio-temporal variations in DPM are explained by</th><th>Ecological hypothesis</th><th>Acoustic population</th><th>Distribution</th><th>Deviance explained (%)</th><th>AIC</th><th>ΔAIC</th><th>Cross Validation (RSpearman)</th></tr></thead><tbody><tr><td rowspan="3">M0</td><td rowspan="3">Chlorophyll-a + sea surface temperature + mean sea ice</td><td rowspan="3">Only depends on environmental variables, does not take into account monthly and annual variations in acoustic presence</td><td>ANT BW</td><td>Tweedie</td><td>11.6</td><td>92019</td><td>3234</td><td>0.38</td></tr><tr><td>SEIO PBW</td><td>Tweedie</td><td>43.5</td><td>34165</td><td>3748</td><td>0.49</td></tr><tr><td>SWIO PBW</td><td>Tweedie</td><td>20.9</td><td>76628</td><td>6755</td><td>0.47</td></tr><tr><td rowspan="3">M1</td><td rowspan="3">M0 + Month and Year effect</td><td rowspan="3">Account for inter-annual and monthly changes in acoustic whale distribution</td><td>ANT BW</td><td>Tweedie</td><td>16.4</td><td>90471</td><td>1686</td><td>0.45</td></tr><tr><td>SEIO PBW</td><td>Tweedie</td><td>52</td><td>31952</td><td>1535</td><td>0.52</td></tr><tr><td>SWIO PBW</td><td>Tweedie</td><td>32.2</td><td>72155</td><td>2282</td><td>0.58</td></tr><tr><td rowspan="7">M2</td><td rowspan="7">M1 + longitude and latitude effects</td><td rowspan="7">Take into account the areas of distribution of species and migration patterns</td><td>ANT BW</td><td>Tweedie</td><td>21.2</td><td>88785</td><td>0*</td><td>0.51</td></tr><tr><td>ANT BW</td><td>Negative binomial</td><td>20.9</td><td>88809</td><td>24</td><td>0.51</td></tr><tr><td>SEIO PBW</td><td>Tweedie</td><td>57</td><td>30417</td><td>0*</td><td>0.54</td></tr><tr><td>SEIO PBW</td><td>Negative binomial</td><td>56.9</td><td>31569</td><td>1152</td><td>0.54</td></tr><tr><td>SWIO PBW</td><td>Tweedie</td><td>37.6</td><td>69873</td><td>0*</td><td>0.63</td></tr><tr><td>SWIO PBW</td><td>Negative binomial</td><td>37.1</td><td>70617</td><td>744</td><td>0.63</td></tr></tbody></table>

EOF $p_1$ captures the highest amount of variability in the spatio-temporal data; the second $p_2$ is orthogonal to the first spatial factor $p_1$ and captures the second highest amount of variability. The same principle applies to the next dimensions.

For each dimension $j$, the spatial factor is related to a set of temporal indices called (temporal) loadings defined as $\alpha_j = (\alpha_j(t_1), \dots, \alpha_j(t_m))$, where $m$ is the number of time steps. If the loading of a time step is positive ($\alpha_j(t) > 0$), then the spatio-temporal variable is distributed according to the related spatial factor $p_j$. If the loading of a time step is negative ($\alpha_j(t) < 0$), the spatio-temporal variable is distributed according to the opposite of the related spatial factor $-p_j$. More details on the interpretation of EOFs can be found in Alglave et al. (2024)$^{55}$.

As with standard PCA, EOFs can be combined with clustering analysis (here Hierarchical Clustering on Principal Components) to regroup locations that have similar temporal trends (clustering on spatial factors) or time steps that have similar spatial patterns (clustering on temporal loadings)$^{53}$. Clusters of locations can be interpreted as distinct functional habitats, while clusters of time steps can be interpreted as ecological seasons. Both can be plotted on a single graph, and locations can be related to seasons if they are in the same direction on the graph (similar to the case in PCA).

Note that both EOF dimensions and clustering can be strongly affected by extreme values. Given the high predictions surrounding Kerguelen Island, which distort the EOF representation of broad-scale spatial patterns, we exclusively retained predictions outside the Kerguelen area to perform the EOFs and clustering for the SEIO PBW and SWIO PBW acoustic population. EOF analysis and clustering were carried out using the 'FactoMineR' library$^{56}$ in R$^{51}$.

# Results

## Acoustic effort

Acoustic effort and detection data for each site and acoustic population are summarized in Table Supplementary 1. The total recording effort exhibited variation between sites, ranging from 1,630 days (NCRO) to 4,750 days (WKER), with a range of recording durations between sites, totaling approximately 2.3 million minutes (NCRO) to 6.8 million minutes (WKER). The proportion of days with acoustic presence exhibited variation across sites and acoustic populations, ranging from 0.0% at MADW for the SEIO PBW acoustic population up to 95.6% at NEAMS for the ANT BW acoustic population. Acoustic effort data reflect continuous and intermittent recording throughout the study period, with the frequency of detections varying over time. Although some sites show nearly constant acoustic presence (e.g., MAD, NEAMS), others exhibit more intermittent detections (e.g., RTJ, ELAN). The acoustic presence of each acoustic population and the recording periods are detailed in the supplementary materials (Fig. Supplementary 1, Supplementary 2, Supplementary 3).

## Spatio-temporal variations in environmental covariates

The thermal gradient in the study area (Fig. 3) ranges from -4°C to more than 25°C, exhibiting a discernible latitudinal pattern, particularly evident during austral summer and autumn, with a SST peak in the northwest. These temperatures remain consistent from year to year and demonstrate minimal variability (Fig. 3a). Conversely, chlorophyll-a concentration exhibits an inverse latitudinal gradient with concentrations ranging from 7.38 mg/m$^3$ in nutrient-rich zones to 0.04 mg/m$^3$ in nutrient-poor zones, reaching its peak near the Kerguelen Plateau in austral summer. In 2014, inter-annual variability was observed, with a notable peak in primary productivity in austral spring, compared to the other years (Fig. 3b). The sea ice cover (Fig. 3c) shows stable seasonal patterns, with slight variability, peaks in austral winter, with ice cover ranging from 0% to almost 50% of the region.
## Model selection

Table 1 delineates the hypotheses, model configurations, and AIC values pertinent to each acoustic population: ANT BW, SEIO PBW and SWIO PBW. For all species examined, model M2, characterized by its incorporation of spatial, monthly and annual effects in conjunction with all environmental covariates such as chlorophyll-a concentration, sea surface temperature, and sea ice cover, exhibited the lowest AIC and the highest explanatory power regarding deviation, although anticipated. Among the two tested distributions, the Tweedie distribution demonstrated superior performance in the context of model M2. In addition to AIC-based model selection, cross-validation using Spearman's rank correlation ($R_{Spearman}$) further supported the superiority of model M2, revealing improved predictive performance for more complex models incorporating temporal and spatial effects. Ultimately, the chosen model represented 21.2% of the deviation for ANT BW, 57% for SEIO PBW, and 37.6% for SWIO PBW, and was validated through both AIC and cross-validation metrics.

## Monthly habitat models and predicted distributions

### ANT BW

ANT BW songs were detected at all recording stations, demonstrating a large representation within the study area (Fig. Supplementary 1). Detection rates are significantly higher at the NEAMS site and comparatively lower at the SSEIR and RTJ sites (Fig. Supplementary 1). The covariates that contribute the most significantly to the variance in our model include sea surface temperature, chlorophyll-a concentration, sea ice extent, GPS coordinates, temporal year parameters, and month. These covariates were identified as statistically significant ($p < 0.001$)) (Fig. Supplementary 4). There is a stable positive correlation between chlorophyll-a concentration and acoustic presence, except at very high chlorophyll-a concentrations, where further increases in acoustic presence are observed (Fig. 4a). In particular, this correlation starts with a logarithmic value of $log(mean Chla) = -1.8$ mg/m³ (Fig. 4a). Seasonal peaks in acoustic events were notably discernible, particularly during austral winter and spring (Fig. 4a). In contrast, there is a negative correlation between temperature and acoustic presence, markedly above 20°C, with a plateau and even a second peak in the range of 10 to 20°C (Fig. 4a). The relationship between sea ice extent and acoustic presence reveals variability at lower percentages of sea ice coverage, with an increase in daily positive minutes (DPM) concurrent with initial ice formation, followed by a secondary positive peak at higher percentages of ice coverage (Fig. 4a). The annual variability in the acoustic presence exhibited significant peaks in the years 2014, 2017, and 2019, with predicted contributions of approximately 0.3, 0.35, and 0.4, respectively, followed by an observable downward trend post-2019 (Fig. 4a). The inter-annual variability in the predicted contributions over the years ranged from around 0.1 in 2010 to a maximum of 0.4 in 2019, subsequently decreasing to approximately 0.2 in 2022.

The predicted areas of acoustic presence extend throughout the study region, with a higher acoustic presence predicted around Kerguelen Island and its plateau, and adjacent to the offshore areas of Madagascar (Fig. 5a). The intra-annual variability in acoustic presence is evident in all four seasons. The highest acoustic presence is observed during austral winter and spring, with predicted contributions ranging between 0.3 and 0.4. In contrast, summer and autumn show significantly lower contributions, ranging between 0 and 0.1 (Fig. 4a). Despite these seasonal fluctuations, the Kerguelen area remains consistently occupied throughout the year, with vocalizations detected in all seasons (Fig. 4a). All predicted DPM maps are available in the supplementary material, from Fig. Supplementary 8 to Fig. Supplementary 20.

### SEIO PBW

SEIO PBW are mainly detected acoustically at locations located east of the observation network, specifically at SWAMS and NEAMS, whereas they are nearly absent from northwestern sites such as MAD and southern locations such as ELAN (Fig. Supplementary 2). The covariates that account for the most significant variance in the model include chlorophyll-a concentration, sea surface temperature, sea ice coverage, GPS coordinates, month, and year. All of these covariates are statistically significant ($p < 0.001$) (Fig. Supplementary 4).

Chlorophyll-a concentration exhibits a notable relationship with acoustic presence. Very low chlorophyll-a concentrations ($log (mean Chla) = -3$ mg/m³) are associated with a very low acoustic presence, while higher chlorophyll-a concentrations ($log (mean Chla) = -0.4$ mg/m³) are associated with a stronger acoustic presence. Between these values, the relationship remains stable (Fig. 4b). Regarding sea surface temperature, a negative relationship with acoustic presence is observed at low temperatures, specifically 0 to 5°C, followed by a rebound and stabilization between 10 and 20°C, and a return to a negative relationship at temperatures exceeding 20°C (Fig. 4b). There is no significant correlation with the percentage of sea ice coverage, except in higher percentages (Fig. 4b). GPS coordinates contribute significantly to spatial variability, with a higher acoustic presence observed in the eastern part of the study area, particularly around the Kerguelen Islands (Fig. 5b).

In terms of temporal variability, there is marked seasonal variation in the acoustic presence. Activity peaks during the austral summer and extends into the autumn (between 0 and 1). The acoustic presence decreases significantly during the winter months, with values dropping to near -2 and remaining low during early spring (-1; Fig. 4b). The inter-annual variability in acoustic presence is evident, with a steady increase observed until 2019 (from approximately -0.7 to 0), followed by a subsequent decline to below -1 in 2022 (Fig. 4b). All predicted DPM maps are available in the supplementary material, from Fig. Supplementary 21 to Fig. Supplementary 33.

### SWIO PBW

SWIO PBW are broadly present within the network, but are quasi-absent from extreme northeastern sites like RTJ and extreme southern sites like ELAN (Fig. Supplementary 3). The highest detection rates occur in the western part of the network, particularly at offshore sites near Madagascar, such as MAD, MADE, and MADW (Fig. Supplementary 3).
<figure><img src="image_9.png" /><figcaption>**Fig. 4.** Functional relationships obtained with generalized additive models between DPM and environmental variables for the three acoustic populations: (a) ANT BW (explained deviance = 21.2%), (b) SEIO PBW (explained deviance = 57%), and (c) SWIO PBW (explained deviance = 37.6%). Solid lines represent the estimated smooth functions and the dashed lines represent the approximate 95% confidence intervals. The x-axis corresponds to the linear predictor and the y-axis to the covariate. Tick marks on the x-axis show the distribution of data points.</figcaption></figure>

The covariates that have the most significant impact on the variance of our model include chlorophyll-a concentration, sea surface temperature, GPS coordinates, month, and year. All of these covariates are statistically significant (p < 0.001)) (Fig. Supplementary 4). Chlorophyll-a concentration shows a positive correlation with acoustic presence only at higher concentrations (log (mean Chla) -1.2 mg/m³) (Fig. 4c). For SST, there is a generally negative relationship with acoustic presence at temperatures near 0°C and a stronger negative
<figure><img src="image_10.png" /><figcaption>**Fig. 5.** Predicted log(DPM) distribution for the three acoustic populations: (a) ANT BW, (b) SEIO PBW, and (c) SWIO PBW for the four seasons, represented by January (summer), April (autumn), July (winter), and November (spring), for the years 2014 and 2022. The red square indicates the Kerguelen Island region, and the blue square corresponds to the area east of Madagascar. The color scale is not the same for all acoustic populations. Predictions for each month and year are provided in **supplementary materials**.</figcaption></figure>
correlation at temperatures exceeding approximately 15°C (Fig. 4c). No significant correlation is observed with the increasing percentage of sea ice coverage (Fig. 4c). The predicted areas of acoustic presence are spread throughout the region, with notably high predictions off Madagascar and around Kerguelen, while the acoustic population is not predicted to be present south of 50°S or in the northeast region (Fig. 5c). Seasonally, the acoustic population exhibits strong peaks in activity during austral summer and autumn (between 0 and 0.5), with a notable peak in 2018 (approximately 0.8). Since 2019, the acoustic population has shown a more stable pattern (around 0.2) (Fig. 4c). All predicted DPM maps are available in the supplementary material, from Fig. Supplementary 34 to Fig. Supplementary 46.

## EOF to identify functional areas

For the three acoustic populations, the first two dimensions of the EOF were selected, which explain more than 60% of the variance for the ANT BW and SWIO PBW acoustic populations, and more than 75% for the SEIO PBW acoustic population. Additional dimensions were excluded from the analysis due to their tendency to represent noise (Fig. Supplementary 6). Each acoustic population exhibits strong spatio-seasonal patterns with succession between a summer/autumn season and a spring/winter season. The cluster analysis highlights the different functional areas predicted for whales during the different seasons for each species (Fig. Supplementary 7).

For the ANT BW subspecies, three temporal clusters and four spatial clusters were predicted (Fig. 6a). Spatial clusters 1 and 4 are predicted for summer and autumn, mainly in the Kerguelen Plateau and higher latitudes. Spatial cluster 2 is predicted for spring and shows limited acoustic activity, concentrated mainly in the southwestern region. Spatial cluster 3 is predicted for winter, predominantly near the coasts or islands within the area.

For SEIO PBW, two temporal clusters and three spatial clusters were predicted (Fig. 6b). Spatial clusters 1 and 2 are predicted for summer and autumn, mainly in the eastern areas surrounding the Kerguelen Plateau. Spatial cluster 3, predicted for winter and spring, suggests that this acoustic population is likely outside the study area.

For SWIO PBW, three temporal clusters and three spatial clusters were predicted (Fig. 6c). In winter, the acoustic population is expected to be mainly located in spatial cluster 3, off Madagascar. In spring, acoustic activity is predicted in clusters 1 and 3, either off Madagascar or beyond the zone. In summer and autumn, the acoustic presence is predicted in clusters 1 and 2, extending from off Madagascar to the Kerguelen Plateau.

## Discussion

Using generalized additive models (GAMs) and empirical orthogonal functions (EOFs), we identified key relationships between acoustic presence and environmental variables. ANT BW acoustic activity was strongly associated with low sea surface temperatures, high chlorophyll-a concentrations, and extensive sea ice cover, primarily during austral winter and spring. In contrast, the acoustic populations of SEIO and SWIO PBW showed peak activity in temperate waters (12-15 °C) with intermediate chlorophyll-a concentrations, revealing different ecological niches. Three critical habitats were identified in the southern Indian Ocean: (1) Madagascar, which serves as a key potential resting area during the austral winter; (2) the Kerguelen Plateau, an essential summer and fall habitat characterized by high primary productivity, which likely supports both feeding and resting activities; and (3) an overlap zone for the SEIO and SWIO PBW acoustic populations in summer, suggesting possible interactions such as hybridization or competition. Our approach combines the flexibility of GAMs with the dimensionality reduction power of EOFs. These methods, coupled with the extensive spatial and temporal scales of the dataset, provide novel insights into the migratory patterns, habitat use, and ecological drivers of blue whale distributions<sup>57,58</sup>.

A principal advantage of our approach is its ability to predict the acoustic presence of whales in unsampled regions, offering a comprehensive understanding of their distribution over vast marine areas. We used the number of positive minutes per day, averaged weekly, metric to effectively capture inter-annual and seasonal variations of blue whale acoustic presence over a dataset with extensive spatial and temporal coverage. The incorporation of EOFs further improved spatio-temporal dimensionality reduction, providing novel insights into critical habitats and seasonal dynamics of blue whale acoustic populations.

The results of our model align with the established knowledge on the ANT BW acoustic population. These cetaceans are predominantly detected during austral winter and spring, with peak acoustic activity correlated with low temperatures, high chlorophyll-a concentrations, and extensive sea ice cover. The spatial predictions of the model, reflecting the temporal variability, corroborate known migratory behaviors, where the ANT BW feed in the southern ocean during the summer and migrate to lower latitudes during the winter months<sup>2,3,5,20</sup>. Furthermore, inter-annual variability in detection patterns is evident, with a notable decline in acoustic presence after 2019. This fluctuation could be linked to changes in environmental conditions, such as variations in sea ice extent and primary productivity, or anthropogenic factors, suggesting dynamic habitat use and potential shifts in blue whale distribution over time.

These findings are consistent with previous research that emphasizes the importance of sea ice and the availability of prey in the distribution of ANT BW. The positive correlation between acoustic detections and sea ice extent supports the results of several studies<sup>29,30</sup>. Furthermore, the prominence of chlorophyll-a as a proxy of prey availability is consistent with observations in similarly productive areas such as the Northeast Pacific and Northern Patagonia<sup>59,60</sup>.

In particular, our results suggest that the eastern waters of Madagascar may serve as a significant resting area for ANT BW during austral winter. These regions, characterized by moderate temperatures and productivity, are mainly visited during this season, indicating their potential role as temporary refuges or resting areas. The Kerguelen Plateau appears to play a prominent role during the austral summer and autumn, where a significant acoustic presence of blue whales is observed. This area may function both as a resting zone before migration
<figure><img src="image_11.png" /><figcaption>**Fig. 6.** Functional areas and associated seasons for the three acoustic populations: (a) ANT BW, (b) SEIO PBW, and (c) SWIO PBW. (Left column) Projection of the standardized temporal loadings (seasons) and the standardized spatial factors (spatial clusters) onto the first two dimensions of the EOF. (Right column) Mapping of spatial clusters using standardized spatial factors to identify functional areas. The same colors are used to represent spatial clusters on the maps (right column) and on the first two dimensions of the EOF (polygons in the left column). In the left column, each season is represented by a square with a specific color. If the loadings of a season are in the same direction as a spatial cluster, then they are correlated and correspond to a couple season x functional area.</figcaption></figure>

to higher latitudes to feed in summer and as a feeding site for non-migrating individuals⁵,⁶¹. Its substantial primary productivity, particularly in krill abundance⁶², supports this hypothesis, highlighting the importance of Kerguelen as a critical habitat for blue whale acoustic populations.

SEIO and SWIO PBW exhibit different spatial distribution patterns throughout the year. However, these acoustic populations were recorded simultaneously during the austral summer and autumn, with maximum activity observed in regions characterized by intermediate temperatures (12–15 °C) and high chlorophyll-a concentrations. This spatial and temporal overlap supports hypotheses from previous studies, which suggest potential interactions between these two acoustic populations, such as hybridization or competition for food resources⁶³. EOF analysis further improves understanding of habitat use by identifying spatio-temporal patterns consistent with those described in other dynamic models⁶⁴,⁶⁵.
In contrast to ANT BW, pygmy blue whales show a minimal association with sea ice, which has limited influence on the predictions in our model. This observation is consistent with studies of SEIO PBW acoustic populations, which do not migrate as far south to feed. Satellite telemetry studies have shown that these whales are not detected south of 42°S, which confirms their tendency to avoid regions with extensive sea ice cover⁶⁶,⁶⁷. The hypothesis of geographical separation within the SEIO PBW acoustic population is supported by recent studies documenting the presence of SEIO PBW in the Kerguelen Plateau and southern Australian waters during the same austral summer and autumn season²¹,²². Acoustic data from Kerguelen and southern Australia overlap seasonally, reinforcing the idea that the Kerguelen Plateau is an important feeding ground for these whales, especially given its high primary productivity. Satellite tracking further supports this, showing that SEIO PBW individuals migrate south from Indonesia to Kerguelen for feeding but also travel to southern Australia, indicating a potential division within the population into two distinct migratory routes⁶⁶,⁶⁷. These observations suggest that the SEIO PBW population may follow different migratory routes, with migrations to Kerguelen driven by foraging rather than reproduction, as these whales are known to reproduce in Indonesian waters. Although Kerguelen's abundant krill⁶² makes it a critical food source, separation into distinct migratory routes and ecological zones may not be fixed.

Similarly, SWIO PBWs disperse as they migrate across the southwestern Indian Ocean. Their range extends from the south of Marion and Prince Edward Islands, where acoustic detections have been reported between December and April⁶⁸ to the south of Amsterdam Island in the east, with detections occurring from December to July⁷. The region surrounding Prince Edward Island has been identified as a significant feeding ground during this period⁶⁸. The concomitant recording of SWIO PBW songs near Kerguelen in our study posits the potential for Kerguelen to also function as a feeding site. This suggests that there may be a spatial dispersion within the SWIO PBW population during summer - mid-fall, with individuals following different migratory pathways while exploiting similar ecological resources. Northward migration begins in April, with individuals traveling toward Madagascar, where the migratory route splits: one branch passes through the Mozambique Channel²⁴, while the other skirts offshore around Réunion and Mauritius⁷. The northernmost acoustic detections have been recorded near Mayotte between May and July²⁵. However, all acoustic monitoring stations across the Indian Ocean report an absence of calls from August to October⁶,⁷,²⁴,⁶⁸, which may indicate a shift in vocal behavior or a seasonal movement of the population to a more northerly, equatorial region that remains unsurveyed. Migration to the central and eastern tropical Indian Ocean is unlikely, as no calls have been detected around the Seychelles or Diego Garcia⁶⁹,⁷⁰. During their return south (October to January), SWIO PBW calls are not detected by offshore acoustic stations in the western Indian Ocean⁷. However, they are consistently detected along the western and eastern coasts of Madagascar⁶,²⁴,²⁵, suggesting that coastal waters provide a more favorable habitat during this part of the migratory cycle.

Our results further advance current knowledge by emphasizing the role of the Kerguelen Plateau for pygmy blue whale acoustic populations as a shared habitat during austral summer and autumn. This area, with its high primary productivity and abundance of krills⁶², is likely to support feeding and resting activities for these acoustic populations. Identifying key habitats underscores the ecological importance of these areas for the conservation of blue whales. The region serves as a biodiversity hotspot, supporting not only blue whales but also other predators such as elephant seals and penguins⁷¹,⁷².

This study has certain limitations to consider. First, ambient noise, which can affect call detection, was not explicitly taken into account in our model. High levels of noise from natural or anthropogenic sources could mask whale calls and lead to underestimations of acoustic activity⁵⁹,⁶⁴. Including ambient noise as a covariate in future studies would improve the accuracy of the model. To address these limitations, we estimate a detection range of 50 km around the hydrophones to account for sound propagation and minimize bias. An estimated detection range of 50 km was determined as a conservative figure, reflecting the substantial variability in the detection ranges documented in the literature. These ranges exhibit substantial variation, ranging from 10 km to more than 700 km, depending on factors such as geographical region, seasonal variations, and source levels between acoustic populations⁶⁸,⁷³,⁷⁴.

Secondly, this analysis intentionally excluded D-calls, which represent short acoustic signals typically associated with foraging activities and are believed to aid in prey detection or coordination during feeding⁹. These vocalizations were omitted because of their inability to be reliably attributed to a distinct acoustic population. However, the inclusion of D-calls in future research could provide deeper insights into the various behavioral contexts of blue whale vocalization, especially in highly productive areas such as the Kerguelen Plateau. This proposition is corroborated by Torterotot et al. (2023)⁷⁵, who reported high occurrence of D-calls near the Kerguelen Plateau. In contrast, the blue whale songs examined in this study are predominantly associated with reproductive activities, similar to those observed in humpback whales¹³,⁷⁶. However, their presence in feeding regions⁷⁷,⁷⁸ implies that they may perform additional functions.

Third, we used chlorophyll-a concentration as a proxy for prey availability, which may not directly reflect krill biomass, the primary food source for blue whales. This limitation stems from the difficulty of obtaining direct measurements of the biomass of krill in such vast regions and time frames. Future work should incorporate direct measurements of krill biomass, such as acoustic surveys or satellite data, to better capture prey dynamics and refine habitat models.

In addition, we focus on broad-scale variables such as SST, chlorophyll-a, and sea ice extent. Although this approach is ecologically relevant on large scales, finer-scale factors such as nutrient upwelling or mesoscale eddies could improve local habitat models and provide deeper ecological insights.

In conclusion, this study sheds light on the complex relationships between the acoustic populations of blue whales and their environment. By identifying functional zones shaped by environmental factors, it establishes a foundation for targeted conservation efforts to preserve these critical marine ecosystems.
## Data availability

The environmental data used in this study are publicly available from dedicated databases: Sea Surface Temperature (SST): Data were obtained from the Copernicus Marine global ocean reanalysis (DOI: 10.48670/moi-00021), covering the period 1992 to 2024 with a spatial resolution of 1/12° and daily temporal resolution. Chlorophyll-a concentration: Data were sourced from Mercator-Ocean (DOI: 10.48670/moi-00019), available from January 1993 to February 2024 at a 1/4° resolution. Sea Ice Concentration: Data were obtained from the NOAA/NSIDC Climate Data Record, Version 4 (DOI: 10.7265/efmz-2t65), based on passive microwave measurements and available at a 25 km × 25 km resolution with daily and monthly temporal scales. These environmental data were used to generate Fig. 3 in the manuscript and were included as explanatory variables in the models used throughout the study. The acoustic data used in this study (DOI: 10.18142/229) were used to generate Fig. 2 and served as the response variable in our model. The raw acoustic data is not publicly available, but the detection data can be obtained from Jean-Yves Royer (jean-yves.royer@univ-brest.fr) upon reasonable request. Additional data sets were not generated or analyzed beyond those specified.

Received: 3 February 2025; Accepted: 16 May 2025

Published online: 02 July 2025

## References

1. IWC. Chairman's report of the thirty-first annual meeting. *Rep. int. Whal. Commn.* vol. 30, 25-41 (1980).

2. Samaran, F., Adam, O. & Guinet, C. Discovery of a mid-latitude sympatric area for two southern hemisphere blue whale subspecies. *Endanger. Species Res.* **12**, 157-165 (2010).

3. Samaran, F. et al. Seasonal and geographic variation of southern blue whale subspecies in the Indian ocean. *PLoS ONE* **8**, e71561 (2013).

4. Stafford, K. M., Chapp, E., Bohnenstiel, D. R. & Tolstoy, M. Seasonal detection of three types of “pygmy” blue whale calls in the Indian ocean. *Mar. Mamm. Sci.* **27**, 828–840 (2011).

5. Leroy, E. C., Samaran, F., Bonnel, J. & Royer, J.-Y. Seasonal and diel vocalization patterns of Antarctic blue whale (*Balaenoptera musculus intermedia*) in the southern Indian ocean: A multi-year and multi-site study. *PLoS ONE* **11**, e0163587 (2016).

6. Dréo, R., Bouffaut, L., Leroy, E., Barruol, G. & Samaran, F. Baleen whale distribution and seasonal occurrence revealed by an ocean bottom seismometer network in the western Indian ocean. *Deep Sea Res. Part II* **161**, 132-144 (2019).

7. Torterotot, M., Samaran, F., Stafford, K. M. & Royer, J.-Y. Distribution of blue whale populations in the southern Indian ocean based on a decade of acoustic monitoring. *Deep Sea Res. Part II* **179**, 104874 (2020).

8. Cummings, W. C. & Thompson, P. O. Underwater sounds from the blue whale, *Balaenoptera musculus*. *J. Acoust. Soc. Am.* **50**, 1193-1198 (1971).

9. Mellinger, D. K. & Clark, C. W. Blue whale (*Balaenoptera musculus*) sounds from the north Atlantic. *J. Acoust. Soc. Am.* **114**, 1108-1119 (2003).

10. Thompson, P. O., Findley, L. T., Vidal, O. & Cummings, W. C. Underwater sounds of blue whales, *Balaenoptera musculus*, in the gulf of California, Mexico. *Mar. Mamm. Sci.* **12**, 288-293 (1996).

11. Thode, A. M., D'Spain, G. L. & Kuperman, W. A. Physical constraints on blue whale sound production mechanisms derived from matched-field processing. *J. Acoust. Soc. Am.* **105**, 993-993 (1999).

12. Lewis, L. A. et al. Context-dependent variability in blue whale acoustic behaviour. *R. Soc. Open Sci.* **5**, 180241 (2018).

13. Oleson, E. M., Calambokidis, J., Barlow, J. & Hildebrand, J. A. Blue whale visual and acoustic encounter rates in the southern California bight. *Mar. Mamm. Sci.* **23**, 574-597 (2007).

14. McDonald, M. A., Mesnick, S. L. & Hildebrand, J. A. Biogeographic characterization of blue whale song worldwide: Using song to identify populations. *J. Cetac. Res. Manag.* **8**, 55-65 (2006).

15. Leroy, E. C., Samaran, F., Stafford, K. M., Bonnel, J. & Royer, J.-Y. Broad-scale study of the seasonal and geographic occurrence of blue and fin whales in the southern Indian ocean. *Endanger. Species Res.* **37**, 289-300 (2018).

16. Leroy, E. C., Royer, J.-Y., Bonnel, J. & Samaran, F. Long-term and seasonal changes of large whale call frequency in the southern Indian ocean. *J. Geophys. Res. Oceans* **123**, 8568-8580 (2018).

17. Miller, B. S. et al. Antarctic sonobuoy surveys for blue whales from 2006-2021 reveal contemporary distribution, changes over time, and paths to further our understanding of their distribution and biology. *Front. Mar. Sci.* **11**, 1324816 (2024).

18. Davis, G. E. et al. Exploring movement patterns and changing distributions of baleen whales in the western north Atlantic using a decade of passive acoustic data. *Glob. Change Biol.* **26**, 4812–4840 (2020).

19. Ichihara, T. The pygmy blue whale, *Balaenoptera musculus brevicauda*, a new subspecies from the Antarctic. *Whales, dolphins, and porpoises* **79–111** (1966).

20. Stafford, K. M. et al. Antarctic-type blue whale calls recorded at low latitudes in the Indian and eastern Pacific oceans. *Deep Sea Res. Part I* **51**, 1337-1346 (2004).

21. Gavrilov, A. N. & McCauley, R. D. Acoustic detection and long-term monitoring of pygmy blue whales over the continental slope in southwest Australia. *J. Acoust. Soc. Am.* **134**, 2505-2513 (2013).

22. McCauley, R. D., Gavrilov, A. N., Jolliffe, C. D., Ward, R. & Gill, P. C. Pygmy blue and Antarctic blue whale presence, distribution and population parameters in southern Australia based on passive acoustics. *Deep Sea Res. Part II* **157**, 154-168 (2018).

23. Ljungblad, D., Clark, C. W. & Shimada, H. Sc/49/sh17 a comparison of sounds attributed to pygmy blue whales (*Balaenoptera musculus brevicauda*) recorded south of the Madagascar plateau and those attributed to 'true' blue whales (*Balaenoptera musculus*) recorded off Antarctica. *Rep.-Int. Whaling Commiss.* **48**, 439-442 (1998).

24. Cerchio, S., Rasoloarijao, T., Mueller-Brennan, B. & Cholewiak, D. Migratory movements and acoustic activity of blue whales (*Balaenoptera musculus*) and fin whales (*Balaenoptera physalus*) in the eastern Mozambique channel. *Technical report* (2024).

25. Dréo, R. et al. Singing around the volcano: Detecting baleen whales in the Mozambique channel based on their song rhythms, from seismic and hydroacoustic data. *J. Acoust. Soc. Am.* (2025).

26. Virgili, A. et al. Combining multiple visual surveys to model the habitat of deep-diving cetaceans at the basin scale: Large-scale modelling of deep-diving cetacean habitats. *Glob. Ecol. Biogeogr.* **28**, 300–314 (2019).

27. Lambert, C., Mannocci, L., Lehodey, P. & Ridoux, V. Predicting cetacean habitats from their energetic needs and the distribution of their prey in two contrasted tropical regions. *PLoS ONE* **9**, e105958 (2014).

28. Nicol, S., Worby, A. & Leaper, R. Changes in the Antarctic sea ice ecosystem: Potential effects on krill and baleen whales. *Mar. Freshw. Res.* **59**, 361-382 (2008).

29. Shabangu, F. W., Yemane, D., Stafford, K. M., Ensor, P. & Findlay, K. P. Modelling the effects of environmental conditions on the acoustic occurrence and behaviour of Antarctic blue whales. *PLoS ONE* **12**, e0172705 (2017).

30. Shabangu, F. W. et al. Seasonal occurrence and diel calling behaviour of Antarctic blue whales and fin whales in relation to environmental conditions off the west coast of South Africa. *J. Mar. Syst.* **190**, 25-39 (2019).
31. Belkin, I. M. & Gordon, A. L. Southern ocean fronts from the Greenwich meridian to Tasmania. *J. Geophys. Res. Oceans* **101**, 3675–3696 (1996).

32. Orsi, A. H., Whitworth, T. III. & Nowlin, W. D. Jr. On the meridional extent and fronts of the Antarctic circumpolar current. *Deep Sea Res. Part I* **42**, 641–673 (1995).

33. Behrenfeld, M. J. et al. Climate-driven trends in contemporary ocean productivity. *Nature* **444**, 752–755 (2006).

34. Smith, W. H. & Sandwell, D. T. Global sea floor topography from satellite altimetry and ship depth soundings. *Science* **277**, 1956–1962 (1997).

35. Moore, J. K. & Abbott, M. R. Surface chlorophyll concentrations in relation to the Antarctic polar front: seasonal and spatial patterns from satellite observations. *J. Mar. Syst.* **37**, 69–86 (2002).

36. Royer, J.-Y. Oha-sis-bio-observatoire hydroacoustique, https://doi.org/10.18142/229 (2009).

37. Lutron, X. An introduction to underwater acoustics (2002).

38. d'Eu, J.-F., Royer, J.-Y. & Perrot, J. Long-term autonomous hydrophones for large-scale hydroacoustics monitoring of the oceans. In *2012 Oceans-Yesosu*, 1–6 (IEEE, 2012).

39. Socheleau, F.-X. & Samaran, F. Detection of mysticete calls: a sparse representation-based approach. *Technical report* (2018).

40. Guilment, T., Socheleau, F.-X., Pastor, D. & Vallez, S. Sparse representation-based classification of mysticete calls. *J. Acoust. Soc. Am.* **144**, 1550–1563 (2018).

41. Torterotot, M., Royer, J.-Y. & Samaran, F. Detection strategy for long-term acoustic monitoring of blue whale stereotyped and non-stereotyped calls in the southern indian ocean. In *OCEANS 2019-Marseille*, 1–10 (IEEE, 2019).

42. Godet, C. et al. Matching zooplankton abundance and environment in the south Indian ocean and southern ocean. *Deep Sea Res. Part I* **163**, 103347 (2020).

43. Širović, A. et al. Seasonality of blue and fin whale calls and the influence of sea ice in the western Antarctic peninsula. *Deep Sea Res. Part II* **51**, 2327–2344 (2004).

44. Meier, W. N. et al. Arctic sea ice in transformation: A review of recent observed changes and impacts on biology and human activity. *Rev. Geophys.* **52**, 185–217 (2014).

45. Meier, W. N., Fetterer, F., Windnagel, A. K. & Stewart, J. S. Noaa/nsidc climate data record of passive microwave sea ice concentration. (g02202, version 4). (National Snow and Ice Data Center, 2020).

46. Wood, S. N. *Generalized Additive Models: An Introduction with R* (Chapman and Hall/CRC, 2017).

47. Arcuti, S. et al. Spatio-temporal modelling of zero-inflated deep-sea shrimp data by tweedie generalized additive. *Statistica (Bologna)* **73**, 87–101 (2013).

48. Cressie, N., Sainsbury-Dale, M. & Zammit-Mangion, A. Basis-function models in spatial statistics. *Annu. Rev. Stat. Appl.* **9**, 373–400 (2022).

49. Akaike, H. A new look at the statistical model identification. *IEEE Trans. Autom. Control* **19**, 716–723 (1974).

50. Wood, S. *Generalized Additive Models: An Introduction with R* 2 edn (Chapman and Hall/CRC, 2017).

51. R Core Team. *R: A Language and Environment for Statistical Computing*. (R Foundation for Statistical Computing, 2021).

52. Thorson, J. T. et al. Empirical orthogonal function regression: Linking population biology to spatial varying environmental conditions using climate projections. *Glob. Change Biol.* **26**, 4638–4649 (2020).

53. Alglave, B. et al. Investigating fish reproduction phenology and essential habitats by identifying the main spatio-temporal patterns of fish distribution. *ICES J. Mar. Sci.* **81**, 1563–1574 (2024).

54. Thorson, J. T., Ciannelli, L. & Litzow, M. A. Defining indices of ecosystem variability using biological samples of fish communities: a generalization of empirical orthogonal functions. *Prog. Oceanogr.* **181**, 102244 (2020).

55. Alglave, B., Dufée, B., Obakrim, S. & Thorson, J. T. Empirical Orthogonal Functions for ecology (2024). Preprint.

56. Lé, S., Josse, J. & Husson, F. FactoMineR: A package for multivariate analysis. *J. Stat. Softw.* **25**, 1–18, https://doi.org/10.18637/jss.v025.i01 (2008).

57. Monnahan, C. C., Branch, T. A., Stafford, K. M., Ivashchenko, Y. V. & Oleson, E. M. Estimating historical eastern north pacific blue whale catches using spatial calling patterns. *PLoS ONE* **9**, e98974 (2014).

58. Anderson, S. C., Ward, E. J., English, P. A. & Barnett, L. A. sdmtmb: an r package for fast, flexible, and user-friendly generalized linear mixed effects models with spatial and spatiotemporal random fields. *BioRxiv* 2022–03 (2022).

59. Burtenshaw, J. C. et al. Acoustic and satellite remote sensing of blue whale seasonality and habitat in the northeast pacific. *Deep Sea Res. Part II* **51**, 967–986 (2004).

60. Buchanan, S. J. & Quiñones, R. A. First insights into the oceanographic characteristics of a blue whale feeding ground in northern Patagonia, Chile. *Mar. Ecol. Prog. Ser.* **554**, 183–199 (2016).

61. Branch, T. A. et al. Past and present distribution, densities and movements of blue whales *Balaenoptera musculus* in the southern hemisphere and northern Indian ocean. *Mammal Rev.* **37**, 116–175 (2007).

62. Hunt, B. & Swadling, K. Macrozooplankton and micronekton community structure and diel vertical migration in the heard island region, central Kerguelen plateau. *J. Mar. Syst.* **221**, 103575 (2021).

63. Attard, C. R. et al. Hybridization of southern hemisphere blue whale subspecies and a sympatric area off Antarctica: Impacts of whaling or climate change. *Mol. Ecol.* **21**, 5715–5727 (2012).

64. Barlow, D. R. & Torres, L. G. Planning ahead: dynamic models forecast blue whale distribution with applications for spatial management. *J. Appl. Ecol.* **58**, 2493–2504 (2021).

65. Abrahms, B. et al. Memory and resource tracking drive blue whale migrations. *Proc. Natl. Acad. Sci.* **116**, 5582–5587 (2019).

66. Double, M. C. et al. Migratory movements of pygmy blue whales (*Balaenoptera musculus brevicauda*) between Australia and Indonesia as revealed by satellite telemetry. *PLoS ONE* **9**, e93578 (2014).

67. Mustika, P. L. K. *et al*. The first record of the southbound movements of satellite-tagged pygmy blue whales (*B. m. brevicauda*) from Savu Sea (Indonesia) to the subantarctic waters. *Mar. Mamm. Sci.* e13167 (2024).

68. Shabangu, F. W. et al. Diverse baleen whale acoustic occurrence around two sub-antarctic islands: A tale of residents and visitors. *Sci. Rep.* **14**, 21663 (2024).

69. Stafford, K. M. et al. Acoustic detections and sightings of blue whales *Balaenoptera musculus* in the Seychelles, western tropical Indian ocean (2020–2022). *Endanger. Species Res.* **52**, 203–208 (2023).

70. Leroy, E. C., Royer, J.-Y., Alling, A., Maslen, B. & Rogers, T. L. Multiple pygmy blue whale acoustic populations in the Indian ocean: Whale song identifies a possible new population. *Sci. Rep.* **11**, 8762 (2021).

71. Toole, M. O., Guinet, C., Lea, M.-A. & Hindell, M. A. Marine predators and phytoplankton: how elephant seals use the recurrent Kerguelen plume. *Mar. Ecol. Prog. Ser.* **581**, 215–227 (2017).

72. Hindell, M. A. *et al*. Foraging habitats of top predators, and areas of ecological significance, on the Kerguelen plateau. *The Kerguelen Plateau Mar. Ecosyst. Fish.* 203–215 (2011).

73. Samaran, F., Guinet, C., Adam, O., Motsch, J.-F. & Cansi, Y. Source level estimation of two blue whale subspecies in southwestern Indian ocean. *J. Acoust. Soc. Am.* **127**, 3800–3808 (2010).

74. Širović, A., Hildebrand, J. A. & Wiggins, S. M. Blue and fin whale call source levels and propagation range in the southern ocean. *J. Acoust. Soc. Am.* **122**, 1208–1215 (2007).

75. Torterotot, M., Samaran, F. & Royer, J.-Y. Long-term acoustic monitoring of nonstereotyped blue whale calls in the southern indian ocean. *Mar. Mamm. Sci.* **39**, 594–610 (2023).

76. McDonald, M. A., Calambokidis, J., Teranishi, A. M. & Hildebrand, J. A. The acoustic calls of blue whales off California with gender data. *J. Acoust. Soc. Am.* **109**, 1728–1735 (2001).
## Acknowledgements

This project was funded by the ANR (grant ESDIR2023). The authors wish to thank the captains and crews of the Marion Dufresne R/V for the successful deployments and recoveries of the instruments of the OHA-SIS-BIO array. The cruises were funded by the French Polar Institute and the French Oceanographic Fleet with additional support from Institut National des Sciences de l'Univers Centre National de la Recherche Scientifique.

## Author contributions

Mathilde Michel: Conceptualization, Methodology, Formal Analysis, Original Draft Writing, Visualization. Baptiste Algave: Conceptualization, Methodology, Validation, Formal Analysis, Writing-Review & Editing. Maxime Olmos: Conceptualization, Methodology, Validation, Formal Analysis, Writing-Review & Editing. Maëlle Torterotot: Validation, Formal Analysis, Resources, Data Curation, Writing-Review & Editing. Auriane Virgili: Validation, Writing-Review & Editing. Salomé Martin-Marin: Formal Analysis, Resources. Jean-Yves Royer: Writing-Review & Editing, Project Administration, Funding Acquisition. Flore Samaran: Conceptualization, Validation, Writing-Review & Editing, Project Administration, Funding Acquisition.

## Declarations

### Competing interests

The authors declare no competing interests.

### Additional information

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-025-02941-9.

**Correspondence** and requests for materials should be addressed to M.M.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2025

<footnote n="77">77. Širović, A., Hildebrand, J. A., Wiggins, S. M. & Thiele, D. Blue and fin whale acoustic presence around Antarctica during 2003 and 2004. *Mar. Mamm. Sci.* **25**, 125–136 (2009).</footnote>

<footnote n="78">78. Thomisch, K., et al. Spatio-temporal patterns in acoustic presence and distribution of Antarctic blue whales *Balaenoptera musculus* intermedia in the Weddell sea. *Endanger. Species Res.* **30**, 239–253 (2016).</footnote>