I Introduction
The recent demand on high quality of service communications for commercial, scientific and military applications of underwater exploration necessitates a high data rate, low latency, and longrange underwater networking solutions [1]. Fulfilling these demands is a formidable challenge for most of the electromagnetic frequencies due to the highly attenuating aquatic medium. Therefore, acoustic systems have received considerable attention in the past decades. As a result of the surfaceinduced pulse spread and frequencydependent attenuation, acoustic communication data rates are restrained to around tens of kbps for ranges of a kilometer, and less than kbps for longer ranges [2]. Moreover, the low propagation speed of acoustic waves (1500 m/s) induces a high latency [3], especially for longrange applications where realtime communication and synchronization are challenging.
Alternatively, underwater optical wireless communication (UOWC) has the virtues of supporting high speed connections in the order of Gbps [4, 5], providing low latency as a result of the high propagation speed of light in the aquatic medium ( m/s) [6], and enhancing security thanks to the pointtopoint links [7]. Nevertheless, the reachable UOWC range is delimited by severe underwater channel impairments of absorption and scattering. The absorption effect refers to the energy dissipation of photons while being converted into other forms (e.g., heat, chemical, etc.) along the propagation path. The scattering is regarded as the deflection of the photons from its default propagation path, which is caused either by water particles of size comparable to the carrier wavelength (i.e., diffraction) or by constituents with different refraction index (i.e., refraction). Therefore, the relation between absorption and scattering primarily characterizes the fundamental tradeoff between range and beam divergence angle. That is, a collimated light beam can reach remote receivers within a tight span whereas a wider light beam can communicate with nearby nodes in a broader span. Accordingly, the rangebeamdwidth tradeoff primarily determines the connectivity and outage in underwater optical wireless networks (UOWNs) [8].
The directed nature of UOWC also necessitates efficient pointingacquisitioningandtracking (PAT) mechanisms, especially for the collimated light beams propagating over a long range. Hence, accurate node location information has a critical importance for effective PAT mechanisms [9]
. However, the estimated node locations may not refer to actual node positions either because of the localization errors or the random movements caused by the hostile underwater environment. These have negative impacts on the link reliability and performance as a result of the poor PAT efficiency. A potential solution for location uncertainty would be employing adaptive optics to provide a degree of robustness by adjusting the bandwidth to cover a wider spatial area around the estimated receiver location.
Limited communication range and hostile underwater channel conditions also entail communicating over multiple hops in order to improve both connectivity and performance of UOWNs. In this respect, there is a dire need for analyzing the endtoend (E2E) performance of multihop UOWCs under prominent relaying techniques. Moreover, multihop UOWCs require effective routing protocols that can tackle the rangebeamwidth tradeoff and adversity of underwater environment. It is therefore of utmost importance to develop robust relaying and routing techniques which manipulate rangebeamwidth tradeoff by means of adaptive optics, which is main focus of this paper.
Ia Related Works
Recent efforts on UOWC can be exemplified as follows: Arnon modeled three types of UOWC links: lineofsight (LoS), modulating retroreflector, and nonLoS (NLoS) [10]. In [11], authors model and evaluate underwater visible light communication (UVLC) vertical links by dividing the UWON into layers and taking account of the inhomogeneous characteristics (e.g., refractive index, attenuation profiles, etc.) of each layer. Experimental evaluation of orthogonal frequency division multiplexing based UVLC is conducted in [12] and performance of automatic repeat request (ARQ) based UVLC system is analyzed in [13]. Considering Poisson point process based spatial distribution of nodes, Vavoulas et. al. analyzed the connectivity of UOWNs under the assumption of omnidirectional communications [14]. To take the directivity of UOWC into account, the connectivity analysis of UOWNs is readdressed in [8] which is further extended to investigate impacts of connectivity on the localization performance [15].
Akhoundi et. al. introduced and investigated a potential adaptation of cellular code division multiple access (CDMA) for UOWNs [16, 17]. In [18], authors characterized the performance of relayassisted underwater optical CDMA system where multihop communication is realized by chip detectandforward method. Assuming identical error performance at each hop, Jamali et. al. consider the performance analysis of multihop UOWC using DF relaying [19]. Since the commonly adopted BeerLambert UOWC channel model assumes that the scattered photons are totally lost, authors of [20] modify this renown model to consider the received scattered photons by a single parameter which is a function of transceiver parameters and water type. Thereafter, they consider a dualhop UVLC system and determine optimal relay placement to minimize the bit error rate (BER) for both DF and AF relaying.
IB Main Contributions
The main contributions of this paper can be summarized as follows:

To model the displacement of the actual locations from the estimates, we consider an uncertainty disk whose radius is a design parameter that depends on the localization accuracy and/or mobility of nodes. Assuming nodes are equipped with adaptive optics, robust UOWC links are provisioned by deriving the divergence angles for scenarios with and without a PAT mechanism. Numerical results show that location accuracy and effective PAT mechanism are crucial to reach a desirable performance.

We investigate and compare two prominent relaying techniques; decode & forward (DF) and optical amplify & forward (AF), each of which is analyzed for important E2E performance metrics. For both relaying methods, closed form expressions are derived for various performance metrics, e.g., E2ERate, E2E bit error rate (BER), transmission power, amplifier gain, etc. For a given path, we also formulate and solve minimum total power consumption problems under both DF and AF relaying schemes.

In light of the adaptive divergence angles and E2E performance analysis, both centralized and distributed routing protocols are developed based on the availability of global or local network state information, respectively. The centralized techniques employ variations of shortest path algorithms which are tailored to achieve various objectives such as maximum achievable E2ERate, minimum E2EBER, and minimum total power consumption. By manipulating the rangebeamwidth tradeoff, we propose Light Path Routing (LiPaR); a distributed routing protocol that combines traveled distance progress and link reliability. Obtained results show that LiPaR can provide a superior performance compared to centralized schemes without a PAT mechanism.
IC Paper Organization
The remainder of the paper is organized as follows: Section II introduces the system model. Section III presents the uncertainty disk model and derives the adaptive divergence angles. Section IV analyzes and compares DF and AF relaying techniques. Section V develops the centralized and distributed routing protocols. Section VI presents the numerical results. Finally, Section VII concludes the paper with a few remarks.
Ii System Model
We consider a twodimensional underwater optical wireless network (UOWN) which consists of a surface stations/sinks (SS) and nodes/sensors as shown in Fig. 1. The SSs are responsible for disseminating the data collected from sensors to mobile or onshore sinks. Nodes are equipped with two optical transceivers to enable bidirectional connections between sensors and SSs [c.f. Fig. 2]. Transmitters are assumed to be capable of adapting their beam divergence angle within a certain range, i.e., . The location estimates of the SS and node are denoted by and respectively. Likewise, actual locations of the SS and node are denoted by and respectively. In addition to , node is aware of .
According to the Beer’s law, aquatic medium can be characterized for wavelength as a combination of absorption and scattering effects, i.e., where , and are absorption, scattering and extinction coefficients respectively. These coefficients vary with water types (e.g., clear, harbor, turbid. etc.) and water depths (e.g., shallow, deep, etc.). The propagation loss factor of the LoS channel between node () and node () is defined as follows
(1) 
where is the perpendicular distance and
is the angle between the receiver plane and the transmitterreceiver trajectory, i.e., the pointing vector of the transmitter. On the other hand, geometric loss of the LoS channel is given as
[21](2) 
where is the receiver aperture area of , is the halfbeamwidth divergence angle of ^{1}^{1}1 is measured at the point where the light intensity drops to of its peak., and is the concentrator gain, that is defined as
(3) 
is the angle of incidence w.r.t. the receiver axis, is the concentrator fieldofview (FoV)^{2}^{2}2 can be and down to for hemisphere and parabolic concentrators, respectively., and is the internal refractive index.
As shown in Fig. 2, UOWC channel characteristics pose a fundamental tradeoff between communication range and divergence angle (i.e., beamwidth). That is, a wide divergence angle can reach nearby nodes whereas a narrow divergence angle is able to reach distant destinations. Therefore, manipulation of this tradeoff has significant impacts on several performance metrics such as the degree of connectivity, distance progress for the next hop, routing efficiency, and E2E performance. For instance, a wide divergence angle provides more path diversity at the expense of a higher number of hops which increases energy consumption and E2E delay. On the other hand, a narrow divergence angle can deliver a desirable performance over long distances, which requires precise and agile PAT mechanisms to keep the transceivers aligned over long ranges. Throughout the paper, this tradeoff is to be discussed in the E2E performance analysis of different relaying techniques and the design of the routing protocols.
Iii Adaptive Beamwidth and PAT Mechanism
For a reliable and wellconnected transmission, the inherent directivity of UOWC requires perfectly aligned pointtopoint links between the transceivers. First and foremost information needed for a proper alignment is node coordinates that can be obtained either by pure optical [22] or hybrid acousticoptical localization algorithms [23]. Nonetheless, localization errors introduce uncertainty over the location estimates, which may be further exacerbated by the random movements of nodes due to the hostile underwater conditions. Thus, the accuracy of location estimates plays an essential role in the E2E performance of the multihop communication. One solution to this problem is adaptively changing the properties of the light beam. Indeed, misalignment caused by the location uncertainty can be alleviated by adapting the divergence angle to cover a larger spatial region at the cost of degraded performance. In the remainder, we therefore probe an initial investigation into the robust divergence angles and its impact on the E2E performance by considering the following three different cases.
Iiia Divergence Angles for PAT with Perfect Location Estimates
In the ideal case, estimated and actual locations are the same (i.e., ), the pointing vector is aligned to , and divergence angle is adjusted to exactly cover the node frame, which is shown by plaid sectors in Fig. 3.b. Accordingly, the full divergence angle of the perfect PAT scenario is given by
(4) 
which follows from the law of sines and the fact that radius is perpendicular to the tangent points.
IiiB Divergence Angles for PAT under Location Uncertainty
The location uncertainty can be caused by estimation errors of localization method and/or the random movements of the transceivers, which results in a radial displacement of the actual location from the estimated location. As shown in Fig. 3.a, we model node locations by an uncertainty disk centered at the location estimates and has a radius of where is the radius of the node frame and is the uncertainty metric. Under the location uncertainty, the divergence angle must be calculated based on the worst case scenarios where the transmitter is located at one of the tangents on the uncertainty disk [c.f. Fig. 3.c]. In this case, the divergence angle is twice of the perfect angle covering the uncertainty disk, i.e,
(5) 
which can be defined as a robust divergence angle as it assures the worst case beamwidth to establish a link between the transceivers. Noting that this is the first step of the PAT procedure, i.e., pointing, acquisitioning and tracking are the next steps to keep nodes aligned via feedbacks, which results in Fig. 3.b. Although developing a PAT mechanism is out of our scope, our purpose is to show its critical role in the multihop UOWN performance.
IiiC Divergence Angles under the Absence of PAT
In the absence of the PAT mechanism, we assume that nodes keep their body frame directed towards the location of the closest SS to align transmitter’s pointing vector with the SS receiver. Although pointing to a certain location is not possible, nodes can still reach nodes within a certain angle around the pointing vector by adjusting the beamwidth.
For the sake of a clear presentation, let us first consider the location uncertainty in the receiver node only [c.f. Fig. 4], where the origin of the local Cartesian coordinate system is set to the estimated transmitter location, i.e., . In this local coordinate system, the locations of and the destination SS (SS) are given by and , respectively. Hence, the angles between the xaxis and vectors pointing () and the SS () are given by
(6) 
respectively. Accordingly, the divergence angle that is centered around covering the uncertainty disk is given by (7) where is the signum function.
(7) 
Fig. 4 illustrates the half divergence angle of the first and second cases in (7) by green/red and blue colors, respectively. The last case is the scenario where passes through the uncertainty disk. Similar to the previous subsection, we consider the wost case scenario when there is uncertainty about both receiver and transmitter locations. Therefore, we consider two additional local coordinates originated at tangent locations [c.f. Fig. 5] which are given by
(8)  
(9) 
Noting that is the same in all local coordinates, for () and () should be obtained as in (6). By substituting and into (7), one can obtain divergence angles for transmitters located at tangent points () and (), respectively. In case of the location uncertainty at both sides, the divergence angle is then given by
(10) 
Iv Performance Analysis of Relaying Techniques
In this section, we consider the E2E performance analysis of potential relaying techniques for an arbitrary multihop path from a source node to one of the sink nodes. Such a path is defined as an ordered set of nodes, i.e., where is the number of hops and the first (last) element represents the source (sink).
Iva Decode & Forward Relaying
In DF transmission, the received optical signal at each hop is converted into electrical signal, then decoded, and finally reencoded before retransmission for the next hop. Thus, the DF relaying requires high speed data converters, decoders, and encoders to sustain the achievable data rates in the order of Gbps. For an arbitrary path, received power at node is given as
(11) 
where is the average optical transmitter power of previous node, is the transmitter efficiency of node , is the receiver efficiency of node , and is the composite channel gain.
The most common detection technique for OWC is intensitymodulation/directdetection (IM/DD) with onoff keying (OOK). Assuming that photon arrivals follows a Poisson process, photon arrival rate of node is given as [24]
(12) 
where is the observed power at node , is the detector counting efficiency, is the transmission rate for the hop , is pulse duration, is Planck’s constant, and is the underwater speed of light. Hence, the photon counts when binary ’0’ transmitted is given by where is the total noise power, is the additive noise power due to dark counts and is the background illumination noise power. Similarly, the photon counts when binary ’1’ transmitted is given by
. Assuming a large number of photon reception, the Poisson distribution can be approximated by a Gaussian distribution as per the Central Limit Theorem. For a given data rate,
, the BER of a single hop is given by [10](13) 
where is the complementary error function. For a certain BER threshold, , data rate of the hop is derived from (13) as
(14) 
which is inversely proportional to the BER. Following from (11)(14), minimum transmit power to ensure a predetermined rate and BER is given by
(15) 
where . For a given error and data rate, communication range between two generic nodes () is given by
(16) 
where is the principal branch of product logarithm function, and . We refer interested readers to Appendix A for derivation details.
Before proceeding into the E2E performance analysis, it is necessary to discuss singlehop performance in light of Section III and the above derivations. For a given transmit power and communication range, increasing the divergence angle deteriorate the and , which are inversely proportional as per (13) and (14). On the contrary, communication range decreases by increasing (decreasing) () and decreasing (increasing) the transmit power (divergence angle). Noting the coupled relation of the divergence angle with all these metrics, singlehop performance analysis clearly shows the need for effective PAT mechanisms to obtain a superior E2E performance, which is analyzed next.
Denoting
as the Bernoulli random variable which represent the erroneous decision of node
, total number of incorrect decision made along the path is given as . Assuming ’s are independent but nonidentically distributed, is a PoissonBinomial random variable. Since a transmitted bit is received correctly at the sink if the number of erroneous detection is even, E2E BER is derived as(17) 
where
is the set of odd numbers and
is the set of all subsets of integers that can be selected from .can be expeditiously calculated from polynomial coefficients of the probability generating function of
in where is the cardinality of [25]. Notice that reduces to a Binomial variable if all hops are assumed to be identical as in [19] , which is hardly the case in practice.The achievable E2E rate is determined by the minimum of the data rates along the path, i.e.,
(18) 
which implies that a predetermined E2E data rate requires . Accordingly, formulates the optimization problem which minimizes the total transmission powers while ensuring a target E2E BER rate, , and E2E data rate by setting .
(19) 
where is the BER vector of the path, is the transpose of vector of ones, and denotes the pairwise inequality for vectors. Based on the mild assumption of ensuring each hop has a BER no more than ,
can be shown to be a convex problem using the logconcavity of the PoissonBinomial distribution. We refer interested readers to Appendix
B for a formal convexity analysis of .IvB Optical Amplify & Forward Relaying
Although DF greatly improves the E2E performance by limiting the background noise propagation, it introduces extra power consumption and signal processing delay. Furthermore, synchronization and clock recovery are additional challenges to be tackled in Gbps links. Alternatively, the AF relaying executes opticaltoelectrical (OEO) conversion at each node, which amplifies the received signal electrically and then retransmits the amplified signal for the next hop. The main drawback of the AF transmission is propagation of noise added at each node, which is amplified and accumulated through the path. As a remedy to costly OEO conversion, all optical AF relaying has the advantages of realizing the high speed transmissions without the need for OEO conversion and sophisticated optoelectronic devices. As illustrated in Fig. 6, amplified and transmitted power of alloptical AF scheme can be modeled for intermediate nodes as follows
(20)  
(21) 
where is the local noise at the receiver, is the amplifier gain, is the amplified spontaneous emission (ASE) noise that is modeled as an additive zeromean white Gaussian, is the power spectral density per polarization, is the frequency, is spontaneous emission parameter of the amplifier, and is the amplifier bandwidth [26]. By assuming the independence of signals, channel gains, and noises, amplified and transmitted powers at node can be respectively written in the following generic form
(22) 
(23) 
which follows from the recursion of (20) and (21). Because of the propagating noise along the path, there is no a unique way of ensuring a target E2E performance. However, a practical and tractable means of analyzing the E2E performance is fixing the transmission power at each hop. Noting that the combination of background light and dark current noise is more dominant than the ASE noise [27], the amplifier gain can be calculated as
(24) 
which is obtained by substituting (20) into (21), equalizing (21) to a fixed transmission power , and then solving for . Let us denote the transmitted signaltonoiseratio (SNR) at the source node as and the received SNR in relay node as . By substituting (24) into (IVB) and neglecting the ASE related last term of (IVB), the received SNR at the sink node can be obtained after some algebraic manipulations as
(25) 
where . Accordingly, the E2EBER of optical AF relaying over an arbitrary route is given by
(26) 
which is derived by rewriting (13) for . For a given E2EBER target , the achievable E2E data rate of AF scheme is derived as
(27) 
which follows from (12) and (26). Accordingly, the minimum total transmit power, that ensures a given E2E data and error rate pair, can be obtained by solving the following problem
(28) 
where is the maximum transmit power. Since is a monotonically increasing function of , the optimal transmit power can be easily obtained by standard line search methods. Accordingly, optimal amplifier gains, , that assures minimum E2E energy cost can be calculated by substituting into (24).
V Centralized and Distributed Routing Schemes
Depending on the available network information, routing protocols can be designed either in distributed or centralized fashion. In this sense, distributed solutions are suitable for scenarios where nodes have local information about the network state of its neighborhood. On the other hand, centralized routing relies upon the availability of the global network topology and thus yields a superior E2E performance. However, it is worth mentioning that collecting a global network state information may yield extra communication overhead and energy cost. Therefore, the remainder of this section addresses centralized and distributed UOWN routing techniques.
Va Centralized Routing Schemes
Let us represent the network by a graph where is the set of sensor nodes, is the set of edges, and is the edge weight matrix that can be designed to achieve different routing objectives. In this respect, the remainder of this section considers three primary routing objectives: minimum E2EBER, maximum E2ERate, and minimal total power consumption.
VA1 Minimum E2EBER Routing
The main goal of this routing scheme is to find a path that minimizes the E2EBER, , while ensuring a predetermined E2E data rate, . As previously shown in (18), we satisfy this constraint by setting the data rate of each hop to . By neglecting the fortunate events of corrections at even number of errors, we make the following mild assumption for the DF relaying; the sink can correctly detect a bit ifandonlyif it is successfully detected at each hop.
Thus, we equivalently consider the maximization of the E2E bit success rate (BSR), , instead of minimizing the . Since the shortest path algorithms are only suitable to additive costs, BSR maximization can be transformed from the multiplication into a summation as follows
(29)  
(30) 
where denotes an arbitrary path. By putting the product into an additive form, the maximum BSR route can be calculated by Dijkstra’s shortest path (DSP) algorithm with a time complexity of . Accordingly, edge weights is set to the nonnegative values of where is the BER between and . Once the route is calculated, can be computed by substituting into (17).
In the case of AF relaying, optimal route must account for the noise propagation along the path. Following from (26) and (27), both E2E data rate and BER can be minimized by maximizing the SNR at the destination node, . Based on the same rationale in (29), maximizing the SNR can be transformed into an additive form as follows
(31)  
(32)  
(33) 
which can also be calculated using the DSP algorithm by setting nonnegative edge weights to where is the SNR between and . Once the route is calculated, the can be calculated by substituting into (26).
VA2 Maximum E2E Data Rate Routing
Our target in this routing scheme is to find a path that maximizes E2E data rate while ensuring a predetermined E2EBER, , which is nontrivial since the is coupled with the hop counts as well as . Therefore, we first relax the problem by fixing the BER at each edge and then find the optimal route with the maximum as follows
(34) 
The problem defined in (34) is known as the widest path or bottleneck shortest path problem and can expeditiously be solved by employing the DSP algorithm via modification of the cost updates [28]. Depending upon the calculated path, it is necessary to adjust the BER of links to maximize which can be formulated as follows
(35) 
which maximizes the minimum rate along the path. In order to have an insight into , let us consider a twohop path with channel gains . If these links are set to the same BER, , we have since due to . Thus, is maximized when by manipulating error rates, i.e., . Also notice at the optimal point that the second constraint of is active (i.e., ) since enhancing decrease due to the fundamental tradeoff between rate and error. In light of these, a line search algorithm can find the maximum data rate by minimizing where denotes the absolute value of . Once the maximum is obtained, BER of each link can be calculated by substituting into (13).
VA3 Minimum Power Consumption
The UOWN nodes have limited energy budget that has a significant impact on the network lifetime. Accounting for the monetary cost and engineering hardship of battery replacement, an energy efficient routing technique is essential to minimize the energy consumption of multihop communications. For a given pair of E2E data rate, , and BER constraint, , the objective of this routing scheme is to find the optimal path which provides the minimum total transmission power. Because of the inextricably interwoven relationship between hop counts, and
finding the most energy efficient route falls within the class of mixedinteger nonlinear programming (MINLP) problems which is known to be NPHard. Therefore, we propose fast yet efficient suboptimal solutions for both DF and AF relaying. For the DF relaying, E2ERate constraint can be satisfied by exactly setting the data rate of each edge weight to
. Then, edge weights are calculated using (15) by fixing the BER at each hop. In this way, the minimum power cost path can be obtained by using the DSP algorithm. Once the route is calculated, optimal BER that minimizes the total energy consumption can obtained by solving .Compared to the DF relaying, finding the optimal route for the AF relaying is more complicated. As can be seen from the constraint of (28), the optimal transmission power varies with the number of hops and channel gains along the routing path, which prevents using the DSP algorithm to find the optimal path. Alternatively, we consider selecting the path with minimum among a number of feasible path that satisfies constraints which can be verified by (26) and (27). In this regard, shortest path (KSP) algorithm can be employed as a generalization of the DSP such that it also finds paths with the highest s in the order of [29]. For feasible paths that satisfy and , we solve the problem in (28) and select the path with minimum as the final route.
VB Light Path Routing (LiPaR): A Distributed Routing Scheme
Since node location information is a prerequisite for pointing to establish links, geographic routing schemes can be regarded naturally potential methods for UOWNs. Even though there exists many successful geographic routing protocols designed for omnidirectional RF communications in terrestrial wireless sensor networks, they cannot be directly applied to UOWNs because of the directed nature of optical wireless communications and hostile underwater environment. Accordingly, we propose a distributed routing algorithm, namely Light Path Routing (LiPaR), that only requires the location information of the sink and neighboring nodes. LiPaR does not rely upon a PAT mechanism and assumes that each node directs its pointing vector to the closest sink. Each forwarder node first determines the nodes within its communication range, which is certainly subject to the fundamental rangebeamwidth tradeoff of UOWCs. Therefore, the set of feasible relays within the neighborhood of can be given by
(36) 
which is defined based on a predetermined data rate and BER pair at each hop. In what follows, the forwarder for the next hop is chosen among as follows
(37) 
which maximizes the distance progress by also taking the link quality into account. Notice that varies both with distance and angle between and pointing vector of , . Let us consider two candidate nodes and with equal Euclidean distance but different angles . In case of is chosen as a next forwarder, the transmitter needs to adjust its divergence angle to a tighter beamwidth, which maximizes the average distance progress since is minimized. Following (37), the current forwarder adapts its beamwidth to cover the next one based on (10) and transmit packets after informing the next node about its decision. This procedure is repeated along the path until one of the sink node is reached.
Par.  Value  Par.  Value  Par.  Value 

mrad  rad  
dBm 
Vi Numerical Results
In this section, we provide the performance evaluations using the default parameters listed in Table I which is mainly drawn from [10]. Simulations are conducted on Matlab and presented results are averaged over random realizations.
Via Validation of Adaptive Beamwidth Calculations
Before presenting the performance evaluations, let us illustrate the validation of the beamwidth calculations provided in Section III. In Fig. 7, we demonstrate the adaptive beamwidths between a transmitter located at and a randomly located receiver for three different cases: Case1) PAT with perfect location estimates [c.f. Fig. 6(a)Fig. 6(c)], Case2) PAT under location uncertainty [c.f. Fig. 6(d)Fig. 6(f)], and Case3) Adaptive beamwidths in the absence of PAT [c.f. Fig. 6(g)Fig. 6(i)]. Thanks to the availability of the actual node locations, Case1 is able to tune the divergence angle for a minimum beamwidth to barely cover the receiver node, which naturally delivers the best performance. Compared to the first case, Case2 requires larger beamwidths to compensate the location uncertainty of the transceivers, which degrade the achievable performance. In the last case, the transmitter is directed toward to the sink node located at and required to adjust its beamwidth to cover a receiver located far way from the transmitterreceiver trajectory. Accordingly, performance degrades as the beamwidth increases to establish a link towards a receiver with a higher . In the remainder of this section, these three cases will be indicated by means of superscripts, e.g., DF refers to the DF relaying in Casei, .
ViB E2E Performance Evaluation of DF and AF Relaying
For the E2E performance evaluation of DF and AF relaying schemes, we consider a routing path with hops to reach a sink node far away from the source node. Although the length of hops are the same, the angle between the relays and sourcesink trajectory is random. Fig. 8 demonstrate the impacts of location uncertainty on rate, power, and BSR under different pointing schemes. In Fig. 8, previous discussions are validated as Case always delivers a superior performance compared to Case, . As the location uncertainty increases, the performance degrades for all cases. For , increasing the hop count serve as a remedy to compensate the distance related losses. However, for , decreasing the hop lengths by increasing the hop count triggers the distancebeamwidth tradeoff. That is, divergence angle becomes significant even for a closer node because perpendicular distance and the distance between receiver and pointing vector are close. Therefore, Fig. 8 exhibits a performance enhancement until and degradation after that. This behavior has also an impact on the location uncertainty such that Case3 is not effected from location uncertainty driven beamwidth enlargement because aforementioned effect makes this additional divergence angle increase less significant.
The E2ERate performance is shown in Fig. 7(a) where DF relaying has better rates AF relaying thanks to its ability to detecting and regenerating the bits at each hop. Unlike the first two cases, rates of Case3 monotonically reduces for due to the low rates as a result of larger beamwidths as explained above. In Fig. 7(b), the DF relaying is also shown to be more energy efficient than the AF relaying, which is expected since the noise propagation along the path necessitates a higher transmission power. We should note that schemes requiring transmission powers more than may not be feasible because of the power budget limitations of underwater sensor nodes. Reminding that we merely count for the transmission power, the AF relaying may be far energy efficient than the DF relaying because it does not need energy and computational power hungry conversion and modulation circuitry. Finally, Fig. 7(c) shows that the AF relaying provides a better BSR performance than the DF relaying, this is indeed caused by main difference between DF and AF relaying. Unlike Case3, existence of a PAT mechanism eliminates this difference and provide a desirable performance starting from .
In Fig. 9, we show how different cases of rate, power, and BSR is affected from various water types; pure water (), ocean water (), and coastal water (). Apparently, the rate, power, and BSR performance degrades as the density of the water particulates increases. We also still observe that Case always delivers a superior performance compared to Case, . In Fig. 8(a), we observe that the DF relaying rate surpasses that of the AF relaying in all cases and water types. Interestingly, the benefits of multihop communications become significant as the hostility of the underwater environment increases. Excluding the coastal water, Fig. 8(b) shows that the DF relaying is more energy efficient than the AF relaying. In Fig. 8(c), we observe that coastal water deteriorates the BSR performance in a great extent.
Overall, Fig. 8 and Fig. 9 clearly shows that UOWC urges multihop communication which is established on PAT mechanism along with accurate node locations. Otherwise, E2E performance can be affected from harsh channel conditions and fundamental tradeoff between range and beamwidth. This also necessitates efficient routing protocols which are examined next.
ViC Comparison of the Proposed Protocols
This section compares the performance of the proposed centralized and distributed routing schemes. We consider a network area of
where nodes are uniformly distributed. At each realization, the source node is located at a random location at the sea bed. On the other hand, sink locations are arranged to be centered at the sea surface with equidistant intervals.
Fig. 10 shows the percentage of failure in finding a route between source and one of the sink nodes. In the axis, we denote the routing schemes with their objective and underlying pointing method, e.g., RATE corresponds to maximum rate routing based on Case pointing method, . It is obvious from Fig. 10 that Case1 and Case2 was always able to find a path thanks to their PAT mechanism. However, Case3 suffers from connectivity of the network which is degraded in a great extent due to the rangebeamwidth tradeoff. Nonetheless, increasing the node density helped to reduce failures as a result of increased connectivity. Finally, LiPaR is shown to deliver a good performance in terms of finding a path towards one of the sinks. Reminding that LiPaR cannot guarantee an E2E performance, having Case3 with higher failures is because of the infeasibility to assure an E2E performance.
Fig. 11 shows the average hop count of the calculated paths. A common behavior for all centralized scheme is that number of hops increase with node density. This is expected since having shorter hop lengths provides a better performance in all terms of the performance. However, we observe a slight decrease in average hop count of LiPaR due to its hopbyhop nature. As the node density increases, it is more probable to reach the sink nodes at lower number of hops.
Fig. 12 demonstrates the E2ERate performance of different routing schemes described in Section VA2. Results show that the DF relaying provides a higher E2ERate in all cases. The rate is also monotonically increases with the node density due to better connectivity and higher performance at individual links. Again, LiPaR delivers a better performance than Case3 and it is always fixed to 1 GGbps because candidate nodes are defined based on this fixed rate.
Fig. 13 depicts the E2E total power consumption of different routing schemes described in Section VA2. Results show that the DF relaying requires less total power than the AF relaying. However, the DF relaying can still be more power hungry because of the OEO conversion and signal processing. The one exception to this relation is the Case3 where AF relaying demands less transmit power. Nonetheless, Case3 requires infeasibly high transmit power, which makes LiPaR a better solution as it requires less power and operates in a distributed manner.
The BSR performance of different routing schemes described in Section VA1 is demonstrated in Fig. 14 where all schemes reaches desirable levels except the AF which suffers from the absence of a PAT mechanism and it is not possible to compensate the loss by detecting and regenerating the signals as in DF.
Finally, we have evaluated the performance metrics for various number of sinks, i.e., . Since the centralized schemes already have the entire network information, we did not observe a significant change in their performance for different . However, LiPaR has shown a considerable performance enhancement in terms of failures. The percentage of failures for is recorded as with a slight change in average number hops. Indeed, high number of nodes make it possible to provide a connection opportunity to relays near the surface such that node clusters who can reach these nodes is granted access to the sinks.
Vii Conclusions
Multihop communication is a promising solution to alleviate the short communication range limitation of UOWCs. In particular, a proper operation of UOWNs relies upon the degree of connectivity which can be improved via multihop communications. In this regard, analyzing the E2E performance of multihop UOWC is necessary to gain insight into the UOWNs. In this paper, we accordingly investigated multihop UOWCs underlying two prominent relaying techniques: the DF and AF relaying. Since pointing is a prerequisite to establish each link on a certain path, we accounted for the location uncertainty and the availability of a PAT mechanism. Numerical results show that a PAT mechanism along with adaptive optics can provide desirable E2E performance. Thereafter, we developed centralized and distributed routing schemes. For the centralized routing, we proposed variations of shortest path algorithms to optimize E2E rate, BSR, and total power consumption. Finally, we developed a distributed routing protocol LiPaR which combines traveled distance progress and link reliability and manipulates the rangebeamwidth tradeoff to find its path. LiPaR is shown to provide better performance than the centralized scheme without a PAT mechanism.
Appendix A Derivations of Communication Ranges
Range expressions of LoS can be derived by substituting (11) into (14). After some algebraic manipulations, (11)(14) can be put in the form of which has the following root . Notice that this root () is a solution for the perpendicular distance, hence, the communication range can be obtained from the Euclidean distance where is the angle between the pointing vector and transmitterreceiver trajectory.
Appendix B Convexity Analysis
Let us start with the convexity analysis of the objective function. Since nonnegative weighted sum of convex functions is convex, it is sufficient to prove the convexity of each term, . is a function of and . By omitting the hop indices, the second derivative test for these terms can be given as
which are always positive due to the assumption of . Noting that coefficients in (15) are also nonnegative, each term of the objective is a convex function. The mild assumption of is also necessary to assure the logconcavity of the PoissonBinomial distribution. Therefore, we take the logarithm of both sides of the first constraint to put the problem in a convex form. We refer interested readers to [30, 31, 32] for a more indepth convexity analysis of Binomial and PoissonBinomial distribution.
References
 [1] A. Celik, N. Saeed, T. Y. AlNaffouri, and M.S. Alouini, “Modeling and performance analysis of multihop underwater optical wireless sensor networks,” in 2018 IEEE Wireless Communications and Networking Conference (WCNC), Apr. 2018, pp. 1–6.
 [2] N. Saeed, A. Celik, T. Y. AlNaffouri, and M.S. Alouini, “Underwater optical wireless communications, networking, and localization: A survey,” arXiv preprint arXiv:1803.02442, 2018.
 [3] I. F. Akyildiz, D. Pompili, and T. Melodia, “Underwater acoustic sensor networks: research challenges,” Ad Hoc Networks, vol. 3, no. 3, pp. 257 – 279, 2005.
 [4] H. M. Oubei, C. Li, K.H. Park, T. K. Ng, M.S. Alouini, and B. S. Ooi, “2.3 Gbit/s underwater wireless optical communications using directly modulated 520 nm laser diode,” Opt. Express, vol. 23, no. 16, pp. 20 743–20 748, Aug 2015.
 [5] C. Shen, Y. Guo, H. M. Oubei, T. K. Ng, G. Liu, K.H. Park, K.T. Ho, M.S. Alouini, and B. S. Ooi, “20meter underwater wireless optical communication link with 1.5 Gbps data rate,” Opt. Express, vol. 24, no. 22, pp. 25 502–25 509, Oct 2016.
 [6] T. Ballal, T. Y. AlNaffouri, and S. F. Ahmed, “Lowcomplexity bayesian estimation of clustersparse channels,” IEEE Trans. Commun., vol. 63, no. 11, pp. 4159–4173, Nov. 2015.
 [7] H. Kaushal and G. Kaddoum, “Underwater optical wireless communication,” IEEE Access, vol. 4, pp. 1518–1547, 2016.
 [8] N. Saeed, A. Celik, T. Y. AlNaffouri, and M.S. Alouini, “Connectivity analysis of underwater optical wireless sensor networks: A graph theoretic approach,” in IEEE Intl. Conf. Commun. Workshops (ICC Workshops), May 2018, pp. 1–6.
 [9] N. Saeed, A. Celik, T. Y. AlNaffouri, and M. Alouini, “Underwater optical sensor networks localization with limited connectivity,” in IEEE Intl. Conf. Acoust, Speech and Signal Process. (ICASSP), Apr. 2018, pp. 3804–3808.
 [10] S. Arnon, “Underwater optical wireless communication network,” Optical Engineering, vol. 49, pp. 49 – 49 – 6, 2010.
 [11] N. Anous, M. Abdallah, M. Uysal, and K. Qaraqe, “Performance evaluation of los and nlos vertical inhomogeneous links in underwater visible light communications,” IEEE Access, vol. 6, pp. 22 408–22 420, 2018.
 [12] S. Hessien, S. C. Tokgöz, N. Anous, A. Boyacı, M. Abdallah, and K. A. Qaraqe, “Experimental evaluation of ofdmbased underwater visible light communication system,” IEEE Photonics Journal, vol. 10, no. 5, pp. 1–13, Oct 2018.
 [13] T. Shafique, O. Amin, M. Abdallah, I. S. Ansari, M.S. Alouini, and K. Qaraqe, “Performance analysis of singlephoton avalanche diode underwater vlc system using arq,” IEEE Photonics Journal, vol. 9, no. 5, pp. 1–11, Oct 2017.
 [14] A. Vavoulas, H. G. Sandalidis, and D. Varoutas, “Underwater optical wireless networks: A connectivity analysis,” IEEE J. Ocean. Eng., vol. 39, no. 4, pp. 801–809, Oct 2014.
 [15] N. Saeed, A. Celik, S. Alouini, and T. Y. AlNaffouri, “Performance analysis of connectivity and localization in multihop underwater optical wireless sensor networks,” To appear in IEEE Transactions on Mobile Computing, pp. 1–1, 2018.
 [16] F. Akhoundi, J. A. Salehi, and A. Tashakori, “Cellular underwater wireless optical CDMA network: Performance analysis and implementation concepts,” IEEE Trans. Commun., vol. 63, no. 3, pp. 882–891, March 2015.
 [17] F. Akhoundi, M. V. Jamali, N. B. Hassan, H. Beyranvand, A. Minoofar, and J. A. Salehi, “Cellular underwater wireless optical CDMA network: Potentials and challenges,” IEEE Access, vol. 4, pp. 4254–4268, 2016.
 [18] M. V. Jamali, F. Akhoundi, and J. A. Salehi, “Performance characterization of relayassisted wireless optical CDMA networks in turbulent underwater channel,” IEEE Trans. Wireless Commun., vol. 15, no. 6, pp. 4104–4116, June 2016.
 [19] M. V. Jamali, A. Chizari, and J. A. Salehi, “Performance analysis of multihop underwater wireless optical communication systems,” IEEE Photon. Technol. Lett., vol. 29, no. 5, pp. 462–465, March 2017.
 [20] M. Elamassie, F. Miramirkhani, and M. Uysal, “Performance characterization of underwater visible light communication,” IEEE Trans. Commun., pp. 1–1, 2018.
 [21] J. M. Kahn and J. R. Barry, “Wireless infrared communications,” Proc. IEEE, vol. 85, no. 2, pp. 265–298, Feb 1997.
 [22] N. Saeed, A. Celik, T. Y. AlNaffouri, and M.S. Alouini, “Robust 3D localization of underwater optical wireless sensor networks via low rank matrix completion,” in IEEE 19th Intl. Workshop Signal Process. Adv. Wireless Commun. (SPAWC), Jun. 2018, pp. 1–5.
 [23] N. Saeed, A. Celik, T. Y. AlNaffouri, and M. Alouini, “Energy harvesting hybrid acousticoptical underwater wireless sensor networks localization,” Sensors, vol. 18, no. 1, 2017.
 [24] S. Arnon and D. Kedar, “Nonlineofsight underwater optical wireless communication network,” J. Opt. Soc. Am. A, vol. 26, no. 3, pp. 530–539, Mar 2009.

[25]
M. Fernandez and S. Williams, “Closedform expression for the poissonbinomial probability density function,”
IEEE Trans. Aerosp. Electron. Syst., vol. 46, no. 2, pp. 803–817, April 2010.  [26] G. P. Agrawal, FiberOptic Communication Systems. John Wiley & Sons, 2012, vol. 222.
 [27] S. Karp, R. Gagliardi, S. Moran, and L. Stotts, “Optical Channels: Fibers,” Clouds, Water, and the Atmosphere: Plenum Press New York, 1988.
 [28] V. Kaibel and M. A. Peinhardt, On the Bottleneck Shortest Path Problem. KonradZuseZentrum für Informationstechnik, 2006.
 [29] D. Eppstein, “Finding the k shortest paths,” SIAM Journal on Computing, vol. 28, no. 2, pp. 652–673, 1998.
 [30] A. Celik and A. E. Kamal, “Multiobjective clustering optimization for multichannel cooperative spectrum sensing in heterogeneous green crns,” IEEE Transactions on Cognitive Communications and Networking, vol. 2, no. 2, pp. 150–161, Jun. 2016.
 [31] A. Celik, A. Alsharoa, and A. E. Kamal, “Hybrid energy harvestingbased cooperative spectrum sensing and access in heterogeneous cognitive radio networks,” IEEE Transactions on Cognitive Communications and Networking, vol. 3, no. 1, pp. 37–48, Mar. 2017.
 [32] A. Celik and A. E. Kamal, “Green cooperative spectrum sensing and scheduling in heterogeneous cognitive radio networks,” IEEE Transactions on Cognitive Communications and Networking, vol. 2, no. 3, pp. 238–248, Sept 2016.
Comments
There are no comments yet.