Extending UWOC System Applications through Photon Transmission Dynamics Study in Harbor Waters

: Underwater wireless optical communication (UWOC) in harbor waters can facilitate real-time monitoring underwater instruments for environmental monitoring, underwater inspection, and maintenance tasks. This study delves into the complex dynamics of UWOC in four distinct harbor waters. The research employs Monte Carlo method incorporated with Fournier–Forand scattering phase function for simulating photon transmission. Key parameters such as the Transmi�ed full divergence angle, received aperture, and Field of View (FOV) are meticulously evaluated for their impact on power loss and time delay spread. Notably, the normalized power loss and time delay spread are found to be more signiﬁcantly aﬀected by communication distance than water quality , and the traditional Beer–Lambert law is ineﬀective in harbor water. The power loss of Harbor II, III, and IV are found to be 14.00 dB, 31.59 dB, and 41.59 dB lower than that of Harbor I, and the time delay spread of Harbor II, III, and IV is 30.56%, 9.67%, and 0.49% times that of the Harbor I under certain conditions. In addition, increasing the received aperture and FOV, particularly over longer distance, make li�le contribution to reduce the power loss and mitigate the time delay spread. Based on the ﬁxed transmi�ed full divergence angle, the most applicable received FOV range is 1–3.2 rad, and the most ideal received aperture is 0.15–0.4 m . Under these conditions, the variation in normalized power loss is less than 2 dB. Additionally, the time delay spread remains within the same order of magnitude with the a�enuation length (AL) held constant. These conclusions hold substantial technical relevance for the engineering design of UWOC systems in harbor waters.


Introduction
Comprehensive ocean exploration has emerged as a critical focus, garnering significant a entions within scientific, industrial, and military sectors.Key areas such as safety, weather monitoring, and marine commerce are heavily dependent on an in-depth understanding of the marine environment [1].Yet, exploring this complex terrain, particularly in deep ocean areas, is fraught with challenges.These include high costs, inherent dangers, and technical difficulties.In response to these challenges, there has been rapid advancement in technologies such as underwater wireless sensing, underwater imaging detection, and development of both manned and unmanned underwater vehicles [2].Consequently, the demand for advanced underwater communication technology, capable of high-rate data transmission, has intensified.Such technology is essential for not only transmi ing observational data, but also for gaining a deeper understanding of the marine environment [3].
Underwater wireless communication is a pivotal technology in ocean exploration and operations.Electromagnetic waves, typically employed for communication, face significant a enuation in underwater environment due to the high electrical conductivity of seawater.As a result, sound waves have become the prevalent medium for underwater wireless communication.Despite their ability to cover distances of kilometers, acoustic systems are limited by their transmission bandwidth, typically only reaching speeds of a few Kbps.Moreover, the slow propagation speed of sound underwater poses substantial challenges for real-time, high-speed communication [4].Further complicating the issue, acoustic communication equipment often suffers from drawbacks such as large size and high energy consumption, which impedes its development and widespread adoption.Finally, certain frequencies of acoustic signals can significantly harm marine organisms, affecting their health, predation behaviors, and migration pa erns.
In contrast, underwater wireless optical communication (UWOC) presents several advantages over traditional radio frequency, acoustic communication, and wired communication methods.UWOC is notable for its higher data transmission rates, lower latency, and greater flexibility [5][6][7].Capable of covering hundreds of meters and potentially achieving Gbps transmission rates, UWOC holds promising application potential in areas such as high-capacity data transmission in deep-sea environments, remote control and monitoring of underwater robotics, and deep-sea oil and gas exploration [8][9][10][11][12].These advantages position UWOC as a technology with significant potential for advancing the field of underwater communication.
Currently, the mechanical shell of UWOC system in practical applications is predominantly constructed by titanium alloy, an inert metal material.The system's optical window typically utilizes materials such as sapphire, acrylic, or quar glass [13].Additionally, the outer layers of the power supply and network transmission cables are designed to be waterproof, corrosion-resistant, and capable of withstanding pressure.In summary, the materials mentioned above are essentially harmless to aquatic organisms and marine environment in short-term exposure scenarios.Moreover, the light source's intensity used for communication is generally not excessively high.However, during prolonged underwater operation, marine microorganisms may adhere to the surface of shell and optical window, significantly impairing the communication system's functionality.
The transmission of light through seawater channels is significantly influenced by dynamic seawater conditions.During transmission, light absorption reduces the total propagation energy, while sca ering caused by water particles leads to increased pulse spread in the time domain of the optical signal [14].Additionally, multiple sca erings result in increased geometric loss in communication links [15].So, the optical signal's spatial and temporal distribution characteristics show strong randomness [16].The high cost of underwater measurement adds the complexity of accurately modeling and analyzing the optical signal variation in underwater channels, a topic that has garnered global research interests [17][18][19][20].Now the most literature about UWOC predominantly concentrates on light transmission properties and communication technology research in several water types, including clear seawater, coastal seawater, Class Harbor II water, pool water, and laboratory sink water [21][22][23].The research on different kinds of harbor waters are limited.There are few studies on the transmission characteristics of light in Harbor-I, Harbor-II, Harbor-III and Harbor-IV water.Studying UWOC in harbor water helps in deep understanding how light behaves in such environments and developing techniques to mitigate these effects.Furthermore, harbors are critical infrastructure with high security needs, and UWOC can facilitate real-time monitoring and control of underwater Autonomous Underwater Vehicles (AUVs), sensors, and instruments for environmental monitoring, un-derwater inspection, and maintenance tasks, which are essential for the sustainable management of harbor areas.Above all, the current engineering application of UWOC systems is mainly under deep sea and ocean water conditions [24][25][26], very li le research on laser signal transmission in various harbor waters has limited its wider application.
Furthermore, harbor waters are typically more turbid than the open ocean, primarily due to sediment presence, runoff, and human activities.The myriad substances found in harbor water-including salt, organic compounds, dissolved materials, and potential pollutants-have a profound effect on the light's absorption and sca ering, thereby affecting the absorption of photons [27].Specifically, higher concentrations of organic ma er and suspended particles in harbor waters tend to increase sca ering and absorption.The scattering and absorption effects of light are significantly influenced by the size, composition, and concentration of these particles.Additionally, the dynamics of bubble flow in 2D bubble columns can further affect the transmission characteristics of light signals [28].This is particularly relevant with methane, which is produced by the anaerobic actions of underwater microorganisms and oxygen, generated through the photosynthesis of aquatic plants in harbor waters.
In recent years, the Monte Carlo method has gained popularity for simulating underwater optical signal transmission [29][30][31][32][33][34][35][36].Gabriel C.'s analysis of underwater channel impulse response and time dispersion across various seawater types, link distances, and received apertures indicates that time dispersion is often negligible, reducing the likelihood of inter-symbol interference [32].The two-term Henyey-Greenstein Function has been employed for more accurate simulation of power distribution on the receiving surface [29,33].Shijian Tang's application of the Monte Carlo method, coupled with the two-term Gamma Function to represent impulse response curves, demonstrates that the two-term Gamma Function model is effective in scenarios with large a enuation lengths.[34,35].
Comparing with other types of water, the sca ering phase function commonly used at present are not suitable for Monte Carlo simulation in Harbor water, which may lead to incorrect estimation of channel characteristics.The Fournier-Forand sca ering phase function has been recognized in recent studies as a be er fit for the characteristics of Pe old Harbor water, particularly when compared to the traditional Henyey-Greenstein, two-term Henyey-Greenstein Function, and Sahu and Shanmugam functions [15].
In this study, we would explore the intricate dynamic characteristics of UWOC systems under four distinct Harbor water conditions, which is vital to develop and extend applications in these water conditions.The Monte Carlo method serves as the core simulation technique for modeling photon transport.The rest of the paper is organized in the following manner: Section 2 details the optical parameters of the underwater channel, the selection of an appropriate sca ering phase function and the process of the Monte Carlo numerical simulation.Section 3 presents the simulation results and provides a discussion on the normalized power loss and time delay spread as functions of the Transmi ed full divergence angle, Received aperture, FOV, and a enuation length (AL).Section 4 concludes the paper with a summary of the findings.This theoretical calculation offers a dependable technical reference for the engineering design of UWOC systems operating in highly turbid waters.

The Optical Parameters of the Underwater Channel
The interaction between photons and particles within seawater is both frequent and intense.When a laser beam travels through water, its light is absorbed by various elements such as water molecules, phytoplankton, dissolved organic ma er, and organic debris.Concurrently, the light undergoes sca ering caused by particles suspended in the seawater [37][38][39].To mathematically describe these effects, two wavelength-dependent optical characteristic parameters are utilized: the absorption coefficient, denoted as ( ), and the sca ering coefficient, denoted as ( ).Consequently, the a enuation coefficient, ( ), which represents the cumulative effect of both absorption and sca ering, can be expressed as [40]: In Equation ( 1), represents the wavelength of light.The ratio of the sca ering coefficient ( ) to the a enuation coefficient ( ) is known as the albedo, symbolized as .This ratio, detailed in Equation ( 2), indicates the proportion of sca ering loss relative to the total loss [41][42][43].For the purpose of this study, which aims to analyze the impact of water turbidity on UWOC, harbor water is categorized into four distinct types: Harbor-I, Harbor-II, Harbor-III, and Harbor-IV.A communication wavelength of 532 nm is chosen for analysis.
Table 1 presents the absorption coefficient ( ), sca ering coefficient ( ), a enuation coefficient ( ), and albedo ( ) for these four types of harbor water.The optical parameters for Harbor II water are sourced from Pe old's report [38].For the remaining three Harbor water types, their optical parameters are extrapolated using the same albedo, , while varying the a enuation coefficients, ( ). and ѱ represent the initial elevation angle and azimuth angle in Figures 1  and 2, respectively.It is important to note that the fixed z-axis is always perpendicular to the Reception Plane.Photons emi ed from the light source undergo multiple sca erings by particles within the water, following various trajectories.Some of these photons, termed 'On-target photons', reach the photosensitive surface of the receiver and are absorbed, as shown by the blue line.Others, known as 'Off-target photons', deviate from the receiver's aperture and FOV after multiple sca erings shown by the green line, as depicted in Figure 1.represent the sca ering angles of the photon at each event, while ѱ , ѱ , …ѱ denote the azimuth angles associated with these events, with indicating the -th photon sca ering occurrence.
A notable feature in this figure is the gray do ed lines, which depict the photon's motion coordinate system.This system undergoes continuous transformations with each sca ering event, reflecting the dynamic nature of the photon's path.The specifics of how this coordinate system is determined by the Equation ( 14).Additionally, it is important to highlight that the x, y, and z axes shown in Figure 2 are part of a fixed coordinate system.This system is established by the positioning of the transmi er and receiver, providing a stable reference frame for understanding the photon's motion and sca ering behavior.
Additionally, the interaction between seawater particles and photons results in two distinct types of light: Forward-sca ering light and Backsca ering light [21,36], as shown in Figure 3.In the simulation process, the receiver is positioned on the z-axis.The distance between the Transmission Plane and the Reception Plane is governed by a variable known as the a enuation length.The Transmission Plane lies within the x-y plane and is parallel to the Reception Plane, maintaining no angle of declination.This configuration ensures that both the Transmission Plane and Reception Plane are perpendicular to the z-axis, as depicted in Figure 3.This perpendicular arrangement is crucial for accurately simulating the propagation of light or signals in the system and analyzing their interactions with various environmental factors.
In underwater environments, the sca ering direction of a photon post-collision with seawater particles is generally random.To evaluate the distribution of light energy following such sca erings, the volume sca ering phase function is employed [44,45].In this study, this function is denoted as ( ), the refers to the newly acquired sca ering angle in the photon motion coordinate system following a sca ering event, it has nothing to do with and the z axis in Figure 1.Instead, is changing in the photon motion coordinate system, which is constantly changing with each sca ering event and is determined by Equation (14).The relationship between the volume sca ering phase function and the sca ering coefficient ( ) is articulated in Equation ( 3) [46].Here, ( , ) represents the sca ering phase function, which plays a crucial role in determining how light is redistributed after interacting with particles in seawater.
To date, several sca ering phase functions have been proposed to study the light transmission characteristics in underwater channels.Notable among these are the Henyey-Greenstein Function (HG) [46], the two-term Henyey-Greenstein Function (TTHG) [47], Sahu and Shanmugam Function (SS) [46], and the Fournier-Forand Function (FF) [48,49].The FF sca ering phase function, in particular, has been shown to align well with the characteristics of Pe old Harbor water, as indicated in the relevant literature [38].
Further supporting the use of the FF function, simulation results have demonstrated its similarity to the measured data of Pe old Harbor water, especially when specific parameters are applied: a slope (µ) of 3.5835 and a seawater refractive index (n) of 1.13 [15].Given these findings, this paper adopts the FF sca ering phase function for its simulation analysis, and the sca ering angle is generated according to the Equation ( 4), ensuring a direct and unbiased determination of the angle.The expression for the FF sca ering phase function is provided in Equation ( 4).In this expression, ѱ represents the azimuth angle of the photons, and denotes their sca ering angle.The variables and are intermediate variables without direct physical significance, with representing the value of when the sca ering angle equals to 180°. (4)

The Initial Parameters of Photon Packets
The Monte Carlo method is a stochastic simulation technique used to model the propagation of light in seawater.This method involves tracking the random propagation paths of numerous photons, thereby facilitating statistical analysis and calculation.In this approach, the light source is modeled as a collection of photon packets, each characterized by specific distribution properties.The simulation process is guided by the optical and water quality parameters pertinent to the communication system.The photon sca ering angle, determined by the volume sca ering phase function, plays a crucial role in this process.By tracking the trajectory of each photon, the method enables a detailed analysis of their motion paths.Ultimately, this simulation allows for the accumulation and analysis of data regarding the number of photons received and their distribution characteristics.Such comprehensive data are crucial for understanding and optimizing underwater communication systems.
In the Monte Carlo simulation, the initial parameters of a photon are defined by its initial coordinates ( , , ) on the Transmission Plane and the directional components ( , , ) of the photon's direction vector along the x, y, and z axes, respectively.These directional components, commonly referred to as Direction Cosines, determine the photon's propagation direction, while the coordinates indicate its current position.As illustrated in Figure 4, we assume that the photon packet travels along the z-axis and lies on the Transmission Plane.
In the depicted setup, the Transmission Plane is situated within the x-y plane.The solid blue arrow in Figure 4 represents the initial direction vector of the photon packet.Additionally, the dashed blue rectangle illustrates the initial photon motion coordinate system, which undergoes changes with each photon sca ering event.It is important to note that, despite the dynamic nature of the photon motion coordinate system, the coordinate system established by the transmi er and receiver remains constant.In other words, the x, y, and z axes, as shown in Figure 4, are fixed.This distinction between the changing photon motion coordinate system and the static system set by the transmi er and receiver is crucial for understanding the dynamics of photon movement and sca ering in the simulation.Here, ( , ) is the polar coordinate for the initial photon position in the Transmission plane, the distance between the photon packet's current position, and the center of the light source is denoted as .Additionally, the initial azimuth angle (ѱ ) and elevation angle ( ) of the photon are specified.The photon's initial weight is set to = 1.
The initial coordinates of the photon ( , , ) and the Direction Cosines ( , , ) are derived using trigonometric functions, as shown in Equations ( 5) and ( 6).
In the Monte Carlo simulation process, the beam of the light source does not necessarily have to be azimuthally isotropic.The initial angle and the position for a photon packet should always be random, irrespective of the symmetrical nature of the beam.and ѱ can be randomly determined, this random selection is guided by Equation (7), where represents a random value within the normalized interval [0, 1].This method ensures that the angles are uniformly distributed, aligning with the radial symmetry of the light source and beam configuration.
The energy distribution of the laser beam in the simulation follows Gaussian distribution characteristics.This distribution is crucial for determining the radius on the Transmission Plane, which represents the distance between a given photon and the center of the light source.The calculation of is based on the Gaussian distribution, as outlined in Equation (8).In this equation, represents the waist radius of the laser beam, and is a random value drawn from the normalized interval [0, 1].
The initial elevation angle of each photon is derived based on the initial divergence half-angle ( ) of the beam, the beam waist radius ( ), and the radius ( ) on the Transmission Plane.The specific relationship among these parameters is articulated in Equation ( 9) [15].This equation ensures that the elevation angle of each photon is accurately aligned with the beam's divergence characteristics.

=
The Equations ( 5)-( 9) discussed previously demonstrate that the initial elevation angle, initial azimuth angle, and initial position of each photon in the Monte Carlo simulation are determined through random calculations.This randomization process is crucial for accurately simulating the inherently stochastic nature of photon behavior in harbor water.

The Movement and Sca ering of the Photon
The Monte Carlo simulation process is comprehensively illustrated in Figure 5, which presents a detailed flowchart.A critical aspect of this simulation involves determining the random sca ering path length ( ), sca ering angle ( ), and azimuth angle (ѱ) for each photon during its movement through the simulated environment.The sca ering path length is defined as the distance a photon travels between two consecutive scattering events.This length is crucial as it influences the photon's trajectory and the overall simulation accuracy.
The sca ering angle of the photon after each collision with the particles in harbor water is randomly calculated by the sca ering phase function ( ) in Equation (3).The new azimuth angle ѱ is calculated randomly by the equation ѱ = 2 , and the random sca ering path length is calculated randomly by the Equation (11). is also a random value drawn from the normalized interval [0, 1] in the calculation of the propagation length and azimuth angle ѱ .The photon's weight is updated by the equation = * , where, represents the -th photon sca ering event.Two critical steps in the simulation are the block «Calculate Euclidian distance to Rx plane» and «Calculate y, z intersection point», which mean that if the Euclidian distance between the photon's position and the Reception Plane's center is smaller than the received aperture, and if the angle between the y axis of the photon moving coordinate system and the z axis of the coordinate system built by the transceiver terminal is smaller than the Received FOV, the photon is thought to be able to reach the Reception Plane, then record the corresponding sca ering angle meanwhile.Russian roule e is also introduced, and its function is: when the photon's weight is greater than the se ing threshold, the photon sca ering process continues to be tracked, and the subsequent flow and algorithm are consistent with the previous.When the photon's weight is less than the se ing threshold, the photon is considered to have been annihilated completely.
During the simulation, the Transmission Plane and the Reception Plane are parallel to each other, and the fixed coordinate axis z is always perpendicular to the Reception Plane, so the orientation of the receiver can be confirmed.The position of the receiver can be confirmed by the se ing values of the AL during simulation.
Equation ( 10) then provides further detail on how this angle is determined, with representing the specific sca ering angle chosen during a sca ering event, and still represents a random value drawn from the normalized interval [0, 1].These steps are integral to accurately simulating the random paths of photons through water, crucial for understanding light behavior in underwater communication channels.When the random sca ering path length , sca ering angle and azimuth angle ѱ of the photon are confirmed, the position of the photon can be updated as [ , , ] in Equation (12), where , , and are the Direction Cosine of the photon in the current state, the sca ering path length is selected randomly in Equation (11) [50,51], where is the sca ering coefficient, and is a random value drawn from the normalized interval [0, 1].

= − ( ) /
(11) The photon packet's new weight, , is obtained further after sca ering events.is the ratio of the sca ering coefficient ( ) to the a enuation coefficient ( ) known as the albedo, 0 < < 1.The initial weight = 1, so with the increase of , its weight decays gradually, as shown in Equation ( 13) and Figure 6.When the photon's new position is determined and is still in front of the Reception Plane, it is still possible to be absorbed by the receiver.The sca ering angle and azimuth angle of the photon are reselected to obtain the new Direction Cosine [ , , ] in Equation ( 14), where = , and ѱ are the new sca ering angle and azimuth angle of the photon, respectively.is one of the moving photon's Direction Cosine, and its square must be more than 0. In addition, in the Equation ( 14), to ensure that 1 − is meaningful and rational, 1 − must be greater than 0, namely, 0 <

Termination of Photon Motion or Reception of Photon
In the simulation, the termination of a photon's motion is defined under two specific conditions.First, a photon is considered absorbed by the particles in the harbor water if its weight decreases to a threshold of 10 −6 or less.The choice of the minimum threshold will affect the calculation time of the simulation, and is also related to the link loss needed to add into the simulation, that is, the AL.The larger the link loss to be added, the smaller the threshold.Second, the photon's transmission is deemed to have ended once it reaches the receiver's photosensitive surface.
Throughout the photon tracking process, key a ributes of each photon are recorded, including its position, received angle, weight, and propagation distance.A photon is classified as 'Received' if it falls within the bounds of the receiver's photosensitive surface and if its received angle is within the specified FOV.If either of these conditions is not met, the photon is considered 'not Received', and its tracking is discontinued.
This methodical approach allows for precise determination of photon reception, facilitating an accurate assessment of the efficiency and effectiveness of the UWOC system in various harbor water conditions.

The Simulation Results and Discussion
In this paper, we primarily investigate the characteristics of UWOC with varying communication distances in four distinct types of harbor water.To facilitate this study, certain parameters of the transmi er and receiver are predetermined and held constant, with the main variable parameters being the transmi ed full divergence angle, the received aperture, and the received FOV.
The transmi ed full divergence angle significantly impacts the initial sca ering angle of photon, as described in Equation ( 9), and the corresponding beam waist radius influences the initial position of photons, as indicated by Equations ( 5) and ( 8).These factors dictate the varied paths photons traverse, thereby affecting the number and time of photons reaching the Reception Plane.Consequently, this variation directly influences the normalized power loss and time delay spread.The received aperture and FOV determine the receiver's capability to collect photons, with the different number of photons received resulting in variations in received optical power.Moreover, the disparate photon counts lead to discrepancies in the statistical timing of photons arriving at the Reception Plane, ultimately affecting both normalized power loss and time delay spread.
Reflecting the focus of our research, we have established a set of initial fixed parameters for use in our simulation model.These parameters, detailed in Table 2, include the laser wavelength, the number of photons transmi ed, the single photon energy, and the total transmi ed light energy.
This approach ensures a consistent baseline for analyzing the impact of the aforementioned variable parameters on the performance of UWOC systems in different harbor water conditions.Table 2. Initial fixed parameters utilized in the simulation.

Wavelength Transmitted Photon Number Single Photon Energy
Transmitted Light Energy 532 nm 10 13  3.74 × 10 −19 J 3.74 × 10 −6 J The specific se ings for the transmi ed full divergence angle, as outlined in Table 3, include values of 0.5 mrad, 5 mrad, 10 mrad, and 100 mrad.The Transmi ed half divergence angle ( ) of the Gaussian beam is defined as when the transmission distance → ∞, the angle between the position with 1/e times the center maximum amplitude of the Far-field spot with z axis, as shown in the Figure 7 below.The Transmi ed full divergence angle is 2 .The Far-field spot of the Gaussian beam is the spot determined by the position with 1/e times the center maximum amplitude, the spot radius ( ) in the propagation direction is shown in Equation (15), where is the beam waist radius.
So, the Transmi ed full divergence angle 2 can be derived: Corresponding to these full divergence angles, the Transmi ed beam waist diameter (2 ) for each se ing is calculated using Equation ( 16).In the analysis of how these varying transmi ed full divergence angles affect the communication system, both the received aperture and the received FOV are held constant.They are fixed at 200 mm for the received aperture and /2 rad for the received FOV.This methodical approach allows for a focused investigation into the impact of the transmi ed full divergence angle on UWOC system performance, ensuring that other variables do not confound the results.

rad rad
The specific se ings for the received FOV and aperture are also detailed in Table 3.The received aperture sizes are set at 25 mm, 50 mm, 100 mm, 200 mm, and 400 mm.The values for the received FOV are varied and include /90 rad, /30 rad, /22.5 rad, /15 rad, /7.5 rad, /2 rad, /1.5 rad, and rad.In the analysis that focuses on the impact of these variable received aperture and FOV on the communication system, the transmi ed full divergence angle is consistently fixed at 0.5 mrad.This approach enables a targeted examination of how changes in the received aperture and FOV affect the performance of the UWOC system, isolating these variables to understand their specific contributions.

Analytical Examination of the Full Divergence Angle in Transmission
In this study, we analyze the variation characteristics of laser power in relation to both the a enuation length and the divergence angle.Figure 8 illustrates the normalized power loss received as a function of the AL and the transmi ed full divergence angle.The AL is defined as the product of the a enuation coefficient ( ) and the communication distance , expressed mathematically as AL = ( ) × .
Simulation results indicate a consistent trend across all four types of harbor water: as AL increases, there is a corresponding increase in the a enuation loss of received power.Notably, the a enuation of received optical power initially occurs at a rapid rate and then stabilizes.This pa ern is a ributed to the gradual divergence of the light beam, which leads to a swift a enuation of light energy over a short distance.As the transmission distance extends, the light beam enters a state of full sca ering.In this phase, the proportion of energy a enuation due to light sca ering gradually diminishes, making absorption the predominant factor affecting light transmission.This trend is particularly evident in Figure 8a, highlighting the dynamic interplay between sca ering and absorption in determining the a enuation characteristics of UWOC.Furthermore, our analysis reveals that a larger transmi ed full divergence angle correlates with greater normalized optical power a enuation loss.This trend holds consistently across the four types of harbor water.
The normalized power loss of four harbor water in Figure 8 is different based on the same AL, which is inconsistent with the traditional Beer-Lambert law.This indicates that the Beer-Lambert law is completely inapplicable in harbor water due to the extreme scattering property, which is also confirmed by the reference [52], and the common approach in many studies of directly deriving communication distance by dividing the a enuation coefficients of two different water types proves to be inapplicable in the context of harbor water.The Monte Carlo algorithm is the valid solution for these highly sca ering communication scenarios.This finding underscores the complexity of the relationship between water quality, communication distance, and light a enuation in UWOC, challenging assumptions that may hold true in other types of aquatic environments.
In addition, with the same AL and divergence angle, it is shown that the normalized power loss is decreasing (↓) as the water quality deteriorates ( ↑) and the communication distance decreases ( ↓) in Figure 8.It can be concluded that in harbor water, the influence of communication distance on normalized power loss is more significant; in other words, the impact of water quality is overshadowed by the sheer magnitude of communication distance -induced losses.This phenomenon is a ributed to that intensified photons scattering in more deteriorating water might increase the diversity of photons reaching the Reception Plane, this diversification may offset some of the power loss, leading to a reduction in normalized power loss.
Our findings also indicate that the differences in normalized power loss among smaller divergence angles, specifically 0.5 mrad, 5 mrad, and 10 mrad, are relatively minor compared to a larger angle like 100 mrad.The primary reason for this is that a smaller transmi ed full divergence angle confines the light spot within a more limited range.This confinement results in reduced geometric loss and diminishes the likelihood of photons being sca ered by particles in the water.Consequently, a smaller divergence angle ensures that a higher number of photons can reach the photosensitive surface of the receiver and be effectively absorbed.This observation underscores the significance of divergence angle control in minimizing power loss in UWOC, particularly in environments with scattering particles.
Differences in the simulation results depicted in Figure 8a-d have been calculated, for example with a full divergence angle of 100 mrad and an AL of 30, the power loss of Harbor II, III, and IV are found to be 14.00 dB, 31.59dB, and 41.59 dB lower, respectively, than that of Harbor I. So, the normalized power loss of Harbor II, III, and IV is 3.98%, 0.069%, and 0.0069% times that of the Harbor I.In our next analysis, we examine the time delay spread in relation to both the AL and the transmi ed full divergence angle across the four types of harbor water, as depicted in Figure 9.We observe that under identical harbor water condition, the time delay spread increases with both the increase in AL and the transmi ed full divergence angle.Conversely, for a given AL and divergence angle, the time delay spread decreases as the quality of harbor water deteriorates.This trend mirrors the observed behavior in normalized power loss, indicating that the influence of communication distance on time delay spread is more significant than water quality; in other words, the impact of water quality is overshadowed by the substantial losses induced by the communication distance .We consider that increase in sca ering events may boost the incidence of forward sca ering, which could reduce the variability in the paths photons travel to the Reception Plane, thus resulting in a decrease in time delay spread.
Additionally, Figure 9d highlights a specific scenario where the AL ranges from 2-5 and the transmi ed full divergence angle is between 0.5 mrad and 5 mrad.In this case, the time delay spread is notably low, ranging only from 1.09 × 10 −11 -1.22 × 10 −11 s.This minimal spread can be a ributed to the relatively short transmission distance, approximately 0.45 m to 1.14 m, which is insufficient to contribute significantly to the time delay spread.This finding emphasizes the importance of considering both the transmission distance and beam divergence characteristics when assessing the time delay in UWOC.
The results in Figure 9a-d have been also numerically analyzed, the time delay spread of Harbor II, III, and IV is 30.56%,9.67%, and 0.49% times that of the Harbor I, when full divergence angle is 100 mrad, AL is 30.The differences of normalized power loss are much smaller than that of time delay spread.This indicates that the influence of absorption and sca ering on power loss is much greater than that of sca ering on time delay spread.

Analytical Evaluation of the Received Aperture and FOV
The study also examines how normalized power loss correlates with AL and the size of the Received aperture at the receiver's end, as illustrated in Figure 10.Our simulations reveal a consistent relationship between these variables across the four types of harbor water.Notably, as the size of the received aperture increases, the normalized power loss decreases.However, the rate of this decline slows and eventually stabilizes.This trend is a ributed to the fact that a larger received aperture can capture more sca ered photons, thus enhancing the absorption efficiency of the receiver.
Additionally, it is evident that larger AL values correspond to greater normalized power losses.This relationship underscores the significant impact of communication distance and beam a enuation characteristics on the efficiency of photon reception in UWOC systems.The relationship between normalized power loss, AL, and received FOV at the receiving end is further explored in Figure 11, focusing on the four types of harbor water.Simulation results indicate that the trends in normalized power loss correspond closely with both AL and the received FOV across these water types.In scenarios where the received FOV is relatively small (<0.2793 rad), there is a noticeable decrease in normalized power loss as the FOV increases.This decrease is a ributed to the fact that a larger FOV allows more photons, particularly those deviating from the optical axis, to become Ontarget photons and be absorbed by the receiver.However, as the received FOV continues to expand beyond this point, the impact on normalized power loss diminishes.This plateau in the effect is due to the large FOV eventually encompassing almost the entire peak intensity region of the Gaussian beam on the Reception Plane.Consequently, further increases in FOV do not significantly affect the proportion of photons captured, thereby stabilizing the normalized power loss.
In our final analysis, we examine the time delay spread at the receiving end in relation to AL, the size of the received aperture, and the received FOV across the four types of harbor water.This analysis is illustrated in Figures 12 and 13, and our findings indicate that, under constant harbor water conditions, the time delay spread increases with the rising AL.In the range of shorter ALs (2-10 AL), as shown in Figure 12, the time delay spread enlarges with an increase in the size of the received aperture.This trend is attributed to the fact that a larger aperture increases the likelihood of capturing multiplesca ered photons.However, in scenarios involving larger ALs (15-30 AL), the time delay spread exhibits minimal variation in response to changes in the received aperture size, remaining within the range of 1.05 × 10 −7 s to 7.39 × 10 −9 s.The primary reason for this is that at longer communication distances, most sca ered photons become Off-target photons due to multiple sca erings, and thus they are unlikely to be absorbed by the receiver.Consequently, in the range of 15-30 AL, the time delay spread remains relatively constant.In the case of a small received FOV (less than 0.8 rad), as depicted in Figure 13, there is a noticeable increase in the time delay spread as the FOV gradually expands.This pattern arises because a larger received FOV enhances the receiver's capability to capture a greater number of multiple-sca ered photons.As a result, the time delay spread grows with the increasing FOV.
However, beyond a certain point, the rate of increase in time delay spread significantly slows down.This slowdown occurs because a larger FOV eventually encompasses most of the multiple-sca ered photons within the Reception Plane.Therefore, further increases in the FOV do not substantially add to the number of photons captured, leading to only marginal changes in the time delay spread.This finding illustrates the intricate balance between FOV and the efficiency of photon capture in determining the time delay characteristics of UWOC systems.
Lastly, the simulation results and theoretical analysis collectively suggest that increasing the received aperture and FOV, particularly over longer communication distance, make li le contribution to reduce the normalized power loss and mitigate the time delay spread, thus not significantly enhancing communication quality.The primary reason is that, although a larger received aperture and FOV can capture more sca ered photons, beyond a certain threshold, these parameters expand to cover nearly the entire peak intensity region of the Gaussian beam or most of the multiple-sca ered photons on the Reception Plane.Based on these findings and the fixed transmi ed full divergence angle at 0.5 mrad, the most suitable received FOV range is about 1-3.2 rad, and the most ideal received aperture is about 0.15-0.4m.The variations in normalized power loss, under threshold conditions, predominantly register at less than 2 dB.Furthermore, statistical analyses concerning the time delay spread, also under these threshold conditions, reveal that it generally maintains the same order of magnitude, given the same AL.
To gain a thorough understanding of the UWOC system's performance in harbor water condition, we meticulously analyze the computational cost, a critical factor in our numerical simulations.The simulations are conducted on a computer equipped with an AMD Ryzen 9 5900X 12-core processor, operating at a main frequency of 3.7 GHz, and 64.0 GB of RAM, selected for its capability to efficiently process computationally intensive tasks.On average, the simulation process for each data point presented in

Conclusions
This study has critically examined the UWOC within four distinct types of harbor water, employing the Monte Carlo method to simulate photon behavior and analyze key communication parameters.Our findings reveal that the Transmi ed full divergence angle, received aperture, and FOV profoundly affect the UWOC system's performance.According to the simulation results, the following three main conclusions can be drawn:


Firstly, increased AL and divergence angle lead to higher received normalized power loss and time delay spread, the attenuation initially occurs at a rapid rate and then stabilizes, with these effects being more pronounced in clearer harbor water types. Secondly, this work reveals that the impact of communication distance on normalized power loss and time delay spread is more pronounced than that of water quality.This insight challenges that the common Beer-Lambert law no longer applies in harbor water.Communication distance, estimated by simply dividing the attenuation coefficients of different water types in many studies, is ineffective in harbor water contexts. Consequently, increasing the received aperture and FOV, particularly over longer communication distances, contributes minimally to reducing normalized power loss or mitigating time delay spread, thus not significantly enhancing communication quality.Based on the fixed transmitted full divergence angle, the most suitable received FOV range is about 1-3.2 rad, and the most ideal received aperture is about 0.15-0.4m.This result holds substantial technical relevance for the engineering design of UWOC across a range of harbor water conditions.
It should be noted that this work simplifies the absorption and sca ering models to facilitate analysis.We assume that the absorption and sca ering characteristics of harbor water remain consistent throughout the photon's travel path.However, real-world water bodies may display more complex sca ering and absorption behaviors due to dynamic water movements.This simplification represents a minor limitation of our study.
All these insights and analysis results pave the way for more efficient design and optimization of UWOC systems, particularly in challenging harbor water environments, by providing a detailed understanding of the interplay between light sca ering, absorption properties, and system configuration.Moreover, this study lays a technical and theoretical foundation that could facilitate the replication and expansion of our research into more turbid inland waters.Future research may develop more comprehensive models to accurately represent the diverse environmental conditions found in natural water bodies.

Figure 1 .
Figure 1.Schematic diagram of On-target photon and Off-target photon in underwater channel.

Figure 2
Figure 2 illustrates the multiple sca ering events that a photon emi ed from the light source undergoes before reaching the Reception Plane.The initial elevation angle and azimuth angle, denoted as and ѱ , respectively, correspond to those presented in Figure 1.The angle , , …represent the sca ering angles of the photon at each event, while ѱ , ѱ , …ѱ denote the azimuth angles associated with these events, with indicating the -th photon sca ering occurrence.A notable feature in this figure is the gray do ed lines, which depict the photon's motion coordinate system.This system undergoes continuous transformations with each sca ering event, reflecting the dynamic nature of the photon's path.The specifics of how this coordinate system is determined by the Equation (14).Additionally, it is important to highlight that the x, y, and z axes shown in Figure2are part of a fixed coordinate system.This system is established by the positioning of the transmi er and receiver, providing a stable reference frame for understanding the photon's motion and sca ering behavior.Additionally, the interaction between seawater particles and photons results in two distinct types of light: Forward-sca ering light and Backsca ering light[21,36], as shown in Figure3.In the simulation process, the receiver is positioned on the z-axis.The distance between the Transmission Plane and the Reception Plane is governed by a variable known

Figure 4 .
Figure 4. Schematic diagram of the initial position and orientation of the photon packet in the Transmission Plane.

Figure 5 .
Photons are still moving ?

Figure 6 .
Figure 6.Schematic illustration of changes in the photon packet's weight.

Figure 7 .
Figure 7. Schematic Transmi ed divergence angle and Far-field light spot.

Figure 8 .
Figure 8. Comparative analysis of normalized power loss relative to a enuation length and transmi ed full divergence angle across four types of harbor water.(a) Harbor I; (b) Harbor II; (c) Harbor III; (d) Harbor IV.

Figure 9 .
Figure 9. Variation of time delay spread in relation to a enuation length and transmi ed full divergence angle across four types of harbor water.(a) Harbor I; (b) Harbor II; (c) Harbor III; (d) Harbor IV.

Figure 10 .
Figure 10.Comparative analysis of normalized power loss relative to a enuation length and received aperture across four types of harbor water.(a) Harbor I; (b) Harbor II; (c) Harbor III; (d) Harbor IV.

Figure 11 .
Figure 11.Analysis of normalized power loss in relation to a enuation length and received FOV across four types of harbor water.(a) Harbor I; (b) Harbor II; (c) Harbor III; (d) Harbor IV.

FiguresFigure 12 .Figure 13 .
Figure 12.Assessment of time delay spread relative to a enuation length and received aperture across four types of harbor water.(a) Harbor I; (b) Harbor II; (c) Harbor III; (d) Harbor IV.

Table 3 .
Se ing values for the transmi ed full divergence angle, transmi ed beam waist diameter, received aperture and received FOV.