Détermination de la base des cumulus

(Redirigé depuis Formule de Hennig)

La détermination de la base des cumulus, en présence de bulles de convection, est le calcul de la hauteur à laquelle commence ce type de nuages lors d'un réchauffement de la surface terrestre et avec une humidité donnée. Elle peut être déterminée avec une grande précision (de l'ordre de quelques pourcents) en application des lois de la thermodynamique.

Cumulus humilis à base bien définie qui dépend de la température et de l'humidité disponible.

Historique

modifier

Explication qualitative

modifier

À la fin du XIXe siècle, la détermination de la base des cumulus était mal comprise. Ainsi Alfred Angot affirma :

« La hauteur des nuages éprouve dans le cours de l'année une variation très nette : une même couche de nuages est en moyenne à une altitude plus grande en été qu'en hiver. Cela se comprend aisément, car l'air en même temps qu'il est plus chaud, est généralement plus sec en été qu'en hiver; il faut donc s'élever davantage pour rencontrer une couche assez froide où commence la condensation[1]. »

La notion d'air sec est ambigüe car en général l'air contient plus de vapeur d'eau en été même si l'humidité relative est plus basse. L'auteur ne connaissait pas la formule de Hennig découverte 4 ans plus tôt et ébauchée par James Espy 60 ans plus tôt. L'explication énonçant « rentrer une couche assez froide » est vague par nature. En fait le calcul de la base des nuages convectifs est extrêmement précis.

Découverte de la formule

modifier

Dans le monde germanique, il est dit que la formule donnant la base des nuages a été publiée pour la première fois par Richard Hennig en 1895[2]. La hauteur de la base des nuages est exprimée comme suit:

 

b est la hauteur de la base des nuages, T est la température et D est le point de rosée. avec k = 122,6 m/K.

Cependant, la détermination de la hauteur de la base des cumulus a été pour la première fois évoquée par James Espy en 1841[3] (soit 54 ans plus tôt) qui affirma que si la différence entre la température et le point de rosée est de 4 degrés Fahrenheit, alors la base des nuages sera à 300 verges de hauteur. Exprimée de manière moderne, la formule d'Espy est la suivante :

 

où à nouveau, b est la hauteur de la base des nuages, T est la température et D est le point de rosée.

Espy affirme donc que  .

On convertit maintenant en unités standards. On a: 1 verge = 0,914 4 m et 1 degré Fahrenheit = 100180 K. Donc,

 

Aussi surprenant que cela puisse paraître, le résultat est très proche de la valeur admise de 125 m/K. Cependant le raisonnement utilisé par l'auteur est hautement discutable. En effet, l'auteur affirme que le gradient adiabatique sec γ est 100 verges par 1,5 degré Fahrenheit. Donc, si l'on convertit en unités standards, l'on obtient:

 

Le gradient adiabatique sec est  g = 9,8065 m/s² est l'accélération de la gravité et cp=1006 J/(kg K) est la capacité calorifique de l'air à pression constante. On constate donc que l'estimation par Espy du gradient adiabatique sec est minorée de 6,5 %. Ensuite, il utilise un argument peu clair où il relie la décroissance moyenne de la température en fonction de la hauteur qui serait de 1 degré Fahrenheit par 100 verges avec le gradient adiabatique sec qui est de 1.5 degré Fahrenheit. Il semble faire une moyenne arithmétique où le coefficient serait 100 verges/1,25 degré Fahrenheit. Il obtiendrait alors 131 m/K. Cependant l'explication n'est pas satisfaisante[Note 1].

Autres formules

modifier

La FAA propose une formule qui sous-estime la hauteur des cumulus de l'ordre de 5%. Cette formule est la suivante[5] :

 k = 400 pieds / ⁰C;

T est la température au sol et D est le point de rosée.

Une formule plus précise (formule de Hennig) est la suivante[6],[7] :

 a = 0,125 km / ⁰C;

Le calcul du coefficient k (ou a) peut être effectué rigoureusement en appliquant les lois de la thermodynamique et en supposant que la bulle d'air ne se mélange pas avec l'air extérieur. On notera que le coefficient k varie de quelques pourcents en fonction de la température au sol.

Calcul rigoureux du coefficient k (ou a)

modifier
 
Émagramme (voir la discussion ci-dessous) qui montre que T soulevé le long de l'adiabatique sèche et à rapport de mélange constant nous permet de trouver le niveau de condensation par convection et donc la hauteur de la base d'un cumulus

On suppose que la parcelle d'air ne se mélange pas avec l'air extérieur. On définit le rapport de mélange r(z) comme le rapport de la pression partielle de vapeur d'eau   sur la pression atmosphérique totale. Soit z l'altitude. On a donc

 

Le point de rosée est défini par la température D telle que

 

  est la pression de vapeur saturante.

On différencie et donc:

 

On élève la parcelle de l'altitude z à l'altitude z + d z.

On a donc:

 

Comme il n'y a pas de mélange, on a r(z+ dz) = r(z). Donc,

 

Donc,

 

Donc,

 

Dans l'article pression atmosphérique, il est démontré que    est la masse volumique de l'air. Donc,

 

On substitue dans (1) et donc :

 

La pression de vapeur saturante de l'eau peut s'exprimer à partir de la formule de Clapeyron comme suit[8] :

 

On considère la pression de référence p0 = 1013.15 hPa et la température de référence T0 = 373.15 K.   est la chaleur latente de vaporisation (à 20 C), R = 8.3144621 J /mol /K est la constante universelle des gaz parfaits, M = 0.01801 kg/mol est la masse molaire de l'eau, T0 = 273.15 K et T est la température de l'atmosphère.

La pression de vapeur saturante de l'eau (en)   est aussi donnée par la formule de Rankine[9] (qui diffère légèrement) :

 

On peut donc écrire

 

avec  .

Dans la formule de Rankine, on a   ce qui n'est pas applicable en l'espèce car on a en général  

On obtient donc:

 

On substitue dans l'équation (2) et donc :

 

Donc,

 

On simplifie un peu et on suppose que  . On obtient alors:

 

On a vu que  ,  ,  . Donc, pour D = 300 K, on obtient

 

L'expression ci-dessus donne le gradient du point de rosée dans une parcelle en ascension en fonction de l'altitude.

Le gradient adiabatique sec est  .

On obtient donc:

 

Donc en mètres,  . Ces conditions correspondent à une journée typique de vol à voile où T = 27 ⁰C. Aux conditions standard de température et de pression, on a   et le facteur 127 est remplacé par 125. On retrouve exactement le coefficient proposé par Stull[6].

Par conséquent, pour une température de 27 ⁰C et un point de rosée de 12 ⁰C la base des nuages sera à 15 × 127 = 1 905 m. La formule approchée ci-dessus donnait une base de nuages à 6 000 pieds (i.e. 1 800 mètres).

Remarque d'utilisation

modifier

Le calcul ci-dessus est fait pour une pression atmosphérique normale de 1 013,15 hPa. Dans la réalité, la pression au sol varie d'heure en heure et de jour en jour avec le passage des systèmes météorologiques. Il faut donc ajouter une correction correspondant à la pression réelle.

Par exemple, dans l'émagramme ci-joint, la pression au sol est de 1 000 hPa et la hauteur déduite à partir de celui-ci est donc la hauteur au-dessus de ce niveau de pression. En appliquant la formule ci-dessus, la base des cumulus se situerait à h = 30 × 127 = 3 809 mètres (ou h = 30 × 125 = 3 750 mètres). Cependant, la lecture de l'émagramme donne une hauteur approximative de 3 500 mètres. L'émagramme ci-joint provient du livre de Stull[10] et est semble-t-il légèrement différent de l'émagramme original. Vu que l'échelle est logarithmique, il est difficile de déterminer la hauteur exacte du cumulus. En outre, en examinant soigneusement l'émagramme original de Stull[10], on constate que l'origine de la courbe de rapport de mélange est à 0.7 ⁰C et non à 0 ⁰C ce qui introduit une erreur systématique de 100 mètres environ. Stull a publié un autre émagramme[11] qui permet de trancher la question. Dans ce graphe, on considère 2 courbes : la droite verticale correspondant à la température potentielle de 30⁰C et la courbe de rapport de mélange constant ayant pour origine -3 ⁰C. Ces courbes se coupent à 4100 mètres d'altitude environ. Donc, si l'on soustrait 3 × 125 = 375 mètres, la hauteur de la base du cumulus sera 3725 mètres environ ce qui est proche de la valeur théorique (3 750 ou 3 809 mètres).

L'exploitation graphique d'un émagramme ou l'utilisation de la formule longue est donc délicate. Il est plus simple d'utiliser la règle des 400 pieds par kelvin de la FAA[5] ou des 125 mètres par kelvin de Hennig[6],[7] qui permettent d'évaluer la hauteur des nuages à moins de 5% près.

Notes et références

modifier
  1. Bohren[4] fait une lecture incorrecte du livre d'Espy où Bohren fait dire à Espy que le coefficient b vaudrait 100 verges / degré Fahrenheit (or celui-ci parle de 300 verges / 4 degrés Fahrenheit). En plus, Bohren affirme que « this rule of thumb seems to not to be widely known » (Traduction: « La formule de Hennig ne semble pas être largement connue »)[4] ce qui est très discutable car c'est ce que la FAA enseigne[5] !

Références

modifier
  1. Alfred Angot, Traité élémentaire de météorologie, Gauthier-Villars, , 417 p. (lire en ligne), p. 203
  2. (de) Richard Hennig, « Eine einfache Formel, die ungefähre Höhe der Wolkenbildung bei adiabatischen Zuständen zu bestimmen. », Meteorologische Zeitschrift, vol. 12,‎ , p. 129 (lire en ligne)
  3. (en) James Espy, Philosophy of Storms, Charles Little &, , 552 p. (lire en ligne), p. 4
  4. a et b (en) Craig F. Bohren et Bruce A. Albrecht, Atmospheric Thermodynamics, Oxford University Press, , 402 p. (ISBN 978-0-19-509904-1), p. 275
  5. a b et c (en) anonyme, Pilot Handbook of Aeronautical Knowledge Chapitre 11, Federal Aviation Administration (lire en ligne), p. 11-14
  6. a b et c Meteorology for Scientists and Engineers, p. 101
  7. a et b (en) Roger P. Frey, Whether the Weather : Aviation Meteorology from A to Z, Books on Demand, , 180 p. (ISBN 978-3-7386-1574-6, lire en ligne), p. 166
  8. Meteorology for Scientists and Engineers, p. 96
  9. Nicolas Vandewalle, « Chapitre 7: Changement d'état », Université de Liège (consulté le )
  10. a et b Meteorology for Scientists and Engineers, p. 121
  11. Meteorology for Scientists and Engineers, p. 122

Bibliographie

modifier
  • (en) M. K. Yau et R. R. Rogers, Short Course in Cloud Physics, Butterworth-Heinemann, , 3e éd. (ISBN 0-7506-3215-1, EAN 9780750632157), p. 39
  • [Meteorology for Scientists and Engineers] (en) Roland B. Stull, Meteorology for Scientists and Engineers Second Edition, Brooks/Cole, , 502 p. (ISBN 978-0-534-37214-9)

Lien externe

modifier
  • Thermodynamique avancée Cours 11 : La stabilité verticale : Stabilité latente ou convective, UQAM (lire en ligne [PDF])