27.06.2013 Views

Cours 3: Gravimétrie

Cours 3: Gravimétrie

Cours 3: Gravimétrie

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

<strong>Gravimétrie</strong>


Pourquoi la gravimétrie?<br />

● Le but est d'étudier les variations du champ de<br />

gravité pour en déduire la répartitions des<br />

masses au sein de la Terre et ainsi sa structure.


Galilée<br />

Vers 1604 Galilée constate que quand on lâche un objet,<br />

la vitesse varie durant la chute ! Galilée propose que :<br />

vitesse = constante×temps écoulé<br />

Il en conclut que la distance parcourue est :<br />

distance = constante×0,5×temps écoulé2<br />

Son idée est confirmée dans une expérience : une<br />

gouttière inclinée le long de laquelle des clochettes sont<br />

disposées pour indiquer le passage de la bille.<br />

La constante sera notée g et sa valeur déterminée<br />

expérimentalement :<br />

g = 9,81 m∙s­2.<br />

m∙s­2


● La démarche de Galilée est­elle correcte ??


Isaac Newton (1643-1727)<br />

● 1ere loi de Newton :<br />

● 2 nd loi de Newton<br />

●<br />

● On peut la réecrire :


Définition du champ de pesanteur<br />

Il s'agit de l’accélération subit par un corps du à<br />

l'attraction de la Terre. Cela inclue la gravité et la<br />

force centrifuge.


Force centrifuge<br />

Tout corps soumis à une rotation a tendance à se déformer sous l ’effet de<br />

l’accélération centrifuge.<br />

● Vitesse de rotation de la Terre : w=7.3. 10^­5 rad.s­1<br />

ω<br />

● => Vitesse de déplacement V=wr=1700 km.s­1 à l'équateur<br />

r<br />

R<br />

φ<br />

g<br />

a c = ω 2 .r = ω 2 R.cosφ


● L'accélération centrifuge est :<br />

● Au poles : ac =0<br />

● A l'équateur : ac = 0.034 m.s­2<br />

ω<br />

r<br />

R<br />

φ<br />

g<br />

a v =a c .cosφ=ω 2 .R.cos 2 φ<br />

a c = ω 2 .r = ω 2 R.cosφ


Que vaut g sur Terre ?<br />

● g=9,8324 m.s­2 aux<br />

poles et<br />

● g=9.7803 m.s­2 à<br />

l'équateur


Définition de la verticale<br />

● L'orientation de g définit la verticale.<br />

● On simplifie le problème en supposant que<br />

l'orientation de g est la meme que ga (ie la<br />

force centrifuge ne change l'orientation de g).


Comment mesure t­on la gravité ?


Mesurer la gravité<br />

● Expérience du pendule


Principe du gravimètre<br />

● On regarde comment se déforme un ressort<br />

auquel est accroché un poids.


Gravimètre à ressort<br />

● Précision : 1/100 mGal<br />

● 1 Gal = 1 cm/s­2 = 0.01 m.s­2


Théorème de Gauss<br />

● Toute masse crée un champ de gravité :<br />

● Le théorème indique que le champ de gravité<br />

crée par une sphère est identique à celui d'une<br />

masse concentrée en son centre


●<br />

Masse de la Terre<br />

● => Quelle est la masse de la Terre (G=6.67e­11)?


Masse de la Terre


Evaluation masse d'une montagne<br />

Chimborazo (Volcan de la cordillère des Andes)<br />

Bouguet (1698­1758)


Chimborazo (Volcan de la cordillère des Andes)


Le champ de gravité est­il le meme<br />

partout ?<br />

Quelle est la forme de la Terre ?


Le champ de gravité est­il le meme<br />

partout ?<br />

Pôle Nord : 90° N 9,83245 ms -2<br />

Paris : 49° N 9,8094 ms -2<br />

Equateur : 0° N 9,7803 ms -2<br />

Java : 6° S 9,7818 ms -2<br />

Melbourne: 38° S 9,7999 ms -2


●<br />

La force centrifuge peut­elle<br />

expliquer ces différences ?<br />

● La force centrifuge n'est pas suffisante pour<br />

expliquer les variations de g.


Que se passe t­il si la distance au centre de<br />

la Terre change avec la latitude ?<br />

● G=6.67e­11 ; R=6370e3 m ; M=2.96e24 Kg;<br />

●<br />

●<br />

● => Une variation de 50 Km de la distance au<br />

centre de la Terre entraine une variation de g de<br />

~0.15 m.s­2


Pôle Nord : 90° N 9,83245 ms -2<br />

Paris : 49° N 9,8094 ms -2<br />

Equateur : 0° N 9,7803 ms -2<br />

Java : 6° S 9,7818 ms -2<br />

Melbourne : 38° S 9,7999 ms -2<br />

● Ces variations ne peuvent s'expliquer<br />

uniquement par la force centrifuge.<br />

● La Terre n'est pas sphérique !<br />

● => La Terre à une forme d'ellipsoide :<br />

● R_pole = 6357 Km R_équateur=6378 Km


Variation de la pesanteur à la<br />

surface du globe<br />

● Le champ de pesanteur vaut en moyenne 981<br />

gals.<br />

● La pesanteur dépend de<br />

● l'altitude.<br />

● La latitude car la Terre n'est pas une sphère<br />

parfaite<br />

● De la position, car les masses ne sont pas réparties<br />

uniformément au sein de la Terre.


Au premier ordre la forme de la Terre peut<br />

Surface ellipsoïdale<br />

de référence<br />

s'approximer par un ellipsoide<br />

Axe de rotation<br />

● si la Terre était homogène et immobile sa surface serait<br />

une sphère parfaite et une équipotentiel de gravité (g=<br />

constant)<br />

● Rotation Terre => Terre est un ellipsoide applati au poles à<br />

cause de la force de centrifuge


En regardant dans le détail ce n'est<br />

Surface ellipsoïdale<br />

de référence<br />

pas tout à fait vrai<br />

Axe de rotation<br />

Le Géoïde<br />

● La répartition des masses au sein de la Terre<br />

est complexe => ce n'est pas exactement un<br />

ellipsoide


Ellipsoide et géoide<br />

● Ellipsoide : si la Terre était homogène et<br />

immobile sa surface serait une sphère parfaite<br />

et une équipotentiel de gravité (g= constant)<br />

● Rotation Terre => Terre est un ellipsoide applati<br />

au poles.


Géoide<br />

● Géoide = surface équipotentielle passant par le<br />

niveau moyen des océans.<br />

● En se promenant sur le Géoide, on garde la<br />

meme énergie potentiel : on ne monte pas, on<br />

ne descend pas.<br />

● Le géoide est en tout point perpendiculaire au<br />

vecteur de gravité.<br />

● L'altitude précise = distance au géoide.


Comment déterminer la forme du<br />

géoide ?


Comment mesurer le Géoide ?<br />

● Les mesures par satellites ne permettent de<br />

voir le géoide dans les grandes longueurs<br />

d'onde.


Mesure par satellite


Échelle : violet = -95m ; rouge = +50m


Comment interpréter le géoide ?<br />

● Hauteur géoide = distance par rapport à<br />

l'ellipsoide.<br />

● L'energie potientielle de gravité est la meme<br />

partout sur le géoide :<br />

●<br />

● Si la masse augmente il faut monter pour avoir<br />

la meme energie potentiell => bosse dans le<br />

géoide.<br />

● Si on enlève de la masse il faut descendre pour<br />

avoir la meme energie potentielle.


● Les variations du géoide peuvent­elle<br />

s'expliquer par la structure de la croute terrestre<br />

? Ou par la topographie ?


Tomographie sismique vs géoide à grande<br />

longueur d'onde<br />

Vitesse des ondes sismiques<br />

Rouge=lent ; bleu = rapide<br />

Géoide :<br />

Rouge =bosse ; bleu = creux<br />

Les ondulation du Géoide sont liées à la géodynamique interne


De quoi dépend le géoide ?<br />

● 1)Densité des matériaux (chimie, température),<br />

● 2) topographie<br />

● 3) l'épaisseur de la croute


Bilan : que voit ­on grace au


● Les ondulations de petite échelle réflètent la<br />

topographie<br />

● Exemple d'un mont sous marin :


Conclusion : Plus on regarde les fluctuation du géoide à<br />

grande échelle, plus on voit profondément dans la Terre<br />

● Les variations de petites échelles dépendent de<br />

la topographie<br />

● A l'échelle du millier de Km, le géoide dépend<br />

du volcanisme, de la convection dans le<br />

manteau superficielle.<br />

● A grande échelle, les mesures pas satellites, le<br />

géoide dépend de la Température et de la<br />

dynamique interne.


On peut étudier la répartition des masses<br />

au sein de la Terre de deux façons<br />

● 1) Géoide : on étudie la forme de la Terre.<br />

Efficace en milieu marin. Donne une vision<br />

globale de la Terre gràce aux satellites. Mais on<br />

ne voit pas les détails.<br />

● 2) On mesure partout g à la surface de la Terre.<br />

=> On peut voir plus de détail sur les<br />

continents.


Aller plus loin dans l'interprétation<br />

du champ de gravité


Mesure par bateau<br />

● On distingue les variations du géoide à petite<br />

échelle. Elles sont dues à la topographie


A l'échelle du millier de km<br />

● Anomalie de 1000 à 2000 km // aux chaines de volcans<br />

● => Anomalie positive correspondent aux volcans<br />

● => Relation avec la convection du manteau?


Les anomalies de gravité<br />

● La gravité dépend de divers facteurs : Altitude,<br />

latitude, topographie en surface, structures<br />

profondes<br />

● => On souhaiterait isoler ces différents facteurs.<br />

● => On introduit différentes corrections


Correction d'altitude (Faye) =<br />

correction à l'air libre<br />

● Quantifie la variation de g due à l'altitude ie la<br />

distance du point de mesure du géoide.<br />

● Sur l'ellipsoide on a :<br />

● A une altitude h :<br />

● Avec un dévelopement limité:<br />

M<br />

h<br />

Référence<br />

= géoïde


Correction plateau<br />

● On quantifie la variation de g due à la matière<br />

se trouvant entre l'ellipsoide et le point de<br />

mesure.<br />

● On suppose une densité de 2.67<br />

● La matière entre l'ellipsoide contribue à la<br />

gravité de la façon suivante :<br />

●<br />

M<br />

p<br />

d<br />

Référence<br />

= géoïde<br />

Cette correction tient compte de la densité,d,<br />

du matériel présent entre la surface de<br />

référence et le point de mesure..


Correction topographique<br />

● Plutot que de supposer qu'entre le point de<br />

mesure et l'ellipsoide il y'a de la matière partout,<br />

on tient compte de la topographie réelle.<br />

● Cette correction peut etre importante dans les<br />

zones à relief contrasté (montagnes).<br />

Loin de la montagne<br />

g o<br />

g o<br />

∆g<br />

g mont<br />

g m


● Il s'agit de la somme<br />

des corrections<br />

d'altitude, de plateau<br />

et de topographie.<br />

Correction Bouger


Comment interpréter les différentes<br />

corrections<br />

● Anomalie à l'air libre : dépend de la masse sous<br />

le point de mesure, pas de son altitude.<br />

● Anomalie Bouguer : information sur ce qu'il y'a<br />

en profondeur.


Comment interpréter les différentes<br />

corrections<br />

● Anomalie à l'air libre : information sur le relief<br />

(faible)<br />

● Anomalie Bouguer : information sur ce qu'il y'a<br />

en profondeur sous le niveau de la mer.<br />

Valeurs plus élevées.<br />

● => On peut pas savoir le volume, la densité, la<br />

profondeur de la masse perturbante à partir<br />

uniquement de l'anomalie de Bouguer.


Anomalie sous une dorsale océanique<br />

● Anomalie à l'air libre faible : 80 mgal<br />

● Elle suit la topographie


● Que doit valoir l'anomalie à l'air libre et<br />

l'anomalie de bouguer dans ces deux cas?


● Si le relief est compensé en profondeur (équilibre<br />

isostasique) l'anomalie air libre est ~ 0 au milieu


● Quelle est l'anomalie de Bouguer dans ce cas<br />

là ?


Anomalie d'une fausse océanique<br />

● Anomalie air libre fortement négative au niveau de la fosse et du mur<br />

interne (­280 mgal) et positive au niveau de l'arc interne (+80 mgal)<br />

● Pas d'équilibre isostasique


Anomalie de Bouguer au dessus<br />

d'un rift ?


Anomalie au dessus d'un rift


D'autres exemples<br />

Topographie Anomalie Air libre


Mars Gravity<br />

● Tharsis<br />

● Large free­air anomaly indicates it is<br />

uncompensated<br />

● But it’s too big and old to last like this<br />

● Mantle plume provides dynamic<br />

compensation?<br />

● Utopia<br />

● Probably compensated initially<br />

● Filled later when lithosphere was thicker<br />

● Crustal thickness<br />

● Assume Bouguer anomalies caused by<br />

thickness variations in a constant density<br />

crust<br />

● Need to choose a mean crustal thickness<br />

● Isidis basin sets a lower limit<br />

Free Air<br />

Zuber et al., 2000


Anomalie de Bouguer en France


Interprétation<br />

● Les chaines de montagnes sont liées à des<br />

anomalies négatives<br />

● => Il existe un défaut de masse en profondeur


Principe de l'isostasie<br />

● La masse en surface est compensée en<br />

profondeur.


Isostasie : principe<br />

● A l'équilibre les forces de pression et de gravité<br />

s'équilibre.<br />

● Si on rajoute ou soustrait une masse cet<br />

équilibre est rompu => ceci entrain des<br />

mouvements jusqu'à ce que l'équilibre est<br />

atteint.<br />

● Il y'a équilibre lorsque la pression est la meme<br />

partout à une profondeur donnée.


Deux visions de l'isostasie :<br />

1 Pratt (1855)<br />

● Chaines de montagne résulte de la dilatation thermique de<br />

la croute dues à diverses sources de chaleur<br />

● => Montagnes = dilatation de la croute terrestre = baisse<br />

densité<br />

ρ m<br />

ρ c3<br />

1<br />

ρ c1<br />

ρ c2 ρ c2 ρ c3<br />

Surface de compensation<br />

La surface de compensation : surface sur laquelle la pression exercée<br />

par chaque colonne 1,2, 3 est égale<br />

2<br />

3


Deux visions de l'isostasie : Airy<br />

● Les reliefs sont compensés par une racine<br />

crustale (et inversement les dépression par une<br />

„anti racine”.)<br />

ρ m<br />

ρ c<br />

1 2<br />

Surface de compensation<br />

La surface de compensation : surface sur laquelle la pression exercée par chaque colonne 1,2, 3 est égale<br />

3


Anomalies isostasiques<br />

● On peut calculer l'altitude d'un structure (chaine<br />

montagne) si l'équilibre isostasique était atteint.<br />

● On comparant cette altitude théorique et<br />

l'altitude réelle nous dit si la structure est en<br />

équilibre ou non<br />

● => On peut ainsi prédire si les structures<br />

géologiques vont monter ou descendre.


● Tectonique des plaques => évolution<br />

perpetuelle du relief<br />

● => La réalisation de l'équilibre isostasique<br />

prend du temps<br />

● => On peut avoir un désiquilibre temporaire<br />

isostasique<br />

● => on alors une anomalie gravimétrique air libre


Réajustement isostasique : scandinavie et<br />

antarctique


Rebond post­glacière


Isostasie et érosion des chaines de<br />

montagnes


Isostasie régionale : flexure de<br />

plaque


Subsidence<br />

● Explique l'épaisseur de certaines formations<br />

sédimentaires


Conclusion<br />

● <strong>Gravimétrie</strong> permet de :<br />

● 1) Estimer la densité de matériaux<br />

● 2) Comprendre la géodynamique globale<br />

(grande longueur d'onde)<br />

● 3) Etude de la bathymétrie (petite longueur<br />

d'onde).<br />

● 4) Comprendre les mouvements verticaux de la<br />

lithosphère


● Gravité: definition …<br />

● De quoi dépend la gravité que peut­on voir et<br />

ne pas voir avec ?<br />

● Mesure du champ de gravité<br />

● Interprétation des mesures (anomalies bouguet)<br />

● Isostasie<br />

● Evolution du relief...


Anomalie air libre sur la lune


Anomalie air libre sur la lune<br />

● Anomalie négative sur le pourtoure des cratères : de la<br />

masse a été enlevé<br />

● Au centre : anomalie positive : sur­compensation<br />

lithospherique


Mars Gravity<br />

● Tharsis<br />

● Large free­air anomaly indicates it is<br />

uncompensated<br />

● But it’s too big and old to last like this<br />

● Mantle plume provides dynamic<br />

compensation?<br />

● Utopia<br />

● Probably compensated initially<br />

● Filled later when lithosphere was thicker<br />

● Crustal thickness<br />

● Assume Bouguer anomalies caused by<br />

thickness variations in a constant density<br />

crust<br />

● Need to choose a mean crustal thickness<br />

● Isidis basin sets a lower limit<br />

Free Air<br />

Zuber et al., 2000


Résumé<br />

● Before we can start interpreting gravity gravity anomalies we need to make sure we’re comparing apples to apples…<br />

● Free­Air correction<br />

■ Assume there’s nothing but vacuum<br />

between observer and reference ellipsoid<br />

■ Just a distance correction<br />

Δg<br />

Δg<br />

FA<br />

FA<br />

δg<br />

δ ⎛ GM ⎞<br />

= Δr<br />

= ⎜ h 2 ⎟<br />

δr<br />

δr<br />

⎝ r ⎠<br />

2gh<br />

= −<br />

r<br />

● Bouguer correction<br />

■ Assume there’s a constant density plate between observer and reference ellipsoid<br />

■<br />

■<br />

Remove the gravitation attraction due to the mass of the plate<br />

If you do a Bouguer correction you must follow up with a free­air correction Δg<br />

B = 2π<br />

Gρh

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!