## Abstract

We analyzed the observations of scintillations in a laser beam (532 nm, $\sim 200\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{mW}$ power) traveling along a 144 km path at an altitude of 2.2–2.4 km above sea level, just above the atmospheric boundary layer, between the islands of La Palma and Tenerife. The observations were performed during nighttime on 18 July 2011, by means of a telescope with an aperture diameter of 1 m. Strong scintillations were observed. The estimates of spatial spectra and correlation functions indicated that the observed intensity fields possess, statistically, a locally isotropic structure, which agrees with the idea of a locally isotropic turbulence. The estimates of spatial autospectra and autocorrelation functions of the intensity field indicated that the characteristic scale of the internal structure of the observed clusters is 6.5–8 mm, while the characteristic size of the clusters is 4–5 cm. The major contribution to the observed scintillations comes from the inhomogeneities of the intensity field with scales from 1–2 cm up to 10–12 cm. The analysis of the cross-spectra indicated that the hypothesis of frozen turbulence introduced by Taylor can be used for the description of spatiotemporal structure of intensity fluctuations of laser beams traveling through long paths in the atmosphere.

© 2012 Optical Society of America

## 1. Introduction

Inhomogeneities of the atmospheric refractivity caused by turbulence destroy the spatial coherence of light propagating through the atmosphere. This phenomenon manifests itself, for example, in stellar scintillations that can be observed by the naked eye. In the optical range, the inhomogeneities responsible for scintillations, to a high degree of accuracy, are proportional to the temperature fluctuations of air. A model of the spatial structure of temperature fluctuations was developed by Obukhov and Corsin [1,2]. The Obukhov–Corsin model is based on the elaboration of the fundamental principles suggested by Kolmogorov for the description of developed, locally homogeneous turbulence. The spatial spectrum of the fluctuations that follows from this model is commonly referred to as the Kolmogorov spectrum.

The description of the temporal evolution of the temperature fluctuations at a fixed spatial point relies upon the hypothesis of the frozen turbulence introduced by J. Taylor [3] for the interpretation of measurements of turbulence in wind tunnels. Based on the fact that the air velocity fluctuations in such a tunnel are significantly smaller than the average velocity of the air flow, the hypothesis states that the statistics of turbulent fluctuations in the coordinate frame moving with the average velocity of the flow do not change with time.

In atmospheric turbulence, however, the relative fluctuations of the velocity are significantly larger than those in wind tunnels. In order to verify the applicability of the Taylor hypothesis to the atmosphere, a comparison was performed between the spectra of temperature fluctuations measured at a 70 m high tower and those obtained at an aircraft flying in the vicinity of the tower [4]. The wind velocity was $9\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}/\mathrm{s}$, while the air speed of the aircraft was 55–$88\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}/\mathrm{s}$. The comparison corroborated the applicability of the hypothesis of frozen turbulence to a high accuracy. Furthermore, experimental studies of the spatiotemporal structure of fluctuations in the surface air [5] indicated that the evolution of inhomogeneities of air temperature and density, drifted by the wind, is small during the time of their displacement over a distance of the order of the scale of the inhomogeneities.

The hypothesis of frozen turbulence forms the basis of the description of the light intensity fluctuations caused by turbulence, as a function of time [6]. Diffraction effects in combination with the frozen turbulence hypothesis [6] explain, in particular, the “running shadows” at the input aperture of a telescope directed at a bright star. This was observed by astronomers as early as in the 17th century, when they removed the ocular and focused their eye on the aperture. In general, the validity of the hypothesis is undoubted for intensity fluctuations in a plane wave propagated through a thin layer of a turbulent medium. For narrow laser beams and large propagation distances, however, where the turbulent spread of the beam exceeds its diffractive divergence, its applicability is far from being that obvious. In a recent study [7], discussing the measurements of the fluctuations of a laser beam over a 149 km long atmospheric path, it was concluded that: “These results challenge the assumption of ‘frozen’ turbulence (Taylor hypothesis) commonly used in analyses.”

Satellite-based observations of stars through the atmosphere in limb sounding geometry [8,9] indicate that the fluctuations of a laser beam along a limb path would become noticeable at perigee heights of 65–75 km. At perigee heights below 25–30 km, the fluctuations will be strong. For the description of strong fluctuations, the Markov approximation [10,11] was introduced. In this approximation, however, the equations that link the observed fluctuations with the atmospheric turbulence are still rather complicated. Their solutions have not yet been investigated sufficiently for the use of this approximation in computer-based modeling of the laser beam propagation along a limb path. A detailed experimental investigation is required to confirm the validity of this approximation.

One of the goals of an experimental campaign performed in the Canary Islands in July 2011, within the framework of a project on long-path greenhouse gas measurements [12] for basic demonstration of the new concept of infrared-laser occultation [13], was a better understanding of the spatiotemporal structure of intensity fluctuations caused by atmospheric turbulence within the observation plane after long-path laser beam propagation. A 144 km long optical path between the Northern Optical Telescope (NOT) observatory on La Palma and the Optical Ground Station (OGS) observatory on Tenerife was employed for the purpose [12]. A similar path has been used in an earlier study [14] of wave front fluctuations. The velocity of the transversal drift along this path with fixed transmitter and receiver locations is two to three orders of magnitude smaller than the satellite velocities in space-borne laser occultation [13]. Therefore, in order to best learn from the experimental results also for satellite observations, it was important to study the spatial spectra and correlation functions of the intensity and test the applicability of the Taylor hypothesis to the analysis of laser beam fluctuations. This study is discussed in this paper.

## 2. Observations of Intensity Fluctuations

The laser beam investigated in this study was emitted by a continuous wave, frequency-doubled, Nd:YAG laser (200 mW, 532 nm) located at the NOT observatory on La Palma (28.76° N, 17.89° W, 2392 m MSL altitude). The outgoing beam was collimated up to a diameter of 15 cm and directed with narrow beam divergence at the OGS telescope, located on Tenerife (28.30° N, 16.51° W, 2410 m MSL altitude), which has an aperture diameter of 1.016 m [12]. Within the OGS telescope optical assembly, the path distance to the coudé focus was 39 m. The optical path length between the transmitter and receiver was 144 km. The major part of the optical path was located above the ocean. The ray perigee height, taking the regular refraction into account, is about 2.2 km altitude. The ray path is, therefore, above the atmospheric boundary layer residing below 2 km altitude. In the absence of turbulence, due to diffraction effects only [15] and without intentional beam divergence, the beam size at the aperture of the OGS telescope would have a diameter of about 40 cm. The turbulent spread of the beam, however, significantly exceeded the diffractive divergence, and the diameter of the spot on the dome, where the telescope was mounted, significantly exceeded the aperture size due to the beam divergence [12]. Turbulence broke the beam into many separate clusters. Because of the vision inertia, the pattern of “running shadow” moving approximately in the horizontal direction could be observed with the naked eye. For a fast registration of the intensity field ($I$) of the received laser radiation, we employed a Prosilica GE-680 CCD camera mounted on an optical bench in the coudé focus room. The camera was equipped with a Fujinon HF75SA-1 lens with a focal length of 75 mm, which produced an image of the aperture of the OGS telescope onto the CCD matrix. The camera produced film records in the AVI format with a depth of 8 bits per pixel at a frame rate of 205 frames per second. The dynamic range limitation determined by the bit depth resulted in a large number of under-/overexposed records that were rejected during the processing of the observations. The correction of the spectra distortion imposed by the finite exposure time is discussed below.

The observations were performed at nighttime in order to exclude any background scattered radiation from the sun. Figure 1 shows an example of the consecutive distributions of intensity $I(\mathit{\rho},t)$ acquired from the aperture of the telescope using an acquisition time of about 5 ms, where $\mathit{\rho}=(x,y)$ is the 2D vector of coordinates in the aperture plane, and $t$ is the time.

It should be noted that the units of intensity are dimensionless throughout this paper. The images are vignetted by parts of the optical system assembly and limited by the visual field of the CCD camera lens. The dark spot in the middle of the image is the shadow of the secondary mirror, which has a diameter of 0.327 m. The knowledge of this diameter allows to transform the discrete coordinates of $\mathit{\rho}$, initially specified by the pixel numbers, into the geometrical coordinates in the aperture plane.

The intensity field is broken into chaotic clusters whose positions change from frame to frame. Upon inspection of Fig. 1, however, one cannot see any noticeable ordered movement associated with “running shadows”. Figure 2 shows an example of light flux $J(t)={\iint}_{S}I(\mathit{\rho},t){\mathrm{d}}^{2}\mathit{\rho}$ through the OGS aperture as a function of time. The observed fluctuations are primarily caused by the fluctuations of the atmospheric refractivity. Slow variations of attenuation from diffusion or absorption are negligible.

Despite significant spatial averaging over the aperture, which has a diameter of 1 m, the light flux indicates strong fluctuations with time. The ratio of the root mean square (rms) deviation $\sqrt{\overline{{J}^{2}}-{\overline{J}}^{2}}$ and the mean $\overline{J}$ estimated over the time interval of 120 s is 0.56. The probability distribution of $J$ differs from the normal distribution because the light flux has short positive spikes. The fluctuations of $J(t)$ are caused by both turbulent displacements of the beam as a whole (beam wandering) and the beam spread.

Figure 3 shows the estimates of the mean intensity $\overline{I(\mathit{\rho},t)}$ (left panel) in the telescope aperture and rms fluctuations of the intensity divided by the mean $\sqrt{\overline{I(\mathit{\rho},t{)}^{2}}-{\overline{I(\mathit{\rho},t)}}^{2}}/\overline{I(\mathit{\rho},t)}$ (right panel). The mean values are evaluated from a representative sequence of 22,542 frames from the two records taken on 18 July 2011, within 00:05:46 to 00:08:23 [$\mathrm{UTC}+1\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{h}$] (termed study time interval hereafter); the accumulation time was thus about 120 s. One pixel in the CCD matrix corresponds to ${l}_{P}=1.18\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{mm}$ in the aperture, which defines the Nyquist frequency of ${\kappa}_{\mathrm{max}}=\pi /{l}_{P}=2.65\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{rad}/\mathrm{mm}$. The white dashed circle indicates the border of the secondary mirror shadow. The coordinate axes are parallel to the sides of the CCD matrix. The white arrow indicates the horizontal direction. This rotation of the image is determined by the altitude-azimuth mounting of the OGS telescope. Throughout this paper, we will use the coordinate system rotated in such a way that the $x$ axis is aligned horizontally. The shadows of the mounting rods of the secondary mirror are clearly visible. The temporal averaging damps the turbulent fluctuations in the intensity, and the measured estimates of its mean values indicate only insignificant fluctuations in comparison with the values of an individual measurement (i.e., a single recorded frame) shown in Fig. 1.

The ratio of the rms deviations and mean values in the right panel is mostly within a range of 1.5–1.7, except narrow bands around the edges and faulty areas of the CCD matrix due, most likely, to particles of dust on the surface of the matrix. Large rms variations are typical for strong scintillations [14,16–18].

In the left panel of Fig. 3, we selected three rectangular regions, 1, 2, and 3, respectively. For each rectangle, we estimated the mean and rms intensities over time, which were then averaged over all the rectangles. The ratio of the rms and the mean fluctuates within a range of 1.59–1.62, while the mean values have a range of 26.38–28.80. This allows for treating the observed fluctuations of the intensity in the approximation of a statistically homogeneous field using the concepts of mean values and spectral densities [19]. The sizes of ${l}_{x,y}$ along the $x$ and $y$ axes of the selected areas used for the analysis define the resolution $\mathrm{\Delta}{\kappa}_{x,y}=2\pi /{l}_{x,y}$ of the spatial spectrum.

## 3. Spectra and Correlation Functions of the Intensity Fluctuations in the Telescope Aperture

Observed distributions of intensity $I(\mathit{\rho};t)$ in the aperture of the OGS telescope are described by the following expression:

By substituting into Eq. (1) the expression of ${I}^{(0)}$ from the Fourier integral, we arrive at the following equation:

Estimates of ${F}^{(I)}(\mathit{\kappa},0)$ were evaluated for the three rectangular regions shown in the right-hand panel of Fig. 3. To this end, for each frame with number $g$ (moment of time ${t}_{g}$) we evaluated, using the Hann window, the discrete Fourier transform ${f}^{(I)}(\mathit{\kappa},{t}_{g})$ of the intensity $I(\mathit{\rho},{t}_{g})$. Values of ${|{f}^{(I)}(\mathit{\kappa},{t}_{g})|}^{2}$ were averaged over a series of frames. As a result, we obtained the following estimates of the autospectrum for each rectangle:

where $G$ is the number of frames in the series. In these results, $G=\mathrm{22}\mathrm{,}\mathrm{542}$. The autospectra were used to estimate the autocorrelation functions [19] ${B}^{(I)}(\delta \mathit{\rho})=\u3008I(\mathit{\rho},t)I(\mathit{\rho}+\delta \mathit{\rho},t)\u3009$ of the intensity fluctuations in the telescope aperture:In order to determine the applicability of the hypothesis of frozen turbulence for the analysis of scintillations, we considered the cross-spectral densities ${F}^{(I)}(\mathit{\kappa};\delta t)$, $\delta t\ne 0$:

The cross-spectral density is, generally speaking, a complex value being equal to $|{F}^{(I)}(\mathit{\kappa},\delta t)|\mathrm{exp}(i{\phi}^{(I)}(\mathit{\kappa},\delta t))$, where ${\phi}^{(I)}(\mathit{\kappa},\delta t)$ is known as the phase spectrum, $-\pi \le {\phi}^{(I)}<\pi $. The coherency spectrum $\mathrm{Coh}(\mathit{\kappa},\delta t)$ is defined as follows:

If the hypothesis of frozen turbulence $I(\mathit{\rho},t+\delta t)=I(\mathit{\rho}-\mathit{u}\delta t,t)$ holds for arbitrary $\mathit{\rho}$ and $t$, the following equation must also hold:

## 4. Estimates of Spectra and Correlation Functions from Observations of Scintillations

Figure 4 presents the estimates of the 2D spectra of fluctuations, normalized to the dispersion, for the three rectangles shown in the left panel of Fig. 3. We estimated the dispersion of scintillations by using Parseval’s formula:

For wave numbers $\kappa >1.0-1.2\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{rad}/\mathrm{mm}$, the estimate of the power spectral density (PSD) is unobtainable due to noise. The noise becomes noticeable for small values of PSD around 0.1. The isolines are circular in nature, slightly oblate in the direction of axis ${\kappa}_{x}$. For all three rectangles, the isolines are close to each other, while the PSD varies by three orders of magnitude. This corroborates the assumption of the statistical homogeneity of the intensity fluctuations in the telescope aperture plane.

Figure 5 shows the estimate of the normalized autocorrelation function ${B}^{(I)}(\delta \mathit{\rho})/{B}^{(I)}(0)$ for rectangle 1. Autocorrelations for rectangles 2 and 3 (not shown) are practically indistinguishable from those for rectangle 1.

## 5. Discussion of Spectra and Correlation Functions

Atmospheric turbulence causes fluctuations in the intensity of a laser beam propagating through the atmosphere [6]. Because of local isotropy of turbulence, the 2D PSD isolines of the intensity field in the perpendicular cross-section are nearly circular. Observations along a near-surface path [20] indicated that turbulent fluctuations of density and temperature are isotropic. Those observations were performed with a wind velocity of several $\mathrm{m}/\mathrm{s}$, the exposure time in photographing the intensity field was $2\tau =0.1\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{ms}$, i.e., 50 times smaller than in the observations discussed in this paper. Therefore, the factor distorting the isotropy of the right-hand part of Eq. (4), which equals $\u3008{(\mathrm{sinc}(\mathit{\kappa}\mathit{u}\tau ))}^{2}\u3009$, was close to unity [20].

The 2D PSD ${F}^{(I)}(\mathit{\kappa},0)$ presented in Fig. 4 rapidly decreases in power with the increase of the wave number $\kappa $. The range of the variation of PSD is more than three orders of magnitude. To estimate the characteristic scale of the inhomogeneities making the major contribution to the dispersion of scintillation, we employed the fact that ${F}^{(I)}(\mathit{\kappa},0)$ can be approximately represented as a function of $\kappa =\sqrt{{\kappa}_{x}^{2}+{\kappa}_{y}^{2}}$. Changing to the polar coordinates in the wave number plane $\mathit{\kappa}$ and integrating ${F}^{(I)}(\mathit{\kappa},0)$ over the angular variable, we arrive at the following expression for the dispersion, ${\sigma}_{I}^{2}$:

Changing to the integration over $\mathrm{ln}\text{\hspace{0.17em}}\kappa $ allows for a detailed study of the contribution of small wave numbers in the dispersion. The dimensionless product ${\kappa}^{2}{V}^{(I)}(\kappa )$ is the density of the contribution of the spectral component to the dispersion.

Figure 6 presents this density, ${\kappa}^{2}{V}^{(I)}(\kappa )$, smoothed with a Gaussian window of width $0.004\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{rad}/\mathrm{mm}$, as a scaled function of $\mathrm{\Lambda}=2\pi /\kappa $. The upper limit of the scale is defined by the sizes of the rectangles shown in Fig. 3. The area under the curves is proportional to the dispersion of scintillations. The maximum of dependence ${\kappa}^{2}{V}^{(I)}(\kappa )$ corresponds to inhomogeneities with scales of about $\mathrm{\Lambda}=3-4\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{cm}$. The major contribution to the dispersion of scintillations comes from clusters with scales from 1–2 to 10–12 cm. These are clearly visible in Fig. 1.

To estimate the scale of correlation of scintillation in the telescope aperture, we utilized the fact that the cross-sections of the surface $b(\delta \rho )={B}^{(I)}(\delta \rho )/{B}^{(I)}(0)$, $\delta \rho =\sqrt{\delta {x}^{2}+\delta {y}^{2}}$ by planes ${B}^{(I)}(\delta \rho )=\text{const}$ are almost circular. Figure 7 presents functions $b(\delta \rho )$ for the three rectangles shown in Fig. 3.

The correlation coefficient decreases down to 0.5 at $\delta \rho \approx 7\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{cm}$, which provides an estimate of the characteristic internal scale of the intensity clusters. At a distance of $\delta \rho \approx 45-50\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{cm}$, the correlation coefficient changes its sign, which provides an estimate of the characteristic external scale of the clusters. Because of the law of conservation of energy, the full integral of the correlation function $\int b(\delta \rho )\delta {\rho}^{2}\mathit{d}\delta \rho $ must equal 0; therefore $b(\delta \rho )$ must be an alternating-sign function [6]. Turbulence results in a redistribution of the intensity in space, the formation of clusters must be accompanied by the appearance of low intensity areas, which is clearly visible in Fig. 1.

The 2D PSD ${F}^{(I)}(\mathit{\kappa},0)$ shown in Fig. 4 and autocorrelation function ${B}^{(I)}(\delta \mathit{\rho})/{B}^{(I)}(0)$ shown in Fig. 5 characterize the intensity fluctuations in space at an arbitrary moment of time, under the assumption of stationary scintillations. The coherence spectra shown in Fig. 8 describe the evolution of intensity clusters in time. Their comparison with autospectra ($\delta t=0$) presented in Fig. 4 indicates that during the time interval of 5 ms, the coherence spectrum significantly decreases, especially for large wave numbers $\kappa >0.4\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{rad}\text{\hspace{0.17em}}{\mathrm{mm}}^{-1}$. Therefore, the correlation of intensity fluctuations decreases with time faster for smaller scale inhomogeneities of the intensity field. An intensity cluster in the telescope aperture plane not only drifts but also changes shape. This is noticeable when examining the distribution of intensity of individual frames, as presented in Fig. 1.

Figure 9 shows the 2D phase spectra. A comparison with Fig. 8 indicates that the phase spectrum can be estimated in the area of wave numbers $({\kappa}_{x},{\kappa}_{y})$ where the coherence spectrum is not too small. Outside this area, the uncertainty of the estimate of complex numbers ${F}^{(I)}(\mathit{\kappa},\delta t)$ [19] does not allow for a reliable estimate of their argument $\phi (\mathit{\kappa},\delta t)$. Isolines here look like random trajectories. Because we defined the phase in the interval from $-\pi $ to $\pi $, its linear increase is equivalent to a saw-tooth function with jumps of $2\pi $. Such jumps are observed at ${\kappa}_{x}=\pm 0.3\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{rad}\text{\hspace{0.17em}}{\mathrm{mm}}^{-1}$ in the left panel and at ${\kappa}_{x}=\pm 0.15$ in the right panel. In order to better visualize the linear dependence of $\phi (\mathit{\kappa},\delta t)$ from ${\kappa}_{x}$ between the jumps, the phase spectra were averaged over ${\kappa}_{y}$ in the region of $-0.425\le {\kappa}_{y}\le 0.425\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{rad}\text{\hspace{0.17em}}{\mathrm{mm}}^{-1}$. The width of the region corresponds to an area where ${F}^{(I)}(\mathit{\kappa},\delta t)$ is large enough for $\delta t=5\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{ms}$. The resulting 1D phase spectra as functions of ${\kappa}_{x}$ are presented in Fig. 10. The averaged phase of the cross-spectrum ${F}^{(I)}(\mathit{\kappa},\delta t)$ between the jumps fits closely a straight line. Therefore, during the time interval of $\delta t=5\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{ms}$, details of the field $I(\mathit{\rho},t)$ with scales of the order of 1 cm are seen with a large drift at a constant velocity ${u}_{0}$ in the aperture plane in the direction of the $x$ axis at a distance of ${u}_{0}\delta t$ without any significant change or deviation, in accordance to the hypothesis of frozen turbulence. The hypothesis holds better for smaller $\delta t$. Denoting the change of ${\kappa}_{x}$ between the phase jumps as $\mathrm{\Delta}{\kappa}_{x}$ we arrive at the following expression for the drift velocity: ${u}_{0}=2\pi /(\mathrm{\Delta}{\kappa}_{x}\delta t)$. Substituting the values of $\mathrm{\Delta}{\kappa}_{x}=0.6\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{rad}\text{\hspace{0.17em}}{\mathrm{mm}}^{-1}$ and $\delta t=5\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{ms}$, we obtain the estimate of ${u}_{0}=2.1\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}/\mathrm{s}$ for the observations in our study time interval. The vector of ${\mathbf{u}}_{0}$ is directed along the Earth’s surface: ${u}_{0}=\{{u}_{x},0\}$.

For the intensity fluctuations in a plane wave with a constant wind velocity, the value of ${u}_{0}$ coincides with the wind velocity component transverse to the direction of wave propagation [6]. In our observations, the estimated wind velocity along the propagation path, extracted from atmospheric analyses of the European Centre for Medium-Range Weather Forecasts (ECMWF), indicated significant variation both in its absolute value and direction (Fig. 11). The transverse component changes its sign. The scintillations are determined by the structural characteristic of turbulence ${C}_{n}^{2}$ [6]. The dependence of ${C}_{n}^{2}$ from the distance along the propagation path is, however, unknown. The evaluation of the weighting function that determines the relative contribution of turbulence in different sections of the ray path in the observed scintillations, under the condition of strong scintillations and a limited laser beam, is a challenge [11]. This makes it difficult to perform a quantitative comparison between the velocity of the “running shadows” and the transverse wind component. Notwithstanding, our estimate of the order of magnitude of ${u}_{0}$ is reasonably close to the independent information about the wind velocity along the propagation path.

A time span $\delta t$ equal to 5 ms and 10 ms corresponds to a shift ${u}_{0}\delta t$ of the clusters by 10 mm and 20 mm, respectively, which is comparable with the estimates of the characteristic external scales of the clusters inferred from the correlation function $b(\delta \rho )$ (Fig. 7). It follows from the coherence spectra (Fig. 8), for the shift ${u}_{0}\delta t$ of 10 mm, that the correlation between values of the 2D PSD for wave number $\kappa =0.1\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{rad}\text{\hspace{0.17em}}{\mathrm{mm}}^{-1}$ is close to 0.5. The value of $\kappa =0.1$ corresponds to a scale of $2\pi /0.1\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{mm}\approx 63$, which, by the order of magnitude, is close to the estimate of the characteristic size of a cluster. Therefore, clusters, when they drift by a distance of the order of their characteristic scale, are not destroyed. This was clearly noticeable also in visual observations of the “running shadows” on the dome of the OGS observatory.

Thus, Taylor’s hypothesis holds for a spatiotemporal description of the intensity fluctuation field of a laser beam for observations through a turbulent medium at a large distance. A similarity between the evolution of the optical clusters and the evolution of temperature inhomogeneities in the atmosphere [5] is readily observed: Their lifetime roughly coincides with the time of their drift by a distance equal to their characteristic scale.

The estimate of velocity ${u}_{0}$ of the movement of intensity clusters in the telescope aperture was used for an estimate of the anisotropy of observed spectra (4) described by factor $\u3008{(\mathrm{sinc}(\mathit{\kappa}\mathit{u}\tau ))}^{2}\u3009$. In our observations, the exposure time of the GE-680 camera can be assumed to be equal to $2\tau =2\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{ms}$. For ${u}_{0}=2.1\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}/\mathrm{s}$ and $\kappa =0.8\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{rad}/\mathrm{mm}$ (Fig. 4), we infer that $0\le |\tau {\mathit{u}}_{0}\mathit{\kappa}|\le 0.17$ and $1\ge {(\mathrm{sinc}(\mathit{\kappa}{\mathit{u}}_{0}\tau ))}^{2}\ge 0.35$. Therefore, the factor of $\u3008{(\mathrm{sinc}(\mathit{\kappa}\mathit{u}\tau ))}^{2}\u3009$ must result in a greater oblateness of the isolines in the direction of the ${\kappa}_{x}$ axis for larger values of ${\kappa}_{x}$.

To quantitatively estimate the influence of this factor arising at a finite exposure time, assume that the drift velocity of the clusters $\mathit{u}=\tilde{\mathit{u}}+\delta \mathit{u}$, where $\tilde{\mathit{u}}$ is the average drift velocity and $\delta \mathit{u}$ is the velocity fluctuation, and $|\tilde{\mathit{u}}|={u}_{0}$. Let us also assume that the orthogonal components of the 2D vector $\delta \mathit{u}$ are statistically independent and have a normal distribution with dispersion ${\sigma}_{u}^{2}$ and zero mean $\tilde{\delta \mathit{u}}=0$. For the evaluation of $\u3008{(\mathrm{sinc}(\mathit{\kappa}\mathit{u}\tau ))}^{2}\u3009$ we took ${\sigma}_{u}=0.15{u}_{0}$.

Figure 12 presents the 2D PSD of intensity ${F}^{(0)}(\mathit{\kappa};0)$ in the telescope aperture obtained from the following formula:

It is noticeable that isoline shapes for ${F}^{(0)}(\mathit{\kappa};0)>1\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\mathrm{mm}}^{4}/{\mathrm{rad}}^{2}$ (for $\kappa <0.8-0.9\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{rad}/\mathrm{mm}$) in Fig. 12 became closer to circles, compared to the ones in Fig. 4. This indicates that the notion of isotropy of the intensity fluctuation field in the telescope aperture for scales exceeding 6–7 mm agrees with the observations, consistent with previous work in the atmospheric surface layer [20]. The spectral components with larger wave numbers, corresponding to smaller scales, are interfered by the noise. Nevertheless, it is reasonable to assume the statistical isotropy of the scintillation field also for smaller scales or larger wave numbers. A noticeable anisotropy for very small wave numbers in Fig. 12 as well as in Fig. 4 relates to scales of 70–80 mm and more. These scales, however, correspond to the resolution limit of the spectral analysis. In this area, experimental estimates of PSD are therefore unreliable.

## 6. Conclusions

We analyzed the observations of scintillations in a laser beam (532 nm, $\sim 200\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{mW}$ power) traveling along a 144 km path at an altitude of 2.2–2.4 km above sea level, just above the atmospheric boundary layer, between the islands of La Palma and Tenerife. The observations were performed during nighttime on 18 July 2011. The observed intensity field of the laser beam in the telescope aperture with a diameter of 1 m at the end of the propagation path was broken by the turbulence into chaotic clusters.

The estimates of spatial spectra and correlation functions indicated that the observed intensity fields possess, statistically, a locally isotropic structure, which agrees with the idea of a locally isotropic turbulence. Rms intensity deviation divided by the mean intensity in the aperture is about 1.6, which signifies strong scintillations. The estimates of spatial auto-spectra and autocorrelation functions of the intensity field indicated that the characteristic scale of the internal structure of the observed clusters is 6.5–8 mm, while the characteristic size of the clusters is 4–5 cm. The major contribution to the observed scintillations comes from the inhomogeneities of the intensity field with scales from 1–2 cm up to 10–12 cm.

The analysis of the cross-spectra indicated that the hypothesis of frozen turbulence introduced by Taylor can be used for the description of spatiotemporal structure of intensity fluctuations of laser beams traveling through long paths in the atmosphere.

The authors are grateful to S. Schweitzer (Wegener Center, CityUniv. of Graz, Austria) for her complementary support to the scintillation observations as part of the July 2011 Canary Islands experiment. The research was funded by the European Space Agency’s Support to Science Element (STSE) and supported by the Earth Observation Future Missions program, including ESA’s OGS facility. ECMWF is acknowledged for access to its atmospheric analysis data.

## References

**1. **A. M. Obukhov, “Structure of the temperature field in a turbulent flow,” Izv. Akad. Nauk. SSSR Ser. Geogr. Geophys. **13**, 58–69 (1949).

**2. **S. Corrsin, “On the spectrum of isotropic temperature fluctuations in an isotropic turbulence,” J. Appl. Phys. **22**, 469–473 (1951). [CrossRef]

**3. **G. I. Taylor, “The spectrum of turbulence,” Proc. R. Soc. Lond. Ser. A **164**, 476–490 (1938).

**4. **L. R. Tsvang, “Some characteristics of the spectra of temperature pulsations in the boundary layer of the atmosphere,” Izv. Akad. Nauk. SSSR Ser. Geophys. **10**, 1594–1600(1963).

**5. **A. S. Gurvich, “Influence of the temporal evolution of turbulent inhomogeneities on frequency spectra,” Izv. Atmos. Ocean. Phys. **16**, 231–237 (1980).

**6. **V. I. Tatarskii, *The Effects of the Turbulent Atmosphere on Wave Propagation* (Keter, 1971).

**7. **M. A. Vorontsov, G. W. Carhart, V. S. Rao Gudimetla, T. Weyrauch, E. Stevenson, S. L. Lachinova, L. A. Beresnev, J. Liu, K. Rehder, and J. F. Riker, “Characterization of atmospheric turbulence effects over 149 km propagation path using multi-wavelength laser beacons,” in *Proceedings of the 2010 AMOS Conference* (2010), p. E18.

**8. **A. S. Gurvich and V. Kan, “Structure of air density irregularities in the stratosphere from spacecraft observations of stellar scintillation: 1. Three-dimensional spectrum model and recovery of its parameters,” Izv. Atmos. Ocean. Phys. **39**, 300–310 (2003).

**9. **V. F. Sofieva, A. S. Gurvich, F. Dalaudier, and V. Kan, “Reconstruction of internal gravity wave and turbulence parameters in the stratosphere using GOMOS scintillation measurements,” J. Geophys. Res. **112**, D12113 (2007). [CrossRef]

**10. **V. I. Tatarskii, “Shortwave propagation in a medium with random heterogeneities in approximate Markovian random process,” Sov. Phys. JETP **29**, 1133–1138 (1969).

**11. **S. M. Rytov, Yu. A. Kravtsov, and V. I. Tatarskii, *Principles of Statistical Radiophysics, Vol. 3: Elements of Random Fields* (Springer-Verlag, 1987).

**12. **J. S. A. Brooke, P. F. Bernath, G. Kirchengast, C. B. Thomas, J.-G. Wang, K. A. Tereszchuk, G. González Abad, R. J. Hargreaves, C. A. Beale, J. J. Harrison, S. Schweitzer, V. Proschek, P. A. Martin, V. L. Kasyutich, C. Gerbig, O. Kolle, and A. Loescher, “Greenhouse gas measurements over a 144 km open path in the Canary Islands,” Atmos. Meas. Tech.5 (in press).

**13. **G. Kirchengast and S. Schweitzer, “Climate benchmark profiling of greenhouse gases and thermodynamic structure and wind from space,” Geophys. Res. Lett. **38**, L13701(2011). [CrossRef]

**14. **N. Perlot, D. Giggenbach, H. Henniger, J. Horwath, M. Knapek, and K. Zettl, “Measurements of the beam-wave fluctuations over a 142 km atmospheric path,” Proc. SPIE **6304**, 63041O (2006). [CrossRef]

**15. **A. Ishimaru, “The beam wave case and remote sensing,” in *Topics in Applied Physics*J. W. Strohbehn, ed. (Springer-Verlag, 1978), pp. 129–170.

**16. **M. E. Gracheva and A. S. Gurvich, “Strong fluctuations in the intensity of light propagated through the atmosphere close to the Earth,” Radiophys. Quantum Electron. **8**, 511–515 (1965).

**17. **G. R. Ochs and R. S. Lawrence, “Saturation of laser-beam scintillation under conditions of strong atmospheric turbulence,” J. Opt. Soc. Am. **59**, 226–227 (1969). [CrossRef]

**18. **G. R. Ochs, R. R. Bergmtan, and J. R. Snyder, “Laser-beam scintillation over horizontal paths from 5.5 to 145 kilometers,” J. Opt. Soc. Am. **59**, 231–234 (1969). [CrossRef]

**19. **A. M. Yaglom, *Correlation Theory of Stationary and Related Random Functions* (Springer-Verlag, 1987), Vols. 1 and 2.

**20. **A. S. Gurvich, V. V. Pakhomov, and A. M. Cheremukhin, “The isotropy of refractive-index fluctuations of small-scale turbulence in the atmospheric surface layer,” Izv. Atmos. Ocean. Phys. **7**, 49–52 (1971).