Climate change is rapidly transforming the coastal margins of Greenland. At the same time, there is increasing recognition that marine-terminating glaciers provide unique and critical habitats to ice-associated top predators. We investigated the connection between a top predator occupying glacial fjord systems in Northwest Greenland and the properties of Atlantic-origin water and marine-terminating glaciers through a multiyear interdisciplinary project. Using passive acoustic monitoring, we quantified the summer presence and autumn departure of narwhals (Monodon monoceros) at glacier fronts in Melville Bay and modeled what glacier fjord physical attributes are associated with narwhal occurrence. We found that narwhals are present at glacier fronts after Greenland Ice Sheet peak summer runoff and they remain there during the period when the water column is becoming colder and fresher. Narwhals occupied glacier fronts when ocean temperatures ranged from –0.6 to 0.8°C and salinities between 33.2 and 34.0 psu at around 200 m depth and they departed on their southbound migration between October and November. Narwhals' departure was approximately 4 weeks later in 2019 than in 2018, after an extreme 2019 summer heatwave event that also delayed sea ice formation by 2 months. Our study provides further support for the niche conservative narwhal's preference for cold ocean temperatures. These results may inform projections about how future changes will impact narwhal subpopulations, especially those occupying Greenland glacial fjords.

Greenland's marine-terminating glaciers connect the ice sheet to the ocean and provide a critical boundary where heat, freshwater, and nutrient exchanges take place. Buoyant freshwater runoff from inland ice sheet melt is discharged at the base of marine-terminating glaciers, forming vigorous upwelling plumes. It is understood that subglacial plumes modify waters near glacier fronts and increase submarine glacier melt by entraining warm ambient waters at depth. However, ocean observations along Greenland's coastal margins remain biased toward summer months which limits accurate estimation of ocean forcing on glacier retreat and acceleration. Here, we fill a key observational gap in northwest Greenland by describing seasonal hydrographic variation at glacier fronts in Melville Bay using in situ observations from moorings deployed year-round, CTDs, and profiling floats. We evaluated local and remote forcing using remote sensing and reanalysis data products alongside a high-resolution ocean model. Analysis of the year-round hydrographic data revealed consistent above-sill seasonality in temperature and salinity. The warmest, saltiest waters occurred in spring (April–May) and primed glaciers for enhanced submarine melt in summer when meltwater plumes entrain deep waters. Waters were coldest and freshest in early winter (November–December) after summer melt from sea ice, glacier ice, and icebergs provided cold freshwater along the shelf. Ocean variability was greatest in the summer and fall, coincident with increased freshwater runoff and large wind events before winter sea ice formation. Results increase our mechanistic understanding of Greenland ice-ocean interactions and enable improvements in ocean model parameterization.

Passive acoustic monitoring has been an effective tool to study cetaceans in remote regions of the Arctic. Here, we advance methods to acoustically identify the only two Arctic toothed whales, the beluga (Delphinapterus leucas) and narwhal (Monodon monoceros), using echolocation clicks. Long-term acoustic recordings collected from moorings in Northwest Greenland were analyzed. Beluga and narwhal echolocation signals were distinguishable using spectrograms where beluga clicks had most energy >30 kHz and narwhal clicks had a sharp lower frequency limit near 20 kHz. Changes in one-third octave levels (TOL) between two pairs of one-third octave bands were compared from over one million click spectra. Narwhal clicks had a steep increase between the 16 and 25 kHz TOL bands that was absent in beluga click spectra. Conversely, beluga clicks had a steep increase between the 25 and 40 kHz TOL bands that was absent in narwhal click spectra. Random Forest classification models built using the 16 to 25 kHz and 25 to 40 kHz TOL ratios accurately predicted the species identity of 100% of acoustic events. Our findings support the use of echolocation TOL ratios in future automated click classifiers for acoustic monitoring of Arctic toothed whales and potentially for other odontocete species.