The Humidity Resource
This resource is compiled from two of my posts (#8 and #10) in a 2006 thread entitled:
Relative Humidity <> Absolute Humidity
It is intended as a resource and tutorial for any "moist air" problems involving water vapor pressure, relative humidity, specific humidity, absolute humidity, mixing ratio, etc. This resource contains a few minor updates from the original posts, and a section on station pressure following two worked examples.
For anyone looking for a chart for solving humidity problems, Google "psychrometric chart." There are free ones out there, I recommend the Carrier chart if you can find it.
For those who want formulas for programs, I will try. Most of this is pretty straightforward with the exception of what equation to use for saturated water vapor pressure.
Relative humidity is a ratio of actual water vapor pressure to saturated water vapor pressure. Saturated vapor pressure is the maximum that can exist at a given temperature; excess will condense out as dew or frost. If unsaturated air is cooled, at some temperature, it will become saturated and dew will form. This dewpoint is another way to describe the actual water vapor pressure.
Let ea = actual water vapor pressure, es = saturated water vapor pressure, RH = relative humidity, T = temperature, and Td = dewpoint temperature, then
RH = ea(T)/es(T) or ea(T) = RH * es(T)
ea(T) = es(Td)
RH = es(Td)/es(T)
If one knows the shape of the saturated vapor pressure curve, actual water vapor pressure and relative humidity can be calculated from temperature and dewpoint. More on the saturated vapor pressure calculation in a moment.
If ea(T) is known, other humidity properties of moist air can be calculated assuming ideal gas law. If total pressure is P, and water vapor pressure is ea, then pressure of dry air is P-ea. Mixing ratio, as kg of water vapor/kg of dry air, is given by
MR = 0.622*ea/(P-ea)
the constant is the ratio of molecular weights, water/air.
Specific humidity, kg of water vapor per kg of moist air mixture, is
MR/(1+MR).
We all know the ideal gas law, PV = nRT. But the number of moles, n, can be written as m/M where m is the mass of gas, and M is the molecular weight in g/mol. Then PV = (m/M)*RT.
Rearranging, the density m/V = P*M/(T*R). The absolute humidity is just the density of water vapor,
ea*M/(T*R)
where M = 18.02 g/mol, and R = 8.314472 Pa·m³/(mol·K), and T is absolute temperature in kelvins.
The key to all this is the es(T) curve. There are MANY equations to calculate or approximate it. See notes below on models and approximations. However, the approximation used by the National Weather Service for humidity calculations in surface observations is one due to Bolton,
es(t) = 611.2 Pa * exp(17.67*t/(243.5 + t))
where 611.2 Pa is the saturated vapor pressure at 0 °C, and
t is temperature in °C, not K.
(this fits steam table data to about 0.4% at 50 °C, and better at lower temperatures)
With this equation:
*Calculate es at T and Td
*Calculate relative humidity
*Calculate mixing ratio, specific humidity, or absolute humidity as desired from above.
Models and Approximations
As mentioned, there are MANY formulas to calculate es(T), and that makes this subject confusing.
I use "models" for formulations that have a thermodynamic basis and fit data well over a wide temperature range. The model formulas tend to be rather complex. Calculation is acceptable for a table (one time), but not for on-the-fly calculations in automated stations. "Approximations" are empirical fits (usually to a model) over a limited temperature range, but are faster to calculate. Also most approximations are much easier to "reverse," that is calculate the temperature at which a given saturation pressure occurs.
Whether models or approximations, there are always different formulas for vapor pressure over water and over ice, and for the models, both fit solid experimental data. Unfortunately, below 0 °C, meteorologists like to calculate relative humidity over supercooled water, not ice; however, scant data exists to justify the extrapolations they make to formulas only validated for water above 0 °C. The models extrapolate quite differently below 0 °C. (I have no clue which, if any, is right.)
I won't list all the formulas, but you can Google terms below, along with water vapor pressure and find them. A quick survey of models:
*Goff-Gratch: Originally 1945, updated in 1946, 1957, 1965. The standard of WMO (WMO No.49) and basis of Smithsonian Meteorological Tables.
*Hyland-Wexler: 1983, used by ASHRAE as the standard in heating/air conditioning calculations for several years (now they use a steam table formulation). Several weather equipment manufacturers use Hyland-Wexler or approximations to it.
*Sonntag, 1994 (can't find out much but a mention, and the formula)
*IAPWS Formulation, 1997: These are the guys who develop official Steam Tables for power plants, etc. Goff-Gratch and Hyland-Wexler claim validity to 100 °C. IAPWS is valid to the critical point of water (374 °C). Only it should be used above 100 °C
NOTE: See update on models in Post #4
A quick survey of approximations:
*The original Magnus equation dates to 1844; the slightly modified form above was by Murray, 1967, with the coefficients by Bolton, 1980. Many people have fit the Magnus form to different temperature ranges by tweaking coefficients. It fits well over about a 50 °C range, generally.
*Bogel equation: Bogel added a fourth parameter to Magnus equation which improves fit and range.
*Flatau, et al, developed 6th and 8th order polynomial fits. The 6th order fits -50 °C to +50 °C to about 0.01%, the 8th order fits -85 °C to +70 °C, but fit is about 0.1%, near the maximum temperature.
Arden Buck has a good 1981 paper on fitting Magnus and Bogel coefficients to different temperature ranges. It is available on the web as a pdf. The Magnus and Bogel forms can explicitly be solved in reverse; given actual water vapor pressure, dewpoint can be calculated. The higher-order Flatau method requires numerical iteration to get dewpoint.
Examples
I was asked to add one or more examples to demonstrate use of the above equations and added two. The first is based on data from Pontiac, MI, a NWS office (and airport) in the northern Detroit suburbs, and closest NWS office to me.
Their 2:00 pm report (2008-03-26) gave temp 43 °F, dew point 25 °F, RH = 49%. Because they are 298 m above sea level, "normal" barometric pressure is about 97.8 kPa, today's is about 98.3 kPa. (That computation needs to be the subject of another post; this one will be complex enough).
NWS uses a dog's breakfast of units in their reports. If you want to get agreement with their figures, you have to know which are source data and which are rounded. For temperature use the °F data. It needs to be converted to °C and you need two or more decimals because water vapor pressure is a strong function of temperature. 43/25 convert to 6.1111 °C/ -3.8889 °C. (The data is NOT that accurate, but it is that precise in calculation.)
Evaluating ea from the dewpoint, 0.4588 kPa.
Evaluating es at temperature, 0.9420 kPa, so RH is 48.7%, which agrees with published number of 49%
The mixing ratio is 0.622*(0.4588 kPa)/(98.3 kPa - 0.4588 kPa) = 0.00292 kg H2O/kg dry air
SH = MR/(1+MR)= 0.00291 kg H2O/kg moist air
AH (remember T is in kelvins) is ea*M/(T*R) = 3.56 g/m³ H2O.
For second example, we'll use some place warm, specifically Miami, where it is drier based on RH but wetter on any absolute measure. Temperature is 76 °F, dew point 54 °F, RH = 46%.
Converting the temperatures, we get 24.444 °C, 12.222 °C
ea(dew point) = 1.422 kPa
es(temp) = 3.0639 kPa, RH = 46.4%, again agrees with published figure.
Miami airport is 2 m above sea level. We ignore altitude correction, but local pressure is 102.4 kPa.(Pressure of the dry air is 102.4 - 1.4 kPa = 101 kPa, approx.)
Mixing ratio, using ea, is 0.00876 kg H2O/kg dry air
SH is 0.00868 kg H2O/kg moist air
AH is 10.36 g/m³
Well, I am over allowed length. The station pressure section will be post #2.
This resource is compiled from two of my posts (#8 and #10) in a 2006 thread entitled:
Relative Humidity <> Absolute Humidity
It is intended as a resource and tutorial for any "moist air" problems involving water vapor pressure, relative humidity, specific humidity, absolute humidity, mixing ratio, etc. This resource contains a few minor updates from the original posts, and a section on station pressure following two worked examples.
For anyone looking for a chart for solving humidity problems, Google "psychrometric chart." There are free ones out there, I recommend the Carrier chart if you can find it.
For those who want formulas for programs, I will try. Most of this is pretty straightforward with the exception of what equation to use for saturated water vapor pressure.
Relative humidity is a ratio of actual water vapor pressure to saturated water vapor pressure. Saturated vapor pressure is the maximum that can exist at a given temperature; excess will condense out as dew or frost. If unsaturated air is cooled, at some temperature, it will become saturated and dew will form. This dewpoint is another way to describe the actual water vapor pressure.
Let ea = actual water vapor pressure, es = saturated water vapor pressure, RH = relative humidity, T = temperature, and Td = dewpoint temperature, then
RH = ea(T)/es(T) or ea(T) = RH * es(T)
ea(T) = es(Td)
RH = es(Td)/es(T)
If one knows the shape of the saturated vapor pressure curve, actual water vapor pressure and relative humidity can be calculated from temperature and dewpoint. More on the saturated vapor pressure calculation in a moment.
If ea(T) is known, other humidity properties of moist air can be calculated assuming ideal gas law. If total pressure is P, and water vapor pressure is ea, then pressure of dry air is P-ea. Mixing ratio, as kg of water vapor/kg of dry air, is given by
MR = 0.622*ea/(P-ea)
the constant is the ratio of molecular weights, water/air.
Specific humidity, kg of water vapor per kg of moist air mixture, is
MR/(1+MR).
We all know the ideal gas law, PV = nRT. But the number of moles, n, can be written as m/M where m is the mass of gas, and M is the molecular weight in g/mol. Then PV = (m/M)*RT.
Rearranging, the density m/V = P*M/(T*R). The absolute humidity is just the density of water vapor,
ea*M/(T*R)
where M = 18.02 g/mol, and R = 8.314472 Pa·m³/(mol·K), and T is absolute temperature in kelvins.
The key to all this is the es(T) curve. There are MANY equations to calculate or approximate it. See notes below on models and approximations. However, the approximation used by the National Weather Service for humidity calculations in surface observations is one due to Bolton,
es(t) = 611.2 Pa * exp(17.67*t/(243.5 + t))
where 611.2 Pa is the saturated vapor pressure at 0 °C, and
t is temperature in °C, not K.
(this fits steam table data to about 0.4% at 50 °C, and better at lower temperatures)
With this equation:
*Calculate es at T and Td
*Calculate relative humidity
*Calculate mixing ratio, specific humidity, or absolute humidity as desired from above.
Models and Approximations
As mentioned, there are MANY formulas to calculate es(T), and that makes this subject confusing.
I use "models" for formulations that have a thermodynamic basis and fit data well over a wide temperature range. The model formulas tend to be rather complex. Calculation is acceptable for a table (one time), but not for on-the-fly calculations in automated stations. "Approximations" are empirical fits (usually to a model) over a limited temperature range, but are faster to calculate. Also most approximations are much easier to "reverse," that is calculate the temperature at which a given saturation pressure occurs.
Whether models or approximations, there are always different formulas for vapor pressure over water and over ice, and for the models, both fit solid experimental data. Unfortunately, below 0 °C, meteorologists like to calculate relative humidity over supercooled water, not ice; however, scant data exists to justify the extrapolations they make to formulas only validated for water above 0 °C. The models extrapolate quite differently below 0 °C. (I have no clue which, if any, is right.)
I won't list all the formulas, but you can Google terms below, along with water vapor pressure and find them. A quick survey of models:
*Goff-Gratch: Originally 1945, updated in 1946, 1957, 1965. The standard of WMO (WMO No.49) and basis of Smithsonian Meteorological Tables.
*Hyland-Wexler: 1983, used by ASHRAE as the standard in heating/air conditioning calculations for several years (now they use a steam table formulation). Several weather equipment manufacturers use Hyland-Wexler or approximations to it.
*Sonntag, 1994 (can't find out much but a mention, and the formula)
*IAPWS Formulation, 1997: These are the guys who develop official Steam Tables for power plants, etc. Goff-Gratch and Hyland-Wexler claim validity to 100 °C. IAPWS is valid to the critical point of water (374 °C). Only it should be used above 100 °C
NOTE: See update on models in Post #4
A quick survey of approximations:
*The original Magnus equation dates to 1844; the slightly modified form above was by Murray, 1967, with the coefficients by Bolton, 1980. Many people have fit the Magnus form to different temperature ranges by tweaking coefficients. It fits well over about a 50 °C range, generally.
*Bogel equation: Bogel added a fourth parameter to Magnus equation which improves fit and range.
*Flatau, et al, developed 6th and 8th order polynomial fits. The 6th order fits -50 °C to +50 °C to about 0.01%, the 8th order fits -85 °C to +70 °C, but fit is about 0.1%, near the maximum temperature.
Arden Buck has a good 1981 paper on fitting Magnus and Bogel coefficients to different temperature ranges. It is available on the web as a pdf. The Magnus and Bogel forms can explicitly be solved in reverse; given actual water vapor pressure, dewpoint can be calculated. The higher-order Flatau method requires numerical iteration to get dewpoint.
Examples
I was asked to add one or more examples to demonstrate use of the above equations and added two. The first is based on data from Pontiac, MI, a NWS office (and airport) in the northern Detroit suburbs, and closest NWS office to me.
Their 2:00 pm report (2008-03-26) gave temp 43 °F, dew point 25 °F, RH = 49%. Because they are 298 m above sea level, "normal" barometric pressure is about 97.8 kPa, today's is about 98.3 kPa. (That computation needs to be the subject of another post; this one will be complex enough).
NWS uses a dog's breakfast of units in their reports. If you want to get agreement with their figures, you have to know which are source data and which are rounded. For temperature use the °F data. It needs to be converted to °C and you need two or more decimals because water vapor pressure is a strong function of temperature. 43/25 convert to 6.1111 °C/ -3.8889 °C. (The data is NOT that accurate, but it is that precise in calculation.)
Evaluating ea from the dewpoint, 0.4588 kPa.
Evaluating es at temperature, 0.9420 kPa, so RH is 48.7%, which agrees with published number of 49%
The mixing ratio is 0.622*(0.4588 kPa)/(98.3 kPa - 0.4588 kPa) = 0.00292 kg H2O/kg dry air
SH = MR/(1+MR)= 0.00291 kg H2O/kg moist air
AH (remember T is in kelvins) is ea*M/(T*R) = 3.56 g/m³ H2O.
For second example, we'll use some place warm, specifically Miami, where it is drier based on RH but wetter on any absolute measure. Temperature is 76 °F, dew point 54 °F, RH = 46%.
Converting the temperatures, we get 24.444 °C, 12.222 °C
ea(dew point) = 1.422 kPa
es(temp) = 3.0639 kPa, RH = 46.4%, again agrees with published figure.
Miami airport is 2 m above sea level. We ignore altitude correction, but local pressure is 102.4 kPa.(Pressure of the dry air is 102.4 - 1.4 kPa = 101 kPa, approx.)
Mixing ratio, using ea, is 0.00876 kg H2O/kg dry air
SH is 0.00868 kg H2O/kg moist air
AH is 10.36 g/m³
Well, I am over allowed length. The station pressure section will be post #2.
Comment