Ocean Acoustic Tomography and Thermometry Statistics (2026): 60+ Data Points on Heard Island and ATOC, Seismic Ocean Thermometry, and Argo Heat Content

Ocean acoustic thermometry statistics 2026: earthquake sound gauged a 3,000 km Indian Ocean path with 0.0060 K error, plus Heard Island, ATOC and Argo data.

A sound wave crossing 3,000 km of the East Indian Ocean was enough to pin down the average temperature of that whole section with a standard error of just 0.0060 kelvin (Wu et al., Science 2020). The sound came from small repeating earthquakes, not from a ship or a loudspeaker, and the result revived an idea first tested in 1991, when coded signals from Heard Island were monitored at 16 sites around the globe (Munk et al., JASA 1994). The stakes are large: ocean warming accounted for 91% of the heating in the climate system (IPCC AR6 WG1 Summary for Policymakers 2021), and the upper 2,000 m gained about 23 ± 8 ZJ in 2025 alone (Pan et al., Advances in Atmospheric Sciences 2026). We aggregated data from Science, the Journal of the Acoustical Society of America (JASA), Journal of Geophysical Research, the Proceedings of the Royal Society A, Caltech, the Argo Program, the IPCC, NOAA NCEI, the CTBTO and the University of Washington Heard Island archive. This roundup is about temperature; for shipping noise, sonar and marine mammals, see our separate underwater acoustics and ocean sonar statistics.

TL;DR

  • The 1991 Heard Island Feasibility Test was monitored at 16 sites across six ocean regions (Munk et al., JASA 1994).
  • Heard Island transmitted at source levels of 209-220 dB re 1 µPa from 175 m depth (Bowles et al., JASA 1994).
  • Only about half of North Pacific seasonal and year-to-year sea level change came from thermal expansion (ATOC Consortium, Science 1998).
  • A 1994 test at 3,250 km range estimated a six-day temperature change with an uncertainty of about 0.006°C (Worcester et al., JASA 1999).
  • Seismic ocean thermometry resolved a 3,000 km section with a standard error of 0.0060 K over 2005-2016 (Wu et al., Science 2020).
  • Warming in that region appeared to be as much as 69 percent greater than had been believed, depending on the comparison dataset (Caltech News 2020).
  • The Kuroshio Extension warmed at 11.8 ± 5.0 mK per year (2σ) from 2008 to 2021 by seismic and in situ data (Peng and Callies, Proceedings of the Royal Society A 2025).
  • Argo runs close to 4,000 active floats and collected its 3 millionth profile in July 2024 (Argo Program status page).
  • Ocean warming accounted for 91% of the heating in the climate system (IPCC AR6 WG1 SPM 2021).
  • The upper 2,000 m ocean gained about 23 ± 8 ZJ in 2025 relative to 2024 (Pan et al., Advances in Atmospheric Sciences 2026).
  • The CTBTO hydroacoustic network has 11 stations: six hydrophone stations and five T-phase stations (CTBTO).

1. The Physics: Why Travel Time Works as a Thermometer

Acoustic thermometry rests on one stubborn fact: a 1°C change in water temperature changes sound speed by about 4 m/s (University of Washington Applied Physics Laboratory, Heard Island Feasibility Test archive). Warm water carries sound faster, so a path that warms delivers its signal a little earlier. The trick is that the clock does the averaging. A single thermometer reads one point, while a pulse crossing thousands of kilometers reports the mean of everything it passed through.

That averaging is the whole point, because the signal hunted for is tiny and the local noise is huge. The Heard Island archive puts the nominal climate signal at about 0.005°C per year at 1,000 m depth, while temperatures at a single point can swing by 1-2°C from one month to the next. On those figures, one month of local “ocean weather” (eddies and internal waves, the ocean’s equivalent of passing storms) is 200 to 400 times larger than a full year of the climate trend (1 / 0.005 = 200; 2 / 0.005 = 400). Point sensors have to average that noise away statistically; a long acoustic path averages it physically.

MetricValueSource
Sound speed change for a 1°C change in water temperatureAbout 4 m/sUniversity of Washington APL, Heard Island archive
Nominal climate warming signal at 1,000 m depthAbout 0.005°C per yearUniversity of Washington APL, Heard Island archive
Depth of the sound speed minimum (SOFAR channel axis) in the deep mid-latitude oceanAbout 800 to 1,000 mDOSITS (University of Rhode Island), Sound Travel in the SOFAR Channel 2023
Launch angle below which rays stay trapped in the deep sound channelLess than about 12° from horizontalDOSITS (University of Rhode Island) 2023
Travel-time perturbation caused by mesoscale eddiesRoughly 100 msHowe et al., Frontiers in Marine Science 2019
Travel-time change caused by a Gulf Stream or Kuroshio meanderAbout 1 secondHowe et al., Frontiers in Marine Science 2019
Receiver timing and positioning accuracy needed for tomography1 ms and 1 mHowe et al., Frontiers in Marine Science 2019

Context: the SOFAR channel (Sound Fixing and Ranging) is a depth band where sound speed is lowest, so sound bends back toward it from above and below instead of escaping. DOSITS places its axis at about 800 to 1,000 m in the deep mid-latitude ocean. Most recent available data for the timing and eddy figures: Howe et al., Frontiers in Marine Science 2019.

2. The Heard Island Feasibility Test (1991)

Heard Island was a proof of concept, and the failure statistics are as telling as the successes. The test made 35 transmissions, and its 10 acoustic sources failed with a half-life of 1.2 days (University of Washington APL, Heard Island archive, The Experiment). By the time 25 transmissions had been made, only three sources were still working, and the last transmission used the final two. The question it answered was not “how warm is the ocean” but “can a coded signal survive a trip around the planet and still be timed precisely”.

It could. The Munk et al. (1994) overview in JASA reports that the transmissions were monitored at 16 sites spread across the North Atlantic, South Atlantic, North Pacific, South Pacific, Indian Ocean and Southern Ocean. A separate marine mammal survey found no drop in the density of whale and dolphin groups near the source during transmissions, though sperm whales, heard during 23% of baseline listening time, went unheard during transmissions (Bowles et al., JASA 1994). ATOC later ran dedicated marine mammal research alongside its thermometry (see section 3).

MetricValueSource
Start date, transmissions made, sources deployed and source failure half-lifeJanuary 26, 1991 (planned for 10 days); 35 transmissions; 10 sources; half-life of 1.2 daysUniversity of Washington APL, Heard Island archive
Number of sites that monitored the transmissions16Munk et al., JASA 1994
Source center frequency, maximum bandwidth and source depth57 Hz, 30 Hz, 175 mBowles et al., JASA 1994
Source level of the transmissions209-220 dB re 1 µPaBowles et al., JASA 1994
Cetacean group density before vs during transmissions0.0157 vs 0.0166 groups per km²Bowles et al., JASA 1994
Sperm whales heard: share of 1,181 baseline minutes vs 1,939 transmission minutes23% vs none heardBowles et al., JASA 1994

Outlier note: the idea was older than the test. A 1960 precursor detonated explosive charges off Perth and listened near Bermuda, about 19,820 km away (University of Washington APL, Heard Island archive). Transmissions ran one hour in every three (Bowles et al., JASA 1994). Most recent available data: 1991 experiment, published 1994.

3. ATOC and a Decade of North Pacific Thermometry (1994-2006)

ATOC turned the Heard Island demonstration into an observing system, and its key number is a precision figure. Acoustic travel times were measured with a precision of about 20-30 milliseconds at 3,000-5,000 km ranges, equivalent to about 0.006°C averaged over range and depth; at 5,000 km the signal takes nearly an hour to arrive (University of Washington APL, ATOC epilogue). That is timing to roughly one part in 100,000 of the travel time (0.03 s / 3,600 s is about 0.0008%).

The science payoff was a correction to how sea level was being read. Comparing gyre-scale heat content from tomography with satellite altimetry and a circulation model, the ATOC Consortium (Science 1998) found that only about half of seasonal and year-to-year sea level change came from thermal expansion. Over 1996-2006 the program fed four independent estimates of the North Pacific temperature field (Dushaw et al., JGR Oceans 2009). It also came with dedicated marine life studies: juvenile elephant seals exposed to the source showed only subtle diving changes, judged likely to have minimal biological significance (Costa et al., JASA 2003). Chronic noise from shipping is a separate exposure question, covered in our ship propeller cavitation noise statistics.

MetricValueSource
ATOC partners at launch11 institutions from 7 nations, 2 sources, 18 receiversMIT News 1994
ATOC source frequency75 HzMIT News 1994
1994 test range and receiver3,250 km, 700 m vertical arrayWorcester et al., JASA 1999
Uncertainty of a six-day temperature change at 3,250 kmAbout 0.006°CWorcester et al., JASA 1999
Uncertainty of range-averaged sound speed at 3,250 kmAbout 0.05 m/sWorcester et al., JASA 1999
Annual cycle of heat flux (peak-to-peak) in the gyre-scale North Pacific150 ± 25 W/m², a 0.2°C cycle in vertically averaged temperatureATOC Consortium, Science 1998
Elephant seal study: ATOC source depth and source level939 m, 195 dB re 1 µPa at 1 mCosta et al., JASA 2003
Seals whose descent rate shifted more than one standard deviation from control during exposure9 of 14Costa et al., JASA 2003

Context: received sound levels at the seals averaged 128 dB, ranging from 118 to 137 dB re 1 µPa in the 60-90 Hz band, against ambient levels of 87-107 dB (Costa et al., JASA 2003). Most recent available data: ATOC ended in 2006; see Dushaw et al., JGR Oceans 2009 for the full decade.

4. Seismic Ocean Thermometry: The 2020 Breakthrough

Dedicated sound sources need ships, moorings and environmental permits. The 2020 idea was to stop building sources and use the ones the planet already provides. Wu et al. used 2,047 pairs of repeating earthquakes off Sumatra to track temperature across a 3,000 km section of the equatorial East Indian Ocean from 2005 to 2016, with a standard error of 0.0060 kelvin (Wu et al., Science 2020; pair count as reported by Smithsonian Magazine 2020). Repeating earthquakes are quakes that rupture the same patch of fault again and again, so each pair gives two near-identical sound sources years apart.

The headline that travelled was a warming trend larger than earlier estimates: Caltech reported that ocean warming in the region appeared to be as much as 69 percent greater than had been believed, depending on which dataset is used. The same release carried the caveat that more data were needed before drawing conclusions. Note the T wave here: it is the part of an earthquake’s energy that converts to sound and travels through the water, not a wave in the solid rock.

MetricValueSource
Length of ocean section averaged and period analyzed3,000 km (equatorial East Indian Ocean), 2005-2016Wu et al., Science 2020
Standard error of the section-mean temperature anomaly0.0060 kelvinWu et al., Science 2020
Time scales of fluctuations detected12 months, 6 months and about 10 daysWu et al., Science 2020
Repeating earthquake pairs used2,047Smithsonian Magazine 2020, reporting the Science paper
Sound travel time vs the part of it caused by water temperatureAbout half an hour vs about one-tenth of a secondCaltech News 2020
Warming in the region vs earlier estimatesAs much as 69 percent greater, depending on the comparison datasetCaltech News 2020
Depths sampled by the earthquake sound wavesBelow 2,000 m, where conventional measurements are fewCaltech News 2020

Context: by October 2026 the paper had 78 citing works in OpenAlex. The 69 percent figure is a regional, dataset-dependent comparison, not a global revision of ocean warming (see Data watch in Methodology).

The follow-up papers turned a single demonstration into a toolkit, and the most consequential number is a trend with an honest error bar. Seismic and in situ data together show steady subsurface warming of 11.8 ± 5.0 mK per year (2σ uncertainty) in the Kuroshio Extension from 2008 to 2021, with no substantial trend between 1997 and 2008 (Peng and Callies, Proceedings of the Royal Society A 2025). Converted to the units climate reports use, 11.8 mK per year is 0.118 K per decade (11.8 × 10 / 1,000).

The same paper shows why seismic data is a complement rather than a rival to floats. Relative to the prior, eddy-driven error variance in area-mean temperature fell by 30% with in situ data, 39% with seismic data and 50% with both; for path-mean estimates the combined reduction was 83% versus 45% from in situ data alone. Meanwhile, sensors keep getting more sensitive: CTBTO hydrophones allow smaller earthquakes (magnitude below 4.0) to be used (Wu et al., JGR Solid Earth 2023), and a fiber-optic cable on the seafloor detected 92 oceanic earthquakes in four days, most missing from existing catalogs (Shen and Wu, JGR Solid Earth 2024).

MetricValueSource
Kuroshio Extension subsurface warming rate, 2008-202111.8 ± 5.0 mK per year (2σ)Peng and Callies, Proceedings of the Royal Society A 2025
Eddy error-variance reduction in area-mean temperature: in situ only vs seismic only vs both30% vs 39% vs 50%Peng and Callies, Proceedings of the Royal Society A 2025
Path-mean error-variance reduction: seismic plus in situ vs in situ alone83% vs 45%Peng and Callies, Proceedings of the Royal Society A 2025
Start of the Wake Island station’s digital record used for the Kuroshio work1997Peng et al., JGR Oceans 2024
Earthquake magnitude usable with CTBTO hydrophonesBelow 4.0Wu et al., JGR Solid Earth 2023
Oceanic earthquakes detected by an ocean-bottom fiber-optic cable in a 4-day experiment off Oregon92Shen and Wu, JGR Solid Earth 2024
Smallest repeating earthquakes usable with the fiber-optic cableM 3.5Shen and Wu, JGR Solid Earth 2024
Frequency bands used in a published step-by-step workflow for passive deep-ocean thermometry3-6 Hz and 6-12 HzSmets et al., JGR Oceans 2022

Context: distributed acoustic sensing (DAS) turns an ordinary seafloor fiber-optic cable into thousands of vibration sensors along its length. The method is also spilling into ocean dynamics: seismic tomography in the East Indian Ocean picked up biweekly temperature signals carried by Yanai waves (a type of wave trapped along the equator), and the variance in that band differed by up to a factor of 2 between ocean models (Peng et al., Journal of Physical Oceanography 2025).

6. Argo and Ocean Heat Content: The Yardstick Acoustic Methods Are Checked Against

Every acoustic result is ultimately compared with Argo, so its coverage defines what acoustics can add. Argo runs close to 4,000 active floats and collected its 3 millionth profile in July 2024 (Argo Program status), at about 13,000 profiles a month from 30 countries (Argo Program, About). But a standard float descends to 2 km and no further. The deep ocean below that is exactly where earthquake sound samples (see section 4), and where the Deep Argo mission is working toward a standing array of 1,228 floats (Argo Program, Deep Argo Mission).

The heat being measured is enormous and accelerating. Ocean warming accounted for 91% of the heating in the climate system, and the observed heating rate rose from 0.50 [0.32 to 0.69] W m⁻² over 1971-2006 to 0.79 [0.52 to 1.06] W m⁻² over 2006-2018 (IPCC AR6 WG1 Summary for Policymakers 2021). In 2025 the upper 2,000 m gained about 23 ± 8 ZJ relative to 2024 by the IAP/CAS estimate, with about 33% of the ocean area in its top three warmest conditions on record (Pan et al., Advances in Atmospheric Sciences 2026).

MetricValueSource
Active Argo floatsClose to 4,000Argo Program status page
Argo profile milestones1 millionth in November 2012, 2 millionth in November 2018, 3 millionth in July 2024Argo Program status page
Deep Argo target array size and maximum float depth1,228 floats; some models capable of 6,000 mArgo Program, Deep Argo Mission
Share of climate-system heating that went into the ocean91%IPCC AR6 WG1 SPM 2021
Global ocean heat content increase, 1971-20180.396 [0.329 to 0.463] yottajoules (likely range)IPCC AR6 WG1 Chapter 9, 2021
Upper 2,000 m ocean heat gain in 2025 vs 2024 (IAP/CAS estimate)About 23 ± 8 ZJPan et al., Advances in Atmospheric Sciences 2026
Share of global ocean area in its top three (and top five) warmest conditions on record in 2025About 33% (about 57%)Pan et al., Advances in Atmospheric Sciences 2026

Context: NOAA NCEI publishes its 0-700 m and 0-2,000 m heat content series from 1955 to the present, updated quarterly (NOAA NCEI Global Ocean Heat Content). Global sea surface temperature in 2025 was 0.49°C above the 1981-2010 baseline (Pan et al. 2026).

7. Arctic Thermometry and the Global Listening Networks

The Arctic is where acoustic thermometry has delivered some of its clearest climate signals, in an ocean where sea ice gets in the way of most other instruments. In 1994, 19.6 Hz signals sent about 2,630 km across the Arctic showed the Atlantic Layer to be about 0.4°C warmer than expected, and a 1999 follow-up found about 0.5°C of further warming (Howe et al., Frontiers in Marine Science 2019). The Atlantic Layer is the warm, salty water from the Atlantic that flows into the Arctic beneath the colder surface layer.

The newest data come from CAATEX, a US-Norwegian campaign whose results were published in 2026. Six moorings transmitted and recorded every three days during 2019-2020, and signals carried over 2,600 km (AIP Scilights 2026). The paper concludes that long-range transmissions can give precise, year-round under-ice measurements, although large mode-2 travel-time uncertainty prevented a definitive comparison with the 1994 data (Dzieciuch et al., JASA 2026). Lower latitudes depend on infrastructure built for other reasons, chiefly the CTBTO hydroacoustic network, the same kind of always-on listening now used to score coral reef soundscape health.

MetricValueSource
1994 Transarctic (TAP) signal frequency and path length19.6 Hz over about 2,630 kmHowe et al., Frontiers in Marine Science 2019
Atlantic Layer temperature in 1994 vs expectation, and further warming by April 1999About 0.4°C warmer; about 0.5°C further warmingHowe et al., Frontiers in Marine Science 2019
CAATEX moorings, center frequency and bandwidth (2019-2020)6 moorings; 35 Hz; about 4 Hz full bandwidthDzieciuch et al., JASA 2026
CAATEX transmission interval and distance coveredEvery three days; over 2,600 kmAIP Scilights 2026
CTBTO hydroacoustic stations11 (6 hydrophone stations, 5 T-phase stations)CTBTO, Hydroacoustic Monitoring
Depth of CTBTO hydrophones and design lifetime600 to 1,200 m; 20 to 25 yearsCTBTO, Hydroacoustic Monitoring
North Atlantic tomography array2 sources and 7 receiversDushaw, Journal of Atmospheric and Oceanic Technology 2018
North Atlantic temperature uncertainty: one data type vs acoustic plus other dataAbout 50% of prior vs about 25% of prior (94% of variance resolved)Dushaw, Journal of Atmospheric and Oceanic Technology 2018

Context: a T-phase station is a seismometer on a steep oceanic island that picks up sound waves converting into ground motion at the shore. The CTBTO record has value for history as well as for monitoring: researchers point to a nearly 20 year record of treaty recordings as a basis for historical ocean thermometry (Dall’Osto and Heyburn, JASA 2019). For the North Atlantic array, see Dushaw, Journal of Atmospheric and Oceanic Technology 2018.

Summary: Ocean Acoustic Tomography and Thermometry by the Numbers

MetricValueSource
Sound speed change per 1°C of water temperatureAbout 4 m/sUniversity of Washington APL, Heard Island archive
Heard Island 1991: sites monitoring the transmissions16Munk et al., JASA 1994
Heard Island source level209-220 dB re 1 µPaBowles et al., JASA 1994
Six-day temperature change uncertainty at 3,250 kmAbout 0.006°CWorcester et al., JASA 1999
Share of North Pacific seasonal and year-to-year sea level change due to thermal expansionOnly about halfATOC Consortium, Science 1998
Seismic thermometry standard error over a 3,000 km section0.0060 kelvinWu et al., Science 2020
Repeating earthquake pairs in the first study2,047Smithsonian Magazine 2020, reporting the Science paper
Regional warming vs earlier estimatesAs much as 69 percent greaterCaltech News 2020
Kuroshio Extension warming, 2008-202111.8 ± 5.0 mK per year (2σ)Peng and Callies, Proceedings of the Royal Society A 2025
Path-mean error-variance reduction, seismic plus in situ vs in situ alone83% vs 45%Peng and Callies, Proceedings of the Royal Society A 2025
Oceanic earthquakes found by a seafloor fiber-optic cable in 4 days92Shen and Wu, JGR Solid Earth 2024
Active Argo floatsClose to 4,000Argo Program status page
Share of climate-system heating absorbed by the ocean91%IPCC AR6 WG1 SPM 2021
Upper 2,000 m ocean heat gain, 2025 vs 2024About 23 ± 8 ZJ (IAP/CAS)Pan et al., Advances in Atmospheric Sciences 2026
Arctic Atlantic Layer, 1994 acoustic measurement vs expectationAbout 0.4°C warmerHowe et al., Frontiers in Marine Science 2019
CTBTO hydroacoustic stations11CTBTO, Hydroacoustic Monitoring

Methodology and Sources

We collected about 75 raw data points from peer-reviewed papers, agency pages and institutional press releases, kept the 60-plus that could be traced to a primary document, and dropped figures we could only see in search snippets (including a widely repeated Heard Island reception distance we could not confirm in a primary page). Abstract-level figures for papers behind paywalls were read from their indexed abstracts.

Last updated: October 25, 2026. We update this roundup quarterly, and the next refresh is expected when the IAP/CAS-led annual ocean heat content update for 2026 is published in Advances in Atmospheric Sciences in January 2027.

Try VoxBooster — 3-day free trial.

Real-time voice cloning, soundboard, and effects — wherever you already talk.

  • No credit card
  • ~30ms latency
  • Discord · Teams · OBS
Try free for 3 days