Equatorial plasma bubbles (EPBs) are large structures consisting of depleted plasma that generally form on the nightside of Earth's ionosphere along magnetic field lines in the upper thermosphere/ionosphere. While referred to as `bubbles', EPBs tend to be longer along magnetic latitudes and narrower along magnetic longitudes which are on the order of thousands and hundreds of kilometers, respectively. EPBs are a well documented occurrence with observations spanning many decades. As such, much is known about their general behavior, seasonal variation of occurrences, increasing/decreasing occurrences with increasing/decreasing solar activity, and their ability to interact and interfere with radio waves such as GPS. This dissertation expands on this understanding by focusing on the detection and tracking of EPBs in the upper thermosphere/ionosphere along equatorial to low latitudes. To do this, far ultraviolet (FUV) emission observations of the recombination of O$^+$ with electrons via the Global-Scale Observations of the Limb and Disk (GOLD) mission are analyzed. GOLD provides consistent data from geostationary orbit with the eastern region of the Americas, Atlantic, and western Africa. The optical data can be used to pick out gradients in brightness along the 135.6 nm wavelength which correlate with the location of EPBs in the nightside ionosphere. The dissertation provides a novel method to look at and analyze 2-dimensional data with inconsistent time-steps for EPB detection and tracking. During development, preprocessing of large scale (multiple years) data proved to be the largest time sync. To that end, this dissertation tests the possibility of using convolution neural networks for detection of EPBs with the end goal of reducing the amount of preprocessing necessary. Further, data from the Ionospheric Connection Explorer's (ICON's) ion velocity meter (IVM) are compared to EPBs detected via GOLD to understand how the ambient plasma around the EPBs behave. Along with the ambient plasma, zonal and meridional thermospheric winds observed by ICON's Michelson Interferometer for Global High-resolution Thermospheric Imaging (MIGHTI) instrument are analyzed in conjunction with the same EPBs to understand how winds coincident with EPBs behave. An analysis of winds before EPBs form is also done to observe the potential for both zonal and meridional winds' ability to suppress and amplify EPB formation. / Doctor of Philosophy / Equatorial plasma bubbles (EPBs) are large structures that generally form during post- sunset along Earth's magnetic equator. While referred to as `bubbles', EPBs tend to be thousands of kilometers from north to south and hundreds of kilometers from east to west and well over a thousands kilometers in altitude. EPBs are a well documented occurrence with observations spanning many decades. This includes their ability to interfere with radar and GPS. This dissertation expands on the scientific community's understanding by focusing on the detection and tracking of EPBs along the magnetic equator. To do this, observations from the NASA Global-Scale Observations of the Limb and Disk (GOLD) mission are analyzed. GOLD provides consistent observations looking over the eastern region of the Americas, Atlantic, and western Africa. A unique method to look at and analyze this data for EPB detection and tracking is developed. This dissertation also tests the possibility of using machine learning for detection of EPBs. Further, data from the NASA Ionospheric Connection Explorer (ICON) mission is compared to EPBs detected via GOLD to understand how the behavior of the upper atmosphere and the conductive region therein, known as the ionosphere, interact with the EBPs themselves.

The L5 signal of the Global Positioning System (GPS) is becoming available on an increasing number of Block IIF satellites. As the third civilian signal, L5 is superior in signal design to the L1 C/A and L2C civilian signals. This new signal has been marked healthy for use on selected satellites since 2010, yet the hardware capable of tracking the L5 signal is still in the early stages of development. This work investigates the characteristics of the new signal and the quality of data produced by L5-tracking receivers. Commonly used receiver models chosen for this study are the Leica GRX1200+GNSS, the Trimble NetR8, and the Javad Delta TRE-G3TH. The metrics used in this analysis to assess the quality of data produced by these receivers are signal strength, receiver phase noise, receiver code noise, and multipath. The data used in these analyses were obtained from the International GNSS Service for the days of the year 275 to 281 in 2011. Metrics averaged over the GPS week 1656 provide a good indication of the overall performance of the receivers. / text

Le PPP (Precise Point Positioning) est une méthode GNSS (Global Navigation Satellite Systems), basée sur le concept SSR (State Space Representation). Grâce aux améliorations récentes des modèles atmosphériques, le PPP en temps réel (RT-PPP) peut être également amélioré. L'objectif principal de ce travail est d'étudier le RT-PPP et l'infrastructure optimisée en termes de coûts et d'avantages pour réaliser la méthode en utilisant des corrections atmosphériques. Pour cela, différentes configurations d'un réseau GNSS dense et régulier existant en France, le réseau Orphéon, sont utilisées. Ce réseau compte environ 160 sites, propriété de Geodata-Diffusion (Hexagon Geosystems). Dans un premier temps, le mode «PPP-RTK flottant» a été évalué, il correspond au RT-PPP avec des améliorations issues des corrections de réseau, mais avec les ambiguïtés flottantes. Ensuite, des corrections de réseau sont appliquées pour améliorer le mode « PPP-RTK » où les ambiguïtés sont fixées à leurs valeurs entières. Pour le PPP-RTK flottant, une version modifiée du package RTKLib 2.4.3 (beta) est utilisée pour prendre en compte les corrections réseau. Les effets ionosphériques de premier ordre ont été éliminés par la combinaison iono-free et le retard troposphérique zénithal est estimé. Les corrections ont été appliquées en introduisant des paramètres troposphériques a priori contraints. Une modélisation adaptative basée sur les OFCs (Optimal Fitting Coefficients) a été mise en place pour décrire le comportement de la troposphère, en utilisant des estimations des retards troposphériques pour les stations Orphéon. Cette solution permet une communication monodirectionnelle entre le serveur et l'utilisateur. Les gains réalisés sur le temps de convergence pour obtenir un positionnement de 10 centimètres de précision ont été quantifiés statistiquement. La topologie du réseau a été évaluée, en réduisant le nombre de stations de référence (jusqu'à 75%), via une configuration de réseau lâche. Dans la deuxième étape, le PPP-RTK est réalisé grâce au logiciel PPP-Wizard 1.3 et avec les produits temps réel CNES (Centre Nacional de Estudes Spatiales) pour les orbites, les horloges et les biais de phase des satellites. Le RT-IPPP (RT-Integer PPP) est réalisé avec estimation des délais troposphériques et ionosphériques. Les corrections ionosphériques et troposphériques sont introduites en tant que paramètres a priori contraints au PPP-RTK. Pour générer des corrections ionosphériques, il a été mis en place un algorithme d'interpolation à distance inversée (IDW–Inverse Distance Weighting). Les améliorations apportées au positionnement horizontal dues aux corrections atmosphériques SSR externes provenant d’un réseau (dense ou lâche) sont prometteuses et peuvent être utiles pour les applications qui dépendent principalement du positionnement horizontal. / PPP (Precise Point Positioning) is a GNSS (Global Navigation Satellite Systems) method, based on SSR (State Space Representation) concept. Thanks to recent improvements in atmospheric models, Real-time PPP (RT-PPP) can also be improved. The main objective of this work is to study the RT-PPP and the optimized infrastructure in terms of costs and benefits to realize the method using atmospheric corrections. Therefore, different configurations of a dense and regular GNSS network existing in France, the Orpheon network, are used. This network has about 160 sites and is owned by Geodata-Diffusion (Hexagon Geosystems). Initially, ‘float PPP-RTK’ was evaluated, it corresponds to RT-PPP with improvements resulting from network corrections, although with ambiguities kept float. Further on, network corrections are applied to improve “PPP-RTK” where ambiguities are fixed to their integer values. For the float PPP-RTK, a modified version of the RTKLib 2.4.3 (beta) package is used to apply network corrections. First-order ionospheric effects were eliminated by the iono-free combination and zenith tropospheric delay estimated. The corrections were applied by introducing a priori constrained tropospheric parameters. Adaptive modeling based on OFCs (Optimal Fitting Coefficients) has been developed to describe the behavior of the troposphere, using estimates of tropospheric delays for Orpheon stations. This solution allows one-way communication between the server and the user. The gains achieved in convergence time to 10 centimeters accuracy were statistically quantified. Network topology was assessed by reducing the number of reference stations (up to 75%) using a sparse network configuration. In the second step, PPP-RTK is realized using the PPP-Wizard 1.3 software and CNES (Centre National d'Etudes Spatiales) real-time products for orbits, clocks and phase biases of satellites. The RT-IPPP (RT-Integer PPP) is performed with estimation of tropospheric and ionospheric delays. Ionospheric and tropospheric corrections are introduced as a priori parameters constrained in PPP-RTK. To generate ionospheric corrections, it was implemented an Inverse Distance Weighting (IDW) algorithm. Improvements achieved in horizontal positioning due to external SSR corrections from a (dense or sparse) network are promising and may be useful for applications that depend primarily on horizontal positioning. / O PPP (Precise Point Positioning) é um método GNSS (Global Navigation Satellite Systems) baseado no conceito SSR (State Space Representation). Graças às melhorias recentes nos modelos atmosféricos, o PPP em tempo real (RT-PPP) também pode ser aprimorado. O objetivo principal deste trabalho é estudar o RT-PPP e a infraestrutura otimizada em termos de custos e benefícios para realizar o método usando correções atmosféricas. Portanto, são utilizadas diferentes configurações de uma rede GNSS densa e regular existente na França, a rede Orphéon. Esta rede tem cerca de 160 estações, sendo propriedade da Geodata-Diffusion (Hexagon Geosystems). Inicialmente, foi avaliado o "float PPP-RTK", que corresponde ao RT-PPP com melhorias resultantes de correções de rede, embora mantendo as ambiguidades como float. Em um segundo momento, as correções de rede são aplicadas para aprimorar o "PPP-RTK", onde ambiguidades são fixadas para seus valores inteiros. Para o float PPP-RTK, uma versão modificada do software RTKLib 2.4.3 (beta) é empregada de modo a levar em consideração as correções de rede. Os efeitos ionosféricos de primeira ordem foram eliminados pela combinação iono-free e o atraso troposférico é estimado. As correções são aplicadas introduzindo parâmetros troposféricos a priori injuncionados. Uma modelagem adaptativa baseada em OFCs (Optimal Fitting Coefficients) foi implementada para descrever o comportamento da troposfera, utilizando estimativas de atraso troposférico para estações da rede Orpheon. Tal solução permite a comunicação unidirecional entre o servidor e o usuário. Os ganhos alcançados no tempo de convergência para acurácia de 10 centímetros foram quantificados estatisticamente. A topologia de rede foi avaliada reduzindo o número de estações de referência (até 75%) usando uma configuração de rede esparsa. Na segunda etapa, o PPP-RTK é realizado usando o software PPP-Wizard 1.3, bem como os produtos para tempo real do CNES (Centre National d’Etudes Spatiales) de órbitas, relógios e biases de fase de satélites. O RT-IPPP (RT-Integer PPP) é realizado com estimativa de atrasos troposféricos e ionosféricos. As correções ionosféricas e troposféricas são introduzidas como parâmetros a priori injuncionados no PPP-RTK. Para gerar correções ionosféricas, foi implementado um algoritmo baseado na ponderação pelo inverso da distância (IDW–Inverse Distance Weighting). As melhorias alcançadas no posicionamento horizontal com o uso das correções SSR externas de uma rede (densa ou esparsa) são promissoras e podem ser úteis para aplicações que dependem principalmente do posicionamento horizontal.

Made available in DSpace on 2016-03-15T19:38:53Z (GMT). No. of bitstreams: 1 ALBERTO FARES AKEL JUNIOR.pdf: 5112998 bytes, checksum: f18fc33d2f9508c3ec265c0efa016b43 (MD5) Previous issue date: 2015-08-21 / Fundação de Amparo a Pesquisa do Estado de São Paulo / We study the behavior of the Earth-ionosphere waveguide through the modeling of the propagation of very low frequency radio waves (VLF). We use the computational model LWPC (Long Wave Propagation Capability) to estimate changes in amplitude and phase of the VLF signals detected by the SAVNET network (South America VLF NETwork), and thus try to understand the behavior of the lower ionosphere under different ionization conditions. The research was divided into two parts. The first part investigates the behavior of the VLF signals in quiescent regimes of ionization. Amplitude and phase simulations for the were carried out, modifying adapting polynomials for the β and h parameters (or Wait s parameters) as a function of the zenithal angle. The second part of this research, uses these polynomials in the study of the lower ionosphere under transient ionization regimes in two distinct conditions: first during of solar flares and second during solar eclipse. For the simulations under solar flare conditions, we calculate the changes in β and ℎ′ parameters during the 25/03/2008 solar explosion. With these values, we calculate the electronic density profile through an exponential model and we find that the electronic density at 75 km is ∼ 104 cm−3, that is twenty times higher than during quiescent conditions. To evaluate our parameter estimates, we calculate the variation of the Wait s parameters for the case of twelve solar events of different classes. We note that the variations Δℎ′ found in this work are larger than that in Muraoka, Murata e Sato (1977) because they consider the variations in the conductivity gradient. For the solar eclipse simulations on 11/07/2011, we investigate its effect on the VLF phase. For this, we use the obscuration coefficient to estimate the guide height variation along the whole path during the eclipse. The simulations reproduce the phase behavior during the eclipse. However, a delay of about twenty four minutes was observed between the simulated and observed measurements. The observed delay is a direct consequence of own estimates of the perturbed ionospheric height and it causal relation with the obscuration during the eclipse. lower ionosphere, VLF, modeling, ionospheric disturbances, solar flares, solar eclipse. / Neste trabalho realizamos o estudo do comportamento do guia de ondas terra-ionosfera através da modelagem da propagação ondas de rádio de frequência muito baixa (VLF). Para isto, utilizamos o modelo computacional LWPC (Long Wave Propagation Capability) para estimar as variações de amplitude e fase de sinais de VLF detectados nos trajetos da rede SAVNET (South America VLF NETwork) e assim compreender o comportamento da baixa ionosfera em diferentes regimes de ionização. A pesquisa foi dividida em duas partes. A primeira parte, investigou o comportamento do sinal VLF em regimes quiescente de ionização, assim realizou-se simulações de amplitude e fase adaptando polinômios que definem os parâmetros β e ℎ′ (ou parâmetros de Wait) em função do ângulo zenital solar. Na segunda parte desta pesquisa, aplicou-se os polinômios no estudo da baixa ionosfera sob regimes transientes de ionização em duas condições distintas. A primeira para o caso de explosões solares e a segunda um para eclipse solar. Nas simulações relativas a explosões solares, calculamos as variações dos parâmetros β e ℎ′ durante o evento do dia 25/03/2008. Com esses valores, calculamos o perfil de densidade eletrônica, através de um modelo exponencial e observamos que a densidade eletrônica em 75 km é ∼ 104 cm−3, ou seja, vinte vezes maiores que antes da explosão. Para avaliar nossas estimativas, calculamos a variação dos parâmetros de Wait para doze eventos de diferentes classes. Observamos que as variações Δℎ′ neste trabalho são sempre maiores do que as descritas em Muraoka, Murata e Sato (1977), devido elas considerarem as variações no gradiente de condutividade. Nas simulações relativa ao eclipse solar do dia 11/07/2011, investigamos seu efeito na fase observada. Para esse estudo, utilizou-se o coeficiente de obscurecimento para realizar as simulações, desta forma foi possível estimar a variação da altura do guia ao longo de todo o trajeto durante o eclipse. As simulações reproduziram o comportamento da fase durante o eclipse. Entretanto, foi observado um atraso entre as medidas calculadas e observadas de aproximadamente ∼ vinte e quatro minutos. O atraso observado é diretamente decorrente da estimativa da altura de referência da ionosfera pertubada e de sua relação causal com o obscurecimento durante o eclipse.

The genetic algorithm (GA) is a popular random search and optimization method inspired by the concepts of crossover, random mutation, and natural selection from evolutionary biology. The real-valued genetic algorithm (RGA) is an improved version of the genetic algorithm designed for direct operation on real-valued variables. In this work, a modified version of a genetic algorithm is introduced, which is called a modified genetic algorithm with micro-movement (MGAM). It implements a particle swarm optimization(PSO)-inspired micro-movement phase that helps to improve the convergence rate, while employing the e'cient GA mechanism for maintaining population diversity. In order to test the capability of the MGAM, we firrst implement it on five generally used test functions. Then we test the MGAM on two typical nonlinear dynamical systems. The performance of the MGAM is compared to a basic RGA on all these applications. Finally, we implement the MGAM on the most important application, which is the plasma physics-based model of the solar wind-driven magnetosphere-ionosphere system (WINDMI). In order to use this model for real-time prediction of geomagnetic activity, the model parameters require up-dating every 6-8 hours. We use the MGAM to train the parameters of the model in order to achieve the lowest mean square error (MSE) against the measured auroral electrojet (AL) and Dst indices. The performance of the MGAM is compared to the RGA on historical geomagnetic storm datasets. While the MGAM performs substantially better than the RGA when evaluating standard test functions, the improvement is about 6-12 percent when used on the 20D nonlinear dynamical WINDMI model.

This thesis addresses solitary waves and their significance for auroral particle acceleration, coronal heating and incoherent scatter radar spectra. Solitary waves are formed due to a balance of nonlinear and dispersive effects. There are several nonlinearities present in ideal magnetohydrodynamics (MHD) and dispersion can be introduced by including theHall termin the generalised Ohm’s law. The resulting system of equations comprise the classical ideal MHD waves, whistlers, drift waves and solitarywave solutions. The latter reside in distinct regions of the phase space spanned by the speed and the angle (to the magnetic field) of the propagating wave. Within each region, qualitatively similar solitary structures are found. In the limit of neglected electron intertia, the solitary wave solutions are confined to two regions of slow and fast waves, respectively. The slow (fast) structures are associated with density compressions (rarefactions) and positive (negative) electric potentials. Such negative potentials are shown to accelerate electrons in the auroral region (solar corona) to tens (hundreds) of keV. The positive electric potentials could accelerate solar wind ions to velocities of 300–800 km/s. The structure widths perpendicular to themagnetic field are in the Earth’s magnetosphere (solar corona) of the order of 1–100 km (m). This thesis also addresses a type of incoherent scatter radar spectra, where the ion line exhibits a spectrally uniform power enhancement with the up- and downshifted shoulder and the spectral region in between enhanced simultaneously and equally. The power enhancements are one order of magnitude above the thermal level and are often localised to an altitude range of less than 20 km at or close to the ionospheric F region peak. The observations are well-described by a model of ion-acoustic solitary waves propagating transversely across the radar beam. Two cases of localised ion line enhancements are shown to occur in conjunction with auroral arcs drifting through the radar beam. The arc passages are associated with large gradients in ion temperature, which are shown to generate sufficiently high velocity shears to give rise to growing Kelvin-Helmholtz (K-H) instabilities. The observed ion line enhancements are interpreted in the light of the low-frequency turbulence associated with these instabilities. / Denna avhandling handlar om solitära vågor och deras roll i norrskensacceleration och koronaupphettning, samt deras signatur i spektra uppmätta med inkoherent spridningsradar. Solitära vågor bildas genom en balans mellan ickelinjära och dispersiva effekter. Ickelinjäriteter finns det gott om i ideal magnetohydrodynamik (MHD) och dispersion kan införas genom att inkludera Halltermen i den generaliserade Ohms lag. Det resulterande ekvationssystemet omfattar de klassiska vågorna inom ideal MHD, visslare, driftvågor och solitära vågor. De sistnämnda återfinns i väldefinierade områden i fasrummet som spänns upp av farten och vinkeln (mot magnetfältet) för den propagerande vågen. Inom varje sådant område återfinns kvalitativt lika solitära våglösningar. Om man försummar elektronernas tröghet begränsas de solitära våglösningarna till två områden med långsamma respektive snabba vågor. De långsamma (snabba) strukturerna är associerade med täthets-kompressioner (förtunningar) och positiva (negativa) elektriska potentialer. De negativa potentialerna visas kunna accelerera elektroner i norrskensområdet (solens korona) till tiotals (hundratals) keV medan de positiva potentialerna accelererar solvindsjoner till hastigheter på 300–800 km/s. Strukturbredderna vinkelrät mot magnetfältet är i jordens magnetosfär (solens korona) av storleksordningen 1–100 km (m). Denna avhandling tar även upp en typ av inkoherent spridningsradarspektra, där jonlinjen uppvisar en spektralt uniform förstärkning. Detta innebär att den upp- och nedskiftade skuldran och spektralbandet däremellan förstärks simultant och i lika hög grad. Effektförstärkningen är en storleksordning över den termiska nivån och är ofta lokaliserad till ett höjd-intervall av mindre än 20 km nära jonosfärens F-skiktstopp. Observationerna beskrivs väl av en modell med solitära vågor som propagerar transversellt genom radarstrålen. Två fall av lokaliserade jonlinjeförstärkningar visas sammanfalla med att norrskensbågar driver genom radarstrålen. I samband med bågarnas passage uppmäts stora gradienter i jontemperatur, vilket visas skapa tillräckligt kraftiga hastighetsskjuvningar för att Kelvin-Helmholtz-instabiliteter ska tillåtas växa. De observerade jonlinjeförstärkningarna tolkas i skenet av den lågfrekventa turbulensen som är kopplad till dessa instabiliteter.