## Abstract

Using a recently proposed model for the refractive index fluctuations in oceanic turbulence, optical beam propagation through seawater is explored. The model provides an accurate depiction of the ocean through the inclusion of both temperature and salinity fluctuations to the refractive index. Several important statistical characteristics are explored including spatial coherence radius, angle-of-arrival fluctuations, and beam wander. Theoretical values of these parameters are found based on weak fluctuation theory using the Rytov method. The results presented serve as a foundation for the study of optical beam propagation in oceanic turbulence, which may provide an important support for further researches in applications for underwater communicating, imaging, and sensing systems.

### Keywords

- oceanic turbulence
- laser beam
- spatial coherence radius
- angle-of-arrival fluctuations
- beam wander

## 1. Introduction

The study of optical wave propagation through random media is a perpetually important topic for its many applications in the atmosphere and the ocean. Random fluctuations in the index of refraction cause beam spreading (beyond that due to pure diffraction), loss of spatial coherence, random wandering of the instantaneous beam center, and random fluctuations in the irradiance and phase [1]. The index of refraction fluctuations, generally referred to as optical turbulence, is one of the most significant quantities in optical wave propagation. For different random media, there are some differentiations among the index of refraction fluctuations. The index of refraction of atmosphere is primarily caused by fluctuating temperature. The refraction index in seawater is induced not only by temperature fluctuations but also by fluctuations of salinity. Changes in the optical signal due to absorption or scattering by molecules or particles are not considered here. Under the assumption of a statistically homogeneous and isotropic ocean, the power spectrum of oceanic turbulence is determined by fluctuations of refraction index.

With the development of underwater optical communications, imaging, sensor, and laser radar, it is indispensable to investigate the propagation behavior of laser beams through water medium. Knowledge of beam spreading is extremely important in a free space optics (FSO) communications link because it determines the loss of power at the receiver. The spatial coherence radius defines the effective receiver aperture size in a heterodyne detection system [1]. To the coherence degradation of laser beams, the spatial coherence radius can also be described as the strength of oceanic turbulence. Angle-of-arrival (AOA) fluctuations of an optical wave in the plane of the receiver aperture are associated with image jitter (dancing) in the focal plane of an imaging system [1] so that it plays a critical role in beam wave propagation applications such as imaging, lasercom, and other related areas. Movement of the short-term beam instantaneous center (or “hot spot”) is commonly called beam wander [1]. Beam wander is an important propagation characteristic of laser beams, which determines their utility for practical applications, such as laser communication [2, 3] and global quantum communication [4].

In this chapter, Section 2 describes a brief introduction of oceanic turbulence including the power spectrum and several significant oceanic parameters. The spatial coherence radius of a plane wave and a spherical wave propagating through oceanic turbulence has been investigated in Section 3, which are valid in both weak and strong fluctuations. Section 4 describes the angle-of-arrival fluctuations for plane- and spherical-wave models of oceanic turbulence. Based on the oceanic power spectrum, the beam wander effect with analytical and numerical methods in weak fluctuation theory is shown in Section 5. These results may provide an inroad for understanding laser beam propagation through oceanic turbulence, and the theoretical findings may provide an important support for further researches in applications for underwater communicating, imaging, and sensing systems.

## 2. Nature of oceanic turbulence

Turbulence is a random, three-dimensional motion with the velocity and vorticity irregularly distributed in time and space [5]. In general, turbulence is accepted to be an energetic, rotational, and eddying state of motion that results in the dispersion of material and the transfer of momentum, heat, and solutes at rates far higher than those of molecular processes alone [6]. It is characterized by an energy transfer from large to small scales where the dissipation of kinetic energy is taking place [5]. Oceanic motions are constrained to flow along density surfaces by the Earth’s rotation and the density stratification. In the upper ocean, microscale turbulence is generated by surface winds, air-sea cooling, or evaporation. In the ocean interior, microscale turbulence develops when internal waves develop strong shears and overturn and break, much like surface gravity waves [7]. These breaking events play a fundamental role in the ocean circulation, because they mix the densest waters at the ocean bottom with the lighter waters above, thereby allowing the densest waters to come back to the surface [7]. Much of the turbulence induced in benthic boundary layers is driven by external processes resulting from the fluxes of buoyancy and momentum through the nearby boundary, such as a tidally driven current, geothermal heat flux, and so on [6]. They are driven by sources of energy outside the benthic boundary layer itself [6].

### 2.1. Power spectrum of oceanic turbulence

Since the power spectrum of oceanic turbulence proposed in 2000 [8], there has been remarkable interest in the study of propagation characteristics using laser beams in seawater. The power spectrum of oceanic turbulence has been simplified for homogeneous and isotropic water media [9], which is applicable for isothermal water [10]. When the eddy thermal diffusivity and the diffusion of salt are assumed to be equal, the power spectrum for homogeneous and isotropic oceanic water is given by the expression [11].

where

### 2.2. Oceanic parameters

In this section, the abovementioned important parameters should be presented in detail that will benefit to more accurately comprehend the oceanic turbulence. In particular, the four significant parameters, such as the rate of dissipation of kinetic energy per unit mass of fluid, the rate of dissipation of mean-squared temperature, the Kolmogorov microscale, and the ratio of temperature and salinity contributions to the refractive index spectrum, will be mainly involved in the following subsections.

#### 2.2.1. The rate of dissipation of kinetic energy per unit mass

The rate of dissipation of the kinetic energy of the turbulent motion per unit mass of fluid through viscosity to heat is usually denoted by

where * v*is the kinematic viscosity, the tensor

*is given by*s

_{ij}

Recorded

#### 2.2.2. The rate of dissipation of mean-squared temperature

The effect of turbulence on the fluid temperature field can be described as the rate of dissipation of mean-squared temperature [6],

where

In isotropic turbulence when the mean-squared gradients of temperature are the same in all directions, so that the rate of dissipation of mean-squared temperature becomes [6].

It is noted that the rate of dissipation of mean-squared salinity,

#### 2.2.3. Kolmogorov microscale

The turbulent flow contains eddies of various sizes, and the energy is transferred from larger eddies to smaller eddies until it is drained out by viscous dissipation. Kolmogorov’s hypothesis asserts that for large Reynolds numbers (i.e., inertial subrange), the small-scale structure of turbulence is statistically steady, isotropic, and locally homogeneous, and independent of the detailed structure of the large-scale components of turbulence [12]. Kolmogorov microscale is the smallest scale in turbulent flow. At the Kolmogorov scale, viscosity dominates and the turbulent kinetic energy is dissipated into heat. The length scale of the turbulent motions at which viscous dissipation becomes important must depend on factors that provide measures of the turbulent motion and of its viscous dissipation [6]. Kolmogorov length scale [13].

where a range of ^{−5} m in very turbulent regions to 0.01 m in the abyssal ocean [6].

#### 2.2.4. The ratio of temperature and salinity contributions to the refractive index spectrum

The parameter * w*is the ratio of temperature and salinity contributions to the refractive index spectrum given by

where

## 3. Spatial coherence radius

When one coherent optical wave propagates through a random medium, various eddies impress a spatial phase fluctuation on the wave front with an imprint of the scale size [1]. The accumulation of such fluctuations on the phase leads to a reduction in the “smoothness” of the wave front [1]. Hence, turbulent eddies further away experience a smoothness of the wave front only on the order of the transverse spatial coherence radius, which Andrews and Phillips denote by

### 3.1. Plane wave

Under Rytov approximation, the wave structure function (WSF) of a plane wave propagating through isotropic and homogeneous turbulence is defined by [1].

where * k*is the optical wave number related to the wavelength

*is the path length,*L

By expanding the zero-order Bessel function in power series, the WSF is written in the form

where the power spectrum given by Eq. (1) is used and the order of summation and integration is interchanged. In addition,

Based on the properties of hypergeometric function and Pochhammer symbol [14], after very tedious calculations [10], the WSF of a plane wave in certain asymptotic regimes is

The separation distance at which the modulus of the complex degree of coherence (DOC) falls to

### 3.2. Spherical wave

Under Rytov approximation, the WSF of a spherical wave is defined by [1].

Similarly, the WSF of a spherical wave is derived in [10].

and the spherical-wave spatial coherence radius as [10].

### 3.3. Discussions

Based on the formula of Eqs. (10), (11), (13), and (14), the WSF of both a plane wave and a spherical wave can be written as

Equation (15) indicates that the spatial coherence radius is the only parameter characterizing the WSF, and under Rytov approximation, the Kolmogorov five-thirds power law of wave structure function is valid for oceanic turbulence in the inertial range if the power spectrum of oceanic turbulence proposed by Nikishov is adopted.

According to Ref. [1], under Rytov approximation, the definitions of wave structure function of a plane wave and a spherical wave are given by Eqs. (8) and (12), respectively. It is known that the expression for wave structure function depends on the mutual coherence function. Rytov approximation is limited to weak fluctuations. However, for the special cases of a plane wave and a spherical wave, it has been shown that mutual coherence function derived by strong fluctuation theories is the same as that derived by Rytov approximation [1]. Only a plane-wave and a spherical-wave case are considered in this section. Thus, the results of the wave structure function and the spatial coherence radius obtained are valid in both weak and strong fluctuations.

## 4. Angle-of-arrival fluctuations

Angle-of-arrival (AOA) fluctuations play an important role in a diverse range of fields including atmospheric turbulence [15, 16], free space optical communication [17], ground-based astronomical observations [18], and so on.

Angle-of-arrival fluctuations of an optical wave in the plane of the receiver aperture are associated with image dancing in the focal plane of an imaging system. Fluctuations in the AOA can be described in terms of the phase structure function [1]. In order to understand it easily, let

Under the geometrical optics method, the AOA is defined by [19].

Based on the homogeneous and isotropic oceanic turbulence, the mean

where

### 4.1. Angle-of-arrival fluctuations of plane wave

The phase structure function associated with an unbounded plane wave is given by [1].

where the normalized distance variable

Based on Eqs. (18) and (19), the AOA fluctuations for a plane wave can be expressed as [1].

the variance of AOA in the

By using

Substituting Eq. (21) into Eq. (20), the analytical expression of AOA fluctuations for a plane wave is

To clarify the physical explanation, we introduce the plane-wave spatial coherence radius

where

### 4.2. Angle-of-arrival fluctuations of a spherical wave

In the case of a spherical wave, the phase structure function is defined by [1].

Based on Eq. (25), the AOA fluctuations for a spherical wave can be written as [1].

Similarly, the phase structure function of a spherical wave is expressed as

and the analytical expression of AOA fluctuations for a spherical wave is

Substituting the spherical-wave spatial coherence radius

where

### 4.3. Discussions

As mentioned in Section 1, the spatial coherence radius plays an important role in a heterodyne detection system. To the best of our knowledge, it is interesting to research the relation between AOA fluctuations and the spatial coherence radius, because spatial coherence radius

## 5. Beam wander

Movement of the short-term beam instantaneous center (or “hot spot”) is commonly called beam wander [1]. This phenomenon can be characterized statistically by the variance of the hot spot displacement along an axis or by the variance of the magnitude of the hot spot displacement [1]. An estimate of the short-term beam radius is obtained by removing beam wander effects from the long-term beam radius [1]. It is much convenient to use the geometrical optics approximation method in the turbulent area. Beam wander is an important characteristic of laser beams, which determines their utility for practical applications, such as ground-to-satellite laser communication [2, 3] and global quantum communication [4].

### 5.1. A general model

The far-field angular spread of a free-space propagating beam of diameter

Based on the introduction of a general model [1], beam wander can be expressed as

where bracket < > denotes an ensemble average,

The large-scale filter function is [1].

where

Because beam wander is caused mostly by a large-scale turbulence near the transmitter, the last term can be dropped in Eq. (32) and the geometrical optics approximation is [1].

Substituting from Eqs. (1), (32), and (33) into Eq. (31), (31) leads to

Equation (34) is applicable for collimated, divergent, or focused Gaussian-beam waves, and it can represent our general expression for the variance of beam wander displacement under weak irradiance fluctuations.

### 5.2. Special cases

In this section, two special cases (i.e., collimated beam and focused beam) are analyzed.

For collimated beam (

where

and

For focused beam (

where

and

To atmospheric turbulence, the focused beam case leads to a greater beam wander variance for the same size beam at the transmitter as that for the collimated beam [1]. However, because of the complexity of oceanic power spectrum, the analytical expressions of collimated beam and focused beam are also less concise than those of atmospheric turbulence. Therefore, it is not simple to distinguish whose variance of beam wander is the larger one directly. In Section 5.4, numerical calculations are used to discuss the abovementioned property of different beams.

### 5.3. Dimensionless quantity B W

In order to obtain influences of beam wander on laser beam propagation through oceanic turbulence, the relation between beam wander and the turbulence-induced beam spot size is investigated in detail by using theoretical and numerical methods. Using the dimensionless quantity

### 5.4. Definition of relative beam wander

To obtain the difference of beam wander among various beam types, the relative beam wander

It is clear that the increment

From Section 5.3, it is impossible to avoid the beam wander effect for laser propagation; therefore, achieving small beam wander is imperative. In this section, the relative beam wander describes the increment of beam wander between focused and collimated beam, and this quantity benefits us to select how to obtain small value of beam wander. Based on Eq. (39), it is feasible to obtain small beam wander as long as the arbitrary beam type is under

## 6. Conclusions

In summary, this chapter has used the same idea to convert atmospheric turbulence concerning beam wave propagation to corresponding oceanic turbulence. The general characteristics of an optical wave propagating through the ocean are greatly affected by small fluctuations in the refractive index that are the direct consequence of small temperature and salinity fluctuations transported by the turbulent motion of the ocean. Therefore, the propagation process of one optical wave suffers beam spreading, loss of spatial coherence, angle-of-arrival fluctuations, beam wander, and so on. In this chapter, three important statistical quantities including spatial coherence radius, angle-of-arrival fluctuations, and beam wander have been investigated.

The analytical formulae for the wave structure function and the spatial coherence radius of a plane wave and a spherical wave propagating through oceanic turbulence have been derived in Section 3, which are valid in both weak and strong fluctuations. It has been shown that under Rytov approximation, the Kolmogorov five-thirds power law of wave structure function is also valid for the oceanic turbulence in the inertial range if the power spectrum of oceanic turbulence proposed by Nikishov is adopted. These results are of considerable theoretical and practical interest for operations in communication, imaging, and sensing systems involving turbulent underwater channels.

Furthermore, spatial coherence radius can be described as the coherence degradation of laser beams propagating through ocean induced by the strength of turbulence; thus, the relation between angle-of-arrival fluctuations and the spatial coherence radius has been researched. Both for a plane wave and for a spherical wave, it is shown that the angle-of-arrival fluctuations are inversely proportional to five-thirds order of spatial coherence radius in the inertial range. In terms of the influences of spatial coherence radius, the angle-of-arrival fluctuations of plane and spherical waves in oceanic turbulence have the similar behavior to that of atmospheric turbulence.

In addition, based on the oceanic power spectrum, the beam wander effect has been studied with analytical and numerical methods in weak fluctuation theory, and the analytical expressions for beam wander of collimated and focused beam in oceanic turbulence have also been derived. For the dimensionless quantity,

In this chapter, the classical treatments of optical wave propagation have been concerned with part of special cases, such as uniform plane wave, spherical wave, collimated beam, and focused beam. The results presented serve as a foundation for the study of optical beam propagation in oceanic turbulence, which may provide an essential support for further researches in applications for underwater communicating, imaging, and sensing systems. Thus, these simple optical wave models are useful in describing certain aspects of wave propagation in oceanic turbulence. However, due to inherent infinite extent, these models are not adequate in describing laser beams when finite size of the transmitted wave and high-order Gaussian wave should be taken into account in the near future.

## Acknowledgments

Lu Lu and Zhiqiang Wang wrote this manuscript. Lu Lu, Zhiqiang Wang, and Chengyu Fan arranged the structure of this manuscript. Lu Lu and Xiaoling Ji performed theoretical calculations and physical analysis in Section 2. Lu Lu, Zhiqiang Wang, Pengfei Zhang, Chunhong Qiao, Jinghui Zhang, and Chengyu Fan performed Sections 3–5. Lu Lu deduced the analytical expressions. Zhiqiang Wang and Pengfei Zhang performed the numerical simulations and analyzed data. Pengfei Zhang, Jinghui Zhang, Chunhong Qiao, and Chengyu Fan revised this manuscript. All authors contributed this work equally and acknowledged the support by the National Natural Science Foundation of China (NSFC) under grants 61405205 and 61475105.

## Conflict of interest

The authors declare no conflict of interest.

## Notes

All authors are very much thankful to the valuable suggestions in the spatial coherence radius section by Prof. Yahya Baykal.