Les ceintures de radiation correspondent à la région de la magnétosphère dans laquelle se trouvent les particules de hautes énergies. Le couplage entre le vent solaire et la magnétosphère donne lieu à des variations des flux de particules sur plusieurs ordres de grandeurs. L’objectif de cette thèse est d’observer et caractériser ces variations de flux d’électrons au passage de différents types d’événements tels que les régions d’interaction en co-rotation (CIRs) et les éjections de masse coronale interplanétaires (ICMEs). Pour cela, nous avons traité et analysé les données de plusieurs types: paramètres du vent solaire, indices géomagnétiques et flux d’électrons dans les ceintures de radiation. Dans les trois premiers chapitres, nous rendons compte de la complexité de l’environnement spatial Terrestre et présentons les différentes données utilisées. Les travaux de thèse sont ensuite organisés en quatre chapitres. Premièrement, nous utilisons les mesures des satellites NOAA-POES afin de caractériser les flux d’électrons dans les ceintures. Nous étudions ensuite les différences de variations de flux causées par les CIRs et les ICMEs en fonction de l’énergie des électrons et du paramètre L*. Après avoir montré le fort lien entre les intensités d’orages magnétiques et les variations de flux, nous nous focalisons sur les ICMEs et la variabilité des orages qu’elles causent. Enfin, nous insistons sur l’importance des enchaînements d’événements. Après avoir quantifié la forte tendance qu’ont les ICMEs à former des séquences, nous réalisons une étude statistique sur les orages qu’elles causent, puis trois études de cas afin d’illustrer leurs effets sur les ceintures.

Depuis le début de l’ère spatiale avec le lancement du satellite Spoutnik 1 en 1957, les ceintures de radiation terrestres n’ont cessé de faire l’objet d’études du fait de leur dangerosité pour les satellites mais aussi pour l’être humain. En effet, lors d’une forte activité solaire, l’injection de particules dans cet environnement radiatif peut induire des flux jusqu’à 1000 fois plus élevés que par temps calme. Par conséquent, il est important d’en comprendre la physique ainsi que la dynamique au cours de ce que l’on appelle un orage géomagnétique. Dans ce but, le Département Physique Instrumentation Environnement et Espace (DPhIEE) de l’ONERA développe depuis maintenant plus de 20 ans la famille de modèles Salammbô reproduisant de façon robuste et en trois dimensions la dynamique des particules piégées dans ces ceintures. Néanmoins, bien que précis au-delà d’environ 100 keV, la physique et les hypothèses prises en compte dans ce modèle restent insuffisantes en deçà. En effet, aux basses énergies, les ceintures de radiation ne peuvent plus être considérées comme homogènes autour de la Terre. L’objectif de cette thèse a donc été de prendre en compte une quatrième dimension, le temps magnétique local (MLT), afin de mieux reproduire l’évolution des structures fines lors d’un orage géomagnétique. La première partie s’est portée sur l’optimisation du schéma numérique. L’ajout d’une quatrième dimension induit, via l’apparition d’un terme d’advection, une forte diffusion numérique qu’il convient de limiter, tout en tenant compte du temps de calcul. L’équation statistique implémentée a alors été discrétisée selon un schéma de type Beam-Warming du second ordre couplé à un limiteur Superbee, garantissant une propagation satisfaisante de la distribution initiale. Une fois les problèmes numériques maitrisés, les différents mécanismes physiques pilotant la dynamique des particules piégées ont été implémentés dans le code, avec une attention toute particulière sur la dépendance en MLT de l’interaction onde-particule. La prise en compte des champs électriques magnétosphériques fut également nécessaire. En effet, ils constituent l’un des moteurs principaux du mouvement des particules de basses énergies. Le modèle Salammbô 4D a ensuite été validé par comparaison avec le modèle 3D déjà existant sur une simulation de l’orage magnétique de Mars 2015. Les résultats ont montré une bonne restitution de la dynamique des ceintures de radiation, avec en plus l’accès à la phase principale de l’orage. Cet évènement a ensuite été modélisé à plus basse énergie pour constater la dynamique asymétrique des électrons piégés avec le rôle prépondérant du champ électrique de convection. La comparaison avec les données du satellite THEMIS a montré une bonne modélisation des différents processus physiques, notamment celui de « dropout » par traversée de la magnétopause. Enfin, la mise en place d’une condition limite dynamique modulée par les paramètres du vent solaire et dépendante du MLT ouvre de nombreuses perspectives.

Cette étude s’inscrit dans le domaine de la description de la dynamique des ceintures de radiation terrestres. Elle consiste à modéliser le phénomène de diffusion radiale en travaillant avec une résolution spatio-temporelle plus fine que celle utilisée pour décrire la dynamique des ceintures par le biais d’une équation de diffusion. La démarche s’est organisée en trois temps. Tout d’abord, l’objectif a été d’étudier le phénomène de diffusion radiale d’un point de vue théorique afin de mettre en lumière les principaux pilotes du processus et d’expliciter une formulation des coefficients de diffusion radiale. Une fois l’expression de ces coefficients établie, l’objectif a ensuite été de les quantifier. Pour cela, nous avons développé des protocoles analytiques et numériques puis des protocoles expérimentaux. Nous avons discuté les résultats obtenus ainsi que les atouts et les limites de ces protocoles. Cette étude met en évidence le rôle central de l’asymétrie des variations du champ électromagnétique et des champs électriques induits dans le processus de diffusion radiale. Elle propose des pistes pour la quantification numérique et expérimentale de ces deux pilotes. Elle apporte également un regard critique sur les travaux de la littérature. Elle ouvre la voie pour une nouvelle quantification des coefficients de diffusion basée sur une modélisation adéquate de la dynamique de l’environnement électromagnétique

Les ceintures de radiation provoquent des dégâts irréversibles sur les satellites les traversant, détériorant ainsi les instruments de mesure embarqués. Les étudier est utile au développement de matériaux adaptés et résistants. Depuis les années 90, l'ONERA-DESP étudie les ceintures de radiations des planètes magnétisées, telle que la Terre ou Jupiter, grâce au modèle Salammbô. Salammbô prend en compte les processus physiques de l'environnement radiatif pour recréer les populations d'électrons peuplant les ceintures. Dans cette étude, il s'agit de développer un modèle des ceintures de radiations internes de Saturne, basé sur les travaux précédents. Avant les années 2000, Pioneer 11 et Voyager 2 ne permettaient pas un développement suffisamment avancé d'un modèle de ceintures de radiations de Saturne. La mission Cassini apporte ensuite quantités d'observations et de données pour mener une étude plus approfondie de ces ceintures. Cette thèse débute par l'analyse de la magnétosphère interne de Saturne : anneaux, satellites, nuages de neutres...L'interaction des particules des ceintures de radiations avec ces différents paramètres se traduit par le calcul de coefficients de diffusion. Ces coefficients sont intégrés à l'équation de transport et permettent de comprendre les mécanismes perturbant la distribution des électrons au sein des ceintures de radiations. Les résultats obtenus sont comparés aux mesures faites par les sondes Pioneer 11, Voyager 2 et Cassini.

Le nombre croissant de charges utiles mises à poste dans l’espace et l’utilisation de systèmes électroniques embarqués de pointe nécessitent une rigoureuse connaissance de l’environnement spatial de notre planète en vue des choix technologiques à poser. Afin de mieux connaître l’environnement radiatif, une étude du transport de plasma dans la région interne de la magnétosphère terrestre est proposée ici, sachant que ce sont ces particules énergétiques qui vont venir peupler les ceintures de radiation. Dans le cadre de la modélisation, on s’intéresse à la source du plasma, son transport et les pertes de particules rencontrées au cours de ce dernier. Une étude paramétrique comparative de modèles de champs magnétique et électrique est ainsi menée, dévoilant l’influence capitale du champ électrique sur la dynamique des particules. Le même type d’étude est entrepris avec des modèles de la densité d’hydrogène dans l’exosphère, qui permettront par la suite d’estimer la perte de protons par échange de charge au cours de leur transport. Ces résultats ont motivé une étude plus approfondie du champ électrique sur la base des mesures proton et électron des détecteurs MEPED des satellites NOAA POES, en montrant la pénétration de particules à très faibles valeurs de L lors d’orages géomagnétiques. L’étude d’événements particuliers à différents niveaux d’activité magnétique et la simulation du phénomène permettent d’arriver à une caractérisation du champ électrique de convection dans la magnétosphère interne.

La composante ionisante des rayonnements spatiaux s'appelle dose cumulée. Elle entraîne l'apparition d'une charge piégée dans les oxydes ainsi que des états d'interface. Ces quantités sont à l'origine des dégradations électriques observées dans les circuits. Un modèle numérique des effets de dose dans la silice est développé. Il décrit la génération, le transport, le piégeage, les phénomènes de guérison ainsi que la génération des états d'interface tout en s'attachant à prendre à compte de manière assez fine l'effet du champ électrique et de la température. Le calage de notre modèle avec des données expérimentales faites sur la gamme des débits de dose de laboratoire avec différentes températures et diverses conditions de polarisations a permis de caractériser les paramètres des différents mécanismes. Sous faible champ électrique, certains composants présentent une sensibilité accrue au faible débit de dose. Notre modèle explique ce phénomène par l'inversion du champ électrique à fort débit de dose. Enfin, l'extrapolation au débit de dose spatial de notre modèle a permis de discuter la représentativité des normes de test en vigueur.

Les ceintures de radiation de la Terre ont été mises en évidence au début de l’ère spatiale en 1958 par J. Van Allen grâce aux données du satellite Explorer 3. Cet environnement peut être très hostile aussi bien pour le corps humain que pour l’électronique embarquée sur les satellites. Les modèles actuels de référence sont ceux de la NASA : AE8 pour les électrons et AP8 pour les protons. Ils ont été élaborés vers la fin des années 70 et le début des années 80 à l’aide de mesures éparses et d’une interpolation entre chacune d’elles. Parallèlement à cela, la modélisation physique des ceintures de radiation s’est accélérée avec l’essor de l’informatique et des mesures satellites in situ. L’ONERA / DESP dispose à la fois d’un modèle physique dynamique performant des ceintures de radiation et d’un large éventail de mesures satellites. L’objectif de ce travail de thèse a été de tirer parti de cette dualité que peu de laboratoires possèdent grâce à l’assimilation de données. Elle peut être considérée comme un processus d’association d’observations ponctuelles à une prédiction globale réalisée par un modèle. Deux méthodes ont été envisagées aussi bien dans le cas électrons que protons. La première a consisté à poursuivre les travaux d’assimilation directe entamés au département DESP afin de l’adapter à l’échelle du cycle solaire et ainsi poser des bases solides à la définition de nouveaux modèles de spécification de l’environnement spatial. La seconde est l’adaptation d’un filtre de Kalman d’Ensemble (EnKF) issu du domaine de l’océanographie aux ceintures de radiation. Technique d’assimilation évoluée, elle permet d’affiner la restitution et ouvre la voie à la prédiction et à une météorologie spatiale objective.