The performance of free-space optical applications can be improved using beams of different wavelengths for the auxiliary actions of pointing/tracking or turbulence correction. Chromatic dispersion owing to the atmosphere is an issue for multiwavelength systems, and the dispersion of electromagnetic signals is typically predicted based on refractive conditions from standard atmospheric models. However, for long near-horizontal paths near the Earth’s surface, substantial refractive index gradients that are associated with features such as inverse temperature layers and ducts can be encountered. These features can significantly alter the ray trajectory, the chromatic divergence, and the angle of arrival of directional beams relative to standard atmosphere predictions. A ray tracing approach was implemented to examine the chromatic divergence and angle of arrival of the rays through various practical and extreme atmospheric conditions involving a temperature inversion layer. Over a distance of 150 km along the ground, a brief encounter with the layer can cause pairs of rays with wavelengths 532 and 1550 nm to diverge up to 4.5 times greater than their standard atmosphere predictions. For a single wavelength, a linear increase of angle of arrival with initial launch angle was found for the standard atmosphere, but this trend was significantly altered in the presence of an inversion layer. Extreme refractive conditions with a large inversion layer were simulated to produce optical ducting over long distances. Chromatic separation of rays as large as 280 m was observed when only one of the two wavelengths remained in the duct.

## 1.

## Introduction

The propagation of an optical signal through the clear atmosphere is affected mainly by atmospheric refraction and turbulence. As light propagates through the Earth’s lower atmosphere in a horizontal or near-horizontal direction, refraction is responsible for some of the most interesting atmospheric optical phenomena such as mirages—the Fata Morgana and Novaya Zemlya effects—the green flash, and ducting. These phenomena are usually associated with the formation of significant refractive index gradients close to the ground.

The value of the atmospheric refractive index is close to unity and it varies with temperature, pressure, humidity, altitude, and wavelength of the signal. For the standard atmosphere, the vertical refractive index gradient is $\sim -3.9\times {10}^{-5}\text{\hspace{0.17em}\hspace{0.17em}}{\mathrm{km}}^{-1}$ near the ground at optical frequencies and its magnitude decreases with increasing altitude ($\sim -3.3\times {10}^{-5}\text{\hspace{0.17em}\hspace{0.17em}}{\mathrm{km}}^{-1}$ at 1-km altitude). Overall, this gradient tends to bend the optical path of the light toward the surface of the Earth. Due to the wavelength dependency of the refractive index and its gradient, the amount of bending of the propagating signal is also wavelength dependent. This leads to chromatic dispersion where the ray paths for different wavelengths separate as they propagate. Dispersion in the atmosphere can generally be ignored for short propagation paths, but it becomes significant for applications such as laser communication or astronomical observation near the horizon, where long propagation paths (tens to hundreds of km) may be involved. Chromatic corrections for these applications are commonly predicted based on the standard atmospheric models. However, the real atmospheric condition may be quite different from the standard models, and unusual features, such as temperature inversion and atmospheric ducts, could significantly affect the chromatic divergence and propagation direction for highly directional beams. Understanding the range of behavior for beams under various dispersion conditions is important when considering correction approaches for applications such as laser communication, LIDAR, target designation, and astronomical observation, where accurate estimation of the propagation path is required for optimum performance.

In this study, we apply a ray tracing algorithm to atmospheric and duct profile models available in the literature that are modified to include a temperature inversion layer. Refractivity of the air as a function of wavelength is modeled with a reduced refractivity coefficient, and a temperature inversion layer is simulated based on a duct profile found in the literature. The chromatic divergence between rays of wavelength 532 nm (green) and 1550 nm (IR) is examined as a function of source height and initial ray trajectory angle relative to the layer position. In certain cases, the encounter with the temperature inversion layer can create significant separation between the rays of different color. The effect of the temperature inversion layer on the angle of arrival (AOA) of the rays at the target location was also studied. Deviation from a linear relationship between the AOA and launch angle was observed as the rays encountered the inversion layer.

## 2.

## Theory and Background

## 2.1.

### Chromatic Dispersion in the Atmosphere

The trajectory of an optical beam propagating through the atmosphere is fundamentally dependent on the dispersion relationship (wavelength dependency of the refractive index) of the air. If the refractive index field for a particular volume of air is known, then the propagation path of an optical beam through the volume can be traced by applying Snell’s law of refraction. The atmospheric refractive index profile for dry air can be expressed as^{1}

^{1}

^{,}

^{2}

## Eq. (2)

$${A}_{D}(\lambda )={10}^{-8}[8342.13+\frac{2406030}{130-\frac{1}{{\lambda}^{2}}}+\frac{15997}{38.9-\frac{1}{{\lambda}^{2}}}]\frac{288.2}{1013.25},$$## 2.2.

### Standard Atmospheric Models

For numeric simulation of optical ray trajectory through the inhomogeneous atmosphere for a specific wavelength, the altitude-dependent (vertical) refractive index profile is required. From Eq. (1), the vertical refractive index can be obtained from altitude-dependent pressure and temperature profiles. We note that at optical frequencies, the humidity-dependent contribution to the refractive index is very small and can be ignored.^{3} The temperature and pressure profiles of the atmosphere are commonly modeled based on geographic reference atmospheric models. These models are described by a sequence of spherical layers in the lower troposphere with each layer having a constant vertical temperature gradient.

The altitude-dependent atmospheric properties of six commonly cited reference atmospheric models are listed in a tabular form in an Air Force Geophysics Laboratory technical report by Anderson et al.,^{4} which is publicly available. These layered atmospheric models assume constant temperature gradients for the vertical temperature profile for each layer. In our simulations, we applied the 1976 U.S. Standard Atmosphere and Subarctic Winter Atmosphere reference models. The Subarctic Winter Atmosphere has the largest temperature gradient among the six reference atmospheric models and hence tends to cause the largest bending of the propagating optical signals. For these atmospheric models, the temperature profile in a given layer in the lower troposphere is modeled as

^{5}where ${P}_{0}$ is the sea-level atmospheric pressure in hPa units (1013.25 hPa), ${M}_{a}$ is the mean molar mass of air molecules in gm/mol, $g=9.81\text{\hspace{0.17em}\hspace{0.17em}}{\mathrm{ms}}^{-2}$ is the gravitational acceleration, and $R=8.3144\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{J}{\mathrm{K}}^{-1}\text{\hspace{0.17em}}{\mathrm{mol}}^{-1}$ is the ideal gas constant.

Inserting Eqs. (3) and (4) into Eq. (1), we get the following expression for the atmospheric refractive index profile:

## Eq. (5)

$$n(\lambda ,h)=1+{A}_{D}(\lambda )\frac{{P}_{0}}{{T}_{0}}{\left(\frac{{T}_{0}-\alpha h}{{T}_{0}}\right)}^{\frac{{M}_{a}g}{R\alpha}-1}.$$Considering the reduced refractivity ${A}_{D}(\lambda )$ to be independent of the altitude, an expression for the refractive index gradient can be obtained analytically, and it is given as

## 2.3.

### Refractivity, Effect of Earth Curvature, Modified Refractivity, and Refractive Conditions

The value of the refractive index of the atmosphere is close to unity, and its value varies slowly in the vertical direction for the reference atmospheric models. Therefore, the refractive index at optical frequencies only changes in the fifth to sixth decimal places over the visible spectra. Hence, it is convenient to use the quantity:

which is known as the refractivity. The refractive index gradient is related to the gradient of the refractivity byFor long horizontal or low-angle optical propagation, the refractivity model should include the curvature of the Earth. Optical signals traveling at low grazing angles tend to move away from the Earth’s surface due to the Earth’s curvature. To account for this effect and simplify the analyses, the modified refractivity ($M$) is used, where it is given as

## Eq. (9)

$$M=N+\left(\frac{h}{{R}_{E}}\right)\times {10}^{6}=N+157\times h\text{\hspace{0.17em}\hspace{0.17em}}(M-\text{units}),$$^{6}As we are primarily concerned with horizontal or nearly horizontal propagation of optical beams, we used the modified refractivity and its gradient in our numeric simulation model.

In numerical modeling, the gradient of the modified refractivity is used to simulate actual ray trajectory for horizontal and low-angle optical signal propagation. Based on the modified refractivity gradient values, the refractive conditions can be divided into four categories (Fig. 1).^{7}^{,}^{8} For radar or optical ray tracing, when $\partial M/\partial h>157$ $M$-units/km, rays will propagate away from the Earth’s surface, and this is referred to as the subrefraction condition. When $157\text{\hspace{0.17em}\hspace{0.17em}}M\text{-}\text{units}/\mathrm{km}\ge \partial M/\partial h>118\text{\hspace{0.17em}\hspace{0.17em}}M\text{-}\text{units}/\mathrm{km}$, standard wave propagation occurs, and this is known as the standard condition. With $118\text{\hspace{0.17em}\hspace{0.17em}}M\text{-}\text{units}/\mathrm{km}\ge \partial M/\partial h>0\text{\hspace{0.17em}\hspace{0.17em}}M\text{-}\text{units}/\mathrm{km}$, rays bend downward, and this condition is termed as super-refraction. In ducting/trapping condition where $\partial M/\partial h\le 0$, the rays (particularly for radar) can become trapped between the ground and the ducting layer. As a point of comparison, for the standard condition, the radio refractivity gradient $dN/dh$ is $\sim -40\text{\hspace{0.17em}}\mathrm{N}\text{-units}/\mathrm{km}$ near the Earth’s surface whereas for the visible spectra the value is $\sim -39\text{\hspace{0.17em}}\mathrm{N}\text{-units}/\mathrm{km}$.

For our study, we are particularly interested in the super-refraction and ducting region as the temperature inversion layer creates strong negative refractivity gradients for optical frequencies. The gradients in the ducting layer or temperature inversion layers can give rise to chromatic dispersion and significantly impact the propagation path of optical signals. The inversion layer parameters for this study were chosen based on experimental results from time-lapse camera measurements that were carried out in Las Cruces, New Mexico.^{9} Generally, the U.S. Standard Atmosphere with the modeled temperature inversion layer produces the super-refraction condition for our simulations. To study extreme atmospheric refractive conditions, we also slightly modified and applied the aforementioned temperature inversion profile to the Subarctic Winter Atmosphere. Under certain conditions, ducting of optical signals in the vicinity of the inversion layer was observed in our subarctic winter results.

## 2.4.

### Temperature Inversion, Atmospheric Ducts and Webster Duct Model

In layered tropospheric models, such as 1976 U.S. Standard or Subarctic Winter models, the temperature is usually highest near the ground and generally decreases with altitude. Temperature inversion occurs when this trend is broken and the temperature increases with altitude. A hot layer of air can become trapped between two relatively cold layers to form a temperature inversion layer, and such a structure can extend over hundreds of kilometers along the ground. Conditions such as rapid cooling of the Earth’s surface and the adjacent air early in the night or adiabatic compression of sinking air over the land due to certain regional geographic features can create these sandwiched temperature layers. The denser (cool) air usually has a higher refractive index than the rarer (hot) layer of air in the troposphere. A sufficiently high gradient can trap optical/electromagnetic signals by creating a propagation channel where the rays swing up and down within the atmosphere between the ducting layer and Earth’s surface. This phenomenon is known as ducting or trapping. The presence of temperature inversion layers over long distances can lead to propagation channels or ducts for optical/electromagnetic signals and could increase the range of radio signals to over thousand miles with very low path loss.

Based on their formation process and modified refractivity profiles, the atmospheric ducts are commonly classified into four types: evaporation ducts, surface-based ducts, elevated ducts, and surface ducts.^{8} Generally, the modified refractivity for evaporation ducts is modeled with log-linear function whereas surface-based and elevated ducts are modeled as a trilinear curve and the surface ducts are modeled with a bilinear curve. Evaporation ducts are particularly important in radio refractivity as radio signals interact with the polarity of water vapor molecules. Thus, radio refractivity is a strong function of relative humidity.^{3}^{,}^{10} For optical signals, humidity has a negligible effect, so evaporation ducts are ignored for our study.^{3} Figure 2 shows representative refractivity profiles of three duct types, excluding the evaporation duct.

There are several detailed analytical models available in the literature describing the form of the inversion profile for the duct types mentioned above. For our study, we chose the Webster duct profile^{11} and applied this profile in the surface-based, elevated, and surface duct types to examine their effect on chromatic dispersion. The Webster duct model, introduced by Webster,^{11} was later modified by Grabner and Kvicera,^{12} where the refractivity profile is defined by the following equation:

## Eq. (10)

$$N(h)={N}_{0}+{G}_{N}h+\frac{\mathrm{\Delta}N}{2}\text{\hspace{0.17em}}\mathrm{tanh}\left[\frac{2.96(h-{h}_{0})}{\mathrm{\Delta}h}\right],$$## Eq. (11)

$$N(\lambda ,h)=[{A}_{D}(\lambda )\frac{{P}_{0}}{{T}_{0}}{\left(\frac{{T}_{0}-\alpha h}{{T}_{0}}\right)}^{\frac{{M}_{a}g}{R\alpha}-1}]\times {10}^{6}+\frac{\mathrm{\Delta}N}{2}\text{\hspace{0.17em}}\mathrm{tanh}\left[\frac{2.96(h-{h}_{0})}{\mathrm{\Delta}h}\right].$$Taking the Earth’s curvature into consideration, the wavelength- and altitude-dependent modified refractivity then becomes

## Eq. (12)

$$M(\lambda ,h)=[{A}_{D}(\lambda )\frac{{P}_{0}}{{T}_{0}}{\left(\frac{{T}_{0}-\alpha h}{{T}_{0}}\right)}^{\frac{{M}_{a}g}{R\alpha}-1}]\times {10}^{6}+\frac{\mathrm{\Delta}N}{2}\text{\hspace{0.17em}}\mathrm{tanh}\left[\frac{2.96(h-{h}_{0})}{\mathrm{\Delta}h}\right]+157\times h.$$For the simulations results presented here, we used a parameterization of the Webster duct model published by Wijerathna et al.^{9} These authors analyzed time-lapse imagery, and the results suggested that a temperature inversion-type of layer was present in the early morning at their location in New Mexico. Parameters for the Webster duct model were obtained by a fitting approach associated with the image shift results. For five consecutive mornings in late November, 2014, these authors found the mean of the calculated duct width to be $\mathrm{\Delta}h=48.5\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{m}$ and the duct height was found to be ${h}_{0}=10.8\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{m}$. The parameters were calculated for the time of the day where the maximum image stretching (distortion) was observed and the authors assumed $\mathrm{\Delta}N$ to be $-2$ N-units. This profile represents a surface-based temperature inversion profile. The three temperature inversion profiles (surface, surface-based, and elevated) were modeled in our simulations by changing the position of the duct height ${h}_{0}$.

## 3.

## Simulation and Results

## 3.1.

### Ray Tracing Method and Algorithm

When refraction effects are much larger than diffractive effects, light propagation through an inhomogeneous media can be simulated using ray tracing methods. This condition generally holds for the propagation path and refractive features that we are studying. The ray trajectory is usually described by the solution of the eikonal equation.^{13}^{–}^{17} A commonly used solution involves the application of a linear transfer equation coupled with a bending equation. The bending equation determines the change in ray trajectory and is a function of change in refractive index between two layers. Starting with Snell’s law of refraction,

## Eq. (13)

$${n}_{1}\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}{\theta}_{1}={n}_{2}\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}{\theta}_{2}$$The small-angle approximation holds for our numeric simulations as the rays travel in a horizontal/near-horizontal path, with a total angular change of $\sim 20\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{mrad}$ over a 150-km path along the ground. The bending equation between layers is obtained from Eq. (14) as

With index $j$ for discrete layers, the bending equation for the numeric simulation is

The rays are assumed to follow a linear path between planes with the linear transfer equation given as

where two-dimensional (2-D) ray propagation is considered and $\mathrm{\Delta}{x}_{j-1}$ is the horizontal step size for the numeric simulation. To obtain accurate ray trajectories with this approach in a medium with significant index gradients, relatively small step sizes are required in both vertical and horizontal planes. Ray height $h$ is considered in the vertical direction. Hence, although simple, this numerical model for ray tracing is not time efficient when long paths (tens to hundreds of kilometers) are considered.Our problem of interest involves refractive index gradients defined in the vertical direction and ray trajectories at low angles relative to the horizontal. For this situation, a second-order ray tracing algorithm^{18} was developed that relaxes the step size constraint and allows for faster but accurate ray trajectory calculation. We implemented this second-order ray tracing method in our simulation, where a quadratic correction term is introduced to the transfer equation. This term comes from the solution of a 2-D form of the eikonal equation describing the ray trajectory in an inhomogeneous media. The eikonal equation is given as

For numeric simulations, the transfer equation [Eq. (19)] and bending equation [Eq. (20)] can be rewritten as follows:

and respectively. The additional term in the transfer equation models a curved ray trajectory between the steps that gives an accurate ray height approximation with fewer iterations, thus reducing the computational time. This is particularly important in our study as the rays propagate through the inhomogeneous atmosphere over a very large distance along the ground.## 3.2.

### Chromatic Dispersion Simulation Results

The chromatic dispersion between rays with wavelength 532 and 1550 nm was simulated using the second-order ray tracing algorithm for several practical atmospheric conditions for a distance of 150 km along the ground. Figure 3 provides an example of how the presence of an elevated inversion layer, described by the Webster duct model, impacts the ray trajectory and chromatic divergence.

Figure 3(a) shows ray trajectories through the U.S. Standard Atmosphere gradient profile. Three pairs of 532- and 1550-nm rays are traced from left to right, and the apparent upward trajectory of the rays beyond the 100-km distance is an effect related to the curvature of the Earth’s surface. Each pair is launched at a slightly different downward initial angle and as the rays traverse the atmosphere, the 532-nm wavelength rays bend more toward the ground than the 1550-nm rays. Because the refractivity gradient is altitude dependent (indicated by the color bar), the overall amount of ray bending is primarily dependent on the source altitude. Figure 3(b) shows the separation between rays within a pair (${h}_{\mathrm{diff}}$) as a function of distance along the Earth’s surface. Regardless of the slight launch angle difference, the rays for this case diverge about 5 m within each pair at the 150-km distance along the Earth’s surface in the standard atmosphere. In general, we find for the standard atmosphere case that ${h}_{\mathrm{diff}}$ is also primarily a function of source altitude, where the separation becomes smaller as the source is moved higher in altitude.

However, if the source is located above the inversion layer [Fig. 3(c)], then the shorter wavelength ray encounters the layer feature earlier (in terms of distance along the Earth’s surface) than the longer wavelength. This can separate the rays with different colors further or even bring them closer to their standard atmospheric prediction, depending on several factors. Figure 3(d) shows chromatic separation between ray pairs launched at the same angles as in Fig. 3(a) but in the presence of the inversion layer. The ray separation depends on the relative source position with respect to the inversion layer, launch angle, and wavelength of the source. For the case shown in Fig. 3, the separation can be as large as 23 m, which is significantly larger than the $\sim 5\text{-}\mathrm{m}$ prediction for the standard atmosphere.

Figure 4 presents ray separation results for packets of 101 pairs of rays, where each pair is aimed at a slightly different angle toward the inversion layer. Results are shown for source locations below, inside, and above the inversion layer. For a specific source height (${h}_{\mathrm{s}}$), ray packets were launched over an angular range of a few milliradians, such that the rays would encounter the inversion layer. Figure 4 shows ray separation as a function of distance along the Earth’s surface for four source altitudes. For the case in Fig. 4(a) with the source below the layer, the result for the ray packet is essentially the same as Fig. 3(b) without the layer. The other cases illustrate that the largest separation effects occur when the source is above the layer and the rays are launched at larger initial angles.

For all of the results in Sec. 3.2, the elevated temperature inversion layer was modeled using the following Webster duct parameter values: $\mathrm{\Delta}N=-2$ $N$ units, $\mathrm{\Delta}h=50\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{m}$, ${h}_{0}=50\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{m}$. The value of the duct width $\mathrm{\Delta}h$ was 50 m (average value is 48.5 m in the work by Wijerathna et al.^{9}) and the duct height value ${h}_{0}=50\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{m}$ (average value is 10.8 m in the work by Wijerathna et al.^{9}) was used to model an elevated inversion layer to ensure maximum possible interaction of the propagating rays with the inversion layer. This profile coupled with U.S. Standard Atmosphere does not produce a strong enough gradient for ray trapping or ducting but rather corresponds to a super-refraction refractive condition. We ran simulations with different ${h}_{0}$ values to model surface and surface-based duct type inversion layers as well. Our simulations show that the ray separation ${h}_{\mathrm{diff}}$ can be considered independent of the duct type but strongly dependent on the relative distance between the source height (${h}_{\mathrm{s}}$), duct height (${h}_{0}$), and launch angle.

For the source located above the inversion layer, we ran simulations with the 532-nm wavelength rays to study the ray AOA at the 150-km distance as a function of launch angle. For the standard atmosphere, the AOA increases essentially linearly with increasing launch angle for a particular source height. However, in the presence of an inversion layer, this general trend is broken and the AOA can decrease with increasing launch angle. Figure 5(a) illustrates the dependence of the AOA on launch angle in the standard atmosphere, and Fig. 5(b) shows the effect of the elevated inversion layer (as defined above) on the AOA with the source located at ${h}_{\mathrm{s}}=200\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{m}$ and rays launched at angles between $-6.5$ and $-5.5\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{mrad}$.

## 3.3.

### Ducting and Dispersion Simulation Results Under Extreme Conditions

Observation of ducting of radio signals goes back as early as World War II, and such anomalous refraction is linked with temperature inversion and the evaporation duct.^{19} Hence, it is a curious question to ask whether ducting of optical signals is a possibility. As humidity has a negligible effect on optical signal propagation, we found that a large temperature gradient along with a large temperature inversion is required for optical ducting. The subarctic winter atmospheric profile exhibits a large temperature gradient with height, and literature suggests that large temperature inversion profiles can exist in the subarctic regions.^{20}^{,}^{21} Thus, the subarctic region could produce atmospheric conditions that support optical ducting. To examine this possibility, we applied our simulation approach to the Subarctic Winter Atmosphere model with an embedded temperature inversion. Currently, we do not have data for the parametrization of duct profiles for the subarctic atmosphere and we did not find suitable analytic subarctic inversion models for use in our approach. Hence, we applied a temperature inversion layer identical to the one described in Sec. 3.2 to the Subarctic Winter Atmosphere with a change in the parameter $\mathrm{\Delta}N=-3\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{N}$ units to simulate a stronger inversion. This inversion layer along with the subarctic winter refractive index profile created conditions that cause actual ducting of the optical signal. Figure 6 shows the ducting effect for 101 green (532 nm) rays launched near the center (${h}_{s}=40\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{m}$) of the ducting layer (duct height = 50 m) with initial launch angles ranging from 0 to 0.4 mrad. The ducting of rays limits the chromatic divergence as long as the rays of different wavelength are all trapped. However, our simulation studies show that for specific source altitudes there can exist an angular launch range over which rays with one color can become trapped by the ducting layer and the rays of the other color may not meet the ducting condition. The chromatic separation between 101 pair of 532- and 1550-nm rays is shown in Fig. 7 for the ducting case described above, with different launch angle ranges and the source located at 40 m above the ground. Chromatic separation is small after the 150-km propagation distance when both colors are ducted [Fig. 7(a)], but for larger launch angles, the 1550-nm rays leave the duct and chromatic separations up to 280 m are demonstrated [Fig. 7(b)]. Although the temperature inversion layer for this example is *ad hoc*, it illustrates the parameters required to produce optical ducting and extreme chromatic separation.

## 4.

## Conclusions

A ray tracing simulation was developed to study the wavelength dependence (dispersion) of optical propagation through atmospheric refractive index features. The dispersion is incorporated through a reduced refractivity coefficient in a modified refractivity expression, and a second-order ray tracing algorithm is applied for improved efficiency and accuracy. Trajectory altitude differences for 532- and 1550-nm wavelength rays were compared using the simulation. We find that for standard atmosphere refraction conditions, the separation between the rays depends largely on the altitude of the source and is nearly independent of the launch angle (for near-horizontal rays with initial angles in milliradians). The chromatic separation becomes smaller as the source is moved higher in altitude.

The presence of a temperature inversion layer introduces significantly higher gradients than the standard atmosphere and can cause much more ray separation (factor of 4.5) than the standard atmosphere depending on the relative source position with respect to the inversion layer, launch angle, and wavelengths of the source. For the inversion layer we studied, we find that the largest ray separation effects occur when the source is above the layer and the rays are launched at relatively larger initial angles. It is important to note that the divergence of chromatic dispersion from the standard atmospheric model in the presence of an inversion layer becomes much more significant when long propagation distances ($\ge 100\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{km}$) are involved. The AOA at the detector location increases linearly with increasing launch angle at the source location for the standard atmosphere. However, in the presence of an inversion layer, this trend is broken. It seems possible that this response can be used as an indicator of the presence of such a layer.

We found that the application of the Subarctic Winter Atmosphere model with the addition of a large temperature inversion structure allows for optical ducting and extreme chromatic dispersion. For a source height near the layer height and over a small launch angle window, ducting of the optical signal of one color (532 nm) can occur while the rays of the other color (1550 nm) do not duct. This leads to significant chromatic dispersion, and the simulation produced ray altitude differences as large as 280 m for the 150-km propagation distance.

## Acknowledgments

Earlier results of this work appeared in the conference proceedings paper: M. Abdullah-Al-Mamun and D. Voelz, “Chromatic refraction in the lower atmosphere associated with refractive index gradient features,” Proc. SPIE 11133, Laser Communication and Propagation through the Atmosphere and Oceans VIII 2019, 11133-12 (13 August 2019). This work was supported by the Directed Energy Joint Transition Office (DE JTO); Award No. N00014-17-1-2535.

## References

**,” Appl. Opt., 42 (3), 354 –366 (2003). https://doi.org/10.1364/AO.42.000354 Google Scholar**

*Ray tracing and refraction in the modified US 1976 atmosphere***,” Metrologia, 2 (2), 71 –80 (1966). https://doi.org/10.1088/0026-1394/2/2/002 MTRGAU 0026-1394 Google Scholar**

*The refractive index of air***,” Johns Hopkins APL Tech. Digest, 4 (1), 12 –26 (1983). JHADDQ 0270-5214 Google Scholar**

*Anomalous microwave propagation through atmospheric ducts***,” (1986). Google Scholar**

*AFGL atmospheric constituent profiles (0–120 km)***,” Am. J. Phys., 65 (5), 404 –412 (1997). https://doi.org/10.1119/1.18555 AJPIAS 0002-9505 Google Scholar**

*On the barometric formula***,” Bull. Géodésique, 54 (3), 395 –405 (1980). https://doi.org/10.1007/BF02521480 Google Scholar**

*Geodetic reference system 1980***,” J. Atmos. Oceanic Technol., 31 2650 –2670 (2014). https://doi.org/10.1175/JTECH-D-13-00152.1 Google Scholar**

*Radar beam tracing methods based on atmospheric refractive index***,” IEEE Trans. Antennas Propag., 63 (6), 2728 –2738 (2015). https://doi.org/10.1109/TAP.2015.2418788 Google Scholar**

*Channel model for the surface ducts: large-scale path-loss, delay spread, and AOA***,” in Proc. OSA Imaging and Appl. Opt., (2017). https://doi.org/10.1364/PCAOP.2017.PTu2D.2 Google Scholar**

*Time lapse image stretch measurements to obtain atmospheric refractivity duct model parameters***,” J. Res. Nat. Bur. Stand., 50 (1), 39 –41 (1953). https://doi.org/10.6028/jres.050.006 Google Scholar**

*The constants in the equation for atmospheric refractive index at radio frequencies***,” IEEE Trans. Antennas Propag., 30 (4), 796 –800 (1982). https://doi.org/10.1109/TAP.1982.1142840 Google Scholar**

*Raypath parameters in tropospheric multipath propagation***,” Fundamentals of Photonics, John Wiley & Sons, Inc., Hoboken, New Jersey (2007). Google Scholar**

*Ray optics***,” Fundamentals of Photonics, John Wiley & Sons, Inc., Hoboken, New Jersey (2007). Google Scholar**

*Wave optics***,” Fundamentals of Photonics, John Wiley & Sons, Inc., Hoboken, New Jersey (2007). Google Scholar**

*Fiber optics***,” Gradient Index Optics, Academic Press, New York (1978). Google Scholar**

*Axial gradients***,” Gradient Index Optics, Academic Press, New York (1978). Google Scholar**

*Radial gradients***,” New Mexico State University, (2016). Google Scholar**

*Refraction in the lower troposphere: higher order image distortion effects due to refractive profile curvature***,” Proc. IEEE, 73 (2), 265 –283 (1985). https://doi.org/10.1109/PROC.1985.13138 IEEPAD 0018-9219 Google Scholar**

*Tropospheric radio propagation assessment***,” J. Phys. Oceanog., 35 2444 –2456 (2005). https://doi.org/10.1175/JPO2829.1 Google Scholar**

*Temperature inversions in the subarctic North Pacific*## Biography

**Mohammad Abdullah-Al-Mamun** is a research assistant in the Klipsch School of Electrical and Computer Engineering, New Mexico State University (NMSU). He received his BS and MS degrees in physics from the University of Dhaka, Bangladesh. He also received an MS degree in physics from NMSU in 2018. His research interests include laser propagation through the atmosphere, atmospheric turbulence, adaptive optics, and optical system design. He is a student member of SPIE.

**David Voelz** is a professor in the Klipsch School of Electrical and Computer Engineering, NMSU. His research interests include optical propagation through turbulence, spectral and polarization sensing, imaging theory, astronomical instrumentation development, and laser communications. He received his PhD in electrical engineering from the University of Illinois in 1987. He is a fellow of SPIE and OSA.