Quantifying Soundscapes in the Ross Sea, Antarctica Using Long-Term Autonomous Hydroacoustic Monitoring Systems

Deployment of long-term, continuously recording passive-acoustic sensors in the ocean can provide insights into sound sources related to ocean dynamics, air–sea interactions, and biologic and human activities, all which contribute to shaping ocean soundscapes. In the polar regions, the changing ocean climate likely contributes to seasonal and long-term variation in cryogenic sounds, adding to the complexity of these soundscapes. The Korea Polar Research Institute and the U.S. National Oceanic and Atmospheric Administration have jointly operated two arrays of autonomous underwater hydrophones in the Southern Ocean, one in the Terra Nova Bay Polynya (TNBP) during December 2015–January 2019 and the other in the Balleny Islands (BI) region during January 2015–March 2016, to monitor changes in ocean soundscapes. In the BI region, we found distinct seasonal variations in the cryogenic signals that were attributed to collisions and thermal/mechanical fracturing of the surface sea ice. This is consistent with sea-ice patterns due to annual freeze–thaw cycles, which are not clearly observed in TNBP, where frequent blowing out of sea ice by katabatic winds and icequakes from nearby ice shelves generate strong noise even in austral winters. Another advantage of passive acoustic recordings is that they provide opportunities to measure biodiversity from classifying spectral characteristics of marine mammals: we identified 1. Leopard seals (Hydrurga leptonyx; 200–400 Hz), most abundant in the BI region and TNBP in December; 2. Antarctic blue whales (Balaenoptera musculus; distinctive vocalization at 18 and 27 Hz), strong signals in austral winter and fall in the BI region and TNBP; 3. Fin whales (B. physalus; fundamental frequency in the 15–28 Hz and overtones at 80 and 90 Hz), maximum presence in the BI region during the austral summer and spring months; 4. Antarctic minke whales (B. bonaerensis; 100–200 Hz), strongest signals from June to August in the BI region; 5. Humpback whales in TNBP; 6. Unidentified whales (long-duration downsweeping from 75 to 62 Hz), detected in TNBP. Long-term soundscape monitoring can help understand the spatiotemporal changes in the Southern Ocean and cryosphere and provide a means of assessing the status and trends of biodiversity in the Ross Sea Region Marine Protected Area.


INTRODUCTION
The ocean soundscape is a combination of underwater sounds from various sources corresponding to a particular location and time, including from marine animals, tectonic activities, ocean dynamics, air-sea interactions, and human activities. The basic concept of underwater sound sources and their acoustic characteristics have been developed in studies on ocean ambient noise (Knudsen et al., 1948;Wenz, 1962). Knudsen et al. (1948) classified the main sources into water motion, marine life, shipping, and other man-made sources, and investigated the ambient noise characteristics (spectra and pressure levels) from these sources. The term "soundscape" was coined by Schafer (1969), who established the World Soundscape Project (WSP) in the late 1960s-early 1970s. These early soundscape studies mainly focused on the ecological field and various interdisciplinary studies (Matsinos et al., 2016).
Listening to sounds and measuring the ocean soundscapes via passive acoustic monitoring (PAM) in the ocean helps understand the physical and biological ocean environment and their interactions. Long-term PAM has revealed that anthropogenic activity (e.g., shipping and seismic surveys) and geophysical noise from stronger winds and icequakes increase ocean ambient noise in the Northern Pacific Ocean (Andrew et al., 2002(Andrew et al., , 2011McDonald et al., 2006;Chapman and Price, 2011;Frisk, 2012), Atlantic and Arctic oceans (Moore et al., 2012), Indian Ocean (Miksis-Olds et al., 2013), and the Southern Ocean (Matsumoto et al., 2014;Haver et al., 2017). Increased ocean ambient noise has a negative influence on marine mammal communication and ecological behaviors (Edds-Walton, 1997;Clark et al., 2009;Castellote et al., 2012).
Despite the difficult approach and risks in the recovery of moorings due to sea ice, thousands of hours of underwater acoustic recordings have been made throughout the Southern Ocean south of 60 • S since 2001 (Miller et al., 2021). One of the main objectives of this study is to examine the distribution and behavior of endangered Antarctic blue whales (Balaenoptera musculus intermedia) and fin whales (B. physalus) Van Opzeeland et al., 2013). Another objective is to monitor tectonic activities (Dziak et al., 2010) and glacial ice-ocean boundary processes (Pettit, 2012;Pettit et al., 2012;Glowacki et al., 2015;Dziak et al., 2019;Glowacki and Deane, 2020).
In this paper, we give an overview the Southern Ocean soundscape in two different regions of the Western Ross Sea, Antarctica-the Balleny Islands (BI) region, and Terra Nova Bay Polynya (TNBP)-located ∼800 km from each other along the eastern Antarctic coastline. We (1) present various biotic and abiotic sound sources in the BI and TNBP regions in the Ross Sea Marine Protected Area (MPA), (2) show influences of geophysical components (wind speed and sea ice concentration) on the ocean ambient noise and marine mammal vocal activity, and (3) discuss how the geographical differences change the soundscapes in these regions.

Autonomous Underwater Acoustic Recording System and Mooring
The Korea Polar Research Institute and the U.S. National Oceanic and Atmospheric Administration have jointly monitored the soundscape of the Southern Ocean, specifically in the BI region (January 2015-March 2016) ( Figure 1A) and the TNBP (December 2015-January 2019) ( Figure 1B) in the Ross Sea.
The BI is a volcanic island chain trending northwestsoutheast, and an array of five hydrophones was deployed parallel and to the east of the islands by icebreaker R/V Araon in January 2015 and recovered in March 2016. Dziak et al. (2017) found various sound sources (volcano-seismic, cryogenic, and biogenic) in the long-term spectrogram of a hydrophone deployed 400 km north off the coast of the Antarctic continent and 150 km away from the BI. In this study, we analyze the same hydrophone recordings under the identifier BAL15.
In the TNBP, the south and north arrays have been operated along the northeastern edge of the Drygalski Ice Tongue and northern sea-ice margin of the polynya, respectively. The TNB south array was first deployed in December 2015 (Dziak et al., 2019), and there were two recovery/redeployments in February 2017 and March 2018. Final recovery was done in January 2019 by icebreaker R/V Araon. The TNB north array was first deployed in February 2017, recovery redeployments were implemented in March 2018, and the recovery was done in January 2019. The locations were changed in every deployment due to the sea ice distribution. Table 1 presents the locations and recording time periods.
The PAM system we used is called Autonomous Underwater Hydrophone (AUH) (Fox et al., 2001). The AUH package deployed in the BI region and TNBP composed of a single ceramic hydrophone, with a filter/amplifier, clock, and a lowpower processor called CF2 from Persistor Instrument, which were all powered by an internal battery pack. The hydrophone model was ITC-1032, which is omnidirectional with a nominal sensitivity of -192 dB re 1 V µPa −1 . The instrument continuously  records at a sampling rate of 1 kHz with 16-bit resolution. The pre-amplifier has an eight-pole anti-aliasing filter at 400 Hz and has a filter curve to equalize the spectrum against typical ocean noise over the pass-band (Dziak et al., 2019). In the BI region, a microprocessor controlled, temperature-correcting crystal oscillator with an average time drift of 1.95 s year −1 was deployed for1 year , and a low-power cesium atomic clock (low power) with an average time drift of ∼0.1 s year −1 was used in the TNB array. The electronics were housed in both titanium (for the BI) and fiberglass composite pressure cases (for TNB), and attached to a standard oceanographic mooring with an anchor, acoustic release, mooring line, and syntactic foam float. The moorings were designed to place the sensor at a depth of 400 m to minimize the chance of iceberg collision. To address the detection ranges of the recorders in the Southern Ocean, we estimate transmission loss of 25 and 400 Hz waves using KRAKEN (Porter, 1992) propagation model for a source at 50 m depth, a receiver at 400 m depth, and sound profiles in the TNBP (Supplementary Figure 1). It presents that some of the detected low-frequency signals (e.g., Antarctic blue and fin whale calls with source level of ∼189 dB over 25-29 Hz and 15-28 Hz, respectively; Širović et al., 2007) could originate from 100s of km away from the recorders, whereas sounds around 400 Hz could have relatively local origins.

Spectral Analysis
All data processing and measurements were performed using MATLAB, 2020b. We calculated the hourly/daily averaged spectrum to (1) characterize the seasonal and regional variations in the long-term spectrograms (Figure 2), (2) quantify the marine mammal contribution to the spectra, and (3) present their diel patterns. Fast Fourier transforms (FFTs) were performed with 1 Hz resolution for a 0% overlapped "Hanning" window of length 1 s (1,000 samples), and the hourly/daily averaged spectral amplitudes were calculated for the entire deployment period. Sources of repeated signals in the long-term spectrogram are identified by comparing with the characteristics of various sounds in the Southern Ocean reported in previous studies(e.g., Dziak et al., 2015;Erbe et al., 2017;Menze et al., 2017). Short-duration spectrograms were obtained to verity the exact sound sources by performing FFTs with 1 Hz resolution for a 50% overlap "Hanning" window of 1 s without time averaging (Supplementary Figure 2). To identify inter-annual and seasonal characteristics, kernel-smoothed (bin width 0.5) spectral probability densities (SPDs) were estimated from the daily spectra at equally spaced sound levels between 50 and 120 dB (Figure 3). Marine mammal contributions (MMC) in the power spectral density (PSD MMC ) are calculated for the hourly averaged sound spectra using the method adopted by Menze et al. (2017). We subtract the estimated spectrum without MMC from the measured spectrum. The spectrum without MMC is interpolated by curve fitting with a power function PSD dB = af b + c (where f is the frequency and a, b, and c are constants determined by fitting). The core frequency bands for the subtracted and fitted frequency bands are chosen from the BI and TNB data that reflect their own frequency response and ambient noise characteristics ( Table 2 and Supplementary Figure 3). The PSD MMC value is estimated only when the signal to the background noise ratio of MMC is higher than 1 dB (Menze et al., 2017).

Sea Ice Concentration and Wind Speed Data
Sea-ice concentration data were obtained from the Arctic Radiation and Turbulence Interaction Study with a grid spacing of 3.125 km resolved from the Advanced Microwave Scanning Radiometer 2 (AMSR-2) dataset (Spreen et al., 2008). For the BI region, we averaged the values within a 50 km radius of the BAL15 mooring location (shaded circle in Figure 1C). This criterion was used in a previous soundscape study of the Southern Ocean (Menze et al., 2017). For TNBP, the daily averaged seaice concentration was calculated in a rectangular area with coordinates 75.5-74.7 • S latitude and 164-166.5 • E longitude (shaded area in Figure 1D). TNBP is a coastal latent heat polynya in the western sector of the Ross Sea, Antarctica. The Drygalski Ice Tongue, located at the southern boundary of the polynya, blocks sea ice moving in from the south, and strong katabatic winds from the Nansen Ice Shelf (NIS) frequently push sea ice away from the coast and create open waters even in the winter season (Kurtz and Bromwich, 1985;Fusco et al., 2002;Sansiviero et al., 2017;Stevens et al., 2017). Considering the high mobility of sea ice in the polynya, we set the masking area to cover all hydrophones and surrounding oceans in the polynya. In addition to sea ice concentration, the monthly sea-ice extent was obtained from the National Snow and Ice Data Center (Fetterer et al., 2002) and was used to estimate the distance between the BAL15 station and the nearest open sea. The daily wind speed in TNB was calculated from meteorological observations obtained using the Automatic Synoptic Observation System (ASOS) at Jang Bogo Station (74 • 37.4 S, 164 • 13.7 E) from December 2015 to December 2018. The daily wind speed in the BI region was extracted from the European Center for Medium Range Weather Forecasts (ECMWF) climate reanalysis dataset ERA5 (Hersbach et al., 2020).

RESULTS AND DISCUSSION
Ambient Sound Characteristics in the Balleny Islands Region and Terra Nova Bay Polynya We can identify the dominant sound sources composing the soundscapes in a region by finding repeated signals in the spectrogram. The broadband (10-400 Hz) vertical lines in all the three spectrograms are icequakes, cryogenic signals originating from the impact and breakup of surface sea ice, ice shelf, and icebergs (MacAyeal et al., 2008;Dziak et al., 2010Dziak et al., , 2015Dziak et al., , 2019. The horizontal peaks in the 15-30 Hz band generally associated with fin and/or Antarctic blue whale calls (Watkins et al., 1987;Ljungblad et al., 1998) are shown in both regions. However, the horizontal lines in the 80-90 Hz band, usually from fin whale calls (Širović et al., 2004), are clearly shown only in the BI region, not in the TNB records. The vocalizations of leopard seals (Rogers et al., 1996) are apparent in continuous 200-400 Hz broadband and strongest in 300-370 Hz from October to December in the BI region and TNBP. Strong anthropogenic sounds continued for several days in May of 2017.
Supplementary Figure 2 shows short-term spectrograms detailing various sources in the BI region and TNBP. Supplementary Figure 2A shows typical earthquake and icequake signals recorded in the BI region. Acoustic waves generated from earthquakes are referred to as T-waves, characterized by low-frequency (<100 Hz) wave trains that continue for tens to hundreds of seconds. Icequakes are impulsive, last for a short duration, and have a broadband frequency content. Supplementary Figure 2B is a spectrogram in 1-100 Hz band to scrutinize whale vocalization in the BI region. Lots of short-duration, narrow-band signals exist at 80 Hz, but the frequency is lower than typical upper part of fin whale calls in other part of the Southern Ocean: 89 Hz in the West Antarctic Peninsula in 2002 (Širović et al., 2004), 98 Hz in the Atlantic sector in 2009-2010 (Menze et al., 2017). and 99 Hz in the East Antarctica in 2003 (Širović et al., 2009) shows "bioduck" sounds from Antarctic minke whale (Risch et al., 2014) in the BI region in the frequency band of 100-250 Hz. In Supplementary Figure 2D, leopard seal vocalizations with  Nevertheless, the ranges of core frequency and fitting frequency bands are differently adapted by reflecting the distorted frequency response of the data (Supplementary Figure 3). main energy in the 280-370 Hz band, the "bioduck" sounds from Antarctic minke whale in the 120-250 Hz band, and Antarctic blue whale "D-calls" (Erbe et al., 2017) in the 60-100 Hz band are recorded in TNB in December 2018. Supplementary Figure 2E shows a "chorus" of Antarctic blue whales (Thomisch et al., 2016) at around 27 Hz and 18 Hz in the TNB records.
Supplementary Figure 2F shows that anthropogenic ship noise continued for several days (April 29-13 May, 2017) and noisy status of the ocean continued for several hours during storms in TNB. Strong winds can enhance noise levels through various physical processes, such as higher ocean waves, fracturing and collision of ice blocks, and vibration and fracturing of icebergs and ice shelves (Chaput et al., 2018;Dziak et al., 2019). Humpback whale song in 20-100 Hz band (Payne and McVay, 1971; Supplementary Figure 2G) and long-duration (200 s) downsweeps from 75 to 62 Hz (Supplementary Figure 2H) were found in TNB in December 2016. Figure 3 shows the SPDs for the daily averaged spectrum of 1-year deployment records from BAL15 and TNBN17. To determine the seasonal variation in the dominant sound sources for each region, we plotted the SPDs for the entire period (first row; Figures 3A,B), austral winter months, June-August (second row; Figures 3C,D), early summer months, November-December (third row; Figures 3E,F), and late summer, January to April (fourth row; Figures 3G,H). The year-round plot in the BI region ( Figure 3A) shows a bimodal pattern with two separated red lines, and the lower lineation is largely the same as the SPD pattern in the austral winter season in the BI region ( Figure 3C). This implies that the bimodal pattern in the BI region is a seasonal difference. However, the bimodal pattern in the year-round SPD plot for TNB ( Figure 3B) is not as clear as that for the BI region, and one of the modes is not shown in any particular seasonal plot. In the SPD of the BI for the austral winter season (Figure 3C), sound levels in the 100-200 Hz band increased due to Antarctic minke whale calls. Antarctic minke whales can breathe in dense sea-ice environments by navigating between ice floes and breaking the ice to make breathing holes, and the BI region in austral winter season would be a safe refuge from killer whales (Orcinus orca) (Filun et al., 2020). In contrast to the BI region, there was no significant increase in the Antarctic minke whale bands, but leopard seal vocalizations contributed to the frequency band of 280-360 Hz in TNBP even in winter season. The leopard seal calls in the BI region and TNBP were strongest in the early summer season (Figures 3E,F). We plotted the SPDs in the late summer season in the frequency band of 10-100 Hz to identify fin and Antarctic blue whale vocalizations.
In the BI region (Figure 3G), there are strong peaks at 18 Hz, 27 Hz, and gentle peaks at 80 and 90 Hz. In the BI records, we could not present a spectrogram with clear fin or Antarctic blue whale calls due to the noise in the bands (Supplementary  Figures 2B, 4). However, the strong peaks at 18 and 27 Hz can be attributed by Antarctic blue whales, as reported in the previous acoustic monitoring in the Ross sea in 2004 (Širović et al., 2009). Meanwhile, the two peaks at 80 and 90 Hz, which are the same frequency bands of the signals shown in Supplementary  Figure 2B, reveal the existence of fin whales in the region. In the TNB late summer SPD (Figure 3H), the amplitude rise is not significant in the overtone fin whale call band (∼90 Hz), and the peaks at 18 Hz, 27 Hz associated with the Antarctic blue whale chorus (Supplementary Figure 2E) are visible. Another research has identified humpback whale signals off the Antarctic coast at about 70 • S (Van Opzeeland et al., 2013), and it is possible that some of the acoustic energy seen in the SPDs of the BI in the 50-300 Hz range is sourced from humpbacks. However, we did not see a definitively identifiable signal in the SPDs for humpback whales; this may be because humpback whale acoustic energy is spread over the spectrum from a few tens of hertz up to several kilohertz rather than being concentrated in a narrow band.

Spectral Amplitude and Its Correlation With Sea Ice and Wind Speed
To address the seasonal variations in the three dominant sound sources (fin and blue whale calls, cryogenic signals, and leopard seal calls) in the BI region and TNBP, mean sound levels were calculated within specific frequency bands: (1) 15-30 Hz for fin and blue whales, (2) 50-65 Hz for cryogenic signals, and (3) 300-350 Hz for leopard seals. Although we selected the 50-65 Hz band for the cryogenic signals where MMC were insignificant in the SPDs (Figure 3), the broadband energy from cryogenic signals can strongly affect other marine mammal bands. To address this issue, we estimate spectral levels by calculating MMC to the hourly averaged sound spectra (PSD MMC , Menze et al., 2017). Despite the fact that a lower Nyquist frequency (500 Hz) and intrinsic background fluctuations in the higher frequency range (BI: > 370 Hz, TNB: > 400 Hz) cause disturbances when estimating confident background noise levels in the core frequency band of leopard seals (310-360 Hz, Table 2), we practically select upper fitting bands to compare relative time variations in the PSD MMC for each mooring ( Table 2 and Supplementary Figures 3B,D). Figure 4 shows (A) the mean amplitudes of the BI recordings in the specific bands, (B) estimated PSD MMC value, and (C) the geophysical components, sea ice concentration, and wind speed. The mean amplitudes in the 15-30 Hz range (baleen whale band; red line in Figure 4A) and PSD MMC value associated with Antarctic blue whales (red line in Figure 4B) increase from February to April and turn downward from mid-April; these changes coincide with the rapid increase in sea-ice concentration (navy solid line in Figure 4C). The highest levels in April may be related to high densities of krill and greater abundance of baleen whales at the ice margins reported in previous studies  Figure 1D) represented as navy lines, and time series of daily wind speed from Jang Bogo station plotted as green lines.
(de la Mare, 1997;Brierley et al., 2002). The PSD MMC value of Antarctic blue whales almost disappeared in July 2015; then they reappeared in February 2016 along with the blue whales' acoustic presence in other Antarctic regions (Širović et al., 2009). The sea ice concentration in the BI region generally remains over 80% from May to November, and the sound levels in most frequency bands decreased during the sea-ice full period. There is an extraordinary event in which the sea ice concentration drops to 50% after a strong storm in late September (Figure 4C), and the sound levels increased in all the frequency bands ( Figure 4A). This implies that sea ice coverage is the major factor determining background noise levels in the BI region. The PSD MMC value of Antarctic minke whales (olive green line in Figure 4B) is generally significant from late May to mid-September during the sea-ice full period; however, the PSD MMC value of leopard seals is the highest from early November to January, which correspond to the breeding season of leopard seals (Southwell et al., 2003). Figure 5 shows plots similar to those in Figure 4, but for the TNBP during 2016-2018. The mean sound levels were estimated by combining the northern and southern records in 2017 and 2018. Although the narrow amplitude gaps between the bands in 2016 compared with the later period (2017-2018, Figure 5A) implies a poor frequency response correction, we can compare the relative amplitude variation within every deployed period. In 2016, noisy summer (January-March) and silent fall-austral winter-early spring season (April-October) were conspicuous features. Dziak et al. (2019) investigated hundreds of icequakes from the NIS and Campbell glaciers prior to the NIS calving event (April 7, 2016) from January to March 2016. The broadband energy from the icequakes peaked in all the frequency bands in February 2016 ( Figure 5A). The seasonal pattern in TNB 2017 was different from that in 2016. The amplitudes in the 50-65 Hz band (blue line in Figure 5A) were high in February but then decreased until mid-April. However, the sound levels rebounded in May and reached high levels in July. The reason for the "noisy winter" in 2017 can be found in the wind speed and sea-ice concentration. In TNBP, the sea ice concentrations frequently dropped below 80% even in the austral winter seasons (Figure 5C), in contrast to that in the BI region ( Figure 4C). This phenomenon is due to the blowing out of sea ice from the polynya by katabatic winds. Yoon et al. (2020) reported that the number of katabatic wind events (the wind direction is westerly, and the wind speed exceeds 25 ms −1 ) in TNBP was the highest in the 2017 austral winter (April-October) during the 2015-2017 period, and the polynya opened for the longest time period in 2017. Figure 5C shows that the sea ice concentration dropped most frequently in 2017. The frequent katabatic winds and high mobility of sea ice in open sea can make for a "noisy winter" in the polynya region in contrast to the silent winters in the BI region and other sea-ice regions, e.g., Beaufort Sea (Lewis and Denner, 1988;Kinda et al., 2013), Antarctic Peninsula (Dziak et al., 2010(Dziak et al., , 2015, and Atlantic sector of the Southern Ocean (Menze et al., 2017). In the TNBP, we estimate the PSD MMC only for Antarctic blue whales and leopard seals ( Figure 5B). Although we detected Antarctic minke whale signals in the TNBP (Supplementary Figure 2D), the PSD MMC estimation could not proceed because upper fitting windows are influenced by vocalizations of leopard seals from late August ( Figure 5B). The leopard seals' calls in the TNBP start earlier than those in the BI region and Atlantic sector of the Southern Ocean (Klinck, 2008;Van Opzeeland et al., 2010;Menze et al., 2017), and the chorus levels are remarkably loud. These results might be attributed to a higher density of leopard seals and year-round existence of pack ice in the polynya region, the inhabitation of leopard seals. Annual PSD MMC values of Antarctic blue whales in TNB are conspicuous in March; however, the amplitude variations do not show any correlations with sea-ice concentrations and vary yearly (Figure 5).
We normalize the hourly PSD MMC values between 0 and 1 for each day and plot them over a day to find diel patterns in the BI (Figure 6) and TNB (Figure 7) data. The PSD MMC value of Antarctic blue whale is slightly lower at midnight than at midday in the BI region and TNBP (Figures 6A, 7A), corresponding to the diel pattern of the calling rate in Southern Indian Ocean (Leroy et al., 2016). The PSD MMC value of Antarctic minke whales in the BI region is increased at midnight (Figure 6B). Minke whales were observed to be more active in nighttime in Southern Ocean (Menze et al., 2017;Shabangu et al., 2020) and Northern Ocean (Risch et al., 2013(Risch et al., , 2019, owing to the diurnal vertical migration of their prey (zooplankton) into deepwater during the day and to the surface at nighttime (Demer and Hewitt, 1995). However, the seasonal variation in the diel patterns was not significant in our study (Supplementary Figure 5). Figures 6C, 7B show the active calling of leopard seals during nighttime, which reduces in daytime. Similar diel patterns have been reported by previous studies (Thomas and DeMaster, 1982;Klinck, 2008;Van Opzeeland et al., 2010), and these patterns most likely arise from the haul-out behavior of the leopard seals that rest on ice during the solar midday.
To investigate the seasonal changes in the sea-ice concentration and their effects on the marine mammal vocalization, we plotted the daily averaged PSD MMC values of Antarctic blue whales, Antarctic minke whales (only for the BI data), and leopard seals over the day of the year with sea ice concentration represented by color (Figure 8). The bimodal pattern of sea-ice concentrations in the BI region is clearly represented by a sudden color change from navy to yellow in  April, while the PSD MMC value of Antarctic blue whales starts decreasing with a sudden change in sea-ice concentrations in the BI region ( Figure 8A). Supplementary Figure 6 shows a scatterplot of the PSD MMC values of Antarctic blue whales, together with the number of days that the sea ice concentration exceeded 80% in the last 100 days; the distance to the closest open sea has also been represented in different colors. Furthermore, whale sound levels diminish with an increase in the duration of the ice-covered period and the distance to the nearest open sea. This may be due to the migration of whales to find reliable access to open water for breathing (Van Opzeeland et al., 2013;Thomisch et al., 2016). The presence period of Antarctic blue whale calls in the TNBP is shorter than that in the BI region, and we found no persistent correlation of this result with the sea-ice concentrations in the TNBP (Figure 8E). The Antarctic minke whale calls are observed in the BI region from May to October when the sea ice concentration is higher than 70% (Figure 8B). The leopard seal PSD MMC increasing from early November in the BI region and from late August in the TNBP, which resulted in an energy peak in late December that ceased to exist in February (Figures 8C,E).
The austral winter season (July-August) in the TNBP is another noisy period attributed to repeated storm noise and cryogenic signals as shown in the long-term spectrogram (Figures 2B,C, 5A). To determine the dominant environmental factor controlling low-frequency cryogenic noise, we created scatterplots of sound levels in the 50-65 Hz band over the wind speed with sea-ice concentration represented by color (Figure 9). The linear fittings to the data of sea-ice concentration greater than 80% and less than 80% are represented by orange and blue lines, respectively. Table 3 lists the slopes, intercepts of the fitting lines, and Pearson's correlation coefficients between the sound levels and the wind speed. In the BI region, the significant difference in the intercepts under the two sea-ice conditions (89 dB for SIC < 80%; 80 dB for SIC > 80%), gentle slopes, and low correlation coefficients indicate that the low-frequency cryogenic noise is dominantly controlled by sea-ice conditions in the BI region. In contrast, sea-ice concentration did not cause significant changes in TNB. The fitting lines for the two sea ice concentrations in TNB were largely the same, and the slopes and correlation coefficients were greater than those of the BI record. This implies that the vibration and fracturing of ice shelves related to strong winds could be a dominant source generating low-frequency noise in the TNBP rather than sea ice.

CONCLUSION
We reviewed soundscapes obtained from long-term hydroacoustic recording in two different regions of the Western Ross Sea, Antarctic: the BI region and TNBP.
The study showed that cryogenic events and marine mammals are the major sound sources in both regions. Tectonic earthquake and anthropogenic noises due to vessels were also found at various times in the spectrograms. The major marine mammal signals in the long-term spectrogram and SPDs of both regions are from Antarctic blue whales in the late summer-fall months and leopard seals in early summer. Typical fin whale signals are not detected in both region despite there are peaks at 80 and 90 Hz in the SPD plots of the BI region. A fin whale chorus was observed between March and May in the BI recordings; however, it could not be identified in the TNB recordings, which is in accordance to the study conducted by Širović et al. (2009) in the Ross Sea. Antarctic minke whales calls are detected in both regions and peaked in austral winter season in the BI region reflecting their preference to dense sea-ice environments (Filun et al., 2020). We found humpback whale vocalization and longduration (>200 s) downsweeps of unknown origin from 75 to 62 Hz in TNBP. Correlation with sea ice in the BI region reflects baleen whale behaviors such as high density in the ice margin and migrations to find reliable access to open water for breathing. The diel patterns associated with the vocalization of Antarctic blue whales (louder in daytime) and Antarctic minke whales (louder in nighttime) may to be related to their feeding ecology. Louder vocalization of leopard seals in midnight results from their haul-out behavior.
The influences of wind and sea ice on ambient noise were different in the two regions. In the BI region, sea ice, often near 100% concentration, reduced ambient sound levels, similarly to other polar regions. However, the ambient noise in TNB increased due to strong winds, regardless of the seaice concentration; hence, the austral winter in 2017 was very noisy because of the frequent katabatic winds coming from the land. The regional differences seem to be attributed to the higher mobility of sea ice by the polynya process and the strong noise generated near ice shelves.
Although monitoring long-term soundscapes in extremely cold waters is challenging because of the logistical difficulties and harsh conditions at sea, it can significantly help understand the spatiotemporal changes in the Southern Ocean in the context of cryosphere changes, and can provide a means of assessing the status and trends of biodiversity in the Ross Sea region Marine Protected Area. Future intensive studies, which quantitatively measure the occurrence of individual species and cryogenic events as well as elucidate acoustic propagation by considering sea-ice and water temperature variations, will provide a better understanding of the soundscapes in Southern Ocean.

ACKNOWLEDGMENTS
We thank the captain and crew of R/V Araon, J. Lee, M. Fowler, and S. Nieukirk for the deployment and recovery of hydrophone arrays. We also thank T. Choi, ECMWF, University of Bremen, NSIDC for sharing their meteorological and sea ice data. We also thank the reviewers CW and SM as well as the topical editor for the review, which greatly improved this manuscript.  (Porter, 1992) at 25 and 400 Hz. The sound speed profiles are calculated from the conductivity, temperature, and depth (CTD) measurements performed by the icebreaker R/V Araon (Dziak et al., 2019). Blue solid lines represent the transmission loss in a 400 m deep receiver when the sound source is at a depth of 50 m in the case of a flat sea surface with a flat bottom.

SUPPLEMENTARY MATERIAL
Supplementary Figure 2 | Spectrograms of the various acoustic sources recorded in the BI region or TNBP as shown above each spectrogram: (A) earthquake (T-phase) and impulsive icequakes; (B) baleen whale calls and noise in the BI region; (C) bioduck sounds from Antarctic minke whale in the BI region; (D) leopard seal sound in the 200-400 Hz band, bioduck sounds of Antarctic minke whale in the 120-250 Hz, and Antarctic blue whale's "D-calls" in 50-100 Hz; (E) a "chorus" of Antarctic blue whales (Thomisch et al., 2016) at around 27 Hz and 18 Hz; (F) broadband storm noise continued for several hours and anthropogenic noise from a ship that continued for several days; (G) leopard seals and humpback whale vocalizations; (H) undefined downsweeps continued hundreds of seconds.
Supplementary Figure 3 | Examples of ambient sound spectrum with marine mammal contributions (MMC) from Antarctic blue whales, Antarctic minke whales, and leopard seals for the BI (BAL15) and TNB (TNBN17) regions. Measured power spectral density (PSD) (black line) and fitted interpolation functions in the frequency bands associated with Antarctic blue whales (red), Antarctic minke whales (olive green), and leopard seals (green). PSD MMC values are estimated in the shaded frequency bands by using the method adopted by Menze et al. (2017).
(A) Averaged PSD MMC recorded in BAL15 on 2015.06.20 from 12:00 to 13:00 and the fitted interpolation functions for bands associated with Antarctic blue whales and minke whales. Note that the PSD MMC values of Antarctic minke whales is estimated only for the BAL15 data from February to September, where the upper fitting band (300-351 Hz) is stable due to the absence of sounds emitted by leopard seals. (B) Averaged PSD MMC recorded in BAL15 on 2015.12.10 from 12:00 to 13:00 and the fitted interpolation functions for the band associated with leopard seals. PSD MMC values of leopard seals in the BI region are estimated from October to January, where the lower fitting band (39-161 Hz) is stable in the absence of Antarctic minke whale calls. (C) Averaged PSD MMC recorded in TNBN17 on 2017.03.15 from 12:00 to 13:00 and fitted interpolation functions for Antarctic blue whales. Owing to the big fluctuation under 18 Hz and the superposed sound energy from various sources in the neighboring frequency band, a narrower fitting band than that applied to the BI data is adopted. (D) Averaged PSD MMC recorded in TNBN17 on 2017.12.10 from 20:00 to 21:00 (solid black line) and on 2017.08.20 from 12:00 to 13:00 (dashed black line), and the fitted interpolation functions for leopard seals. We could not calculate the MMC of Antarctic minke whales in the TNBP, as leopard seals produced sounds from late August.