FR 2933203 A1 20100101 FR 0803634 A 20080627 PROCEDE DE CARACTERISATION DE L'INTENSITE DE CONVECTION D'UN NUAGE, PAR UN RADAR METEOROLOGIQUE La présente invention concerne un procédé de caractérisation de l'intensité de convection d'un nuage par un radar météorologique. Elle s'applique notamment à des radars météorologiques embarqués sur des avions. En cours de vol, le pilote d'un avion utilise généralement un système d'observations météorologiques issues d'un radar lui donnant des informations sur l'état de l'atmosphère dans le but de détecter les risques inhérents aux situations météorologiques permettant entre autre de restituer un champ de vitesse Doppler et un champ de réflectivité radar en trois dimensions, 3D. Ces champs sont obtenus en effectuant avec le radar des balayages successifs suivant différents angles de site et d'azimut. En chaque point d'une grille 3D couvrant le domaine d'observation, on dispose ainsi en particulier d'une information de réflectivité radar dépendant de la longueur d'onde utilisée et des propriétés de réflexion, d'absorption et de diffusion des cibles présentes dans le volume de résolution radar. Sur l'écran de contrôle, dans le poste de pilotage, le pilote voit alors s'afficher une information météorologique analysée généralement selon quatre niveaux de couleurs correspondant soit à l'absence de signal, soit à un facteur de réflectivité faible, modéré ou fort. Le pilote interprète ces informations en termes de risque pour décider s'il peut conserver la trajectoire prévue pour son avion ou s'il doit la modifier pour éviter une zone présentant des risques météorologiques, par exemple une zone où un orage se produit. Pour simplifier la démarche du pilote dans sa prise de décision, il est utile de 25 traduire les champs météorologiques observés par le radar en un champ décrivant le niveau de risque météorologique. Généralement, les radars météorologiques embarqués sur les avions fonctionnent en bande X. La réflectivité radar mesurée par un radar météorologique en bande X utilisant une seule longueur d'onde, présente 30 une ambiguïté en ce qui concerne sa relation avec la nature et les caractéristiques des hydrométéores, telles que des gouttes d'eau, de neige ou de grêle notamment, qui l'ont générée. En effet, cette relation est fortement non bijective. Par conséquent un même niveau de réflectivité peut être généré par des hydrométéores de natures très différentes, provenant soit de la neige ou de la grêle notamment, correspondant à des situations météorologiques également très différentes du point de vue du risque aéronautique, par exemple des nuages stratiformes et des nuages de convection épaisse tels que les cumulonimbus. Ce défaut de bijectivité concerne notamment la grêle composée de grêlons de taille supérieure à environ un centimètre de diamètre sphérique équivalent. Ces grêlons sont situés hors du domaine de diffusion de Rayleigh, dans lequel travaillent usuellement les radars météorologiques. Ils se trouvent alors dans le domaine de diffusion de Mie qui est tel que lorsque la taille du grêlon augmente, sa réflectivité diminue. II est donc nécessaire de lever cette ambiguïté et d'obtenir avec un radar météorologique à la fois des informations sur la nature dynamique des nuages générateurs du champ de réflectivité, nuages stratiformes ou nuages convectifs intenses de type cumulonimbus d'orage, et sur la nature microphysique des hydrométéores, issus de pluie, de neige ou de gros grêlons. Des solutions sont connues pour lever cette ambiguïté. En particulier on sait déterminer le niveau de risque en chaque point de grille en utilisant de simples seuils sur le facteur de réflectivité radar observé en ces points. Si celui-ci est supérieur à 40 dBZ, le risque est considéré comme fort tandis que s'il est inférieur à 30 dBZ le risque est considéré comme faible. Entre ces deux valeurs, le risque est considéré comme modéré. Ces seuils peuvent éventuellement être différents selon les régions survolés conformément notamment à ce qui est décrit dans le brevet US 7 129 885. Cette méthode d'estimation du risque par la seule valeur ponctuelle de la réflectivité radar a notamment pour inconvénient qu'elle peut conduire à une sous-estimation ou une surestimation du risque dans certaines situations météorologiques. Par exemple, une forte averse de neige dans un système stratiforme peut être vue par un radar avec des facteurs de réflectivité élevés, supérieurs à 40 dBZ, alors que la situation ne présente pas de risque pour l'avion en vol. Un but de l'invention est notamment d'éviter cet inconvénient et de lever 35 l'ambiguïté précitée de façon fiable. A cet effet, l'invention a pour objet un procédé de caractérisation de l'intensité de convection d'un nuage, la réflectivité dudit nuage à une onde électromagnétique étant distribuée dans l'espace, la distribution de la réflectivité étant discrétisée selon un réseau de points (i, j, k) de l'espace en trois dimensions, au moins un profil est défini comme une fonction d'un paramètre, lui-même fonction numérique donnée à deux dimensions (i, j) de la distribution de la réflectivité en chaque point du réseau, ladite fonction variant uniformément entre une valeur constante minimale et une valeur constante maximale, la fonction étant égale à la valeur minimale lorsque le paramètre est inférieur à un seuil bas (Seuil min) et étant égale à la valeur maximale lorsque le paramètre est supérieur à un seuil haut (Seuil max), le nuage étant caractérisé comme convectif lorsque le profil est égal à l'une des valeurs constante et comme stratiforme lorsqu'il est égal à l'autre valeur constante. Avantageusement, la fonction du paramètre est par exemple normalisée. La fonction numérique à deux dimensions (i, j), définissant un paramètre, est par exemple calculée dans le plan horizontal. Un nuage convectif correspond par exemple à la valeur constante haute. Le profil peut être un profil global combinant plusieurs profils, ledit profil global étant défini par une fonction normalisée d'un paramètre combinant les différents paramètres numériques liés aux dits profils. La combinaison des différents paramètres est par exemple une moyenne pondérée. Un paramètre est par exemple l'intégrale de la réflectivité selon la dimension verticale (z, k). Un paramètre peut être une mesure de la hauteur de colonne, selon la dimension verticale, pour laquelle les points pris en compte sont ceux dont la valeur de réflectivité est supérieure à un seuil donné. Un paramètre peut aussi être l'altitude maximum où les réflectivités sont supérieures à un seuil donné. Un paramètre est par exemple la réflectivité maximum calculée dans une verticale donnée. Un paramètre est par exemple le gradient vertical de réflectivité entre une zone de réflectivité maximum et l'altitude maximum où les réflectivités sont supérieures à un seuil donné. Un paramètre peut aussi être le gradient horizontal de la réflectivité à une altitude donnée. Le paramètre est par exemple calculé dans un volume déterminé. D'autres caractéristiques et avantages de l'invention apparaîtront à l'aide de la figure qui suit faite en regard de dessins annexés qui représentent : la figure 1, une illustration d'un réseau de réflectivité radar à l'intérieur d'un nuage ; la figure 2, un profil caractérisant un nuage, utilisé par un procédé selon l'invention ; - la figure 3, un exemple de combinaison possible de profils caractérisant un nuage, dans un procédé selon l'invention. La figure 1 illustre une grille de réflectivité dans l'espace. L'espace est discrétisé selon un réseau de points de l'espace, en trois dimensions. Ce réseau forme une structure appelée grille par la suite. Chaque point de la grille, repéré par son indice i selon une ligne I, son indice j selon une ligne J et son indice k selon une ligne K présente une réflectivité Z(i, j, k) qui est détectée par un radar météorologique. Les indices i, j définissent la position d'un point dans le plan horizontal, et l'indice k défini la position verticale de ce point. L'espace ainsi discrétisé est mémorisé dans le calculateur du radar avec les réflectivités associées Z(i, j, k). Une matrice à trois dimensions est ainsi stockée que l'on confondra par la suite avec la grille de l'espace. On pourra ainsi de façon indifférente parler de grille de réflectivité Z(i, j, k) ou de matrice de réflectivité Z(i, j, k). La matrice de réflectivité peut être renouvelée ou mise à jour au fur et à mesure des périodes d'observation du radar. La matrice de réflectivité correspond à un domaine d'observation couvert par le radar. Le procédé selon l'invention permet de caractériser l'intensité de convection d'un nuage à partir des la distribution spatiale de sa réflectivité, plus particulièrement en fonction de la distribution spatiale discrétisée, selon une grille par exemple. En fonction de cette caractérisation, il est ensuite possible d'estimer, par un pilote ou de façon automatique au moyen d'un programme d'analyse, les risques liés à ce nuage. La figure 2 illustre par une courbe 1 une fonction de normalisation définissant un type de profil selon la valeur d'un paramètre considéré, ce profil défini par une courbe fonction de ce paramètre est du type à deux dimensions 2D. L'invention part du fait que pour estimer le niveau de risque en un point de la grille de réflectivité, la valeur ponctuelle du facteur de réflectivité apporte une réponse ambigüe donc insuffisante. Pour compléter l'information relative au seuil de réflectivité, l'invention considère un ensemble de paramètres relatifs à la distribution spatiale de réflectivité permettant de déterminer le caractère convectif épais ou stratiforme d'un nuage. La connaissance de ce caractère permet non seulement de connaître la nature physique des hydrométéores responsables d'un signal reçu par un radar mais aussi de pondérer correctement le risque lié à la dynamique des nuages par l'identification du caractère convectif ou stratiforme. En ce qui concerne l'identification de profils, plusieurs critères doivent être pris en compte pour obtenir ce profil. Ces critères sont ensuite exploités, par exemple à l'aide d'un arbre de décision, ou encore en intervenant comme les termes d'une somme pondérée. Dans ce dernier cas, les critères étant de natures différentes, ils sont par exemple normalisés, par exemple de façon que chacun d'eux varie entre la valeur 0 et la valeur 1 comme illustré sur la figure 2. Il s'agit donc de choisir un seuil haut , Seuil max sur la figure 1, un seuil bas, Seuil min sur la figure 2, et une fonction de normalisation pour chaque paramètre. La fonction peut être linéaire comme dans la figure 1, variant donc linéairement entre Seuil min et Seuil max. La fonction peut néanmoins présenter une autre forme. Ces seuils sont choisis tels qu'ils correspondent, suivant le critère adopté, pour l'un à une situation convective intense, cas du seuil haut, et pou l'autre à une situation stratiforme, cas du seuil bas. Ainsi, après normalisation du critère, un profil convectif intense obtient le score 1 tandis qu'un profil stratiforme obtient le score O. Si le profil est ambigu ou de convection peu intense, il obtient un score intermédiaire défini par la partie linéaire de la courbe. Les différents critères normalisés sont pondérés en fonction de la confiance que l'on accorde à chacun. Le score final est comparé à un seuil afin de déterminer la nature du profil. La figure 3 illustre par un schéma le principe d'identification de profils stratiformes ou convectifs utilisé par un procédé selon l'invention pour estimer des risques météorologiques. A partir du champ de réflectivité 3D 20, et pour chaque critère d'identification on calcule le paramètre correspondant 21 auquel on applique une fonction de normalisation du type de celle illustrée par la figure 2. On obtient ainsi autant de champs de profils 22, en 2D, que de critères choisis. L'exemple de la figure 2 comporte n critères. Dans un cas limite, il est possible de n'utiliser qu'un seul critère, n = 1 dans ce cas. Les champs de profils sont par exemple par la suite combinés pour une fonction 23 qui peut être une moyenne pondérée afin d'obtenir un seul champ de profil 24 utilisable pour la détermination du type microphysique des hydrométéores et du niveau de risque encouru par l'avion. Plusieurs critères d'identification peuvent donc être utilisés. Comme indiqué précédemment, un radar météorologique permet de restituer un champ de vitesse Doppler et champs de réflectivité en trois dimensions, 3D. En chaque point d'une grille 3D, préalablement définie et telle qu'illustré par exemple par la figure 1, et couvrant un domaine d'observation, on dispose ainsi d'une information de réflectivité. Les caractéristiques microphysiques des hydrométéores sont à déterminer pour chaque point de la grille mais le profil stratiforme/convectif en peut être identifié en ne considérant qu'un seul point. L'invention propose plusieurs paramètres qui permettent l'identification des profils. Ils sont calculés à partir des réflectivités mesurées en tous points de grille à la verticale du point considéré. A partir du champ 3D de réflectivité fourni par le radar, et éventuellement interpolé dans un repère approprié, par exemple le repère géocentré latitude ù longitude ù altitude, le procédé selon l'invention calcule un champ 2D de valeurs pour chaque paramètre d'identification. Un premier critère possible, appelé VIZ, est l'intégrale de la réflectivité le long de la verticale, en un point donné selon l'horizontale. II est représentatif de la quantité des précipitations en ce point donné. VIZ est déterminé en calculant l'intégrale de la réflectivité Z selon la verticale :35 =m~ VIZ(x, y) = J Z(x, y, z)dz 0 où x et y désignent les coordonnées cartésiennes du point selon l'horizontale et z la coordonnée verticale, l'intégrale étant calculée jusqu'à une 5 coordonnée maximale zmax le long de la verticale. Etant donné que les valeurs de réflectivité selon la verticale sont données dans la grille 3D de manière discrète, le calcul numérique du champ 2D de VIZ est donc donné par la relation suivante : kmax 10 VIZJ = 1Zy-k (2) k=0 où i et j représentent les indices selon les deux dimensions horizontales et k l'indice vertical, Z;~k est la réflectivité au point d'indices i, j et k de la grille de réflectivité. 15 Cette somme définie par la relation (2) peut être remplacée par la somme des intensités de précipitation selon la verticale. Un autre paramètre possible est la hauteur de colonne, notée Hco,onne. L'une des caractéristiques principales des zones convectives est leur extension 20 verticale importante. Le critère de hauteur de colonne représente la hauteur cumulée d'une colonne de valeurs de réflectivités supérieures à un seuil donné. Le critère de hauteur de colonne, en un point d'indices i, j selon l'horizontale, est ainsi donnée par la relation suivante : kmax Hcolonne = az;k 4A~ 25 k=0 1 1 si Z > seuil avec 8 = Z 0 sinon Dans cette relation (3), AA"k représente le pas en distance suivant l'axe vertical au niveau du point de la grille d'indices i, j et k, k étant l'indice selon la verticale, et kmax étant l'indice du point le plus haut. Z est la réflectivité au 30 point i, j, k.. 7 (1) (3) Un autre critère qui peut être utilisé par l'invention est complémentaire de la hauteur de colonne. Il sera appelé par la suite écho top et noté Etop. II permet d'améliorer la probabilité de détection dans certaines configurations, notamment pour les cumulonimbus en forme d'enclume. Il correspond à l'altitude de l'écho le plus haut sur une colonne d'observation détecté par le radar. Pour chaque colonne repérée dans le plan horizontale par un couple d'indices (i, j) dans la grille de réflectivité, on recherche l'élément d'indice vertical k le plus élevé pour lequel la réflectivité est supérieure à un seuil donné. En chaque couple (i, j), la valeur du champ d'écho top est l'altitude de l'élément trouvé, ainsi le critère Etop est donné par la relation suivante : Etopu = maxk (Auk , Zuk > 0) (4) où Auk désigne l'altitude d'un point de la grille de réflectivité et max k le maximum suivant k. Un autre critère possible est le facteur de réflectivité maximum rencontré dans une colonne, Zmax. Il permet de localiser rapidement les zones de précipitation intense, telle que la grêle, les fortes pluies ou la neige. Zmax, selon une colonne définie par ses indices i, j dans le plan horizontal, est obtenu selon la relation suivante : Zmax ~ =maxk (Z,~k ) (5) Un autre critère possible est le gradient vertical de réflectivité, gradä (Z). II s'agit du gradient vertical de réflectivité entre l'altitude correspondant au facteur Zmax et celle du sommet des échos, c'est-à-dire entre l'altitude zM où l'on observe le maximum de réflectivité et l'écho top Etop. Si le maximum de réflectivité s'étend sur plusieurs points de grille, l'altitude moyenne est alors considérée. Le gradient gradä (Z) est donné par la relation (7) issue de la relation (6) suivante : grade Z(x, y, z) _ Z(x, y, zM) ù Z(x, y, Etop) zM ù Etop (6) Dans le domaine numérique, en considérant les points discrétisés de la grille, en un point d'indices i, j, k : grade Z = Z(Auk,Zjk = Zmax) ù Z(A;, , A,jk = Etop) (7) Un autre critère possible est le gradient horizontal de réflectivité, gradH(Z). Il s'agit du gradient de réflectivité selon un plan horizontal à une altitude z choisie. En raison du caractère bidirectionnel de ce gradient, on considère par exemple sa norme, donnée par la relation suivante : IlgradH Z(x, y, z)II = z ùv2 ( aZ\2 'ôZ + ù \.ôx,Z ôyiZ (9) Du point de vue numérique, en un point i, j, k de la grille de réflectivité, le 15 gradient horizontal s'obtient selon la relation suivante : 2 2 Z(i+1) j,k=k, ù Zij,k=k, \ + Zi(j+I ),k=k ù Zij.k=k \ d. d . IJ,k=k, J )l,k=k, 20 en notant d; et di les pas selon les dimensions horizontales x et y respectivement. Un autre critère qui peut être utilisé pour l'identification d'un profil est la réflectivité Z1500 au de dessus de l'isotherme 0° C, à une altitude d'environ 25 1500 mètres. Z1500 = Z(i, j, k15oo ) où k15o0 est tel que Ajkisä =Ajkoc +1500, Alko°c étant l'altitude de 1/2 IlgradH Zuk Il ij,kk (10) 30 l'isotherme 0°C. Les critères ou paramètres précédents sont fondés sur l'analyse d'une colonne de réflectivité. Il est aussi possible selon l'invention d'utiliser des critères volumiques. En particulier, d'autres critères fondés sur des volumes différents peuvent être utilisés. Le champ de profil calculé peut ainsi devenir un champ 3D. Dans l'exemple d'une sphère, pour chaque point i, j, k de l'espace où l'on veut calculer un type de profil, on considère une sphère de rayon fixé centré sur ce point. Certains des paramètres précédents acquièrent alors une définition différente. Le critère VIZ devient un critère d'intégration dans l'espace qui représente alors la somme des réflectivités dans le volume de référence, la sphère dans le présent exemple. Le critère Zmax devient un facteur de réflectivité maximum observé dans ce volume de référence. D'autres critères volumiques peuvent être définis : Le volume Vzseuil total pour lequel le facteur de réflectivité Z est supérieur à un seuil ZSeU11 déterminé. Il s'agit de sommer les volumes élémentaires pour lesquels Z dépasse une valeur seuil Zseuii (par exemple 35 dBZ). Ce paramètre indique le volume ou la convection est parvenue à son terme de développement. Plusieurs seuils peuvent être envisagés ; Le volume continu Vcontinu, zseuil pour lequel le facteur de réflectivité Z est supérieur à un seuil déterminé Zseu;,. Le procédé est identique au calcul de Vzseuu, à ceci près que l'on doit assurer la continuité du volume calculé. Ainsi, pour être comptabilisé, chaque point de grille doit posséder au moins un certain nombre de voisins tel que Z > Zseuil. Ce nombre est fixé initialement et il est inférieur ou égal à 18 qui est le nombre total d'arêtes et de faces d'un cube. On peut, en outre, ajouter une condition sur la hauteur minimale à considérer et prendre par exemple l'altitude de l'isotherme 0 °C. Ce paramètre donne le volume de la zone génératrice de la précipitation convective ; La réflectivité moyenne Zmoy dans Vcontinu, zseuil, ce paramètre s'applique aussi aux autres volumes comme une colonne ou une sphère par exemple ; Le rapport des axes d'un ellipsoïde pour lequel Z est supérieur à un seuil Zseuil fixé. II s'agit de représenter au mieux, la figure géométrique 35 correspondant à Z > Zseui1 par un ellipsoïde dont les axes sont av selon la verticale et aH selon le plan horizontal. On forme alors le rapport av/aH et on fixe un seuil permettant de séparer un profil de type convectif d'un profil de type stratiforme. Il est également possible d'utiliser des grandeurs statistiques sur des 5 volumes, comme par exemple la variance spatiale de réflectivité. Cela indique si l'on est en présence d'un système homogène, composé d'une seule cellule, ou de plusieurs cellules. Pour chacun des paramètres d'identification précédemment définis, deux 10 seuils, haut et bas, sont nécessaires dans la fonction de normalisation pour calculer le profil, comme illustré notamment par la figure 2. Ces seuils sont par exemple choisis d'après l'analyse d'observations radar d'événements types, par exemple une tempête de neige ou de grêle caractéristique relevée en un point donné du globe. 15 Les caractéristiques des distributions des intensités de précipitation et donc de réflectivité peuvent différer selon la climatologie du lieu d'observation. II peut alors être nécessaire d'adapter les seuils suivant des considérations climatologiques. Par ailleurs, les seuils d'altitude, pour les critères écho-top et hauteur de 20 colonne, peuvent nécessiter une adaptation en fonction de la latitude de la zone observée. En effet, l'étagement des nuages en altitude est conditionné par la tropopause, dont l'altitude est fonction de la latitude du lieu. Comme l'illustre la figure 3, il est possible d'effectuer une identification 25 multicritères. L'identification d'un profil de convection peut être faite selon un seul critère mais on augmente significativement la fiabilité de cette identification en utilisant plusieurs critères. Un champ de profil 22 est calculé pour chacun des paramètres choisis 21. Les différents champs ainsi obtenus sont combinés pour obtenir un seul champ de profil 24. La combinaison de 30 ces champs peut être faite par le calcul d'une moyenne 23 avec des poids, éventuellement différents, associés aux critères d'identification. Ces poids peuvent également être adaptés suivant des considérations climatologiques comme pour les seuils des paramètres L'invention a notamment comme avantage de permettre l'extraction d'information d'intensité de convection à partir des informations de réflectivité, de manière efficace. L'invention s'applique aux observations de réflectivité effectuées par un radar météorologique, notamment aéroporté, fonctionnant en bande X. Elle peut aussi s'appliquer pour des radars fonctionnant dans d'autres bandes de fréquences, par exemple en bande W, Ka, C ou S, de type polarimétrique ou non polarimétrique, à effet Doppler ou non. 15 20 REVENDICATIONS1. Procédé de caractérisation de l'intensité de convection d'un nuage, la réflectivité dudit nuage à une onde électromagnétique étant distribuée dans l'espace, caractérisé en ce que la distribution de la réflectivité étant discrétisée selon un réseau de points (i, j, k) de l'espace en trois dimensions, au moins un profil (22) est défini comme une fonction (1) d'un paramètre (21), lui-même fonction numérique donnée à deux dimensions (i, j) de la distribution de la réflectivité en chaque point du réseau, ladite fonction variant uniformément entre une valeur constante minimale et une valeur constante maximale, la fonction étant égale à la valeur minimale lorsque le paramètre est inférieur à un seuil bas (Seuil min) et étant égale à la valeur maximale lorsque le paramètre est supérieur à un seuil haut (Seuil max), le nuage étant caractérisé comme convectif lorsque le profil est égal à l'une des valeurs constante et comme stratiforme lorsqu'il est égal à l'autre valeur constante. 2. Procédé selon la revendication 1, caractérisé en ce que la fonction (1) du paramètre (21) est normalisée. 3. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que la fonction numérique à deux dimensions (i, j), 20 définissant un paramètre, est calculée dans le plan horizontal. 4. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce qu'un nuage convectif correspond à la valeur constante haute. 25 5. Procédé selon l'une quelconque des revendications précédente, caractérisé en ce que le profil est un profil global (24) combinant plusieurs profils (22), ledit profil global étant défini par une fonction normalisée d'un paramètre combinant les différents paramètres numériques liés aux dits 30 profils (22). 6. Procédé selon la revendication 5, caractérisé en ce que la combinaison des différents paramètres (21, 22) est une moyenne pondérée. 7. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce qu'un paramètre est l'intégrale de la réflectivité selon la dimension verticale (z, k). 8. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce qu'un paramètre est une mesure de la hauteur de colonne, selon la dimension verticale, pour laquelle les points pris en compte sont ceux dont la valeur de réflectivité est supérieure à un seuil donné. 9. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce qu'un paramètre est l'altitude maximum où les réflectivités sont supérieures à un seuil donné. 10. Procédé selon l'une quelconque des revendications précédentes, 15 caractérisé en ce qu'un paramètre est la réflectivité maximum calculée dans une verticale donnée. 11. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce qu'un paramètre est le gradient vertical de réflectivité entre 20 une zone de réflectivité maximum et l'altitude maximum où les réflectivités sont supérieures à un seuil donné. 12. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce qu'un paramètre est le gradient horizontal de la réflectivité 25 à une altitude donnée. 13. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que le paramètre est calculé dans un volume déterminé. 30