Comprendre

Auteurs: Loïc Rossi, Emmanuel Marcq

Transfert radiatif

Auteurs: Loïc Rossi, Emmanuel Marcq

Processus radiatifs

L'interaction entre le rayonnement (ondes du champ électromagnétique) et la matière (chargée au niveau microscopique : électrons, noyaux atomiques) se caractérise par trois types de processus radiatifs : l'absorption, l'émission et la diffusion.

Émission

Un processus d'émission intervient lorsque la matière cède localement de l'énergie au champ électromagnétique (création d'un photon). Cette émission est qualifiée de thermique lorsque le champ et la matière sont à l'équilibre thermodynamique à une même température T (on peut voir alors leur interaction comme celle d'un gaz de photons perpétuellement absorbés et émis par la matière environnante, et se mettant en équilibre de la sorte). Il existe également des phénomènes d'émission non thermiques (ex: fluorescence, émission laser), mais ils sortent du cadre de ce cours ; certains de ces processus sont vus dans ce chapitre.

Illustration d'un processus d'émission spontanée
EmissionSpontanee.png
Émission spontanée d'un photon par un système matériel possédant au moins deux niveaux d'énergie. L'énergie totale est conservée dans ce processus.

Absorption

Un processus d'absorption a lieu lorsqu'inversement, le champ électromagnétique cède de l'énergie à la matière environnante (destruction d'un photon).

Illustration d'un processus d'absorption
Absorption.png
Absorption d'un photon par un système matériel possédant au moins deux niveaux d'énergie. L'énergie totale est conservée dans ce processus.

Diffusion

Lorsque la présence de matière interagit avec l'onde électromagnétique en absorbant et réémettant un photon dans un temps très court, donc sans lui soutirer d'énergie. La réémission s'effectuant dans tout l'espace, il en résulte une déviation des ondes électromagnétiques hors de la direction d'incidence. On appelle ce phénomène diffusion.

Illustration d'un processus de diffusion
circle10_z.gif
Illustration de la diffusion d'une onde électromagnétique par une sphère diélectrique de rayon comparable à la longueur d'onde et d'indice de réfraction relatif n_r = 1,78. Le champ électrique total figure à gauche, et la composante diffusée à droite.
Crédit : R. Hogan (Univ. of Reading)

On regroupe les phénomènes d'absorption et de diffusion sous le nom d'extinction : ces deux phénomènes ont comme point commun de diminuer l'intensité d'un faisceau lumineux dans la direction d'incidence, soit en l'absorbant, soit en en déviant tout ou partie.


Épaisseur optique / Loi de Beer-Lambert

Profondeur optique

On définit alors la transmittance t comme le rapport entre l'intensité lumineuse I après traversée et celle I_0 avant la traversée du milieu : t(\lambda) = \frac{I(\lambda)}{I_0(\lambda)}. Cette transmittance vaut 1 si le milieu traversé est transparent, 0 si le milieu est parfaitement opaque, et prend une valeur intermédiaire si le milieu est partiellement opaque. Elle dépend a priori de la longueur d'onde comme la formule l'indique.

La transmittance a l'avantage d'être une grandeur aisément mesurable, mais son interprétation physique directe est malcommode. Il est préférable pour cela d'introduire la grandeur appelée profondeur optique \tau(\lambda) qui se déduit de la transmittance comme suit : \tau(\lambda) = - \ln \left[ t(\lambda) \right], soit I(\lambda) = I_0(\lambda) e^{-\tau(\lambda)}. Il est alors facile de déduire des propriétés mathématiques du logarithme que \tau(\lambda) est une grandeur additive le long d'un même rayon (voir figure ci-contre), et qu'elle offre donc un paramétrage naturel de l'abscisse curviligne le long de ce rayon. Il est en outre possible de relier physiquement la variation locale de \tau(\lambda) le long d'un rayon à la présence de matière traversée : c'est la loi de Beer-Lambert .

Additivité de la profondeur optique
tau_add.png
Crédit : Emmanuel Marcq CC-BY-SA

definitionLoi de Beer-Lambert

Elle s'exprime comme suit le long d'un rayon lumineux : d\tau(\lambda) = n(s) \sigma_{\mathrm{ext}}(\lambda) \, dss désigne l'abscisse curviligne le long du rayon (croissant dans le sens de propagation), n(s) la densité volumique de diffusants et/ou absorbants et \sigma_{\mathrm{ext}}(\lambda) la section efficace d'extinction des diffusants et/ou absorbants. En présence de plusieurs types de diffuseurs et/ou d'absorbeurs, il suffit d'ajouter leurs contributions individuelles. L'évolution de la profondeur optique totale \tau(\lambda) = \int_{s=-\infty}^{+\infty} d\tau(\lambda) selon la longueur d'onde au cours d'une occultation ou d'un transit offre donc la possibilité à l'observateur de déduire les propriétés matérielles du milieu traversé (densité volumique et nature physique des particules responsables de l'extinction).


Diffusions Rayleigh et Mie

Régimes de Rayleigh et de Mie

Le régime de diffusion Rayleigh concerne la diffusion de la lumière par des particules petites devant la longueur d'onde (rayon r \ll \lambda), ce quelque soit le type d'aérosol (sphérique ou non) : molécules isolées de gaz, micro gouttes, etc.

La théorie de Mie est quant à elle valable pour toutes les valeurs de r, mais uniquement pour des particules sphériques (gouttes liquides par exemple). Elle tend vers le cas Rayleigh lorsque r \ll \lambda. Pour des particules non sphériques (cristaux, aérosols de Titan, poussières), cette théorie n'est valable qu'en première approximation, et les méthodes utilisées pour aller plus loin sont hors-programme.

Section efficace de diffusion

Cette section efficace évolue de façon très différente selon le régime. Dans le régime Rayleigh, la section efficace est une fonction rapidement croissante de la longueur d'onde (une illustration courante étant la diffusion bien plus importante du bleu que de rouge dans l'atmosphère terrestre, donnant sa couleur à notre ciel). En régime de Mie, cette dépendance est bien plus faible, si bien que les diffuseurs apparaissent peu colorés en l'absence d'absorption (sur Terre, les nuages d'eau liquide sont un bon exemple).

Section efficace de diffusion
Qext_log_log_v3.png
Représentation log-log de la section efficace de diffusion en fonction de la longueur d'onde pour des particules sphériques d'un rayon r \approx 1\,\mathrm{\mu m}. Aux grandes longueurs d'onde, r \ll \lambda et on se trouve dans le régime Rayleigh (\sigma_{\mathrm{ext}} \propto \lambda^{-4}, apparaissant comme une asymptote de pente -4). Pour r \sim \lambda et r > \lambda, on se trouve dans le régime de Mie, où la section efficace est presque indépendante de \lambda et voisine de la section géométrique des particules.
Crédit : Loïc Rossi CC-BY-SA

Fonction de phase

La fonction de phase définit la répartition angulaire du rayonnement diffusé. C'est en quelque sorte la probabilité pour un rayon incident d'être diffusé dans une direction donnée.

Fonctions de phase de Rayleigh et de Mie
mie_scattering_Sharayanan_ccbysa.png
Figure montrant les indicatrices des différents processus de diffusion. La lumière incidente est supposée aller de gauche à droite. Plus la flèche est longue, plus la diffusion sera importante dans cette direction. À gauche, le régime de Rayleigh avec une diffusion relativement isotrope. Au centre, la diffusion de Mie, qui diffuse fortement vers l'avant. À droite, la diffusion de Mie pour des particules encore plus grandes.
Crédit : Sharayanan CC-BY-SA

Albédo de simple diffusion

Quand un rayonnement est diffusé et absorbé, on définit l'albédo de simple diffusion \varpi_0 comme la proportion du rayonnement qui est diffusée. Si \varpi_0 = 1, toute la lumière est diffusée, il n'y a pas d'absorption. En revanche si \varpi_0 = 0, tout la lumière est absorbée, il n'y a plus de diffusion. En général, \varpi_0 dépend de la longueur d'onde : \varpi_0(\lambda).


Surfaces

Types de réflexion

Quand du rayonnement parvient à une surface, il peut être réfléchi de plusieurs façons :

Réflexion diffuse et spéculaire
earth_specular.jpg
Photographie de la Terre prise depuis la station spatiale internationale. On distingue nettement la réflexion spéculaire sur la surface de l'océan et celle diffuse due aux nuages ou à la surface des continents.

Spectroscopie

Selon la longueur d'onde, le rayonnement incident sur une surface peut-être réfléchi mais aussi absorbé. On pourra alors étudier la réflectance d'une surface en fonction de la longueur d'onde (la réflectance valant 1 si tout la lumière est réfléchie ; 0 si tout est absorbé.).

L'étude de la réflectance d'une surface en fonction de la longueur d'onde peut être utile pour caractériser celle-ci. Ainsi certains minéraux ont un spectre avec des bandes d'absoption caractéristiques. C'est ce qui permet de caractériser des surfaces comme celles de Mars (pour détecter des minéraux hydratés par exemple), de la Lune ou bien celle d'astéroïdes et de comètes.


Température de brillance

definitionDéfinition

La température de brillance T_B(\lambda) d'un objet est la température du corps noir qui émettrait la même intensité que celle émise par l'objet à la longueur d'onde \lambda : I(\lambda) = B\left[\lambda, T_B(\lambda) \right]B\left[\lambda, T \right] désigne la fonction de Planck à la longueur d'onde \lambda pour un corps noir de température T. Comme les courbes représentant les fonctions de Planck ne se croisent jamais, il y a correspondance unique entre I(\lambda) et T_B(\lambda) : la température de brillance n'est qu'une autre façon de décrire un spectre. En particulier, elle ne représente pas forcément la température d'un objet physique.

La formule de Planck permet d'exprimer analytiquement T_B(\lambda) selon la formule T_B(\lambda) = \frac{hc}{k\lambda} \ln^{-1} \left( 1 + \frac{2hc^2}{\lambda^5 I(\lambda)} \right).

Lecture graphique de la température de brillance
nimbus-earth.png
Spectre thermique terrestre observé par l'instrument NIMBUS. Les courbes en pointillé représentent les spectres des corps noits aux températures indiquées en Kelvin. Ceci permet une lecture directe de ce spectre en température de brillance : un minimum situé autour de 215 K vers 15 µm, et un maximum un peu en dessous de 300 K entre 10 et 12 µm.
Crédit : Adapté de Hanel et al. (1970)

Intérêt

Son intérêt principal réside dans l'interprétation des spectres thermiques issus d'un objet (ici une planète et son atmosphère). En effet, sous réserve de certaines hypothèses :

La température de brillance correspond alors à la température du milieu situé à une profondeur optique égale à 1 : T_B(\lambda) \approx T\left[\tau(\lambda) = 1\right] avec \tau comptée depuis l'observateur et le long du rayon. Intuitivement, cela résulte d'un compromis dans l'intensité intégrée le long du rayon et reçue par l'observateur :


Méthodes

Auteurs: Loïc Rossi, Emmanuel Marcq

Nadir

Auteurs: Loïc Rossi, Emmanuel Marcq

Sondages thermiques

Raies en absorption vs. émission

Nous sommes à présent en mesure d'interpréter pourquoi les raies caractéristiques de certains gaz apparaissent tantôt en absorption, tantôt en émission dans les spectres thermiques observés. En effet, le centre d'une raie spectrale présente une absorption massique plus grande que les ailes de cette même raie. En conséquence, du point de vue d'un observateur extérieur à l'atmosphère, la profondeur optique unité est atteinte à des altitudes plus élevées aux centres des raies que sur leurs ailes. À partir de là, on a les cas de figure suivants :

Raies en absorption ou en émission
emis_abs.png
Profils d'une raie spectrale vue en émission (haut) ou en absorption (bas). L'altitude \tau(\lambda)=1 est atteinte plus haut au coeur de la raie (longueur d'onde \lambda_c) que dans les ailes lointaines de la raie (longueur d'onde \lambda_a) : z_c < z_a avec \tau(\lambda_c,z_c) = 1 et \tau(\lambda_a,z_a)=1.
Crédit : Emmanuel Marcq CC-BY-SA

Application

Lorsque le composé responsable d'une raie d'absorption est bien connu spectroscopiquement (son absorption en fonction de la longueur d'onde notamment) et que son profil d'abondance vertical est supposé connu au sein de l'atmosphère étudiée, il est possible de savoir à quelles altitudes sont atteintes les profondeurs optiques unité pour les différentes longueurs d'onde (ainsi, comme expliqué ci-dessus, le coeur des raies va sonder plus haut que les ailes lointaines des raies, où les rayons pourront même provenir de la surface si l'atmosphère est transparente par ailleurs). L'analyse du flux spectral reçu par l'observateur selon la longueur d'onde va alors lui permettre de reconstituer le profil thermique dans les plages d'altitudes associées aux différentes raies spectrales.

De façon semi-quantitative, lorsque la section efficace \sigma_{\mathrm{abs}}(\lambda) et le profil vertical en densité de l'absorbant n(z) sont connus, il est alors possible de trouver à quelle altitude z(\lambda) on a \tau \left[ z(\lambda) \right] = \int_{z(\lambda)}^{+\infty} \sigma_{\mathrm{abs}}(\lambda) n(z')\,dz' = 1. On sait alors qu'il règne une température T\left[ z(\lambda) \right] \approx T_B(\lambda) à cette altitude. On a supposé dans ce calcul que la visée de l'observateur pointe verticalement vers le centre de la planète (visée nadir).


Sondages de composition

Principe

Lorsque le profil vertical de température est supposé connu dans une atmosphère planétaire (par exemple au moyen de la méthode décrite précédemment), il est possible d'utiliser les raies spectrales de composés identifiés (et bien connus spectroscopiquement) pour en déduire leur profil vertical d'abondance sur une plage d'altitudes donnée.

De façon semi-quantitative et en visée nadir, on peut procéder comme suit : la connaissance du profil thermique T(z) permet, pour chaque longueur d'onde \lambda, de retrouver l'altitude z(\lambda) correspondant à la température de brillance observée — en d'autres termes telle que T_B(\lambda) = T\left[z(\lambda)\right]. Or, on sait que cette altitude correspond à la profondeur optique unité à la longueur d'onde considérée, si bien que l'on a \tau \left[ z(\lambda) \right] = \sigma_{\mathrm{abs}}(\lambda) \times N\left[z(\lambda) \right] = 1\sigma_{\mathrm{abs}}(\lambda) désigne la section efficace d'absorption du composé mesuré (supposée ici indépendante de la pression et de la température) et N(z) la densité de colonne du composé recherché au-dessus de l'altitude z, c'est-à-dire N(z) = \int_{z}^{\infty} n(z')\,dz' avec n(z') la densité volumique de l'absorbant à l'altitude z'. On peut alors en déduire par différentiation sur les altitudes sondées le profil de densité volumique n(z) du composé mesuré sur l'intervalle sondé.

Cependant, afin de moins dépendre des diverses approximations déjà évoquées, il est préférable d'utiliser un modèle de transfert de rayonnement afin de modéliser les spectres thermiques attendus et de les comparer ensuite aux observations en ajustant les différents paramètres (opération appelée fit en anglais de laboratoire). Notons que la procédure décrite ci-dessus ne fonctionne en pratique que lorsque les raies spectrales sont résolues, c'est-à-dire à une très haute résolution spectrale.

Il n'est en revanche pas possible d'inverser simultanément le profil thermique et le profil d'abondance du composé par la seule observation de ses raies d'absorption. On parle de problème dégénéré.

Absorption de quelques gaz dans l'IR thermique
GEISA2011.png
Absorption de la vapeur d'eau (rouge), du dioxyde de carbone (vert), de l'ozone (bleu) et du méthane (violet) dans les conditions standard de température et de pression. L'échelle verticale est logarithmique.
Crédit : Emmanuel Marcq (à partir de la base de données GEISA 2011) CC-BY-SA

Spectre réfléchis

Réflectance

La lumière stellaire réfléchie et reçue par un observateur extérieur a subi au moins un événement de diffusion, soit au sein de l'atmosphère (molécules de gaz, aérosols), soit par la surface pour une atmosphère peu opaque autour d'une planète tellurique. Ces propriétés de réflexion sont résumées par la fonction R(\theta_0,\varphi_0,\theta_1,\varphi_1) appelée réflectance bidirectionnelle. C'est une fonction à quatre variables, deux pour caractériser la direction du rayon incident (d'angle zénital d'incidence \theta_0 et d'azimuth \varphi_0) et deux pour celle du rayon émergent (d'angle d'émergence zénital \theta_1 et d'azimuth \varphi_1) . C'est une fonction positive, dont la valeur minimale est nulle dans les directions où l'intensité réfléchie est nulle. Elle ne se déduit pas simplement des fonctions de phase calculables en régime de Rayleigh et de Mie, notamment à cause des phénomènes de diffusion multiple et de la possible contibution de la surface.

Transmission atmosphérique

Dans le cas d'une couche réfléchissante (surface, sommet des nuages) surmontée d'une atmosphère purement absorbante (non diffusante en elle-même), on a alors une relation anayltique entre l'intensité réfléchie I(\theta_1,\varphi_1) et celle de la source I(\theta_0,\varphi_0) : I(\theta_1,\varphi_1) = I(\theta_0,\varphi_0) R(\theta_0,\varphi_0,\theta_1,\varphi_1) \exp \left[ - \tau \left(\frac{1}{\cos \theta_0} + \frac{1}{\cos \theta_1} \right) \right]\tau désigne l'épaisseur optique de l'atmosphère au-dessus de la couche réfléchissante (comptée selon la verticale). Le terme en exponentielle provient de la loi de Beer-Lambert pour le rayon incident et le rayon émergent, en tenant compte de leur inclinaison éventuelle qui les amène à traverser des épaisseurs atmosphériques plus importantes que des rayons purement verticaux.

Hormis l'étude la couche réfléchissante dans des longueurs d'onde où \tau est négligeable, cette méthode permet de mesurer le spectre de l'absorption atmosphérique décrite par \tau(\lambda) : on a ainsi accès à la densité de colonne intégrée au-dessus de la couche réfléchissante en examinant la profondeur des raies d'absorption des composés atmosphériques. Notons qu'au contraire de l'émission thermique, la présence de composés absorbants se traduit toujours par des raies en absorption et non plus en émission.

Illustration de l'intensité réfléchie
refl.png
Crédit : Emmanuel Marcq CC-BY-SA

Occultations et transits

Auteurs: Loïc Rossi, Emmanuel Marcq

Occultations

definitionLa méthode des pelures d'oignon

L'observation d'une occultation au cours du lever (ou du coucher) de la source à travers l'atmosphère d'une planète permet d'inverser le profil vertical d'extinction k_{\mathrm{ext}}(z;\lambda) = n(z) \sigma_{\mathrm{ext}}(\lambda) à la longueur d'onde \lambda à partir de l'observation des transmissions t(z_i,\lambda) selon différentes altitudes tangentesz_i (altitude minimale du rayon lors de sa traversée de l'atmosphère). Un des algorithmes permettant cette inversion est appelé méthode des pelures d'oignon : on découpe l'atmosphère en coquilles sphériques concentriques homogènes à la manière d'un oignon. L'altitude tangente z_0 la plus élevée pour laquelle on mesure \tau(z_0,\lambda) > 0 permet de déduire l'extinction locale k_{\mathrm{ext}}(z_0,\lambda) au sein de cette couche la plus extérieure. La profondeur optique observée juste en-dessous, à l'altitude tangente z_1 < z_0, est alors un peu plus grande. Comme on connaît déjà k_{\mathrm{ext}}(z_0,\lambda), on en déduit k_{\mathrm{ext}}(z_1,\lambda) pour la couche immédiatement intérieure. Et ainsi de suite jusqu'au rayon le plus bas pour lequel on puisse mesurer une transmission t(z_N,\lambda) > 0.

La mesure des profils d'extinction k_{\mathrrm{ext}}(z,\lambda) dans l'atmosphère permet alors d'en déduire plusieurs paramètres importants :

Ces occultations sont donc une méthode précieuse de sondage atmosphérique, mais assez délicate à mettre en oeuvre. Il faut en effet pouvoir mesurer de nombreuses transmissions précisément, tout en connaissant parfaitement chacune des altitudes tangentes. Ceci n'est en général possible que pour un satellite en orbite autour de la planète, et nous restreint donc au système solaire.

Méthode dite "par pelure d'oignon"
onion.png
Visualisation géométrique de la méthode d'inversion surnommée "par pelure d'oignon" : la transmission observée t(z_1) permet, connaissant l'épaisseur traversée, de remonter à k_{\mathrm{ext}}(z_1). À partir de là, la transmission t(z_2) permet, connaissant k_{\mathrm{ext}}(z_1), d'en déduire k_{\mathrm{ext}}(z_2) et ainsi de suite.
Crédit : Emmanuel Marcq CC-BY-SA

Transits primaires

Transits primaires

L'observation d'un transit primaire se traduit par une baisse relative de luminosité en provenance de la source partiellement occultée, appelée profondeur du transit, notée \delta et habituellement mesurée en ppm (parties par million). Pour une planète sans atmosphère, la profondeur est indépendante de la longueur d'onde (la réciproque n'est pas vraie : ainsi un plafond nuageux élevé peut conduire à la même observation, par exemple). En revanche, si l'atmosphère présente une extinction variable selon la longueur d'onde (que ce soit par absorption ou par diffusion Rayleigh par exemple), alors le rayon apparent de la planète (et donc la profondeur du transit) sembleront légèrement plus grands dans les longueurs d'onde où l'atmosphère cause davantage d'extinction.

On peut ainsi montrer, sous certaines hypothèses simplificatrices, que la variation relative de la profondeur du transit selon la longueur d'onde est liée à celle de la section efficace d'extinction (moyennée verticalement) au sein de l'atmosphère, selon la loi : \frac{d \delta / d\lambda}{\delta} = \frac{2H}{R} \frac{d\sigma_{\mathrm{ext}} / d \lambda}{\sigma_{\mathrm{ext}}}H désigne l'échelle de hauteur de l'atmosphère et R le rayon de la planète. Ces variations sont donc d'autant plus notables que l'échelle de hauteur est grande, ce qui favorise donc les atmosphères très dilatées. C'est par exemple le cas pour les Jupiter chauds (atmosphère légère constituée principalement de H2 et à température élevée).

Spectroscopie d'un transit primaire
transit_prim.png
Visualisation de la profondeur d'un transit primaire en fonction de la longueur d'onde. Ici, l'atmosphère est transparente aux grandes longueurs d'onde, et de plus en plus opaque aux courtes longueurs d'onde, ce qui augmente la profondeur du transit.
Crédit : Emmanuel Marcq CC-BY-SA

Transits secondaires

L'intérêt d'un transit secondaire réside dans la possibilité d'effectuer une spectroscopie différentielle de la planète : si l'on mesure un spectre S_{\mathrm{star + planet}}(\lambda) juste avant ou juste après ce transit secondaire, ainsi qu'un spectre S_{\mathrm{star}}(\lambda) de l'étoile seule pendant ce transit secondaire, on peut en déduire le spectre S_{\mathrm{planet}}(\lambda) = S_{\mathrm{star + planet}}(\lambda) - S_{\mathrm{star}}(\lambda) émis ou réfléchi par la planète seule vue à angle de phase nul.

Le spectre S_{\mathrm{planet}}(\lambda) peut alors être analysé comme le serait n'importe quel spectre (par exemple en termes de température ou de composition pour un spectre thermique, ou bien de propriétés des diffuseurs pour un spectre réfléchi). C'est du moins le cas en théorie, car en pratique le faible rapport signal à bruit de tels spectres interdit à ce jour toute analyse trop poussée de ces spectres, de résolution spectrale souvent médiocre.

Courbe de lumière
stevenson.png
Courbe de lumière lors du transit secondaire de WASP-43b. Le flux en provenance de l'étoile est normalisé à 1. Le surcroît autour du transit secondaire est causé par l'émission thermique de la planète, ce qui a permis de reconstituer son profil thermique. L'insert en haut à droite représente le transit primaire.
Crédit : Figure tirée de Stevenson et al. (2014)