World Ocean Simulation System (WOSS) library
World Ocean Simulation System (WOSS) Library
Author
Federico Guerra - feder.nosp@m.ico@.nosp@m.guerr.nosp@m.a-tl.nosp@m.c.com
Version
1.13.0

This document provides a short techical description of the World Ocean Simulation System (WOSS) library and of its integration into Multi InteRfAce Cross Layer Extension (NS-Miracle).

Library overview

The WOSS library aims to model a channel power delay or frequency attenuation profiles according to a more accurate representation of propagation phenomena as described by the laws of acoustical physics [3]. WOSS is a multi-threaded framework that permits the integration of any existing underwater channel simulator that expects environmental data as input and provides as output a channel realization represented using the channel profiles just mentioned.

From version 1.3.5 the simplified acoustic model based on empirical equations for the calculation of delay, attenuation and noise power has been moved into NS-Miracle. These equations can be found in [1], [2].

Currently WOSS integrates the Bellhop ray-tracing program [4], while the NS-Miracle integration library retains the previously mentioned formula for the noise power.

From version 1.3.0 WOSS supports time evolution, for a better representation of the channel variability. Bellhop calculations are performed across the whole bandwidth with custom resolution, while noise calculations are performed at the geometric average of the highest and the lowest frequencies of the band being used, as pointed in [5]. To calculate the solution of the propagation equations between a transmitter and a receiver, Bellhop requires the knowledge of the Sound Speed Profile (SSP) (The Sound Speed Profile is the propagation speed of sound considered as a function of water depth. Different profiles lead to potentially very different propagation effects, including surface sound channels, deep sound channels, convergence zones, shadow zones, and so on, see [3]), the Bathymetric Profile (BP) and the type of Bottom Sediments (BS, required to model acoustic power losses due to bottom reflections). In this respect, WOSS offers a technology independent database API and a number of classes to manipulate SSP, BP and BS data. For the SSP, it employs a custom implementation of the World Ocean Database [7-11], a collection of SSPs calculated with the TEOS-10 [6] exact formula based on the depth, temperature and salinity data measured during a number of experiments around the world; the measurements are divided by location and day or season of the year when the measurement was performed (recall that sound propagation is affected by water temperature, which in turn undergoes seasonal changes, especially in the superficial layer). The bathymetry data have been taken from the General Bathymetric Chart of the Oceans [12], a public database offering samples of the depth of the sea bottom with an angular spacing of 30 seconds of arc. Finally, the type of bottom sediments is provided by a geo-acoustic analysis of the National Geophysical Data Center's Deck41 data-base [13].

Thanks to this automation the user only has to specify the location in the world and the time where the simulation should take place. This is done by setting the simulated date and the wanted latitude and longitude of every node involved. The simulator automatically handles the rest. In more detail, the simulator picks the location (i.e., latitude, longitude and depth) of the transmitter and the receiver and queries the database manager for samples of bottom sediments, measured SSPs (for simplicity, the SSP can also be assumed to be constant, on average, throughout the network area) and for bathymetry data along the path. Full customization of surficial sediment, bathymetry and SSP data is also possible if more accurate data is available. The power delay profile obtained with Bellhop can be used to model multiple receptions: in detail the user can choose to coherently combine the complex channel gains with a custom resolution time window and obtain one or more complex channel taps.

Figs. 1.1 and 1.2 show the coherent sum of all channel taps, obtained with WOSS for an example scenario. Fig. 1.1, represents the attenuation incurred by an acoustic wave at 4 kHz transmitted in August, approximately 20 km offshore the harbor of the city of Taranto, Southern Italy, with the transmitter located at 40.32N, 17.12E, 300 m depth. Darker shades of gray represent stronger signal power. The figure shows that the signal reaches the top right corner surface, about 11 km from the transmitter, bearing sufficient power to allow correct reception. In contrast, Fig 1.2 shows the same scenario in March: in this case, the average temperature of the water is lower, changing the way the sound is refracted and reflected by the sea bottom. Here, any acoustic receiver deployed between 9 and 11 km from the transmitter may be unable to acquire the transmitted signal. Furthermore, a shadow zone of about 2 km appears in front of the transmitter, causing failure in a hypothetical short range communication scenario.

Fig. 1.1
Fig. 1.2

Technical Overview

In this section, a short techical overview of the WOSS library and its APIs is provided.

The Foundation Classes

WOSS foundation is a complete set of classes for:

  • geographical coordinates [woss::Coord, woss::CoordZ]: measured in decimal degrees, this classes offers methods for distance, bearing, and marsden square calculations and other mathematical formulas;
  • date time [woss::Time]: date time arithmetics;
  • simulation time reference [woss::TimeReference]: importing a simulation time reference from the external network simulator;
  • random number generation [woss::RandomGenerator]: importing a custom random generator from an external network simulator or library;
  • sound speed profile [woss::SSP]: sound speed calculations and manipulations with equations taken from the scientific literature;
  • surficial sediment [woss::Sediment]: geoacoustic sediment parameter definitions and manipulations;
  • channel power delay profile [woss::TimeArr]: profile sampling, averaging etc.
  • pressure [woss::Pressure]: arithmetic and conversion operators;
  • transducer [woss::Transducer]: accurate power consumption calculations and beam pattern capabilities of real transducers;
  • altimetry [woss::Altimetry]: a base class that should be used for implementing and simulating different sea surface wave spectra.
  • a definitions handler class for "plugging-in" at run-time any class that inherits from the ones described above [woss::DefHandler].
  • test suite [woss::WossTest]: a base class that should be used for writing tests of specific WOSS framework APIs.

The Database APIs

The woss::WossDb API separates database technology implementation (SQL, NetCDF, textual…) from data behaviour implementation, permitting maximum code modularity and re-usability. Data behaviour are defined with WossBathymetryDb, WossSedimentDb, WossSSPDb.

The woss::WossDbCreator API provides interface for creation and initializiation of WossDb objects, providing the user a ready-made database instance.

Custom defined bathymetry, SSP and sediment databases are accessed through the woss::WossDbManager API, to provide a further refinement of data (eg. averaging, arithmetics…).

The woss::WossResPressDb and woss::WossResTimeArrDb classes provide the possibility to store and retrieve channel simulation results in any user defined fashion.

The Channel Simulator APIs

The woss::Woss API has the task of integrating a channel simulator. It gives interface for:

  • setting up the enviroment from databases or from custom data;
  • configuring, running and reading results of any channel simulator;
  • averaging and giving output results over any number of runs and frequencies;
  • optionally reading channel simulator's results [woss::ResReader, woss::WossResReader].

The woss::WossCreator API provides interface for creation and initializiation of woss::Woss objects, relieving the user from this task.

The woss::WossManager API has the logic to control and manipulate and create woss::Woss objects. It offers the abilities to:

  • carefully create and utilize Woss objects to minimize CPU load;
  • plan a strategy for time-varying and/or multi-frequency channel simulations;
  • use on-the-fly custom bathymetry, sediment or SSP data for any tx-rx pair.

Finally the woss::WossController API has the task of inter-connecting all major APIs involved.

WOSS block diagram

WOSS block diagram

Current capabilities

The current version provides:

See also
Database Descriptions for further info on custom database, Underwater Acoustic Transducers for further info on transducers, Altimetry models for further info on altimetry models, Time evolution modeling for further info on time evolution modeling.

References

  1. R. Urick, Principles of Underwater Sound. McGraw-Hill, 1983.
  2. M. Stojanovic, “On the relationship between capacity and distance in an underwater acoustic communication channel,” ACM Mobile Comput. and Commun. Review, vol. 11, no. 4, pp. 34–43, 2007.
  3. F. Jensen, W. Kuperman, M. Porter, and H. Schmidt, Computational Ocean Acoustics. Springer-Verlag, 2000.
  4. Bellhop code, http://oalib.hlsresearch.com/Rays/index.html .
  5. P. C. Etter, Underwater acoustic modeling and simulation. Spon Press, Taylor & Francis group, 2003.
  6. TEOS-10 source code http://www.TEOS-10.org .
  7. World Ocean Atlas 2001, https://www.nodc.noaa.gov/OC5/WOA01/pr_woa01.html .
  8. World Ocean Atlas 2005, https://www.nodc.noaa.gov/OC5/WOA05/pr_woa05.html .
  9. World Ocean Atlas 2009, https://www.nodc.noaa.gov/OC5/WOA09/pr_woa09.html .
  10. World Ocean Atlas 2013, https://www.nodc.noaa.gov/OC5/woa13/ .
  11. World Ocean Atlas 2018, https://www.nodc.noaa.gov/OC5/woa18/ .
  12. World Ocean Atlas 2023, https://www.nodc.noaa.gov/OC5/woa23/ .
  13. General bathymetric chart of the oceans, http://www.gebco.net .
  14. National geophysical data center, seafloor surficial sediment descriptions, http://www.ngdc.noaa.gov/mgg/geology/deck41.html .