Connectivity on Underwater MI-Assisted Acoustic Cooperative MIMO Networks

In traditional underwater wireless sensor networks (UWSNs), it is difficult to establish reliable communication links as the acoustic wave experiences severe multipath effect, channel fading, and ambient noise. Recently, with the assistance of magnetic induction (MI) technique, cooperative multi-input-multi-output (MIMO) is utilized in UWSNs to enable the reliable long range underwater communication. Compared with the acoustic-based UWSNs, the UWSNs adopting MI-assisted acoustic cooperative MIMO are referred to as heterogeneous UWSNs, which are able to significantly improve the effective cover space and network throughput. Due to the complex channel characteristics and the heterogeneous architecture, the connectivity of underwater MI-assisted acoustic cooperative MIMO networks is much more complicated than that of acoustic-based UWSNs. In this paper, a mathematical model is proposed to analyze the connectivity of the networks, which considers the effects of channel characteristics, system parameters, and synchronization errors. The lower and upper bounds of the connectivity probability are also derived, which provide guidelines for the design and deployment of underwater MI-assisted acoustic cooperative MIMO networks. Monte Carlo simulations were performed, and the results validate the accuracy of the proposed model.

The performance of the underwater MI-assisted acoustic cooperative MIMO networks is obtained in [21], in which the performance gain of underwater MI-assisted acoustic cooperative MIMO networks is validated in terms of synchronization errors, signal-to-noise ratio (SNR), effective communication time, and the upper bound of the throughput in physical layer. In this paper, we focus on the system level performance of underwater MI-assisted acoustic cooperative MIMO networks in terms connectivity and analyze the impact of clock synchronization on connectivity. Because of the complex channel characteristics and the heterogeneous network architecture, the connectivity analysis of MI-assisted acoustic cooperative MIMO networks is much more complicated than that of UWSNs. In particular, the acoustic cooperative MIMO and MI communications are performed in inter-cluster and intra-cluster manners, respectively. Therefore, the underwater acoustic channel and MI channel should both be considered, which increases the complexity of connectivity analysis of the networks. Besides the complex channel characteristics, the heterogeneous architecture of MI-assisted acoustic cooperative MIMO networks also poses a great challenge for connectivity analysis. Compared with the traditional UWSNs [7,22], the inter-cluster and intra-cluster connectivity need to be calculated. Meanwhile, the synchronization error is inevitable in the connectivity analysis of the networks.
Without loss of generality, the networks are supposed to be deployed in the bounded three-dimensional space R 3 . In this underwater MI-assisted acoustic cooperative MIMO networks, the clusters are distributed according to a homogeneous Poisson point process and the nodes in a cluster are distributed according to another homogeneous Poisson point process. When considering the probability that a certain cluster can be directly connected with surface BS, it is not only related to its position in space R 3 , but also related to the maximum transmission range of the cluster, which depends on the number of nodes in the cluster [23]. When considering the multi-hop fashion, some work has been well explored to analyze the connectivity of wireless sensor networks assuming two-dimensional coordinates [24]. We extend it to the three-dimensional scene in this network model and give the upper and lower bounds of connectivity.
The major contributions of this paper can be summarized as follows: • A mathematical model is developed to analyze the connectivity of underwater MI-assisted acoustic cooperative MIMO networks, which considers the three-dimensional underwater bounded environment. Connectivity probabilities of direct and multi-hop fashions are derived.

•
Due to the characteristics of heterogeneous networks, inter-cluster and intra-cluster connectivities need to be calculated, and the influence of synchronization errors on connectivity is also considered. Monte Carlo simulations were carried out, and the results agree well with the analytical model that validate the accuracy of the model.

•
Based on the framework, we present numerical and simulation results that quantitatively analyze the effects of various system parameters on the connectivity in this networks, considering the effects of channel characteristics, operating frequency, operating range, and density of nodes. The results provide the guideline for the deployment of underwater MI-assisted acoustic cooperative MIMO networks according to the practical requirement.
The remainder of this paper is organized as follows. In Section 2, a brief introduction to the network architecture of underwater MI-assisted acoustic cooperative MIMO networks is given. Based on the system architecture, the connectivity problem is formulated in Section 3. In Section 4, the transmission ranges of node and cluster are derived. Next, connectivity probabilities of direct and multi-hop fashions are further analyzed in Section 5. In Section 6, the connectivity performance evaluation by Monte Carlo simulations is presented. Finally, the conclusions of this paper are drawn in Section 7.

System Model
The complex system architecture constitutes one of the major challenges in the connectivity analysis. In this section, we describe the system architecture and channel model, and the synchronization accuracy is analyzed.

System Architecture
The system architecture of underwater MI-assisted acoustic cooperative MIMO networks is illustrated in Figure 1. It consists a surface base station (BS) and a great number of underwater sensor nodes. The surface BS aggregates the data from underwater nodes for environment sensing, which includes water velocity, temperature, salinity, etc. Underwater sensor nodes are managed in the form of clusters with assigned IDs. In each cluster, the sensor selected as the master node collects the monitoring data from slave nodes and coordinatea the intra-cluster synchronization. The acoustic cooperative MIMO is utilized between the surface BS and clusters as well as different clusters for information exchange and command interaction. As the cooperative MIMO technique collects all the intra-cluster nodes power, the long range underwater acoustic communication is enabled, which expands the effective cover space. Compared with traditional UWSNs [1], the cluster of the networks can be considered as an enhanced node due to the adoption of MI-assisted acoustic cooperative MIMO.
The process diagram of MI-assisted acoustic cooperative MIMO networks is shown in Figure 2. It includes five phases: Data request, Aggregation, Broadcast, Synchronization, and Communication. During the Data request period, the BS sends the data request to clusters. As the cluster receives the request, the monitoring data are transmitted from slave nodes to the master node by MI communication during the Aggregation period. The master node aggregates the monitoring data from different slave nodes and then broadcasts them to all the slave nodes during the Broadcast period. As a result, all the intra-cluster slave nodes share the same monitoring data. Timing-sync Protocol for Sensor Networks (TPSN) is adopted to adjust the clocks of slave nodes to that of the master node, which is operated in the Synchronization period [25]. After the Synchronization period, all the nodes within a cluster have the same monitoring data as well as the same clock. Finally, the monitoring information converged to the surface BS via the direct link or the multi-hop fashion in the Communication period.

Synchronization Aggregation
Broadcast Communication Data request

Channel Characteristics
There are two types of channel in the underwater MI-assisted acoustic cooperative MIMO networks: the underwater acoustic channel and the magnetic induction channel.
In underwater MI-assisted acoustic cooperative MIMO networks, underwater acoustic communication is adopted for communication between clusters in terms of cooperative MIMO. In this paper, we consider Thorp's underwater acoustic channel model [26], which is based on empirical data. Thorp's model provides the attenuation for the underwater acoustic channel, which is influenced by the frequency and the communication distance between the sender node and the receiver node. The transmission power model is formulated as [26] where P r,ac is received power, P t is the transmitted power, and A(d, f ) is the channel attenuation depending on the distance d and the frequency f (kHz) where A 0 is the normalizing constant and θ is the spreading factor. The spreading factor has a value between 1 and 2 depending on the depth. a( f ) is the absorption coefficient and is expressed in decibels per kilometer 10 log a( f ) = 0.11 In underwater MI-assisted acoustic cooperative MIMO networks, MI is adopted for information aggregation and synchronization within a cluster. As a near-filed communication technique, the transmission power model is formulated as [27] where P r,MI is received power, P t is the transmitted power, and r is the distance of the master node and one slave node; ρ = Q T Q R η T η R r 3 T r 3 R π 2 G where Q T = ωL T /R LT and Q R = ωL R /R LR are quality factors of transmitter and receiver, respectively; L T and L R are the inductances of the transmitting and receiving coils; R LT and R LR are the transmitting and receiving coil resistances; η T = R S /(R S + R LT ) and η R = R L /(R L + R LR ) are efficiencies of the transmitter and receiver coils; R S and R L are source and load resistances; r T and r R are the radius of the transmitting and the receiving coil; ω is resonant frequency; G = e −r/δ is a parameter related to the channel state; and δ is conductivity. For sea water, fresh water, and dry oil, δ is δ = 4.5 S/m, δ = 0.07 S/m, and δ = 0.0048 S/m, respectively.
Highly conductive sea water induces significant eddy current that incurs very high path loss. Therefore, in seawater, the path loss of MI communication also needs to take into account the influence of medium loss [28]. The attenuation coefficient is where µ 0 is vacuum permeability. Therefore, the path loss caused by the influence of seawater medium PL α (dB) is Finally, the path loss of MI communication in seawater is

Synchronization Accuracy Analysis
For the TPSN synchronization algorithm [25], the clock drift of two nodes is assumed to be unchanged in this small span of time. As shown in Figure 2, k − 1 slave nodes adjust their clocks to that of the master node in turn. The node z is taken as an example to illustrate the TPSN synchronization process. The node z sends synchronization-pulse packet to the master node at T 1 , and this package contains time T 1 . The master node receives the packet at T 2 , where T 2 is equal to T 1 + ∆ 1 + ∆ 2 . Here, ∆ 1 and ∆ 2 represent the clock drift between the two nodes and propagation delay, respectively. The master node returns an acknowledgement packet to the node z at T 3 , which contains the values of T 1 , T 2 , and T 3 , as well as beacon signal f s generated by the master node. The node z receives the packet at T 4 . The clock drift and propagation delay of the node z can be calculated as As the drift is indicated in Equations (8) and (9), the node z is able to correct its clock aligned with that of the master node.
However, the node's clock is driven by crystal oscillation. Doe to the existence of impure crystals, environmental changes, and component aging, frequency and synchronization errors are unavoidable.
First, the frequency error of the synchronization is analyzed herein [21]. The average relative clock drift in ∆T of the node z to the master node can be expressed as where f c,z (t) is the oscillator frequency of the node z and f c,0 (t) is the oscillator frequency of the master node. For node z, its working frequency f s,z (t) is proportional to its crystal oscillator frequency, thus where K is the amplification factor of the frequency multiplier. After receiving the acknowledgement packet at the node z, the frequency is estimated according to the local oscillator of the nodef where f s is the working frequency of the master node and ε s,z is the error of the frequency estimation.
The frequency offset at the node z for the frequency synchronization can be determined as According to Equation (13), the optimal frequency offset is (a z − 1) f s . Due to the estimated error ε s,z , the frequency cannot be perfectly synchronized and the error can be defined as For the node z, the duration of the time slot is ∆t z , which mainly includes the sending time and transmission time of synchronization-pulse packet and acknowledgement packet where L z,0 and L 0,z are the size of synchronization-pulse packet and acknowledgement packet, B is the bandwidth, r 0,z is the distance between master node and node z, and v is the propagation speed of signal. Due to the frequency error of crystal oscillator at each node, the total drifting time of each node will gradually increase along with time. For the node z, the time synchronization error to the beginning of beamforming can be calculated by

Problem Formulation
Based on the system architecture and operation mechanism of underwater MI-assisted acoustic cooperative MIMO networks given in Section 2, the problem of connectivity analysis is formulated in this section.
Without loss of generality, the network is supposed to be deployed in a bounded 3D space R 3 . In this area, the clusters {N j , j = 1, 2, · · · } are regarded as independent units, which are distributed according to a homogeneous Poisson point process with the density of ρ c . The nodes {N i , i = 1, 2, · · · } in one certain cluster are subject to another homogeneous Poisson distribution with the density of ρ n .
In MI-assisted acoustic cooperative MIMO networks, a cluster with k nodes and the surface BS can successfully establish a connected link if and only if the following two conditions are satisfied: (1) k − 1 slave nodes can connect topologically with the master node; and (2) the cluster with k nodes participating in cooperative MIMO and the surface BS can connect topologically via a direct link or a multi-hop fashion. Specifically, the definitions of node connected and cluster connected are provided as follow [24]: When there is a master node within the maximum transmission range of the slave node, it is considered that it can connect to the master node, and the slave node is defined as connected.

Definition 2.
When the surface BS is within the maximum transmission range of the cluster, it is considered that it can connect to the surface BS directly, and the cluster is defined as connected.

Definition 3.
When there are relays that can connect topologically with each other until the surface BS, it is considered that it can connect to the surface BS in multi-hop fashion, and the cluster is defined as connected.
From the above definitions, the connectivity of the jth cluster in the networks can be defined as where P C,k,j is probability that the jth cluster including k nodes connects to the surface BS, denoted as Equation (18), and P N,k is the probability that there are k nodes participating in the cooperative MIMO in the jth cluster, indicated as Equation (19). Specifically, where x j is the position of the jth cluster, and P(x j ↔ surfaceBS) is the probability that the jth cluster at x j is connected to the surface BS, which depends on the distribution and maximum transmission range of the cluster. As the maximum transmission range of the cluster depends on the number of nodes participating in the cooperative MIMO, the maximum transmission range of the cluster with k nodes is recorded as R ac,k . According to Definition 1, we have where S = πR 2 MI and R MI is the maximum transmission range of the node. Substituting Equations (18) and (19) into Equation (17), we have From Equation (20), to derive the connectivity probability P C,j , the following terms are requested: the maximum transmission range of the node R MI , the maximum transmission range of the cluster R ac,k , and the probability that the jth cluster at x j connects to the surface BS P(x j ↔ surfaceBS). These are analyzed in the following sections.

Inter-Cluster Acoustic Transmission Range
It is assumed that there are k nodes distributed in the jth cluster. Let τ i denote the time delay of the ith sensor, where i = 1, 2, · · · , k. The signal from the k nodes arrives at the surface BS can be written as where H i is the amplitude of the ith transmitted signal, h i = 1 A(d i , f ) is the underwater acoustic channel envelope [26], f is the signal frequency, and n(t) is the noise. Similar to the assumption in [21,29], as k nodes share the same transmitted signal, the amplitude H i can be recorded as H.
Considering the surface BS deployed in the far-field, all the channel envelops are the same.
Therefore, the signal at the surface BS can be rewritten as We consider the perfect CSI information, the channel envelope h i and the delay τ i are known. From Equation (22), the received signals by surface BS cannot be enhanced as the phase difference caused by the time delay. To maintain the reliable communication, the clusters should be close to the surface BS that leads to a limited effective cover space. To tackle this problem, pre-compensation phases are added to the nodes to eliminate the time delays, and this technique is referred as cooperative MIMO. After the phase compensation, the received signal at the surface BS can be expressed as where w i ∈ C is the pre-compensation phase and λ is the wave length. Define w = [w 1 · · · ω k ] T , a φ = [ e −j2π f τ 1 . . . e −j2π f τ i ] T , and σ 2 = E{|n(t)| 2 }. The received signal y(t) can be expressed in a compact form [29]: Then, the SNR of received signal at the surface BS can be calculated as The optimal compensation vector w can be obtained by solving the following optimization problem: The optimization problem in Equation (26) can be solved by the Lagrangian multiplier method, and the optimal compensation vector is given as below [29] The above analysis does not consider the synchronization error. Due to the existence of the impure crystals, component aging, and time drifting, frequency and time estimation errors are unavoidable. Taking account of the frequency and time estimation errors, the signals received at the surface BS can be expressed as where h = [Hhe j2π( f +ε 1, f )(t+ε 1,t ) · · · Hhe j2π( f +ε k, f )(t+ε k,t ) ], • is Hadamard product. ε i, f is frequency estimation error of the ith node and ε i,t is the time estimation error of the ith node, which are explained in Section 2. Therefore, the SNR of signals received at the surface BS is given as [21] SNR ac, Based on Equation (29) the maximum transmission range of a cluster is determined as where SNR ac,th indicates the minimum SNR that the surface BS is able to receive and recover the transmitted signal from the cluster.

Intra-Cluster MI Transmission Range
In underwater MI-assisted acoustic cooperative MIMO networks, within a cluster MI is adopted for information aggregation and synchronization. To obtain the intra-cluster transmission range, the SNR of MI is needed. From the transmission power model in Equation (4) in Section 2, the SNR of the signal received by the master node can be expressed as Furthermore, the maximum transmission range of intra-cluster communication can be provided as where SNR MI,th indicates the minimum SNR that the master node is able to receive and recover the transmitted signal from slave nodes.

Connectivity Analysis
From the problem formulation in Section 3, the probability P C,k,j participates in the connectivity analysis. This section provides the derivations of P C,k,j in direct and multi-hop fashions. As the analysis form of P C,k,j in multi-hop is very complicated, the lower and upper bounds are developed, which can provide the guidelines for the networks design in practice.

Directly Connected
Without loss of generality, it is supposed that the R 3 is a R 3 cube. There is a surface BS and m clusters (m = 1, 2, · · · ) distributed in the R 3 . According to the homogeneous distribution of the clusters, P C,k,j is the probability that the jth cluster falls into the hemisphere with surface BS as the center of the sphere and the radius of R ac,k . We define the area Z 3 where clusters can establish a direct link with surface BS. Therefore, P C,k,j can be derived as where V R 3 = R 3 . However, when considering the boundary effect, the transmission range R ac,k gets larger with the increase of the number of intra-cluster nodes. Based on the values of R ac,k and R, V Z 3 can be calculated as: Case 1: When 0 ≤ R ac,k < R 2 , Z 3 is a hemisphere, as shown in Figure 3a. We have

Case 2: When
2 , the hemisphere partially extends beyond the sides of the cube to form four hemispheric crowns, as shown in Figure 3b. We have

Case 3: When
√ 2R 2 ≤ R ac,k < R, the hemisphere extends beyond the sides of the cube and does not reach the bottom of the cube, as shown in Figure 3c. V Z 3 is calculated by calculus. We have Case 4: When R ≤ R ac,k < √ 5R 2 , the hemisphere extends beyond the sides of the cube and extends beyond the bottom of the cube to form a spherical crown, as shown in Figure 3d. We have

Case 5: When
2 , R ac,k continues to increase, the hemisphere does not completely contain the entire cube, as shown in Figure 3e. We have Case 6: When R ac,k ≥ √ 6R 2 , the cluster at any point in R 3 can be connected to BS, and then V Z 3 = R 3 .
Based on the above analysis, P C,k,j can be calculated as Equation (39). (39)

Lower Bound of Connectivity Probability in Multi-Hop Fashion
Assuming there is a connected link with a relay cluster B between the cluster A and surface BS, there are a nodes in cluster A and b nodes in relay cluster B. Then, the connectivity probability can be calculated as where P a is the probability that the cluster A has a nodes, P b is the probability that the relay B has b nodes, P(A ↔ B) is the probability that the cluster A is connected with relay B, and P(B ↔ surfaceBS) is the probability that the relay B is connected with surface BS. Furthermore, From the above analysis, the lower bound of the connectivity probability can be obtained by assuming that all clusters have the same number of nodes, which means that the maximum communication range of all clusters is the same.Therefore, we can obtain where P C,k,j is the connectivity probability that the cluster with k nodes is connected with surface BS in multi-hop fashion, and it is worth noting that the number of nodes in the cluster as a relay is also k.
To derive the lower bound of connectivity probability of jth cluster with k nodes P C,k,j , we first map the networks on a discrete cube, as shown in Figure 4 The straight line cb connecting the jth cluster and the surface BS is set as the horizontal line. The jth cluster is located at the center of the cube and sides of the cube are parallel to cb. The length of each side is d = √ 6 6 R ac,k , which means that clusters in two adjacent cubes can be connected with each other. The clusters are distributed according to a homogeneous Poisson point process; as a result, they can be randomly divided into different cubes.
If the jth cluster can connect with the surface BS through multi-hop, there could be one or more paths that are distinguished by different relay choices. The set of these paths is denoted as P c = {P 1 c , P 2 c , · · · }, where P 1 c , P 2 c , · · · denote all possible paths. When the cluster density is high, we use the maximum probability that a certain path exists as the lower bound of P c , i.e., max i {P(P i c )}. When the cluster density is low, we calculate the lower bound of the probability that there is at least one open path, i.e., 1 − P(|P c | = 0). The lower bound is P c , where |P c | is the number of the existing paths and |P c | = 0 means that there is no paths between the jth cluster and surface BS. To get the valid lower bound of connectivity probability P(P c ), the larger of the two bounds is selected as the lower bound of P(P c ). Therefore, For the different paths that the jth cluster can connect with surface BS, the shorter path can yield the higher connectivity probability. Therefore, the shortest path of the relay nodes that exist in the cubes for which a straight line cb passes can yield the maximum existing probability. The number of cubes W can be derived as where D is the distance between the jth cluster and the surface BS and a means rounding a to the nearest integer no less than a.
When the cluster density is high, the lower bound of connectivity probability is max i P(P i c ) = P W (There exists at least one cluster in a cube) where q = e −ρ c d 3 is the probability that there is no cluster in a cube. When the cluster density is low, we need to calculate the probability P(|P c | = 0). When the jth cluster cannot connect with the surface BS, there must be a hull C that contains the jth cluster and excludes the surface BS, as shown in Figure 4. The cubes inside the hull contain at least one cluster, while those cubes outside the hull have no cluster. Hence, P(|P c | = 0) can be evaluated by counting the number of hulls and the hulls must be passed by the straight line cb.
Firstly, the number of faces of a hull is denoted as n, the number of hulls that begin at a face and has n faces is denoted as σ(n). Obviously, σ(n) is less than the number of self-avoiding walks beginning at a face, i.e., σ(n) ≤ 8 · 7 n−1 . The total number of such closed hulls is denoted as CN There are relays on one side of the hull and no relays on the other side. Hence, the probability that the jth cluster cannot connect with the surface BS P(|P c | = 0) is Substituting Equations (45) and (47) into Equation (43), we have (48)

Upper Bound of Connectivity Probability in Multi-Hop Fashion
To get the upper bound of connectivity probability, we consider the probability in terms of hops. Therefore, the connectivity probability of the jth cluster is where P 0 is probability that the jth cluster is connected with surface BS directly, P 1 is probability that the jth cluster is connected with surface BS through one relay, and P 2 is probability that the jth cluster is connected with surface BS through two relays. From Equation (49), we can derive that the more paths there are, the greater is the probability of connectivity. The connectivity probability increases with the number of interconnected nodes until there are no isolated nodes. The probability that there are no isolated clusters is the upper bound for the connectivity probability: These clusters are distributed according to a homogeneous Poisson point process with density ρ c , where the probability of the cluster is isolated is Substituting Equation (52) into Equation (51), we have

Simulation Results
To evaluate the accuracy of developed mathematical models, some simulations were carried out. The simulation conditions were set as follows. The R 3 was a 20 km × 20 km × 20 km cube. The clusters were deployed according to a homogeneous Poisson point process with the density of ρ c , while the intra-cluster nodes were distributed according to another homogeneous Poisson point process with the density of ρ n . The transmitting power of each node was 10 dBm, and the underwater noise level was −83 dBm. The MI communication was operated at the frequency of 10 MHz with the bandwidth of 20 kHz. The underwater acoustic communication was performed at the frequency of 10 kHz with the bandwidth of 10 kHz. The frequency of crystal oscillator used for nodes' clock synchronization was 100 kHz. The packet lengths of L 0,z and L z,0 were set to 100 bits. The propagation speeds of MI and sound were 3.33 × 10 7 m/s and 1500 m/s, respectively.

Transmission Range
The maximum transmission ranges of the cluster are plotted in Figure 5. With the increase of SNR threshold, the maximum transmission range decreases subsequently, which is consistent with the analysis. At the same time, connectivity is affected by communication range of the sensor node, which is enlarged by cooperative MIMO compared with single node. We show that the cooperative MIMO using MI-assisted synchronization and pure acoustic-based synchronization both have a great improvement in maximum transmission range compared with a single node in Figure 5. Therefore, the underwater MI-assisted acoustic cooperative MIMO networks can leverage the gains of MIMO in a distributed fashion to enhance connectivity of UWSNs. Meanwhile, MI communication works at ultra-high frequency and the signal transmission delay can be neglected. Thus, the influence of frequency and time estimation errors on its SNR is almost 0, and the maximum transmission range is almost the same as with perfect synchronization. However, in a pure acoustic-based synchronization system, due to the existence of error, the maximum transmission range is obviously less than the perfect synchronization. To explore the transmission range versus the variety of coil radius and SNR threshold, a simulation experiment was conducted, and the results are indicated in Figure 6. As shown in Figure 6, the range of the cluster can be extended by increasing the coil radius. Compared with underwater acoustic communication, MI communication has a smaller transmission range, so it is suitable for short distance inter-cluster communication. At the same time, it is worth noting that there is great difference in the transmission range of the MI communication between sea water and lake water. Highly conductive sea water induces significant Eddy current that incurs very high path loss, which results in the inevitable decrease of the transmission range. However, the transmission range of underwater MI communications can be increased by using the optimal operating frequency and larger coil antennas. Lake water, r=5m Lake water, r=2m Lake water, r=1m Sea water, r=5m Sea water, r=2m Sea water, r=1m Figure 6. The intra-cluster maximum transmission range.

The Directly Connected Probability of the jth Cluster Consisting of k Nodes P C,k,j
To get the directly connected probability of the jth cluster with consisting of k nodes, i.e., P C,k,j , it is necessary to calculate the probability that the actual distance R j between the surface BS and the jth cluster is no more than the maximum transmission range R ac,k , i.e., P C,k,j = P(R j ≤ R ac,k ). To obtain the statistical performance, 10, 000 Monte Carlo simulations were performed. In a certain simulation, an underwater MI-assisted acoustic cooperative MIMO network was generated. If R j ≤ R ac,k , the result was marked as 1; if R j > R ac,k , the result was marked as 0. Finally, the probability P C,k,j was obtained by the 10,000 simulation results. As shown in Figure 7, the simulation results verify the correctness of the theoretical analysis in Equation (39).

Directly Connected Probability of the jth Cluster with Different SNR Threshold
The directly connected probability of the jth cluster with different SNR thresholds is investigated in Figure 8. To evaluate the performances of the proposed approach, 5000 Monte Carlo simulations were carried out in each case. Figure 8 shows that the higher SNR threshold leads to the smaller connectivity probability under the same node density. As SNR th increases from 5 to 20 dB, the transmission range of sensor gets shorter; as a result, fewer sensors participate in the cooperative MIMO, which means the smaller transmission range. Therefore, the probability of connectivity is lower. Note that, when node density is low, there is notable gap between analysis and simulation results. That is because we can gat that the transmission range of a cluster has a great influence on the theoretical analysis of its connectivity from Equation (39). While the cluster communicates with the BS in the form of cooperative MIMO, the number of nodes participating in the cooperative MIMO determines its transmission range. When the connectivity of the jth cluster is theoretically analyzed, the number of nodes may be non-integer, which is calculated by ρ n πR 2 MI . However, the number of nodes in the cluster is determined in simulation. When the node density is low, there are few nodes participating in cooperative MIMO, and the non-integer number of nodes causes notable gap between analysis and simulation results, as shown in Figure 8.

Theoretical Bounds of the jth Cluster
In Figure 9, the theoretical upper and lower bounds of underwater MI-assisted acoustic cooperative MIMO networks are compared with the simulation results with a variety of sensor densities. Each simulated connectivity probability was calculated based on 500 simulation iterations. The SNR threshold of the received signal was set to 20 dB. The connectivity probability increases as the cluster density increases. This result is consistent with the connectivity analysis and provide the guideline for the deployment of the networks according to the practical requirement.
Cluster density/V

Directly Connected Probability of the jth Cluster
The jth cluster's directly connected probabilities in the cases of perfect synchronization, pure acoustic-based synchronization, and MI-assisted synchronization were investigated under the same simulation conditions, and the results are indicated in Figure 10. The connectivity probabilities increase as the node density increases. As the node density in the cluster becomes higher, it demonstrates that more sensors within a cluster participate in the cooperative MIMO, resulting in the long-range transmission. Considering the frequency and time estimation errors in the synchronization process, the maximum transmission range of the cluster is smaller than that of the perfect synchronization. Therefore, for the cluster with the same sensor density, the connectivity probability becomes smaller.

Conclusions
Recently, underwater MI-assisted acoustic cooperative MIMO networks have attracted much attention as they significantly improve the effective cover space and network throughput. However, the complex channel characteristics and the heterogeneous architecture make the connectivity more complex than that of traditional UWSNs. In this paper, we investigate the connectivity of the networks considering the influences of channel characteristics, system parameters, and synchronization error. To provide the theoretical guide for the design of underwater MI-assisted acoustic cooperative MIMO networks, the upper and lower bounds are both derived. Finally, Monte Carlo simulations were performed to address the correctness of the proposed mathematical models.

Conflicts of Interest:
The authors declare no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.