Détermination de la base des cumulus

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

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.

Historique

Explication qualitative

À la fin du Modèle:S, la détermination de la base des cumulus était mal comprise. Ainsi Alfred Angot affirma : Modèle:Citation bloc 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

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. La hauteur de la base des nuages est exprimée comme suit:

$b=k\times (T-D)$

où 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 (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 :

$b=k\times (T-D)$

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 $k=\frac{300\mathrm{v}\mathrm{e}\mathrm{r}\mathrm{g}\mathrm{e}\mathrm{s}}{4\mathrm{d}\mathrm{e}\mathrm{g}\mathrm{r}\acute{\mathrm{e}}\mathrm{s}\mathrm{F}\mathrm{a}\mathrm{h}\mathrm{r}\mathrm{e}\mathrm{n}\mathrm{h}\mathrm{e}\mathrm{i}\mathrm{t}}$.

On convertit maintenant en unités standards. On a: 1 verge = Modèle:Unité et 1 degré Fahrenheit = Modèle:Unité. Donc,

$k=\frac{300\times 0,9144}{4\times \frac{5}{9}}=123,4\mathrm{m}/\mathrm{K}$

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:

$\gamma =\frac{\frac{3}{2}\times \frac{5}{9}}{100\times 0,9144}=9,11\mathrm{K}/\mathrm{k}\mathrm{m}$

Le gradient adiabatique sec est $\gamma =g/c_p=9,75\mathrm{K}/\mathrm{k}\mathrm{m}$ où g = 9,8065 m/s² est l'accélération de la gravité et c_p=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.

Autres formules

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

$b=k\times (T-D)$ où k = Modèle:Unité / ⁰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 suivanteModèle:, :

$b=a\times (T-D)$ où a = Modèle:Unité / ⁰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)

Fichier:Emagramme altitude cumulus.svg

É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 $p_{H_{2}O}$ sur la pression atmosphérique totale. Soit z l'altitude. On a donc

$r(z)=\frac{p_{H_{2}O}(T)}{p_a(z)}$

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

$p_{H_{2}O}(T)=p_s(D)$

où $p_s(D)$ est la pression de vapeur saturante.

On différencie et donc:

$\frac{d{p}_{H_{2}O}(T)}{dz}dz=\frac{d{p}_s(D)}{dD}dD(1)$

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

On a donc:

$r(z+dz)=\frac{p_{H_{2}O}(z+dz)}{p_a(z+dz)}$

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

$\frac{p_{H_{2}O}(z+dz)}{p_a(z+dz)}=\frac{p_{H_{2}O}(z)}{p_a(z)}$

Donc,

$(p_{H_{2}O}(z)+\frac{d{p}_{H_{2}O}(z)}{dz}dz))p_a(z)=p_{H_{2}O}(z)(p_a(z)+\frac{d{p}_a}{dz}dz)$

Donc,

$\frac{d{p}_{H_{2}O}(z)}{dz}=\frac{p_{H_{2}O}(z)}{p_a(z)}\times \frac{d{p}_a}{dz}$

Dans l'article pression atmosphérique, il est démontré que $\frac{d{p}_a}{dz}=-\rho g$ où $\rho$ est la masse volumique de l'air. Donc,

$\frac{d{p}_{H_{2}O}(z)}{dz}=-\rho g\frac{p_{H_{2}O}(z)}{p_a(z)}$

On substitue dans (1) et donc :

$-\rho g\frac{p_{H_{2}O}(z)}{p_a(z)}dz=\frac{d{p}_s(D)}{dD}dD(2)$

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

$P=p_0\exp \left[\frac{L}{R/M}(\frac{1}{T_0}-\frac{1}{T})\right]$

On considère la pression de référence p₀ = 1013.15 hPa et la température de référence T₀ = 373.15 K. $L=2.470\times 10^6$ 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, T₀ = 273.15 K et T est la température de l'atmosphère.

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

$p_s(D)=\exp (13.7-5120/D)$

On peut donc écrire

$p_s(D)=\hat{p}{e}^{-\frac{\hat{D}}{D}}$

avec $\hat{D}=\frac{LM}{R}=5350K$.

Dans la formule de Rankine, on a $\hat{D}=5120$ ce qui n'est pas applicable en l'espèce car on a en général $D\approx 300$

On obtient donc:

$\frac{d{p}_s(D)}{dD}=\frac{\hat{D}}{D^2}\hat{p}{e}^{-\frac{\hat{D}}{D}}=\frac{\hat{D}}{D^2}{p}_s$

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

$-\rho g\frac{p_{H_{2}O}(z)}{p_a(z)}dz=\frac{\hat{D}}{D^2}{p}_{s}dD(2)$

Donc,

$\frac{dD}{dz}=-\rho g\frac{p_{H_{2}O}(z)}{p_s(D)p_a(z)}\times \frac{D^2}{\hat{D}}$

On simplifie un peu et on suppose que $p_{H_{2}O}(T)\approx p_s(D)$. On obtient alors:

$\frac{dD}{dz}=-\rho g\frac{1}{p_a(z)}\times \frac{D^2}{\hat{D}}$

On a vu que $\hat{D}=5350K$, $p_a=1.01\times 10^{5}Pa$, $\rho g=11.77Pa/m$. Donc, pour D = 300 K, on obtient

$\frac{dD}{dz}=-\frac{300\times 300}{5350}\times \frac{11.77}{1.01\times 10^5}=-1.86\times 10^{-3}K/m$

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 $\frac{dT}{dz}=-\frac{g}{C_p}=-\frac{9.81}{1006}=-0.00975K/m$.

On obtient donc:

$\frac{d(T-D)}{dz}=-0.00186+0.00975$

Donc en mètres, $z=127\times (T-D)$. 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 $dD/dz=0.00174$ et le facteur 127 est remplacé par 125. On retrouve exactement le coefficient proposé par Stull.

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 = Modèle:Unité. La formule approchée ci-dessus donnait une base de nuages à Modèle:Unité (i.e. Modèle:Unité).

Remarque d'utilisation

Le calcul ci-dessus est fait pour une pression atmosphérique normale de Modèle:Unité. 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 Modèle:Unité 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 = Modèle:Unité (ou h = 30 × 125 = Modèle:Unité). Cependant, la lecture de l'émagramme donne une hauteur approximative de Modèle:Unité. L'émagramme ci-joint provient du livre de Stull 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, 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 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 Modèle:Unité).

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 ou des 125 mètres par kelvin de HennigModèle:, qui permettent d'évaluer la hauteur des nuages à moins de 5% près.

Bibliographie

• Modèle:Ouvrage
• Modèle:Ouvrage

Lien externe

• Modèle:Ouvrage