TY - THES A1 - Cervantes Villa, Juan Sebastian T1 - Understanding the dynamics of radiation belt electrons by means of data assimilation T1 - Verständnis der Dynamik von Strahlungsgürtel-Elektronen durch Datenassimilation N2 - The Earth's electron radiation belts exhibit a two-zone structure, with the outer belt being highly dynamic due to the constant competition between a number of physical processes, including acceleration, loss, and transport. The flux of electrons in the outer belt can vary over several orders of magnitude, reaching levels that may disrupt satellite operations. Therefore, understanding the mechanisms that drive these variations is of high interest to the scientific community. In particular, the important role played by loss mechanisms in controlling relativistic electron dynamics has become increasingly clear in recent years. It is now widely accepted that radiation belt electrons can be lost either by precipitation into the atmosphere or by transport across the magnetopause, called magnetopause shadowing. Precipitation of electrons occurs due to pitch-angle scattering by resonant interaction with various types of waves, including whistler mode chorus, plasmaspheric hiss, and electromagnetic ion cyclotron waves. In addition, the compression of the magnetopause due to increases in solar wind dynamic pressure can substantially deplete electrons at high L shells where they find themselves in open drift paths, whereas electrons at low L shells can be lost through outward radial diffusion. Nevertheless, the role played by each physical process during electron flux dropouts still remains a fundamental puzzle. Differentiation between these processes and quantification of their relative contributions to the evolution of radiation belt electrons requires high-resolution profiles of phase space density (PSD). However, such profiles of PSD are difficult to obtain due to restrictions of spacecraft observations to a single measurement in space and time, which is also compounded by the inaccuracy of instruments. Data assimilation techniques aim to blend incomplete and inaccurate spaceborne data with physics-based models in an optimal way. In the Earth's radiation belts, it is used to reconstruct the entire radial profile of electron PSD, and it has become an increasingly important tool in validating our current understanding of radiation belt dynamics, identifying new physical processes, and predicting the near-Earth hazardous radiation environment. In this study, sparse measurements from Van Allen Probes A and B and Geostationary Operational Environmental Satellites (GOES) 13 and 15 are assimilated into the three-dimensional Versatile Electron Radiation Belt (VERB-3D) diffusion model, by means of a split-operator Kalman filter over a four-year period from 01 October 2012 to 01 October 2016. In comparison to previous works, the 3D model accounts for more physical processes, namely mixed pitch angle-energy diffusion, scattering by EMIC waves, and magnetopause shadowing. It is shown how data assimilation, by means of the innovation vector (the residual between observations and model forecast), can be used to account for missing physics in the model. This method is used to identify the radial distances from the Earth and the geomagnetic conditions where the model is inconsistent with the measured PSD for different values of the adiabatic invariants mu and K. As a result, the Kalman filter adjusts the predictions in order to match the observations, and this is interpreted as evidence of where and when additional source or loss processes are active. Furthermore, two distinct loss mechanisms responsible for the rapid dropouts of radiation belt electrons are investigated: EMIC wave-induced scattering and magnetopause shadowing. The innovation vector is inspected for values of the invariant mu ranging from 300 to 3000 MeV/G, and a statistical analysis is performed to quantitatively assess the effect of both processes as a function of various geomagnetic indices, solar wind parameters, and radial distance from the Earth. The results of this work are in agreement with previous studies that demonstrated the energy dependence of these two mechanisms. EMIC wave scattering dominates loss at lower L shells and it may amount to between 10%/hr to 30%/hr of the maximum value of PSD over all L shells for fixed first and second adiabatic invariants. On the other hand, magnetopause shadowing is found to deplete electrons across all energies, mostly at higher L shells, resulting in loss from 50%/hr to 70%/hr of the maximum PSD. Nevertheless, during times of enhanced geomagnetic activity, both processes can operate beyond such location and encompass the entire outer radiation belt. The results of this study are two-fold. Firstly, it demonstrates that the 3D data assimilative code provides a comprehensive picture of the radiation belts and is an important step toward performing reanalysis using observations from current and future missions. Secondly, it achieves a better understanding and provides critical clues of the dominant loss mechanisms responsible for the rapid dropouts of electrons at different locations over the outer radiation belt. N2 - Die Elektronenstrahlungsgürtel der Erde weisen eine Zwei-Zonen-Struktur auf, wobei der äußere Gürtel aufgrund des ständigen Zusammenspiels zwischen einer Reihe von physikalischen Prozessen, einschließlich Beschleunigung, Verlust und Transport, eine hohe Dynamik aufweist. Der Elektronenfluss im äußeren Gürtel kann über mehrere Größenordnungen variieren und Werte erreichen, die den Satellitenbetrieb stören können. Daher ist das Verständnis der Mechanismen, die diese Variabilität bewirken, von hohem Interesse für die wissenschaftliche Gemeinschaft. Insbesondere die wichtige Rolle die Verlustmechanismen bei der Kontrolle der relativistischen Elektronendynamik spielen ist in den letzten Jahren immer deutlicher geworden. Es ist inzwischen weithin anerkannt, dass Strahlungsgürtelelektronen entweder durch Interaktion mit der Atmosphäre oder durch Transport über die Magnetopause, das so genannte Magnetopauseshadowing, verloren gehen können. Der Verlust von Elektronen in der Atmosphäre erfolgt aufgrund von Pitchwinkelstreuung durch resonante Wechselwirkung mit verschiedenen Arten von magnetosphärischen Wellen, einschließlich plasmasphärischem Hiss, Whistler-Mode-Chorus, und elektromagnetischen Ionenzyklotron-Wellen (EMIC). Darüber hinaus kann die Komprimierung der Magnetopause aufgrund der Erhöhungen des dynamischen Drucks des Sonnenwindes dazu führen, dass Elektronen an hohen L-Shells, wo sie sich in offenen Driftpfaden befinden, erheblich in ihrer Dichte reduziert werden, während Elektronen an niedrigen L-Shells durch radiale Diffusion nach außen verloren gehen können. Nichtsdestotrotz bleibt die Rolle, die jeder physikalische Prozess bei der schnellen Reduktion des Elektronenflusses spielt, nach wie vor ein grundlegendes Rätsel. Die Unterscheidung zwischen diesen Prozessen und die Quantifizierung ihrer relativen Beiträge zur Entwicklung der Strahlungsgürtelelektronen erfordert hochauflösende Profile der Phasenraumdichte (PSD). Solche Profile der PSD sind jedoch schwierig zu bestimmen, da die Beobachtungen von Raumfahrzeugen auf eine einzige Messung in Raum und Zeit beschränkt sind, was auch durch die Ungenauigkeit der Instrumente erschwert wird. Datenassimilationstechniken zielen darauf ab, unvollständige und ungenaue raumgestützte Daten mit physikalisch basierten Modellen auf optimale Weise zu kombinieren. In den Strahlungsgürteln der Erde werden sie verwendet, um das gesamte radiale Profil der Elektronen-PSD zu rekonstruieren, und sie sind zu einem immer wichtigeren Werkzeug geworden, um unser derzeitiges Verständnis der Dynamik der Strahlungsgürtel zu validieren, neue physikalische Prozesse zu identifizieren und die erdnahe gefährliche Strahlungsumgebung vorherzusagen. In dieser Studie werden Messungen der Van-Allen-Probes A und B und der Geostationary-Operational-Environmental-Satellites (GOES) 13 und 15 mit Hilfe eines Split-Operator-Kalman-Filters über einen Zeitraum von vier Jahren vom 01. Oktober 2012 bis zum 01. Oktober 2016 in das dreidimensionale Versatile Electron Radiation Belt-3D-Diffusionsmodell (VERB-3D) integriert. Im Vergleich zu früheren Arbeiten berücksichtigt das 3D-Modell mehr physikalische Prozesse, nämlich gemischte Diffusion, Streuung durch EMIC-Wellen und Magnetopausenverluste. Es wird gezeigt, wie die Datenassimilation mit Hilfe des Innovationsvektors (des Residuums zwischen Beobachtungen und Modellprognose), genutzt werden kann, um fehlende physikalische Prozesse im Modell zu berücksichtigen. Diese Methode wird verwendet, um die radialen Entfernungen von der Erde und die geomagnetischen Bedingungen zu identifizieren, bei denen unser Modell für verschiedene Werte der adiabatischen Invarianten mu und K nicht mit der gemessenen PSD übereinstimmt. Infolgedessen passt der Kalman-Filter die Vorhersagen an die Beobachtungen an, und dies wird als Nachweis dafür interpretiert, wo und wann zusätzliche Quellen- oder Verlustprozesse aktiv sind. Darüber hinaus werden zwei unterschiedliche Verlustmechanismen untersucht, die für die schnellen Verluste von Strahlungsgürtelelektronen verantwortlich sind: EMIC-Wellen-induzierte Streuung und Magnetopausenverluste. Der Innovationsvektor wird bei Werten der Invariante mu im Bereich von 300 bis 3000 MeV/G untersucht, und es wird eine statistische Analyse durchgeführt, um die Wirkung beider Prozesse in Abhängigkeit von verschiedenen geomagnetischen Indizes, Sonnenwindparametern und der radialen Entfernung von der Erde quantitativ zu bewerten. Die Ergebnisse dieser Arbeit stehen in Übereinstimmung mit früheren Studien, die die Energieabhängigkeit dieser beiden Mechanismen nachgewiesen haben. Die EMIC-Wellenstreuung dominiert den Verlust bei niedrigen L-Shells und kann zwi-schen 10%/hr bis 30%/hr des Maximalwertes der PSD über alle L-Shells für feste Werte der ersten und zweiten adiabatische Invarianten betragen. Andererseits wird festgestellt, dass bei den Magnetopausenverlusten über alle Energien hinweg, meist bei höheren L-Shells, Elektronen Verluste zeigen, was zu einer Verstärkung des Verlustes von 50%/hr auf 70%/hr der maximalen PSD führt. Nichtsdestotrotz können beide Prozesse in Zeiten erhöhter geomagnetischer Aktivität über diese L-Shells hinaus wirken und den gesamten äußeren Strahlungsgürtel umfassen. Die Ergebnisse dieser Studie sind zweifacher Art. Erstens zeigt sie, dass der 3D-Daten-Assimilationscode ein umfassendes Bild der Strahlungsgürtel liefert und ein wichtiger Schritt zur Durchführung einer Reanalyse unter Verwendung von Beobachtungen aus aktuellen und zukünftigen Missionen ist. Zweitens erreicht er ein besseres Verständnis und liefert entscheidende Hinweise auf die vorherrschenden Verlustmechanismen, die für die schnellen Verluste von Elektronen an verschiedenen Orten im äußeren Strahlungsgürtel verantwortlich sind. KW - radiation belts KW - Strahlungsgürtel KW - data assimilation KW - Datenassimilation KW - phase space density KW - Phasenraumdichte KW - magnetospheric waves KW - magnetosphärischen Wellen KW - Kalman filter KW - Kalman-Filter Y1 - 2021 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-519827 ER - TY - THES A1 - Aseev, Nikita T1 - Modeling and understanding dynamics of charged particles in the Earth's inner magnetosphere T1 - Modellierung und Untersuchung der Dynamik geladener Teilchen in der inneren Magnetosphäre der Erde N2 - The Earth's inner magnetosphere is a very dynamic system, mostly driven by the external solar wind forcing exerted upon the magnetic field of our planet. Disturbances in the solar wind, such as coronal mass ejections and co-rotating interaction regions, cause geomagnetic storms, which lead to prominent changes in charged particle populations of the inner magnetosphere - the plasmasphere, ring current, and radiation belts. Satellites operating in the regions of elevated energetic and relativistic electron fluxes can be damaged by deep dielectric or surface charging during severe space weather events. Predicting the dynamics of the charged particles and mitigating their effects on the infrastructure is of particular importance, due to our increasing reliance on space technologies. The dynamics of particles in the plasmasphere, ring current, and radiation belts are strongly coupled by means of collisions and collisionless interactions with electromagnetic fields induced by the motion of charged particles. Multidimensional numerical models simplify the treatment of transport, acceleration, and loss processes of these particles, and allow us to predict how the near-Earth space environment responds to solar storms. The models inevitably rely on a number of simplifications and assumptions that affect model accuracy and complicate the interpretation of the results. In this dissertation, we quantify the processes that control electron dynamics in the inner magnetosphere, paying particular attention to the uncertainties of the employed numerical codes and tools. We use a set of convenient analytical solutions for advection and diffusion equations to test the accuracy and stability of the four-dimensional Versatile Electron Radiation Belt (VERB-4D) code. We show that numerical schemes implemented in the code converge to the analytical solutions and that the VERB-4D code demonstrates stable behavior independent of the assumed time step. The order of the numerical scheme for the convection equation is demonstrated to affect results of ring current and radiation belt simulations, and it is crucially important to use high-order numerical schemes to decrease numerical errors in the model. Using the thoroughly tested VERB-4D code, we model the dynamics of the ring current electrons during the 17 March 2013 storm. The discrepancies between the model and observations above 4.5 Earth's radii can be explained by uncertainties in the outer boundary conditions. Simulation results indicate that the electrons were transported from the geostationary orbit towards the Earth by the global-scale electric and magnetic fields. We investigate how simulation results depend on the input models and parameters. The model is shown to be particularly sensitive to the global electric field and electron lifetimes below 4.5 Earth's radii. The effects of radial diffusion and subauroral polarization streams are also quantified. We developed a data-assimilative code that blends together a convection model of energetic electron transport and loss and Van Allen Probes satellite data by means of the Kalman filter. We show that the Kalman filter can correct model uncertainties in the convection electric field, electron lifetimes, and boundary conditions. It is also demonstrated how the innovation vector - the difference between observations and model prediction - can be used to identify physical processes missing in the model of energetic electron dynamics. We computed radial profiles of phase space density of ultrarelativistic electrons, using Van Allen Probes measurements. We analyze the shape of the profiles during geomagnetically quiet and disturbed times and show that the formation of new local minimums in the radial profiles coincides with the ground observations of electromagnetic ion-cyclotron (EMIC) waves. This correlation indicates that EMIC waves are responsible for the loss of ultrarelativistic electrons from the heart of the outer radiation belt into the Earth's atmosphere. N2 - Die innere Magnetosphäre der Erde ist ein sehr dynamisches System, das hauptsächlich vom äußeren Sonnenwind beeinflusst wird, der auf das Magnetfeld unseres Planeten einwirkt. Störungen im Sonnenwind, wie z.B. koronale Massenauswürfe und sogenannte Korotierende Wechselwirkungsbereiche, verursachen geomagnetische Stürme, die zu deutlichen Veränderungen der Populationen geladener Teilchen in der inneren Magnetosphäre führen - Plasmasphäre, Ringstrom und Strahlungsgürtel. Satelliten, die in Regionen mit erhöhten energetischen und relativistischen Elektronenflüssen betrieben werden, können durch tiefe dielektrische Ladung oder Oberflächenladungen bei schweren Weltraumwetterereignissen beschädigt werden. Die Vorhersage der Dynamik der geladenen Teilchen und die Abschwächung ihrer Auswirkungen auf die Infrastruktur sind heutzutage von besonderer Bedeutung, insbesondere aufgrund unserer zunehmenden Abhängigkeit von Weltraumtechnologien. Die Dynamik von Teilchen in der Plasmasphäre, des Ringstrom und in den Strahlungsgürteln sind durch Kollisionen und kollisionsfreie Wechselwirkungen mit elektromagnetischen Feldern, die durch die Bewegung geladener Teilchen induziert werden, stark gekoppelt. Mehrdimensionale numerische Modelle vereinfachen die Betrachtung von Transport-, Beschleunigungs- und Verlustprozessen dieser Partikel und ermöglichen es uns, vorherzusagen, wie die erdnahe Weltraumumgebung auf Sonnenstürme reagiert. Die Modelle beruhen zwangsläufig auf einer Reihe von Vereinfachungen und Voraussetzungen, die sich auf die Modellgenauigkeit auswirken und die Interpretation der Ergebnisse erschweren. In dieser Dissertation quantifizieren wir die Prozesse, die die Dynamik der Elektronen in der inneren Magnetosphäre steuern. Dabei richten wir den Fokus insbesondere auch auf die Unsicherheiten der verwendeten numerischen Codes. Wir verwenden eine Reihe praktischer analytischer Lösungen für Advektions- und Diffusionsgleichungen, um die Genauigkeit und Stabilität des 4-dimensionalen ''Versatile Electron Radiation Belt'' Codes (VERB-4D Code) zu testen. Wir zeigen, dass die im Code implementierten numerischen Schemata zu den analytischen Lösungen konvergieren und der Code sich unabhängig vom angenommenen Zeitschritt stabil verhält. Wir demonstrieren, wie die Genauigkeit des numerischen Schemas für die Konvektionsgleichung die Ergebnisse von Ringstrom- und Strahlungsgürtelsimulationen beeinflussen kann, und dass es von entscheidender Beteutung ist, numerische Schemata höherer Ordnung zu verwenden, um numerische Fehler im Modell zu reduzieren. Mit dem ausführlich getesteten VERB-4D Code modellieren wir die Dynamik der Ringstromelektronen während des Sturms vom 17. März 2013. Wir zeigen, dass die Diskrepanzen zwischen dem Modell und Beobachtungen oberhalb von 4.5 Erdradien durch Unsicherheiten in den äußeren Randbedingungen erklärt werden können und dass die Elektronen durch die globalen elektrischen und magnetischen Felder von der geostationäre Umlaufbahn zur Erde transportiert wurden. Wir untersuchen weiterhin, wie die Simulationsergebnisse von den Eingabemodellen und Parametern abhängen. Wir zeigen, dass das Modell besonders empfindlich für das globale elektrische Feld und die Lebensdauer der Elektronen unterhalb von 4.5 Erdradien ist. Außerdem quantifizieren wir auch die Auswirkungen von radialer Diffusion und subauroralen Polarisationsströmen. Wir haben einen datenassimilativen Code entwickelt, der mithilfe des Kalman-Filters ein Konvektionsmodell für den Transport und den Verlust energetischer Elektronen mit den Satellitendaten der Van Allen Probes kombiniert. Wir zeigen, dass die Verwendung eines Kalman-Filters Modellunsicherheiten im elektrischen Konvektionsfeld, in der Lebensdauer der Elektronen und in den Randbedingungen korrigieren kann. Weiterhin zeigen wir, wie der Innovationsvektor - die Differenz zwischen Beobachtungen und Modellvorhersagen - verwendet werden kann, um physikalische Prozesse zu identifizieren, die im Modell der Dynamik der energetischen Elektronen fehlen. Außerdem berechnen wir radiale Profile der Phasenraumdichte ultrarelativistischer Elektronen mithilfe von Van Allen Probes-Messungen. Wir analysieren die Form der Profile und zeigen, dass die Entstehung neuer lokaler Minima in den radialen Profilen mit den Bodenbeobachtungen von EMIC-Wellen übereinstimmt. Diese Korrelation legt nahe, dass EMIC-Wellen für den Verlust ultrarelativistischer Elektronen vom Herzen des äußeren Strahlungsgürtels in die Erdatmosphäre verantwortlich sind. KW - ring current electrons KW - radiation belts KW - mathematical modeling KW - wave-particle interactions KW - data assimilation KW - Ringstromelektronen KW - Strahlungsgürtel KW - mathematische Modellierung KW - Wellen-Teilchen Wechselwirkungen KW - Datenassimilation Y1 - 2020 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-479211 ER -