RSS-Based Target Localization in Underwater Acoustic Sensor Networks via Convex Relaxation

The received signal strength (RSS) based target localization problem in underwater acoustic wireless sensor networks (UWSNs) is considered. Two cases with respect to target transmit power are considered. For the first case, under the assumption that the reference of the target transmit power is known, we derive a novel weighted least squares (WLS) estimator by using an approximation to the RSS expressions, and then transform the originally non-convex problem into a mixed semi-definite programming/second-order cone programming (SD/SOCP) problem for reaching an efficient solution. For the second case, there is no knowledge on the target transmit power, and we treat the reference power as an additional unknown parameter. In this case, we formulate a WLS estimator by using a further approximation, and present an iterative ML and mixed SD/SOCP algorithm for solving the derived WLS problem. For both cases, we also derive the closed form expressions of the Cramer–Rao Lower Bounds (CRLBs) on root mean square error (RMSE). Computer simulation results show the superior performance of the proposed methods over the existing ones in the underwater acoustic environment.


Background
The target localization techniques based on an underwater acoustic wireless sensor networks (UWSNs) have recently attracted much attention due to their wide applications in many areas, including data collection, pollution monitoring, offshore exploration, disaster prevention, target tracking, and assisted navigation [1]. In such UWSNs, the sensor nodes are usually classified into anchor nodes and target nodes, where the locations of anchor nodes are known, and the locations of target nodes are unknown and need to be determined [2]. The target localization in UWSNs is often classified as four types based on the measurements: time-of-arrival (ToA), time-difference-of-arrival (TDoA), angle-of-arrival (AoA), and received signal strength (RSS). However, the measurements based on the first three types of target localization methods require complicated timing, synchronization, and line of sight, thus are often difficult to be obtained in underwater acoustic environments [3,4].
Currently, most works are focused on two typical RSS-based target localization methods: the maximum likelihood (ML) [5] and least squares (LS) [6]. While the ML estimator can achieve high accuracy and approximate the Cramer-Rao lower bound (CRLB), its performance depends highly on

Contributions
In this paper, we propose a new approach to the RSS-based underwater acoustic localization problem based on the convex relaxation technique in UWSNs. Unlike the SDP method in [13] and WLS methods in [14], we reformulate the original underwater acoustic path loss measurement model as the pseudolinear equation, derive a minimum optimization problem, and then directly transform this optimization problem into a mixed SD/SOCP localization problem for estimating an accurate location of target. Although the WLS methods in [14] can simply solve the underwater RSS localization problem with lower complexity, its accuracy cannot be perfect. This paper is to improve the underwater acoustic RSS localization accuracy. Two cases are considered with respect to the target transmit power. For the first case, where the reference power is assumed to be known as a measure of the target transmit power, we derive a novel WLS estimator by using an approximation to the RSS expressions, and then transform the originally non-convex problem into a mixed SD/SOCP problem for reaching an efficient solution. For the second case, there is no knowledge on the target transmit power, and we treat the reference power as an additional parameter to estimate. In this case, we formulate a WLS estimator by using a further approximation, and present an iterative ML and mixed SD/SOCP algorithm for solving the WLS problem. For both cases, we also derive the Cramer-Rao Lower Rounds (CRLBs) on root mean square error (RMSE) to evaluate the performance of the two proposed WLS estimators by using computer simulations.
The main contributions of our work are summarized as follows: (1) Based on the convex relaxation, we propose a new approach to the RSS-based underwater acoustic localization problem in UWSNs. (2) The correlation between the target node location and transmit power can be removed by using the separation constant technique which will be helpful for solving the proposed localization method.
(3) In the unknown target transmit power case, an auxiliary constant is used to establish the semi-definite constraints., and then an iterative ML and mixed SD/SOCP algorithm is proposed to estimate both the target location and the transmit power. (4) The closed form CRLBs of the proposed scheme are derived.
The following notations are adopted throughout the paper. Bold face lower case letters and bold face upper case letters denote vectors and matrices, respectively. R n denotes the set of n-dimensional real column vectors. l i denotes the ith entry of the vector l. In addition, · denotes the 2 -norm.
The remainder of the paper is organized as follows. In Section 2, we first give the ray trajectory in an underwater medium, and then the RSS models for underwater acoustic node localization are discussed. In Section 3, we present the proposed localization method. In Section 4, the Cramer-Rao Lower Rounds are derived for the proposed method. Computer simulation results are presented in Section 5. Finally, the main conclusions are drawn in Section 6.

System Model
In this section, we first given the arc length of tracing a ray between a target node (denoted as T) and the anchor nodes (denoted as A) in an underwater environment [15,16]. Then, the RSS-based localization problem in underwater acoustics environment is formulated.

Ray Trajectory
We consider a target localization problem in a 3D underwater acoustic environment. The main considered in this paper are the stratification effect of water medium. Aiming at the stratification effect of underwater environment, in this paper, we propose a new approach to the RSS-based underwater acoustic localization problem based on the convex relaxation technique in UWSNs. In this approach, the underwater sound speed profile (SSP) is assumed only linearly depth dependent and can be approximated as [15][16][17] v where z denotes the depth, b indicates the sound speed at the water surface, and a is the steepness of SSP depending on the environment of the stratification effect of water medium. Let [x A , y A , z A ] and [x T , y T , z T ] denote [r A , z A ] and [r T , z T ], respectively, in 3D space. In Figure 1, we can seen that the UWSNs is 3D, while the ray equations are established on a 2D plane. This is because of the cylindrical symmetry around the z-axis when the z-axis crosses anchor node A. In this case, we can transfer the target localization problem to the 2D plane, which includes both nodes and z-axis. When the z-axis does not cross anchor node A, we still consider it to be a general 3D space. Equation (1) reveals that the SSP is only a linear function of the depth, but the SSP and horizontal distance r are nonlinear relations. Ray tracing method is guided by Snell's law [17]: where θ A , θ T are the ray angles at anchor node and the target node locations, respectively; z A , z T are the depth of the anchor node and target node, respectively; and c is constant along a ray traveling between the two nodes. Then, the arc length l of acoustic propagation path between the anchor node and target node has the following form [15,16] Figure 1. The stratification effect description of a ray between a target node and an anchor node though two adjacent layers.
In the next section, we suppose that the speed in each layer is assumed to be a constant. In this case, the arc length l is assumed to be approximated as a straight line, which is used in the following discussion.

Underwater Acoustic RSS Model
The main characteristics of underwater acoustic channel are frequency-dependent attenuation, time-varying multipath propagation, and low speed of sound. Up to now, the underwater acoustic channel is considered the most difficult communication media. Underwater acoustic transmission loss experienced by a narrow band acoustic signal centered at frequency f (kHz) traveling over distance l is given by Urick propagation model, and then it can be transformed into RSS model. In the following, we describe the signal model considered in this paper.
Consider an UWSNs with N anchor nodes and one target node. For a k-dimensional (k = 2 or 3) localization scenario, we suppose that the target node is located at an unknown position x ∈ R k , and the locations of anchor nodes are known as s i ∈ R k , where i ∈ N {1, 2, . . . , N}. Under a centralized processing mode, the anchor nodes convey their RSS measurements to the central processor for estimating x.
By the underwater acoustic propagation log-normal shadowing model, the RSS in dBm at each anchor node i is given by [13,14,18] where P 0 is the reference power, d 0 is the reference distance, and l i is the ray length between the target node and ith anchor node, which can be calculated by Equation (3). For the same layer of acoustic speed, the arc length between the target node and ith anchor node is approximated to l i = x − s i , γ is the path-loss exponent between 2 and 4 depending upon the propagation environment, α f is the absorption coefficient, and n i is the log-normal shadowing effects [19], which results from degradation of the acoustic intensity caused by multipath propagation, refraction, diffraction, and scattering of sound. Note that α f can be obtained in dB/km as a function of signal frequency f by Thorp's formula [20,21] Let P = (P 1 , P 2 , . . . , P N ) T and n = (n 1 , n 2 , . . . , n N ) T be the collections of the RSS measurements and the associated log-normal shadowing effects, respectively. In most studies, the log-normal shadowing is generally modeled as independent and identical distribution (i.i.d.)zero-mean Gaussian random variable. Unfortunately, we do not have such a guarantee about n i in underwater acoustic environment. In this paper, we suppose that the N × N covariance matrix Q of n is symmetric with elements given by where 0 ≤ ρ < 1 is the common correlation coefficient between any two different entries. Given the target location x, the conditional probability density function (PDF) of P is given by [18,22] whereP(x) = (P 1 (x),P 2 (x), . . . ,P N (x)) T , and for all i ∈ N, Then, the ML estimator of x can be formulated as [22] x = arg min Clearly, the ML estimation, whether in Equation (9), yields a non-convex problem. In the next section, we propose a mixed SD/SOCP method for efficiently solving the localization problem by using appropriate approximations for the RSS expressions.

The WLS Estimation by Approximation
Suppose that the noises n i are quite small when compared with the RSS measurements P i . Let l 0 = 1 m without loss of generality. From Equation (4), we have the following approximation expression be the exponent of the second factor term on the right hand side of Equation (10). It is shown in [13] that 0 < |u| 1 since α f is relatively minor in UWSNs, as shown in Equation (5). Therefore, it is reasonable to approximate the factor term 10 u to be its first-order Taylor expansion at the point u = 0, i.e., where the higher-order terms are omitted. Substituting Equation (11) into Equation (10) yields a pseudolinear approximation where η = 10 10γ η. Based on Equation (12), the WLS estimation for the target location x can be formulated as where {η, λ i } N i=1 are dependent upon both the reference power P 0 and the RSS parameters Clearly, the proposed WLS estimation is non-convex.

Case of Known Target Transmit Power
In this case, we suppose that the transmit power of the target node is known in terms of the reference power P 0 . It follows that the parameters {η, λ i } N i=1 in Equation (13) can be determined. For the sake of presentation, by introducing auxiliary variables l, l i , the WLS estimation (Equation (13)) can be equivalently written as where Then, by further introducing auxiliary variables {L, r}, and applying the SDP technique, the WLS estimation (Equation (14)) can be rewritten in the constricted form as Furthermore, according to the Cauchy-Schwarz inequality, we have for all i, j ∈ N, i > j.
Adding Equation (17) into Equation (15) as a constraint, dropping rank-1 constraint, and applying the SOCP technique, the Equation (15) can be transformed to be where the constraint in Equation (18a) for all i, j ∈ N, i > j.
In the sequel, the mixed SD/SOCP problem in Equation (19) is referred to as "SD/SOCP-K", which can be solved efficiently by using the standard interior method [24], for example, the MATLAB CVX package [25]. Since there exist K 1 = 5 optimization variables in Equation (19), the worst-case complexity has order of O((N + K 1 ) 3.5 ), which is less than that for solving a SDP problem [13].

Case of Unknown Target Transmit Power
In practice, it is impractical to acquire the knowledge on the transmit power of the target node in UWSNs. It is thus worth investigating the localization problem with unknown P 0 . To cope with the difficulty caused by the unknown parameters, we explore a further approximation in addition to Equation (12).
By the definitions of β i and η, using the separation constant technique, the approximation in Equation (12) can be reorganized as Let v = α f ln 10 10γ x − s i . Following the explanations beneath Equation (10), it also holds that 0 < |v| 1. It is thus reasonable to approximate the factor in parentheses on the left hand side by using the first-order Taylor expansion about v = 0 as follows Substituting Equation (21) into Equation (20) yields a simplified approximation Hence, instead of solving Equation (13), we obtain a WLS formulation for estimating both x and P 0 in the form where again η = 10 P 0 10γ , and β i = 10 10γ . Clearly, the WLS problem in Equation (23) remains non-convex. To cope with the objective function, we introduce an auxiliary vector g = (g 1 , g 2 , . . . , g N , η + c 0 ) T , where g i = x − s i , g N+1 = η + c 0 , c 0 is an arbitrary constant, such that the objective function is written as ( Bg − V ) T Q −1 ( Bg − V ). Then, we transform the WLS problem in Equation (23) into the following constricted form G g g T 1 0, rank{G} = 1, (25f) Furthermore, according to the Cauchy-Schwarz inequality, we also have Adding Equation (27) into Equation (25) as a constraint, and dropping rank-1 constraint, by the reformulations above, the non-convex optimization problem is casted into a mixed SD/SOCP with the optimization variables {x, η, r, G, g} as follows min x,η tr C G g g T 1 .
for all i, j ∈ N, i > j.
The mixed SD/SOCP problem in Equation (28) is referred to as "SD/SOCP-U", which can be solved by calling the MATLAB CVX package at a computational complexity of order O((N + K 2 ) 3.5 ), where K 2 = 6.
In Algorithm 1, we also propose a hybrid ML-SD/SOCP method by iteratively estimating P 0 and calling the standard CVX package for solving Equation (28), where the maximum iteration number N max can be determined empirically with aid of numerical results on the convergence.

CRLB Analysis
For the target location x, the WLS estimatorx in either Equation (19) or Equation (28) is an unbiased estimator, i.e., E(x) = x, as they are essentially based on Equation (9). Then, the covariance matrix ofx is subject to the CRLB as VAR(x) F −1 , where F is the Fisher information matrix (FIM).
To evaluate the accuracy performance of the unbiased estimation, the RMSE is defined as wherex i is the estimate of the randomly generated target location x i in the ith simulation, and M is the number of independent simulation rounds. Accordingly, we define the CRLB on RMSE by computing the root trace of F −1 . As proved in the Appendix A, the CRLBs on RMSE for both SD/SOCP methods are given, respectively, by where G = ∂P(x) ∂x , H = (G, 1 N ) N×3 , and which is derived from Equation (8).

Simulation Results
In this section, computer simulation results are provided to evaluate the performance of the proposed mixed SD/SOCP methods, i.e., SD/SOCP-K and SD/SOCP-U. The RSS values were generated according to Equation (4) under the underwater acoustic propagation model, where l 0 = 1 m, P 0 = −40 dBm, and γ = 2, which is applicable to the scenario of underwater acoustic spherical spreading. We assumed that the depth of the target node could be obtained through pressure sensors, and, therefore, only the x and y coordinates of target location were estimated, which reduced the computational complexity of the proposed methods. The same hypothesis is also seen in [26]. We also assumed that only one layer of SSP was considered in the simulations. Then, the anchor and target nodes were randomly located within a square region of 100 × 100 m 2 , and the correlation coefficient was set to be ρ = 0.8. A total of Mc = 3000 Monte Carlo simulations were carried out. Here, target location x was a random variable of Monte Carlo estimation. The estimation of random variables x (target location) in each Monte Carlo simulation was independent of each other. For reader's convenience, the simulation parameters are listed in Table 1. Note that SD/SOCP-K was simulated by directly calling the CVX package [25,27], whereas SD/SOCP-U was simulated by using Algorithm 1. For comparison, simulation results on the SDP method in [13] and the WLS methods (including WLS-K and WLS-U) in [14] are provided. For convenience, the discussed estimators are listed in Table 2. Table 2. Summary of the compared methods.

SDP SDP method in [13] WLS-K
WLS-K method in [14] for known transmit power WLS-U WLS-K method in [14] for unknown transmit power SD/SOCP-K Proposed the new SD/SOCP-K method for known transmit power (Equation (19)) SD/SOCP-U Proposed the new SD/SOCP-K method for unknown transmit power (Equation (28)   Note that the acoustic signal attenuates dramatically in underwater acoustic channels, which depends upon both transmission range and frequency. In addition, underwater acoustic signal frequencies are available from tens of hertz to hundreds of kilohertz due to the severe frequency-dependent attenuation [28][29][30][31]. In simulations, we chose several frequency components for f within the range from 34 kHz to 454 kHz, and the associated absorption coefficient α f varies from 0.01 dB/m to 0.1 dB/m by Equation (5). We considered four localization scenarios to evaluate the performance of the proposed methods in underwater acoustic environment. Figure 3 shows the simulation results on the effects of the noise standard derivation σ on the RMSE for σ varying from 1 dB to 5 dB when P 0 is known (Figure 3a) and unknown (Figure 3b), where the number of anchor nodes is fixed at N = 10, and the absorption coefficient is fixed at α f = 0.01 dB/m. By the simulation setups, we had u = −0.05 and v = 0.1151 for the approximations in Equations (11) and (21), respectively. In these figures, it is naturally observed that the RMSE increased with the noise level in all methods. It was also observed that both the proposed mixed SD/SOCP methods yielded smaller RMSE values than the SDP method and WLS method. Furthermore, Figure 3 shows the superior RMSE performance of the proposed methods SD/SOCP-K and SD/SOCP-U when σ is small, while the gap between the SD/SOCP methods and WLS methods become smaller when σ is large.   Figure 4 shows the simulation results on the effects of the number of anchor nodes N on the RMSE for 8 ≤ N ≤ 22, σ = 4 dB, and α f = 0.01 dB/m when P 0 is known (Figure 4a) and unknown (Figure 4b). In this case, it remained that u = −0.05 and v = 0.1151. It was observed that the RMSE appears a decreasing function of N for both the proposed SD/SOCP methods and WLS methods. In contrast, the RMSE of the SDP method approximately increased with N, and it was considerably larger than that of the proposed SD/SOCP methods over the entire range of considered N values. The main reason for this lies in adopting more approximations and replacing the 2 norm with the 1 to solve the original localization problem of the SDP method. On the other hand, as N increased, the gap between SD/SOCP and WLS becomes smaller, while the gap between the achieved RMSE and the CRLB becomes larger. Figure 4 shows the superior RMSE performance of the proposed methods for all chosen N.   To further verify the effectiveness of the proposed methods, it was compared with other discussed methods. We also considered the localization scenario of three-dimensional underwater acoustic wireless sensor networks, where the anchor and target nodes were randomly located within a cube region of 100 m×100 m×100 m. Figure 2b shows an example of a 3D network configuration and estimation accuracy results for the proposed methods. Figure 6 shows the simulation results on the effects of the noise standard derivation σ on the RMSE for σ varying from 1 dB to 5 dB in three-dimensional underwater acoustic wireless sensor networks, where the number of anchor nodes is also fixed at N = 10, and the absorption coefficient was fixed at α f = 0.01 dB/m. By the simulation setups, we had u = −0.05 and v = 0.1151 for the approximations in Equations (11) and (21), respectively. As in Figure 3, in these figures, it can be seen that the proposed methods also showed excellent performance in three-dimensional scenario. It may be concluded that, although both SD/SOCP-K and SD/SOCP-U are derived under the assumptions of small values for both noise and absorption coefficient, the simulations showed that excellent performance could still be achieved, even for high noise level and large absorption coefficient.

Conclusions
A mixed SD/SOCP approach has been proposed for the RSS-based target localization problem in UWSNs with and without the knowledge on the target transmit power. In the case of known transmit power, we treat the reference power as a constant and derive a novel non-convex WLS estimator based on appropriate approximations. Then, the WLS estimation problem is transformed into a mixed SD/SOCP for reaching efficiently an optimized solution. In the case of unknown transmit power, we treat the reference power as an additional unknown parameter and propose a combined ML and SD/SOCP iterative algorithm to jointly estimate both the reference power and the target location. In this paper, we only consider a simplified model, and the depth is assumed to be known in advance. Furthermore, Only simulation results are provided to verify the performance of the proposed methods. We hope to extend the present work in a general model, and for the results to be verified by real data in our future work.

Conflicts of Interest:
The authors declare no conflict of interest.