Underwater Acoustic Communication Range Calculator

This underwater acoustic communication range calculator closes a one-way modem link budget before a sea trial. It combines source level, receiver noise, bandwidth, directivity, detection threshold, spreading geometry and seawater absorption. The result is the estimated marginal range, together with transmission-loss components, sound speed and one-way acoustic latency.

Introduction to underwater acoustic communication range

Underwater radio and optical links usually have limited practical reach, so sound remains the standard carrier for AUV telemetry, acoustic releases, subsea instruments and positioning systems. The central design question is how much propagation loss a modem can tolerate before its received signal-to-noise ratio falls below the demodulator’s requirement.

All acoustic levels entered here must use underwater conventions. Source level is expressed in dB re 1 µPa at 1 m, and received noise is expressed in dB re 1 µPa or dB re 1 µPa²/Hz. An airborne level referenced to 20 µPa is not interchangeable with an underwater value. Confirm both the pressure reference and the 1 m source-distance reference before using a datasheet figure.

How to use the underwater modem range estimator

Start with the modem’s source level and carrier frequency. Frequency strongly affects absorption: lowering it generally increases range, although it may also reduce available bandwidth and data rate. Directivity index represents receive-array gain against isotropic noise; use 0 dB for an omnidirectional hydrophone. Detection threshold is the in-band SNR needed for the target packet or bit error rate. A practical field design can include a fading reserve by adding it to this threshold.

Noise may be entered as an integrated band level or as a 1 Hz spectrum level. When spectrum mode is selected, the calculator adds the bandwidth correction automatically. Enter the receiver’s effective noise bandwidth, not merely a nominal channel spacing.

Choose spherical spreading for an unbounded wavefront, cylindrical spreading only when a waveguide is credible, or the intermediate exponent for a simplified shallow-water estimate. Temperature, salinity, depth and pH feed the Francois-Garrison model. A CTD cast or local climatology is preferable to a guess because environmental values can shift both absorption and sound speed.

Formulas for the sonar equation and seawater absorption

The one-way communication form of the sonar equation is

Here $\mathrm{SL}$ is source level, $\mathrm{TL}$ is transmission loss, $\mathrm{DI}$ is receiver directivity index and $\mathrm{NL}$ is in-band noise. At the marginal range, available SNR equals detection threshold $\mathrm{DT}$, giving

If noise is supplied as spectrum level $\mathrm{NSL}$, it is integrated over bandwidth $B$ in hertz:

The propagation model adds geometric spreading and frequency-dependent absorption:

The spreading exponent is $N=20$ for spherical spreading, $N=10$ for cylindrical spreading and $N=15$ for the practical intermediate law.

Range inside the logarithm must be measured from the 1 m source reference, while absorption is conventionally in dB/km. The implementation therefore evaluates

Francois-Garrison separates absorption into boric acid, magnesium sulphate and pure-water terms:

With $f$ in kHz, $T$ in °C, $S$ in parts per thousand, $D$ in metres and $\theta =T+273$ in kelvin, the boric acid coefficients are

The magnesium sulphate coefficients are

The pure-water coefficient has two temperature branches:

for $T\le 20$ °C, and

for $T>20$ °C, with depth correction

The sound speed $c$ used inside $A_1$ and $A_2$ is

For the displayed sound speed and latency, the calculator uses Mackenzie’s nine-term equation:

Solving $\mathrm{TL}\left(r\right)={\mathrm{TL}}_{\max }$ has a unique positive root because loss increases with range. The script brackets that root and uses bisection. It avoids the tempting update $r\to {10}^{({\mathrm{TL}}_{\max }-\alpha r)/N}$, which can diverge when absorption is important. Its derivative magnitude is $\alpha r\hspace{0.17em}\ln 10/N$, with instability above $N/\ln 10$ dB of absorption.

At the solved range, received level and propagation time are

The remaining link margin at any candidate distance is useful when evaluating a route rather than only the marginal endpoint:

A positive margin means the simplified link closes at that distance, while a negative margin indicates a shortfall. The percentage of the available transmission-loss budget consumed by absorption is

When the optional Thorp comparison is selected, the absorption coefficient in dB/km is calculated from frequency in kHz as

If absorption is neglected entirely, the geometric-spreading-only range provides a useful diagnostic ceiling:

Finally, one decibel of absorption-only allowance corresponds approximately to the following additional path length when the coefficient is locally constant:

Worked example: a 10 kHz AUV telemetry link

Consider a modem with $\mathrm{SL}=190$ dB re 1 µPa at 1 m, a 10 kHz carrier, 70 dB band noise and an omnidirectional receiver with $\mathrm{DI}=0$ dB. The demodulator needs $\mathrm{DT}=10$ dB. Use 10 °C, 35 ppt, 100 m depth, pH 8.0 and spherical spreading with $N=20$.

Francois-Garrison gives about $\alpha =0.950$ dB/km. The allowable transmission loss is

Solving $20{\log }_{10}\left(1000r_{\text{km}}\right)+0.950r_{\text{km}}=110$ gives $r=23.69$ km. Spreading contributes about 87.49 dB and absorption 22.51 dB. The marginal received level is 80 dB re 1 µPa, exactly 10 dB above the 70 dB noise. At roughly 1491 m/s, one-way latency is about 15.9 seconds.

Interpreting the calculated underwater range

The reported distance is where SNR exactly equals the selected threshold. It is not a guaranteed operating radius. At low frequency, geometry consumes most of the budget and added source level or receive gain can produce a substantial range increase. At high frequency, absorption grows linearly with distance, so extra power buys less reach. The absorption-share result helps identify which regime controls the link.

Read latency with range. Acoustic propagation is slow enough that acknowledgements, retries and control loops may take seconds or minutes. A link that closes acoustically may still be unsuitable for an interactive protocol.

Limitations of the underwater acoustic range model

This range-independent model assumes one spreading law, a uniform noise field and a representative water column. It does not calculate refraction, ducts, shadow zones, surface or seabed reflection loss, multipath fading, Doppler, transducer bandwidth, packet overhead or modem implementation loss. Real channels can therefore be much worse—or occasionally better—than this estimate.

Use the result as an engineering upper bound. For deployment planning, add a fading margin of roughly 6–12 dB or design for a conservative fraction of the estimated range. Important links should be checked with measured sound-speed profiles and a propagation model such as Bellhop or RAM.

Francois-Garrison is generally applied from 0.2 kHz to 1 MHz, but accuracy is not uniform across that span. The Thorp option is included for comparison with older low-frequency work and should not be trusted at high frequencies. Mackenzie sound speed is intended for approximately −2 to 30 °C, 30 to 40 ppt and depths to 8000 m; the calculator warns when it extrapolates beyond that envelope.

Questions engineers ask about acoustic modem range

Sources for the acoustic range equations

Sources. The sonar equation and spreading laws follow R. J. Urick, Principles of Underwater Sound, 3rd edition. Absorption follows R. E. Francois and G. R. Garrison, “Sound absorption based on ocean measurements,” Parts I and II, Journal of the Acoustical Society of America 72 (1982), including doi:10.1121/1.388673. Sound speed follows K. V. Mackenzie, “Nine-term equation for sound speed in the oceans,” JASA 70 (1981). The implementation was also compared with measurements reported by G. J. Macaulay, D. Chu and E. Ona, available from the NOAA Institutional Repository.

Transcription note. This implementation uses 8.86/c for the boric acid amplitude. Some secondary implementations use 8.696/c. For sub-decibel work below 1 kHz, consult the original paper and the conventions required by the project.