A multi-instrumental approach to the study of equatorial ionosphere over South America

An extensive study of the morphology and the dynamics of the equatorial ionosphere over South America is presented here. A multi parametric approach is used to describe the physical characteristics of the ionosphere in the regions where the combination of the thermospheric electric field and the horizontal geomagnetic field creates the so-called Equatorial Ionization Anomalies. Ground based measurements from GNSS receivers are used to link the Total Electron Content (TEC), its spatial gradients and the phenomenon known as scintillation that can lead to a GNSS signal degradation or even to a GNSS signal ‘loss of lock’. A new algorithm to highlight the features characterizing the TEC distribution is developed in the framework of this thesis and the results obtained are validated and used to improve the performance of a GNSS positioning technique (long baseline RTK). In addition, the correlation between scintillation and dynamics of the ionospheric irregularities is investigated. By means of a software, here implemented, the velocity of the ionospheric irregularities is evaluated using high sampling rate GNSS measurements. The results highlight the parallel behaviour of the amplitude scintillation index (S4) occurrence and the zonal velocity of the ionospheric irregularities at least during severe scintillations conditions (post-sunset hours). This suggests that scintillations are driven by TEC gradients as well as by the dynamics of the ionospheric plasma. Finally, given the importance of such studies for technological applications (e.g. GNSS high-precision applications), a validation of the NeQuick model (i.e. the model used in the new GALILEO satellites for TEC modelling) is performed. The NeQuick performance dramatically improves when data from HF radar sounding (ionograms) are ingested. A custom designed algorithm, based on the image recognition technique, is developed to properly select the ingested data, leading to further improvement of the NeQuick performance.

Tackling ionospheric scintillation threat to GNSS in Latin America

Scintillations are rapid fluctuations in the phase and amplitude of transionospheric radio signals which are caused by small-scale plasma density irregularities in the ionosphere. In the case of the Global Navigation Satellite System (GNSS) receivers, scintillation can cause cycle slips, degrade the positioning accuracy and, when severe enough, can even lead to a complete loss of signal lock. Thus, the required levels of availability, accuracy, integrity and reliability for the GNSS applications may not be met during scintillation occurrence; this poses a major threat to a large number of modern-day GNSS-based applications. The whole of Latin America, Brazil in particular, is located in one of the regions most affected by scintillations. These effects will be exacerbated during solar maxima, the next predicted for 2013. This paper presents initial results from a research work aimed to tackle ionospheric scintillation effects for GNSS users in Latin America. This research is a part of the CIGALA (Concept for Ionospheric Scintillation Mitigation for Professional GNSS in Latin America) project, co-funded by the EC Seventh Framework Program and supervised by the GNSS Supervisory Authority (GSA), which aims to develop and test ionospheric scintillation countermeasures to be implemented in multi-frequency, multi-constellation GNSS receivers.

Diagnostics of equatorial and low latitude ionosphere by TEC mapping over Brazil

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Dependence of convective flows and particle precipitation in the high-latitude dayside ionosphere on theXandYcomponents of the interplanetary magnetic field

The asymmetries in the convective flows, current systems, and particle precipitation in the high-latitude dayside ionosphere which are related to the equatorial plane components of the interplanetary magnetic field (IMF) are discussed in relation to the results of several recent observational studies. It is argued that all of the effects reported to date which are ascribed to the y component of the IMF can be understood, at least qualitatively, in terms of a simple theoretical picture in which the effects result from the stresses exerted on the magnetosphere consequent on the interconnection of terrestrial and interplanetary fields. In particular, relaxation under the action of these stresses allows, in effect, a partial penetration of the IMF into the magnetospheric cavity, such that the sense of the expected asymmetry effects on closed field lines can be understood, to zeroth order, in terms of the “dipole plus uniform field” model. In particular, in response to IMF By, the dayside cusp should be displaced in longitude about noon in the same sense as By in the northern hemisphere, and in the opposite sense to By in the southern hemisphere, while simultaneously the auroral oval as a whole should be shifted in the dawn-dusk direction in the opposite sense with respect to By. These expected displacements are found to be consistent with recently published observations. Similar considerations lead to the suggestion that the auroral oval may also undergo displacements in the noon-midnight direction which are associated with the x component of the IMF. We show that a previously published study of the position of the auroral oval contains strong initial evidence for the existence of this effect. However, recent results on variations in the latitude of the cusp are more ambiguous. This topic therefore requires further study before definitive conclusions can be drawn.

Application of ionospheric corrections in the equatorial region for L1 GPS users

In the absence of the selective availability, which was turned off on May 1, 2000, the ionosphere can be the largest source of error in GPS positioning and navigation. Its effects on GPS observable cause a code delays and phase advances. The magnitude of this error is affected by the local time of the day, season, solar cycle, geographical location of the receiver and Earth's magnetic field. As it is well known, the ionosphere is the main drawback for high accuracy positioning, when using single frequency receivers, either for point positioning or relative positioning of medium and long baselines. The ionosphere effects were investigated in the determination of point positioning and relative positioning using single frequency data. A model represented by a Fourier series type was implemented and the parameters were estimated from data collected at the active stations of RBMC (Brazilian Network for Continuous Monitoring of GPS satellites). The data input were the pseudorange observables filtered by the carrier phase. Quality control was implemented in order to analyse the adjustment and to validate the significance of the estimated parameters. Experiments were carried out in the equatorial region, using data collected from dual frequency receivers. In order to validate the model, the estimated values were compared with ground truth. For point and relative positioning of baselines of approximately 100 km, the values of the discrepancies indicated an error reduction better than 80% and 50% respectively, compared to the processing without the ionospheric model. These results give an indication that more research has to be done in order to provide support to the L1 GPS users in the Equatorial region.

Unusual early morning development of the equatorial anomaly in the Brazilian sector during the Halloween magnetic storm

The solar events that occurred at the end of October 2003 gave rise to very strong geomagnetic disturbances that peaked twice with Dst values reaching -345 nT around 0000 UT on 30 October and -400 nT around 2300 UT, on the same day. Disturbances in several ionospheric parameters were observed over Brazil. This work will focus on the ionospheric response to the initial westward prompt penetration electric field and on the strong intensification of the equatorial ionization anomaly that occurred because of the electric field polarity reversal that followed in the early morning hours of 29 October. The F layer peak height over the equator first decreased under the strong prompt penetration westward electric field, which was followed by significant height increase under eastward electric field. We have used Sheffield University Plasmasphere Ionosphere Model (SUPIM) with an intensified westward disturbed electric field in the presunrise hours, presumably due to prompt penetration from the magnetosphere, in order to study the effect of such a field in the ionosphere. The simulation results showed that prompt penetration of magnetospheric electric fields of westward polarity to the nightside equatorial region seems to be the most probable cause of the initial F layer height decreases. The intensification of the equatorial ionization anomaly and the unusual enhancement on F layer peak density, which was not modeled by the SUPIM, are explained as caused by the strong eastward electric field that followed the initial phase in combination with a highly variable disturbed meridional/transequatorial wind system as inferred from the F2 layer peak height variations. The highly dynamic wind pattern, with a short-term response (2-4 hours), is compatible with the predictions of some previous theoretical model calculations reported in the literature.

Posicionamento relativo na região equatorial em diversas condições ionosféricas

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Influência de estruturas geológicas bidimensionais no campo geoeletromagnético na presença do eletrojato equatorial

A Terra atua como um grande magneto esférico, cujo campo assemelha-se àquele gerado por um dipolo magnético. Este campo apresenta mudanças de intensidade que variam com a localização e a hora local. A parte principal do campo geomagnético se origina no interior da Terra através de processos eletromagnéticos. Extensivos estudos mostraram ainda que existem contribuições de origem externa ao planeta, principalmente de origem solar. Dentre estas fontes há anomalias do campo magnético que surgem a partir de um aumento diurno da corrente elétrica em uma estreita faixa da ionosfera, de direção leste-oeste, centrada no equador magnético e denominada Eletrojato Equatorial (EEJ). Ocasionalmente estas correntes podem apresentar reversões de fluxo, sendo denominadas Contra-Eletrojato (CEJ). Vários autores têm estudado os efeitos do EEJ e CEJ sobre as observações geoeletromagnéticas. Eles estão interessados no efeito combinado do EEJ e estruturas geológicas condutivas 1-D e 2-D. Nestes trabalhos a estrutura 2-D sempre se apresentava paralela ao eletrojato, o que é uma hipótese bastante restritiva ao se modelar ambientes geológicos mais realistas, em que corpos bidimensionais podem ter qualquer strike em relação ao EEJ. Neste trabalho apresentamos a solução deste problema sem esta restrição. Assim, mostramos os campos geoeletromagnéticos devidos a estruturas bidimensionais que possuam strike oblíquo em relação ao EEJ, através de perfis dos campos elétrico e magnéticos calculados na superfície e formando direção arbitrária à heterogeneidade condutiva 2-D. Com esta resposta avaliamos ainda qual a influência que estruturas bidimensionais exercem sobre a resposta magnetotelúrica, sob influência do Eletrojato Equatorial. Durante o desenvolvimento deste trabalho, utilizamos o método de elementos finitos, tendo por fonte eletromagnética o EEJ e o CEJ, que por sua vez foram representados por uma combinação de distribuições gaussianas de densidade de corrente. Estas fontes foram decompostas nas direções paralela e perpendicular à estrutura 2-D, resultando nos modos de propagação TE₁ e TE₂ e TM acoplados, respectivamente. Resolvemos o modo acoplado aplicando uma Transformada de Fourier nas equações de Maxwell e uma Transformada Inversa de Fourier na solução encontrada. De acordo com os experimentos numéricos realizados em um modelo interpretativo da Anomalia Condutiva da Bacia do Parnaíba, formado por uma enorme estrutura de 3000 ohm-m dentro de um corpo externo condutivo (1 ohm-m), concluímos que a presença do CEJ causa uma inversão na anomalia, se compararmos com o resultado do EEJ. Concluímos também que para as frequências mais altas as componentes do campo elétrico apresentam menor influência da parte interna do corpo 2-D do que da parte externa. Já para frequências mais baixas este comportamento se observa com as componentes do campo magnético. Com relação à frequência, vimos os efeitos do “skin-depth”, principalmente nas respostas magnéticas. Além disso, quando a estrutura 2-D está paralela ao eletrojato, o campo elétrico é insensível à estrutura interna do modelo para todos os valores de frequência utilizados. Com respeito ao ângulo θ_h entre a heterogeneidade e a fonte, vimos que o modo TM se manifesta naturalmente quando θ_h é diferente de 0°. Neste caso, o modo TE é composto por uma parte devido à componente da fonte paralela à heterogeneidade e a outra devido à componente da fonte perpendicular, que é acoplada ao modo TM. Assim, os campos calculados têm relação direta com o valor de θ_h. Analisando a influência do ângulo entre a direção do perfil dos campos e o strike da heterogeneidade verificamos que, à medida que θ_h se aproxima de 90°, os campos primários tornam-se variáveis para valores de θ_p diferentes de 90°. Estas variações causam uma assimetria na anomalia e dão uma idéia da inclinação da direção do perfil em relação aos corpos. Finalmente, concluímos que uma das influências que a distância entre o centro do EEJ e o centro da estrutura 2-D, causa sobre as componentes dos campos está relacionado às correntes reversas do EEJ e CEJ, pois a 500 km do centro da fonte estas correntes têm máxima intensidade. No entanto, com o aumento da distância, as anomalias diminuem de intensidade. Nas sondagens MT, nós também usamos o EEJ e o CEJ como fonte primária e comparamos nossos resultados com a resposta da onda plana. Deste modo observamos que as componentes do campo geoeletromagnético, usadas para calcular a impedância, têm influência do fator de acoplamento entre os modos TE₂ e TM. Além disso, esta influência se torna maior em meios resistivos e nas frequências mais baixas. No entanto, o fator de acoplamento não afeta os dados magnetotelúricos em frequências maiores de 10^-2 Hz. Para frequências da ordem de 10^-4 Hz os dados MT apresentam duas fontes de perturbação: a primeira e mais evidente é devido à presença fonte 2-D (EEJ e CEJ), que viola a hipótese da onda plana no método MT; e a segunda é causada pelo acoplamento entre os modos TE₂ e TM, pois quando a estrutura bidimensional está obliqua à fonte 2-D temos correntes elétricas adicionais ao longo da heterogeneidade. Concluimos assim, que o strike de uma grande estrutura condutiva bidimensional relativamente à direção do EEJ ou CEJ tem de fato influência sobre o campo geomagnético. Por outro lado, para estudos magnetotelúricos rasos (frequências maiores de 10^-3 Hz) o efeito do ângulo entre a estrutura geológica 2-D e a direção do EEJ não é tão importante. Contudo, em estudos de litosfera frequências menores de 10^-3 Hz) o acoplamento entre os modos TE₂ e TM não pode ser ignorado.

Avaliação da influência do eletrojato equatorial na interpretação de dados magnetotelúrico bidimensionais

O campo magnetotelúrico em regiões equatoriais viola a condição de ondas planas por causa de uma fonte fortemente concentrada na direção E-W na ionosfera, denominada eletrojato equatorial. No presente trabalho, procurou-se analisar a resposta magnetotelúrica de fontes que simulam o efeito do eletrojato equatorial. Foram considerados dois tipos de fontes para simular o eletrojato: uma linha infinita de corrente e uma distribuição gaussiana de densidade de corrente em relação a uma das coordenadas horizontais. A resistividade aparente foi obtida através da relação de Cagniard e comparada com os resultados de ondas planas. É mostrada também a comparação entre a fase da impedância na superfície, para os três tipos de fontes (ondas planas, eletrojato gaussiano e linha de corrente). O problema de meios com heterogeneidades laterais foi resolvido em termos de campos secundários, sendo as equações diferenciais solucionadas através da técnica de elementos finitos bidimensionais. Os resultados mostram que o eletrojato tem pouca influência nas respostas (resistividade aparente e fase) de estruturas geológicas rasas. Entretanto, a influência pode ser considerável nas estruturas profundas (maior que 5000 m), principalmente se suas resistividades são altas (maior que 100 Ω.m). Portanto, a influência do eletrojato equatorial deve ser considerada na interpretação de dados magnetotelúricos de bacias sedimentares profundas ou no estudo da crosta terrestre.

Estudos de distúrbios ionosféricos propagantes de larga escala (LSTID’S) na região equatorial durante grandes tempestades magnéticas: estudo de casos

The Ionospheric Disturbances – TIDs – are irregularities on the ionospheric plasma propagating in speeds in the order of tens to a few hundreds of meters per second. This present study detected and characterized the TIDs of LSTIDs (Large Scale Travelling Ionospheric Disturbance) type at low latitudes during intense geomagnetic storms and its propagation over the Brazilian sector. This work also shows as being the first to report systematically propagation of gravity waves over Natal. For this purpose, we used ionospheric records obtained from type of digisonde CADI (Canadiam Advanced Digital Ionosonde) located in Natal and the type DSP (Digisonde Portable Souder) located in Cachoeira Paulista, Fortaleza and São Luis, whereupon we used a dataset of 12 years collected by INPE (National Institute of Space Research). In this study, both calm days, that preceded the storms, and the geomagnetically disturbed days were related during the years 2000 and 1012, which cover a period of maximum and minimum solar activity. And it is presented the variations that happened in the electron density from region F of the ionosphere over the Brazilian sector, especially near the Equator (Natal, Fortaleza and São Luis), caused by ionospheric disturbances in the equatorial region during intense geomagnetic storms, because, as we know of the literature in this area, this phenomenon contributes positively to the emergence of LSTIDs in the auroral region, which may move to the equatorial region where a few cases have been documented and studied systematically. From the observation of signatures if TIDs in ionogram records, a study of the morphology of these events was performed and compared with the main characteristics of the wave of this phenomenon during great magnetic storms, i.e., DST <(-200 nT) and KP > 6. Thus, we obtained the main characteristics of TIDs over our region, i.e., period, vertical wavelength, phase and propagation speed, as well as the delay of these disturbances compared to the beginning of the magnetic storms to the Brazilian Sector.

Equatorial modes observed in atmospheric variables

Wavenumber-frequency spectral analysis of different atmospheric variables has been carried Out using 25 years of data. The area considered is the tropical belt 25 degrees S-25 degrees N. A combined FFT wavelet analysis method has been used for this purpose. Variables considered are outgoing long wave radiation (OLR), 850 hPa divergence, zonal and meridional winds at 850, 500 and 200 hPa levels, sea level pressure and 850 hPa geopotential height. It is shown that the spectra of different variables have some common properties, but each variable also has few features diffe:rent from the rest. While Kelvin mode is prominent in OLR, and zonal winds, it is not clearly observed in pressure and geopotential height fields; the latter two have a dominant wavenumber zero mode not seen in other variables except in meridional wind at 200 hPa and 850 hPa divergences. Different dominant modes in the tropics show significant variations on sub-seasonal time scales.

The effect of solar activity on temperatures in the equatorial mesosphere

A study of the response of neutral temperatures in the equatorial mesosphere to variations in solar activity has been carried out by investigating the correlation between the 10.7 cm solar radio flux and temperatures obtained from a series of 51 rocket soundings conducted over Thumba, India (8°N, 77°E) during the period December 1970–December 1971. A strong positive correlation between these two parameters has been obtained, indicating mesospheric heating effects caused by day-to-day variations in solar EUV emission. The correlation analysis indicates that this response persists over several days and that the peak correlation between the temperatures and the F10.7 index occurs with a time lag of less than 24 hr.