Ces chapitres présentent la structure des planètes et autres objets des systèmes planétaire,
Plusieurs processus sont à l'oeuvre à la surface des planètes et petits corps du Système Solaire. Leur nature, intensité et conséquence varient principalement en fonction de la composition de la surface (silicatée/glacée), de sa distance à l’étoile et la présence ou non d’une atmosphère. Si le corps planétaire est actif, il peut régulièrement faire « peau neuve » (par volcanisme notamment) alors qu’un corps mort est condamné à voir sa surface se dégrader, parfois lentement mais toujours sûrement.
Dans ce qui suit, nous décrirons les principaux processus qui façonnent les surfaces planétaires, à savoir :
Les cratères d’impact sont des dépressions de surface, généralement circulaires, résultant de l’impact de fragments solides d’origine météoritique ou cométaire. Ils sont présents sur quasiment toutes les surfaces planétaires solides du système solaire et la cratérisation représente la principale cause d’altération des surfaces des corps dépourvus d’atmosphère (excepté pour Io et Europe). Les impacts sont aussi le principal mécanisme d’apport et d’excavation de matériel planétaire ; certaines planètes ou satellites se sont formés par agrégation après collision (la Lune, par exemple) ; les volatiles, océans et atmosphères, présents sur certains corps ont sans doute été apportés par des objects impacteurs (en l’occurrence des comètes).
L’étude des cratères est intéressante à plus d’un titre. Leur forme et leur taille renseignent sur la nature (composition, résistance, stratifications, porosité) des surfaces planétaires cibles, de leurs impacteurs et les propriétés d’une éventuelle atmosphère. Leur distribution informe sur l’âge des surfaces planétaires et apporte des clés dans la compréhension de leur histoire d’autant que les impacts peuvent être à l’origine d’évènements importants voire catastrophiques (formation de la Lune, extinction des dinosaures sur Terre il y a 65 millions d’années…), les plus larges ayant même pu modifier les paramètres orbitaux de certains corps. Les études statistiques des cratères fournissent aussi des informations sur la population des corps impacteurs du système solaire qui sont ce qui reste des planétésimaux de l’accrétion planétaire et donc, à ce titre, des objets très primitifs, témoins privilégiés de la jeunesse du système solaire.
Ce bestiaire de cratères illustre la variété en termes de forme et d’échelle des cratères planétaires. La morphologie des cratères dépend avant tout de leur taille. Les plus petits d’entre eux, les cratères simples, présentent une forme en bol, avec des bords surélevés. La majorité des cratères lunaires ayant un diamètre inférieur à 15 km sont de ce type. Au delà de ce diamètre, ils présentent un pic central avec éventuellement des terrasses et des dépôts et sont dits complexes. Les cratères complexes sont moins profonds que les cratères simples. Le diamètre de transition entre ces deux types de cratères varie de façon inversement proportionnelle à la gravité de la planète : un pic central apparaît dans les cratères martiens dont la taille est supérieure à 10 km alors que sur Terre, les pics apparaissent dès que les cratères ont un diamètre supérieur à 2-3 km. Dans les cas d’impacts plus gros, le pic central est remplacé par un anneau montagneux voire par des anneaux multiples si le diamètre augmente encore. Une partie des matériaux excavés par l’impact peut, après avoir parcouru une trajectoire balistique, retomber et provoquer la création de cratères secondaires, à proximité du cratère principal. Enfin, si l’impact est suffisamment puissant pour percer la croûte et provoquer des épanchements, on parle de bassins d'impact. Le bassin d'Hellas sur Mars (plus de 2000 km de diamètre) est sans doute le plus grand bassin d'impact du système solaire. Avec le temps, l’érosion ou encore les mouvements du sol altèrent les cratères dont nous n’observons finalement qu’une forme dégradée.
Les corps qui ne sont pas protégés par une atmosphère (Mercure, Lune, astéroïdes, satellites glacés du système extérieur…) sont directement et constamment exposés aux rayons d’origine solaire et galactique ainsi qu’au bombardement micro-météoritique. Ils subissent, de plus, de grandes variations de température. Ces différents phénomènes dégradent progressivement leur surface, généralement caractérisée par une couche superficielle très poreuse appelée régolithe.
La notion de régolithe est large ; elle désigne « tout matériel d’origine continentale, quel qu’en soit l’âge, recouvrant les roches saines et dures ». Cependant, en planétologie, on l’utilise généralement pour évoquer la couche de poussière, rocheuse ou glacée, recouvrant la surface de corps solides dépourvu d’atmosphère protectrice (Mercure, la Lune, les satellites de Mars, tous les satellites du système extérieur sauf Titan, les astéroïdes et, dans une moindre mesure, Mars et les comètes).
Les régolithes planétaires sont d’abord créés par impacts météoritiques (voir chapitre Cratérisation des surfaces) puis évoluent sous l’effet de l’érosion spatiale (ou « space weathering » en anglais) c’est-à-dire des effets combinés du bombardement micro-météoritique, de la collision des rayons cosmiques d’origine solaire ou galactique ou encore de l’irradiation et de la pulvérisation cathodique (« sputtering ») par les particules du vent solaire. Les chocs thermiques, auxquels sont particulièrement soumis les surfaces sans atmosphère, contribuent également à la désagrégation physique (ou thermoclastie) des roches et donc au développement du régolithe.
Les propriétés et le degré de maturité d’un régolithe varient en fonction de la composition et de la position dans le Système Solaire de l’objet planétaire. Le régolithe lunaire est de loin le mieux connu. Il recouvre l'ensemble de la surface du satellite sur une profondeur de 2 à 8 mètres dans les mers et pouvant même excéder 15 m dans les terres les plus anciennes (4.4 Ga). Il repose sur plusieurs mètres d’un méga-régolithe constitué de gros blocs rocheux, éjectas d’anciens grands impacts. Le régolithe de Mercure est probablement très semblable à celui de la Lune quoique peut être légèrement plus développé car le flux micro-météoritique y est plus important et le contraste thermique entre le jour et la nuit accru. Le développement d'un régolite mature est, en revanche, nettement plus lent sur les astéroïdes en raison de leur faible gravité. Sur Mars, qui possède une atmosphère tenue, l’érosion spatiale s’est combinée à d’autres formes d’érosion (hydrique, éolienne… voir chapitre Erosion et sédimendation des surfaces avec atmosphère) pour former un épais manteau de poussière et de débris. Io est aussi un cas à part car le volcanisme qui y sévit efface immédiatement les traces d’impact. Enfin, la volatilité de la glace d’eau, ainsi que des glaces de CO2 (dioxyde de carbone) ou de CH4 (méthane), rend les surfaces glacées du système solaire particulièrement vulnérables à l’érosion spatiale.
Les surfaces des corps dotés d’une atmosphère (Vénus, Terre, Mars, Titan, éventuellement Pluton) sont protégées de l’érosion spatiale et subissent modérément le bombardement météoritique et les effets thermiques. En conséquence, elles n’ont pas ou peu de régolithe. En revanche, elles sont soumises à l’action conjuguée de l’air et de solvants liquides lorsqu’ils existent. Partout où il y a une atmosphère, même extrêmement dense (Vénus) ou, au contraire, tenue (Mars), l'activité éolienne transporte les sédiments les organisant notamment en champs de dunes. Sur Pluton de telles formations n'ont pas encore été observées sans doute pour des raisons de résolution. Sur Terre, aujourd’hui ou sur Mars, hier, le cycle de l’eau a, en outre, façonné la surface via l’érosion pluviale, fluviatile ou glaciaire créant, transformant et distribuant la matière sédimentaire. Sur Titan, c’est le méthane et l’éthane qui modifient les paysages.
Les corps solides du système solaire présentent des visages multiples et en même temps étrangement familiers. C’est que des processus semblables à ceux que l’on observe sur Terre y sont à l’œuvre. La planétologie comparée consiste à s’appuyer sur la connaissance de notre planète pour comprendre comment ont évolué d’autres mondes. Au passage, nous observons des processus a priori connus se développer dans des environnements radicalement différents et du même coup enrichissant la compréhension que nous en avons.
Ce principe s’applique particulièrement bien lorsqu’on compare les paysages terrestres à ceux des trois autres corps du système solaire possédant une atmosphère : Vénus, Mars et Titan. Tous présentent à leur surface des dunes, preuves de l’activité éolienne qui y sévit (voir figure ci-contre). Cependant la taille, la forme de ces dunes varient d’une planète à l’autre. La connaissance des dunes terrestres permet de comprendre l’origine de ces divergences et donne des clefs pour déduire de l’observation des dunes extraterrestres les régimes de vents qui les ont sculptées. On reconnait également la signature de l’érosion fluviale à la surface de Mars, pourtant sèche aujourd’hui. L’étude des paléo-réseaux fluviaux ou des vallées de débâcle permet de reconstituer une partie de l’histoire géologique de la planète rouge. Sur Titan, par analogie avec le cycle de l’eau sur Terre, les phénomènes météorologiques liés au cycle du méthane et de l’éthane peuvent aussi être mieux compris et prédits. Enfin, même Pluton, récemment observé par la sonde américaine New Horizons, présente des paysages connus, en particulier des sols polygonaux typiques des régions glaciaires et périglaciaires des hautes latitudes sur la Terre et sur Mars (voir figure ci-contre).
Sur les corps les plus actifs tels que la Terre, Io ou Europe, le renouvellement de la surface s’opère principalement par l’activité volcanique et tectonique. Il s’agit de processus endogènes c’est-à-dire ayant une cause interne, par opposition aux processus exogènes décrits précédemment.
Le volcanisme est lié au transfert de matière (magma, éléments volatiles et matériaux cristallisés) de l’intérieur vers la surface. Il est l’une des expressions les plus spectaculaires de l’activité interne d’un corps. Il participe au renouvellement des surfaces en recouvrant de ses épanchements (laves) les cicatrices du passé et peut aussi contribuer à la création ou à l’enrichissement d’une atmosphère. De nombreuses planètes ou satellites du système solaire portent sur leur surface les traces d’une activité volcanique passée (Mars, Lune), récente (Vénus) et même présente (Terre, Io, Europe, Encelade).
L’activité tectonique est l’ensemble des mécanismes de mouvements de surface responsables de déformations à grande échelle de la croûte d’une planète. Beaucoup de corps présentent à leur surface des failles, des chaînes dorsales et escarpements témoignant d’une activité tectonique passée. Cependant, à ce jour, le mouvement de plaques tectoniques n’a été observé que sur la Terre.
Le processus d’impact et ses conséquences varient avec la vitesse du corps impacteur et la nature du sol impacté et de son impacteur. Si la planète possède une atmosphère, le projectile est freiné et chauffé ce qui peut entrainer sa vaporisation partielle voire totale ou sa fragmentation. Les météorites de moins de 10 m de diamètre parviennent rarement jusqu’au sol terrestre. Les modèles d’ablation atmosphérique prédisent un nombre réduit de cratère de moins de 20 km de diamètre sur Titan (seul satellite du système solaire possédant une atmosphère substantielle, voir chapitre Erosion et sédimendation des surfaces avec atmosphère), ce que semble confirmer les observations de la sonde Cassini.
Lorsque le projectile, ou ce qu’il en reste, atteint la surface, le processus d’impact commence ; on le décompose classiquement en trois phases qui, en réalité, se chevauchent dans le temps : la phase de contact et compression, la phase d'excavation et la phase de modification et relaxation. C'est ce qu'illustre la figure ci-contre.
L'appliquette "Cratérisation" permet d'appréhender les effets d’un impact météoritique sur Mercure, la Terre (avec ou sans atmosphère), la Lune et Mars en fonction des caractéristiques de l’impacteur (vitesse, angle d’arrivée, taille, densité) et de la surface impactée (densité).
En l’absence d’échantillon du sol, la datation par comptage de cratères est la seule méthode pour estimer l’âge relatif des surfaces planétaires. Celle-ci s’appuie sur deux règles simples:
Ces règles reposent sur l’idée que la population des impacteurs a évolué au cours du temps ; la taille des projectiles et le taux de cratérisation étaient nettement plus importants dans la jeunesse du système solaire, à une époque où les débris étaient abondants. Ces derniers ont progressivement été mobilisés pour former les planètes, les plus gros planétésimaux disparaissant en premier jusqu’à ce que les impacteurs moyens puis petits se fassent rares aussi. Si le bombardement météoritique a affecté de façon uniforme la surface d’une planète donnée, certaines régions en ont gardé toutes les cicatrices alors que d’autres ont connu depuis des épisodes de rajeunissement (par volcanisme par exemple).
L’échelle ci-contre renseigne sur le niveau de cratérisation des principales surfaces solides du système solaire, les surfaces les plus jeunes étant les moins cratérisées.
Ainsi, l’étude de la distribution des cratères (nombre de cratères en fonction de leur taille) permet-elle de donner un âge relatif à différentes unités de surface. Pour déterminer leur âge absolu il faudrait connaître précisément l’histoire de l’évolution du flux d’impacteurs dans le système solaire. Une partie de cette histoire a pu être retracée grâce à la datation radiogéniques d’échantillons lunaires collectés lors des missions Apollo. Ces datations précises, comparées à la distribution des cratères lunaires, ont permis de dresser des courbes d’évolution dans le temps de la densité et de la taille des impacteurs, révélant notamment un pic d’impacts il y a environ 4 milliard d’années lors d’une phase appelée le Grand Bombardement Tardif (ou Late Heavy Bombardement, LHB). Les surfaces du système solaire ayant atteint le niveau de saturation sont sans doute vieilles de 4 milliards d’années.
En tenant compte du fait que le flux des impacteurs devait varier avec la distance au Soleil et, lorsque cela est nécessaire, de la présence d’une atmosphère, les enseignements du cas lunaire peuvent être extrapolés afin de dater les autres surfaces planétaires. Cependant il est important de garder à l’esprit que cette extrapolation est sujette à caution ; la position des planètes et la densité des atmosphères ont pu, en effet, varier au cours de l’histoire du système solaire.
Les surfaces des corps sans atmosphère sont de véritables champs de bataille subissant en permanence:
Tous ces processus sont exogènes, c’est-à-dire ayant une cause externe à l’objet qu’ils affectent. Lentement mais sûrement ils érodent les surfaces mais sont aussi à l’origine de la formation des atmosphères extrêmement tenues, appelées exosphères, que l’on trouve autour de Mercure, de la Lune, de la plupart des satellites glacés et même des anneaux de Saturne.
Les surfaces sans atmosphère sont constamment bombardées par des grains météoritique de diamètre <1 mm. Ce bombardement micro-météoritique ne contribue pas à augmenter le volume de régolithe (le régolithe lunaire ne grandit que de 1 mm/106 ans et, en raison de leur faible gravité, les astéroïdes perdent même constamment une partie de leur régolithe) mais il en modifie durablement les propriétés et la distribution.
Les impacts micro-météoritiques pulvérisent progressivement les premiers millimètres du sol, réduisant la taille des particules à la surface. Ce phénomène de fragmentation (ou «comminution» en anglais) est, en partie, compensé par un phénomène d’agglutination: lorsque les impacts sont suffisamment rapides, certains matériaux du sol fondent et, en refroidissant, se soudent (formant des sphérules de verre sur la Lune par exemple) ou soudent entre eux des fragments de roches et de minéraux donnant naissance à des particules plus grosses. Certains matériaux sont même vaporisés sous l’effet des micro-impacts avant d’être redéposés à la surface. Le régolithe lunaire est constitué d’environ 30% d'agglutinates, agrégats dont la taille varie de quelques micromètres à quelques millimètres et présentant à leur surface des nanoparticules de fer intégrées lors de la vaporisation puis re-condensation de minéraux ferrifères (olivine et pyroxène notamment). L’érosion spatiale sur la Lune est donc synonyme d’un obscurcissement (l’albédo diminue) et d’un rougissement de la surface avec le temps. Le régolithe lunaire est dit mature lorsque les processus de fragmentation et d’agglutination se compensent ; la taille des grains est alors ~60 μm. Un régolithe immature est constitué de grains plus gros et d’une proportion réduite d’agglutinates.
Dans le même ordre d'idée, sur les surfaces glacées, le bombardement micro-météoritique participe à la recristallisation de la glace lorsqu’elle est amorphe (c’est-à-dire sans arrangement précis, par opposition à la glace cristalline qui présente une structure héxagonale) à la surface par un processus de recuit (« annealing » en anglais) et lutte donc contre le travail d’amorphisation mené par les rayons solaires UV et les particules ionisées énergétiques (voir Radiations d’origine solaire et cosmique).
Enfin, sur les astéroïdes où la vitesse d’échappement est faible, le bombardement micro-météoritique, aidé par d’autres processus tels que le « sputtering» (voir Radiations d’origine solaire et cosmique), contribue à l’éjection et à la perte des particules les plus petites. Ainsi s’attend-on à trouver un sol plus grossier à la surface des plus petits astéroïdes.
Le bombardement micro-météoritique modifie également la distribution des composés des régolithes. Les premiers millimètres du sol lunaire sont en permanence « labourés » par des micro-impacts ce qui a pour effet d’homogénéiser la composition verticale (en profondeur) du régolite. On parle d’"impact gardening" (de l’anglais « garden », jardiner). Ce processus est néanmoins très lent – il faut au moins 100 000 ans pour entièrement retourner et mélanger le premier centimètre du sol lunaire. Les couches plus profondes du régolithe ne sont retournées qu’à l’occasion d’impacts plus importants donc plus rarement.
La distribution horizontale des composés du régolite est, quant à elle, contrôlée par les lois de retombée balistique des éjectas autour du cratère principal (voir Processus de cratérisation des surfaces) et varie peu sous l’effet du bombardement micro-météoritique. Les micro-impacts peuvent néanmoins, localement, apporter de nouveaux éléments à la composition de surface.
Les surfaces sans atmosphère sont également soumises à un bombardement permanent par des particules plus ou moins énergétiques en premier lieu desquelles des photons X et ultra-violet (UV) solaires, des ions issus du vent solaire et des rayons cosmiques provenant de notre Galaxie ou d’au-delà. Ces radiations modifient chimiquement, physiquement et structurellement les surfaces sur une profondeur allant de quelques micromètres à quelques mètres, en fonction de l’énergie des particules.
Le vent solaire est un flux de plasma essentiellement composé de particules d’hydrogène et d’hélium ionisées dont l’énergie est modérée (0.3-3 keV/nucléon). Ce flux varie, en température et en vitesse, avec l’activité du Soleil. Lors d’éruptions solaires et d’éjection de masse coronale, des rafales de particules solaires particulièrement énergétiques (1-100 MeV/nucléon) balayent notre système stellaire.
Les corps pourvus d’un champ magnétique propre (Mercure, Terre, Ganymède) sont protégés en grande partie des radiations, leur magnétosphère déviant les particules chargées le long des lignes de champ et agissant ainsi comme un bouclier. A l’inverse, les magnétosphères des géantes gazeuses, en piégeant et accélérant les particules chargées, produisent d'intenses ceintures de rayonnement et soumettent les satellites qui leur sont les plus proches à de grandes doses de radiations. En particulier, Io et Europe, autour de Jupiter, reçoivent des doses 100 à 1000 fois plus élevées que la Lune.
Les principaux effets du bombardement par les particules solaires et cosmiques sur les surfaces sont les suivants :
Sur ce dernier point, notons que les surfaces glacées sont particulièrement sensibles aux radiations car elles sont trois fois moins résistantes que les surfaces silicatées et plus volatiles (c’est-à-dire susceptibles de changer de phase). Rappelons que la glace d’eau peut exister sous plusieurs formes: différents états cristallins (en fonction essentiellement de la température) ou amorphes. A basses températures, le bombardement par les particules UV et les ions peut modifier la structure cristalline de la glace en surface, voire même entrainer son amorphisation. Europe, qui baigne dans magnétosphère jovienne et est, par conséquent, soumise à des taux de radiation particulièrement élevés, présente une surface jeune mais largement amorphisée alors que la phase cristalline est stable à la surface de Callisto, satellite près de 3 fois plus éloignée de Jupiter. Ganymède, qui se trouve entre Europe et Callisto et est, de surcroît, protégé par un champ magnétique propre, présente de la glace amorphe aux pôles (là où les lignes de champs sont ouvertes) et cristalline ailleurs.
Les corps sans atmosphère peuvent connaitre des variations de températures considérables au cours d’une journée. Plus l’inertie de la surface (c’est-à-dire sa capacité à stocker la chaleur) est faible et plus la rotation du corps est lente, plus le contraste jour/huit est important. Mercure, en particulier, subit les chocs thermiques les plus violents du système solaire : la température à sa surface peut varier de -170°C à 430°C.
Sur les surfaces silicatés du Système Solaire, la différence de réponse (dilatation/contraction) des minéraux des roches à l’alternance jour/nuit induit des contraintes mécaniques (surpression) pouvant déboucher sur la fissuration progressive voire l’éclatement de certaines roches. Plus les changements de température sont prononcés et rapides, plus ce processus de désagrégation, appelé thermoclastie, est efficace. En outre, les roches dont la taille excède la profondeur de peau diurne (la profondeur du sol qui subit les fluctuations diurnes du flux solaire - en général quelques centimètres) sont soumises à un fort gradient de température qui peut les fragiliser à long terme.
Sur les surfaces glacées du système solaire, en raison de la grande volatilité de la glace d’eau (c’est-à-dire de sa capacité à changer de phase), le cycle jour/nuit peut s’accompagner d’un phénomène de migration/ségrégation thermique.
Les surfaces des satellites glacés sont généralement constituées d’un mélange, aux proportions variables, de glace et d’un composé optiquement sombre (matière organique, soufrée ou silicatée). Les régions les plus riches en glace étant aussi les plus brillantes, elles sont moins efficacement chauffées par le Soleil (elles réfléchissent une grande partie du flux solaire) et le taux de sublimation de la glace y est bas. Inversement, dans les régions les plus sombres, le taux de sublimation de la glace peut être élevé. Ce déséquilibre permet un transfert de la glace des régions sombres et chaudes vers régions brillantes et froides. Ce transfert prend fin quand toute la glace du sol des régions sombres a disparu (laissant un sol encore plus sombre) et s’est redéposé dans les régions plus claires (accentuant alors leur brillance). Au passage, il est fréquent qu’une partie des volatiles se perde dans l’espace ou vienne enrichir une exosphère.
Le phénomène de migration/ségrégation thermique a pour effet de renforcer les contrastes d’albédo à la surface et, en séparant la glace de la matière sombre, va à l’encontre des processus de bombardements (météoritiques ou par des particules de haute énergie) qui tendent à homogénéiser le régolithe. Ce phénomène peut être local (exemple de Callisto) ou global (exemple de Japet). Sur Japet, même si l’origine de la matière sombre est vraisemblablement exogénique (en provenance de Phoebe), il est fort probable que le phénomène de migration/ségrégation participe à accentuer le contraste entre les basses et moyennes latitudes, très sombres, de la face avant du satellite et les pôles, particulièrement brillants : le jour, la glace des régions équatoriales se sublime et migre vers les pôles, plus froids, ou elle se re-condense. Sur Callisto, les crêtes des cratères des régions équatoriales sont recouvertes d’un manteau blanc résultant sans doute de la migration de la glace du fond des dépressions, généralement plus chaudes car doublement chauffées, à la fois par le flux solaire direct et par le flux solaire réfléchi sur les parois. Ce processus de modification du paysage par sublimation est aussi à l’œuvre sur Mars où ni l’eau ni le CO2 ne sont stables à la surface.
L’appliquette "Migration" vise à évaluer l’efficacité du phénomène de migration/ségrégation par rapport à d’autres processus d’érosion spatiale sur les principaux satellites glacés du Système Solaire.
L'appliquette Migration
Les surfaces des corps avec atmosphère sont certes protégées de l’érosion spatiale mais elles subissent d’autres formes d’érosion, souvent plus efficaces, liées à l’activité atmosphérique. Ces formes d’érosion se traduisent généralement par une perte graduelle de substance et notamment de relief. Elles sont aussi parfois à l’origine de paysages spectaculaires.
Il existe plusieurs agents et types d’érosion. Dans ce qui suit nous aborderons :
Et nous mentionnons seulement ici:
Quelque soit son moteur, l’érosion comporte trois phases étroitement liées :
Les gaz (par le biais du vent), l’eau, le méthane, l’éthane liquides sont des fluides. Avant de décrire leur action sur les surfaces planétaires, quelques rappels sur l’écoulement des fluides s’imposent.
Un fluide est un milieu matériel déformable (il change de forme sous l'action d'une force extérieure). Il est constitué de particules libres de se mouvoir les unes par rapport aux autres. Les liquides sont des fluides peu compressibles ; ils conservent le même volume quelque soit leur forme et présentent une surface propre. Les gaz, en revanche, sont des fluides compressibles ; ils tendent à occuper tout l'espace disponible et n'ont pas de surface propre.
Qu’ils soient gazeux ou liquides, les fluides sont caractérisés par leur densité ρ et viscosité η. La densité ou masse volumique est la masse du fluide par unité de volume (en kg/m3). Elle est une mesure du nombre de molécules par unité de volume. La viscosité est une mesure de la résistance d’un fluide au changement de forme (en kg/(m⋅s) ou Pa.s); elle détermine la vitesse de mouvement du fluide. Lorsque la viscosité augmente, la capacité du fluide à s'écouler diminue. Les liquides ont une densité et une viscosité supérieures à celles des gaz: les molécules sont plus rapprochées, des liaisons (forces de van der Waals, interactions dipolaires) s'établissent entre elles qui augmentent la cohésion de l'ensemble. Le tableau ci-dessous donne les caratéristiques des atmosphères et agents liquides des planètes solides du Système Solaire.
Fluide | Densité ρ () | Viscosité η ( Pa s) | Vitesse typique (m/s) |
Eau liquide | 1000 | 1540 | 5 |
Méthane liquide | 450 | 184 | |
Atmosphère terrestre | 1.27 | 17.1 | 40 |
Atmosphère martienne | 0.027 | 10.8 | |
Atmosphère de Vénus | 71.92 | 33.0 | |
Atmosphère de Titan | 5.3 | 6.3 | 0.5-1 |
Glace | 992 |
Pour un fluide s'écoulant sur une paroi (le vent ou un agent liquide au dessus d’un sol), la viscosité décrit la contrainte de cisaillement, c’est-à-dire la force tangentielle (par opposition aux forces normales, perpendiculaires à la paroi) qui s’applique par unité de surface sur la paroi. Cette contrainte de cisaillement s’accompagne de l'existence d'un gradient de vitesse d'écoulement dans le fluide. En effet, il existe une couche limite contiguë à la paroi, dans laquelle la vitesse du fluide passe de zéro, au niveau de la paroi, à sa la valeur maximale correspondant à celle d'un fluide libre. L'épaisseur de cette couche limite varie suivant l'état de la surface (plus la surface est lisse, plus la couche est mince). Plus précisément, on peut montrer que la vitesse de l’écoulement croit avec la hauteur au dessus de la paroi de manière logarithmique comme l'illustre la figure ci-contre.
Un fluide peut s’écouler de différentes façons. Lorsque les lignes de courant (c’est-à-dire les tangentes en tous points parallèles à la direction de l'écoulement) restent parallèles entre elles et rectilignes, l’écoulement est dit laminaire. Au contraire, quand l’écoulement est désorganisé et le siège de tourbillons, on parle de régime turbulent (voir la figure ci-contre).
Afin de déterminer le régime en vigueur autour d’un "obstacle" (une roche dans l’eau ou un grain de sable dans le vent par exemple), on définit le nombre de Reynolds: où est le diamètre de l’obstacle et la vitesse terminale ou de sédimentation du fluide (voir ici et là). Si le nombre de Reynolds est grand, les forces inertielles l’emportent sur les forces de frottement liées à la viscosité du fluide ; le régime est turbulent et des tourbillons se développent autour de l’obstacle. Au contraire, si le nombre de Reynolds est petit, en faisant « coller » le fluide à l’obstacle, les forces de viscosité tendent à faire disparaître les tourbillons ; le régime est laminaire.
Les vents sont provoqués par les contrastes de températures à la surface d’une planète et par la rotation de la planète. Les effets de l’activité éolienne sur un paysage sont multiples. Le vent déplace la matière sédimentaire, la redistribue, l’organise parfois en dunes et en modifie les propriétés physiques par abrasion mécanique.
Si l’action destructrice du vent est limitée par rapport à celle d’un agent liquide, les roches se trouvant sur le chemin de saltation (processus de transport des grains par sauts successifs, voir ici) des grains déplacés par le vent sont néanmoins progressivement polies, taillées, striées. L’efficacité de l’abrasion par le vent dépendant de la vitesse d’impact des grains, il est peu probable qu’elle soit élevée sur Titan et Vénus. En revanche, elle devrait l’être sur Mars comme semble le confirmer les images de la surface.
Initier le mouvement d’une particule requiert d’avantage d’énergie que d’entretenir ce mouvement. Dans cet effort, l’idée que plus la particule est petite, plus elle est facile à mettre en mouvement est erronée ; il n’est pas si aisé d’arracher des particules fines à la surface. Ralph Bagnold (1896-1990), grand explorateur des déserts terrestres, a mis en évidence l’existence d’une taille de particule pour laquelle le mouvement à la surface était le plus facile à initier. Ce diamètre seuil peut êtres estimé par la relation empirique suivante : où est la densité du grain, est l’accélération de pesanteur de la planète, la densité de l’air et η sa viscosité. On a pu vérifier sur Terre que la répartition en taille des grains constituant les dunes était globalement distribuée autour de .
La vitesse de cisaillement (qui caractérise la force que le vent exerce sur le sol) seuil de mise en mouvement d’un grain de diamètre est alors : .
La vitesse de l’écoulement (ici du vent), qui varie de manière logarithmique avec la hauteur au dessus du sol (voir la figure), est liée à la vitesse de cisaillement par la formule empirique suivante : où est un facteur lié à la rugosité de surface, environ égal à 1/30 de la taille des grains lorsque ces derniers sont compactés à la surface. Typiquement, sur la Terre et sur Mars, =0.2-0.3 mm. C’est le gradient de vitesse d’écoulement au dessus du sol qui crée un cisaillement et permet de transférer l’énergie du vent vers les grains.
En tout logique, plus la gravité de la planète est faible, plus la taille des particules faciles à déplacer est grande et plus la vitesse seuil de mise en mouvement est petite. Sur Terre comme sur Titan, où les particules sont respectivement constituées de quartz et d’hydrocarbures solides, ∼0.2 mm. Cependant, sur Titan, qui 44 fois moins massif que la Terre, un vent à 1 m de la surface de 0.5 m/s suffit à initier leur mouvement alors que sur Terre un vent minimum de 4.5 m/s est requis. Sur Mars, la gravité est réduite aussi mais la faible densité atmosphérique élève la vitesse de cisaillement seuil de mise en mouvement. L’atmosphère dense de Vénus, au contraire, facilite la mise en mouvement de la matière sédimentaire. L’appliquette Erosion permet d’apprécier la facilité d’entrainement de la matière sédimentaire sur différents corps planétaires.
En dessous de la vitesse seuil , aucune particule, même petite, ne décolle.
En effet, pour les particules dont le diamètre , la vitesse de cisaillement seuil de mise en mouvement décroit en mais reste supérieure à car le déplacement des petites particules est gêné par la présence d’une couche visqueuse tout près de la surface et/ou l’existence de forces cohésives (forces électrostatique ou de van der Waals). Il a fallu plusieurs tempêtes (« dust devils ») sur Mars pour balayer la poussière installée sur les panneaux solaires de Spirit et d’Opportunity.
Pour les grains dont le diamètre , la viscosité de l’air n’est plus un obstacle mais la vitesse seuil de cisaillement nécessaire pour initier leur mouvement croit en selon : . Evidemment, plus le grain est grand et/ou lourd, plus il est difficile à mettre en mouvement. Un air dense, cependant, facilitera son déplacement.
Un bon moyen de prédire comment le vent déplacera un grain de diamètre et de densité , est de comparer la vitesse de cisaillement de l’écoulement à la vitesse limite de chute du grain dans l’atmosphère , aussi appelée vitesse terminale: où est un facteur caratéristique lié à la taille et à la forme du grain. pour les grains de forme sphérique caractérisés par un nombre de Reynold suffisamment grand pour considérer que l’écoulement autour d’eux est de nature turbulente et pour les grains plus petits pour lesquels la viscosité de l’air joue un rôle non négligeable (petit nombre de Reynolds). Pour ces cas là on peut aisément montrer : .
Lorsque la vitesse de cisaillement de l’écoulement est nettement supérieure à , le grain mis en mouvement entre en suspension ; il peut alors s’élever très haut dans l’atmosphère et traverser des distances intercontinentales. La condition étant plus facilement vérifiée pour des grains petits et/ou légers, ce type de mouvement concerne essentiellement les poussières.
Lorsque la vitesse de cisaillement de l’écoulement est comprise entre la vitesse seuil de mise en mouvement et la vitesse limite de chute (), le grain arraché à la surface est entrainé par le vent et se déplace alors par sauts successifs. C’est le phénomène de saltation. Les grains sont soulevées par le vent sur une longueur (qui correspond à une hauteur de quelques centimètres ou dizaines de centimètres sur Terre) et retombent sous l'effet de leur propre poids, en rebondissant et en éjectant d'autres particules par impact. Une fois le mouvement de saltation initié, il requiert un vent moins fort pour être entretenu dans la mesure où, à chacun de ses sauts, le grain transmet une partie de son énergie cinétique à de nouveaux grains. Notons néanmoins que le mouvement de saltation peut être contrarié par le passage sur une surface rugueuse (ou la végétation sur Terre) et la présence d’une force de cohésion entre les grains (notamment la tension superficielle si les grains sont humides).
Enfin, les grains de plus grande dimension et/ou plus massifs, plus difficiles à soulever (), roulent ou glissent à la surface. Ils se déplacent ainsi de proche en proche sans jamais perdre le contact avec le sol. En réalité, leur mouvement est davantage déclenché par l'impact des particules en saltation plutôt que par l'action du vent. C’est la reptation.
Des vents excédant la vitesse seuil de déplacement ont été directement mesurés sur la Terre, Vénus et Titan. Mais cela ne suffit pas à expliquer pourquoi au moins 13% de la surface de Titan, 1.5% des continents terrestres, 0.7% du sol de Mars et 0.004% de la surface de Vénus observée par le Radar de Magellan sont couverts de dunes. Le mécanisme de formation de ces accumulations sableuses est complexe et nous n’en donnerons que le principe ici.
Les premiers grains mis en saltation, en retombant au sol, transmettent une partie de leur énergie à d’autres grains qui peuvent alors plus facilement se mettre en mouvement. La quantité de grains en saltation augmente ainsi progressivement jusqu’à ce que le flux soit saturé c’est-à-dire jusqu’à ce que le fluide ne puisse plus se charger en sable. Commence alors l’accumulation ; le sable érodé en amont se dépose sur le sol. Si la quantité de sable est suffisante, une dune apparaitra.
La morphologie des dunes est très variable (en étoile, barkhane, linéaire, transverse, … etc.); elle dépend non seulement du régime de vent en vigueur (unidirectionnel ou multidirectionnel, permanent ou oscillant) mais aussi de la disponibilité et de la nature des sédiments sur place. Dans les régions où la direction du vent s’inverse périodiquement, à l’échelle d’un jour ou d’une saison, les dunes sont souvent de type linéaire. L’essentiel des dunes observées sur Titan sont de ce type. Elles sont très semblables aux structures observées sur Terre en Namibie ou encore dans le désert d’Arabie Saoudite.
Il est important de souligner que le flux en masse de sédiment transporté par le vent étant proportionnel au cube de la vitesse de cisaillement , l’orientation des dunes est donc plus volontairement contrôlées par des vents rares mais exceptionnellement puissants plutôt que par des vents présents tout au long de l’année mais dont l’intensité est faible. En particulier, il a été récemment montré que les dunes de Titan était sculptées lors de rares tempêtes tropicales, ce qui a permis de résoudre le mystère de leur orientation, contraire à celle des vents moyens prédits par les modèles climatiques mais en accord avec la direction de propagation des vents forts soufflant aux moments des Equinoxes.
Enfin, soulignons qu’il existe d’autres formes édifices éoliens, notamment des yardangs, (crêtes rocheuses sculptée par le vent) et trainées (« wind streaks ») sans expression topographique.
Si la surface de Mars a pu voir l’eau couler dans le passé, aujourd’hui seules les surface de la Terre et de Titan sont soumises à l’érosion pluviale et fluviale. Il pleut, en effet, sur Terre partout où la température l’autorise. Sur Titan, ce sont les cycles du méthane et de l’éthane qui régissent la météo. Lorsque le sol est perméable, les « eaux » des pluies s’y infiltrent, venant parfois enrichir des aquifères - ou « alcanofères » sur Titan- souterrains, jusqu’à le saturer avant de s’écouler à la surface.
Au fur et à mesure, l’écoulement liquide décompose et désagrège le socle rocheux en place (aussi appelé substrat), par altération physique, mécanique mais aussi chimique, et participe ainsi à la production d’une masse sédimentaire ensuite transportée à l’état de grains ou de manière dissoute vers de plus basses altitudes
Altération physique: Plusieurs processus physiques provoquent la fragmentation mécanique du matériel rocheux sans en affecter la composition. En particulier, sur Terre, l’eau, en s’infiltrant, dans les fissures des roches les fragilise et contribue à leur désagrégation par cryoclastie (processus de fragmentation lié au cycle de gel/dégel) ou haloclastie (processus de fragmentation lié à la formation de cristaux de sels suite à l’évaporation de l’eau sur Terre). Ces processus étant liés à des changements de phase de l’agent liquide, ils ne sont sans doute pas très efficaces sur Titan où les variations de température à la surface sont très limitées (<2° pendant la journée, <4° d’une saison à l’autre).
Abrasion mécanique: Les débris solides transportés dans les écoulements sont aussi de puissants agents d’érosion; ils entaillent le substrat rocheux pour creuser des vallées et saper des berges. Un liquide étant plus dense qu’un gaz, l’activité fluviale est un agent d’érosion nettement plus efficace que l’activité éolienne : elle exerce une pression plus forte sur les sols, est capable de transporter des débris plus gros et est davantage aidée dans ses attaques mécaniques par la gravité. Pendant leur transport prolongé, les débris voient généralement leur taille se réduire et leur forme s’arrondir ; 10 km suffisent à façonner un galet de calcaire sur Terre - 300 km pour un galet siliceux. L’ampleur de ce travail d’érosion dépend de la vitesse et de la viscosité du fluide en mouvement ainsi que de la nature (notamment de la dureté) des sédiments et du socle rocheux. Est-ce l’érosion fluviale qui a façonné les pierres arrondies de 2 à 20 cm de diamètre photographiées par la sonde Huygens à la surface de Titan (voir la figure ci-contre)? Et si oui, combien de kilomètres on été nécessaires pour leur donner leur forme ? Des travaux préliminaires sur l’érodabilité de la glace d’eau à -180°C (la matière probable de la croûte de Titan) suggère que l’érosion par les rivières d’hydrocarbure sur Titan est aussi efficace que l’érosion fluviale sur Terre.
Altération chimique:Les activités pluviale et fluviale peuvent aussi modifier la nature chimique du socle rocheux, notamment par :
On l’a vu avec l’activité éolienne, la capacité d’un fluide à mobiliser un sédiment de diamètre d et de densité σ dépend de la densité ρ du fluide, de sa viscosité η et, bien sûr, de sa vitesse d’écoulement. En utilisant une loi d’échelle, l'appliquette "Erosion" permet de comprendre comment les différents paramètres en jeu agissent sur la mise en mouvement de la matière sédimentaire et de comparer l’action d’entrainement du vent sur différents corps planétaires à celle d’un solvant liquide (eau ou méthane liquide). Cependant, sans entrer dans le détail, soulignons que le modèle d’initiation de mouvement présenté dans la section "Activité éolienne: transport (initiation)" s’applique mal au mouvement des petites particules dans un fluide liquide, notamment parce que les forces de cohésion sont moindres et de natures différentes que celles en jeu dans un fluide gazeux. La figure ci-contre montre que des vitesses d’écoulement moins importantes sont requises pour transporter des petites particules dans l’eau ou le méthane liquide.
Une fois mis en mouvement, les sédiments peuvent être transportés. Les modes de transport des sédiments par un agent liquide sont sensiblement les mêmes que ceux vu pour l’activité éolienne: suspension (les particules ne sont jamais en contact avec le fond de la rivière), saltation (les grains se déplace en rebondissant sur le fond) et traction (les grains se déplacent en roulant ou en glissant au fond sans jamais perdre le contact avec le sol). Le transport des sédiments dans un liquide par saltation ou traction est aussi appelé charriage. A cette liste vient néanmoins s’ajouter la possibilité de transporter certains composés sous une forme dissoute (voir ici). Soulignons que l’addition d'une faible quantité de substance en suspension ou en solution peut augmenter grandement la viscosité du liquide.
Pour prédire le mouvement d’un grain dans une rivière, on peut comparer la vitesse de cisaillement du fluide à ce qu’on appelle la vitesse de sédimentation c’est-à-dire la vitesse minimale qu'un flot doit avoir pour transporter, plutôt que déposer, un sédiment de diamètre d et de densité σ. La vitesse de sédimentation est à l’activité fluviale ce que la vitesse terminale est à l’activité éolienne. Elle dépend de la pesanteur, de la taille de la particule, de sa densité et de celle du fluide, et, pour les plus petites particules (celles dont le nombre de Reynolds est inférieur à 1), de la viscosité du fluide.
A l’instar de ce que l’on a vu pour l’érosion éolienne, la vitesse de sédimentation peut s’écrire: pour les petites particules sphériques autour desquels l’écoulement est laminaire et pour les particules sphériques plus grosses dont autour desquels l’écoulement est de nature turbulente.
Si l’écoulement est gravitaire (uniquement produit par l’action de la pesanteur), la vitesse de cisaillement basal (au fond d’un lit de rivière) est lié, en première approximation, à la profondeur du flot h (ou la hauteur des « eaux »), à la pente du lit () et à la gravité par la relation : . Dans la pratique, on peut considérer que si , les grains sont suspendus dans le liquide et si , ils sont charriées (par saltation ou traction) ; les sédiments restent alors confinés dans une zone proche du fond. Plus la pente locale est forte ou plus le niveau des «eaux» est haut (notamment en période de crue), plus la matière sédimentaire sera facilement et abondamment transportée.
Notons que parce que la vitesse de l’écoulement n'est pas constante sur une section de cours d'eau (elle est maximale un peu en-dessous de la surface et dans l'axe du cours d'eau et minimale sur le fond et près des berges), une rivière profonde aura peu d'action sur le fond au contraire d’un écoulement très superficiel (quelques décimètres). Sur Terre, à vitesse égale en surface, la force érosive des wadi (lits de rivières généralement asséchées, en milieu aride) est en effet bien plus forte que celle des rivières des pays tempérés.
La figure ci-contre illustre les modes de transport des gros grains sédimentaires sur la Terre, Titan et Mars. Elle montre notamment que, pour une même vitesse d’écoulement, des sédiments plus gros peuvent être charriés sur Titan par rapport au cas terrestre. D’autre part, du fait de leur densité plus faible, les sédiments composés de glace sont a priori plus faciles à transporter que des sédiments de nature organique.
Les « eaux » fluviales sont généralement collectées au sein de bassins de réception et viennent alimenter un réseau fluviatile hiérarchisé (rigoles, ruisseaux, rivières et fleuves sur Terre). Notons, une différence notable avec le cas éolien : l’agent liquide, contrairement au vent, ne se répartit pas sur toute la surface mais suit la ligne de plus grande pente en restant confiné dans un lit. Les filets d' « eau » confluent et fusionnent en chenaux de taille croissante. La mise en place de ces réseaux fluviatiles dépend de la pente régionale, du débit de liquide et de la nature du substrat (notamment de sa perméabilité). Sur pente forte, les chenaux sont multiples et confluent: le réseau fluviatile est dit en tresse: c'est le cas des portions amonts des cours d'eau (torrents de montagne). Quand la pente devient faible les différents cours d’ « eau » se rejoignent généralement en un unique chenal d'écoulement, souvent sinueux: le réseau à méandres caractérise la plaine alluviale proche de l'embouchure.
La matière sédimentaire mobilisée est déposée à l’endroit où la pente diminue formant ce que l’on appelle des cônes de déjection ou cônes alluviaux. La consolidation des sédiments est à l'origine de la formation des couches sédimentaires mais cette matière peut aussi être re-mobilisée lors de nouveaux épisodes pluvieux. Sur Terre, les sédiments finissent souvent leur voyage dans les océans. C’est peut être aussi le cas sur Titan dont les plus grands lacs sont connectés à des réseaux fluviatiles complexes (voir la figure ci-contre). Les amas de dépôts à l’embouchure des fleuves sont appelés deltas. De part et d’autre du lit limité par les berges, il est aussi fréquent de trouver des levées alluviales, topographies bombées formées par les dépôts des crues. Quelques exemples de paysages fluviaux sur la Terre, Titan et Mars sont présentés sur la figure ci-contre.
Les dessins que forment les réseaux fluviatiles à la surface renseignent sur la nature du substrat rocheux, la pente locale et une éventuelle activité tectonique. Ils fournissent également des informations précieuses sur les climats présents ou passés. Cependant il faut garder à l’esprit que le flux de sédiments transportés variant en , la forme des chenaux d’écoulement est généralement représentative d’événements catastrophiques et notamment d’épisodes de crue. Enfin, soulignons que tous les chenaux ne sont pas d’origine fluviale : les chenaux visibles à la surface de Vénus ont été creusés par des coulées de lave.
Quel point du cours la figure ci-contre illustre-t-elle ?
Une météorite (sphérique) composée de roches denses ayant une densité et un rayon de 1 km frappe la Lune avec une vitesse de .
Calculer l’énergie cinétique de ce projectile ?
Quelle est l’amplitude équivalente M de cet impact sur l’échelle de Richter ? On utilisera la formule: . Qu’en serait-il si la météorite était composée de fer (donc de densité de )?
3) Quelle fraction de cette énergie est nécessaire à la vaporisation totale du projectile ? On considérera que l’énergie de vaporisation est égale à .
Quel est le diamètre du cratère transitoire crée par cet impact ? On supposera que la Lune a la même densité que la météorite et que le projectile arrive sur la surface lunaire avec un angle de 30°.
D’après votre analyse de l’exercice "Catégorie de cratères", s’agit-il d’un cratère simple ou complexe ?
Quel est le volume de matière déplacé sachant que la profondeur du cratère transitoire vaut le tiers de son diamètre ?
Sachant qu’un tiers du volume déplacé est éjecté et redéposé hors du cratère, combien d’impacts de ce type faudrait-il pour recouvrir d’éjectas toute la surface de la Lune sur une épaisseur moyenne de 5 m ?
Cet exercice s'inspire d'un exercice proposé dans "Planétologie" de C. Sotin, O. Grasset et G. Tobie, Edition Dunod, Paris, 2009.
Une météorite de fer ayant une densité de et un diamètre de frappe la Lune avec un angle de et une vitesse de .
Estimer la taille du cratère formé par cet impact.
Des roches sont excavées du cratère avec une vitesse de . A quelle distance du cratère principal se formeront les cratères secondaires ?
Reprendre les questions 1) et 2) pour un impact sur Mercure. Comparer avec le cas lunaire et conclure.
Commentez cette figure extraite de Shkuratov & Bondarenko (Icarus 149, 329, 2001) donnant l’épaisseur du régolithe de différentes régions de la face visible de la Lune en fonction de l’âge de la surface . A votre avis, quel type d’observation a permis d’estimer ?
Triton, le plus grand satellite de Neptune, possède une atmosphère tenue, composée presque uniquement d’azote. Cette atmosphère a probablement comme origine des geysers dont les traces (en l’occurrence des traînées sombres orientées dans le sens du vent dominant résultant de l’éjection puis de la retombée à la surface de panaches de poussières de 2 à 8 km de haut, cf. figure) ont été observées sur la calotte polaire australe du satellite par la sonde Voyager 2. La densité de l’atmosphère de Triton à la pression de surface (∼5 Pa) est de et la viscosité de l’azote à la température de surface (∼38 K) est . L'accélération de pesanteur à la surface de Triton est 0.78 .
Supposons que des grains de glace d’eau soient présents à la surface de Triton.
Quelle taille ont les grains susceptibles d’être déplacé ?
Quelle vitesse doit avoir le vent à 1 m au dessus de la surface pour les déplacer?
Comparez cette vitesse à la vitesse du son dans l’atmosphère de Triton (environ 127 m/s) et concluez sur la probabilité qu’une future mission, dotée de l’instrumentation adéquate, trouve à la surface du satellite des dunes.
A votre avis, quels processus pourraient être à l’origine de la matière sédimentaire sur Triton ?
Cet exercice s'inspire d'un exercice de "Planetary Surface Processes" de H. Jay Melosh, Cambridge University Press, 2011.
Considérons un cours d’eau particulièrement calme s’écoulant à 1 m/s. Le fond de ce cours d’eau est à 1 m de la surface.
Combien de temps faut-il à un grain de sable de 2 µm de diamètre pour atteindre le fond du cours d’eau ?
Même question pour une particule fine d’argile de 0.2 µm de diamètre?
A votre avis, les particules d’argile fines participent-elles à la sédimentation au fond du cours d'eau ?
La sonde Voyager 1 a détecté 9 volcans actifs à la surface de Io. En supposant que chacun de ses volcans a un taux d’éruption de /, déterminer :
La vitesse moyenne de renouvellement de la surface sur Io en cm/an.
Le temps nécessaire au renouvellement total de la surface de Io sur une épaisseur d’1 km.
Cette appliquette illustre les effets d’un impact météoritique sur Mercure, la Terre (avec ou sans atmosphère), la Lune et Mars en fonction des caractéristiques de l’impacteur (vitesse, angle d’arrivée, taille, densité) et de la surface impactée (densité).
Pour le détail des formules à partir desquelles a été construire cette appliquette, se référer à : G. S. Collins, H. J. Melosh, R. A. Marcus: Earth Impact Effects Program: A Web-based computer program for calculating the regional environmental consequences of a meteoroid impact on Earth, Meteoritics & Planetary Science 40, Nr 6, 817–840 (2005).
A vous de jouer en répondant notamment aux questions ci-dessous!
On rappelle que les astéroïdes sont composés de roches et de métaux ; leur densité varie entre 2000 et 8000 et leur vitesse à l’entrée de l’atmosphère terrestre est généralement comprise entre 11 et 21 km/s. Les comètes, quant à elle, sont essentiellement composées de glace ; leur densité est comprise entre 500 et 1500 et leur vitesse à l’entrée de l’atmosphère terrestre est généralement comprise entre 30 et 72 km/s.
En utilisant l'appliquette, vérifiez, quand cela est possible, les résultats des exercices Création du régolithe lunaire par cratérisation et Cratères secondaires.
Retrouvez le diamètre de transition entre cratère complexe et cratère simple pour chaque corps planétaire. Comment évolue-t-il avec la gravité ?
En comparant les sorties de l’appliquette pour les cas « Terre » et « Terre sans atmosphère », déduisez le principal effet de l’atmosphère.
Retrouvez le coefficient de la loi en puissance qui lie le diamètre final d’un cratère simple à l’accélération de pesanteur du corps sur lequel il se trouve.
Estimez le diamètre de la météorite à l’origine de Meteor Crater, en Arizona, sachant que ce cratère a un diamètre d’environ 1.2 km, que l’impacteur était très probablement riche en fer et en nickel, et qu’il a frappé la Terre avec un angle d’environ 80°. Vérifiez que l’ordre de grandeur théorique de la profondeur finale du cratère est en accord avec la réalité.
Mêmes questions pour le cratère lunaire complexe Tycho dont le diamètre est de 85 km sachant que l’impacteur qui lui a donné naissance avait une trajectoire basse au dessus de l’horizon (c’est-à-dire avec un angle d’arrivée d’au moins 45°) et que sa densité était proche de celle de la Lune.
En raison de la grande volatilité de la glace, les surfaces glacées sans atmosphère sont soumises à un phénomène de migration/ségrégation thermique. L’appliquette ci-dessous vise à évaluer l’efficacité de ce phénomène par rapport à d’autres processus d’érosion spatiale sur les principaux satellites glacés du Système Solaire. Elle s’inspire du travail de thèse de J.R. Spencer: The surfaces of Europa, Ganymède, and Callisto- An investigation using Voyager IRIS Thermal Infrared Specta, Ph.D dissertation by John R. Spencer, 1999.
Lisez l’essentiel à savoir ci-dessous et essayez de répondre aux questions.
Le taux instantané de sublimation de la glace peut, en première approximation, s’exprimer de la façon suivante : où est la densité volumique de la glace (), est la masse molaire de l'eau (), est la constante universelle des gaz parfaits () et , la température instantanée (en K).
s’obtient en égalisant le flux solaire (entrant) et le flux émis par la surface (sortant) : où est la constante de Stefan-Boltzmann (), est la constante solaire (c’est-à-dire la puissance reçue du Soleil par unité de surface normale aux rayons solaires à la distance héliocentrique de 1 UA) (), est la distance héliocentrique en UA du corps glacé, est est l’albédo de la surface et est l’angle d’illumination du Soleil à la surface (l’angle entre la normale à la surface et la direction de l’ensoleillement). Il dépend de la latitude, de l’heure locale et éventuellement de la saison. Ici on considère que .
est la pression de vapeur saturante de la glace en Pascal, c’est-à-dire la pression à laquelle la phase gazeuse de l’eau est à l’équilibre avec sa phase solide à la température . Dans la gamme de températures des satellites glacés du système solaire (130-150 K), il a été établi semi-empiriquement que : .
Mise en jambe : par une analyse dimensionnelle, retrouvez la dimension de ?
Comparez l’amplitude du phénomène de ségrégation thermique entre les satellites galiléens (), les satellites saturniens () et ceux d’Uranus (). Vous vous placerez à l’Equateur, à midi, en été et prendrez un albédo de pour la glace équatoriale.
Tracez le taux de sublimation lié au phénomène de ségrégation thermique en fonction de la latitude et de l’albédo de la glace pour un satellite galiléen. Comparez son intensité sur Europa (), Ganymède () et Callisto ().
À quelle(s) latitude(s) ce phénomène est-il le plus actif ?
Dans le système de Jupiter, la vitesse de « laboure » des régolithes par impacts micro-météoritiques est de quelques mm/an. L’intensité du phénomène de « sputtering », quant à elle, décroit avec la distance à Jupiter : mm/an sur Europe, mm/an sur Ganymède et mm/an sur Callisto. Que peut-on en déduire sur l’efficacité du phénomène de ségrégation thermique sur Europe, Ganymède et Callisto ? Discutez.
L’activité éolienne et, lorsqu’elle existe (sur Terre et Titan), l’activité fluviale participent efficacement au transport des sédiments des surfaces planétaires dotées d’une atmosphère. L’initiation de ce mouvement se fait plus ou moins facilement en fonction de la nature des sédiments (densité), des caractéristiques du fluide (densité et viscosité de l’air ou de l’agent liquide) et de la pesanteur. C’est ce qu’illustre cette appliquette.
Le modèle simplifié sur lequel s’appuie cette appliquette est décrit en partie dans la section activité éolienne. Pour plus de détails se référer à : « Planetary Surface Processes » de H. Jay Melosh, Cambridge University Press, 2011.
A vous de jouer en essayant de répondre aux questions ci-dessous.
Mise en mouvement des sédiments sur Vénus, la Terre, Mars et Titan
Les caractéristiques (densité, viscosité) des fluides (atmosphères, liquides) présents à la surface de Vénus, la Terre, Mars et Titan sont donnés ici. Rappelons que la composition des sédiments varie d’un corps planétaire à l’autre : du quartz (2650 ) sur Terre, du basalte (2900 ) sur Vénus et Mars et de la glace d’eau (992 à 92 K) et/ou de la matière organique (1500 ) sur Titan.
Testez l’appliquette pour un cas quelconque et expliquez la courbe obtenue et notamment les 2 régimes qui se dégagent.
Sur quel objet planétaire le transport des particules sédimentaires par le vent est-il le plus facile ? Classez les planètes par ordre de facilité du transport éolien et expliquez.
Toutes ces planètes présentent des dunes à leur surface, sur laquelle les grains constituant ces dunes sont sans doute les plus fins ?
Trouvez la combinaison Planète-Atmosphère-Sédiments pour laquelle l’entrainement de la matière sédimentaire à la surface requiert les vents les plus faibles.
La taille typique des grains constituant les dunes sur Terre est de l’ordre de 200 μm. Commentez.
Comparez la mise en mouvement d’un grain à la surface par le vent à celle d’un grain situé au fond d’une rivière sur Terre. Expliquez.
Même question pour Titan. Notez que la composition des sédiments sur Titan n’est pas encore bien identifiée. Il pourrait s’agir de matière organique, de glace d’eau ou d’un mélange de ces composés. L’agent liquide est le méthane liquide.
Comparez le transport fluvial sur Titan à celui sur Terre.
Dans le passé, de l’eau liquide coulait sans doute à la sa surface de Mars. Comparez l’activité paléo-fluviale de Mars à celle présente de la Terre. Que se passe-t-il pour les grains les plus petits ?
Quelle doit être la vitesse minimum du vent à 1 m du sol pour la mise en mouvement de grains à la surface de Mars. Qu’en déduisez vous ?
-------------FIN DU CHAPITRE -------------FIN DU CHAPITRE -------------FIN DU CHAPITRE -------------FIN DU CHAPITRE -------------
Présentation du module
Prérequis
Dans cette partie, nous verrons comment s'établit la structure thermique dans les couches externes fluides (principalement gazeuses, ce que l'on appelle les atmosphères) des planètes, ainsi que les conséquences de l'existence d'une telle structure.
Il est possible de parcourir la partie Découvrir avec un simple bagage de Terminale scientifique ou d'amateur de vulgarisation scientifique. En revanche, la bonne compréhension des phénomènes en jeu et la capacité à calculer même approximativement les conditions moyennes au sein d'une atmosphère planétaire exigent un bagage en physique générale niveau licence, à savoir plus précisément :
L'observation de la seule biosphère connue à jour (celle de la Terre) conduit les exobiologistes à poser comme nécessaire la présence d'eau liquide (ou au moins d'un liquide aux propriétés analogues comme l'ammoniac) à la surface d'une planète pour qu'une chimie prébiotique complexe, puis une activité biologique au sens propre, puisse s'y développer. Si bien que la notion d'habitabilité planétaire est de nos jours quasiment devenue un synonyme de présence possible d'eau liquide.
Or, si la disponibilité de l'eau dans l'Univers ne fait guère de doutes (la molécule H2O étant l'une des plus répandues), la question de son apport sur les planètes telluriques fait encore l'objet de débats. Surtout, la permanence de son état liquide est encore plus difficile à obtenir, et nécessite une fourchette de conditions de pression et de température bien particulières (ainsi, à la pression atmosphérique terrestre, doit-on se trouver entre 0°C et 100°C pour que l'eau puisse demeurer liquide). Les conditions de pression et de température au sein des atmosphères planétaires de leur sommet jusqu'à l'éventuelle surface constituent donc l'un des facteurs essentiels conditionnant les phénomènes pouvant s'y dérouler (qu'ils soient de nature biologique, ou plus simplement chimique ou météorologique).
Une autre question cruciale est celle de la détectabilité de telles planètes dans notre voisinage galactique. Le seul moyen envisageable pour caractériser ces planètes consiste en l'étude spectroscopique (c'est-à-dire, décomposé selon ses différentes "couleurs") du rayonnement qui nous parvient. Ce rayonnement peut nous parvenir principalement par deux processus physiques distincts :
Il existe trois modes de transport de la chaleur au sein des atmosphères planétaires, qui déterminent les variations de température au sein de ces atmosphères :
Une des plus importantes caractéristiques des atmosphères planétaires est leur épaisseur verticale. En toute rigueur, leur densité décroît continûment avec l'altitude jusqu'à rejoindre celle du milieu interplanétaire, si bien qu'il est difficile de leur attribuer une épaisseur bien définie. On peut néanmoins caractériser la rapidité avec laquelle cette densité décroît avec l'altitude (atmosphère plus ou moins bien "tassée" verticalement). Cela définit ce que l'on appelle l'échelle de hauteur atmosphérique, qui représente la différence d'altitude entraînant une division de la pression atmosphérique (liée à la densité) par un facteur constant (). Le lecteur intéressé par une définition quantitative pourra se reporter ici.
Cette échelle de hauteur est le résultat du compromis entre deux phénomènes physiques : la gravitation qui tend à tasser les molécules de l'atmosphère vers le bas, et l'agitation thermique des molécules qui tend à les disperser dans tout l'espace, y compris vers le haut. À ce titre, et toutes choses égales par ailleurs, l'échelle de hauteur atmosphérique est :
Planète (ou satellite) | Vénus | Terre | Mars | Jupiter | Io | Saturne | Titan | Uranus | Neptune | Triton | Pluton |
---|---|---|---|---|---|---|---|---|---|---|---|
Échelle de hauteur (en km) | 16 | 8,4 | 11 | 25 | 7,9 | 48 | 21 | 27 | 22 | 14 | 18 |
Dans le système solaire, les valeurs extrêmes vont de pour la Terre à environ pour Saturne. Ces valeurs sont en général très petites devant le rayon de la planète, si bien que l'on peut négliger la courbure de la planète et considérer l'atmosphère comme une succession de couches planes empilées de bas en haut. C'est ce que l'on appelle l'approximation plan-parallèle.
Là où la convection est le mode de transport dominant d'énergie dans une atmosphère, on constate une décroissance régulière de la température avec l'altitude selon un coefficient (en °C/km ou K/km) appelé gradient adiabatique. En effet, si l'on considère une masse de gaz au cours de son transport dans un courant de convection vertical, celle-ci devra lutter contre la pesanteur et donc fournir de l'énergie pour ce faire. Or, le seul "réservoir" d'énergie d'un gaz parfait réside dans sa capacité calorifique. Il y aura donc une conversion partielle de son énergie thermique (en fait, de son enthalpie puisqu'on y inclut le travail des forces de pression) vers de l'énergie potentielle de pesanteur, et donc une baisse de la température de la parcelle d'air d'autant plus grande que celle-ci aura acquis davantage d'altitude (voir ici pour la démonstration). Si une région de l'atmosphère est soumise à cette circulation et en négligeant les autres modes de transport d'énergie, la température y décroît alors avec l'altitude en suivant ce gradient adiabatique.
En pratique cependant, les atmosphères planétaires ne sont pas constituées que de gaz parfaits, mais comportent des gaz en équilibre avec leur propre phase condensée (liquide ou solide). C'est le cas par exemple sur Terre de la vapeur d'eau qui constitue une proportion variable de l'atmosphère terrestre et se trouve parfois en équilibre avec des gouttes d'eau liquide ou des cristaux de glace d'eau. Ou encore de Titan où c'est le méthane gazeux qui se trouve parfois au contact de gouttes de méthane liquide. En ce cas, il existe un réservoir d'énergie supplémentaire pour une parcelle d'atmosphère en mouvement ascendant, à savoir l'énergie libérée par le gaz condensable lorsqu'il se convertit en gouttelettes liquides ou en cristaux solides, ce que l'on appelle la chaleur latente de condensation. Ce réservoir supplémentaire d'énergie limite alors la baisse de température avec l'altitude vers une valeur plus faible. On parle alors de gradient adiabatique humide, que l'on distingue du gradient adiabatique sec en l'absence de condensation.
La couche atmosphérique où la convection est le mode dominant de transport d'énergie s'appelle la troposphère, caractérisée par la décroissance en température décrite ci-dessus. C'est la couche atmosphérique la plus profonde, au contact de la surface pour les planètes telluriques. Au-dessus de la troposphère, les densités plus faibles rendent le transport d'énergie par rayonnement comparativement plus efficace que la convection, car le milieu dilué devient davantage transparent au rayonnement thermique.
Le corps noir est un objet physique idéal qui absorbe tout le rayonnement électromagnétique qu'il reçoit (sa réflectivité est donc nulle à toutes les longueurs d'onde).
La propriété fondamentale du corps noir est que l'intégralité du rayonnement électromagnétique en provenance de cet objet est d'origine thermique. Le spectre de ce rayonnement ne dépend alors que de la température du corps noir en question. En particulier :
Certains objets réels sont de bonnes approximations du corps noir idéal, du moins sur certains intervalles de longueur d'onde et dès que le rayonnement réfléchi y est négligeable devant l'émission thermique et en l'absence de processus d'émission autres que thermiques. C'est par exemple le cas de la plupart des objets du quotidien dans le domaine infrarouge moyen (pour les longueurs d'onde autour de .), ou encore des étoiles dans le domaine visible.
Il est d'usage de distinguer deux intervalles spectraux différents lorsque les planètes ont une température notablement plus faible que leur étoile (ce qui est toujours le cas dans le système solaire, mais pas toujours pour les planètes extrasolaires !).
La température d'équilibre d'une planète est la température théorique de sa surface (si on suppose cette température uniforme) en l'absence d'atmosphère. C'est une grandeur théorique qui n'a pas vocation à être mesurée, contrairement à la température effective.
La température d'équilibre se détermine à partir d'un simple bilan de rayonnement (visible et thermique). Cela revient à négliger les autres sources d'énergie que le rayonnement de l'étoile hôte (le Soleil pour la Terre par exemple) : géothermie, réactions chimiques ou nucléaires, transitions de phase, etc. Sont pris en compte :
La température de surface influe ici sur le rayonnement thermique. Elle est égale à la température d'équilibre lorsque le bilan est équilibré, à savoir : Puissance lumineuse reçue = Puissance lumineuse réfléchie + Puissance rayonnée thermiquement, ce qui est équivalent à Puissance lumineuse absorbée = Puissance rayonnée thermiquement.
Une remarque importante est que cette définition repose sur l'hypothèse irréaliste d'une température de surface homogène sur l'ensemble de la planète, donc avec une redistribution parfaite de l'énergie. Cette température est donc un outil théorique plus qu'une température physiquement mesurable. Le lecteur intéressé par une approche plus quantitative (mais identique conceptuellement) pourra se reporter ici.
La température effective est une mesure de la puissance émise thermiquement par une planète. Elle se définit comme la température du corps noir (idéal) émettant la même puissance totale (en comptant toutes les longueurs d'onde) que la planète par unité de surface. Contrairement à la température d'équilibre, c'est une grandeur expérimentalement mesurable.
Nous connaissons assez bien le système solaire pour pouvoir mesurer les températures effectives des planètes et les comparer aux températures d'équilibre théoriques. Le résultat est résumé sur le tableau ci-dessous :
Planète (ou satellite) | Mercure | Vénus | Terre | Lune | Mars | Jupiter | Saturne | Titan | Uranus | Neptune |
---|---|---|---|---|---|---|---|---|---|---|
Température d'équilibre (°C) | 161 | -42 | -19 | -2 | -63 | -163 | -192 | -191 | -215 | -227 |
Température effective (°C) | 161 | -42 | -19 | -2 | -63 | -149 | -178 | -191 | -214 | -214 |
Température moyenne de surface (°C) | 161 | 462 | 15 | -2 | -58 | N/A | N/A | -179 | N/A | N/A |
Pour la plupart des planètes extrasolaires (hormis les plus grosses et les plus chaudes), seule la température d'équilibre peut être estimée (en supposant un albédo donné par un modèle théorique). Les ordres de grandeur de ces deux températures sont comparables lorsque la source d'énergie principale de l'atmosphère est le rayonnement de l'étoile hôte, comme c'est le cas dans le système solaire. Pour les planètes telluriques (et le satellite de Saturne Titan), ces deux températures sont mêmes égales car les sources d'énergie interne à la planète ont un effet négligeable sur l'atmosphère, ce qui n'est pas le cas pour les géantes gazeuses.
On constate également sur le tableau précédent que pour les corps possédant une surface solide, la température moyenne de la surface est toujours au moins égale à la température effective (égale pour un corps sans atmosphère comme la Lune ou bien Mercure, supérieure pour ceux possédant une atmosphère). Ce phénomène est appelé effet de serre et sera expliqué plus en détail à la page suivante.
Le phénomène essentiel à l'origine de l'effet de serre au sein d'une atmosphère réside dans la différence d'absorption des rayonnements infrarouge thermique (en provenance de la planète) et visible/UV (en provenance de l'étoile) par les constituants de l'atmosphère. Les constituants gazeux d'une atmosphère (en excluant les particules solides ou liquide en suspension comme les poussières ou les cristaux et gouttelettes des nuages) sont en général transparents pour la lumière visible émise par leur étoile. En revanche, certains de ces gaz (comme la vapeur d'eau H2O, le dioxyde de carbone CO2 ou encore le méthane CH4) absorbent très bien le rayonnement infrarouge d'origine thermique émis par la planète.
Cette différence d'absorption entre les rayonnements conduit à une séparation entre les régions :
Or, le bilan d'énergie de la planète impose que ce soit la couche rayonnant vers l'espace qui soit à la température effective permettant un équilibre entre la puissance reçue et celle absorbée. Il faut donc que l'énergie absorbée en profondeur puisse être transportée jusqu'à cette altitude de rayonnement. Comme l'atmosphère profonde située entre ces deux niveaux est opaque aux infrarouges, le rayonnement n'est pas le mode le plus efficace de transport, et c'est la convection qui prend le relais. Cette atmosphère profonde, s'étendant depuis l'altitude d'émission infrarouge jusqu'en bas (surface ou intérieur planétaire pour les géantes) n'est autre que la troposphère définie précédemment. Afin que ce transport d'énergie par convection puisse avoir lieu, il faut que la température de surface soit plus élevée que celle au sommet de la troposphère selon le gradient adiabatique. La température au sommet de la troposphère étant égale à la température effective, la température de surface est en conséquence plus élevée, ce qui est la définition même de l'effet de serre.
Ce sont les couches atmosphériques situées au-dessus de la troposphère, où la convection joue un rôle négligeable.
La couche atmosphérique située au-dessus de la troposphère est (en général, voir page suivante) appelée mésosphère. Le transport d'énergie s'y fait exclusivement par rayonnement. La température y décroît avec l'altitude en tendant vers une valeur appelée température de peau atmosphérique. Cette décroissance s'y effectue de façon beaucoup plus modérée que dans la troposphère située en dessous et soumise au gradient adiabatique.
Au sommet de la mésosphère, l'atmosphère devient complètement transparente à tous les rayonnements (les rayonnements visible ou IR thermique ne peuvent donc y déposer leur énergie) et extrêmement ténue (la convection est donc inefficace). Le transport d'énergie y est donc assuré faute de mieux par des processus de conduction qui sont eux-mêmes très inefficaces à grande distance. Cette zone connaît donc d'énormes contrastes de température verticaux et horizontaux car l'énergie qui y est déposée par les particules énergétiques de l'espace interplanétaire ou les rayonnements X et γ de l'étoile s'évacue très difficilement, ce qui conduit à l'appellation de thermosphère. La température y croît avec l'altitude, comme montré plus en détail ici.
Certaines atmosphères planétaires possèdent une couche supplémentaire appelée stratosphère située entre la troposphère et la mésosphère. Cette couche est une couche radiative (la convection n'y joue aucun rôle dans le transport vertical de la chaleur) et connaît une inversion de température : la température y croît avec l'altitude ! Cette inversion est causée par une absorption partielle de la lumière et/ou des UV stellaires assez haut dans l'atmosphère, si bien que cette énergie ne peut pas s'évacuer par convection et seulement difficilement par radiation. Il se crée alors une anomalie chaude qui déforme le profil de température, allant jusqu'à l'inversion de température.
Dans le système solaire, Vénus et Mars ne possèdent pas de stratosphère (ces atmosphères principalement constituées de CO2 rayonnent très efficacement en infrarouge le peu de puissance absorbé à haute altitude, si bien que les anomalies de températures n'altèrent pas la forme du profil thermique). La Terre en possède une, causée par l'absorption des UV solaires par l'ozone (O3), sous-produit du dioxygène (O2) d'origine biologique. Les planètes géantes en possèdent également (causée par des composés hydrocarbonés absorbant les UV) ainsi que Titan (par absorption des UV solaires sur les particules du brouillard photochimique produit dans la haute atmosphère).
Nous allons à présent aborder les lois quantitatives permettant de modéliser simplement les profils verticaux de température au sein des atmosphères planétaires. Cela nécessite quelques rappels sur le rayonnement thermique, dit de "corps noir".
L'intensité lumineuse , définie comme la puissance émise par unité de surface émettrice, par angle solide autour de la direction du rayon et par unité de longueur d'onde émise par tout corps noir idéal de température , est donnée par la loi de Planck :
où , et désignent respectivement les constantes fondamentales de Planck, de la vitesse de la lumière et de Maxwell-Boltzmann. Cette fonction possède des propriétés mathématiques aux conséquences importantes pour la suite du cours.
Elle donne la position du maximum en de à température donnée, comme illustré précédemment.
Autrement dit, plus le corps est chaud, et plus il émet principalement à des longueurs d'ondes courtes et ce de façon inversement proportionnelle. Cela justifie la séparation du spectre lumineux en :
La séparation entre les deux domaines est prise de façon conventionnelle autour de . Dans le contexte exoplanétaire, une remarque importante s'impose dès maintenant : la plupart des exoplanètes actuellement connues sont extrêmement chaudes, avec des températures excédant souvent , si bien que la limite entre infrarouge thermique et lumière stellaire est décalée vers de plus courtes longueurs d'onde, voire devient complètement dénuée de sens. Cela empêche notamment d'appliquer tels quels les modèles atmosphériques conçus dans le système solaire qui distinguent ces deux catégories.
Lorsque l'on ne s'intéresse pas au détail du spectre émis par le corps noir, il est souvent intéressant de calculer le flux (c'est à dire la puissance par unité de surface émettrice) total émis par le corps noir dans un demi-espace (par exemple, pour une surface planétaire, vers le haut). Pour cela, il suffit d'intégrer la loi de Planck sur sa variable spectrale , et sur les d'angle solide en question. Le calcul donne alors le résultat suivant, connu sous le nom de loi de Stefan-Boltzmann :
où est connu sous le nom de constante de Stefan-Boltzmann. La puissance émise par un corps noir dépend donc énormément de sa température (une augmentation relative de de sa température entraîne ainsi une augmentation d'environ du flux émis).
Le corps noir est un modèle abstrait que l'on ne rencontre pas dans la vie courante. Le spectre thermique émis par un corps donné se trouvant à l'équilibre thermodynamique à la température peut alors s'exprimer comme où est une grandeur sans dimension appelée émissivité (qui dépend de la température, mais de façon moins marquée que la fonction de Planck si bien que par souci d'alléger les notations, on ne la note pas en général comme on le devrait en toute rigueur).
Considérons un corps noir en contact radiatif avec un corps réel à travers un filtre laissant seulement passer les radiations à la longueur d'onde . On sait qu'une fois l'équilibre thermodynamique atteint, ces deux corps en contact radiatif auront la même température . Si l'on note la fraction du rayonnement incident absorbée par le corps réel, que l'on appelle absorbance, il en renvoie la fraction complémentaire . Un bilan net des flux (nul à l'équilibre) à travers le filtre donne alors la relation , ce qui se simplifie en . C'est la loi de Kirchhoff, que L'on résume souvent en "les bons absorbeurs sont les bons émetteurs".
Cette page développe de façon quantitative les notions vues de façon qualitative ici.
Le bilan radiatif à l'équilibre imposant l'égalité entre la puissance rayonnée par la planète et la puissance absorbée par la planète, on obtient alors l'équation suivante :
qui se résout directement, après simplification du rayon de la planète (ce qui signifie qu'en première approximation, la température d'une planète ne dépend pas de sa taille) en :
ce qui permet de constater que cette température décroît avec la distance à l'étoile, et est proportionnelle à celle de l'étoile. Ainsi, toutes choses égales par ailleurs, pour une étoile naine rouge d'une température moitié de celle du Soleil, il faut pour conserver une température d'équilibre donnée se rapprocher de l'étoile d'un facteur quatre : on peut d'ores et déjà affirmer que les zones habitables autour des petites étoiles de faible température (naines rouges) sont très proches de ces dernières. Notons au passage que la température d'équilibre d'une planète est bornée par celle de son étoile, plus précisément comprise entre (à très grande distance) et à la limite où l'orbite de la planète est tangente à son étoile (et la planète de rayon négligeable devant l'étoile).
Cette équation relie l'augmentation de la pression en descendant avec la masse volumique locale (autrement dit, elle exprime le fait que l'origine physique de la pression au sein des atmosphères est le poids de la colonne de gaz située à la verticale). La différence de pression entre le haut et le bas d'une couche d'épaisseur (la direction verticale étant bien définie en géométrie plan-parallèle) dépend donc de la masse contenue dans un volume de section horizontale et d'épaisseur , d'où, par équilibre des forces verticales s'exerçant sur ce volume
Une simplification par fait donc apparaître : la pression décroît bien avec l'altitude, selon la masse volumique et la gravité locales.
La forme habituelle de cette équation , où désigne la constante des gaz parfaits, P la pression, V le volume occupé, n le nombre de moles et T la température n'est pas vraiment adaptée à une formulation locale (intensive, dirait-on en thermodynamique). Il vaut mieux la présenter sous la forme , où l'on voit apparaître la densité molaire (homogène à des ) locale. Or, cette grandeur est proportionnelle à la masse volumique selon la relation où désigne la masse molaire. Il est alors possible d'exprimer la masse volumique du gaz parfait en fonction des conditions de pression et température locales, ainsi que de la masse molaire du gaz constituant : .
On suppose ici que l'atmosphère est constituée d'un gaz parfait de masse molaire , et que l'atmosphère est de surcroît isotherme à la température selon l'altitude. L'utilisation de l'équation d'état du gaz parfait au sein de l'équilibre hydrostatique donne, par substitution de , avec désignant une grandeur homogène à une altitude. On l'appelle échelle de hauteur, et son interprétation est plus claire en intégrant l'équation différentielle où elle apparaît, avec la condition à la limite inférieure :
L'échelle de hauteur représente donc la hauteur caractéristique avec laquelle la pression décroît avec l'altitude pour tendre vers dans l'espace interplanétaire à grande distance de la planète (mais l'approximation plan-parallèle, ainsi que la thermodynamique usuelle à l'équilibre cessent d'être valides à quelques dizaines d'échelles de hauteur au-dessus de la surface).
Dans le cas d'une atmosphère non isotherme, la résolution formelle est un peu plus complexe, mais l'idée générale d'une décroissance localement exponentielle selon une échelle de hauteur locale dépendant de la température locale reste valable.
On peut reformuler la constante des gaz parfaits selon où désigne la constante de Maxwell-Boltzmann et la constante d'Avogadro, puis simplifier dans l'expression de . On obtient alors où désigne la masse d'une molécule de gaz : une molécule de gaz à la hauteur caractéristique possède donc une énergie potentielle de pesanteur du même ordre que son énergie cinétique microscopique (thermique) moyenne. On comprend donc bien pourquoi représente le compromis entre l'agitation thermique qui tend à disperser les atmosphères ( est croissant avec ), et le poids qui a tendance à tasser les atmosphères vers le bas : décroît avec (atmosphère dense) et (gravité forte).
Dans un modèle purement radiatif d'une colonne d'atmosphère (sans convection ni conduction), il est relativement facile d'estimer l'effet de serre causé par une atmosphère (transparente en lumière visible et partiellement opaque au rayonnement infrarouge thermique) entourant une planète tellurique.
On supposera que la surface possède une émissivité égale à en infrarouge thermique, et que celle de l'atmosphère (directement reliée à son absorbance via la loi de Kirchhoff) est prise constante et égale à dans tout le domaine infrarouge thermique (c'est ce que l'on appelle l'approximation grise). L'atmosphère est considérée ici isotherme à la température . On négligera aussi les flux d'énergie éventuels provenant de l'intérieur de la planète, et on supposera que l'étoile émet de façon négligeable dans l'infrarouge thermique, situé loin de son maximum d'émission dans le visible (ou le proche IR pour les plus froides d'entre elles).
La situation est très simple pour les flux stellaires. désigne le flux moyen à la surface de la planète, qui se déduit du flux à incidence normale appelé constante solaire (ou stellaire) par l'égalité des puissances : (voir le raisonnement définissant la température d'équilibre pour plus de détails, désigne ici le rayon planétaire). On en déduit immédiatement : un facteur 2 s'explique aisément par le fait que seul un hémisphère est éclairé, et l'autre facteur 2 par la moyenne du cosinus de l'angle d'incidence intervenant dans le calcul local du flux.
En vertu de la définition de l'émissivité, l'atmosphère rayonne donc dans chacun des demi-espaces inférieur (vers la surface) et supérieur (vers l'espace). En vertu de la loi de Kirchhoff, cette émissivité est égale à son absorbance, si bien que la fraction complémentaire du rayonnement en provenance de la surface (considérée comme un corps noir) réussit à la traverser, le reste étant absorbé (on néglige les processus de diffusion ici ; seules les émissions et absorptions sont prises en compte).
Le bilan des flux à la surface donne alors à l'équilibre radiatif (synonyme d'égalité entre la somme des flux entrants et la somme des flux sortants) : , tandis que celui au niveau de la couche atmosphérique donne . Nous avons donc deux équations pour les deux inconnues et , et la résolution du système donne alors : et où l'on aura reconnu la température d'équilibre définie précédemment.
Le modèle vu précédemment a l'inconvénient de ne pas pouvoir excéder une augmentation de température à la surface de . Ceci est insuffisant dans le cas des atmosphères très épaisses comme celle de Vénus, où le rapport excède ! Cela signifie que de telles atmosphères ne peuvent se modéliser par une unique couche isotherme, même totalement absorbante aux rayons infrarouges. Il existe différents modèles plus complexes permettant de mieux rendre compte des effets de serre intenses.
Une première idée est d'ajouter, au-dessus de la première couche atmosphérique complètement opaque au rayonnement thermique de la planète, une ou plusieurs couches (la dernière couche immédiatement avant l'espace pouvant être partiellement transparente). Ces différentes couches atmosphériques peuvent alors chacune adopter des températures différentes, et former ainsi un profil de température décroissant avec l'altitude. Il faut ainsi environ une centaine de couches opaques pour rendre compte de la température de surface de Vénus.
L'étude d'un modèle à deux couches atmosphériques fait l'objet du mini-projet associé à ce chapitre.
Une vision plus réaliste mais ne faisant toujours intervenir que des échanges d'énergie par rayonnement consiste à découper l'atmosphère en un mille-feuille constitué d'une infinité de couches atmosphériques infiniment fines (d'un point de vue radiatif). En restant dans l'approximation grise en infrarouge thermique et transparente en lumière visible, il est même possible (mais hors-programme au niveau licence) de démontrer l'expression du profil de température en fonction de la profondeur optique en infrarouge thermique : . Notons que dans ce modèle, on obtient : le seul équilibre radiatif tend à créer une discontinuité de température au niveau de la surface, ce qui déclencherait alors des processus de convection pour y remédier. Un tel contraste thermique est néanmoins observable à la surface des planètes telluriques éclairées par le Soleil, comme une plage sur Terre par beau temps (le sable peut alors être brûlant et l'air frais), ou mieux encore dans les déserts de Mars.
Néanmoins, dans les atmosphères épaisses ou pour expliquer l'existence des stratosphères, l'absorption de la lumière stellaire par l'atmosphère doit être prise en compte (par exemple, seuls quelques pourcents de la lumière solaire atteint directement la surface de Vénus). Des expressions analytiques deviennent alors délicates à trouver, mais des modèles numériques peuvent être utilisés pour déterminer les profils de température dans une colonne d'atmosphère (ce que l'on appelle un modèle 1D radiatif). On peut également profiter de la puissance de calcul des ordinateurs pour abandonner d'autres approximations : il est par exemple indispensable d'abandonner l'approximation grise en infrarouge thermique si l'on veut simuler le spectre du rayonnement thermique émis par la planète.
Lorsqu'une parcelle de gaz se déplace verticalement de façon adiabatique, sa température varie sous l'effet des variations de pression. Si l'on considère un déplacement élémentaire entre l'altitude et d'une masse d'un gaz de capacité calorifique à pression constante Cp, d'entropie S et de volume V à la pression P et à la température T, un bilan de son enthalpie H donne puisque (déplacement adiabatique). La variation de pression étant reliée au déplacement vertical selon la loi hydrostatique , on obtient alors en posant , appelé gradient adiabatique (sec). La détente adiabatique d'une parcelle de gaz ascendante conduit donc à un refroidissement proportionnel à la différence d'altitude selon le gradient adiabatique.
Une autre façon, peut-être plus intuitive, de considérer ce phénomène est d'interpréter la relation intermédiaire obtenue : la variation d'enthalpie du gaz (son "énergie thermique" en tenant compte des forces de pression) est directement reliée à sa variation d'énergie potentielle. Faire monter une parcelle de gaz lui coûte de l'énergie potentielle, ce qui est prélevé sur l'énergie thermique interne de ce gaz en l'absence de chaleur communiquée depuis l'extérieur.
Certaines atmosphères comportent des espèces chimiques condensables. L'exemple par excellence est la vapeur d'eau sur Terre, qui peut se condenser en glace ou un eau liquide. On rencontre aussi ce cas de figure sur Titan avec cette fois le méthane, ou encore dans les atmosphères des géantes gazeuses au niveau de leurs couches nuageuses. Le bilan précédent doit alors être modifié pour tenir compte de la libération de chaleur latente causée par le changement d'état qui peut arriver lorsque l'espèce condensable est saturée.
Le nouveau bilan d'enthalpie est alors donné par où désigne l'enthalpie massique de condensation et la masse d'espèce volatile qui se condense au cours du déplacement au sein de la parcelle de gaz. On arrive alors à l'expression suivante pour le nouveau gradient adiabatique (dit gradient adiabatique humide) : où désigne la fraction massique du volatil au sein de la parcelle de gaz. On constate alors que est plus faible que en valeur absolue : la libération de chaleur latente par liquéfaction ou condensation compense partiellement le refroidissement dû à l'ascension.
Vénus | Terre | Mars | Jupiter | Saturne | Uranus | Neptune | Titan | |
---|---|---|---|---|---|---|---|---|
-10.5 | -9.8 | -4.5 | -2 | -0.71 | -0.67 | -0.85 | -1.3 | |
-5 | -0.5 |
La comparaison entre le profil thermique à un instant donné et le gradient adiabatique au même endroit permet de connaître la stabilité de l'atmosphère vis-à-vis des phénomènes de convection. Supposons pour bien comprendre un profil thermique isotherme. Si un mouvement local amène une parcelle de gaz à un niveau plus élevé de façon assez rapide pour qu'aucun échange thermique n'ait lieu (par conduction ou rayonnement), celle-ci va se refroidir en suivant le gradient adiabatique, et sera donc plus froide et plus dense que ses environs immédiats. Cette parcelle aura donc tendance à retomber jusqu'à son niveau de départ, puisqu'un gaz plus froid est également plus dense toutes choses égales par ailleurs : par exemple, pour un gaz parfait, . On est donc en présence d'une atmosphère stable.
À l'inverse, si le profil thermique décroît plus fortement avec l'altitude que ce qu'indique le gradient adiabatique, cette parcelle de gaz sera certes refroidie si elle est soumise à un déplacement ascendant adiabatique, mais elle se retrouvera tout de même légèrement plus chaude que l'atmosphère environnante, et donc moins dense. Elle pourra donc continuer son mouvement ascendant jusqu'à ce qu'elle rencontre une zone stable où le profil thermique décroît moins vite que le gradient adiabatique. Une telle zone où des mouvements de convection à grande échelle peuvent se développer à partir d'une petite perturbation est dite instable. L'effet à long terme de ces mouvements de convection va conduire à un mélange qui homogénéisera le profil vertical de température jusqu'à retrouver une situation marginalement stable, c'est-à-dire avec un profil thermique suivant exactement le gradient adiabatique.
Les profils thermiques purement radiatifs tels que ceux modélisés ici ont tendance à voir leur pente croître en valeur absolue à mesure que la profondeur optique infrarouge croît en s'enfonçant dans l'atmosphère profonde. Sous couvert d'hypothèses raisonnables concernant la composition du gaz considéré parfait (pour ) et la croissance de selon le niveau de pression dans l'atmosphère, il est possible (mais hors-programme) de montrer que la pente du profil radiatif excède, en valeur absolue, le gradient adiabatique pour voisin de l'unité. Les régions atmosphériques situées en dessous () deviennent donc instables vis-à-vis de la convection qui s'y développe, et le profil thermique se met alors à suivre non plus la valeur donnée par le seul équilibre radiatif, mais le gradient adiabatique. On appelle cette couche atmosphérique troposphère. Les couches situées au-dessus () sont quant à elles stables vis-à-vis de la convection, et l'équilibre radiatif y est valable : on se trouve alors dans la stratosphère ou la mésosphère, selon l'existence ou non d'une inversion de température.
Notons qu'il existe quand même une troposphère dans les atmosphères des planètes telluriques trop peu opaques au rayonnement infrarouge thermique pour avoir (par exemple Mars, et dans une moindre mesure la Terre). En ce cas, l'instabilité de départ est causée par la discontinuité de température au niveau de la surface planétaire (voir ici), qui donne naissance à des mouvements de convection s'étendant jusqu'à une altitude équivalente à une échelle de hauteur environ.
Les profils thermiques les plus simples ne comportent qu'une troposphère surmontée d'une mésosphère, et le profil thermique y décroît toujours avec l'altitude. Mais il existe parfois au sein de la zone purement radiative une anomalie, une zone où la température croît avec l'altitude. Une telle zone est appelée stratosphère. Pour qu'une telle couche existe au sein d'une atmosphère, il faut qu'elle absorbe elle-même une partie du flux stellaire (dans le domaine visible, UV ou proche IR) et qu'elle soit relativement mauvaise émettrice en infrarouge thermique afin que l'énergie reçue par absorption du flux stellaire ne soit pas immédiatement perdue par rayonnement infrarouge thermique. Si l'on néglige les processus de diffusion lumineuse (ce qui est une hypothèse souvent vérifiée dans le domaine infrarouge thermique en l'absence de nuages, mais assez inexacte pour la lumière stellaire à plus courte longueur d'onde), le critère quantitatif pour l'existence d'une stratosphère est d'avoir une région verticale d'épaisseur optique en lumière stellaire et en infrarouge thermique tels que .
Dans le système solaire, la Terre possède une stratosphère due à la présence d'ozone, qui est un très bon absorbant de la lumière UV du Soleil. Comme, à l'altitude où cette absorption a lieu, l'atmosphère est froide et sèche, et que l'atmosphère terrestre est pauvre en , il y a peu d'absorption du rayonnement infrarouge, et donc aussi une faible émissivité infrarouge ( et étant les gaz à effet de serre principaux au sein des atmosphères telluriques). Les conditions d'existence d'une stratosphère sont donc réunies. En revanche, les atmosphères de Vénus et de Mars, constituées principalement de qui est un excellent émetteur infrarouge, ne possèdent pas de stratosphère. Dans le système solaire extérieur, on trouve également des stratosphères, dues à la présence de méthane () au sein de ces atmosphères qui absorbe bien dans l'infrarouge proche émis par le Soleil. Dans le cas de Titan, la stratosphère est due non seulement au méthane, mais aussi à l'absorption de la lumière solaire par les particules du brouillard photochimique qui l'entoure à haute altitude.
Au sommet de la mésosphère, vers un niveau de pression de , l'atmosphère devient trop peu dense pour être efficacement absorbante au rayonnement infrarouge et ainsi échanger de l'énergie de façon radiative. Le seul phénomène encore capable de transporter l'énergie devient alors la conduction thermique, obéissant à la loi de Fourier : où désigne la conductivité thermique du milieu et le flux de chaleur ainsi transporté. La structure thermique dans cette couche est alors dictée par la position des sources et des puits de chaleur :
Les positions respectives de ces puits et de ces sources causent un profil thermique croissant avec l'altitude, et pouvant atteindre des températures très élevées la journée car la conductivité thermique d'un tel milieu dilué est très faible, la chaleur peut donc y être piégée de façon très efficace. On nomme donc cette couche thermosphère. Sur Terre, la dissociation des molécules de par les UV solaires est une source de chaleur intense (ces molécules très fragiles vis-à-vis des rayonnements dissociants et/ou ionisants sont nombreuses dans l'atmosphère terrestre), si bien que les températures thermosphériques peuvent atteindre des valeurs très élevées, supérieures à . Pour les planètes géantes du système solaire, la source d'énergie est principalement due au chauffage par effet Joule dans l'ionosphère (friction des électrons libres). En revanche, dans les atmosphères telluriques riches en comme celles de Vénus et Mars, la dissociation des molécules est relativement difficile et le dioxyde de carbone est un radiateur efficace même à faible pression, ce qui entraîne des maxima de température diurne bien plus faible, pouvant même disparaître complètement pendant la nuit. On appelle alors parfois cette couche cryosphère lorsque ce phénomène se produit.
Voici quelques questions sur les définitions des grandeurs employées
Vous venez d'être embauché par un célèbre réalisateur Hollywoodien en tant que conseiller scientifique pour son prochain film de science-fiction. L'action se déroulera sur une lune tellurique nommée Pandore d'une planète géante appelée Polyphème en orbite autour de l'étoile . Toutes les données numériques pertinentes se trouvent ci-dessous.
Étoile (α Centauri) :
Polyphème :
Pandore :
On se propose à présent d'estimer l'opacité infrarouge de l'atmosphère de Pandore à l'aide d'un modèle simple. L'atmosphère est supposée parfaitement transparente en lumière visible et absorbe la totalité des rayonnements infrarouges thermiques. La température de l'atmosphère, supposée uniforme, sera notée .
Afin d'améliorer ce modèle, on ajoute une seconde couche atmosphérique partiellement opaque aux IR thermiques au-dessus de la première couche (qui reste complètement opaque à ces mêmes IR thermiques). On note la température de la couche supérieure et celle de la couche profonde . La couche 1 absorbe une fraction du rayonnement IR thermique. Ces deux couches sont toujours considérées parfaitement transparentes en lumière visible.
Compte tenu de la composition atmosphérique, s'attend-on à trouver une stratosphère sur Pandore ? Si oui, quelle serait l'espèce chimique responsable ?
Calculer le gradient adiabatique sec . Le gradient adiabatique humide sera-t-il inférieur ou supérieur en valeur absolue ?
Représenter l'allure du profil thermique moyen de Pandore. On considérera que la troposphère s'étend sur une échelle de hauteur , et on fera figurer l'échelle de hauteur, les différentes couches atmosphériques et les températures à leurs limites quand cela est possible.
Le but de ce mini-projet est de gagner une compréhension plus intuitive de l'équilibre radiatif au sein d'une atmosphère tellurique, et notamment du phénomène d'effet de serre, au moyen d'un modèle simplifié basé sur des hypothèses simplificatrices vues dans ce cours. Ceci permettra de vérifier par le calcul les résultats obtenus en manipulant ce modèle.
Le modèle est donc plan-parallèle, et traite le spectre électromagnétique en deux domaines distincts : le visible/UV/IR proche (shortwave, SW) et l'infrarouge thermique (longwave, LW ou IR). Ces deux domaines sont traités chacun de façon grise. L'atmosphère est constituée de deux couches considérées isothermes, pouvant absorber et émettre tout ou partie des rayonnements IR thermiques. Le rayonnement SW n'est quant à lui absorbé que par la surface, l'atmosphère étant donc parfaitement transparente. Les échanges autres que radiatifs sont négligés par ce modèle.
Voici la signification des différents champs à remplir :
Calculatrice
Modèle à deux couches
Dans cette partie, on fixe l'opacité de la couche n°2 à proximité de la surface à une valeur très élevée, et on ajuste seulement l'opacité de la couche n°1.
-------------FIN DU CHAPITRE -------------FIN DU CHAPITRE -------------FIN DU CHAPITRE -------------FIN DU CHAPITRE -------------
Les objectifs de ce chapitre sont, d'une part, d'acquérir quelques notions de base de dynamique des fluides et, d'autre part, d'explorer les différents régimes dynamiques qu'on peut observer sur les planètes de notre Système Solaire.
Les concepts abordés dans ce chapitre font appels à des notions de:
On connait très peu de choses sur les atmosphères des exoplanètes connues. Or cette question est fondamentale, tant pour les planètes géantes gazeuses dont l'atmosphère constitue la majeure partie de la planète, que pour les planètes à surface solide où le rôle de l'atmosphère dans leur habitabilité est primordial. Une atmosphère étant une enveloppe fluide en mouvement, son étude passe par sa dynamique et la compréhension de sa circulation.
Il existe deux méthodes complémentaires pour mieux comprendre la circulation atmosphérique : l'observation et la théorie. Notre connaissance de la circulation des atmosphères des exoplanètes est limitée car les observations sont à l'heure actuelle très difficiles. C'est pourquoi le recours à la théorie de la circulation atmosphérique, en se basant sur la mécanique des fluides, sert de socle à l'étude des circulations atmosphériques des exoplanètes. L'utilisation de modèles informatiques pour la simulation d'une atmosphère permet d'étendre cette théorie à des cas plus complexes tenant compte des nuages, aérosols, surface, etc ...
Les observations de la circulation atmosphèrique d'exoplanètes existent mais sont très rares. Par exemple, on a pu mesurer pour la planète HD 189733b son spectre d'émission et en tirer une carte de température. On a aussi pu mesurer les vents en haute altitude (d'au plus 10 000 km/h) de la planète HD 209458b par effet Doppler des lignes d'absorption du monoxyde de carbone présent dans son atmosphère. Toutefois, ces exemples restent rares si bien qu'en comparaison les observations des atmosphères planètaires de notre système solaire semblent exister à profusion. La connaissance de la dynamique atmosphérique des exoplanètes passe donc aussi et surtout par l'étude comparée des atmosphères du sytème solaire. Toutefois, notre système solaire ne présente pas toute la gamme des types d'exoplanètes connues, et par exemple l'étude d'une planète gazeuse géante chaude en rotation synchrone avec son étoile sur une orbite de 3 jours, cas assez exotique au regard du sytème solaire, passera également par la théorie et les moyens de simulation ...
Détecter et connaître les grandes caratéristiques physiques d'une exoplanète (orbite, taille, masse) est la première étape dans la compréhension d'une exoplanète. Connaître son atmosphère est la suite logique et la science actuelle en est à ses balbutiements dans ce domaine. L'étude des planètes du système solaire, tellement plus accessibles, et l'utilisation de modèles de climats s'avèrent nos meilleurs atouts pour comparer et extrapoler notre savoir à toutes les planètes en général.
On a vu que la notion de température d'équilibre ou effective requiert plusieurs hypothèses. L'une d'entre elles est que l'énergie rayonnée par unité de surface de la planète vers l'espace est homogène, c'est-à-dire qu'elle est la même en tous points, depuis l'équateur jusqu'aux pôles.
Or, la figure 1 nous donne un exemple pour la Terre qui montre la limite de cette hypothèse. On peut y voir que l'énergie absorbée par la Terre provenant du Soleil (et donc principalement dans la bande visible du spectre lumineux) dépend de la latitude, comme attendu d'après la figure 2. On y remarque surtout que l'énergie émise par la Terre vers l'espace, principalement de la bande infrarouge, dépend aussi de la latitude et est différente de la courbe d'énergie reçue. Ceci est un simple constat, issu d'observations, qu'il va s'agir d'interpréter.
En comparant les valeurs des puissances émises et reçues dans la figure 1, on constate qu'il y a une perte d'énergie vers les pôles et un gain d'énergie vers l'équateur. En vertu du principe de conservation de l'énergie, et en supposant que la planète est dans un état stationnaire, on ne peut que conclure que de l'énergie est transférée de l'équateur vers les pôles. Ce transfert d'énergie est causé par le mouvement global de l'atmosphère de la planète (et des océans s'ils sont présents).
On verra plus loin dans le cours que cette différence en énergie est également le moteur de la circulation atmosphérique.
Plus généralement, ce cours sur la dynamique atmosphérique va de pair avec celui sur la structure thermique des atmosphères. On y aborde ici les mouvements fluides à grande échelle, tandis que le cours sur la structure thermique aborde l'équilibre radiatif et les mouvements convectifs verticaux uniquement.
Les échanges d'énergie à l'échelle globale d'une planète sont fondamentaux pour comprendre la circulation à grande échelle. Historiquement, la première description et compréhension de ces mouvements a été faite sur Terre au 18ème siècle.
L'exemple le plus célèbre de la preuve d'une circulation atmosphérique à grande échelle provient de l'observation des vents permettant la navigation sur les océans terrestres. Les vents dominants soufflent d'Est en Ouest aux latitudes tropicales (environ 30°S à 30°N, et qu'on appelle alizé). Une explication satisfaisante de l'origine de ces vents fut donnée par George Hadley en 1735 : il existe une circulation à l'échelle de la planète qui transporte des masses d'air depuis l'équateur jusqu'aux tropiques en formant une cellule, qu'on appelle cellule de Hadley tel qu'on peut le voir sur la Figure 1.
L'origine des alizés est la suivante : quand une masse d'air tombe vers la surface aux tropiques, elle retourne à l'équateur. Ce faisant elle subit la force de Coriolis : elle est donc déviée vers sa droite dans l'hémisphère Nord, ou vers sa gauche dans l'hémisphère Sud, c'est-à dire vers l'Ouest dans les deux cas. Inversement, une masse d'air qui part de l'équateur vers les tropiques subira une déviation vers l'Est et formera un courant jet en altitude au niveau des tropiques et orienté vers l'Est.
L'origine des alizés est ainsi comprise, mais quid de l'origine de la cellule de Hadley en elle-même ? Elle provient de la différence d'énergie reçue par la Terre en fonction de la latitude. La hauteur d'échelle atmosphérique est plus grande à l'équateur (ou, autrement dit, l'air est plus dilaté), ce qui cause un gradient horizontal de pression entre l'équateur et les pôles, de plus en plus marqué à mesure que l'on s'élève du sol (voir Figure 2). Ce gradient de pression provoque un mouvement d'air en altitude, de l'équateur vers les tropiques. Le reste de la cellule de Hadley se met en place pour "fermer" la circulation d'air, par simple conservation de la masse.
Le mécanisme de la cellule de Hadley, connu depuis longtemps sur Terre, semble être universel dans les atmosphères des planètes, mais avec des caractéristiques différentes (extension en latitude, nombre de cellules) selon les propriétés physiques de l'atmosphère. Il est essentiel à la compréhension des mouvements à grande échelle d'une atmosphère et permet d'expliquer le transfert d'énergie de l'équateur vers les tropiques. Il est en réalité très complexe à comprendre dans les détails, et on trouve sur Terre d'autres mécanismes de transfert d'énergie au-delà des tropiques (voir Figure 1).
La circulation de Hadley n'est pas une cellule de convection. Son origine est une différence de chauffage en latitude, suivant l'horizontale, et non un gradient thermique vertical qui génère une force de flotabilité qui s'oppose à la gravité comme dans le cas de la convection. Ainsi, pour déterminer les dimensions de la cellule de Hadley il faut connaître la vitesse de rotation de la planète et le gradient horizontal de pression, des phénomènes qui ne rentrent pas en jeu dans les causes de la convection.
La rotation de la planète est la clé de la circulation atmosphérique, comme on vient de le voir sur Terre avec le rôle de la force de Coriolis dans la formation de la cellule de Hadley. Ainsi, il est fondamental de garder en tête les périodes de rotation pour les planètes du système solaire.
L'obliquité, qui est l'angle entre l'axe de rotation et la perpendiculaire au plan de l'orbite autour du Soleil, nous renseigne sur les saisons. Dans le cas de Vénus, l'obliquité proche de 180° signifie que la rotation est rétrograde : Vénus tourne de l'Est vers l'Ouest, dans le sens opposé aux autres planètes.
Corps | Période de révolution (jours) | Période de rotation | Obliquité |
---|---|---|---|
Vénus | 224.7 | 243.0 jours | 177.4° |
Terre | 365.2 | 23h56 | 23.4° |
Mars | 687.0 | 24h37 | 25.2° |
Jupiter | 4335.3 | 9h50 | 3.1° |
Saturne | 10757.7 | 14h14 | 26.7° |
Titan | 10757.7 | 15.95 jours | 26.7° |
Les informations que nous possédons sur ces planètes proviennent d'observations, de résultats de modèles et des équations régissant l'atmosphère.
Les méthodes d'étude des atmosphères planétaires comprennent des moyens différents tels que la télédétection (en utilisant des telescopes sur la Terre ou des vaisseaux spatiaux); des mesures directes in situ obtenues par sondes, ballons, atterrisseurs; des études de laboratoire; et des modèles de circulation atmosphèrique. Les mesures des mouvements atmosphériques sont accomplies principalement par les moyens suivants :
La Terre et Mars présentent des circulations à grande échelle très similaires typiques d'une planète en rotation rapide - les deux planètes ont des périodes de rotation similaires (Table). La différence principale est due à la présence des océans sur la Terre. En définitive, les deux atmosphères sont mises en mouvement par la différence d'énergie reçue sur la surface en fonction de la latitude.
La circulation atmosphérique à grande échelle se caractérise par des cellules de circulation entre différentes latitudes. Si la cellule de Hadley rejoint l'équateur aux latitudes moyennes à environ 30 ° Nord et Sud comme vu auparavant, les cellules de Ferrel ne fonctionnent pas du tout sur le même mécanisme que les cellules de Hadley. Elles mettent en jeu des ondes planétaires baroclines qui dominent les mécanismes de transfert. Les cellules de Ferrel se situent entre les latitudes 30 et 60 ° dans chaque hémisphère, dans lesquelles on trouve un courant-jet. On trouve également une cellule polaire de latitude 60° au pôle dans chaque hémisphère, délimitée par le vortex polaire.
Cette structure en cellule est un comportement moyen de l'atmosphère, qui varie avec les saisons et surtout avec les conditions météorologiques locales, dépendant de nombreaux phénomènes transitoires et chaotiques. Par exemple le centre de la cellule a tendance à se déplacer de part et d'autre de l'équateur lorsque le point subsolaire varie en latitude selon la période de l'année. De plus, il n'existe pas non plus de courant ascendant ou descendant continu et clairement mesurable dans les branches des cellules pour un observateur momentané en un endroit donné. Il s'agit là d'un comportement moyen qui peut être dominé par des effets locaux ou temporaires (cycle jour/nuit, topographie, perturbation météorologique, etc ...). La manifestation la plus tangible de ces cellules à tout un chacun demeure néanmoins la présence de courants jets.
Tout comme la Terre, Mars possède également une cellule de Hadley. Du fait de l'atmosphère plus ténue de Mars (la pression au sol y est de seulement 6 mbar en moyenne, par rapport à 1000 mbar sur la Terre), les gradients de pression relatifs des pôles à l'équateur sont plus grands et la cellule de Hadley y est plus étendue. Pendant l'équinoxe on trouve deux cellules de Hadley centrées sur l'équateur comme sur Terre, mais durant le solstice, on ne trouve plus qu'un seule grande cellule de Hadley des moyennes latitudes de l'hémisphère d'été vers les moyennes latitudes de l'hémisphére d'hiver. Il faut noter que la différence notable de position de l’ascendance au moment du solstice entre la Terre et Mars vient de la faible inertie de la surface martienne (contrairement aux océans terrestres), qui fait que le maximum de température se déplace complètement dans l’hémisphère d’été alors qu’il reste dans les tropiques sur Terre.
L'atmosphère de Vénus est caractérisée par une dynamique complexe: une super-rotation rétrograde domine dans la troposphère/basse mésosphère, tandis qu'une circulation solaire-antisolaire peut être observée dans la thermosphère. La super-rotation s'étend depuis la surface jusqu'au sommet des nuages avec des vents de seulement quelques mètres par seconde près de la surface et atteignant une valeur maximale de 100 au sommet des nuages (70 km), ce qui correspond à un periode de rotation de 4 jours terrestres. Les processus responsables du maintien de la super-rotation zonale dans la basse atmosphère et sa transition vers la circulation solaire-antisolaire dans la haute atmosphère sont encore méconnus.
En plus de la super-rotation zonale, une cellule de Hadley s'écoulant de l'équateur aux pôles avec des vitesses de moins de 10 a été observée au sommet des nuages. La circulation de Hadley ou méridienne sur Vénus joue un rôle important dans le transport d'air chaud vers les pôles et d'air froid vers l'équateur. Il s'agit d'une cellule dans chaque hémisphère et elle n'est pas limitée aux régions proches de l'équateur comme sur la Terre où la cellule de Hadley est limitée par la force de Coriolis, mais elle s’étend jusqu’aux pôles.
Il y a très peu mesures directes de la circulation atmosphérique de Titan, un satellite de Saturne, avec une période de rotation de 16 jours. Comme sur Vénus, une super-rotation est prévue dans la haute atmosphère; néanmoins, au contraire de Vénus, les saisons sur Titan jouent un rôle important dans la circulation générale. Les observations obtenues par la sonde Cassini-Huygens montrent la présence d'un seul jet très intense direct vers l'est avec des vitesses de 190 entre les latitudes 30-50 dans l'hémisphère nord. L'existence d'une circulation de Hadley de l'équateur au pôle a été également prédite. L’hémisphère Nord est l’hémisphère d’hiver au moment des observations de Cassini. Les cellules de circulation prédites vont d’un pôle (été, ascendance) à l’autre (hiver, subsidence) pendant la majeure partie de l’année, avec une transition autour de l’équinoxe au cours de laquelle l’ascendance se déplace pour changer d’hémisphère, la cellule se réduisant sur l’hémisphère de printemps pendant qu’une autre se développe sur l’hémisphère d’automne, puis prend sa place.
Sur les planètes telluriques les vitesses du vent sont mesurées par rapport à la surface de la planète. Comme les géantes gazeuses et les géantes glacées ne disposent pas d'une surface solide, les vents sont généralement mesurés par rapport à la période de rotation de leur champ magnétique. Chacune des quatre planètes géantes (Jupiter, Saturne, Uranus, et Neptune) montrent une circulation zonale intense formée par un système de courants-jets qui alternent leur circulation (Est-Ouest) avec la latitude.
Les planètes géantes gazeuses de notre Système solaire, Jupiter et Saturne, sont composées essentiellement de gaz légers, comme l'hydrogène et l'hélium. On observe aussi la formation de nuages constitués de méthane, d'ammoniac, et du sulfure d'hydrogène.
Au niveau des nuages d'ammoniac, Jupiter possède environ 6-8 courants-jets par hémisphère, avec une vitesse maximale à l'équateur où le jet dirigé vers l'est atteint une vitesse zonale de m s-1, et à 23° Nord les vents atteignent jusqu'à 180 m s-1. Saturne possède 4-6 courants-jets par hémisphère, avec un très fort et large jet équatorial direct vers l'est à une vitesse maximale de m s-1. Jupiter affiche également de nombreux tourbillons de courte durée et des tempêtes de longue durée comme la Grande Tache rouge.
Les planètes géantes glacées de notre Système solaire, Uranus et Neptune, sont constituées principalement d'eau, méthane ou d'ammoniac.
Les deux planètes glacées montrent un jet équatorial dirigé vers l'ouest - avec une vitesse de m s-1 et de m s-1 respectivement pour Uranus et Neptune - et deux jets intenses aux latitudes moyennes dirigés vers l'est avec une vitesse maximale de 200 m s-1. Comme sur Jupiter, forts tourbillons, convection et tempêtes ont été observés sur Neptune. Au contraire, sur Uranus aucune tempête n'a été observée.
Comme toute équation de la physique, les équations régissant la dynamique atmosphérique doivent s'exprimer dans un système de coordonnées et un référentiel choisis arbitrairement. Un tel système naturellement adapté à une sphère, et donc à une planète, sont les coordonnées sphériques , où est la distance au centre de la sphère, est la l'angle de la longitude, et est l'angle de la latitude.
On définit également un repère local pour tout point de l'espace de coordonnées , avec comme base le triplet où est dirigé vers l'Est, dirigé vers le Nord, et selon la verticale locale vers le haut. Le référentiel d'étude est ce référentiel local, lié à la rotation de la planète, il s'agit donc d'un référentiel tournant, donc non galiléen.
Pour résumer, on travaille dans deux référentiels différents, ce qui donne trois systèmes de coordonnées différents :
L'exercice suivant permet de se familiariser avec la manipulation mathématique des coordonnées et des repères. Les resultats serviront à établir l'équation fondamentale de la dynamique.
Montrer que dans le référentiel de la planète, un point de coordonnées sphériques a pour coordonnées cartésiennes
Exprimer , et dans le repère de la planète en fonction des angles et
Montrer qu'une vitesse dans le référentiel local au point de coordonnées sphériques s'exprime par .
Exprimer , et en fonction de , et .
Pour résoudre les équations régissant une atmosphère, on se place dans le référentiel local, lui-même dans un référentiel tournant avec la planète. Ce référentiel n'est pas inertiel, c'est-à-dire qu'il est en accélération par rapport à un référentiel inertiel. Afin de poser le principe fondamental de la dynamique, il est essentiel de tenir de compte de l'accélération apparente du référentiel d'étude, sous la forme de pseudo-forces.
Traditionnellement et pour des raisons pratiques on parle de force, ou pseudo-force, en multipliant l'accélération apparente par la masse de l'objet. On parlera ici plutôt de l'accélération d'une force apparente afin de s'affranchir du terme de masse, qui disparaitra dans les équations finales.
Les deux forces apparentes à considérer dans le cas d'une atmosphère sont la force centrifuge et la force de Coriolis.
L'expression mathématique des forces apparentes requiert l'emploi du produit vectoriel, défini ainsi:
Le produit vectoriel des vecteurs et s'écrit et correspond au vecteur orthogonal à la fois à et tel que le triplet soit de sens direct. et que le module du produit soit où est l'angle direct de vers . Ainsi, si et sont colinéaires, leur produit vectoriel sera le vecteur nul.
Le produit vectoriel s'écrit avec des coordonnées dans un système cartésien uniquement :
Cette relation n'est pas valable en coordonnées sphériques. Il faut faire la transformartion en coordonnées cartésiennes pour pouvoir utiliser cette relation.
La force centrifuge est la force apparente due au fait que le référentiel d'étude est en rotation, donc en accélération. En effet, l'orientation de la vitesse d'un point lié à la planète varie, mais pas son module. Ainsi, une particule au repos dans un référentiel galiléen aura une force apparente dans le référentiel de la planète.
On définit le vecteur rotation comme étant le vecteur orienté selon l'axe de rotation de la planète et de module avec la période de rotation de la planète. L'accélération de la force centrifuge s'exprime avec le vecteur position depuis l'axe de rotation où la force s'applique. Dans un système de coordonnées sphériques, on a , avec le vecteur position depuis le centre de la planète.
Une autre manière de l'exprimer est de dire qu'il s'agit d'une accélération perpendiculaire à l'axe de rotation, orientée vers l'extérieur et de valeur
La force centrifuge est regroupée avec la force de gravité, dont l'accélération vaut , l'indice faisant référence à la masse de la planète. On obtient une force dont l'accélération totale est . Cette force dérive d'un potentiel, qu'on appelle le géopotentiel, somme de l'action de la gravité et de la force centrifuge.
Démontrer par l'analyse ou par un schéma en 3D que correspond bien à une accélération perpendiculaire à l'axe de rotation, orientée vers l'extérieur et de valeur . Comment s'exprime cette accélération dans le référentiel local ?
Par la suite on fait l'approximation que pour l'accélération due au géopotentiel, on peut se contenter de considérer sa composante uniquement suivant l'axe pour l'additionner à la pesanteur, ceci afin de simplifier les équations du problème. Afin de s'assurer que cette approximation reste acceptable, quelle est l'erreur maximale faite dans le cas de la Terre si on ne considère pas les autres composantes de l'accélération due au géopotentiel ? Quelle est cette erreur relativement à l'accélération de pesanteur ? Et pour Jupiter ?
Dans un référentiel en rotation, une autre force apparente est à prendre en compte lors d'un déplacement. Il s'agit de la force de Coriolis, qui tient compte du fait que le déplacement d'une particule génère une accélération apparente supplémentaire. Par exemple, le mouvement rectiligne d'une particule est apparement dévié pour un observateur situé dans un référentiel tournant.
L'accélération de la force de Coriolis s'exprime avec le vecteur vitesse de la parcelle d'air considérée.
La vitesse d'une parcelle d'air dans l'atmosphère est généralement orientée parallèlement à la surface locale, c'est-à-dire que sa composante radiale (c'est-à-dire sa composante verticale locale) est en générale petit par rapport à au moins une des deux autres. En négligeant la composante radiale, on constate les choses suivantes :
Ceci explique pourquoi certaines structures atmosphériques, tels les ouragans ou les anticyclones, tournent toujours dans le sens des aiguilles d'une montre dans l'hémisphère Nord, et dans le sens contraire dans l'autre hémisphère.
Démontrer par l'analyse ou par un schéma en 3D les 3 points ci-dessus.
Les écoulements atmosphériques peuvent être décrits en utilisant deux points de vue classiques, appelés eulérien ou lagrangien :
L'écoulement est suivi par un observateur depuis une position fixe. C'est le cas, par exemple, d'un atterrisseur sur Mars fixé au sol qui mesure la vitesse du vent, la température, ou la pression. Cette description est souvent préférée car elle est la plus pratique.
Dans ce cas, les particules fluides sont suivies le long de leurs trajectoires. C'est la description la plus intuitive.
Les descriptions lagrangienne et eulérienne sont liées à travers la dérivée particulaire, encore appelée dérivée totale, et qui s'écrit . Soit une grandeur physique vectorielle de l'écoulement, dépendant du point d'observation et du temps . La variation de la grandeur s'écrit :
Où est la vitesse du fluide avec composants: , , .
Dans la cas où la grandeur est un champ scalaire , la relation est la même.
La dérivée particulaire décrit la variation avec le temps en suivant la particule en mouvement (point de vue lagrangien), en revanche décrit la variation locale avec le temps en un point d'observation fixé (point de vue eulérien) .
Le mouvement d'une particule dans un fluide est décrit par la deuxième loi de Newton (conservation de la quantité de mouvement) qui lorsqu'elle est appliquée à la mécanique des fluides donne l'équation de Navier-Stokes. Dans un système en rotation l'équation du mouvement d'une parcelle de fluide est:
avec est la dérivée particulaire qui s'écrit , la vitesse du fluide, la somme des forces s'appliquant sur la parcelle et la densité du fluide.
Soit en détaillant les forces :
avec l'accélération de la force de Coriolis, les forces dues au gradient de pression, l'accélération du géopotentiel, et qui désigne l'accélération dues à la viscosité. Ce qui donne :
(1)
où est le vecteur de rotation de la planète est le gradient de pression.
Comment s'exprime l'accélération de la force de Coriolis dans le repère local ?
On a . Or il a été vu en exercice les expressions des dérivées temporelles des vecteurs , et . Ceci nous permet d'établir les équations de Navier-Stokes dans le référentiel local, en notant que et :
Ce système d'équations décrit tous les types de mouvements atmosphériques à toutes les échelles. Ces équations sont compliquées à résoudre, mais dans bien des cas utiliser une approximation est suffisante pour modéliser de nombreux phénomènes atmosphériques dynamiques.
Déduire les équations de Navier-Stokes en coordonnées sphériques données ci-dessus à partir de l'équation fondamentale de la dynamique (Equation 1).
Les équations de la dynamique sont très compliquées car elles forment un système non linéaire. Ceci signifie que la somme de deux solutions n'est pas forcément solution du problème, ce qui rend la résolution de ces équations très ardue, et à ce jour encore source de recherches. Cependant, en fonction des phénomènes étudiés et des caractéristiques de l'atmosphère planétaire, certains termes de ces équations peuvent en dominer d'autres. Pour estimer les différents termes dans les équations, on utilise la méthode de l'analyse d'échelle. Les ordres de grandeur des différents termes en jeu dans les équations fondamentales de la dynamique seront très différents selon l'échelle des écoulements que l'on souhaite étudier. Dans le tableau ci-dessous on compare les termes dominants sur les planètes à rotation rapide (la Terre) avec ceux sur les planètes à rotation lente (Vénus):
Terre | 10-5 | 10-5 | 10-8 | 10-3 | 10-6 | 10-12 |
Vénus | 10-3 | 10-5 | 10-5 | 10-5 | 10-7 | 10-12 |
avec le rayon de la planète. On a , où est l'altitude depuis la surface.
On peut alors appliquer les approximations suivantes:
On obtient alors les équations primitives de la météorologie :
À ce système d'équations on ajoute l'équation des gaz parfaits:
Avec la constante universelle des gaz parfaits et la masse molaire du gaz qui constitue l'atmosphère, et dépend donc de sa composition. Pour l'air terrestre, on a
ainsi que l'équation de conservation de la masse:
Enfin, le premier principe de la thermodynamique:
Avec le forçage diabatique et la température potentielle : , où , la chaleur spécifique à pression constante et une pression de référence.
On obtient ainsi 6 équations avec 6 inconnues ().
Ce système d'équations primitives est le plus complet utilisé pour l'étude de la circulation générale de l'atmosphère. C'est notamment celui utilisé par les modèles de circulation générale.
Les équilibres géostrophique et cyclostrophique sont deux approximations des équations primitives. Ils sont purement diagnostiques : ils ne contiennent pas de dérivées dans le temps, d'où l'impossibilité de faire des prédictions. Néanmoins, ils sont des outils puissants pour décrire différents écoulements observés dans les planètes.
L'approximation géostrophique est un développement des équations primitives utilisée aux moyennes latitudes sur les planètes à rotation rapide (Terre, Mars). On suppose l'équilibre entre la force de Coriolis et la force due au gradient horizontal de pression. La force centrifuge est négligée.
D'après ces équations, lorsque cet équilibre est valide, la vitesse du vent est directement proportionnelle au gradient horizontal de pression. Notez que l'équilibre géostrophique cesse d'être valide autour des latitudes équatoriales.
En combinant les deux composantes de la vitesse, on peut introduir le vent géostrophique comme :
La circulation générale des planètes à rotation lente (Vénus, Titan), aussi bien que les vortex et les tourbillons sur toutes les planètes, peut être approximée par l'équilibre cyclostrophique. Cela suppose l'égalité entre la composante dirigée vers l'équateur de la force centrifuge et le gradient méridional de la pression. La force de Coriolis est négligée.
L'équation du vent cyclostrophique peut alors être écrite comme :
Montrer à partir de l'équation que l'équation du vent cyclostrophique peut être écrite comme :
où : J kg-1 K-1, et est la coordonnée de pression logarithmique, avec la pression au niveau de référence.
Le nombre de Rossby est un nombre sans dimension qui permet de caractériser les mouvements atmosphériques. Il est défini par:
où est une vitesse caractéristique du système, est la vitesse angulaire de rotation de la planète, et est une longueur caractéristique du système.
Quelle est la dimension du nombre de Rossby ?
Une valeur de nombre de Rossby très supérieure à l'unité indique que la force de Coriolis due à la rotation de la planète est négligeable par rapport à la force d'inertie, dans ce cas on parle d'équilibre cyclostrophique. Dans le cas contraire d'un nombre de Rossby inférieur à l'unité, l'équilibre est dit géostrophique.
Equilibre géostrophique | [Terre, Mars] | |
Equilibre cyclostrophique | [Vénus, Titan, Ouragans] |
Les valeurs du nombre de Rossby pour différents systèmes sont comparées sur le tableau ci-dessous :
Vénus | ∼103 |
Terre | ∼1 |
Mars | |
Titan | |
Tourbillons de poussière | ∼102-103 |
Tornades | ∼103 |
Ouragans |
La Terre et Mars présentent une circulation atmosphérique aux grandes échelles très similaire et typique des planètes à rotation rapide : les deux planètes ont en fait une période de rotation similaire (Voir Tableau). La principale différence entre eux vient de:
Un régime géostrophique des vents zonaux domine la circulation dans les deux planètes en dehors des latitudes tropicales. Aux latitudes moyennes, à la fois sur la Terre et sur Mars, la circulation est caractérisée par deux courant-jets, un dans chaque hémisphère.
Sur la Terre, ces jets sont des vents zonaux qui circulent de l'ouest vers l'est et leur vitesse augmente avec l'altitude jusqu'à la tropopause. Au-dessus de la tropopause ces jets affaiblissent, et puis ils augmentent encore avec l'altitude au-dessus de km jusqu'à km, où ils atteignent une vitesse de 60 m s-1.
Sur Mars, en raison de l'atmosphère très ténue et de l'absence des océans, l'atmosphère réagit presque instantanément au chauffage solaire. C'est aussi la raison pour laquelle les courant-jets dépendent des variations saisonnières. Au solstice d'hiver dans l'hémisphère nord, le courant-jet est centré entre de latitude atteignant une vitesse maximale de 40 m s-1 à 5 km. Il augmente encore à 35 km, où il atteint une vitesse de 110 m s-1. A cette même époque, le courant-jet dans l'hémisphère sud est beaucoup plus faible.
Il faut noter que ces configurations des vents sont une moyenne temporelle et spatiale et ils sont vus rarement sur des journées individuelles. Les configurations des vents d'un jour à l'autre dévient considérablement de cette circulation globale.
Vénus et Titan sont deux planètes à rotation lente, caractérisées respectivement par une période de rotation de 244 jours (Vénus) et 16 jours (Titan). Une description détaillée de leur circulation peut être trouvée ici. Les deux planètes sont caractérisées par des forts vents zonaux dans l'ensemble de l'atmosphère, une caractéristique appelée super-rotation. Sur Vénus, la super-rotation atteint une vitesse supérieure à 100 m s-1 au sommet des nuages (vers 70 km d'altitude), correspondant à une période de rotation de 4 jours terrestres (60 fois plus rapide que Vénus elle-même). Différentes études ont montré que sur les planètes qui tournent lentement, comme Vénus et Titan, les forts vents zonaux au sommet des nuages peuvent être décrits par l’équilibre cyclostrophique. Ce qui donne une possibilité de reconstruire le vent zonal si le champ de température est connu.
D'autres systèmes cyclostrophiques à petite échelle sont les ouragans et les tornades (Figure 1); ils sont caractérisés par un centre de basse pression et de forts vents. En raison d'un nombre de Rossby élevé (Tableau 1), la force de Coriolis peut être négligée pour les tornades et les ouragans, et on suppose l'équilibre entre la force centrifuge et le gradient de pression.
Selon le glossaire de météorologie (AMS 2000), une tornade est une colonne d'air tournant violemment, en contact avec le sol et la base des nuages, et souvent (mais pas toujours) visible comme un nuage en forme d'entonnoir (Figure 1). La plupart des tornades ont des vitesses de vent entre 18 m s-1 et 135 m s-1. Son vortex a un diamètre typique de quelques centaines de mètres et tourne, en général, dans le sens contraire des aiguilles d'une montre dans l'hémisphère Nord. Les tornades se produisent sur tous les continents, mais sont plus fréquentes dans les États-Unis, où le nombre moyen de tornades enregistrées est d'environ 1000 par an, avec la majorité d'entre eux sur les plaines centrales et dans les états du sud. Les tornades sont associées à des orages violents et sont alimentés par l'afflux d'air chaud et humide. En général, ils sont le résultat de l'instabilité produit par la différence de température et d'humidité entre la surface, où l'air est plus chaud, et les niveaux supérieurs de l'orage, où l'air est plus froid.
Ouragan est le nom utilisé pour indiquer les cyclones tropicaux qui se produisent dans l'Atlantique ou le Pacifique Est. Les ouragans sont marqués par une région centrale d'air descendants, l'oeil, enfermé par des orages forts associés à des vents et des pluies intenses. Comme pour le cas des tornades, l'énergie des ouragans est fournie principalement par libération de chaleur latente dans l'air humide. Les ouragans sur la Terre se forment dans les régions tropicales au-dessus des océans chauds, et ils s'affaiblissent lorsqu'ils arrivent sur terre, où la source d'énergie disparaît. Dans l'oeil d'un ouragan, le nombre de Rossby locale est toujours et peut arriver jusqu'à 100. Dans ce cas, l'équilibre devient cyclostrophique.
Les cyclones tropicaux sur la Terre et le vortex polaire de Vénus présentent des similitudes morphologiques et dynamiques, comme on peut voir dans la Figure 2. Le vortex de Vénus et les ouragans sont caractérisés par différentes échelles horizontales et durée de vie: le vortex de Vénus à un diamètre de 12000 km et il semble être permanent; les plus grands cyclones tropicaux observés sur la Terre ont un rayon de moins de 1000 km et durent environ une à deux semaines dans leur phase de maturité. La source d'énergie est aussi différente pour le vortex sur Vénus et les ouragans terrestres: la source d'énergie pour les ouragans est la libération de chaleur latente; le vortex polaire sur Vénus reçoit un apport d'énergie par le dépôt du rayonnement solaire au niveau des nuages et par l'émission thermique dans la basse atmosphère. Malgré leurs différences, les circulations du vortex de Vénus et des ouragans est très similaires: à partir du leur comparaison une meilleure compréhension de la dynamique de Vénus peut être atteinte.
Les tourbillons de poussière, présents à la fois sur la Terre et Mars, sont caractérisés par des vitesses de vent de rotation élevées, des champs électrostatiques importants et sont rendus visibles par la présence de poussière et de sable soulevé. Ils sont distincts des tornades car les tornades sont associées à des orages tandis que les tourbillons de poussière se forment sous un ciel clair.
Sur la Terre (Figure 3), l'étude de tourbillons de poussière est fondamentale pour comprendre leur rôle dans la convection et l'érosion des zones arides. Les tourbillons de poussière se produisent généralement en été dans les régions désertiques plus chaudes. Ils sont des événements transitoires et la plupart ne durent que quelques minutes.
Sur Mars (Figure 4), les tourbillons de poussière peuvent avoir un effet important sur le cycle global de poussière. Les tourbillons de poussière sur Mars ont d'abord été identifiés dans les images prises par l'orbiteur Viking comme des petits nuages avec des longues ombres coniques. En plus, des traces laissées au sol par des tourbillons de poussière ont été détectées dans des images de la Mars Orbiter Camera et ont été également observées au sol par des atterrisseurs. Les tourbillons de poussière martiens et terrestres semblent avoir une morphologie similaire. Cependant, les tourbillons de poussière martiens sont un ordre de grandeur plus grand que ceux terrestres, atteignant souvent quelques kilomètres d'altitude et des centaines de mètres de diamètre avec des bases étroites et des larges sommets.
Pour caractériser les différents types des instabilités atmosphériques on utilise le nombre de Richardson, un nombre sans dimension défini par:
Où est le cisaillement vertical du vent.
est nommée fréquence de Brunt-Väisälä définie comme la différence entre le gradient vertical de température et le gradient adiabatique :
est la fréquence d'oscillation d'une particule soumise à un déplacement vertical. Pour l'atmosphère est instable et une particule déplacée de son état initial s'éloignera irréversiblement. Si , la stabilité est "neutre", la particule déplacée demeura à sa nouvelle altitude. Enfin, pour se produit une oscillation de la particule autour de son état initial.
Dans ce cas, correspondant à un valeur négative de la fréquence de Brunt-Väisälä , la couche atmosphérique est instable et la turbulence est soutenue par la convection.
Parfois on observe de la turbulence dans des couches atmosphèriques thermiquement stables. Cette turbulence, dite de Kelvin-Helmholtz, est crée dans des régions où il y a du cisaillement du vent. Un valeur positive de nombre de Richardson au-dessous d'une valeur critique est une condition nécessaire afin que l'instabilité de Kelvin-Helmholtz se puisse produire. L'écoulement de Kelvin-Helmholtz est le résultat du cisaillement de vitesse entre deux fluides glissant l'un par rapport à l'autre.
Condition suffisante pour une écoulement stable.
Les ondes atmosphériques sont des perturbations des champs atmosphériques qui se propagent dans l'espace et/ou le temps. C'est un mécanisme important dans la dynamique des atmosphères car les ondes permettent de transporter des perturbations, transporter de l'énergie et de la quantité de mouvement d'une région à une autre.
On peut représenter de manière simplifiée une onde atmosphérique par une fonction sinusoïdale, en fonction d'une dimension spatiale de coordonnée et d'une dimension temporelle de coordonnée :
Où est l'amplitude de l'onde; est le nombre d'onde; est la longueur d'onde (en mètres); est la pulsation; est la période (en secondes). est la phase de l'onde, c'est-à-dire la valeur de la perturbation lorsque et . La longueur d'onde est définie comme étant la distance séparant deux crêtes consécutives d'une onde. Si (en mètres par seconde) est la vitesse de propagation de l'onde, on définit la fréquence (en hertz) par : .
Le champ physique représenté par est une variable atmosphérique. Il peut s'agir de la température, pression, le vent, etc ... La dimension de l'amplitude est donc la même que celle de la variable représentée par la perturbation .
Une manière plus compacte et efficace pour représenter une onde est la notation exponentielle. On écrit la perturbation sous sa forme complexe de la manière suivante :
Avec le nombre imaginaire . La perturbation réelle est définie comme étant la partie réelle de sa forme complexe :
En utilisant la relation trigonométrique bien connue , on obtient que l'amplitude complexe vaut :
L'amplitude complexe contient ainsi l'information à la fois sur l'amplitude et la phase de l'onde. Cette notation est très pratique car elle permet notamment de dériver ou d'intégrer une onde par rapport à l'espace ou au temps. Par exemple :
Ainsi, dériver par rapport à la coordonnée spatiale revient à multiplier l'onde complexe par . De même, une dérivation temporelle revient à multiplier par .
De la même manière, on peut montrer que trouver une primitive de l'onde complexe revient à diviser par , donc à multiplier par . De même, multiplier par permet de trouver une primitive par rapport à .
À partir des équations primitives, il est possible de trouver les ondes susceptibles de se propager dans l'atmopshère en utilisant la méthode des perturbations. Il s'agit d'écrire chaque champ (par exemple avec la pression ) comme étant la somme d'une valeur fixe solution des équations et d'une petite perturbation : . Ceci permet de linéariser les équations primitives en obtenant une équation pour les petites pertubations. Ces pertubations correspondent à des ondes que l'on peut ainsi étudier au moyen d'un cadre formel.
L'onde la plus évidente est l'onde sonore dont le calcul va être détaillé ci-dessous. On définit les petites perturbations comme étant des ondes se propageant horizontalement et verticalement :
avec et les nombres d'ondes horizontaux et verticaux respectivement. On fait l'approximation que la rotation et la gravité sont négligeables dans le cas qui nous intéresse. De plus, on suppose un fluide au repos, où la dérivée lagrangienne est égale à la dérivée eulérienne. Ainsi dans les équations primitives la force de Coriolis s'annule, tout comme la force centrifuge. En posant avec la vitesse du son, les équations du mouvement horizontal et de l'équation de continuité s'écrivent ainsi :
Soit, en utilisant les propriétés de la notation exponentielle :
Par identification, on obtient :
qui nous donne la relation entre longueur d'onde et période d'une onde sonore.
Le traitement des ondes atmosphériques est un sujet complexe; dans la page qui suit, nous allons donner un aperçu général des principaux types d'ondes en comparant différentes planètes, mais sans entrer dans le détail.
Une analyse plus détaillée des phénomènes ondulatoires peut être trouvé dans la liste suivante des livres:
Les ondes atmosphériques peuvent se manifester de diverses manières: comme oscillations de la température, de la densité et de la vitesse du vent, ou à travers des structures régulières de nuages. Ils peuvent être classés sur la base de facteurs différents: (1) mécanismes de restauration; (2) échelles de temps et d'espace; (3) ondes stationnaires ou qui se déplacent.
Le mécanisme de restauration des ondes de gravité est la poussée d'Archimède. Elles sont créées: (1) par la topographie (ondes orographiques), lors du passage d'une masse d'air au-dessus d'un relief montagneux; ou (2) par des instabilités dues à la présence d'orages ou de fronts d'air. Ils sont très communs dans les mésosphères des planètes telluriques. Sur Terre, les ondes de gravité révèlent souvent leur présence à travers des formations nuageuses, comme dans le cas des ondes orographiques. À des altitudes plus élevées, des images des nuages polaires mésosphériques, aussi connus sous le nom de nuages noctulescents, montrent souvent des structures ondulatoires probablement dues à la propagation vers le haut des ondes de gravité. Dans l'atmosphère martienne, des ondes de gravité ont été observées dans des images de nuages mésosphériques de glace de CO2 obtenues par la High Resolution Stereo Camera à bord de la mission européenne Mars Express et elles sont produites par la topographie (montagnes et cratères). Plusieurs missions spatiales ont détecté des ondes de gravité dans l'atmosphère de Vénus à la fois comme des oscillations sur le champ de température et comme des structures sur la couche de nuages. Elles sont probablement causées par des mouvements convectifs ou par des instabilité de de Kelvin-Helmholtz. Dans la haute atmosphère des planètes géantes (Jupiter, Uranus, et Neptune), des ondes de gravité ont été détectées par des oscillations dans les profils verticaux de température. Elles sont probablement la manifestation d'ondes de gravité générées près de la tropopause se propageant vers le haut.
Ces ondes sont un type d'ondes de gravité produites par l'instabilité de Kelvin-Helmholtz (Figure 2).
Les ondes de Rossby ou ondes planétaires sont caracterisées par une grande longueur d'onde. Elles sont dues à la variation de la force de Coriolis selon la latitude, quand une masse d'air se déplace à latitudes plus ou moins élevées. Sur la Terre, les ondes de Rossby sont des phénomènes permanents dans les latitudes moyennes et sous-polaires. Souvent, des cyclones et des anticyclones se forment dans les crêtes de l'onde. De la même façon, des ondes de Rossby se forment dans l'atmosphère vers les régions polaires de Mars. Dans les planètes géantes Jupiter et Saturne, il y a une grande variété d'ondes atmosphériques à l'extérieur de la région équatoriale qu'on suppose être des ondes de Rossby. Sur Saturne, un système ondulatoire hexagonal existe autour du pôle Nord, qui pourrait être une onde de Rossby. On peut voir divers exemples d'ondes de Rossby sur différentes planètes dans la figure ci-contre.
Les ondes de marées thermiques sont des ondes d'échelle planétaire excitées par les variations de l'insolation du sol dues au cycle jour-nuit. Ces ondes se manifestent sur le champ de pression et sur les composants du vent et elles évoluent avec le temps solaire local. Les ondes de marée thermique dans l'atmosphère de Mars ont une amplitude beaucoup plus élevée que sur la Terre car le forçage thermique sur Mars est très fort à cause de l'absorption dans le proche infra-rouge du atmosphérique, l'absorption du rayonnement infra-rouge émis par la surface, la présence de la poussière dans l'atmosphère et le fait que l'atmosphère y soit plus ténue. L'effet des marées thermiques sur la circulation zonale et méridienne moyenne est donc très important dans l'atmosphère martienne. Sur Vénus, l'absorption du rayonnement solaire se produit essentiellement dans la couche de nuages et des études ont proposé que les marées solaires jouent un rôle important dans la super-rotation de Vénus à ce niveau.
Dans la figure suivante, on voit que l'énergie absorbée dans le visible dépend de la latitude (courbe bleue). Sachant que la puissance émise par le Soleil et reçue par la Terre est de 1361 W/m2 au total, quelle est l'équation régissant la relation entre la puissance reçue sur Terre et la latitude ?
Quelle serait l'allure de cette figure pour une planète ne possédant pas d'atmosphère ?
Dans la figure suivante, on a fait figurer le vent à deux altitudes différentes (300 et 925 Pa). Comment interprétez-vous l'orientation et l'intensité des vents ? Quel lien faites-vous avec les zones où les nuages sont absents ou présents ? Les zones arides et boisées sur Terre ?
Quel est le lien entre variation particulaire et la variation totale d'un grandeur ? Comment interpréter le cas où le fluide est au repos (vitesse du fluide nulle) ?
Un ami vous affirme que le sens de rotation d'un vortex créé par de l'eau s'écoulant d'un lavabo dépend de l'hémisphère dans lequel on se trouve. Il en veut pour preuve que ce sens est toujours le même dans sa salle de bains et que ses nombreuses expériences de voyage de par le monde ne permettent pas de mettre ce fait en doute. Qu'en pensez-vous ?
Vous vous trouvez sur un manège tournant et souhaitez lancer un balle à votre ami situé de l'autre côté du manège. Si le manège tourne dans le sens des aiguilles d'une montre, Devez-vous lancer la balle à droite, à gauche, ou dans la direction de votre ami pour qu'elle lui arrive directement dans les bras ?
Vue de haut : vous êtes un bonhomme bleu avec un ami orange sur un manège gris.
Contrairement à vous, votre ami a eu mal au coeur et est descendu du manège pour s'asseoir sur un banc. Vous souhaitez néanmoins toujours lui lancer une balle, ce que vous tentez au moment où vous êtes le plus proche de lui. Devez-vous lui lancer la balle à droite, à gauche ou dans sa direction.
Vue de haut : vous êtes toujours un bonhomme bleu sur un manège gris mais votre ami orange n'est plus sur le manège. .
Il est possible de recréer une gravité artificielle dans un vaisseau spatial en le mettant en rotation autour d'un axe central. Les astronautes peuvent ainsi profiter d'une gravité telle que ressentie à la surface de la Terre dans un anneau qui tourne autour de son centre. Afin qu'ils ne ressentent pas de gêne lorsqu'ils se déplacent dans l'anneau, quelle doit être le diamètre minimal de l'anneau ?
Ce mini-projet propose de visualiser des données calculées par un GCM afin d'aborder différents points vus dans le cours sur la dynamique atmosphérique.
Il est préférable d'avoir lu le cours sur le GCM et effectué le mini-projet de ce même cours pour prendre en main l'outil de visualisation et avoir une compréhension des données à manipuler.
Afficher la structure zonale de la température pour une planète similaire à la Terre, en fonction de la latitude et de la longitude. Qu'observez-vous ? Comment interpréter cette structure ?
Faire de même en faisant varier la vitesse de rotation de la planète. Qu'observez vous ? Comment l'interpréter ?
Faire de même en affichant le vent méridien, la circulation méridienne, et la super-rotation. Conclure.
Faire le lien avec les planètes connues du système solaire: Vénus, la Terre, Mars, Jupiter et Saturne.
Afficher et identifier diverses ondes grâce à la pression de surface, les vents, la température
On se propose de vérifier la relation suivante de Lewis et consorts sur le lien entre rotation et taille sur la circulation grande échelle, etc ...
-------------FIN DU CHAPITRE -------------FIN DU CHAPITRE -------------FIN DU CHAPITRE -------------FIN DU CHAPITRE -------------
Plusieurs milliers d’exoplanètes, en orbite autour d'étoiles autres que le Soleil, ont été détectées à ce jour. La plupart d'entre elles ont été découvertes en étudiant les variations de luminosité liées à leur passage devant leur étoile (méthode des transits) et grâce à l'observation du mouvement de l'étoile autour de laquelle les planètes sont en orbite (méthode des vitesses radiales). Ces deux méthodes reposent sur l'exploitation d'informations liées à l'interaction gravitationnelle entre la planète et son étoile. Ce cours est essentiellement un cours de mécanique. Il a pour objectif de modéliser le mouvement d'une planète autour de son étoile et d'en identifier les paramètres caractéristiques. La partie Décrire présente les phénomènes de façon aussi simple que possible et ne requiert pas de compétences en mathématique. La partie Comprendre s'attache à démontrer les résultats énoncés dans la partie précédente. Le problème à deux corps (interaction entre une planète et son étoile) est au coeur de ce module. Le problème à plus de deux corps, qui est mathématiquement nettement plus compliqué, sera uniquement évoqué et permettra d'introduire les points dits de Lagrange. La migration des planètes sera également évoquée. Les illustrations seront systématiquement empruntées au monde des exoplanètes.
Les exoplanètes, comme les planètes et autres corps de notre système solaire sont soumises à la loi de la gravitation ou attraction universelle. Cette loi a été historiquement formalisée par Isaac Newton en 1687. Elle permet de modéliser, de façon générale, l'attraction entre des corps ayant une masse et, par conséquent, le mouvement des corps célestes.
Bien avant Newton, le mouvement des planètes dans le ciel a retenu l'attention de nombreux scientifiques (Aristote, Galilée, Copernic,...) qui ont proposé des modèles pour expliquer ou, au moins, modéliser leurs observations. Les travaux de Copernic font date avec sa proposition en 1543 d'un modèle héliocentrique qui présente, entre autres, le grand avantage de simplifier les calculs par rapport au modèle géocentrique. Cependant, la cause du mouvement n'est pas identifiée et les orbites des planètes restent des cercles.
Johannes Kepler établit, à partir des observations minutieuses du mouvement de Vénus, Mars, Jupiter et Saturne faites par Tycho Brahe, le fait que les orbites des planètes de notre système solaire ne sont pas des cercles mais des ellipses dont le Soleil est un foyer. Il publie, en 1609 et 1619, des lois empiriques qui prendront le nom de Lois de Kepler. Ces lois sont particulièrement intéressantes et utiles parce qu'elles établissent des relations entre les différents paramètres des orbites. En outre, il faut noter que Kepler a certainement ouvert la voie à Newton en affirmant que "deux corps voisins et hors de la sphère d'attraction d'un troisième corps s'attireraient en raison directe de leur masse."
Newton a finalement modélisé la gravitation est ainsi permis une explication des phénomènes observés ainsi que la possibilité de prédire le mouvement des planètes en orbite autour de leur soleil. Les lois de la gravitation de Newton permettront d'ailleurs d'expliquer les résultats empiriques énoncés dans les lois de Kepler. Notons que la mécanique classique issue des lois de Newton n'est elle-même qu'une approximation, qui ne rend plus bien compte des mouvements des corps aux vitesses très élevées ou au voisinage immédiat des masses. Ainsi c'est la relativité générale d'Einstein qui a permis d'expliquer le mouvement du périhélie de Mercure ce qui n'était pas possible à partir des lois de Newton.
Les trois lois de Kepler ont été établies à partir de l'observation du mouvement des planètes au sein de notre système solaire. Elles ne reposent sur aucune modélisation mécanique de l'interaction entre la planète et le soleil. Néanmoins, ces lois seront confirmées, a posteriori, par la théorie de la gravitation universelle. Dans cette partie, ces lois sont énoncées sans démonstration (les démonstrations seront faites plus loin) et illustrées.
Loi des orbites : Les planètes autour de leur soleil ont un mouvement périodique dans un plan. Plus précisément, ces orbites sont des ellipses dont le soleil est un foyer.
L'ellipse est caractérisée géométriquement dans son plan par deux paramètres : La longueur de son demi-grand axe a ( 2*a mesure la taille de l'ellipse dans sa plus grande dimension) et son excentricité ( est un nombre sans dimension ) qui décrit le degré d'allongement de l'ellipse. Le cercle est le cas particulier qui correspond à . Le cas particulier correspond à une parabole (qui ne convient pas pour une orbite fermée mais décrit la chute des corps ). La trajectoire de la planète sur son orbite est décrite en coordonnées polaires par la relation . Sous cette forme, la distance r est mesurée par rapport au centre de l'étoile (qui est sur l'un des foyers) et l'angle polaire est appelé l'anomalie vraie.
La distance planète-étoile est souvent exprimée en Unités Astronomiques (UA ou AU en anglais). 1 UA correspond approximativement à la distance moyenne entre la Terre et le Soleil, soit 149 597 870 000 m.
La valeur correspond à la direction du point de l'ellipse le plus proche de l'étoile. Cette position qui correspond à la distance planète-étoile minimale est appelée périastre (périhélie pour la Terre qui tourne autour du Soleil). L' apoastre (aphélie dans le cas de la Terre ) est la position qui correspond à la distance planète-étoile maximale.
Exprimer la distance étoile-planète lorsque la planète est à son apoastre. Même question lorsqu'elle est située à son périastre .
Calculer, en utilisant les données du catalogue des exoplanètes, la valeur munérique de ces deux distances pour l'exoplanète détectée à ce jour, ayant la valeur d'excentricité la plus élevée du catalogue.
Utiliser les données du catalogue des exoplanètes et l'outil histogrammes pour étudier la variabilité des excentricités des exoplanètes.
Estimez l'excentricité moyenne des exoplanètes détectées à ce jour
L'objectif de cet exercice est d'étudier les tailles des ellipses des exoplanètes qui ont été détectées à ce jour en utilisant les données du catalogue et les outils qui permettent de les visualiser.
Utilisez la boite à outil histogrammes pour étudier la taille des orbites des exoplanètes détectées à ce jour.
Comparez ces valeurs à celle de l'orbite de la Terre.
Utilisez la boite à outils diagrammes pour rechercher un lien entre taille et forme des ellipses.
Loi des aires : En un temps donné, le segment qui joint le centre du soleil au centre de la planète en orbite balaie une surface (aire) égale quelle que soit la position de la planète sur l'orbite.
Loi des périodes : Le carré de la période de révolution T d'une planète sur son orbite est proportionnel au cube du demi-grand axe de l'ellipse.
Pour un système planétaire donné, plus une planète est éloignée de son étoile, plus sa période est grande. Il est important de ne pas tirer de conclusions fausses sur les vitesses relatives de ces planètes. En effet, lorsque les orbites en question ne sont pas des cercles, la vitesse de déplacement de la planète n'est pas constante en accord avec la loi des aires.
Cette loi ayant été établie dans notre système solaire, il pouvait sembler à Kepler que la constante était universelle. Nous verrons par la suite que les lois de la mécanique permettent de calculer cette constante et que la valeur de la constante dépend de la masse de l'étoile. On montre que
où est la constante de gravitation universelle et la masse de l'étoile et celle de la planète.
Dans le cas très fréquent où la masse de l'étoile est très supérieure à celle de la planète , la relation devient et seule la masse de l'étoile intervient.
Utilisez la boite à outils diagrammes pour représenter pour l'ensemble des exoplanètes du catalogue leur période de révolution autour de son étoile en fonction de la valeur du demi-grand axe de leur ellipse.
En échelle logarithmique pour les deux axes montrez que les points sont tous sur des droites de même pente. Expliquez pourquoi, sur le diagramme obtenu, les points ne sont pas tous alignés sur la même droite.
L'attraction gravitationnelle, dans le cas d'un système planète-étoile isolé, permet de démontrer les lois de Kepler et de calculer le mouvement des deux corps. (Ces calculs sont présentés dans la partie comprendre).
Les animations ci-dessous illustrent l'effet de cette force d'attraction en fonction du rapport des deux masses en jeu. La croix rouge imobile correspond au centre de gravité (barycentre) de l'ensemble.
Au cours de la formation et de l'évolution d'un système planétaire, plusieurs types de migrations planétaires peuvent avoir lieu. On distingue trois types de migration de planètes. Les deux premiers types de migrations ont lieu lorsque la proto-étoile est toujours entourée d'un disque soit de gaz (migration de type I) soit de planétésimaux (migration de type II).
Après la dispersion du disque de gaz et de planétésimaux, les planètes peuvent toujours intéragir avec les petits corps rescapés de l'accrétion planétaire. Étant donnée la composition actuelle des planètes Uranus et Neptune (elles contiennent, toutes deux, un coeur solide : signe qu'elles ont du se former plus proche du Soleil), il y a de forts soupçons sur le fait que la migration tardive des planètes du Système solaire soit responsable de l'orbite actuelle d'Uranus et Neptune à, respectivement, 20 et 30 UA du Soleil.
Ces migrations planétaires sont notamment une explication plausible pour l'existence des " Jupiters chauds " : planètes trés massives et trés proches de leur étoile.
Pour en savoir plus sur les migrations planétaires : voir la page en suivant le lien.
Les exoplanètes, comme les planètes et autres corps de notre système solaire sont soumises à loi de la gravitation ou attraction universelle. Cette loi a été historiquement formalisée par Isaac Newton en 1687. Elle permet de modéliser, de façon générale, l'attraction entre des corps ayant une masse et, par conséquent, décrire le mouvement des corps massiques soumis aux forces de gravitation. Dans les pages qui suivent, le principe fondamental de la dynamique est utilisé pour étudié le mouvement de deux corps massiques isolés et démontrer les lois de Kepler.
On cherche à étudier le mouvement de deux corps de masses et assimilés à deux points matériels localisés aux points et qui sont leur centre de gravité. Le système de ces deux corps étant isolé on fait l'étude dans un référenciel supposé Galiléen d'origine arbitraire.
Les forces en présence
Lorsque deux corps massiques sont en présence l'un de l'autre, l'effet de la gravitation qui agit sur ces corps se traduit par une force d'attraction. Si ces deux corps sont assimilés à des points matériels localisés en leur centre de gravité, cette force est proportionnelle aux deux masses en jeu et inversement proportionnelle à la distance au carré entre les deux points.
Cette force explique aussi bien la chute des corps sur Terre que le mouvement d'une planète autour de son soleil ou d'une lune autour de sa planète.
On a la force exercée par sur , et, par symétrie, la force exercée par sur ,
La constante de gravitation universelle est
On remarque que en accord avec la loi de l'action et de la réaction pour un système isolé.
Dans ce référentiel (), le principe fondamental de la dynamique appliqué aux deux corps donne donc deux équations couplées
L'objectif est de connaître la position des centres de gravité et de en fonction du temps.
Le choix du référentiel de travail - Réduction à un problème à un corps
Si on considère l'ensemble des deux corps, le centre de gravité (ou barycentre) est défini de la façon suivante :
On dérive deux fois par rapport au temps et on utilise le fait que , on obtient alors
Ce qui traduit le fait que le barycentre est en mouvement rectligne et uniforme dans le référentiel (). Le repère barycentrique () dont l'origine est le centre de gravité des deux corps est donc lui aussi Galiléen.
On choisit de travailler dans le repère () ce qui permet de découpler le mouvement du barycentre des mouvements relatifs des deux corps.
Les équations du mouvement dans le repère barycentrique ()
On note ,
et on introduit
soit
Avec ces notations, les équations du mouvement dans le repère barycentrique deviennent :
L'expression des forces d'attraction gravitationnelle permet de réécrire ce système sous la forme :
et d'obtenir par différence
Cette équation peut être interprétée comme l'équation du mouvement d'un corps ponctuel fictif de masse (appeléee masse réduite du système) soumis à la force , soit . Dans la suite, ce point fictif sera noté .
Le problème à deux corps se réduit donc à un problème à un corps fictif unique. On aboutit à une équation unique pour le mouvement de dans laquelle n'apparait que l'inconnue .Cette équation est valable dans le repère baycentrique ().
L'expression de la force de gravité permet a priori de modéliser l'interaction entre plusieurs corps (étoiles, planètes, lunes, petits corps ...). Cependant, seul le problème à deux corps peut être mathématiquement résolu sans approximation.
On suppose ici que
Dans le repère barycentrique, le point qui correspond au corps le plus massif est immobile. En effet, la position du centre de gravité est confondue avec celle du corps le plus massif, comme le montre le calcul suivant.
devient
La masse du point est très voisine de celle du corps le plus léger. En effet,
Le corps 2 est donc en orbite autour du corps 1 dont le mouvement est négligeable. Cette situation particulière se présente fréquemment dans le cas d'une planète en interaction avec son étoile qui est souvent nettement plus massive qu'elle.
On part de l'équation obtenue précédemment dans le repère barycentrique ().
Le vecteur position est défini à partir du centre de gravité qui l'origine du repère.
est l'inconnue dont il faut déterminer l'évolution en fonction du temps.
On démontre en utilisant le théorème du moment cinétique que la trajectoire du point est plane. C'est la première loi de Kepler.
Le moment cinétique d'un point en mouvement est défini, pour un repère particulier et par rapport à l'origine du repère, de la façon suivante . Le théorème du moment cinétique (qui se déduit des lois de la mécanique) exprime le fait que la dérivée de ce moment cinétique est le produit vectoriel du rayon vecteur et de la force .
Pour le cas qui nous intéresse, dans le repère barycentrique () au point , on a donc
puisque
Le produit vectoriel de deux vecteurs colinéaires étant nul, on obtient . est donc un vecteur constant.
Or, du fait du produit vectoriel, le vecteur doit être orthogonal à chacun des vecteurs et . On montre ainsi que le plan défini en tout instant par les deux vecteurs et est invariant (orthogonal à l'axe défini par ce vecteur constant ). C'est le plan dans lequel le point évolue.
En conclusion, la trajectoire du point est bien plane.
A partir du moment cinétique, il est aussi possible de démontrer la loi des aires.
Nous avons démontré précédemment que la trajectoire du point restait dans un plan fixe qui est perpendiculaire à la direction du moment cinétique constant. Dans ce plan, il est donc possible de décrire le mouvement en coordonnées polaires. l'angle polaire sera noté ν (il s'agit de l'anomalie vraie définie dans La première loi de Kepler )
Dans ce système de coordonnées, on a , et
Il est possible dès à présent de simplifier l'experssion de l'accélération puisque la force n'a pas de composante que selon l'axe . On a
Avec ces notation, le moment cynétique constant s'écrit .
On démontre ainsi, puisque est un vecteur constant, que la quantité (appelée vitesse aréolaire) reste constante lors du mouvement du corps autour du point .
Or est l'aire balayée pendant la durée dt par le vecteur .
On en déduit donc que la surface balayée (en rouge sur la figure) pendant un temps donné est constante le long de la trajectoire et égale à . C'est la loi des aires de Kepler.
Remarque 1 : On peut également noter que si est une constante, alors le signe de est constant, ce qui traduit le fait que la rotation s'effectue toujours dans le même sens.
Remarque 2 : On peut utiliser la relation pour exprimer les dérivées de et de par rapport au temps et réecrire l'accélération du point en introduisant la fonction . en notant en notant
L'accélération devient alors . Cette expression sera utilisée par la suite pour trouver l'équation de la trajectoire du point .
Les solutions de l'équation du mouvementpermettent de déterminer la façon dont la position du point déterminée par les deux variables et évolue au cours du temps. Cependant, la résolution de cette équation différentielle, n'est pas possible analytiquement. Nous allons donc nous focaliser sur la relation entre la distance r et l'angle polaire qui modélise la trajectoire du point dans le plan de l'orbite.
Le résultat recherché est obtenu en utilisant la fonction introduite précédemment. Avec cette nouvelle variable, l'équation devient
soit après simplifications :
Remarque : Cette équation peut également être obtenue en utilisant la conservation de l'énergie mécanique (cinétique+potentielle) du point
L'équation à résoudre est une équation différentielle du second ordre (présence de ) à coefficients constants avec second membre constant .
La solution de cette équation est la somme d'une solution particulière de l'équation complète avec second membre et de la solution générale de l'équation sans second membre. On choisit comme solution particulière la solution constante égale à .
La solution générale de l'équation sans second membre est une fonction sinusoïdale de phase à l'origine et d'amplitude qui dépendent des conditions initiales du problème.
Au final, on obtient donc une solution de la forme .
Il est maintenant possible de repasser à la fonction et on obtient l'équation en polaire d'une conique :
Par identification aux paramètres de l'équation classique de l'ellipse , on a
et
Difficulté : ☆☆
Vous vous intéressez à la trajectoire du corps fictif de masse qui est en orbite autour du centre de gravité du système isolé planète-étoile, ceci dans le cas particulier d'une orbite circulaire.
Utilisez la relation pour montrer que, dans le cas d'une trajectoire circulaire, la vitesse est constante en module sur toute la trajectoire circulaire suivie par le point fictif autour du centre de masse des deux corps. Exprimez cette vitesse en fonction du rayon du cercle suivi par .
Exprimer la période de l'orbite (temps mis par le corps pour parcourir une fois le cercle).Vérifier, toujours dans le cas d'une trajectoire circulaire, la deuxième loi de Kepler.
Dans le cas particulier où ce point est confondu avec le centre de gravité de la planète ( ), simplifier les expressions obtenues précédemment.
Les animations ci-dessous montrent l'effet de cette force d'attraction en fonction du rapport des deux masses en jeu. La croix rouge imobile correspond au centre de gravité (barycentre) de l'ensemble.
Pour faire les calculs du mouvement d'une planète et de son étoile nous nous avons choisi un repère approprié pour simplifier les calculs (centre du repère au centre de masse, deux axes dans le plan de l'ellipse, ...), il faut maintenant tenir compte du fait que cette ellipse possède une certaine orientation dans l'espace et compléter la liste de paramètres que nous avons pour l'instant. Une orbite elliptique est décrite au moyen de deux plans (le plan de l'orbite et le plan de référence) et de six paramètres
Description de l'ellipse et du mouvement du corps dans le plan orbital
Description du plan orbital
L'orientation du plan orbital est donnée par rapport à un plan de référence.
Orientation du plan orbital par rapport au plan de référence
le transit désigne le passage d'une planète entre son étoile et nous. L'observation consiste à mesurer la variation du flux stellaire lors de ce passage de façon à obtenir des informations sur la planète étudiée et son orbite. La première planète ainsi découverte est HD209458b en 2000 (Charbonneau et al. 2000). Le principe est illustré sur la figure ci-dessous. Cette méthode permet de détecter la présence d'une exoplanète en orbite autout de son étoile et d'avoir accès à certains paramètres de l'ellipse. Pour en savoir plus, voir cours sur la méthode de détection des exoplanètes par la méthode des transits.
Dans notre système solaire, on peut observer depuis la Terre le transit des planètes Mercure et Vénus qui sont sur des orbites plus proches du soleil que celle de la Terre. Johannes Kepler a été le premier à prédire et pouvoir observer le transit de Mercure en novembre 1631, ainsi que celui de Vénus un mois plus tard.
Principe : On mesure par effet Doppler la vitesse d'éloignement ou de rapprochement de l'étoile, on peut détecter ainsi qu'il y a une planète en orbite et estimer la période de révolution. Pour en savoir plus,voir cours sur la méthode de détection des exoplanètes par la méthode des vitesses radiales.
Le cas d'un problème à deux corps, qui a été traité précédemment et qui permet de démontrer les lois de Kepler, est une approximation valable lorsque l'on peut négliger les forces de gravitation dues aux autres corps.
Le problème à N (N>2) corps se pose lorsque N corps massifs interagissent sans que l'on puisse a priori négliger certaines de ces interactions. Dans ce cas, on a un système de N équations à N inconnues qui sont les positions des centres de gravité des N corps de masse .
avec
Trouver analytiquement les solutions de ce système d'équations est impossible dans le cas général. Il faut recourir à des méthodes de résolutions approchées (perturbatives ou numériques).
Dans le cas particulier du problème à trois corps, on s'intéresse ici au mouvement d'un corps 'test' de masse négligeable qui subit l'attraction de deux corps plus massifs et. Le fait que la masse du corps soit négligeable permet de considérer que les mouvements de et ne sont pas perturbés par la présence de .
Pour simplifier la présentation du problème, nous allons nous restreindre au cas où est l'étoile et une planète, beaucoup moins massive, est en orbite circulaire autour de son étoile.
Le mathématicien Joseph-Louis Lagrange (1772) étudie ce problème. Il montre qu'il existe 5 points, dits de Lagrange (notés à ), pour lesquels les forces d'attraction de et se combinent de façon à ce que le corps "test" de masse négligeable ait la même période de révolution que les deux autres corps et les suive donc dans leur mouvement autour du centre de gravité de et.
En contradiction apparente avec les résultats obtenus dans le cadre du problème à deux corps, on peut trouver des corps de masse négligeable qui ont donc une période de révolution égale à celle de la planète mais qui ne sont pas sur la même orbite.
On montre que les points , et (parfois appelés Points d'Euler) correspondent à des positions instables alors que les points et correspondent à des positions stables. Les positions de ces deux derniers points ne dépendent pas des masses des points et. Dans le cas, du système Soleil/Jupiter, ce sont au voisinage de ces points que se trouvent les nombreux astéroïdes troyens qui suivent (ou précèdent) la révolution de la Terre autour du Soleil. D'autres planètes du système solaire sont suivies ou précédées également par des petits corps troyens (la liste des troyens détectés à ce jour dans notre système solaire est disponible sur le site du Minor Planet Center )
Appliquette Système Solaire
Vous allez utiliser l'appliquette pour visualiser les points de Lagrange (astéroïdes Troyens) stables et instables d'un système étoile-planète.
Attention : Cette appliquette utilise un système d'unités arbitraire pour les distances, vitesses, masses et temps de façon à ce que les valeurs numériques restent inférieures à un millier.
Choisissez pour commencer le système 'astéroïdes Troyens' proposé par l'appliquette
Laissez évoluer ce système, sans modifier les conditions initiales, jusqu'à 100 unités de temps pour vérifier que les deux petits corps positionnés aux points de Lagrange et ont bien la même période de révolution autour de l'étoile que la planète.
Vérifiez qu'initialement les trois corps ont des vitesses très voisines.
Modifiez la masse de l'un des astéroides (de 0.001 à 1 par exemple) et observer le changement qui apparait après un temps suffisamment long d'environ 50 unités. Expliquez ce qui se passe et proposez une explication.
Vous allez maintenant utiliser l'appliquette afin d'étudier les mouvements d'un corps de masse négligeable au voisinage des points de Lagrange instables et (situés de part et d'autre de la planète sur l'axe entre la planète et l'étoile).
Créez pour commencer un système à deux corps avec les valeurs numériques suivantes pour l'étoile et la planète : , , , . Trouver (en tâtonnant) une valeur de la vitesse initiale qui permet d'avoir une orbite quasi-circulaire pour la planète. Vérifiez que l'orbite de autour de est voisine d'une orbite circulaire.
Etude des points de Lagrange et
Vous allez maintenant ajouter à ce système un troisième corps de masse très faible aux positions qui correspondent aux points de Lagrange et .
On peut montrer par méthode perturbative que les positions de et par rapport au centre de masse du système étoile-planète sont données par les développements limités suivants:
et
avec
Calculer la valeur de et pour la configuration que vous allez simuler.
Ajoutez dans le système à deux corps, un point de masse négligeable sur l'axe étoile-planète à la position que vous avez précédemment calculée. Donnez à ce petit corps une vitesse initiale quelconque et observez ce que donne la simulation. Etudiez, en prenant quelques valeurs de vitesse différentes, l'impact sur la trajectoire du petit corps. Trouvez une valeur de vitesse initiale qui permet à ce petit corps d'avoir une vitesse angulaire proche de celle de la planète en révolution autour du soleil (au moins pendant une courte durée)
Reprennez les questions précédentes pour la position que vous avez calculée.
Pour vérifier que vous avez compris et retenu les notions de base de ce module
L’observation du centre de notre Galaxie a révélé la présence d’étoiles en orbite autour d’une masse invisible. L’observation de l’étoile S2 autour du centre galactique a été menée sur une dizaine d’années et a ainsi permis de mesurer la masse du corps central invisible. La concentration de masse associée à l’absence de rayonnements visible ou même infrarouge, laisse suspecter la présence d’un trou noir super massif. Dans cette première partie du mini projet, nous vous proposons d’étudier l’orbite de l’étoile S2 et de pouvoir ainsi déterminer la masse du trou noir central.
Appliquette pour le mini-projet
À l'aide de l'appliquette représentant l’orbite projetée dans le plan du ciel de l’étoile ainsi que le trou noir central SgrA*, repérer géométriquement le centre de l’ellipse.
Tracer la projection du grand-axe de l’orbite de S2 et évaluer l’excentricité de l’orbite. (L’excentricité évaluée par le rapport de la distance centre/foyer sur le demi grand axe reste préservée par projection par application du théorème de Thalès)
Lorsque l'étoile S2 est au périastre, elle se situe à une distance angulaire de 0,015". Notre Système solaire étant situé à 8000pc du centre galactique, estimer la distance du périastre au trou noir en unités astronomiques. En déduire le demi-grand axe de l'orbite réelle de l'étoile autour du trou noir en unités astronomiques.
À partir de la figure, déterminer la période de révolution de l'étoile. Grâce aux lois de Kepler, en déduire la masse du trou noir central de notre galaxie (nous sommes dans les conditions ou l'astre central est beaucoup plus massif que l'étoile S2).
De la même façon que pour estimer la masse du trou noir central de notre galaxie, il est possible d'estimer la masse de toutes les planètes possèdant un ou plusieurs satellites. Dans cet exercice, nous allons estimer la masse de Jupiter grâce aux orbites des satellites galiléens (Io, Europe, Ganymède et Callisto). Pour déterminer les positions des satellites par rapport à la planète, nous allons faire appel au serveur d'éphémérides de l'IMCCE.
Déterminez la distance Io-Jupiter pour 6 dates prises à 6h d'intervalle. Pour cela, grâce au serveur d'éphémérides, on se place dans un repére héliocentrique et en coordonnées rectangulaires. On peut ainsi obtenir la position de Io et de Jupiter. Que peut-on en conclure sur la forme de l'orbite du satellite autour de la planète.
Dans la suite de l'exercice, nous ferons l'approximation que l'orbite du satellite est circulaire. On se place à présent dans un repére géocentrique en coordonnées sphériques. Déterminez la distance apparente Jupiter-Io pour une vingtaine de date prise toute les 3h. (Vous pouvez vous servir d'un tableur afin de réaliser les calculs.) La distance Jupiter-Io est comptée positement vers l'ouest et négativement vers l'est.
Tracez cette distance apparente en fonction du temps et déterminer la période de révolution du satellite et le rayon de son orbite. En déduire la masse de Jupiter.
Recommencez l'étude avec les autres satellites galiléens (pensez à échantilloner différemment les dates car les périodes de révolution des satellites sont de plus en plus grandes).
Les masses déterminées avec chacun des satellites sont-elles toujours égales ? D'où viennent ces différences ?
Cet exercice peut aussi être réalisé en observant Jupiter et ses satellites avec un petit télescope ou une lunette et en prenant des clichés de la position du système à intervalle de temps régulier. Pour simuler ces observations, vous pouvez vous appuyer sur le logiciel libre stellarium et effectuer ces mesures à partir de capture d'écran à intervalle de temps régulier.
-------------FIN DU CHAPITRE -------------FIN DU CHAPITRE -------------FIN DU CHAPITRE -------------FIN DU CHAPITRE -------------
Dans cette partie, nous verrons les caractéristiques essentielles des petits corps du Système Solaire principalement les astéroïdes, les objets transneptuniens et les comètes.
La partie Décrire présente ces caractéristiques de façon aussi simple que possible et peut être lue avec des connaissances de niveau de terminale scientifique. En revanche, la partie Comprendre a pour objectif de démontrer quelques uns des résultats présentés précédemment et requiert un bagage en physique générale d’un niveau de licence.
En 2006, l’Union astronomique internationale (UAI) a défini les petits corps du système solaire comme étant tous les objets orbitant autour du Soleil qui ne sont ni une planète, ni une planète naine, ni un satellite. Les petits corps sont donc définis par opposition aux autres objets du système solaire. En pratique, les petits corps du système solaire sont les objets orbitant autour du Soleil mais n’ayant pas une masse suffisante pour avoir une forme presque sphérique et qui n’ont pas fait place nette dans leur voisinage.
La résolution de l'UAI définissant les différents objets du système solaire peut être trouvée ici : https://www.iau.org/static/resolutions/Resolution_GA26-5-6.pdf
L’usage distingue différents types de petits corps et en particulier, les astéroïdes, les objets de Kuiper et les comètes. Néanmoins, l’UAI n’a pas défini ces différentes appellations. Nous appellerons « astéroïdes » les petits corps dont l’ensemble de l’orbite est contenue à l’intérieur de celle de Neptune et « objets de Kuiper » ceux dont au moins une partie de l’orbite est située au-delà de celle de Neptune. D’autre part, nous définirons par « comète » les petits corps qui présentent une activité, c’est-à-dire l’émission de gaz et de poussières, sur une partie de leur orbite.
La majorité des astéroïdes est située entre Mars et Jupiter. Cet ensemble forme la ceinture principale d’astéroïdes. Néanmoins, l’orbite de certains astéroïdes peut croiser celle de la Terre, ceux-ci sont parfois appelés les géo-croiseurs. D'autres astéroïdes, appelés les troyens, circulent sur les mêmes orbites que les planètes et ceux dont des orbites sont comprises entre celles des planètes géantes sont appélés centaures. Les objets de Kuiper sont quant à eux majoritairement situés entre 30 et 50 unités astronomiques et forment la ceinture de Kuiper. Les comètes peuvent provenir d'au moins deux réservoirs : la ceinture de Kuiper et le nuage de Oort qui pourrait s'étendre jusqu'à 100 000 unités astronomiques du Soleil.
D’autre part, les poussières micrométriques gravitant sur le plan de l’écliptique et formant la lumière zodiacale peuvent aussi être considérées comme des petits corps même si leurs dimensions sont extrêmement faibles.
Le premier astéroïde a été découvert par Giuseppe Piazzi dans la nuit du 1er janvier 1801. Il s’agit de Cérès. Les calculs orbitaux montrèrent rapidement que ce nouvel objet se situait entre les planètes Mars et Jupiter. Entre 1802 et 1807, trois autres astéroïdes : Pallas, Junon et Vesta furent découverts sur des orbites très similaires. La recherche des petits corps du système solaire ne présentant pas d’activité cométaire, commença donc au début du 19ème siècle.
L’étude des astéroïdes présente un intérêt notable car elle permet de mieux comprendre les origines et la formation du système solaire dans son ensemble (voir notamment la page sur la composition des astéroïdes).
Afin de pouvoir connaître le plus grand nombre possible d'astéroïdes, des programmes de recherche automatique d’astéroïdes ont été conçus. Il s’agit de caméras CCD couplées à des télescopes dont les images prises à intervalles réguliers sont comparées automatiquement. Ces programmes permettent de découvrir de nombreux astéroïdes tous les mois. La page d’accueil du Minor Planet Center (http://minorplanetcenter.net/) donnant le nombre de découvertes lors du mois et de l’année en cours est particulièrement éloquente.
Dès qu’un nouveau petit corps est observé, il reçoit une désignation provisoire codant la date de la 1ère observation. Une fois que l’orbite d’un petit corps est clairement établie chaque corps reçoit un numéro permanent. Au 10 Févier 2015, 680 035 petits corps avaient été observés dont 427 393 ont une orbite bien définie (http://www.minorplanetcenter.net/mpc/summary).
La majorité des astéroïdes est située entre Mars et Jupiter et plus précisément entre 2.1 et 3.3 unités astronomiques du Soleil. D’autre part, l’inclinaison des orbites des astéroïdes et généralement très faible, inférieure à 4°, les astéroïdes gravitent donc autour du Soleil dans une zone proche du plan de l’écliptique. La distribution des astéroïdes en fonction du demi-grand axe de leurs orbites présente des discontinuités appelées « lacunes de Kirkwood » qui sont des zones de la ceinture d’astéroïdes peu peuplée. En effet, l’orbite des astéroïdes présents dans ces zones n’est pas stable sur de longues périodes à cause des effets de résonances orbitales avec Jupiter.
Sur la figure ci-contre, on remarquera une faible population d’astéroïdes dont le demi-grand axe est compris entre 0.8 et 1.8 unités astronomiques. Parmi ceux-ci, certains peuvent avoir une orbite croisant celle de la Terre, il s’agit des astéroïdes géo-croiseurs qui sont aussi appelés NEO (Near Earth Objects). Les astéroïdes ayant un demi-grand axe de 5.2 unités astronomiques ont une orbite stable et identique à celle de Jupiter et ils se répartissent aux alentours de certains points de Lagrange de Jupiter qui correspondent à des zones de stabilité https://media4.obspm.fr/public/FSU/pages_points-lagrange/impression.html.
Les astéroïdes ont des tailles extrêmement variables. Le plus gros et le plus massif, Cérès, présente un diamètre équatorial de 975 km. D’autre part, Cérès ayant une forme proche de la sphère il est aujourd’hui considéré comme une planète naine. Les grands astéroïdes, du type de Cérès, sont les plus rares, néanmoins ils représentent une fraction notable de la masse contenue dans la ceinture principale. Ainsi, avec une masse de 9.5 × 1020 kg, Cérès représente à lui seul environ un tiers de la masse totale de la ceinture principale. Mais, d’une manière générale, le nombre d’astéroïdes croît lorsque leur taille diminue. Aussi, la très grande majorité des astéroïdes ont des dimensions beaucoup plus petites et peuvent avoir des formes très éloignées d’une sphère. C’est en particulier les cas de l’astéroïde Eros qui a été survolé en 2000 par la sonde spatiale NEAR (Near Earth Asteroid Rendezvous) et qui possède une forme très allongée. Globalement, la distribution en taille des astéroïdes semble suivre une loi en puissance avec une surabondance d’astéroïdes présentant des diamètres d’environ 100 km et 5 km. Cette distribution fournit une contrainte forte sur les modèles concernant l’histoire collisionnelle de la ceinture principale d’astéroïdes.
Les observations spectroscopiques permettent de distinguer de nombreuses classes classes d’astéroïdes. Les 3 principales classes sont les astéroïdes de type C (Carbonés), de type S (Silicatés) et de type M (Métalliques). Les astéroïdes carbonés sont les plus nombreux et représentent 75% des astéroïdes observés. La distribution de ces différents types d’astéroïdes en fonction de la distance au Soleil est différente et est représentée sur la partie basse de la figure ci-contre. En effet, les astéroïdes carbonés sont majoritairement présents dans les zones externes de la ceinture d’astéroïdes alors que les astéroïdes de type S et M se concentrent dans les zones plus internes (voir figure ci-contre). Ces différentes distributions spatiales reflètent à la fois les processus liés à la formation de ces différents petits corps et l'évolution dynamique de la ceinture d'astéroïdes qui a pu conduire au mélange partiel des différentes populations d'astéroïdes.
Néanmoins, il faut aussi garder à l’esprit qu’une grande partie de notre connaissance sur la composition des astéroïdes provient de l’étude des météorites. En effet, la très grande majorité des météorités récoltées sur Terre provient de la ceinture d'astéroïdes. L'étude en laboratoire des météorites permet en particulier de préciser l'abondance des différents métaux (principalement du fer et du nickel), la nature des minéraux ainsi que d'obtenir des informations structuales sur la phase organique complexe et l'identification chimique des certains molécules simples. D’autre part, à quelques exceptions près, les astéroïdes ne contiennent pas de glace et ne présentent donc pas d’activité cométaire.
Les premières sondes spatiales lancées dans les années 1970 et 1980 n’ont pas survolé d’astéroïdes. Le premier survol d’un astéroïde, 951 Gaspra, a été effectué le 29 Octobre 1991 par la sonde Galileo. Depuis, de nombreux astéroïdes ont été survolés par des sondes spatiales. Les images acquises par les télescopes situés à la surface de la Terre ou en orbite autour de la Terre ne révèlent que très peu de détails dela surface des astéroïdes même sur Cérès qui est pourtant l’astéroïde le plus grand. Les survols d’astéroïdes par des sondes spatiales sont donc importants pour révéler des détails de la surface des astéroïdes, mais aussi leur forme précise ainsi que d’autres paramètres physico-chimiques comme la composition et la minéralogie de la surface ou la distribution massique interne.
Les principales missions d’études des astéroïdes sont NEAR (Near Earth Asteroid Rendezvous) qui s’est placée en orbite autour de l’astéroïde Eros en 2000 et Hayabusa qui malgré quelques difficultés techniques a pu prélever, puis ramener sur Terre, quelques poussières à la surface de l’astéroïde Itokawa. D’autre part, la mission DAWN s’est placée en orbite autour de l’astéroïde Vesta en 2011 puis en orbite autour de Cérès à partir du 6 mars 2015.
On désigne par objet trans-neptunien tout objet du système solaire dont l’orbite est entièrement ou en majeure partie située au-delà de celle de Neptune. Le premier objet trans-neptunien connu est Pluton qui a été découvert en 1930. Il fut considéré comme la 9ème planète du système solaire, avant d’être classifiée comme planète naine en 2006. L’existence de nombreux objets au-delà de Neptune avait été postulée par K.E. Edgeworth en 1949 puis par G.P Kuiper en 1951. Jusqu’en 1992, Pluton et son plus gros satellite Charon étaient les seuls objets trans-neptuniens connus. La découverte d’un nouvel objet de Kuiper en 1992 marque le début d’une recherche systématique des objets trans-neptuniens. Début 2015, plus de 1300 objets trans-neptuniens avaient été détectés (http://www.minorplanetcenter.net/iau/lists/t_tnos.html). Les orbites de ces objets correspondent à plusieurs groupes dynamiques. Nous verrons dans la partie suivante qu'on distingue en particulier les objets appartenant à la ceinture de Kuiper de ceux appartenant au disque épars.
Environ la moitié des objets trans-neptuniens présente des orbites au-delà de Neptune ayant des faibles excentricités (e < 0.2) et des demi-grand axes généralement compris entre 37 and 48 unités astronomiques. Ces objets trans-neptuniens sont appelés les objets de Kuiper classiques et forment la ceinture de Kuiper. Cette ceinture a une structure similaire à la ceinture principale d’astéroïdes et est structurée par des résonances avec l’orbite de Neptune. Ainsi de nombreux objets de Kuiper, dont Pluton, sont en résonance 2:3 avec Neptune.
D’autre part, un nombre croissant d’objets trans-neptuniens présentent des orbites de forte excentricité (e > 0.2), des périhélies au-delà de l’orbite de Neptune et des demi-grand axes supérieurs à 48 unités astronomiques. Ces objets font partie du disque épars. L’objet le plus massif de ce disque épars est Eris. On peut noter que l’orbite de Sedna avec un périhélie de 76 UA et un aphélie d’environ 900 UA est tout à fait exceptionnelle. Aussi Sedna peut être considéré comme un objet du nuage de Oort interne. La strcuture du nuage est discuté dans la partie traitant de l'orbite des comètes.
Les objets se trouvant aujourd'hui dans la ceinture de Kuiper et le nuage de Oort se sont très vraisemblablement formés dans la zone actuelle des planètes géantes. Ils ont ensuite été repoussé à des distance héliocentriques plus élevés suite aux interactions gravitationnelles avec les planètes géantes.
Les objets trans-neptuniens étant éloignés et de faible luminosité, il n’est pas aisé d’en déduire leurs caractéristiques physiques comme la taille ou la masse. Néanmoins, certains objets trans-neptuniens possèdent un ou plusieurs satellites. C’est, par exemple, le cas d'Eris qui possède un satellite, Dysnomie, et de Pluton autour duquel gravitent au moins 5 satellites. Si la période de révolution et le rayon de l’orbite des satellites sont mesurés, alors la masse du système peut être calculée grâce à la troisième loi de Kepler. Ainsi les masses d'Eris et de Pluton sont de 1.67 × 1022 kg et de 1.31 × 1022 kg. On remarquera que Eris est légèrement plus massif que Pluton.
Comme pour la masse, la mesure des diamètres des objets trans-neptuniens peut être délicate. Pour des objets dont l’orbite est bien connue, la méthode la plus précise pour mesurer la taille est l’occultation stellaire. Lors d’une occultation stellaire, l’objet étudié passe devant l’étoile et projette une ombre à la surface de la Terre. La mesure de la durée de l’occultation en différents points à la surface de la Terre permet alors de calculer le diamètre de l’objet occultant de manière relativement précise. Les rayons de Pluton et d'Eris sont actuellement estimés à 1153 ± 10 km et 1163 ± 6 km.
Connaissant la taille et la masse, la densité moyenne de ces objets est immédiatement déduite. La valeur de la densité moyenne renseigne sur la composition globale de ces objets. La densité volumique d'Eris étant d’environ 2.5 g.cm-3, celui-ci doit être composé majoritairement de roches. D’autre part, connaissant la taille et la magnitude, l’albédo peut aussi être estimé rapidement. Ainsi Eris possède un albédo de 0.96, ce qui fait d'Eris un des objets les plus réfléchissants du système solaire et qui indique la présence possible de glaces à la surface.
On notera que les occultations stellaires peuvent aussi révéler la présence d’atmosphère. Ainsi, Pluton est entouré d’une faible atmosphère dont la pression à la surface est estimée entre 6 et 24 microbars.
Malgré l’éloignement et la forte magnitude des objets trans-neptuniens, des observations photométriques et spectroscopiques dans le visible et le proche infrarouge permettent de contraindre la nature chimique des surfaces des objets de Kuiper. Sans surprise, l’objet dont la nature chimique de la surface est la mieux connue est le plus brillant, c’est-à-dire Pluton. En effet, les observations spectroscopiques dans le proche infrarouge ont montré que le composé principal est l’azote moléculaire (N2) avec des quantités plus faibles de méthane (CH4) et de monoxyde de carbone (CO). Dans le cas de Pluton, les observations actuelles permettent aussi de mettre en évidence des hétérogénéités à la surface. Le méthane a été détecté sur d’autres objets trans-neptuniens de grande taille, néanmoins l’azote et le monoxyde de carbone n’ont été détecté que sur Pluton (et Triton, le satellite principal de Neptune). En effet, même si le méthane n’est pas le composé le plus abondant, c'est le plus facile à détecter dans le proche infrarouge grâce à ces bandes d’absorption intenses. On notera que ces composés, l’azote, le méthane et le monoxyde de carbone, sont extrêmement volatiles et ne peuvent être présents en phase solide que sur des objets à la fois suffisamment froids pour les condenser et suffisamment massifs pour éviter l’échappement gravitationnel.
Néanmoins, tous les objets trans-neptuniens n’ont pas la même composition de surface. En effet, de la glace d’eau ainsi que des hydrates d’ammoniac ont été détectés à la surface de certains d’entre eux.
A ce jour (Mai 2015), aucun objet trans-neptunien n’a été survolé par une sonde spatiale. Néanmoins, la sonde spatiale New Horizons lancée le 19 janvier 2006 doit survoler Pluton le 14 Juillet 2015 à une distance d’environ 10 000 km. Il s’agira donc du 1er survol d’un objet trans-neptunien par une sonde spatiale. Après le survol de Pluton, la sonde New Horizons continuera sa trajectoire dans la ceinture de Kuiper et devrait encore survoler un ou deux autres objets. Néanmoins, cette phase de la mission et les objets à survoler ne sont pas encore définis.
Les comètes sont parfois visibles à l'œil nu, aussi sont-elles connues depuis l'Antiquité et ont parfois été représentées sur certaines œuvres d'art anciennes. Dans la croyance populaire, l'apparition d'une comète était associée à un bon ou à un mauvais présage suivant les cas. D'après Aristote, les comètes étaient un phénomène atmosphérique. La nature extraterrestre des comètes ne sera démontrée qu'en 1577 par Tycho Brahé qui effectua des mesures de parallaxes montrant que les comètes pouvaient être plus lointaines que la Lune. Certaines comètes ont eu une importance historique. En particulier, le retour de la comète 1P/Halley en 1758 avait été prédit par Edmond Halley. Cette prédiction permit alors de démontrer la théorie de la gravitation newtonienne. D'autre part, la comète 1P/Halley a été également la première comète survolée par des sondes spatiales en 1986.
Certaines comètes ont une orbite très elliptique, leur distance au Soleil peut varier de manière très significative au cours de l’orbite. Aussi lorsqu’elles se rapprochent, elles développent une activité : leur noyau est alors entouré d’une atmosphère constituée de gaz et de poussières. Bien que de faible densité, cette atmosphère peut être très étendue et brillante. C’est cette activité qui distingue les comètes des autres petits corps du système solaire.
De nos jours plus de 3800 comètes ont été observées (http://www.minorplanetcenter.net/). Durant les dernières décennies les comètes les plus remarquables ont été C/1996 B2 Hyakutake, C/1995 O1 Hale-Bopp, C/2006 P1 McNaught. D'autre part, depuis 1P/Halley, 6 autres comètes ont été survolées par des sondes spatiales : 26P/Grigg-Skjellerup, 19P/Borelly, 9P/Tempel, 81P/Wild 2, 103P/Hartley 2 et 67P/Churyumov-Gerasimenko qui est accompagnée par la sonde Rosetta depuis août 2014.
Les comètes ont généralement des orbites très excentriques. Leur distance au Soleil varie donc de manière très significative durant leur orbite. Elles sont généralement classifiées en fonction de leur période de révolution : les comètes ayant une période inférieure à 200 ans sont appelées comètes à courtes périodes alors que celles ayant une période supérieure à 200 ans sont appelées comètes à longue période. L’inclinaison des orbites des comètes à courtes périodes est généralement proche du plan de l’écliptique alors que les comètes à longues périodes ont une inclinaison aléatoire. Il est possible que les comètes à longues périodes puissent évoluer en comètes à courtes périodes suite à des interactions gravitationnelles avec les planètes géantes. On distingue aussi la famille des comètes de Jupiter qui possédent des périodes orbitales inférieures à 20 ans. Ces comètes sont ainsi nommées car leurs orbites sont déterminées par influence gravitationnelle de Jupiter. De plus, étant qu'elles ont des inclinaisons faibles et un sens de révolution autour du Soleil identique à celui des planètes, ces comètes trouvent très vraisemblablement leur origine dans le cainture de Kuiper.
Les comètes dont les orbites présentent de très longues périodes et des inclinaisons aléatoires proviendraient d’un réservoir appelé « nuage de Oort ». Celui-ci se situerait entre 50 000 et 100 000 unités astronomiques et aurait la forme d'une coquille sphérique entourant l'ensemble du système solaire. Il faut noter qu'aucun corps n'a jamais été observé à de telles distances du Soleil ; aussi l'existence du nuage de Oort n’est déduite que des observations des comètes à longues périodes. La densité de matière étant très faible à de telles distances du Soleil, les comètes du nuage de Oort ne se sont pas formées in-situ. Elles se seraient agrégées dans la zone des planètes géantes et en auraient été éjectées suite à des interactions gravitationnelles.
Les comètes dont les orbites présentent des périodes inférieures à 200 ans et des inclinaisons proches du plan de l’écliptique proviendraient quant à elles du nuage de Kuiper.
L’activité cométaire, c’est-à-dire le développement d’une faible atmosphère entourant le noyau solide est la caractéristique qui distingue les comètes des autres petits corps du système solaire. Le développement de cette faible atmosphère est dû à la sublimation des glaces contenues dans le noyau. Cette sublimation est induite par le réchauffement du noyau lorque celui-ci se rapproche du Soleil. Aussi lorsque les comètes sont situées à quelques unités astronomiques du Soleil, différentes grandes structures sont observables :
Ces différentes structures, ainsi que la petite taille et le faible albédo des noyaux cométaires, rendent très difficile l’observation directe de ces noyaux cométaires depuis le sol ou l’orbite terrestre. Aussi notre connaissance sur la composition des noyaux cométaires provient majoritairement de l’étude de leurs atmosphères.
La composition des comètes peut être déterminée par des observations depuis le sol et l’orbite terrestre dans tous les domaines de longueur d’onde. Ces observations ont permis de détecter dans la coma plus de 20 molécules gazeuses. L’eau est la molécule gazeuse la plus abondante ; le dioxyde de carbone (CO2) et le monoxyde de carbone (CO) présentent de fortes abondances de l’ordre de 10%. Ces molécules gazeuses présentent une grande diversité chimique et certaines sont relativement complexes comme le cyano-acétylène (HC3N) ou l’éthylène glycol (HO-CH2-CH2-OH). D’autre part, la nature de certains minéraux peut aussi être déterminée par des observations à distance. Ainsi les grains cométaires sont majoritairement constitués de silicates. Néanmoins, ces grains contiennent aussi une matière organique réfractaire. Cette composante organique présente en phase solide dans les grains cométaires est très difficile à observer depuis la Terre et n’a été mise en évidence de manière directe que par des mesures in-situ.
Notre connaissance sur la composition des noyaux cométaires provient donc majoritairement de l’étude de leurs atmosphères. Les observations de l’environnement cométaire permettent de déduire que les noyaux cométaires sont constitués d’un mélange de minéraux, de composés organiques réfractaires et de glaces de composés volatils.
Il faut noter qu’une partie des poussières extra-terrestres collectées dans la haute atmosphère terrestre ou au sol ont très vraisemblablement une origine cométaire. D’autre part, la sonde spatiale Stardust a collecté des grains dans l’environnement de la comète 81P/Wild2 puis les a ramenés sur Terre en janvier 2006. Ces échantillons disponibles au laboratoire permettent de déterminer de manière relativement précise la composition minéralogique des comètes.
La comète de 1P/Halley est la première à avoir été survolée par une sonde spatiale en 1986. A cette occasion, 5 sondes spatiales ont approché la comète 1P/Halley : Sakigake, Suisei, Vega-1, Vega-2 et Giotto. Le survol le plus proche a été réalisé par la sonde Giotto en mars 1986 à une distance de 596 km. Le noyau de la comète 1P/Halley a donc été le premier à avoir été photographié, ses dimensions ont alors été estimées à 16 × 8 × 7 km. La sonde Giotto a poursuivi sa course et a ensuite survolé la comète 26P/Grigg-Skjellerup en juillet 1992. Depuis 5 autres comètes ont été survolées :
Depuis août 2014, la sonde Rosetta accompagne la comète 67P/Churyumov-Gerasimenko durant une partie de son orbite autour du Soleil. Pour la première fois, il ne s’agit pas d’un survol du noyau, mais d’une mise en orbite autour du noyau ce qui permet de faire une étude plus détaillée de la comète mais aussi de suivre le développement de son activité lors de son approche du Soleil. D’autre part, un module nommé Philae a atterri sur la surface du noyau de 67P/Churyumov-Gerasimenko le 12 novembre 2014. Cette mission devait durer jusqu’en décembre 2015, mais a été prologée jusqu'en septembre 2016..
Les météorites sont des objets d’origine extra-terrestre retrouvés à la surface terrestre et donc disponibles en laboratoire pour des analyses physico-chimiques détaillées. Les météorites ne sont pas des petits corps du système solaire puisqu’elles n’orbitent plus autour du Soleil. Avant d’atteindre le sol terrestre, ces objets extra-terrestres traversent l’atmosphère terrestre. Lors de cette traversée atmosphérique, ils subissent un très fort échauffement produisant une traînée lumineuse appelée météore. Concernant les météorites, les chutes sont distinguées des trouvailles. Les chutes sont les météorites pour lesquelles la traversée atmosphérique a été observée alors que les trouvailles sont celles qui ont été retrouvées au sol sans que la date d’impact soit connue.
Les corps parents de la plupart des météorites sont des astéroïdes. L’étude des météorites en laboratoire permet d’obtenir des informations inaccessibles aux observations astronomiques classiques et de mieux comprendre l’évolution des astéroïdes mais révèlent aussi de précieuses informations concernant la formation du système solaire et les processus physico-chimiques présents dans le disque proto-planétaire. En particulier, certaines météorites contiennent une très grande diversité de composés organiques ou des inclusions qui constituent les tout premiers solides condensés dans la nébuleuse primodiale et qui permettent de dater l’origine du système solaire dans son ensemble.
On peut noter que le flux annuel de météorites à la surface de la Terre est estimé à 1010 – 1011 g/an et que ce flux est dominé par les micrométéorites ayant des dimensions de l’ordre de quelques micromètres.
L’étude des météorites a montré une très grande diversité de ces objets et elles sont subdivisées en deux grandes familles : les météorites différenciées et les météorites non-différenciées. Les météorites non-différenciées contiennent un matériau qui n’a pas été profondément altéré dans le corps parents et sont aussi appelées chondrites car elles contiennent généralement des chondrules qui sont des inclusions de forme sphérique ayant des tailles de l’ordre de quelques micromètres jusqu’au centimètre. Le matériau des météorites différenciées a subi de profondes modifications dans le corps parents.
Les météorites non différenciées, aussi appelées chondrites, sont elles même subdivisées en trois classes principales en fonction de leur composition : i.) les chondrites carbonées qui contiennent plusieurs pourcents en masse de carbone, ii.) les chondrites à enstatite qui sont riches en enstatite (minéral de formule MgSiO3) et iii.) les chondrites ordinaires qui sont les plus courantes. Bien que non-différenciées, les chondrites ont pu subir quelques processus d’altération depuis leur formation : métamorphisme thermique, altération aqueuse ou chocs. Néanmoins, certaines chondrites primitives possèdent une composition élémentaire très proche de celle du Soleil ce qui montre une formation très ancienne et quasi contemporaine de celle du Soleil. Ces météorites font partie des matériaux les plus primitifs que nous possédions sur Terre et peuvent alors fournir des contraintes concernant les conditions régnant dans le disque protoplanétaire ainsi qu’une datation du système solaire grâce à des méthodes de datation isotopique.
Les météorites différenciées sont elles aussi subdivisées en trois classes principales : i.) les achondrites qui sont principalement composées de minéraux, principalement des silicates ii.) les météorites métalliques constituées principalement de fer et de nickel et qui présentent des structures de cristallisation particulières appelées les figures de Widmanstätten et iii.) les lithosidérites constituées d’un mélange de minéraux et de métal. Celles-ci proviennent de corps parents suffisamment massifs pour avoir subi un processus de différentiation. Les météorites métalliques proviendraient du cœur métallique de ces corps parents alors que les achondrites proviendraient des régions externes de ces corps parents.
Certaines chondrites carbonées possèdent une composition élémentaire très proche de celle du Soleil malgré un appauvrissement pour les éléments les plus volatils (H, He, C, N,...). Ces météCes météorites se sont donc formées en même temps que le système solaire et n'ont pas subi de changement de composition important depuis. D’autre part, des études par spectrométrie de masse en laboratoire permettent de mesurer l’abondance d’éléments à l’état de traces. Il est donc possible de mesurer l’abondance d’éléments radioactifs à très longue durée de vie ainsi que l’abondance des éléments produits par ces processus radioactifs. Ces méthodes de datation radiométrique ont permis de déterminer l’âge de formation des météorites les plus primitives : 4.56 milliards d’années. Cette date correspond à la condensation des premiers solides dans le système solaire. On pourra retenir que sans l’étude en laboratoire des météorites, l’âge de notre système solaire ne serait pas connu avec autant de précision.
La datation absolue de roches anciennes est généralement réalisée grâce à la mesure de l'abondance d'éléments naturellement radioactifs. En effet, au cours du temps, l'abondance des noyaux parents décroît alors que celle des noyaux fils augmente. A un instant t, l'abondance des noyaux parents dépend de son abondance initiale et de son temps caractéristique de décroissance : . Il faut noter que le temps de demi-vie , qui correspond au temps nécessaire pour que l’abondance initiale soit divisée par deux, est souvent utilisé à la place du temps caractéristique de décroissance. Ces deux temps sont liés par la relation .
A l’instant t, l'abondance des noyaux fils est égale à si tous les noyaux parents se désintègrent pour donner le même noyau fils. En combinant, les deux équations précédentes, peut être exprimée en fonction de qui sont les deux quantités mesurables. Néanmoins,l'abondance initiale de noyaux fils est généralement inconnue. La roche ayant pu subir une différentiation chimique lors de sa formation, elle peut être inhomogène. Les différentes mesures réalisées sur la même roche seront donc normalisées à l'abondance d'un isotope non radioactif de l'élément fils. En mesurant plusieurs échantillons de la même roche, l’âge de celle-ci ainsi que l'abondance initiale de l'élément fils radioactif pourront être déterminés (voir exercice dans la partie "Se Tester".).
Un des couples d'éléments couramment utilisé est le couple 87Rb (élément parent) et 87Sr (élément fils), les abondances de ces deux éléments sont normalisées à l'abondance de 86Sr. Ces abondances sont mesurées par spectrométrie de masse en laboratoire. Les différentes mesures seront représentées dans un diagramme 87Sr / 86Sr en fonction de 87Rb / 86Sr. La pente de la droite obtenue est liée à l’âge de formation de la roche alors que l’ordonnée à l’origine donne la rapport initial 87Sr / 86Sr. En normalisant à un isotope stable et en effectuant différentes mesures sur une roche hétérogène, les deux inconnues précédentes peuvent donc être déterminées et en particulier l’âge de la roche. Cette méthode n’est valable que si le système étudié, la roche ou la météorité, est clos, c’est-à-dire si aucun des éléments étudiés n’a pu diffuser à l’extérieur du système.
La datation absolue des météorites, et en particulier des inclusions les plus anciennes qu’elles contiennent, nous permet de dater la formation du système à 4.56 milliards d'années. Sans l'étude des météorites et des méthodes de datations liées aux isotopes naturellement radioactifs, l'obtention d'un chiffre aussi précis ne serait pas possible. L’étude des météorites et de l’ensemble des échantillons d’origine extra-terrestre disponibles au laboratoire permet donc de dater de manière absolue la formation du système solaire dans son ensemble.
L’atmosphère des comètes est appelée coma. La principale source des espèces gazeuses présentes autour des noyaux cométaires est la sublimation des glaces contenues dans celui-ci. Les molécules gazeuses produites directement depuis la surface ou la sous-surface du noyau sont appelées « molécules mères ». Une fois dans l’environnement cométaire, ces molécules gazeuses sont soumises au flux ultraviolet du Soleil et peuvent se photo-dissocier en de nouvelles espèces gazeuses plus petites appelée « molécules filles ». Par exemple, la molécule d’eau peut se photo-dissocier en radicaux hydroxyle (OH) et en atomes d’hydrogène. La distribution spatiale autour du noyau est l'élément essentiel pour comprendre l’origine des espèces gazeuses. Afin de déterminer les mécanismes de production des différentes espèces gazeuses présentes dans la coma puis remonter aux abondances des molécules présentes dans le noyau ; les distributions spatiales mesurées grâce aux observations astronomiques doivent être comparées à un modèle. Le modèle le plus simple pour décrire la distribution spatiale des molécules est le modèle de Haser.
Dans le cadre du modèle de Haser, la densité des molécules mères dans la coma n’est régie que par l’expansion générale de la coma et par leur photolyse sous l’effet du rayonnement UV solaire. La densité volumique (en molécules.m-3) en fonction de la distance r au noyau, vérifie donc l’équation de conservation de la masse : où est la durée de vie de la molécule mère étudiée dans l’environnement cométaire et la vitesse d’expansion des gaz dans la coma.
En supposant un état stationnaire, une vitesse constante dans la coma et une symétrie sphérique, l’équation précédente se simplifie : où est appelée « longueur d’échelle parent », elle correspond à la longueur caractéristique de photolyse de la molécule mère.
La condition initiale permettant d’intégrer l’équation différentielle précédente est donnée par le taux de production Q (en molécules.s-1) qui correspond au nombre de molécules mère émises depuis la surface du noyau en 1 seconde. L’intégration de l’équation précédente conduit à : .
Pour les molécules filles, la densité volumique dans la coma est régie par les mêmes processus auxquels il faut rajouter la production directement dans la coma par photolyse de la molécule mère. Dans ce cas, l’équation de conservation de la masse s’écrit avec deux termes dans le membre de droite ; i.) un terme de production qui correspond à la photodissociation de la molécule mère et ii.) un second terme de destruction correspond à la photodissociation de l’espèce fille considérée. Soit et les densités volumiques des espèces mère et fille, ainsi que et leurs temps de vie respectifs dans l’environnement cométaire. L’équation de conservation de la masse pour une espèce fille s’écrit alors : En utilisant les mêmes hypothèses géométriques et dynamiques que précédemment et en imposant la condition initiale n_F (0)=0, on trouve : . Dans cette dernière équation, Q est le taux de production de la molécule mère dont est issue la molécule fille, la longueur d’échelle parent et la longueur d’échelle fille, qui correspondent respectivement aux longueurs caractéristiques de destruction par photolyse des espèces mère et fille.
Ce modèle est très simple, voir simpliste au regard de nos connaissances actuelles sur les comètes. En particulier, la production de gaz à la surface du noyau est inhomogène et les environnements cométaire n’ont pas une symétrie sphérique. Néanmoins, ce modèle permet d’obtenir des ordres de grandeurs pertinents en particulier pour les longueurs d’échelles et des versions modifiées de ce modèle continuent à être utilisées pour calculer les taux de production des espèces gazeuses à partir des observations.
Après la découverte de Dysnomie, l’unique satellite de Eris, son orbite a pu être déterminée (Brown et al., 2007, Science, 316, 1585). Son demi-grand axe et sa période sont 37400 ± 180 km et 15.773 ± 0.002 jours, respectivement.
A partir de la 3ème loi de Kepler, calculer la masse du système Eris-Dysnomie.
Le rayon de Eris a été déterminé à 1163 ± 6 km par occultation stellaire (Sicardy et al., 2011, Nature, 478, 493). Calculer la densité de Eris. Que peut-on en conclure ?
Difficulté : ☆
Le flux de micrométéorites à la surface de la Terre est estimé à (40 ± 20) × 106 kg.an-1 pour des grains ayant des masses comprises entre 10-9 et 10-4 grammes.
Les glaces cométaires sont majoritairement constituées de glace d’eau. Dans la coma, la vitesse d’expansion des molécules d’eau est d’environ 0.5 km.s-1 et la photodissociation de l’eau conduit à la formation de radicaux hydroxyles (OH) et d’atomes d’hydrogène. Les temps de vie typiques de l’eau, du radical hydroxyle et de l’atome d’hydrogène sont de 6 × 104 s, 2 x 105 s et de 106 s, respectivement. En supposant une vitesse d’expansion des radicaux hydroxyles égale à celle des molécules d’eau et de 12 km.s-1 pour l’atome d’hydrogène, calculer les dimensions caractéristiques des coma de H2O, OH et H.
Difficulté : ☆☆ Temps : 15 min
Sachant que le temps de demi-vie du Rubidium 87 est de 48.8 milliards d'années et connaissant les abondances suivantes, calculer l'âge de la météorite dans laquelle ces abondances ont été mesurées.
0.758 | 0.74864 |
0.7255 | 0.7465 |
1.52 | 0.79891 |
1.49 | 0.79692 |
1.555 | 0.80152 |
1.685 | 0.80952 |
0.1542 | 0.7091 |
0.1533 | 0.70895 |
Le but de ce projet est de retrouver les principales caractéristiques dynamiques des astéroïdes et des objets trans-neptuniens à partir des données orbitales compilées par le « Minor Planet Center » ou le « Jet Propulsion Laboratory ».
Les lacunes de Kirkwood correspondent à des minimas de la distribution des astéroïdes en fonction de leur demi-grand axe (ou en fonction de leur période orbitale). Ces lacunes sont dues à des résonances avec Jupiter.
A partir des éléments orbitaux des astéroïdes numérotés, dont l’orbite est connue, qui sont disponibles sur la base de données du JPL (http://ssd.jpl.nasa.gov/?sb_elem), tracer l’histogramme montrant la distribution des astéroïdes en fonction de leur demi-grand axe. Pour construire cet histogramme, on utilisera préférentiellement un pas de 0.01 UA. Afin de mettre en évidence le lien entre la position des lacunes et les résonances avec Jupiter, on représentera sur cet histogramme, la localisation des résonances 4:1, 3:1, 5:2 7:3 et 2:1. On rappelle que le demi-grand axe de Jupiter est de 5.2 UA.
Afin de mettre en évidence, les différentes familles dynamiques d’astéroïdes, on pourra aussi représenter l’ensemble des astéroïdes dans un graphique montrant l’inclinaison en fonction du demi-grand axe.
Les objets trans-neptuniens sont contenus dans deux grands réservoirs : la ceinture de Kuiper classique et le disque épars. Le principal critère permettant de distinguer ces deux réservoirs est l’excentricité de l'orbite des objets qu’ils contiennent. De plus, la ceinture de Kuiper est structurée par des résonances avec Neptune, dont le demi-grand axe est de 30.1 UA.
A partir des éléments orbitaux des objets de la ceinture de Kuiper classique (http://www.minorplanetcenter.net/iau/lists/TNOs.html) et de ceux du disque épars (1) (http://www.minorplanetcenter.net/iau/lists/Centaurs.html) on pourra représenter l’ensemble de ces objets dans deux graphiques montrant l’excentricité et l’inclinaison en fonction du demi-grand axe de ces objets. Afin de mettre en évidence la structuration de la ceinture de Kuiper classique par Neptune, on pourra indiquer sur ces graphiques la positon des résonances 4:3, 3:2, 5:3, 7:4 et 2:1. De plus, sur le graphique de l’excentricité en fonction du demi-grand axe, on pourra aussi tracer les courbes représentant un périhélie constant. Afin de bien distinguer la ceinture de Kuiper classique du disque épars, on pourra pour chacun des deux graphiques proposés utiliser deux axes des abscisses ; un premier en échelle linéaire allant de 28 à 60 UA puis un second en échelle logarithmique s’étendant de 25 à 1 000 UA.
(1) Ce fichier issu du Minor Planet Center contient à la fois les éléments orbitaux des centaures et des objets du disque épars.
-------------FIN DU CHAPITRE -------------FIN DU CHAPITRE -------------FIN DU CHAPITRE -------------FIN DU CHAPITRE -------------
pages_fluide-temperature-setester-qcm/definitions.html
pages_fluide-temperature-setester-qcm/effetdeserre-exo.html
Revoir la définition de l'effet de serre : l'absorption du rayonnement thermique n'est pas le seul critère !
pages_fluide-temperature-setester-qcm/profils.html
Revoir les définitions.
pages_orbite-setester/exercices.html
La masse de l'étoile est bien supérieure à celle de la planète mais n'est pas infinie
Quand la planète est au périastre, elle est au plus proche de l'étoile.
pages_se-tester/qcm.html
pages_se-tester/qcm.html
pages_se-tester/qcm.html
pages_se-tester/qcm.html
pages_se-tester/qcm.html
pages_se-tester/qcm.html
pages_se-tester/qcm.html
pages_se-tester/qcm.html
pages_se-tester/qcm.html
pages_sphere-surface/exo1.html
Ce graphe montre la relation inverse entre la gravité et la taille de cratère marquant la transition entre cratère simple et cratère complexe sur Terre, la Lune, Mars et Mercure. Sur Terre, la transition entre cratère simple et cratère complexe de fait pour un diamètre transitoire d'environ 3 km.
pages_sphere-surface/exo2.html
Attention aux unités !
avec ( en , en , en et en ) AN:
Le choc libère une énergie équivalente à un séisme de magnitude sur l’échelle de Richter. Si la météorite était composée de fer :
sont nécessaires à la vaporisation totale du projectile soit 25% de l’énergie cinétique.
La formule est dans le cours !
Le diamètre du cratère transitoire est de 19.4 km.
Il s’agit d’un cratère complexe dont le diamètre final sera d’environ 23 km et la profondeur de 2-3 km.
Exprimer le volume d’un bol de diamètre et de profondeur .
Le volume d’un bol de diamètre et de profondeur est égal à avec . Le volume de matière déplacé est environ soit .
Le rayon de la Lune est 1737 km.
Si une couche d’éjectas de 5 m d’épaisseur recouvre toute la Lune, le volume total de dépôts est : . Le nombre d’impacts nécessaires pour créer cette couche est / .
pages_sphere-surface/exo3.html
La formule est dans le cours. On rappelle que la densité de la surface de la Lune est . Son accélération de pesanteur est 1.62 .
Le cratère transitoire aura une taille d’environ 5.7 km. Il sera donc de type simple et le diamètre final du cratère sera d’environ 7 km.
La trajectoire d'un éjecta est balistique. On peut aisément montrer que la portée de cette excavation est où est la vitesse initiale de l'éjecta et l’angle d'éjection. En prennant °, on trouve que les roches excavées sont éjectées à environ 150 km du cratère principal.
L'accélération de pesanteur sur Mercure est 3.7 .
Sur Mercure, le cratère principal aura un diamètre final de près de 6 km et la portée de l'excavation sera d'environ 70 km (pour un angle d'ejection de 45°). Plus la gravité est forte, plus les cratères secondaires sont proches du cratère principal.
pages_sphere-surface/exo4.html
L’épaisseur du régolite lunaire augmente avec son âge. De plus, le taux de croissance de ce régolite était beaucoup plus important dans la jeunesse du Système solaire, avant 3.5 milliard d’années, quand le flux météoritique était plus élevé. Ces résultats ont été obtenus à partir d’observations radar de la Lune acquises à la longueur d’onde de 70 cm depuis l’Observatoire d’Arecibo (Porto Rico).
pages_sphere-surface/exo5.html
La formule est dans le cours.
En utilisant la formule (page sur l'érosion éolienne) , la taille optimale des grains susceptibles d'être déplacés est de 4 mm.
Pour déplacer un grain de 4 mm de diamètre sur Triton, il faut une vitesse de cisaillement minimum de 12.7 m/s (obtenue via la formule ) et donc une vitesse de vent à 1 mètre de la surface d’environ 330 m/s (obtenue via la formule ) avec égal à 1/30 de la taille des grains).
Pour former des dunes il faut que le vent à 1 m du sol soit significativement plus fort que la vitesse du son. De tels vents existent sans doute mais à plusieurs dizaines de mètres au-dessus de la surface. Il est donc peu probable de trouver des dunes sur Triton sauf si ce satellite avait une atmosphère plus dense dans le passé.
Elle est produite par impact ou cryo-volcanisme.
pages_sphere-surface/exo32.html
On rappelle: viscosité de l’eau : 1540 10-6 Pas, densité de l’eau : 1000 kg/m3, densité du quartz : 2650 kg/m3.
Environ 5 jours.
On rappelle: densité de l’argile: 1 700
Plus de 1000 jours !
Toute agitation remettant en suspension les particules très fines, comme les argiles, ces dernières ne devraient jamais sédimenter. En fait elles floculent, forment des agrégats plus gros qui finissent par décanter rapidement.
pages_sphere-surface/exo6.html
Le rayon de Io est 1822 km
Les 9 volcans produisent 9*50 km3 par an ce qui correspond à une couche d’épaisseur km/an soit un renouvellement de la surface d’environ 1 cm/an.
Avec un taux de renouvellement d’environ 10-2 m/an il faut environ 105 soit 100 000 ans pour renouveler la surface sur 1 km.
pages_sphere-surface/pageapplicrater.html
10 km pour Mercure, 3.2 km pour la Terre, 17.5 km pour la Lune et 8 km pour Mars. Voir aussi l'exercice Catégorie de cratères.
La présence d’une atmosphère freine l’impacteur ce qui a pour effet de réduire sa vitesse d’arrivée à la surface et donc les effets de l’impact (notamment la taille du cratère et l’amplitude du séisme qu’il génère).
La réponse est -0.22. Voir la page Processus de cratérisation.
La densité de la Terre est .
L'impacteur qui a donné naissance à Meteor Crater faisait environ 50 m de diamètre. La profondeur du cratère est de l’ordre de 180 m.
La densité de la Lune est .
L'impacteur qui a donné naissance au cratère Tycho faisait probablement entre 8 et 10 km de diamètre. Cet impact a eu lieu il y a seulement 109 millions d'années. Tycho a une profondeur d’environ 4.8 km.
pages_sphere-surface/applisegregation.html
où est la dimension d'une distance et celle d'un temps.
Le taux de sublimation est si rapide que ni le bombardement micro-météoritique, ni le « ion sputtering » ne peuvent rivaliser aux basses latitudes de Callisto, Ganymède et, peut être, d’Europa. Aux hautes latitudes d’Europa, en revanche, ce phénomène est a priori négligeable.
pages_sphere-surface/applierosion.html
Lorsque les grains sont suffisamment grands, la vitesse seuil nécessaire à leur mise en mouvement augmente en sqrt(d). Dans ce cas, le vent s’oppose simplement à la gravité (c’est-à-dire au poids des grains). Lorsque les grains sont petits, les forces de cohésion (entre particules ou avec le sol) ne sont pas négligeables et prennent même le dessus sur la gravité ; la vitesse seuil nécessaire à la mise en mouvement des petites particules varie alors en 1/d (approximation). A la transition entre ces deux régimes, il existe une taille limite de grain pour laquelle la vitesse seuil de mise en mouvement est minimum.
Dans l'ordre: Vénus, Titan, Terre, Mars. Les vitesses seuil de mise en mouvement des grains sur Vénus et Titan sont nettement plus faibles que sur Mars et la Terre. Ceci est lié au fait que les densités atmosphériques de Vénus et Titan sont significativement plus grandes. Sur Vénus, de plus, la viscosité de l’atmosphère est élevée et sur Titan la pesanteur est faible (7 fois plus faible que sur Terre).
Vénus. Sur Vénus, la taille des grains les plus facilement transportables (c’est-à-dire requérant la vitesse d’écoulement la plus faible pour être mis en mouvement) est d’environ 94 μm (le minimum de la courbe) contre 160 μm sur Titan, 220 μm sur la Terre et 1100 μm sur Mars.
Titan - Atmosphère Vénus – Glace d’eau.
Cette taille correspond à la taille des grains les plus facilement transportable par le vent sur Terre (le minimum de la courbe est atteint pour un diamètre de grain de 220 μm).
En raison de la plus grande densité de l’eau par rapport à l’air, l’activité fluviale sur Terre déplace plus facilement (seuil d’entrainement plus petit) des grains en moyenne plus gros (le minimum de la courbe est atteint pour un diamètre de grain de 560 μm) que l’activité éolienne.
L’activité fluviale sur Titan déplace plus facilement (seuil d’entrainement plus petit) des grains en moyenne plus gros, que ceux-ci soient composés de glace d’eau ou de matière organique (le minimum de la courbe est atteint pour un diamètre de grain de 400-500 μm) que l’activité éolienne.
Selon ce modèle, le transport fluvial est plus facile à initier sur Titan.
La faible pesanteur de Mars rend le mouvement des particules les plus grandes plus facile. Pour les petits grains, en revanche, les forces de cohésion qui s’opposent au transport fluvial sont les mêmes sur la Terre et sur Mars.
Les seuils de vitesse d’écoulement requis pour la mise en mouvement des grains sédimentaires sur Mars sont particulièrement élevés. On rappelle que la vitesse de l’écoulement (lorsqu’il est turbulent) à 1 m du sol est liée à la vitesse de cisaillement par la relation (approximative) : Il faut donc, à 1 m du sol, un vent minimum de 37 m/s soit de plus de 130 km/h pour déplacer des particules sur Mars! Sur la planète rouge, les sédiments sont donc transportés essentiellement pendant des épisodes de tempête.
pages_fluide-temperature/calculs.html
Ce sont des applications directes des définitions.
Chercher parmi les espèces présentes celles qui absorbent les UV, et comparer leurs abondances à celle des espèces qui rayonnent bien dans l'IR thermique.
Les bons absorbants en UV sont l'ozone et le dioxyde de soufre, tandis que le dioxyde de carbone, la vapeur d'eau, le méthane et le monoxyde de carbone sont de bons gaz à effet de serre. L'existence d'une stratosphère est donc possible, même si probablement moins marquée que sur Terre du fait du plus grand effet de serre.
Applications directes des formules du cours.
-6,8 K/km. Le gradient adiabatique humide sera inférieur en valeur absolue.
Températures importantes : au sommet de la troposphère, au sommet de l'atmosphère (voir le modèle à deux couches).
pages_fluide-dynamique/coord-sphere.html
Projeter les vecteurs sur les axes X,Y et Z.
Calculer en gardant en tête que est un vecteur local qui dépend de la position, donc varie avec les coordonnées .
Exprimer , et en fonction des angles et .
Utiliser les résultats des deux questions précédentes qui donnent les expressions de , , , , et .
pages_fluide-dynamique/equations-dynamique.html
C'est le même exercice que celui sur la force de Coriolis, mais sans négliger la composante .
pages_fluide-dynamique/equations-dynamique.html
Aussi dans ce cas, c'est un exercice similaire à celui de la force de Coriolis.
pages_fluide-dynamique/equilibres-dynamiques.html
pages_fluide-dynamique/setester-bases.html
Un vortex de lavabo est -il un système géostrophique, cyclostrophique ou aucun des deux ?
pages_fluide-dynamique/exo-coriolis.html
La gêne ressentie serait due à une force de Coriolis trop grande lorsqu'ils se déplacent sur la surface interne de l'anneau ...
pages_mecanique-orbite/decriretrois.html
On a la relation qui donne la distance en fonction de l'angle :
Au périastre, et c'est la distance minimale A l'apoastre, et c'est la distance maximale
Dans le catalogue, classer les planètes par value d'excentricité décroissante pour identifier l'exoplanète recherchée.
L'exoplanète qui à l'orbite la plus excentrique est HD 20782 b avec une valeur de 0.97. La valeur de son demi-grand axe a est 1.381 UA.
Tracez l'histogramme des excentricités et pensez à diminuer le nombre de bins de l'histogramme pour obtenir simplement une valeur approchée de la moyenne des excentricités.
L'excentricité moyenne est de 0.18.
pages_mecanique-orbite/decriretrois.html
La taille d'une ellipse est donnée par la valeur de son demi-grand axe a. Pensez à utiliser l'échelle logarithmique.
La Terre se situe à 1 AU de son étoile sur une orbite proche d'un cercle. a=1 AU
Il n'y a pas de corrélation évidente entre ces deux paramètres.
pages_mecanique-orbite/decrirecinq.html
Pensez à utiliser les échelles logarithmiques pour mettre en évidence une relation du type
Utiliser le code couleur de l'outil diagramme pour visualiser la masse de l'étoile et vérifier que ce paramètre explique largement le décalage des droites qui sont toutes parallèles.
On voit que de façon très globale (malgré quelques points aberrants très évidents en haut à gauche du diagramme), les points correspondant à des planètes autour d'étoiles plus massives (en rouge) ont des periodes de révolution plus grandes que ceux qui (en vert) correspondent à des planètes autour d'étoile moins massives. En selectionnant des points particuliers on peut vérifier que pour un système donné, les points sont parfaitement alignés sur une droite (exemple du systeme de Kepler 42 avec ses planètes b, c et d en vert)
pages_mecanique-orbite/exercicecirculaire.html
L'accélération en coordonnées polaires s'écrit de façon générale . Lorsque le mouvement est circulaire est une constante.
L'accélération subie par le point est colinéaire à du fait de la force centrale et ne dépend que de . Pour un mouvement circulaire, on a de plus , la force est constante en module. .est une constante, le mouvement circulaire est donc obligatoirement uniforme. La vitesse est égale à , d'où D'ou la relation . Il est également possible de faire le calcul de cette vitesse est d'utiliser la conservation de l'énergie mécanique : cinétique + potentielle de gravitation.
La circonférence d'un cercle est égale à .
Le temps mis pour faire un tour est . . C'est la troisième loi de Kepler qui dit que
On trouve et .
pages_mecanique-orbite/lagrange.html
La position des points de Lagrange est obtenue en considérant que la masse des astéroïdes est négligeable. Si on l'augmente, l'approximation n'est plus valable.
convient.
On trouve et
par exemple
par exemple
pages_disque-petitscorps/exo-eris-masse-densite.html
m = (1.66 ± 0.02) × 1022 kg
ρ = 2.53 ± 0.07 g.cm-3, Eris est vraisemblablement constitué d'un mélange de glace (env. 1 g.cm-3) et de minéraux (env. 3-5 g.cm-3)
pages_disque-petitscorps/exo-flux-micrometeorites.html
pages_disque-petitscorps/exo-coma.html
L(H2O) = 3 × 104 km L(OH) = 105 km L(H) = 12 × 106 km
pages_disque-petitscorps/exo-datation.html
Dans le graphique représentant 87Sr/86Sr en fonction de 87Rb/86Sr, les points forment une droite de pente 0.0659 ± 0.0001. On en déduit un âge de 4.49 ± 0.01 milliards d'années.
pages_disque-petitscorps/projet.html
pages_disque-petitscorps/projet.html