Next Article in Journal
High Spatial Resolution Detector System Based on Reconfigurable Dual-FPGA Approach for Coincidence Measurements
Next Article in Special Issue
Inter-Mode Crosstalk Estimation between Cores for LPmn Modes in Weakly Coupled Few-Mode Multicore Fiber with Perturbations
Previous Article in Journal
Enhanced Hand Gesture Recognition with Surface Electromyogram and Machine Learning
Previous Article in Special Issue
Novel Optical Modulator Photonic Device Based on TiN/Ti3C2 Heterojunction
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Performance Evaluation of UOWC Systems from an Empirical Channel Model Approach for Air Bubble-Induced Scattering

by
Pedro Salcedo-Serrano
*,
Rubén Boluda-Ruiz
,
José María Garrido-Balsells
,
Beatriz Castillo-Vázquez
,
Antonio Puerta-Notario
and
Antonio García-Zambrana
Telecommunication Research Institute (TELMA), Universidad de Málaga, E-29010 Málaga, Spain
*
Author to whom correspondence should be addressed.
Sensors 2024, 24(16), 5232; https://doi.org/10.3390/s24165232
Submission received: 4 July 2024 / Revised: 29 July 2024 / Accepted: 11 August 2024 / Published: 13 August 2024
(This article belongs to the Special Issue Novel Technology in Optical Communications)

Abstract

:
Underwater optical wireless communication (UOWC) systems provide the potential to establish secure high-data-rate communication links in underwater environments. The uniqueness of oceanic impairments, such as absorption, scattering, oceanic turbulence, and air bubbles demands accurate statistical channel models based on empirical measurements for the development of UOWC systems adapted to different types of water and link conditions. Recently, generalized Gamma and a mixture of two generalized Gamma probability density functions (PDF) were proposed to describe the statistical behavior of small and large air bubbles, respectively, when considering several levels of particle-induced scattering. In this paper, we derive novel closed-form analytic expressions to compute the bit error rate (BER) and outage performance using both proposed PDFs for various scattering conditions. Furthermore, simple asymptotic expressions are obtained to determine the diversity order of each scenario. Monte Carlo simulation results verify the obtained theoretical expressions. Our results also reveal that UOWC systems present lower BER and outage performance under more turbid water cases with respect to the tap water case due to the higher diversity order and despite the significant increases in pathloss at short link distances. Particle-induced scattering provides an inherent mechanism of turbid waters to mitigate air bubble-induced fluctuations and light blockages.

1. Introduction

The underwater medium presents significant challenges for wireless communication systems that aim to extend the telecommunication networks to marine environments [1,2]. Recent bandwidth-hungry applications such as high-definition and real-time video transmission or remote control of underwater vehicles require robust links with low delay latency. At the same time, some of these applications require high transmission security due to the sensitivity of the information. In this context, underwater optical wireless communication (UOWC) systems offer a promising alternative to acoustic wireless communication systems, especially over short and moderate link distances [3]. However, the underwater optical channel presents physical challenges that can degrade the performance of UOWC systems that are not adequately designed for underwater-specific conditions. Traditional theoretical models, often based on Monte Carlo simulations, are limited by assumptions such as idealizing scattering particles as spherical or treating air bubbles under specific conditions [4]. These assumptions may not accurately reflect the true complexity of the marine environment. Some of these assumptions are related to the absorption and scattering effect, oceanic turbulence, and the presence of air bubbles. Each of these effects exhibits different characteristics and affects the performance of the UOWC system in different ways. However, the empirical approach provides a more accurate and realistic depiction of the underwater optical channel, integrating realistic environmental conditions rather than relying on simplified theoretical assumptions. Therefore, reliable UOWC system designs for specific oceanic environments require a performance analysis based on realistic and experimental underwater channel models.
Given the difficulty of conducting experiments in the open ocean, underwater optical channel models are typically developed in laboratory water tanks under controlled conditions. In this respect, several studies have focused on the experimental modeling of UOWC channels under different phenomena, such as random temperature and salinity variations or air bubbles [5,6,7,8,9,10,11] (and references therein). Numerous works have demonstrated the importance of these models by obtaining new closed-form expressions for evaluating the UOWC system performance under these experimental channel models. The vast majority of these works have concentrated on analyzing UOWC systems under oceanic turbulence by means of probability density functions (PDF) such as the Weibull distribution [12,13], the Exponential-Generalized Gamma distribution [14,15], and the Gamma–Gamma distribution [16]. However, although these models comply with current data of irradiance fluctuations due to the joint effect of thermohaline gradient and air bubbles, neither of them have been evaluated for optical power fluctuations due exclusively to air bubbles. In this regard, the thermohaline gradient can mask the real impact of air bubbles on the perturbed wavefront at the receiver. Furthermore, oceanic turbulence due to the inhomogeneous thermohaline distribution represents a minor drawback compared with air bubbles or the absorption and scattering effect in deep waters, as the temperature and salinity of the underwater medium tend to stabilize at constant values [17]. However, the origin of air bubbles is highly heterogeneous (cavitation of vehicle propellers, phytoplankton photosynthesis, or zooplankton respiration [18,19]), and their presence can significantly distort the received footprint. Indeed, a few investigations have demonstrated that the performance degradation in UOWC links is due solely to air bubbles [20,21]. In [20], the dependence of the bit error rate (BER) of a UOWC link using different orders of pulse position modulation on bubble density and bubble size is experimentally analyzed in a water tank. In [21], closed-form expressions for the ergodic capacity and the BER were obtained over a composite channel model which considers the bubble-obstruction and turbulence effect. However, the proposed statistical model is purely mathematical and numerically obtained by using statistics of the generation, size, and horizontal movement of air bubbles. Unlike empirical measurements-based channel models which consider the complicated nature of real underwater environments, statistical channel models derived from theoretical simulations often require simplifying and assuming specific channel parameters. This can result in significant discrepancies between the simulated results and the actual behavior of underwater environments. Furthermore, it should be noted that the aforementioned works neglect the dispersive nature of the water, by ignoring different types of water in the performance analysis.
Recently, novel research has proposed a statistical UOWC channel model based on experimental measurements that considers the combined effect of the scintillation and obstruction of light due to air bubbles and the water turbidity [22]. Then, small air bubble-induced fluctuations are described by a generalized Gamma distribution, and large air bubble-induced fluctuations and light blockage are described by a mixture of two generalized Gamma distributions. Channel measurements suggest that the scattering effect reduces both the scintillation effect and the obstruction of light caused by air bubbles, by inducing a more resilient link at the expense of a higher pathloss value. However, the current research has ignored the scattering phenomenon because it makes use of statistical channel models that are not fitted to realistic underwater environments. In other words, the effectiveness of UOWC systems under air bubbles in turbid waters has not yet been investigated. Therefore, there remains a need for mathematical tools and expressions that can complete the UOWC system performance of submarine links when air bubbles of different sizes perturb the optical wavefront under different types of water.
In this paper, we analyze the UOWC system performance in terms of bit error rate (BER) and outage probability over underwater fading channels in the presence of small and large air bubbles and different levels of scattering following the mentioned empirical-based UOWC channel models [22]. Unlike [5,7,11], where the estimated PDF and the corresponding channel parameters are fitted when considering a joint effect of thermohaline gradient and air bubbles in tap water, the UOWC system is evaluated under UOWC statistical channel models where fitted PDFs and its parameters are based on measurements with optical power fluctuations due to uniquely from air bubbles of two different sizes and a wide range of in-suspension particle concentrations.
Although previous research has derived analytical expressions for BER using a generalized Gamma distribution, it is essential to note that these studies do not address the performance of optical systems under air bubble-induced fading and particle-induced scattering [23,24,25]. Furthermore, some of these papers utilize fitted parameters derived from numerical simulation, which do not accurately capture the complexity of the underwater environment. Consequently, their findings may have limitations when applied to optimizing and designing novel UOWC systems in realistic conditions, such as environments with air bubbles in turbid water. This work extends the knowledge of existing UOWC system performance based on empirical UOWC channel models because we analyze the impact of water turbidity, which considers the oceanic pathloss and the beamwidth expansion due to the scattering effect. In particular, we provide exact and asymptotic closed-form expressions for BER and outage probability over air bubble-induced fading channels. In fact, to our knowledge, there is no asymptotic analysis of BER and outage probability presented in the literature that can obtain novel insights into how channel model parameters affect the overall performance at high signal-to-noise (SNR) ratio. Therefore, by using the asymptotic behavior, a novel diversity order analysis for the different scattering scenarios is included. The presented results show that UOWC links under higher levels of scattering exhibit a lower average BER and outage probability than tap water scenarios. It follows, therefore, that particle-induced scattering can mitigate the random irradiance fluctuations and light blockages due to both small and large air bubbles. Monte Carlo simulations are included to verify the presented results.
The remainder of this article is organized as follows. In Section 2, the system and channel models adopted in this work are described for small and large air bubbles from experimental measurements. In Section 3, we derive new closed-form analytic and asymptotic expressions for the average BER and outage probability when adopting the proposed channel models. Numerical results and discussions are presented and analyzed in Section 4. Finally, we summarize the main conclusions and some future lines of this work in Section 5.

2. System and Channel Models

Figure 1 illustrates a UOWC link between two autonomous underwater vehicles (AUV). In a practical underwater scenario, the transmitted beam is influenced by the underwater particles in terms of absorption and scattering, as well as air bubbles, which are generated by the cavitation effect of AUV propellers, marine flora and fauna respiration, and oxygen systems of scuba divers. A comprehensive model that contemplates the relation of these impairments in realistic scenarios is essential for accurately characterizing the UOWC channel and assessing the performance of UOWC systems. The following section presents the system and channel models adopted in this work, considering the environmental conditions prevalent in underwater settings from experimental measurements presented in [22].

2.1. System Model

Consider a single-input/single-output (SISO) UOWC system with an intensity modulation and direct detection (IM/DD) scheme due to the lower complexity and low cost with respect to coherent schemes [26]. The received electrical signal for the proposed UOWC system is given by
y = L · h b · R · x + n ,
where L is the pathloss due to the water turbidity and air bubbles, h b is the fading coefficient of the UOWC channel due to air bubbles, R is the detector responsivity, assumed hereinafter to be the unity, x is the transmitted optical power, and n is additive white Gaussian noise with zero mean and variance N 0 / 2 . Hence, the received average electrical SNR in the absence of fading can be defined as γ = P o p t 2 T b / N 0 , where P o p t is the average transmitted optical power, and T b is the bit period.

2.2. Air Bubble-Induced Fading Model

Air bubble-induced fading is modeled using the generalized Gamma distribution and a mixture of the sum of two generalized Gamma distributions when considering small air bubbles and large air bubbles, respectively [22]. The generalized Gamma distribution is a versatile statistical model that includes the Gamma distribution, the Chi and Chi-squared distribution, the Exponential distribution, and the Weibull distribution, which are usually employed in modeling power fluctuations in radio-frequency (RF) and optical wireless systems [7,23,24]. In the context of underwater environments, the Exponential-Generalized Gamma mixture has been used to analyze the multihop UOWC systems performance under oceanic turbulence [25,27]. Furthermore, in [22], the generalized Gamma distribution and the mixture of the sum of two generalized Gamma distributions provide an excellent agreement for a wide range of water turbidity scenarios, as well as different air bubble populations with experimental measurements under controlled laboratory conditions [22]. It should be noted that the small and large terms are defined in relation to the beam width of the transmitted beam. When the diameter of an air bubble is comparable to or bigger than such a beam width, it can result in a partial or complete light blockage. In such instances, this constitutes a scenario of large air bubbles. Concerning the PDF, large air bubbles exhibit a bimodal behavior due to the combined effects of light attenuation caused by large bubbles and power fluctuations induced by particles and small bubbles formed through bubble collisions. Conversely, scenarios involving smaller bubbles are assumed to represent instances of small air bubbles, which induce power fluctuations around a given received mean power. All details regarding fitting the PDF to the experimental measurements can be found in [22].
The corresponding PDF for small air bubble-induced fading is given by ([22], Equation (7)).
f h small ( h ; a , d , p ) = p a d Γ ( d / p ) h d 1 e ( h / a ) p ,
where h is the random variable, Γ ( · ) is the Gamma function, a is a positive scale parameter, and d and p are positive shape parameters [28]. In the same way, the considered PDF of large air bubble-induced fading is given by ([22], Equation (8))
f h large ( h ; a 1 , d 1 , p 1 , a 2 , d 2 , p 2 , w ) = w · f h small ( h ; a 1 , d 1 , p 1 ) + ( 1 w ) · f h small ( h ; a 2 , d 2 , p 2 ) ,
where w is the proportion between the blockage effect, the first term, and the power fluctuation effect, the second term, such that w [ 0 , 1 ] .
By fitting f h small ( h ) and f h large ( h ) channel models to experimental UOWC measurements, the specific values of the distribution parameters, as well as some parameters related to the UOWC channel such as the experimental propagation losses, L T , the extinction coefficient of the UOWC link, c, the scintillation index, σ 2 , and the coefficient of determination, R 2 , are obtained with a 520 nm laser diode at a link distance of 3 m. For detailed information on the experimental equipment and the setup used in the measurements, including all relevant parameters and configurations, please refer to the comprehensive description provided in [22]. A summary of the parameters of each distribution and the remaining parameters derived from the measurement of the optical power fluctuations for different levels of commercial antacid is presented in Table 1 and Table 2 when considering small and large air bubble scenarios, respectively. As reported, both proposed models accurately fit the random behavior of the received optical power for both small and large air bubbles, with a coefficient of determination R 2 exceeding 0.95 for all evaluated scattering conditions.

3. Performance Analysis of UOWC Links under Air Bubble-Induced Fading

In this section, the performance of a SISO UOWC system under the proposed underwater channel in the presence of air bubbles and different water turbidity levels is analyzed using the average BER and the outage probability. Moreover, asymptotic expressions are obtained to shed light on the impact of the channel parameters on the system performance by revealing the diversity order and the coding gain.

3.1. Bit Error Rate Analysis

As a case study, we analyze the BER performance of a SISO UOWC system when adopting On-Off Keying (OOK) signaling for the case of perfectly known channel state information at the receiver. As stated in [29], the conditional average BER at the receiver for the case of equally likely transmitted symbols is given by
P b ( E | h ) = Q 2 γ · L · h b ,
where Q ( · ) represents the Gaussian-Q function. Hence, P b in the small air bubble scenario is obtained by averaging over the PDF of small air bubbles defined in Equation (2) as follows:
P b small = 0 Q 2 γ · L · h b · f h small ( h b ) d h b .
By substituting Equation (2) into Equation (5), and making use of ([30], Equation (07.34.21.0012.01)), the exact closed-form expression for the average BER in the presence of small air bubbles can be obtained as follows
P b small ( L , a , d , p ) = p 4 π ( L a ) d Γ d p L d γ d 2 H 2 , 2 1 , 2 1 ( L a ) p γ p 2 | ( 1 d 2 , p 2 ) , ( 1 d 2 , p 2 ) ( 0 , 1 ) , ( d 2 , p 2 ) , ,
where H m , n p , q · is the H-Fox function ([31], Equation (1.1.1)).
Similarly, the average BER when considering large air bubbles is obtained by averaging over the PDF of large air bubbles defined in Equation (3) as follows:
P b large = 0 Q 2 γ · L · h b · f h large ( h b ) d h b .
From Equations (3) and (7) can be formulated as a weighted sum of integrals as follows:
P b large = 0 Q 2 γ · L · h b · w · f h small ( h b ; a 1 , d 1 , p 1 ) d h b + 0 Q 2 γ · L · h b · ( 1 w ) · f h small ( h b ; a 2 , d 2 , p 2 ) d h b .
Thus, P b large can be easily derived from Equation (5) as a weighted sum of P b small as follows:
P b large ( L , a 1 , d 1 , p 1 , a 2 , d 2 , p 2 , w ) = w · P b small ( L , a 1 , d 1 , p 1 ) + ( 1 w ) · P b small ( L , a 2 , d 2 , p 2 ) .
Since the obtained exact closed-form expressions could obscure the impact of PDF channel parameters on the UOWC system performance, more mathematically tractable expressions are obtained based on the asymptotic behavior at high SNR. According to ([32], Proposition 1), the asymptotic BER at high SNR can be derived with the behavior of the PDF of the considered random variable ( h b ) near the origin. Hence, the asymptotic BER at high SNR tends to P b ( G c γ ) G d , where G c is the coding gain, G d is the diversity order, and ≐ is the asymptotic equality sign [32]. On a log–log scale G c and G d specify a relative horizontal shift and the slope of the average BER and outage curves in the asymptotic regime, respectively. Both PDFs expressed in Equations (2) and (3) can be approximated by a single polynomial term for h b 0 , i.e., near the origin, as follows
f h small ( h b ) p a d Γ ( d / p ) h b d 1 ,
and
f h large ( h b ) w p 1 a 1 d 1 Γ ( d 1 / p 1 ) h b d 1 1 + ( 1 w ) p 2 a 2 d 2 Γ ( d 2 / p 2 ) h b d 2 1 .
Then, by substituting Equation (10) into Equation (5), and making use of ([30], Equation (07.34.21.0009.01)), the asymptotic closed-form expression for the average BER in the presence of small air bubbles can be solved as follows:
P b small p Γ 1 + d 2 2 π d ( L a ) d Γ d p γ d 2 .
Since the P b large is a weighted sum of P b small , the asymptotic closed-form expression for the average BER can be solved as follows:
P b large w 2 π p 1 Γ 1 + d 1 2 ( L a 1 ) d 1 d 1 Γ d 1 p 1 γ d 1 2 + 1 w 2 π p 2 Γ 1 + d 2 2 ( L a 2 ) d 2 d 2 Γ d 2 p 2 γ d 2 2 .

3.2. Outage Probability Analysis

An outage probability analysis for a SISO UOWC system under the UOWC channel model in the presence of air bubbles proposed in [22] is presented here to assess the completeness of the UOWC system performance. The outage probability is defined as the instantaneous SNR falls below a predefined threshold, γ t h as [29]
OP = P ( γ T γ t h ) ,
where γ T = 4 · γ · h b · L 2 . Therefore, Equation (14) can be expressed as
OP = P ( 4 · γ · h b · L 2 γ t h ) = 0 γ t h / 4 L 2 γ f h b ( h ) d h = F h b 1 4 L 2 γ ¯ ,
where γ ¯ = γ γ t h is the normalized SNR, and F h b ( · ) is the cumulative distribution function (CDF). The definition of the normalized SNR allows for the extrapolation of outage probability results to any predefined threshold value. The corresponding CDF when considering small air bubbles is obtained as follows ([28], Equation (2)).
F h small ( h ; a , d , p ) = ν ( d p , ( h / a ) p ) Γ d p ,
where ν ( · , · ) is the lower incomplete Gamma function [33]. Hence, by substituting Equation (16) into Equation (2.3) of [33], F h small ( · ) can be easily derived, as follows:
F h small ( h ; a , d , p ) = 1 Γ ( d p , ( h / a ) p ) Γ d p ,
where Γ ( · , · ) is the upper incomplete Gamma function [33]. Therefore, the corresponding outage probability when small air bubbles are considered can be easily derived by substituting Equation (17) into (15), as follows:
OP small = 1 Γ d p , 1 ( a L ) p γ ¯ p 2 Γ d p .
In the same way that the average BER, the outage probability at high SNR also tends to OP ( O c γ ) O d , where O c is the coding gain, and O d is the diversity order. Hence, by replacing Equation (10) into Equation (15), we obtain the asymptotic behavior of the outage probability at high SNR for small air bubbles as follows:
OP small p ( a L ) d d Γ d p γ ¯ d 2 .
Based on Equation (3), the CDF of large air bubble-induced fading can be expressed as
F h large h ; a 1 , d 1 , p 1 , a 2 , d 2 , p 2 , w = w · F h small ( h ; a 1 , d 1 , p 1 ) + ( 1 w ) F h small ( h ; a 2 , d 2 , p 2 ) .
Hence, the outage probability in the presence of large air bubbles is derived by substituting Equation (17) into Equation (20) as follows:
OP large = w 1 Γ d 1 p 1 , 1 ( a 1 L ) p 1 γ ¯ p 1 2 Γ d 1 p 1 + ( 1 w ) 1 Γ d 2 p 2 , 1 ( a 2 L ) p 2 γ ¯ p 2 2 Γ d 2 p 2 .
Finally, by applying Equation (11) into Equation (15), we obtain the asymptotic of the outage probability at high SNR for large air bubbles as follows:
O P large w 1 a 1 L d 1 Γ d 1 + p 1 p 1 γ ¯ d 1 2 + ( 1 w ) 1 a 2 L d 2 Γ d 2 + p 2 p 2 γ ¯ d 2 2 .

3.3. Diversity Order Analysis

As mentioned before, the average BER and the outage probability tend to P b G c γ G d and OP O c γ ¯ O d . As stated in [32], the BER and outage performance exhibit identical diversity orders for sufficiently large SNR, hence G d = O d . The resulting diversity order is useful for evaluating in a unifying manner and building insights about the impact of the underwater channel parameters on the UOWC system performance. From Equations (12) and (19), it can be shown that the diversity order, G d , when considering small air bubbles is given by
G d s m a l l = d 2 .
In the case of large bubbles, both the average BER and outage probability tend to P b G c 1 γ G d 1 + G c 2 γ G d 2 and OP O c 1 γ ¯ O d 1 + O c 2 γ ¯ O d 2 , where G d i = O d i . Therefore, the diversity order refers to the exponent of γ and γ ¯ in Equations (13) and (22), which determine the slope of the average BER and outage probability at high SNR, respectively. Mathematically, the smallest exponent dominates the behavior of these expressions. Thus, the diversity order can be expressed as the minimum of G d 1 and G d 2 . Consequently, the diversity order can be obtained from Equations (13) and (22) as follows:
G d l a r g e = min ( d 1 , d 2 ) 2 .

4. Numerical Results and Discussion

In this section, the average BER and outage probability of UOWC links under air bubble-induced scattering with different levels of water turbidity are evaluated. Furthermore, we provide Monte Carlo simulation results to verify the proposed analytical and asymptotic closed-form expressions.
Regarding the numerical simulation, a quasi-analytical Monte Carlo simulation approach with an equivalent noise source to improve computational efficiency has been implemented. Due to the long simulation time involved, simulation results only up to 10 9 are included in this manuscript. Therefore, the number of bits used in our simulations is 10 10 to ensure precision [34]. As mentioned above, OOK modulation is employed with a rectangular pulse shape. The output waveform takes predefined values corresponding to a 1 or 0, with each bit occurring with a probability of 0.5. It should be noted that no channel coding scheme was incorporated. The transmitted signal is detected using a matched filter at the receiver, with an impulse response matching the rectangular pulse shape. This filter maximizes the sampled SNR, and a maximum likelihood (ML) detector minimizes symbol error probability. Simulations are conducted in baseband, without considering carrier frequency and frequency band.
Concerning the empirical-based channel model, it must be noted that different channel conditions and fitted parameters are obtained from empirical measurements at a link distance of 3 m and a wavelength of 520 nm, which are summarized in Table 1 and Table 2 [22]. The measurement-based channel model allows for a more accurate estimation of the average BER and outage probability, accurately reflecting the complex relationships among absorption, particle-induced scattering, and air bubbles. Additionally, different levels of absorption and scattering were selected to represent a wide range of water types. For instance, for 0 mg/L of antacid, the extinction coefficient is c = 0.16   m 1 , which can be compared to the extinction coefficient of clear ocean water or the Jerlov IB water type [3,35]. In the case of the most turbid scenario, for 14.5 mg/L of antacid, c = 1.24   m 1 , the corresponding extinction coefficient is near the Jerlov 5C water type [35].

4.1. Diversity Order Performance

In Figure 2, we show the diversity order of each level of water turbidity when considering small air bubbles and large air bubbles when considering the same UOWC system. Firstly, in the presence of small air bubbles, the diversity order exhibits an exponentially increasing trend as the concentration of antacid grows, demonstrating a proportional relationship between the two variables. For large air bubbles, we observe that at the lowest level of antacid, i.e., 3.6 mg/L, the diversity order decreases slightly compared with tap water, i.e., 0 mg/L. It may be attributed to the unique and singular impact of partial and total light blockages caused by large bubbles. This effect diminishes at higher antacid concentrations with a more significant collection of scattered photons, as detailed in [22]. For this reason, after the initial singularity in the diversity order of the large air bubble scenario, both graphs provide a coherent and valuable insight into the exponential growth of diversity order in relation to water turbidity.

4.2. Small Air Bubble Scenario

In Figure 3, the average BER of a SISO UOWC system under small air bubbles for different levels of scattering are compared. As can be readily observed, the analytical and asymptotic results show an excellent agreement with Monte Carlo simulation results for all considered cases by validating the accuracy of the exact and asymptotic closed-form expressions in Equations (6) and (12), respectively. Firstly, it is interesting to note that for scenarios with a low antacid concentration, the asymptotic expression approaches Monte Carlo and analytic results even at low SNR faster than higher antacid concentrations. Secondly, as one would expect, the high pathloss due to increased antacid concentration deteriorates the average BER performance compared with the tap water case, as evident in the 3.6 mg/L scenario, which shows a higher average BER than the 0 mg/L scenario. However, the increase in diversity order due to higher water turbidity, as shown in Figure 2, reverses this trend at high SNR regime, especially for more turbid scenarios. While the cases with 3.6 mg/L and 7.3 mg/L of antacid exhibit higher average BER than the 0 mg/L case within the presented SNR range, the scenarios with the highest water turbidity, i.e., 11 mg/L and 14.5 mg/L of antacid concentration, significantly improve the performance of the UOWC system at high SNR under the same population of air bubbles. Thus, antacid-induced scattering mitigates the impact of air bubbles on the performance of UOWC systems. For instance, at 40 dB, the average BER for 14.5 mg/L is 1.35 × 10 1 , while for 0 mg/L, it is 1.7 × 10 3 . However, at 70 dB, the average BER for 14.5 mg/L is 6 × 10 6 , while for 0 mg/L, it is 1.2 × 10 5 . In practice, when the system is assured of having an SNR greater than 58 dB, the scenario with the highest turbidity will exhibit better performance in terms of average BER compared with the scenario with tap water. To illustrate this point, let us consider the natural optical beam spreading and geometric losses due to scattering. Hence, this effect can offer a great degree of robustness to air bubble-induced scattering and light blocking. These findings are in line with previous simulation and experimental results reported in [12,22], where the impact of pointing error and scintillation index also decreases as the water gets more turbid, respectively. Therefore, we conclude that this change in trend can be perfectly attributable to an increase in the severity of scattering. However, this interpretation is limited to the 3 m experimental water tank from which fitted channel model parameters have been obtained. The impact of scattering over longer link distances may not be accurately reflected. For shorter distances, the attenuation effect of turbid waters might be less significant than the beam spreading effect, potentially improving BER performance, as mentioned above. Nevertheless, over longer distances, the pathloss calculated from the Beer–Lambert law as e c d , where d is the link distance, becomes significantly greater in turbid waters, potentially resulting in higher BER performance. Therefore, beam spreading due to scattering can benefit short-range communication links by mitigating air bubble blockage, but overall performance can degrade over longer link distances due to increased pathloss.
The outage probability for small air bubbles is represented in Figure 4. As expected, the Monte Carlo simulation results agree with the analytical and asymptotic expression described in Equations (18) and (19), respectively. Note that the outage probability results reveal similar insights as the average BER into the impact of water turbidity on UOWC system performance.

4.3. Large Air Bubble Scenario

In Figure 5, we consider the average BER in the presence of large air bubble-induced fading for different water turbidity levels. As can be observed, the considered antacid concentration levels are similar to the small air bubble scenarios. Both analytic and asymptotic results demonstrate behavior similar to Monte Carlo simulation results, indicating the accuracy of the exact and asymptotic closed-form expressions obtained in Equations (9) and (13), respectively. In addition, Monte Carlo and analytic results support the convergence of the asymptotic solution, which is again faster for less turbid water cases. It can be observed that the asymptotic behavior of the scenario with a concentration of 14 mg/L of antacid does not fit as fast and accurately as the rest of the scattering scenarios. This is because this scenario exhibits a higher diversity order, i.e., a greater slope in the BER curve. Nonetheless, asymptotic results can be used as a tight upper bound across the presented SNR range.
As discussed in Section 4.2, the behavior of the average BER can be categorized into two regions. Before 60 dB, the average BER is lower when considering a tap water scenario, where the scattering effect is insignificant. However, after 60 dB, low average BER results are achieved in most turbid water scenarios. At 40 dB, the average BER for 14.5 mg/L is 9.1 × 10 2 , while for 0 mg/L, it is 2.1 × 10 2 . However, at 70 dB, the average BER for 14.5 mg/L is 9.4 × 10 7 , while for 0 mg/L, it is 4.2 × 10 5 . Although the scenarios with 3.6 mg/L and 7.3 mg/L of antacid show a higher average BER than the tap water scenario in all the evaluated SNR ranges, after 60 dB, the 7.3 mg/L scenario shows a lower average BER than the 3.6 mg/L scenario. This should be attributable to an increase in the absorption coefficient, i.e., a higher pathloss, in 3.6 mg/L and 7.3 mg/L scenarios with respect to the clear water case. However, it seems that the scattering coefficient remains constant for both cases.
In Figure 6, the outage performance results of the UOWC system under large air bubbles for different levels of water turbidity are plotted. Monte Carlo simulations results validate the analytic and asymptotic results obtained from Equations (21) and (22). The outage probability confirms the conclusions obtained from the above-average BER performance.

5. Conclusions

In this paper, novel closed-form expressions have been developed in order to compute the average BER and outage probability of UOWC systems under empirical underwater channels in the presence of small and large air bubbles when different scattering levels are contemplated. Furthermore, simple asymptotic expressions are also obtained, which allows us to analyze the diversity order of the UOWC system of each scenario. Monte Carlo simulations verify the obtained analytic and asymptotic results.
For the first time, the proposed system performance contemplates an empirical underwater optical channel model obtained from experimental data, which considers the impact of scattering in terms of pathloss and beam spreading. The considered UOWC measurement-based channel model is able to capture the inherent complexities and variability of absorption, scattering, and power fluctuations due to air bubbles with greater accuracy than theoretical models based on simulations. Hence, the accuracy of the empirical modeling improves the average BER and outage estimation under realistic environments. In particular, our results prove that the relative impact of air bubbles on BER performance becomes less significant as water turbidity increases at short link distances. The reason behind this is that the diversity order highly depends on the scattering effect, increasing for higher antacid concentrations, i.e., more turbid waters. Hence, our performance analysis provides novel insights into the existing UOWC systems in realistic water bodies, since the proposed fading statistical models in [22] consider the air bubble-induced fading and light blockage, as well as the water turbidity and the scattering effect.
Since the performance analysis of UOWC systems under a realistic oceanic medium is a critical investigation in current UOWC networks, our technical contribution is particularly important to the engineers working on UOWC systems optimization and design. Moreover, the use of experimental channel models allows for the validation and adjustment of theoretical models, thereby improving their robustness and practical applicability. Some interesting extensions of our work, including the performance evaluation of more sophisticated schemes such as coherent modulation for UOWC systems, as well as Multiple-Input/Multiple-Output (MIMO) UOWC systems, will be considered as the subject of our future research. Furthermore, the presented UOWC channel performance could be enhanced by exploring the impact of scattering and attenuation over extended link distances to provide a more comprehensive understanding of UOWC system performance.

Author Contributions

Conceptualization, methodology, and formal analysis, P.S.-S., R.B.-R., J.M.G.-B. and A.G.-Z.; writing—original draft preparation, P.S.-S.; writing—review and editing, P.S.-S., R.B.-R., J.M.G.-B., B.C.-V., A.P.-N. and A.G.-Z.; supervision, P.S.-S., R.B.-R., J.M.G.-B., B.C.-V., A.P.-N. and A.G.-Z.; project administration and funding acquisition, R.B.-R., J.M.G.-B. and A.G.-Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded in part by the Universidad de Málaga; in part by the Spanish Ministerio de Ciencia, Innovación y Universidades under Grant FPU22/02183; and in part by the Plan Andaluz de Investigación, Desarrollo e Innovación (PAIDI) 2021, under Grant P21_00390.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data will be made available on request.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

The following abbreviations are used in this manuscript:
AUVAutonomous Underwater Vehicles
BERBit Error Rate
CDFCumulative Distribution Function
IM/DDIntensity Modulation/Direct Detection
OOKOn-Off Keying
PDFProbability Density Function
SISOSingle-Input/Single-Output
SNRSignal-to-Noise Ratio
UOWCUnderwater Optical Wireless Communications

References

  1. Stojanovic, M. Underwater wireless communications: Current achievements and research challenges. IEEE Ocean. Eng. Soc. Newsl. 2006, 41, 1–5. [Google Scholar]
  2. Sun, X.; Kang, C.H.; Kong, M.; Alkhazragi, O.; Guo, Y.; Ouhssain, M.; Weng, Y.; Jones, B.H.; Ng, T.K.; Ooi, B.S. A Review on Practical Considerations and Solutions in Underwater Wireless Optical Communication. J. Light. Technol. 2020, 38, 421–431. [Google Scholar] [CrossRef]
  3. Zeng, Z.; Fu, S.; Zhang, H.; Dong, Y.; Cheng, J. A Survey of Underwater Optical Wireless Communications. IEEE Commun. Surv. Tutorials 2017, 19, 204–238. [Google Scholar] [CrossRef]
  4. Gabriel, C.; Khalighi, M.A.; Bourennane, S.; Léon, P.; Rigaud, V. Monte-Carlo-based channel characterization for underwater optical communication systems. J. Opt. Commun. Netw. 2013, 5, 1–12. [Google Scholar] [CrossRef]
  5. Oubei, H.M.; Zedini, E.; ElAfandy, R.T.; Kammoun, A.; Abdallah, M.; Ng, T.K.; Hamdi, M.; Alouini, M.S.; Ooi, B.S. Simple statistical channel model for weak temperature-induced turbulence in underwater wireless optical communication systems. Opt. Lett. 2017, 42, 2455–2458. [Google Scholar] [CrossRef] [PubMed]
  6. Oubei, H.M.; ElAfandy, R.T.; Park, K.H.; Ng, T.K.; Alouini, M.S.; Ooi, B.S. Performance Evaluation of Underwater Wireless Optical Communications Links in the Presence of Different Air Bubble Populations. IEEE Photonics J. 2017, 9, 1–9. [Google Scholar] [CrossRef]
  7. Jamali, M.V.; Mirani, A.; Parsay, A.; Abolhassani, B.; Nabavi, P.; Chizari, A.; Khorramshahi, P.; Abdollahramezani, S.; Salehi, J.A. Statistical Studies of Fading in Underwater Wireless Optical Channels in the Presence of Air Bubble, Temperature, and Salinity Random Variations. IEEE Trans. Commun. 2018, 66, 4706–4723. [Google Scholar] [CrossRef]
  8. Zedini, E.; Oubei, H.M.; Kammoun, A.; Hamdi, M.; Ooi, B.S.; Alouini, M.S. Unified statistical channel model for turbulence-induced fading in underwater wireless optical communication systems. IEEE Trans. Commun. 2019, 67, 2893–2907. [Google Scholar] [CrossRef]
  9. Yousefi, M.; Kashani, F.D.; Aghajani, A.; Rad, M.R.H. Experimental study of the effects of thermally induced optical turbulence on underwater wireless optical communication link parameters. J. Opt. 2020, 22, 025702. [Google Scholar] [CrossRef]
  10. Ke, X.; Bao, J.; Liang, J. Experimental study of underwater wireless optical channel model. Opt. Eng. 2023, 62, 124103. [Google Scholar] [CrossRef]
  11. Qiu, H.; Huang, Z.; Xu, J.; Zhai, W.; Gao, Y.; Ji, Y. Unified statistical thermocline channel model for underwater wireless optical communication. Opt. Lett. 2023, 48, 636–639. [Google Scholar] [CrossRef] [PubMed]
  12. Boluda-Ruiz, R.; García-Zambrana, A.; Castillo-Vázquez, B.; Hranilovic, S. Impact of angular pointing error on BER performance of underwater optical wireless links. Opt. Express 2020, 28, 34606–34622. [Google Scholar] [CrossRef] [PubMed]
  13. Boluda-Ruiz, R.; Salcedo-Serrano, P.; Castillo-Vázquez, B.; García-Zambrana, A.; Garrido-Balsells, J.M. Capacity of underwater optical wireless communication systems over salinity-induced oceanic turbulence channels with ISI. Opt. Express 2021, 29, 23142–23158. [Google Scholar] [CrossRef] [PubMed]
  14. Samir, A.; Elsayed, M.; El-Banna, A.A.A.; Shafique Ansari, I.; Rabie, K.; ElHalawany, B.M. Performance Analysis of Dual-Hop Hybrid RF-UOWC NOMA Systems. Sensors 2022, 22, 4521. [Google Scholar] [CrossRef] [PubMed]
  15. Rahman, Z.; Hassan, M.Z.; Kaddoum, G. Ergodic Capacity Optimization for RSMA-Based UOWC Systems Over EGG Turbulence Channel. IEEE Commun. Lett. 2024, 28, 587–591. [Google Scholar] [CrossRef]
  16. Ding, G.; Du, X.; Du, H.; Wang, S.; Feng, H.; Xu, G.; Xiong, Z.; Jia, Z.; Li, Y. Performance analysis of OSSK-UWOC systems considering pointing errors and channel estimation errors. Opt. Express 2024, 32, 3606–3618. [Google Scholar] [CrossRef] [PubMed]
  17. Weng, Y.; Guo, Y.; Alkhazragi, O.; Ng, T.K.; Guo, J.H.; Ooi, B.S. Impact of turbulent-flow-induced scintillation on deep-ocean wireless optical communication. J. Light. Technol. 2019, 37, 5083–5090. [Google Scholar] [CrossRef]
  18. Vahaji, S.; Chen, L.; Cheung, S.C.; Tu, J. Numerical investigation on bubble size distribution around an underwater vehicle. Appl. Ocean. Res. 2018, 78, 254–266. [Google Scholar] [CrossRef]
  19. Waaland, J.R.; Branton, D. Gas vacuole development in a blue-green alga. Science 1969, 163, 1339–1341. [Google Scholar] [CrossRef]
  20. Chen, D.; Wang, J.; Li, S.; Xu, Z. Effects of air bubbles on underwater optical wireless communication. Invited. Chin. Opt. Lett. 2019, 17, 100008. [Google Scholar] [CrossRef]
  21. Shin, M.; Park, K.H.; Alouini, M.S. Statistical Modeling of the Impact of Underwater Bubbles on an Optical Wireless Channel. IEEE Open J. Commun. Soc. 2020, 1, 808–818. [Google Scholar] [CrossRef]
  22. Salcedo-Serrano, P.; Gómez-García, C.; Iamaguti-Debessa, J.; Boluda-Ruiz, R.; Garrido-Balsells, J.M.; Castillo-Vázquez, B.; Puerta-Notario, A.; García-Zambrana, A. On the effect of air bubbles-induced scattering on turbid waters: An experimental UOWC channel modeling approach. IEEE Access 2024, 12, 75888–75898. [Google Scholar] [CrossRef]
  23. Aalo, V.A.; Piboongungon, T.; Iskander, C.D. Bit-error rate of binary digital modulation schemes in generalized gamma fading channels. IEEE Commun. Lett. 2005, 9, 139–141. [Google Scholar] [CrossRef]
  24. AlQuwaiee, H.; Ansari, I.S.; Alouini, M.S. On the performance of free-space optical communication systems over double generalized gamma channel. IEEE J. Sel. Areas Commun. 2015, 33, 1829–1840. [Google Scholar] [CrossRef]
  25. Le-Tran, M.; Kim, S. Performance analysis of multi-hop underwater wireless optical communication systems over exponential-generalized gamma turbulence channels. IEEE Trans. Veh. Technol. 2022, 71, 6214–6227. [Google Scholar] [CrossRef]
  26. Barry, J.R. Wireless Infrared Communications; Springer Science & Business Media: Berlin/Heidelberg, Germany, 1994; Volume 280. [Google Scholar]
  27. Zedini, E.; Kammoun, A.; Soury, H.; Hamdi, M.; Alouini, M.S. Performance analysis of dual-hop underwater wireless optical communication systems over mixture exponential-generalized gamma turbulence channels. IEEE Trans. Commun. 2020, 68, 5718–5731. [Google Scholar] [CrossRef]
  28. Stacy, E.W. A Generalization of the Gamma Distribution. Ann. Math. Stat. 1962, 33, 1187–1192. [Google Scholar] [CrossRef]
  29. Goldsmith, A. Wireless Communications; Cambridge University: Cambridge, UK, 2005. [Google Scholar]
  30. Wolfram Research, Inc. The Wolfram Functions Site. Available online: https://functions.wolfram.com/ (accessed on 1 July 2024).
  31. Kilbas, A.A. H-Transforms: Theory and Applications; CRC Press: Boca Raton, FL, USA, 2004. [Google Scholar]
  32. Wang, Z.; Giannakis, G.B. A simple and general parameterization quantifying performance in fading channels. IEEE Trans. Commun. 2003, 51, 1389–1398. [Google Scholar] [CrossRef]
  33. Chaudhry, M.A.; Zubair, S.M. On a Class of Incomplete Gamma Functions with Applications; Chapman and Hall/CRC: Boca Raton, FL, USA, 2001. [Google Scholar]
  34. Jeruchim, M.C.; Balaban, P.; Shanmugan, K.S.; Jeruchim, M.C.; Balaban, P.; Shanmugan, K.S. Modeling of Communication Systems. In Simulation of Communication Systems. Applications of Communications Theory; Springer: Boston, MA, USA, 1992; pp. 303–462. [Google Scholar]
  35. Williamson, C.A.; Hollins, R.C. Measured IOPs of Jerlov water types. Appl. Opt. 2022, 61, 9951–9961. [Google Scholar] [CrossRef]
Figure 1. Schematic of a UOWC link application between two AUVs under air bubble scenario due to oxygen systems of scuba divers and AUV propellers.
Figure 1. Schematic of a UOWC link application between two AUVs under air bubble scenario due to oxygen systems of scuba divers and AUV propellers.
Sensors 24 05232 g001
Figure 2. Diversity order when considering small air bubbles and large air bubbles for different levels of antacid in the experimental UOWC link.
Figure 2. Diversity order when considering small air bubbles and large air bubbles for different levels of antacid in the experimental UOWC link.
Sensors 24 05232 g002
Figure 3. Average BER performance when considering an empirical UOWC channel in the presence of small air bubble-induced fading under different levels of antacid concentration at a link distance of 3 m and a wavelength of 520 nm.
Figure 3. Average BER performance when considering an empirical UOWC channel in the presence of small air bubble-induced fading under different levels of antacid concentration at a link distance of 3 m and a wavelength of 520 nm.
Sensors 24 05232 g003
Figure 4. Outage performance when considering an empirical UOWC channel in the presence of small air bubble-induced fading under different levels of antacid concentration at a link distance of 3 m and a wavelength of 520 nm.
Figure 4. Outage performance when considering an empirical UOWC channel in the presence of small air bubble-induced fading under different levels of antacid concentration at a link distance of 3 m and a wavelength of 520 nm.
Sensors 24 05232 g004
Figure 5. Average BER performance when considering an empirical UOWC channel in the presence of large air bubble-induced fading under different levels of antacid concentration at a link distance of 3 m and a wavelength of 520 nm.
Figure 5. Average BER performance when considering an empirical UOWC channel in the presence of large air bubble-induced fading under different levels of antacid concentration at a link distance of 3 m and a wavelength of 520 nm.
Sensors 24 05232 g005
Figure 6. Outage performance when considering an empirical UOWC channel in the presence of large air bubble-induced fading under different levels of antacid concentration at a link distance of 3 m and a wavelength of 520 nm.
Figure 6. Outage performance when considering an empirical UOWC channel in the presence of large air bubble-induced fading under different levels of antacid concentration at a link distance of 3 m and a wavelength of 520 nm.
Sensors 24 05232 g006
Table 1. Experiment results of fitting parameters for generalized Gamma, total losses, the attenuation coefficient, and the scintillation index under small air bubble scenario for different turbidity levels.
Table 1. Experiment results of fitting parameters for generalized Gamma, total losses, the attenuation coefficient, and the scintillation index under small air bubble scenario for different turbidity levels.
Antacid [mg/L] ( a , d , p ) L T [dB]c [ m 1 ] σ 2 R 2
0(1.73, 1.43, 7.95)6.80.160.220.99
3.6(1.72, 1.44, 6.28)11.70.540.230.99
7.3(1.59, 1.68, 4)14.70.770.190.98
11(1.19, 2.75, 3.05)17.40.980.140.99
14.5(0.79, 5.09, 2.41)20.91.240.090.98
Table 2. Experiment results of fitting parameters for generalized Gamma, total losses, the attenuation coefficient, and the scintillation index under large air bubble scenario for different turbidity levels.
Table 2. Experiment results of fitting parameters for generalized Gamma, total losses, the attenuation coefficient, and the scintillation index under large air bubble scenario for different turbidity levels.
Antacid [mg/L] ( a 1 , d 1 , p 1 ) ( a 2 , d 2 , p 2 ) w L T [dB]c [ m 1 ] σ 2 R 2
0(0.028, 2, 1.2)(1.38, 29.38, 22.7)0.25.00.160.320.95
3.6(0.05, 1.110, 0.85)(1.38, 13.95, 12.8)0.239.90.540.300.98
7.3(0.06, 1.320, 0.85)(1.39, 11.93, 11.7)0.2312.90.770.290.98
11(0.047, 3.77, 1.060)(1.36, 10.87, 10.55)0.2615.60.980.270.98
14.5(0.037, 7.28, 0.98)(1.33, 10.87, 10.2)0.2519.11.240.210.98
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Salcedo-Serrano, P.; Boluda-Ruiz, R.; Garrido-Balsells, J.M.; Castillo-Vázquez, B.; Puerta-Notario, A.; García-Zambrana, A. Performance Evaluation of UOWC Systems from an Empirical Channel Model Approach for Air Bubble-Induced Scattering. Sensors 2024, 24, 5232. https://doi.org/10.3390/s24165232

AMA Style

Salcedo-Serrano P, Boluda-Ruiz R, Garrido-Balsells JM, Castillo-Vázquez B, Puerta-Notario A, García-Zambrana A. Performance Evaluation of UOWC Systems from an Empirical Channel Model Approach for Air Bubble-Induced Scattering. Sensors. 2024; 24(16):5232. https://doi.org/10.3390/s24165232

Chicago/Turabian Style

Salcedo-Serrano, Pedro, Rubén Boluda-Ruiz, José María Garrido-Balsells, Beatriz Castillo-Vázquez, Antonio Puerta-Notario, and Antonio García-Zambrana. 2024. "Performance Evaluation of UOWC Systems from an Empirical Channel Model Approach for Air Bubble-Induced Scattering" Sensors 24, no. 16: 5232. https://doi.org/10.3390/s24165232

APA Style

Salcedo-Serrano, P., Boluda-Ruiz, R., Garrido-Balsells, J. M., Castillo-Vázquez, B., Puerta-Notario, A., & García-Zambrana, A. (2024). Performance Evaluation of UOWC Systems from an Empirical Channel Model Approach for Air Bubble-Induced Scattering. Sensors, 24(16), 5232. https://doi.org/10.3390/s24165232

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop