# Deriving Equations for Atmospheric Refraction

On this page I explain how terrestrial atmospheric refraction is calculated from basic physical laws and how the equations and constants used in refraction equations are derived. The same principle can be applied to derive astronomical refraction, although this is not the topic of this page.

Atmospheric Refraction is the bending of light due to a gradient of the refractive index caused by a density gradient of the air. This light bending causes observed objects in the distance to be displaced (mainly up) or distorted and inverted (mirage). In a layer above the surface refraction can vary considerably. Above this layer refraction varies only little and we can derive a mean standard value, see Refraction Coefficient as a Function of Altitude.

In survey it is important to know how much light gets bent to be able to correct the observed height of a distant object for Refraction effects.

## Refraction Coefficient Globe

The density gradient along a light ray is difficult to measure in practice. Consequently, the path of a light ray can not be known exactly. But as long as the air is not too turbulent, we can use an average gradient. In this case a light ray can be approximated by an arc with the radius r.

Refraction can be expressed by the Refraction Coefficient k, which is the ratio of the radius of the earth to the radius of a bent light ray:

(1)

per Definition

source

Terrestrial Refraction; Wikipedia

where'
 $k$ ' =' 'Refraction Coefficient, see (7) $c$ ' =' '1/r = curvature of light ray, see (5) $r$ ' =' 'radius of curvature of the light ray $R$ ' =' '6371 km = radius of the earth

The Refraction Coefficient is defined so, that k = 0 is no refraction, k = 1 is light bending down with radius R and k = −1 is light bending up with radius R.

Note: The phenomenon of refraction is not dependent on the shape of the earth. We could define a refraction coefficient that is not related to the radius of the earth. In fact we could simply use the curvature c = 1/r of the light ray. We can always get this curvature from k by dividing all equations for k through R: c = k/R.

## Refraction Coefficient Flat Earth

We can define a similar refraction coefficient as for the Globe for a Flat Earth. Like for the Globe the FE Refraction Coefficient is defined so, that k = 0 is no refraction, k = 1 is light bending down with radius R and k = −1 is light bending up with radius R.

(2)
where'
 $k_\mathrm{FE}$ ' =' 'Refraction Coefficient for Flat Earth $k$ ' =' 'Refraction Coefficient for Globe, see (7) $c_\mathrm{FE}$ ' =' '1/rFE = curvature of light rays on flat earth, see (5) $r_\mathrm{FE}$ ' =' 'radius of curvature of light rays on flat earth $R$ ' =' '6371 km = Radius of the Earth
k kFE Description
0 -1 Earth appears as a Globe with radius R
0.17 -0.83 Earth appears as a Globe with an apparent radius of 6/5 R, see Refraction Factor
1 0 Earth appears flat

As the curvature of light can be calculated from the Refraction Coefficient: $c = k_\mathrm{FE} / R$, we get that for a flat earth to appear as a Globe, kFE has to be negative which implies a negative curvature of light, light must bend upwards.

## Refraction Factor, Apparent Radius of Earth

Many calculations for observations that involve the radius of the earth R, like for hidden height, the distance to the horizon ect., assume a straight light ray. But there is a trick to keep the same equations for bent light, by increasing the radius of the earth in the equations by the Refraction Factor a and making the light ray straight. The geometry of the two situations is the same. So we can simply replace R in this equations by R':

(3)
where'
 $R^{\,\prime}$ ' =' 'Refracted Apparent Radius of the Earth $R$ ' =' '6371 km = geometrical radius of the earth $a$ ' =' 'Refraction Factor

Note: because our eyes assume all light rays are straight and interprets the scene accordingly, the earth appears to us scaled by the Refraction Factor.

The connection between the Refraction Coefficient k and the Refraction Factor a is:

(4)
where'
 $a$ ' =' 'Refraction Factor (in some literature called K-value) $k$ ' =' 'Refraction Coefficient, see (1)
Refraction Coeff k Refraction Factor a Description
0 1 no Refraction
0.143 7/6 Often used Standard Refraction in Geodesy
0.17 6/5 Standard Refraction at sea level
0.25 4/3 Standard Refraction for Microwaves

If refraction is such that light follows the curvature of the earth (k = 1), then the Refraction Factor is infinite. That means the apparent earth radius is infinity, i.e. the earth appears flat.

## Calculating Curvature of Light

If we have measured atmospheric parameters like the pressure, temperature and temperature gradient, we can use the following formula to calculate the the curvature of light:

(5)
see
where'
 $c$ ' =' 'Curvature of light ray $r$ ' =' '1/c radius of curvature of the light ray $P$ ' =' 'Air pressure at the observer in mbar or hPa or 1/100 Pa, Standard = 1013.25 mbar $T$ ' =' 'Temperature at the observer in Kelvin, Standard = 288.15 K = 15°C $\mathrm{d} T/\mathrm{d} h$ ' =' 'Temperature Gradient in K/m or °C/m, Standard = −0.0065°C/m

Note: this formula is independent on the shape of the earth.

If we know the apparent size of the earth R' from measurements and assuming a Globe with radius R, or if we know the Refraction Coefficient k, we can calculate the curvature of the light rays as follows:

(6)
where'
 $c$ ' =' 'Curvature of light ray on the Globe $c_\mathrm{FE}$ ' =' 'Curvature of light ray on the Flat Earth $r$ ' =' '1/c = radius of curvature of the light ray $R^{\,\prime}$ ' =' 'Apparent Radius of the Earth $R$ ' =' '6371 km = Real Radius of the Earth $k$ ' =' 'Refraction Coefficient on the Globe, see (7) $k_\mathrm{FE}$ ' =' 'Refraction Coefficient on the Flat Earth

From the curvature of the light c or cFE we can calculate the prevailing temperature gradient using (8) for either model respectively.

## Calculating Refraction Coefficient

If we know the curvature of the light rays (5) or if we know the pressure, temperature and temperature gradient at the observer and assume the same conditions along the line of sight, we can use the following formula to calculate the Refraction Coefficient [1]:

(7)
see
where'
 $k$ ' =' 'Refraction Coefficient $P$ ' =' 'Air pressure at the observer in mbar or hPa or 1/100 Pa, Standard = 1013.25 mbar $T$ ' =' 'Temperature at the observer in Kelvin, Standard = 288.15 K = 15°C $\mathrm{d} T/\mathrm{d} h$ ' =' 'Temperature Gradient in K/m or °C/m, Standard = −0.0065°C/m

If we know the curvature of light, e.g. by measuring the apparent size of the earth from observerations and using formula (6), or the refraction coefficient, we can solve the equation (5) for the temperature gradient that is causing the measured refraction:

(8)
where'
 $\mathrm{d} T / \mathrm{d} h$ ' =' 'Temperature Gradient in °C/m or K/m $c$ ' =' 'Curvature of light rays, see (6) $r$ ' =' '1/c = Radius of curvature of the light ray $T$ ' =' 'Absolute Temperature in Kelvin (standard 288.15 K) $P$ ' =' 'Pressure in mbar (standard 1013.25 mbar)

This formula is not dependent on the shape of the earth, so it works for the Globe and the Flat Earth model. But the calculated curvature of light from observations dependends on the assumed model.

So we can use the different calculations of the curvature of light (6) for each model and then using the formula above calculate the required temperature gradient for both models. It turns out that, to get actual observations above the ground layer, the required temperature gradient on a Flat Earth is physically impossible, about −170°C/km at standard conditions, all the way up in every altitude.

As the curvature of light rays is proportional to the negative Refractivity Gradient (20) wen can use (5) to convert between temperature gradients and refractivity gradients with the following equations:

(9)
(10)
where'
 $\mathrm{d} T/\mathrm{d} h$ ' =' 'Temperature-Gradient in K/m or °C/m, Standard = −0.0065°C/m $\mathrm{d} N/\mathrm{d} h$ ' =' 'refractivity gradient in 1/m, standard is −22.4/km for k = 0.143 and a = 7/6 $P$ ' =' 'Air pressure in mbar or hPa or 1/100 Pa, Standard = 1013.25 mbar $T$ ' =' 'Temperature in Kelvin, Standard = 288.15 K = 15°C $N$ ' =' '(n − 1) · 106 = refractivity $n$ ' =' 'refraction index $h$ ' =' 'altitude

## How does Refraction work?

Before we can derive equations for atmospheric refraction we first have to investigate how refraction works. We know that the speed of light in a medium is slower than in a vacuum. This is expressed by the refractive index:

(11)

Definition Index of Refraction

where'
 $n$ ' =' 'Refractive Index or Index of Refraction, a dimensionless number ≥ 1; for vacuum n = 1 $c$ ' =' '299,792,458 m/s = speed of light in vaccum [2] $v$ ' =' 'speed of light in a medium

The refractive index depends on the density and composition of the medium and the wavelength of the light. In a homogeneous solid like glas the refractive index is constant. As soon as light, or any electromagnetic wave, encounters a boundary with a sudden change in the refractive index, the light gets bent by an angle that can be calculated by Snell's law if we know the refractive index of both sides of the boundary. Light gets always bent towards the denser medium, that is the medium with a greater refractive index.

Explanation of how Refraction works; Wikipedia
https://en.wikipedia.org/wiki/Refraction

Rather dark Anatidae and the basics of refraction
Excellent derivation of Snell's law and the refractive index by AB science

A typical value for the refractive index of air at sea level under standard atmospheric conditions is no = 1.000278. Since it is so close to 1, the difference between the refractive index of air and that in a vacuum (n = 1) is of most interest. Since n−1 is very small, it is useful to define Refractivity, which is equal to:

(12)

Definition Refractivity

where'
 $N$ ' =' 'Refractivity $n$ ' =' 'Refractive Index or Index of Refraction

So Refractivity of air at sea level is typically N = 278 N-units.

## Refraction in the Atmosphere

Fig. 1: Ray bundle traversing the atmosphere

Refraction in the atmosphere is more complicated than in a solid with constant refractive index. The air density decreases continuously with increasing altitude. This is called a density gradient. Consequently the refractive index also decreases continuously with increasing altitude until it reaches 1 in the vaccum of space. How is light affected by a refractive index gradient? We can't apply Snell's law here because we have no boundaries.

But based on Fermat's principle we can use calculus to derive an equation that relates curvature of a light ray to any refractive index gradient.

Consider a bundle of light rays propagating through the air with a refractive index gradient as shown in Fig. 1. The radius of curvature of the light ray is denoted as r. All rays in the bundle must traverse the air in the same time dt, if they are part of the same bundle. Lets consider two rays in that bundle.

Path 2: The wave on this path has a radius of r and a velocity v and traverses the path in a time dt. The length of the path that is traversed in time dt is r·dθ. Because length is velocity times time, we can write this as:

 (13)

Path 1: The refractive index n is different since the path is higher in the atmosphere. We denote the refractive index here as n + dn where dn is the difference in the indices. The radius of this path is r + dr. The associated velocity is v + dv and the time interval dt is the same for both rays. Similar to (13) we can write for this ray:

 (14)

We can get rid of dθ and dt by dividing equations (13) and (14):

 (15)

Subtracting 1 from both sides we obtain:

 (16)

The definition of the refractive index is n = c/v, where c is the speed of light in vacuum and v the speed of light in the medium. Differentiating with respect to v gives:

 (17)

We have now a relation between the refractive index and the speed of light in the medium. But we need a relation between the refractive index and the radius of curvature of the light ray. We can get this by replacing the green term in (17) with the blue term of (16):

 (18)

The inverse of the radius of curvature r is called curvature. Hence, the relation (18) describes the curvature of the ray as a function of the refractive index n and the refractive index gradient dn/dh.

We can make some simplifications: if we assume that the ray propagates close to parallel to the ground, so that the angle β in the figure is small, then cos(β) ≈ 1. The refractive index n of air is very close to 1, so we can neglegt the term 1/n.

So the curvature of the light ray can be expressed as:

(19)
where'
 $c$ ' =' 'curvature of the bent light ray (0 = straight) $r$ ' =' 'radius of curvature of the bent light ray $\beta$ ' =' 'inclination angle with respect to horizontal $n$ ' =' 'refractive index $\mathrm{d} n/\mathrm{d} h$ ' =' 'refractive index gradient

The curvature of light in air is proportional to the negative refractive index gradient.

Note: As the refractive index decreases with increasing altitude, the refractive index gradient is negative, but due to the minus sign the curvature is positive. From our derivation follows, that a positive curvature means that light gets bent down toward the surface. This causes objects in the distance to appear higher than they are.

This is true even if the light ray starts initialliy horizontal. This is due to the gradient of the refractive index, not the curvature of the atmosphere. So even on a flat earth with a vertical gradient, horizontal light gets bent down towards the denser part of the atmosphere.

Note: If the refractive index gradient is 0, the curvature 1/r of the light ray is zero, which means the light ray is straight.

Now we can calculate the curvature 1/r of a light ray at any point from the refractive index gradient dn/dh at that point. And as the Refraction Coefficient k is a measure of this curvature, using the relation (1), the Refraction Coefficient can be calculated from this gradient:

(20)
and
where'
 $c$ ' =' '$1 \over r$ = curvature of light rays $r$ ' =' 'radius of curvature of light rays $k$ ' =' 'Refraction Coefficient $R$ ' =' '6371 km = Radius of the earth $\mathrm{d} N/\mathrm{d} h$ ' =' 'refractivity gradient in 1/km, standard is −22.4/km for k = 0.143 $\mathrm{d} n/\mathrm{d} h$ ' =' 'refractive index gradient in 1/km $N$ ' =' '(n − 1) · 106 = refractivity $n$ ' =' 'refractive index $h$ ' =' 'altitude

I found this derivation on the following PDF:

Note: The end result in the linked document is expressed as 1/K, where K is called the K-value. I call this value the Refraction Factor a.

## Calculating Refractivity of Air

Before we can derive the equation for the Rrefraction Coefficient k as a function of atmospheric pressure, temperature and temperature gradient, we have to calculate the Refractivity of the air N dependent on some atmospheric parameters.

The refractivity of the atmosphere at a certain point is a function of many things, including the wavelength of light, the air pressure P, the absolute temperature T, and the humidity e at that point. A commonly used expression is [3] [4]:

(21)
where'
 $N$ ' =' 'Refractivity $P$ ' =' 'Air pressure in mbar or hPa or 1/100 Pa $T$ ' =' 'Absolute Temperature in Kelvin $e$ ' =' 'Humidity expressed as the partial pressure of water vapor, in mbar $K_1$ ' =' '79.0 K/mbar for light (λ = 550 nm) $K_2$ ' =' '4810 K

For an accuracy of less than 0.5% some simplifying assumptions may be made if we limit the range of certain variables [4]:

• Temperature: −50°C ≤ T+40°C
• Pressure: 200 mbar ≤ P1100 mbar
• Humidity and CO2 content can be neglegted

So equation (21) can be simplified to:

(22)
where'
 $N$ ' =' 'Refractivity $P$ ' =' 'Air pressure in mbar or hPa or 1/100 Pa $T$ ' =' 'Absolute Temperature in Kelvin $K_1$ ' =' '79.0 K/mbar for light (λ = 550 nm), see (24)

The constant K1 is calculated as follows: Since according to the gas law the density of air is proportional to P/T and Refractivity N according to (22) is also proportional to P/T, Refractivity is proportional to the density of the air. So we can write:

(23)
gas law
where'
 $N(h)$ ' =' 'refractivity at altitude h $N_o$ ' =' 'refractivity at sea level at International Standard Atmosphere, see below $N_o/\rho_o$ ' =' 'proportionality factor $\rho_o$ ' =' '1.225 kg/m3 = density of air at sea level and standard atmosphere $\rho(h)$ ' =' 'density of air at altitude h $P(h)$ ' =' 'air pressure at altitude h in Pa $P'(h)$ ' =' 'air pressure at altitude h in mbar $T(h)$ ' =' 'absolute temperature at altitude h $R_\mathrm{S}$ ' =' '287.058 J/(kg·K) = Specific Gas Constant of dry air

Note: The factor 100 is present so we can use the pressure P in mbar instead of Pa (Pascal).

The blue term is our K1. Now we only have to find the refractivity No for light of wavelength λ = 550 nm, standard pressure at sea level Po = 1013.25 mbar and standard absolute temperature at sea level To = 288.15 K = 15°C. Humidity is set to 0 and CO2 content is set to standard 450 ppm, which is irrelevant for an accuracy of less than 0.5% as can be explored with the Calculator for Refractivity based on Ciddor Equation.

We can use the Calculator for Refractivity based on Ciddor Equation (or even the older Edlén Equation) to find No = 278. The Ciddor equation was derived empirically by accurate measurements of the refractive index of air under many different conditions [5].

Now we can calculate K1:

(24)

So we finally have refractivity of the atmosphere as a function of altitude:

(25)
where'
 $N(h)$ ' =' 'refractivity at altitude h $P(h)$ ' =' 'air pressure at altitude h in mbar $T(h)$ ' =' 'absolute temperature at altitude h

The Refractivity of air is proportional to air density ρ and P/T.

We can calculate the factor K1 also with the following equation, where the refractivity N can be calculated using the Calculator for Refractivity based on Ciddor Equation, for dry air with wavelength 550 nm and CO2 concentration 450 ppm at standard atmospheric conditions at sea level:

(26)

## Deriving Equation for Refraction

As we know how the Refraction Coefficient k can be calculated from the refractive index gradient dn/dh or the refractivity gradient dN/dh (20), we can derive the equation for the Refraction Coefficient from atmospheric parameters like pressure P, temperature T and temperature gradient dT/dh as shown in (7).

To be able to calculate the refractivity gradient dN/dh, we need equations to calculate pressure P and temperature T as functions of altitude h. This equations are provided for the model of the International Standard Atmosphere. For the troposphere with a linear temperature gradient α:

(27)
with
see
where'
 $P(h)$ ' =' 'Air pressure at altitude h in mbar or hPa or 1/100 Pa $P_{\mathrm{ref}}$ ' =' 'Air pressure at altitude href $T_{\mathrm{ref}}$ ' =' 'Absolute temperatur at href in K $h_{\mathrm{ref}}$ ' =' 'Reference altitude $\alpha$ ' =' 'Temperature Gradient in K/m $h$ ' =' 'Altitude above href $g$ ' =' '9.80665 m/s2 gravitational acceleration $R_\mathrm{S}$ ' =' 'Specific Gas Constant; dry air = 287.058 J/(kg·K)
(28)
where'
 $T(h)$ ' =' 'Absolute Temperature at altitude h in Kelvin $T_\mathrm{ref}$ ' =' 'Absolute Temperatur at href $\alpha$ ' =' 'Temperature Gradient $h$ ' =' 'Altitude above href $h_\mathrm{ref}$ ' =' 'Reference altitude

We can insert this equations in (22) and get refractivity as a function of altitude:

(29)
where'
 $\Delta h$ ' =' 'h − href

We need the refractivity gradient, so we have to calculate the derivative with respect to h:

 (30)

We can split the blue term into 2 terms:

 (31)

Applying this split and simplifying we get:

 (32)

Resubstituting the expressions for P and T:

 (33)

gives:

 (34)

α is simply the temperature gradient α = dT/dh, so we finally get:

(35)

We can insert some values: g/RS = 0.0343 K/m, K1 = 79.0 K/mbar to get finally the curvature of light rays. Note: I used a value g = 9.83 m/s2 to get the published constant 0.0343. This is slightly more than average gravitational acceleration. Gravitational acceleration on earth at sea level is in the range of g = 9.7639...9.8337 m/s2 [6].

#### Curvature of Light

(36)
where'
 $c$ ' =' 'Curvature of light rays $r$ ' =' '1/c = radius of curvature of light rays $P$ ' =' 'Air pressure at the observer in mbar or hPa or 1/100 Pa, Standard = 1013.25 mbar $T$ ' =' 'Temperature at the observer in Kelvin, Standard = 288.15 K = 15°C $\mathrm{d} T/\mathrm{d} h$ ' =' 'Temperature Gradient at the observer in K/m or °C/m, Standard = −0.0065°C/m

Note: In the famous formula for the refraction coefficient k (37) a value of R · K1 = 503 is used. To get this 503 I use a slightly smaller value K1 = 78.95. This is justified, because we dropped a 1/n = 1/1.000278 factor in (19). This way we get a better match between k and the curvature of light c when converting between them using R = 6,371,000.

The curvature of light is:

• proportional to the negative Refractive Index Gradient
• proportional to the negative Refractivity Gradient
• proportional to the negative Density Gradient
• most sensible to the Temperature Gradient
• independent on the shape of the earth

Note: the Temperature Gradient dT/dh is the most important factor that determines the curvature of light. At the same pressure and temperature you can get light bending up or down depending on the Temperature Gradient!

#### Refraction Coefficient Globe

We can also calculate the Refraction coefficient for the Globe model:

(37)
where'
 $k$ ' =' 'Refraction Coefficient $R$ ' =' '6371 km = Radius of the earth $c$ ' =' 'Curvature of light rays (36) $P$ ' =' 'Air pressure at the observer in mbar or hPa or 1/100 Pa, Standard = 1013.25 mbar $T$ ' =' 'Temperature at the observer in Kelvin, Standard = 288.15 K = 15°C $\mathrm{d} T/\mathrm{d} h$ ' =' 'Temperature Gradient at the observer in K/m or °C/m, Standard = −0.0065°C/m

## Influence of Water Vapor

The absolute value of the refractive index is irrelevant for the bending of light, only the gradient, the amount of change per height, is relevant for how much light gets bent.

(38)
where'
 $c$ ' =' 'curvature of light ray $r$ ' =' 'radius of curvature of light ray $\mathrm{d} n / \mathrm{d} h$ ' =' 'refractive index gradient, i.e. the change per altitude

If the refractive index gradient is zero, i.e. no change in refractive index with increasing altitude, then the curvature of the light ray is zero, no matter what the absolute value of the refractive index is.

The refractive index is proportional to the density of the air. So if the density is constant, then the refractive index is constant and we have no light bending in a constant density. Because at constant air temperature the density decreases with increasing altitude, the refractive index decreases proportionally and light gets bent down towards the denser part of the atmosphere, i.e. objects appear loomed up at standard conditions.

The density of moist air can be calculated as follows: [7]

(39)
where'
 $\rho_\mathrm{ma}$ ' =' 'density of moist air $\rho_\mathrm{da}$ ' =' 'density of dry air $K_\mathrm{w}(x)$ ' =' 'correction factor depending on humidity ratio $x$ ' =' 'humidity ratio $R_\mathrm{a}$ ' =' '287.1 J/kg/K = specific gas constant of dry air $R_\mathrm{w}$ ' =' '461.5 J/kg/K = specific gas constant of water vapor

The maximal humidity ratio for air at 15°C is x = 0.01062 [8]. The corresponding correction factor Kw due to maximum saturation humidity at 15°C is then:

(40)

So saturated moist air density is 0.9936 times the density of dry air. Moist air is maximal 0.64% less dense than dry air at 15°C. Can this small difference in density between dry and saturated air have any visible influence at all?

Lets investigate how water vapor can change the density gradient of air and what effect that has on observations. Some Flat Earthers claim, because the air above water contains more moisture than higher up, which means it has less density over water than higher up, and light bends from less dense to denser regions, the light will bend upwards over water and that's the reason why the flat earth appears like a Globe. Is that a physically possible scenario?

The following graphic shows how water vapor can change the density gradient or refractive index gradient respectively from the red curve (dry air) to the blue curve (moist to dry air). Remember, the slope of the curve is responsible for the bending of light (refraction). To get upwards bending, the curve must raise from the left to the right.

The air cannot have less humidity than 0. This limits the influence of water vapor on the refractive index as follows:

1. If humidity stays constant over great height ranges, then it does reduce the magnitude of density very slightly evenly, but does not change the density gradient, i.e. no change in the slope of the curve. So there is no change in the bending of light between dry and moist air.
2. If humidity reduces slowly over great height ranges, then there is a very, very slight change in the density gradient, the slope of the curve. The corresponding change in light bending is there but can be neglegted. Light still bends down.
3. To have any noticable influence on refraction at all, humidity has to change very quick. It has to go from maximum humidity to dry air in a small layer directly above the water to get a strong enough gradient change that has any effect on refraction at all. How fast moist air dries is limited by physics.
4. Normally when water is cooler than the air above light gets bent down stronger than at standard conditions, where the water and air have the same temperature. But when the sun is heating up the water, a small strong water vapor gradient can build, which can create a very small layer with inverted density. In this layer light gets bent upwards. This can cause inferior mirages even if the water is cooler than the air above. See also the excellent video Mirages require a steep temperature gradient... or do they? by ScienceItOut.

Generally the density gradient depends mostly on the temperature gradient, then on pressure and on absolute temperature and only very little on the humidity gradient. A constant humidity has no influence on refraction at all, a decreasing humidity with altitude makes refraction a tiny bit less, an increasing humidity a tiny bit stronger, very marginal, negligible in practice. Only a very strong vapor gradient over water can invert the density gradient enough to cause inferior mirages.

The surface layer is the only place where refraction can vary considerably due to steep temperature gradients due to heat exchange between ground and air. Above the surface layer density always decreases slowly according to the equations of the International Standard Atmosphere. See also Refraction Coefficient as a Function of Altitude.

So a small water vapor layer directly over water can explain some inferior mirages in cases where they are not present without water vapor. But water vapor can not bend light upwards in such a way that it could explain why the earth appears globe shaped if the earth would be flat. Strong enough water vapor layers dry out very quickly and are therefore restricted to small layers over water. Above this layer light will always bend down.

## Correcting for Refraction

If we know the Refraction Coefficient, height measurements using a theodolite can be corrected with the following approximation equations, as long as the distance d is much less than the radius of curvature of the light ray:

(41)
see
where'
 $\rho$ ' =' 'Refraction Angle in radian. Multiply with 180°/π to get degrees $c$ ' =' 'Curvature of light ray, see (5) $r$ ' =' 'radius of curvature of the light ray, see (5) $k$ ' =' 'Refraction Coefficient, see (7) $d$ ' =' 'target distance in m $R$ ' =' '6,371,000 m = radius of the earth

A positive Refraction Angle means the light ray is bent down and the object appears too high. In a zenith angle measurement we have to add the Refraction Angle to the zenith angle to get the correct angle as measured without refraction.

From the Refraction Angle we can calculate the magnitude l of how much an object at distance d appears to be raised due to Refraction:

(42)
where'
 $l$ ' =' 'apparent lift of object due to refraction $\rho$ ' =' 'Refraction Angle in radian, see (41) $c$ ' =' 'Curvature of light ray, see (5) $r$ ' =' 'radius of curvature of the light ray, see (5) $k$ ' =' 'Refraction Coefficient, see (7) $d$ ' =' 'target distance in m $R$ ' =' '6,371,000 m = radius of the earth

As Refraction Angles in radian are numbers much less than 1, this equations are very good approximations. The distance d can be the line of sight distance or the distance along the light ray or the distance along the surface of the earth. They are all essentially the same as long as d is much smaller than R or r respectively.

## References

Terrestrial refraction; Wikipedia
https://en.wikipedia.org/wiki/Atmospheric_refraction#Terrestrial_refraction
Speed of Light; Wikipedia
https://en.wikipedia.org/wiki/Speed_of_light
Atmospheric Refraction; Radio and Microwave Wireless Systems; Prof. Sean Victor Hum
Very good derivation of atmospheric refraction due to refractive index gradients
http://www.waves.utoronto.ca/prof/svhum/ece422/notes/20a-atmospheric-refr.pdf
[4]
The Constants in the Equation for Atmospheric Refractive Index at Radio Frequencies; Ernest K. Smith, Jr., and Stanley Weintraub
http://web.archive.org/web/20180721151000/https://nvlpubs.nist.gov/nistpubs/jres/50/jresv50n1p39_A1b.pdf
Refractivity of Air; An Introduction to Green Flashes; Andrew T. Young
Chronoligy of early measurements of the refractive index of air
https://aty.sdsu.edu/explain/atmos_refr/air_refr.html
Gravity of Earth; Wikipedia
https://en.wikipedia.org/wiki/Gravity%5Fof%5FEarth
Density of Moist Humid Air; The Engineering ToolBox
https://www.engineeringtoolbox.com/density-air-d_680.html
Humidity Ratio of Air; The Engineering ToolBox
https://www.engineeringtoolbox.com/humidity-ratio-air-d_686.html
[9]
An Accurate Method for Computing Atmospheric Refraction
Stone, R. C.; Journal: Publications of the Astronomical Society of the Pacific, v.108, p.1051-1058; Bibliographic Code: 1996PASP..108.1051S