The effect of atmospheric refraction on light has been known and studied since antiquity.1 Considerable work has been accomplished on refractive effects related to astronomical observations.23.–4 For terrestrial propagation scenarios, refractive effects can range from relatively subtle effects, such as an apparent shift in object position, to things as spectacular as mirages and the green flash.5,6 Extreme events such as mirages and green flash have been studied extensively and their physics is well known.7 However, very few efforts have been directed at measuring and quantifying the amount of refractive bending one might expect on an ordinary day. With the development of the Global Positioning System (GPS) satellite network, refractivity profiles have been developed using radio occultation of GPS signals.8 In 2013, a patent was awarded to the Boeing Company for a method to measure the refractivity profile of a parcel of atmosphere using airborne imaging sensors.9
An imaging approach to measure the variations in the path-averaged refractive index gradient is described here. The sea level vertical refractive index gradient in the U.S. Standard Atmosphere model is at a wavelength of 500 nm. At any location, the actual gradient varies due to turbulence, local weather conditions, and time of day. It is important to understand the nature and extent of these variations and how they correlate to changing weather conditions. The imaging approach can be applied to any near-horizontal path through the atmosphere and provides a low-cost solution for ground truth measurement of refractive bending. From the refractive bending measurements, estimates of the refractive index gradient can be obtained. Observations can be done over several days, or months at a time, with the camera in time-lapse mode. Besides being low-cost, the ability to measure refractive bending at visible or infrared wavelengths is an added advantage. The imaging-based quantification of refractive bending can help to validate and improve fine-scale numerical weather prediction models by providing quantification of the vertical temperature lapse rate in the lower atmosphere. The vertical temperature lapse rate in the lower atmosphere is strongly influenced by the radiative properties of the surface, thus it can change significantly with cloudiness, terrain, and changes in ground cover (e.g., snow) yet it is only directly measured approximately twice a day with radiosonde balloon launches that are widely spaced geographically. It is currently difficult to capture the radiative nuances of partial cloudiness, terrain, and ground cover with numerical weather models. Thus if refractive bending observations were conducted with low-cost digital photography along with standard surface weather observations, the deduced temperature lapse rate could be used to increase forecast accuracy in the lower atmosphere.
An experiment was conducted at the Air Force Institute of Technology (AFIT) to demonstrate the time-lapse imaging approach. The methodology is described in Sec. 2. Section 3 discusses the results from the experiment and the correlation to meteorological conditions. A summary of the work and future research directions are presented in Sec. 4.
Relationship Between Refractive Index Gradient and Apparent Change in Object Position
Figure 1 shows an example ray trajectory through the atmosphere. Due to refractive bending of the ray, to a remote observer, the object appears to be displaced from its actual position. From Snell’s law, it follows that, for near-horizontal paths through the atmosphere,10
The refractive index is related to the air density, , via the following equation:11 for the dry air index of refraction normalized for an air density of . Edlen’s equation,
By using Eqs. (1), (2), and (4) and , where is the acceleration due to gravity, a relationship between the curvature of an atmospheric ray path and vertical gradients of air density and temperature can be derived:12
Referring back to Fig. 1, if the slope of the ray is close to zero, which is true for near-horizontal paths,
Time-Lapse Imaging Experiment
The experiment was conducted over 10 days in July–August 2014. An imaging system, comprised of a Canon 40D digital camera (pixel pitch ) and a telephoto lens with focal length 300 mm () was used to image the Good Samaritan Hospital building located 8 miles (12.8 km) away. The camera was mounted on a sturdy tripod and was pointed through a window on the ground floor of the laboratory building. The ground floor laboratory setting provided a temperature controlled environment for the camera system and minimized the base disturbances. A cardboard shield was used to shade the tripod from direct sunlight. The camera and the building were at almost the same elevation; hence the path through turbulence was nearly horizontal and for the most part, about 60 m above ground level.10 An image was captured every minute and written to a laptop computer. A cross-correlation algorithm was used to estimate the image shifts. The images were first cropped to roughly isolate the building, and frames that were not clear enough to provide a good correlation estimate were discarded. A Gaussian window was applied to the images to reduce the effects of the frame edges on the correlation result. The cross-correlation was performed between each image and one of the two reference images. A parabolic fit was applied to the correlation peak to provide an estimate of the shift between the images with subpixel resolution.
Stitching Image Shifts from Different Reference Images
Figure 2 shows example images taken during the afternoon and night of July 25, 2014. Due to different illumination conditions, the daytime images correlated very poorly with the night time images, thus making it difficult to estimate the relative shifts between daytime and nighttime images. Hence, two reference images were chosen: a midday image and a twilight image. The midday image correlated well with all the daytime images, but correlated poorly with the evening and nighttime images. The twilight image, on the other hand, showed good correlation with the late afternoon, evening, and nighttime images, but poor correlation with the morning and early afternoon images. In order to obtain the relative image motion, a cross-correlation algorithm was run between the reference images and every other image in the image bank. The relative shifts with respect to the two reference images are shown in Fig. 3. The jumps/discontinuities in the shift profiles suggest poor correlation. It is apparent from Fig. 3 that during certain periods, especially at daybreak and in the afternoons, the two profiles overlap well. While one profile starts getting noisy after the overlap zones, the other profile picks up right at these zones. The two profiles were stitched at the overlap areas, thus replacing the noisy information from one profile with the useful information obtained from the other profile.
Figure 4 shows changes in the apparent hospital position from the evening of July 24, 2014, to the evening of August 3, 2014. The profile was obtained using the stitching technique described in Sec. 2.3. Gaps in the profile indicate cloudy/foggy/rainy periods when visibility was poor. Two components of image motion are evident in the profile. The faster, random motion is due to turbulence while the slow, vertical drift is due to changes in average refractive index gradient along the path. The thermal inversion that built up during the night caused the index gradient to decrease more with height and hence increase in the curvature of the rays coming from the hospital. During the day, as the ground and the surrounding air warmed up, the thermal inversion disappeared and the index gradient relaxed from the nighttime values and resulted in reduced ray curvature. The hospital building thus appeared to move up during the night and move down during the day. The image motion over the course of a day varied from about a meter on some days to about six times as much on certain days. Using Eq. (10), this corresponds to diurnal variations of to in the average refractive index gradient. The gradient changes were more pronounced during the mornings and nights and less pronounced during the afternoons.
One interesting aspect of the stitched plot in Fig. 4 is that the hospital appears to be almost at the same observed vertical position every afternoon, suggesting the refractive index gradient values are nearly the same each afternoon throughout the observation period. This can be expected as during the afternoon, the boundary layer is well mixed and the conditions are near adiabatic ). A curious dip in the profile beyond the usual floor () is seen in the early evening of July 27, 2014. In addition to satellite imagery that shows periods of cloudiness on that day, radar and observation data indicate brief periods of rain and a thunderstorm in the afternoon followed by clearing; this led to a strong moisture gradient near the surface. Evaporating rainfall created a virtual temperature near the ground a few Kelvin higher than the dry air temperature; with clearing conditions the moisture near the surface quickly decreased with height as the drier air moved in—this caused the virtual temperature to converge to the dry air temperature aloft and thus decrease with height at a lapse rate slightly less than . A “superadiabatic” condition was created, or in other words a layer in the lower atmosphere where the density does not decrease with height as much as it typically does under normal adiabatic conditions. As per Eq. (6), the superadiabatic conditions caused a decrease in curvature for ray paths from the hospital, which caused a drop in the hospital image beyond the level expected for adiabatic conditions. The nighttime excursions in Fig. 4 also show good correlation with meteorological observations. The sharp rise in the profile on the night of July 31 is due to cloud free conditions and weak wind resulting in a very strong thermal inversion. On this night, the temperature increased with height above the ground for the lowest 100 to 200 m (a positive lapse rate). On the contrary, presence of clouds and strong winds did not allow a substantial inversion to develop on the nights of July 26, 27, and 28 which explains the short range of motion seen on those nights. The boundary layer collapses to about 100 m, just above the experimental path in the night. Local disturbances create ripples on top of the stratified boundary layer and this rippling behavior is seen around the nighttime peaks in Fig. 4.
To get an estimate of the refractive index gradient values at different times of the day, observation data for surface temperature and humidity at 2 PM on July 25, 2014, was considered. Assuming the afternoon lapse rate is , (for 500 nm), and , the vertical refractive index gradient computed using Eq. (6) is , which is a smaller gradient in the vertical than the standard atmosphere value (). This is expected, as the lapse rate in the standard atmosphere model is , and it implies there is more refractive bending in a standard atmosphere than in a neutrally buoyant atmosphere (). Since the gradients were observed to be nearly the same every afternoon, this value of was taken as a representative value for the afternoon refractive index gradient during the observation period. Thus, the maximum nighttime values of the average refractive index gradient varied from to .
The U.S. Standard Atmosphere model provides a single value for the atmospheric refractive index gradient. The actual refractive index gradient in any location on any given day varies in correlation with solar heating and other weather conditions. A time-lapse imaging experiment was conducted to estimate the temporal variability of the refractive index gradient. The diurnal variation of the average refractive index gradient was found to be more than twice the standard atmosphere value for this gradient on certain days. The gradient changes were more pronounced during the nights and early mornings. The average refractive index gradient was nearly the same each afternoon throughout the observation period when the atmospheric conditions were nearly adiabatic. Comparisons with observational data showed the gradient variations estimated from the time-lapse imagery correlated well with changing weather conditions too.
The imaging approach can be used as a ground truth measurement of refractive bending for comparison to derived estimation methods, such as those obtained from temperature and moisture gradients. It can thus be used for validation and improvement of fine-scale numerical weather models. Future work will include estimating atmospheric turbulence parameters from the individual and differential random motions of different features in the images.
This work was supported by the Air Force Office of Scientific Research (AFOSR) Multidisciplinary Research Program of the University Research Initiative (MURI) Grant. This research was also supported in part by an appointment at the Air Force Institute of Technology administered by the Oak Ridge Institute for Science and Education through an interagency agreement between the U.S. Department of Energy and AFIT. The views expressed in this paper are those of the authors and do not necessarily reflect the official policy or position of the Air Force, the Department of Defense, or the U.S. Government.