4.4.1 Fonction de poids pour la diffusion par une surface rugueuse

Bien que l’intégrale des coefficients bistatiques donnée en (4.104) et (4.108) soit une solution moins coûteuse en temps de calcul qu’une solution rigoureuse des équations de Maxwell, car elle est issue d’une méthode d’approximation, elle n’en reste pas moins très lourde numériquement. Johnson et Zhang ([45 ]) ont proposé une réécriture de (4.104) et (4.108) qui a l’avantage d’une part d’accélérer considérablement le calcul et d’améliorer la précision, et d’autre part de permettre de mieux apréhender l’influence relative des différentes longueurs d’onde du spectre des vagues dans le processus de diffusion. Le détail du calcule est explicité dans [45], je ne donne dans cette section que les grandes lignes de ce calcul et les résultats importants. En premier lieu, les expressions (4.104) et (4.108) sont rassemblées sous une intégrale unique et RPE s’exprime désormais sous la forme

   |_      _|                           |_             '  '   _| 
    Rv      integral   oo        integral  2p           gv(er,he,fe,kr',f')
RPE=  |_  Rh  _| +    k20k'rdk'r    Y(k'r,f'). |_  gh(er,he,fe,kr',f') _|  df',
     0      kd          0             g3(er,he,fe,kr',f')
     0                                g4(er,he,fe,kr,f )
(4.137)

où les fonctions gv,gh,g3 et g4 sont données par 4.133. Les coefficients de reflexion Ri et Rc donnés en (4.104) et (4.108) sont très proches en valeur absolue et sont de signe opposé. Ri est positif (augmente la reflectivité) et Rc est négative. La somme SR = Ri + Rc, qui caractérise l’effet total de la diffusion, résulte donc de la différence DR = |Ri|-|Rc| |Ri| -~ |Rc|. La différence |Ri|-|Rc| va ainsi être faible devant |Ri | et |Rc|, c’est à dire

DR   «   |Ri| ,                             (4.138)
DR   «   |Rc| .                             (4.139)

Ainsi, fixer une précision relative donnée sur l’intégrale numérique de DR est plus efficace que de fixer cette même précision relative sur les deux intégrales de Ri et Rc effectuée séparement avant de faire la différence |Ri | -|Rc|.

En developpant les fonctions gv,gh,g3 et g4 en série de fourier de l’angle f', on peut démontrer ([45] que R PE s’exprime sous la forme d’une série de Fourier

RPE = RPE,0 + RPE,2 cos(2f0)
(4.140)

R PE,0 est la réflectivité omnidirectionnelle et RPE,2 est l’amplitude de la variation azimutale de la réflectivité. Ces coefficient sont donnés par les relations suivantes

                     |_               _| 
         integral              g'v,0(er,he,q)
R     =    oo  C (k').  g'h,0(er,he,q),  dq,
 PE,0    kd/k0  0 r   |_  g'3,0(er,he,q)  _| 
                       g'4,0(er,he,q)
(4.141)

et

                     |_               _| 
         integral              g'v,2(er,he,q)
R     =    oo  C (k').  g'h,2(er,he,q),  dq,
 PE,2    kd/k0  2 r   |_  g'3,2(er,he,q)  _| 
                       g'4,2(er,he,q)
(4.142)

q = kr '/k0, est le rapport entre la longueur d’onde d’une vague et la longueur d’onde du radiomètre, C0 (kr') et C2(kr') sont respectivement le spectre de courbure (C(kr',f')) omnidirectionnel et l’amplitude de sa variation azimutale. Ainsi, la composante omnidirectionnelle de la reflectivité des petites échelles (RPE,0) est obtenue en intégrant sur les échelles comprises entre kd et l’infini le spectre de courbure omnidirectionnel pondéré par la fonction g'x,0 (où x = v,h,3 ou 4 selon la polarisation). De même, l’amplitude de la composante azimutale de la reflectivité des petites échelles (R PE,2 ) est obtenue en intégrant sur les échelles comprises entre kd et l’infini la seconde harmonique du spectre de courbure pondérée par la fonction g'x,2 (où x = v,h,3 ou 4 selon la polarisation).

Les fonctions de poids g'x,0,, sont illustrées sur les figures 4.22, 4.23 et 4.24 respectivement pour des incidence de 0o, 30o et 60o. On peut déduire de ces figures qu’au nadir (figure 4.22), la réflectivité d’une surface rugueuse (PE) est essentiellement affectée par les vagues dont la longueur d’onde est de l’ordre de c0. Lorsque l’incidence augmente (figures 4.22, 4.23 et 4.24), la gamme des longueurs d’ondes affectant sensiblement la réflectivité s’élargie autour de c0 .


PIC 
PIC
(a)  
(b)
PIC 
PIC
(c)  
(d)
PIC 
PIC
(e)  
(f)
 

FIG. 4.22: Fonctions de poids (harmoniques 0 et 2 en trait plein et tirets) en valeur absolue (à gauche) et signes des fonctions de poids (à droite) à h = 0o, SST = 15oC et SSS = 35 psu (modèle KS pour la constante diélectrique).



PIC 
PIC
(a)  
(b)
PIC 
PIC
(c)  
(d)
PIC 
PIC
(e)  
(f)
PIC 
PIC
(g)  
(h)
 

FIG. 4.23: Fonctions de poids (harmoniques 0 et 2 en trait plein et tirets) en valeur absolue (à gauche) et signes des fonctions de poids (à droite) à h = 30o, SST = 15oC et SSS = 35 psu (modèle KS pour la constante diélectrique).



PIC 
PIC
(a)  
(b)
PIC 
PIC
(c)  
(d)
PIC 
PIC
(e)  
(f)
PIC 
PIC
(g)  
(h)
 

FIG. 4.24: Fonctions de poids (harmoniques 0 et 2 en trait plein et tirets) en valeur absolue (à gauche) et signes des fonctions de poids (à droite) à h = 60o, SST = 15oC et SSS = 35 psu (modèle KS pour la constante diélectrique).