Abstract
A routing-benefited deployment algorithm combining static and dynamic layouts is proposed, and its comprehensive performance evaluation is given in this article. The proposed routing-benefited deployment algorithm is intended to provide a suitable network deployment and subsequent data transmission approach for underwater optical networking and communication. Static nodes are anchored for long-term monitoring, and movable nodes can adjust their depths based on the virtual force and move with the variation of area-of-interest changing. Then, nodes begin to collect data that they can monitor and transmit to sink nodes. Here, the underwater wireless optical communication model is described to actualize the real environment, and the vector-based forwarding protocol is particularly considered to compare the impact of different deployment algorithms on routing. It is shown by simulation experiment results that routing-benefited deployment algorithm outperforms several existing traditional virtual force deployment algorithms in terms of coverage, lifetime, energy consumption balance, packet-loss rate, and time-delay.
Keywords
Introduction
Underwater wireless sensor networks (UWSNs) are composed of homogeneous and/or heterogeneous sensor nodes. These nodes have the functions of data acquisition, processing, and wireless communication, which can be increasingly available for marine environment monitoring, disaster alarm and pollution detection, ocean resource exploration and exploitation, and military supervision.1,2
UWSNs have many differences from terrestrial sensor networks. First, the UWSNs are three-dimensional (3D) rather than two-dimensional (2D) networks. 3 Second, the underwater environment is complicated and fickle, and there are various disturbances such as ocean currents, marine life, seawater salinity, and temperature/pressure, which will seriously affect the performance of communication and networking. 4 Third, underwater is an untouchable and rugged environment. In this case, energy becomes a vital resource because changing batteries seems impossible. 5 In addition, the cost of underwater nodes is extremely high, which poses huge challenges to the deployment of UWSNs. 6
The deployment approaches not only assign the spatial positions of varying types of sensor nodes but also design the topology, even architecture of networks. Therefore, network deployment is a determinant of energy consumption, and the capacity and reliability of networks. 7
The deployment of underwater network nodes is divided into static deployment and dynamic deployment. Static deployment is simple and efficient, but it is only suitable for monitoring the fixed areas. The dynamically deployed nodes can adjust positions spontaneously to improve coverage ratio and connectivity ratio, but the cost of this kind of deployment is high and the mobile energy consumption is huge. 8 Therefore, the combination of static and dynamic deployments is more suitable for the actual situation of UWSNs, which can not only save energy consumption, prolong the network lifetime but also economize costs and improve robustness.
Regarding communication methods, the channels in UWSNs mainly include electromagnetic wave–based, acoustic wave–based, and optical wave–based ones. Since seawater contains salt which is a conductive transmission medium, the electromagnetic signal can only propagate for a short distance. 9 Besides, underwater wireless optical communications can implement several orders of magnitude higher propagation speed, lower time-delay, lower data-rate, and higher bandwidth than underwater acoustic communications.10–12 Therefore, optical channels are more suitable for high-speed, multi-service UWSNs. 13
Moreover, for the performance evaluation of network deployment methods, the current papers mainly consider the coverage and energy consumption performance of the network. Network deployment also greatly affects routing performance such as packet-loss rate (PLR) and time-delay.
Hence, this article proposes a routing-benefited underwater wireless sensor network deployment algorithm (RBDA). Some nodes are static nodes for long-term monitoring. The other nodes are movable, and they move in depth with the altering of area-of-interest according to the virtual force. The main contributions in this article include the following:
A deployment algorithm based on underwater wireless optical communication mode is proposed, which takes into account sunlight background noise and optical transmission loss model to make it more suitable for practical application scenarios.
Considering the network routing process, the VBF (vector-based forwarding) protocol is introduced, and routing performance is employed as one of the evaluation indicators of the deployment method.
The performances of the proposed RBDA, including network coverage (including its convergence speed), lifetime, energy consumption balance, PLR, and time-delay are evaluated and compared with some existing deployment algorithms.
The rest of this article is organized as follows. Section “Related work” provides a brief review of existing underwater deployment algorithms based on different methodologies. Some related preliminaries are outlined in section “Preliminaries.” A hybrid deployment algorithm with static and movable sensor nodes is proposed in section “Algorithm description.” The performance analysis and simulation of our deployment approach and the comparison with the existing several other methods are shown in section “Simulation results and discussion.” The summary and prospects are conducted in section “Conclusion.”
Related work
Regarding the current deployment approaches of UWSNs, they can be categorized into static and dynamic ones according to whether the nodes can move.
The core idea of the static deployment scenarios is using convex polyhedron to fill 3D space. Alam and Haas 14 used four polyhedrons, namely, cube, hexagonal prism, rhombic dodecahedron, and truncated octahedron to fill the 3D area. They proposed several useful definitions, such as volumetric quotient and sensing range. Then, they proposed the concept of cell and duty-cycled scheme to extend the lifetime of the networks in Alam and Haas. 15 Liu 16 further proposed a cluster-based deployment algorithm based on the aforementioned method.
The deployment methods of dynamic nodes can be designed based on graph theory, virtual force, and artificial intelligence, etc. 17 Among them, data centers are required to instruct other nodes in graph theory–based deployments, so it is difficult to achieve complete distributed deployment. Also, artificial intelligence methods tend to be complex in calculation. In contrast, the virtual force–based methods can ensure the distributed deployment of nodes, and the complexity is low. Therefore, they are more suitable for UWSNs composed of battery-powered nodes.
Zou and Chakrabarty 18 first introduced the theory of virtual force to the deployment to 3D environments and verified the feasibility. Li et al. 19 rewrote the virtual force formula, and introduced central gravity and balance force to jointly optimize the movement of nodes. Boufares et al. 20 gave a feasible method to derive the distance threshold for 3D deployments based on virtual force. Liu et al. 21 continued to improve the formula of virtual force and introduced the mobile benefit model to balance the energy consumption.
Most of the deployment algorithms for UWSNs are based on acoustic networks; there are only a few pieces of research works on the deployment of underwater optical wireless sensor networks. Xing et al. 22 and Wang et al. 23 proposed two types of underwater acoustic-optical hybrid wireless networks’ architectures, where the deployment of surface gateway nodes is focused in Xing et al., 22 and the network transmission protocols are concerned in Wang et al., 23 respectively. Yadav and Kumar 24 designed a clustering approach for UWSNs which combined the acoustic, optical, and electromagnetic wave–based communication techniques, but the placement of nodes merely obeys a 2D Gaussian distribution.
Preliminaries
Underwater optical wireless communication model
The underwater wireless optical communication system is briefly shown in Figure 1. The electrical signal to be transmitted is modulated onto a light source (LD or LED) to output an optical signal. Next, the optical signal is transmitted through the underwater channel to the receiver, then converted into an electrical signal by photo-detector (PD), at last, recovered by the demodulator.25,26

Schematic diagram of underwater wireless optical communication system.
Assuming that the power of the optical signal emitted by the light source of the optical transmitter is Pt and the optical power of the signal received by the PD is Pr. Their functional relationship can be expressed as
where d is the distance between the transmitter and the receiver. θ (half-angle beam width) and φ (angle between the optical axis of the receiver and the connection between the transmitter and receiver) are shown in Figure 1. ηt and ηr denote the optical efficiency of transmitter and receiver, respectively. Ar indicates the aperture area of the receiver, and k(λ) symbolizes the attenuation coefficient, which is influenced by the wavelength of the optical signal λ, as well as the general turbidity of the water.
The PD converts the optical power Pr to photocurrent i
where ℜ is the responsivity of PD, and it can be defined as
where ηPD denotes the quantum efficiency of PD, q indicates the electron charge, ħ is the Planck constant, and c represents the propagation velocity of light in vacuum.
Substituting equations (2) and (3) into equation (1), then
Also, the signal-to-noise ratio (SNR) for the underwater optical wireless communication system can be expressed as
where N0 is the noise equivalent power, composed of solar background shot noise power (Nb), signal shot noise power (Ns), dark current noise power (Nd), and preamplifier noise power (Np), that is
where the value of Nb changes with the depth of the receiver in water D, the wavelength, and power of the emitted optical signal. The value of Ns also changes with the wavelength and power of the transmitted optical signal. Nd and Np can be regarded as constants in the simulation for the sake of simplicity since they are not affected by the transmitter, receiver, and the transmitting signal. 27
Applying the Shannon theorem based on known system bandwidth (W) and SNR, the data transmission bit rate (Rb) can be attained. Then, the optical energy Et consumed by each transmitter for transmitting each symbol (bit) is
where Pt can be calculated by equation (4).
3D virtual force model
The basic mechanism of virtual force is that a repulsive force will be produced when two neighbor nodes are close to each other, or an attractive force will be produced when two neighbor nodes are far away from each other. In this case, nodes can adjust the distances between their neighbors and themselves to optimize spatial distribution and coverage performance of the whole network. The preinstall threshold for judging whether the distance between neighbor nodes is close or far is recorded as dth.
The coordinates of ni and nj are (xi, yi, zi) and (xj, yj, zj), respectively, and the virtual force exerted by sensor nj on sensor ni is expressed as
where rc is the communication radius of the node, Kr and Ka represent the repulsive coefficient and the attractive coefficient, respectively. dij represents the Euclidean distance between the two nodes.
In practical terms, realizing the 3D free movement of the nodes is too costly; it is more of engineering significance to assume that the nodes can only adjust their depths. 28 Hence, the moving distance is directly related to the projection of the virtual force in the vertical direction. αij is introduced as
and then, the projection of virtual force in vertical direction exerted by sensor nj on sensor ni is
Besides the virtual force between nodes, when the node is close to the upper or lower boundary, a repulsive force can be also produced to prevent it from moving out of the monitoring space. The direction of this type of repulsive force is away from the boundary, and its magnitude is
where Kb is the repulsive coefficient of the boundary, dbth is the threshold distance determining whether the repulsive force will be exerted, dbi is the depth difference between node ni and the closest boundary of the monitoring area.
According to equations (11) and (12), the total virtual force exerted on ni is
where N is the total number of nodes.
The node moves upward and downward under the virtual force, the distance and direction of movement are
where lmax is the preinstall maximum moving distance with each iteration to avoid high moving energy consumption, and β (unit: s2/kg m) indicates the adjustment coefficient. ω is the mobile benefit coefficient designed for the sake of node mobility effectiveness and energy balance of networks. 21
Underwater sensor network routing
The network layer protocol of UWSNs is one of the research hotspots for the past few years. The commonly used VBF protocol is introduced in this article to evaluate the routing performances of the deployment algorithms.
The VBF is a routing protocol based on geographic location, and its basic principle is shown in Figure 2. 29 nso is the source node and nsi is the sink node. Assuming that the water surface gateway nodes are sink nodes; they are evenly distributed on the water surface. The data collection nodes are source nodes, and every source node matches a sink node. The node ni is the transmitting node, nj is one of the upper-level neighbor nodes of ni, and the forwarding factor µ(i, j) of nj relative to ni can be expressed as
where L is the distance from node nj to the connection between the source node and the sink node, and δ is the angle between the connection from ni to nj and the connection between ni to sink node. dij represents the distance from ni to nj, and κ is the adjustment factor.

Schematic diagram of the VBF protocol.
Before transmitting the packet, the transmitting node calculates the forwarding factors of its upper-level neighbor nodes, and selects the upper-level neighbor node with the smallest forwarding factor as the receiving node.
PLR and average transmission time-delay are the two important routing performance parameters. The PLR is defined as
where Nc is the number of packets successfully collected and transmitted, and Nl is the number of packets unsuccessfully collected or transmitted.
Time-delay (t) refers to the time required for the network to successfully collect and transmit a data packet to the sink node. It is the sum of transmission delay (tt), propagation delay (tp), processing delay (tpr), and queuing delay (tq), that is
Assuming that the size of each data packet collected is equal to lp (bit), and the queuing delay is negligible. Thus
where Nr is the number of relay nodes. The 0th relay node is the source node, and the (Nr + 1)th relay node is the sink node. τt and τpr denote the transmission delay and the processing delay of the node to send a single data packet, respectively. dk(k + 1) is the Euclidean distance between the kth relay node and the (k + 1)th relay node. It should be mentioned that the propagation speed of light is approximately equal to 0.75c because of the underwater environment.
The average time-delay is defined as the ratio of the total delay (T = Σt) to the number of packets collected (Nc) during the lifetime of networks.
Algorithm description
The proposed algorithm can be divided into four main stages: deployment, redeployment, data collection, and data transmission.
Deployment and redeployment approach
3D network topology of underwater optical sensor network is depicted in Figure 3. Beneath the water surface, nodes transmit data through optical links with other nodes, AUV 30 (autonomous underwater vehicle) or submarine, and surface gateway nodes. The data transmission between nodes is mostly from bottom to top, and finally, to the surface gateway nodes for the sake of convenient and safe data storage and processing. 31 Over the water surface, surface gateway nodes can directly or indirectly communicate with buoys, ships, ground stations, satellites, and so on, via radio (including radio frequency links and satellite communication links).

3D network topology of underwater optical sensor network.
The underwater nodes are divided into static nodes and movable nodes. After filling the monitoring space with rhombic dodecahedrons, static nodes are anchored at the center of each rhombic dodecahedron for long-term monitoring. The movable nodes are randomly placed on the water surface. This process is called deployment.
In order to simplify the analysis, it is assumed that the gravity and buoyancy received by movable nodes are the same in size and the opposite in direction. The movable nodes adjust their depths on the basis of their exerted virtual force and move with the change of the area-of-interest. This process is called redeployment.
Besides different mobility of the two types of nodes, the rest of the functions and parameter setting are homogeneous, including as follows:
Data collection with sensing radius rs (types of data include environmental parameters such as pressure and temperature, and services such as voices, images, and videos).
Data transmission with communication radius rc.
Data storage, processing, and computing.
Initial energy E0.
The movement of movable nodes consumes energy Em
where em is the energy consumption for moving per unit distance, and |
Algorithm 1 is the pseudocode of the detailed process redeployment.
It should be noticed that Line_5 shows the process of information exchange between neighbor nodes: each node broadcasts a “Hello” message that includes its ID and coordinate. At the same time, each node receives its neighbors’“Hello” messages, and then calculates the distance between them, and adds this distance and neighbors’ IDs and coordinates into its storage. This process follows the Request To Send/Clear To Send (RTS/CTS) handshake protocol in medium access control (MAC) layer. 32
After the redeployment, the network coverage ratio C is calculated. It is defined as the ratio of the space that can be monitored by all nodes to the total monitoring area
where V is the total volume of the monitoring area that is usually set as a rectangular space. Vn(i) is the volume of the area that can be monitored by the node ni. Assuming that the nodes are equipped with omnidirectional antennas, so their monitoring area is a sphere area with the node as the sphere center and rs as the radius.
If the coverage of the network is not convergent after adjustment, the next redeployment process continues. After the coverage ratio is converged (the numerical difference in coverage ratio between before-adjusted and after-adjusted network is less than or equal to 0.2%), the movable node will stay stable and the network starts to collect and transmit data.
Data collection and transmission
Only the data collection and transmission under the water is considered in this article. It is supposed that data randomly appears in the monitoring space, which will be collected by the nearest node. After simple processing and calculation, the data are transmitted from the bottom to the surface gateway node based on VBF protocol.
Algorithm 2 is the pseudocode of the detailed process of data collection and transmission.
First, check whether the lifetime ends. When too many packets are lost or too many nodes stop working, the network’s lifetime ends (Lines_2∼11). Next, the data packets will be collected (Lines_12∼20). But if the data packets appear outside the coverage of the network or the residual energy of the collecting node is not enough, the data packet will lose (Lines_14∼17). It should be noted that the signal for data collection is sent and received by the same sensor, so the distance of optical transmission is twice the distance from the node to the location of the packet (Line_18). After the data being successfully collected, data transmission will be performed according to the VBF protocol (Lines_21∼38). If there is no upper-level neighbor node (Lines_27∼29) or the residual energy of relay node is not enough (Lines_32∼34) during the transmission process, the packet-loss happens.
Simulation results and discussion
As a calculation example, the monitoring area is set to a cube area of 200 m × 200 m × 200 m. Specific parameter settings are shown in Table 1. The performance of RBDA is evaluated in five aspects: coverage ratio, lifetime, residual energy variance, PLR, and average time-delay, which are also compared with several existing classic deployment methods.18,19,21
Parameter setting.
Coverage ratio
The virtual force model is introduced to increase network coverage ratio. A total of 200 nodes (100 static nodes and 100 movable nodes) are deployed in the monitoring area. The network coverage ratio of different algorithms during the redeployment process with the number of iterations is shown in Figure 4.

Network coverage ratio versus the number of iterations.
RBDA is significantly better than other algorithms in terms of coverage ratio, as well as the rate of convergence. The coverage ratios of the four algorithms when having been converged are 57.35% for virtual force algorithm (VFA), 58.18% for improved-VFA, 58.24% for distributed node deployment algorithm based on virtual force (DABVF), and 60.89% for RBDA, respectively. Besides, the numbers of iterations are 23 times for VFA, 20 times for improved-VFA, 18 times for DABVF, and only 7 times for RBDA when the coverage ratio of the network is converged. In this case, the RBDA can achieve the same coverage ratio as the purely dynamic deployment method with fewer nodes, and therefore, it is more cost-effective.
Lifetime and residual energy variance
The number of data packets collected within normally working state of the network is called lifetime, measuring the utilization rate of network energy. After the network stops working (lifetime ends), the variance of the residual energy of all nodes is calculated, which measures the balance of energy consumption of the network. 33
Figure 5 shows the lifetime of different algorithms with the initial number of nodes obtained by simulation. In the early deployment process of the RBDA, the moving distance of the nodes is shorter, the energy consumption is lower, and the coverage rate of the network is higher, so the lifetime of RBDA nearly doubles the other algorithms.

Network lifetime versus the number of nodes.
Figure 6 shows the variance of the residual energy of different algorithms with the number of nodes. In general, the residual energy variance decreases with the number of nodes. RBDA and DABVF both introduce mobile benefit coefficient, so the residual energy variances are much lower than the other two algorithms, indicating that network energy consumptions of these two algorithms are more balanced.

Residual energy variance versus the number of nodes.
Package-loss rate and average time-delay
The PLR and the time-delay are typical network routing performances. The routing performances between the RBDA and several existing deployment algorithms are simulated and compared using the same VBF protocol.
The PLRs of different algorithms with the number of nodes are shown in Figure 7. It can be seen that the PLR of RBDA is much lower than (nearly half of) the other three algorithms when the number of initial nodes is the same.

Packet-loss rate versus the number of nodes.
The average time-delay of different algorithms with the number of nodes is shown in Figure 8. It can be seen that the average time-delay gradually decreases with the number of nodes. In the case of the same number of nodes, the RBDA has the lowest average packet-delay.

Average time-delay versus the number of nodes.
From Figures 4–8, regarding the common criteria of the deployment algorithms, RBDA can achieve a higher and more fast-converged coverage ratio of networks, a longer lifetime, and a more balanced energy consumption. Besides, RBDA is a routing-benefited deployment algorithm, as two important performances of routing, namely, package-loss rate and time-delay, can be greatly improved by RBDA.
Conclusion
Network deployment is important for efficient networking and reliable communication. We proposed RBDA for the underwater wireless optical sensor network using the virtual force model and VBF protocol. Also, RBDA is compared with relevant algorithms. The numerical results show that it outperforms in coverage, energy, and routing performances.
In the future, more environmental factors should be considered in the deployments, such as water currents and vortices. Localization algorithms and transmitting protocols combined with specific deployment approaches also have research value. Besides, research works for underwater wireless optical communication technologies should be more in-depth to further accelerate the practical application.
Footnotes
Handling Editor: Yanjiao Chen
Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by the National Natural Science Foundation of China (NSFC) under grant nos 61871418 and 61801079, China.
