# Calculator for RHi and Contrail Persistence Criteria

Discussion in 'Tools for Investigating and Debunking' started by Mick West, Jan 12, 2016.

Tags:

I've started on a calculator to convert RHw (relative humidity with respect to water) to RHi (relative humidity with respect to ice).

https://www.metabunk.org/rh

Currently it just takes RHw, temperature and pressure, and then gives you the RHi, and the value of RHw at which RHi is 100% (i.e. the value at which contrails will persist)

Note the pressure has very little effect, the relationship is more one of temperature.

I will expand this calculator to make it more useful

Here's the code I use to convert RHW to RHI based on pressure and temp:
Code:
```# convert RHw to RHi, as per
# ta = temperature of atmosphere, in Celsius
# p = pressure in hPa (hectopascals), which are the same as mb (millibars)
# note this is defined as accurate from -0C to -50C, but it's unclear how inaccurate it is
# at temperatures below -50C
def rhw2rhi(rh, ta, p)
e_sat_w_Ta = (1.0007 + 0.00000346 * p)*6.1121*Math.exp((17.966*ta)/(247.15+ta));
e_sat_i_Ta = (1.0003 + 0.00000418 * p)*6.1115*Math.exp((22.452*ta)/(272.55+ta));
return rh * (e_sat_w_Ta/e_sat_i_Ta);
end
```
From:

Note there seem to be several different ways of calculating the saturation vapor pressure.

The equations above are known as the Arden Buck equations, and come from Buck 1981 - New Equations for Computing Vapor Pressure and Enhancement Factor. They are what as know as "empirical" equations, in that they are simply designed to most closely fit the numbers observed in experiments, and the equations themselves do not necessarily have any really relationship to what is going on in the real physical world. The paper details how the experiments were carried out, and how the equations were worked out.

Essentially they are an approximation (over particular ranges of T) of formulae given by Wexler,

(Where θ is the temperature in Kelvin and exp and ln are natural exponent and logarithms)

Wexler's formulae in turn are derived from a variety of measurement made in the 1970s, after the older (and less accurate) Goff-Gratch equations were derived.

https://en.wikipedia.org/wiki/Goff–Gratch_equation

Several other different approximations are found here:
http://faculty.eas.ualberta.ca/jdwilson/EAS372_13/Vomel_CIRES_satvpformulae.html

But ultimately these are all just approximation to actual real world data. You could equally take the value from a table (interpolating for intermediate values, fractions of degrees), and get the same, or more accurate result, and computationally it would probably be faster.

#### Attached Files:

• ###### Buck 1981 - New Equations for Computing Vapor Pressure and Enhancement Factor.pdf
File size:
350.4 KB
Views:
532
Last edited: Jan 19, 2016
• Like x 11
• Agree x 1
• Useful x 1

I've made it embeddable with [rhcalc].[/rhcalc], and added an interactive graph. You can drag the (Temperatre,RHw) around and see the RHi value move. The two lines are the ice supersaturation lines for RHi (the striaght line, just RHi = 100) and RHw (the curve, which is RHw when RHi = 100).

See OP for embedded calculator.

3. ### JFDeeSenior Member

For descriptive reasons, it would be great to display altitude (standard atmosphere) related to the selected pressure - probably as a separate text label or even an editable text field (with a disclaimer attached).

But the thing is, pressure does not make any significant different to the relationship between RHw and RHi. If you put in 0 or 1000 it only makes a different of about 0.05% in the calculated RHi values.

I believe this is to do with the very low compressibility of liquid water. Note the very small multipliers for pressure

Last edited: Jan 19, 2016
5. ### JFDeeSenior Member

Even so, there is a scale and an entry field for pressure values. It would just be a tad less abstract when it's illustrated that pressure is directly related to altitude.

Oops, that scale should be RH%, not mb. Fixed.

7. ### JFDeeSenior Member

Ah, that got me confused ...

Given the slope of the line is pretty flat from -50 to -40°C, we can derive a simple formula for converting RHw to RHi that works for most contrail conditions.

RHi = RHw * (0.89 - (0.0148 * t))

This is accurate to about 0.1% between -50 and -40, and 1% between -70 and -30.

While using the full calculator would be best, this still might be useful if you want to provide a simple answer to the question "how do you convert RHw to RHi?"

• Like x 2

And perhaps an even more useful question that we can derived a heuristic for is "What relative humidity is required for contrail persistence". Normally this is given as either "above 70%", or the somewhat more accurate "between 60% and 70%". But if we rearrange the above simple formula

RHi = RHw * (0.89 - (0.0148 * t))

to express RHw in terms of RHi:

RHw = RHi / ( 0.89 - (0.0148 * t) )

The add in the condition for persistence (RHi >= 100)

We get

RHw >= 100/( 0.89 - (0.0148 * t))
(Accurate within 1% for -70°C < t < -30°C)

• Like x 3

I modified the calculator to use the more accurate Wexler formula

Only makes about 0.5% difference over useful range, but better to have the correct formula. Code might be useful for someone:
Code:
```	// These are the Wexler functions as described in Buck 1981
function e_sat_w_wex(t,p) {
// Wexler uses t in Kelvin, so convert from C
// p is ignored
t += 273.15;

return 0.01 * exp( -2991.2729*pow(t,-2)
- 6017.0128*pow(t,-1) + 18.87643854
-0.028354721*t + 0.17838301*0.0001*pow(t,2)
-0.84150417 * 0.000000001*pow(t,3) + 0.44412543*0.000000000001*pow(t,4)
+ 2.858487 * Math.log(t));
}

function e_sat_i_wex(t,p) {
// Wexler uses t in Kelvin, so convert from C
// p is ignored
t += 273.15;
return 0.01 * exp(-5865.369*pow(t,-1)
+ 22.241033 + 0.013749042*t
- 0.34031775 * 0.0001 * pow(t,2) + 0.26967687*0.0000001*pow(t,3)
+ 0.6918651 * Math.log(t));
}
```

The above is focussed on calculating the point at which air is ice-supersaturated. i.e. the point at which contrails persist because the ice crystals cannot sublimate back into vapor.

I think it would be useful to develop interactive versions of all the graphs that are commonly used when explaining contrail formation and persistence. In particular the Appleman chart and the mixing line chart.

@Ross Marsden, you generated mixing line charts before:
Do you still have the code and/or parameters used for this? In particular the "Jet Engine" point (off the chart). The blue lines are the equilibrium vapor pressure for water and ice, for which I can use the Wexler or Buck equations.

@skephu, in the video you posted, there are some nice illustrations of the mixing line.

Do you know how these illustrations were generated? Is there code?

And then the Appleman chart.

Does anyone happen to have the actual equation for the curves there and how the contrail factor is included.

Then there's newer versions, as in this post:
Unfortunately the links are broken, @CbIncus, I know you've not been back since 2013, but do you still have the code?

12. ### skephuSenior Member

Last edited: Jan 20, 2016
• Useful x 1

Images sometimes take a while to upload after they are pasted, so if you post too quick they don't show up.

Were the values you use for the initial temperature and mixing ratio (mass ratio of water to air) just for illustration, or were any of them real numbers from somewhere?

14. ### skephuSenior Member

Temperatures for a high-bypass engine were taken from here:
Kärcher B, Hirschberg MM, Fabian P. Small‐scale chemical evolution of aircraft exhaust species at cruising altitudes. Journal of Geophysical Research: Atmospheres (1984–2012). 1996 Jun 27;101(D10):15169-90.

Table 3 contains temperatures, fuel flow, air flow, etc. (sorry, pasting images just doesn't work for me).
But for the graph you included above, the temperatures are increased a bit because it shows the case where the efficiencies of a high-bypass and a non-bypass engine are identical.
For the non-bypass engine, the temperature was calculated from the high-bypass data assuming that all the heat heats the core flow.
I don't remember now where the initial mixing ratio came from but that is also based on real data. Basically it can be calculated from the fuel flow, emission index of water, and air flow.

• Useful x 1
15. ### C_JNew Member

Thanks for the calculator. I just want to make sure I'm understanding, and sorry if this seems ignorant or is a stupid question: when I enter certain values and then hit "calculate", what's the green line with the red dot? Is that the persistent contrail line? In other words, it has to be above that in order for there to be persistent contrails?

Essentially yes, it's the RHw value at that temperature for which RHi=100.

The dashed curved line uses the full Wexler formula. The green line uses a simple derived function:

RHi = RHw * (0.89 - (0.0148 * t))

Intended to be accurate over the range -70C to -30C (and more so from -60 to -40), and useful for manual calculation.

Last edited: Apr 24, 2016

Note you can drag around the blue "RH Water" point, which might help clarify. Drag it along the green line, and notice that RHi stays around 100%

18. ### C_JNew Member

Ok, thanks. So when I plug in my numbers from today: RHw 41, pressure 300, temp -42.5, it says the Rhw has to be above 65.8% for persistent contrails. However, there were persistent contrails all afternoon. Are there other factors that are involved? I don't know how there were persistent contrails today.

The problem is that it's very hard to get an accurate RHw for a particular point in space in time. The weather balloons are widely spaced, and not very accurate. See discussion here:

20. ### C_JNew Member

Well the one I found was close to where I am, and even if the humidity and temperature varied quite a bit, it should still be below the persistent chemtrail line.

Nonetheless, how DO you get an accurate RHw?

Go up there with a cooled mirror hygrometer at the same time the contrail forms.

There are more accurate forecast models than single nearby sounding, but they are still just computed values. Unless you go measure it, you don't know for sure what it is.

22. ### C_JNew Member

Ok, so basically there's no way to figure out the persistent contrail issue? That's frustrating that we can't confirm one way or another. It's interesting, though, that yesterday the readings were almost exactly the same, however, there were NO visible contrails in the sky all day when just a few days ago with the same readings, there were long huge persistent ones. Just doesn't make sense.

What issue?

24. ### TrailblazerModeratorStaff Member

Well, the best indicator that there is RHi above 100% is... persistent contrails! They're a far more accurate indicator than sounding balloons.

And yes, I am fully aware of the circular logic there.

• Agree x 3
25. ### TrailspotterSenior Member

Try earth.nullschool.net server to see the predicted RH values at the altitude of 250hPa in your area at the time. They correlate well with the occurrences of persistent contrails at least in Europe and North America. Here, for example, the "prediction" for California on April 20, 2016 at 21:00 UTC, which was a good day for contrails:

It correlates well with the satellite images of the day:

You can use earth.nullschool.net for real predictions, checking the predicted RH values for the next day or two, it generally works. I did it here a few times, the predicted values tend to change upon recalculation with the most recent data, but these changes are small. As in the weather forecast, the predicted persistent contrails may appear a bit earlier, or a bit later than the initial prediction.

Last edited: Apr 26, 2016
26. ### C_JNew Member

Yes, circular logic, haha. When I see hundreds of persistent contrails on one day and then literally none the next day, and the weather was the same and the humidity was the same, it's curious, that's all.

27. ### LandruModeratorStaff Member

How did you determine the weather and humidity were the same at altitude for both days?

28. ### TrailblazerModeratorStaff Member

Rather than using soundings, which are rather inaccurate and widely spaced (both temporally and geographically), have you tried using weather model data? For instance over the past couple of days there have been no contrails at all visible here in the southern UK, because the air is extremely dry up there:

However, by tomorrow there will be much more humidity around, so I would expect to see the contrails back again.

Maps like this are a much easier way to visualise the changing conditions than individual soundings are.

29. ### C_JNew Member

Does it show the temp and humidity for altitudes that jets fly?

31. ### TrailblazerModeratorStaff Member

Yes. The maps that I posted there are from publically available GFS data (the weather model run by NOAA). These ones give relative humidity at the 300 millibar pressure level, which equates to about 30,000 feet. Temperature is also available. See instantweathermaps.com for charts.

Obviously this is not foolproof, as contrails can also form lower than that, and much higher than that. For example I recently saw contrails forming at about 26,000 feet, even though the 300mb humidity was seemingly too low. In fact the tropopause was below 30,000ft on that occasion, so the air above was much drier, and all persistent contrails were below 30,000ft.

See https://www.metabunk.org/photos-of-...t-you-took-yourself.t1487/page-21#post-179910 and the following post.

32. ### C_JNew Member

Thanks, but I can't seem to find any maps like you're talking about unless I become a paid subscriber. I'm clicking MAPS, then GFS then US REGIONAL, and it takes me to a page to pay.

33. ### skephuSenior Member

Choose Maps, GFS, CONUS, Upper Air, then Relative Humidity at 300mb, at 0h.

34. ### C_JNew Member

Thanks. And what's the 0z, 6z, 12z, 18z?
And what's th 0h, 3h, 6h, etc.?

GFS is a forecast based on data taken at a certain time. the 0z, 6z etc. are Zulu Time for the base data of the forecast.

0h, 3h etc are forecasts relative to that base time.

Example: the most recent one is 2016-05-08 12Z - that's 12 noon (12PM) in Zulu/UTC/GMT, or 5AM PDT, May 9th

Here I've chosen the 24h forecast, so it's giving me the forecast for 12Z (5AM PDT) 24 hours after the base time, i.e. on May 9th (see the two date/times above the map)

36. ### funghi2New Member

guys can you expand on RHi vs RHw? how would ice even have an RH lower than 100 if its ice? little confused about this

37. ### TrailblazerModeratorStaff Member

Both RHi and RHw are a measure of the relative humidity of the air.

They values differ because they are measured relative to different things.

RHi is the humidity relative to the saturation vapour pressure over ice.

RHw is the humidity relative to the saturation vapour pressure over water.

The saturation vapour pressure basically means the humidity that you will find, at equilibrium, over a sheet of ice (for RHi) or a pool of water (for RHw) at the given temperature and pressure.

Water has a higher vapour pressure than ice (i.e. water vapour will come off water more easily than off ice), so it stands to reason that the humidity of a given volume of air will be lower when measured relative to water (RHw) than when measured relative to ice (RHi).

So a given humidity can be expressed as greater than 100% when measured as RHi, but less than 100% when measured as RHw.

There's a rather more technical explanation here: http://www.rhs.com/papers/RH_WMO.pdf

• Like x 1
• Informative x 1

And the "over water/ice" can be read quite literally - i.e. it's related to the situation where the water vapor (in the air) is right next to some liquid or solid water.

In a lab this would be "over" something like a beaker of water.

In the air it's over (around) other liquid micro drops or ice crystals, or an equivalent nucleation seed.

If it's greater than 100% for that substance, then it will condense on that substance.

• Like x 1
39. ### TrailblazerModeratorStaff Member

Yes, another way of thinking about it is that water vapour in the air is more prone to condense onto ice than it is onto water. So as the air gets more humid, the humidity will hit 100% with respect to ice before it hits 100% with respect to water.

40. ### funghi2New Member

i see someone asked but i still dont get the instant weathermap....if i just want to know what the map is at this moment, what z and what h do i select?