Introduction aux méthodes de
télédétection des signaux
géophysiques
Quelles signaux
• Position
• Forme de la Terre ( altitude)
• Champs de potentiels
– Champ magnétique
– Champ de gravité
• Autres grandeurs
– Tsunami, Signaux sismiques
Principe de positionnement radio
• 3 inconnues (position)
• Chaque noeud transmet sa position et un bip synchronisé pour tous les noeud
• Réception simultanée des 3 signaux:
– mesure du temps de retard entre bip
• Résolution du système d’inconnue permettant d’expliquer ces temps de retard.
• Difficulté: 1m = 3.3 ns
– Vitesse de propagation de la lumière
– Synchronisation des nœuds difficile ( 400 m d’erreur)
Le système GPS : fonctionnement et applications
Le système NAVSTAR (NAVigation System by Timing And Ranging) - GPS (Global
Positionning System) est un système de positionnement par satellite qui a été développé et
est contrôlé par le département de la défense américain.
Conçu à l ’origine pour fournir aux forces armées américaines un système de repérage en tout
point du globe avec une très haute précision, le système GPS a depuis lors, trouvé des
applications dans le domaine civil.
Le système GPS répond aux questions suivantes :
Quelle est l’heure, quelle est la position, quelle est la vitesse ?
• Agriculture raisonnée
– Détermination du stress des cultures par télédétection
– Cartographie des zones à traiter ( engrais, irrigation)
– Pilotage par GPS des instruments d’épendage
Positionnement par satellites
Principe de GPS
• 4 inconnues (temps et position)
• Chaque satellite transmet
– Sa position, son temps local
• Réception simultanée des 4 signaux avec le temps de retard entre
satellites et station
• Résolution du système d’inconnue permettant d’expliquer ces
temps de retard.
• Difficulté: 1m = 3.3 ns
– Vitesse de propagation de la lumière
– Horloges atomiques de grande stabilité dans les satellites
LES TROIS SEGMENTS DU SYSTEME GPS
Le segment spatial comporte plus de 25 satellites répartis de façon à assurer depuis la Terre une
visibilité simultanée de 4 à 8 satellites. Les satellites sont répartis selon 6 plans orbitaux, ayant tous une
inclinaison d ’environ 55° par rapport à l ’équateur. Les satellites sont sur des orbites quasi-circulaires, à
une altitude d ’ environ 20184 km et leur période est de douze heures.
Chaque satellite est muni d ’un émetteur-récepteur, d ’une horloge atomique, d ’un ordinateur et
d ’équipements auxiliaires destinés au fonctionnement du système. L ’énergie pour l ’équipement est
fournie par 7 m2 de capteurs solaires tandis qu ’un système de fusées d ’appoint permet de réajuster la
position du satellite sur son orbite et de contrôler sa stabilité.
Le segment de contrôle est constitué de plusieurs stations au sol qui vérifient en permanence le bon
fonctionnement des satellites.
Le segment utilisateur est constitué de récepteurs qui permettent de convertir les données provenant
des satellites en informations sur la position, la vitesse et la date.
1st Generation: Test GPS Satellites (BLOCK I) 1978-85
• The first GPS ----! February 22, 1978
2nd Generation: Operating GPS Sat. (BLOCK II, IIA) 1989-94
• The first GPS ------! February 14, 1989
3rd Generation: Improved GPS Sat. (BLOCK IIR) 1990s-
• 24 satellites
GPS Type Number Year Characters
1st Gen BLOCK I 11 1978-85 Test
2nd Gen BLOCK II, IIA 28 1989-94 Operating
3rd Gen BLOCK IIR 20 1990s- Advanced
EXO
Trouver la durée de visibilité d’un satellite GPS?
Détail du codage 1/2
• f0 = 10.23 Mhz (signal horloge atomique)
• f1= 154 f0 = 1575.42 Mhz
• f2= 120 f0 = 1227.60 Mhz
• P-code (f0) , C/A-code (f0 /10), W-code (f0 /20), D-code
(50 Hz)
L1(t) = a1P(t)W (t)D(t) cos( f1t) + a1CA(t)D(t) sin( f1t)
L2 (t) = a1P(t)W (t)D(t) cos( f2 t)
Remarque sur
l’effet Doppler:
!
De nombreux effets relativistes
EXO
EXO
Equations de Maxwell
r
r r r %E
".B = 0 "#B = µ0 ( j + $ 0 )
%t
r r
# r %B
".E = "#E =&
$0 %t
• Plasma neutre !=0
• Equation du mouvement des électrons
r r r r
"v
me! e = #e( E + v $ B)
! "t
r r r r
"v i
mi = +e( E + v $ B)
"t
r r r
j = #N e ev e + N e ev i
Equations de Maxwell
r
r r r %E
".B = 0 "#B = µ0 ( j + $ 0 )
%t
r r
".E = 0 r %B
"#E =&
%t
• Plasma neutre !=0
• Equation du mouvement des électrons
! r r
me i"v e =!#eE
v i << v e
r
r N ee r
2 2
r N e e $E
j = #N e ev e = E =#
me i" me" 2 $t
Equations de Maxwell
r
r r %E
".B = 0 "#B = µ0$r
r %t
r
".E = 0 r %B
"#E =&
%t
2 2
N ee $p
! "r = "o # ! 2 = "o (1 # 2 )
me$ $
Propagation dans le plasma
• Propagation dans un plasma
• Masse des électrons << Masse des ions
• Fréquence plasma entre 1 et 10 Mhz
Applications
e=1.602 x 10-19 C
m=9.91 x 10-31 kg
"p=8.97 Ne1/2
Propagation entre 10 et 15 Mhz
km
degrés
Termes associés à l’ionosphère
s1
"Liono = Le # L = $ (n(s)# 1).ds ne
N = 40.3 2
s0
f
s1
= $ N(s).ds
s0
• Ne est la densité d’électrons ( en e/m3)
• !
En pratique, deux fréquences sont émises par le GPS (f1= 1227,60 et f2=1575,42!MHz)
La distance mesurée est différente sur les deux fréquences du fait de la réfraction
ionosphérique. On peut ainsi calculer le retard ionosphérique et corriger les signaux
40, 3 40,3
! " g r = + 2 # Ne ds ! "p h = # 2 $ Ne ds + ambiguïte
f f
f12 f2 2
2 (#$g r " #$g r )
1 2
TEC = ! Ne ds =
40, 3( f2 " f1 )
2
Evolution diurne et long terme du TEC
• Sur 2 mois:
– Variabilité d'un jour sur l'autre du maximum d'ionisation
– Corrélation avec l'activité magnétique
TEC sur 2 mois en Californie
Indice Kp à Boulder
Correction ionosphérique
• Correction en temps réel en utilisant les deux
fréquences ( récepteur bi-fréquence)
• Correction en temps quasi-réel en utilisant un
modèle ionosphèrique + meilleurs éphémérides
Source d’erreur Valeurs absolues (m) Remarques
Orbite 5-10 Connaissance de la position du satellite et
5 cm (IGS) autres
Horloge 10 Stabilité en fréquence du satellite par rapport au
Satellite (c.dt) temps GPS, variations d’accélérations et autres
Horloge récepteur 10 Décalage d’horloge du récepteur par rapport au
(c.dT) temps GPS
Ionosphère 2-50 Dépend de l’angle d’élévation du satellite et de
l’activité solaire
Troposphère 2-30 Dépend surtout du contenu en vapeur d’eau de
la basse atmosphère
Multitrajet (code) 0.2-3 Les signaux GPS peuvent êtres réfléchis par
(!P) des objets situés à proximité
Bruit récepteur 0.1-3 Dépend de l’électronique du récepteur et de
(code) (!P) son environnement
Multitrajet (phase) 0.001-0.03 Les signaux GPS peuvent êtres réfléchis par
(!") des objets situés à proximité
Bruit récepteur 0.0002-0.002 Dépend de l’électronique du récepteur et de
(phase) (!") son environnement
S.A (Selective 5-80 Dégradations volontaires des éléments d’orbite
Availability) et d’horloge dans les messages de navigation et
des performances des oscillateurs. S.A a été
éliminée à partir de mai 2000.
EGNOS
Termes associés à la troposphère
s1
"Ltrop = Le # L = $ (n(s)# 1).ds P 5 e
s0 N = 77.6 + 3.73.10 2
s1 T T
= 10#6 $ N(s).ds
s0
• où P est la pression de l’air (hPa), T la température de l’air
!
(K) et e la pression partielle de la vapeur d’eau (hPa).
! • Le délai #L est donc la somme de deux termes!: le
premier, appelé composante hydrostatique, représente la
contribution de l’atmosphère totale et retarde le signal
zénithal GPS de plus de 2 mètres!; le second, appelé
composante humide, représente la contribution de la
vapeur d’eau et retarde le signal zénithal GPS d’une
dizaine de centimètres.
Effet de la troposphère
• La présence de divers gaz et de la vapeur d'eau, ainsi que la variation de la
pression et de la température avec l'altitude, entraînent un écart de l'indice de
réfraction de la troposphère par rapport à l'unité.
• Cette source d'erreur est en
particulier critique pour
l'étude des petites déformations,
associées à des petits séismes
par exemple, ou à des
déformations lentes,
volcaniques ou tectoniques.
Izmit, août 1999
• Des données météorologiques
brutes, prises sur le terrain ou par
imagerie satellitaire peuvent en
revanche être utilisées pour
identifier la partie du signal SAR
cohérente avec ces données
environnementales, et donc a priori
associée à la troposphère.
• Un autre effet indirect est associé à
la topographie du terrain, en
particulier à un effet de baromètre
inversé.
Effet de cellules convectives sur la
mesure interférométrique de la phase.
Jusqu'à deux franges sont visibles.
D'après Hanssen et Feijt, 1996.
EXO
Altimétrie
Principe de l’altimétrie
• Fréquences de l’altimètre jason
– 13.5 GHz ( Bande K) et 5.3 GHz (Bande C)
– GPS (f1= 1227,60 et f2=1575,42!MHz, bande L)
– Quel est l’effet de l’ionosphère ?
Bandes de fréquence micro-ondes
Désignation Gamme de fréquences Gamme de longueur d'onde
Bande L de 1 à 2 GHz 15 à 30cm
Bande S de 2 à 4 GHz 7.5 à 15 cm
Bande C de 4 à 8 GHz 3.75 à 7.5 cm
Bande X de 8 à 12 GHz 2.5 à 3.75 cm
Bande Ku de 12 à 18 GHz
Bande K de 18 à 26 GHz
Bande Ka de 26 à 40 GHz
Bande Q de 30 à 50 GHz
Bande U de 40 à 60 GHz
Bande V de 46 à 56 GHz
Bande W de 56 à 100 GHz
Autres méthodes d’altimétrie
Altimétrie et planétologie:
• Première mesure radar depuis la
Terre (Goldstone, Aricebo, etc) 200
m de précision
• Mesure d’occultation radio des
sondes Mariner 9 et Viking
– Structure de grande échelle
s<8
• MGS: expérience MOLA ( Zuber et al.)
– Résolution latérale (taille du spot): 130-330 m
– Résolution absolue: 10 m
– Résolution relative: 35 cm
MOLA:
Mars Orbiter Laser Altimeter
Range Precision: ~ 37 cm
Absolute Vertical Resolution: <10 m
Surface Spot Size: 130 m
Along Track Shot Spacing: ~330 m
Cross-Track Spacing: ~1-30 km (equator)
Mass: 25 kg
Total number of Shots: ~640 Million
Currently acting as a radiometer
• La topographie dynamique est l'analogue du champ de pression en météorologie: les
courants océaniques de surface (flèches rouges) suivent les courbes de niveau avec une
vitesse proportionnelle à sa pente locale. Dans l'hémisphère Nord, les courants tournent
dans le sens des aiguilles d'une montre autour des bosses et en sens inverse autour des
creux, à l'opposé de ce qui se passe dans l'hémisphère Sud.La topographie dynamique
cartographié a été obtenue à partir des observations du satellite altimérique
Topex/Poséïon . L’échelle va de -110 cm (bleu foncé) à +110 cm (bleu clair). Elle permet
d'identifier la circulation océanique globale avec ses gyres et les courants de bord Ouest
(Gulf Stream dans l'Atlantique, Kuroshio au large du Japon). Notez aussi les deux
niveaux entre le Pacifique et l'Atlantique dus aux différences de salinité.
Moyens d’observations et
systèmes d’alertes
tsunami: moyens spatiaux
Les moyens d’observation
Moyens au sol:
• les sismomètres (mouvement du sol)
• les marégraphes (mouvement de la surface de l’eau)
• les capteurs de pression (mouvement de la hauteur d’eau)
• les mesures sur le terrain
Outils satellitaires:
• les satellites imageurs (photos avant/après)
• les satellites altimétriques (mesure de la hauteur d’eau de l’océan)
Outils de demain:
• télédétection spatiale et radars
Courbes de temps de parcours
Temps des ondes P et S
Distance
Marégraphes et capteurs de pression
Enregistrement de la hauteur d’eau: Capteur de pression en fond de
Mesure des marées et des tsunamis mer et système de communication
Marégraphe sur l’île de Kerguelen
$ Amélioration des modèles de prédiction de la marée
$ Localisation des sources de tsunami
$ Détermination de l’intensité des tsunamis
Les systèmes d’alerte
• Un seul système opérationnel avant décembre 2004:
Le système d’alerte du Pacifique
• 40 nations contribuent dont la France pour la Polynésie Française
• Schéma de principe du dispositif d’alerte:
Réseaux de Marégraphes et
sismomètres capteurs de pression
Localisation et Détection de
intensité du séisme l’onde de tsunami
Diffusion de
Emission d’une Validation de l’alerte l’alerte par
alerte tsunami des réseaux
locaux
La dimension spatiale
Tsunami en Mer vu par Topex et Jason
Elévation de l’océan
(m)
Satellites imageurs : Quickbird
Enregistrement d’images de l’inondation:
Carte des longueurs d’inondation
Plage de Kalutara au Sri Lanka
Juste avant l’arrivée du tsunami Juste après le tsunami: retrait la mer
$ Estimation des dégats et des longueurs d’inondation
Satellites imageurs : Quickbird
Juste avant l’arrivée de la première vague:
la mer se retire sur 400 mètres
Plage de Kalutara au Sri Lanka
400 m
Banda Aceh avant
et après le passage
du tsunami
Mécanisme des ondes sismiques atmosphériques
350 km ionosphere
150 km
A proximité:
Pulse
atmosphérique
A distance:
Séisme
Ondes de surface
Ondes de gravité
TEC….
TEC…. Et Jason Topex
