Simultaneous observations of some unusual whistlers and VLF hiss emission at a low L-value (L = 1.17)

The purpose of this paper is to describe some unusual whistlers (doublets and triplets) and VLF hiss emission recorded simultaneously on February 18, 1998 during nighttime at a low latitude Indian ground station Jammu (geomag. lat., 22° 26′ N; L = 1.17), and to make some discussions about their origin. The detailed structures of the observed VLF hiss emission clearly show that these emissions confined to a narrow continuous frequency band. Some times the frequency of hiss band oscillates and subsequently touches the upper edge of the first whistler component of the doublet. Detailed structures of the dynamic spectra of whistler and VLF hiss emission are briefly presented. From the dispersion analysis of the whistler doublets and triplets, it is found that the individual whistlers of the doublets and triplets simultaneously observed on the same day are one-hop whistlers having propagation path along higher and closely spaced L-values. Our result also shows that VLF hiss has been generated in the equatorial region of higher L-values. Generation and propagation mechanisms are briefly discussed.     

An explanation of the observation of pulsing hiss at low latitude

In this paper we report pulsing hiss emissions observed at the low latitude station, Jammu (geomag. lat. 22°26′N, L = 1.17) in which intensity decreases with the increase in frequency. The entire dynamic spectra contain somewhat irregular structure. To explain these we propose that the hiss emissions are generated through Doppler-shifted cyclotron interactions near the equator and propagate to the earth in the whistler-mode. Further, ULF waves present in the generation region modulate the intensity of the emission resulting in the pulsing nature. The growth rates are computed and discussed in the light of recent works.     

Response of low latitude D-region ionosphere to the total solar eclipse of 22 July 2009 deduced from ELF/VLF analysis

Response of the D-region of the ionosphere to the total solar eclipse of 22 July 2009 at low latitude, Varanasi (Geog. lat., 25.27° N; Geog. long., 82.98° E; Geomag. lat. = 14° 55’ N) was investigated using ELF/VLF radio signal. Tweeks, a naturally occurring VLF signal and radio signals from various VLF navigational transmitters are first time used simultaneously to study the effect of total solar eclipse (TSE). Tweeks occurrence is a nighttime phenomena but the obscuration of solar disc during TSE in early morning leads to tweek occurrence. The changes in D-region ionospheric VLF reflection heights (h) and electron density (n_e: 22.6–24.6 cm⁻³) during eclipse have been estimated from tweek analysis. The reflection height increased from ∼89 km from the first occurrence of tweek to about ∼93 km at the totality and then decreased to ∼88 km at the end of the eclipse, suggesting significant increase in tweek reflection height of about 5.5 km during the eclipse. The reflection heights at the time of totality during TSE are found to be less by 2–3 km as compared to the usual nighttime tweek reflection heights. This is due to partial nighttime condition created by TSE. A significant increase of 3 dB in the strength of the amplitude of VLF signal of 22.2 kHz transmitted from JJI-Japan is observed around the time of the total solar eclipse (TSE) as compared to a normal day. The modeled electron density height profile of the lower ionosphere depicts linear variation in the electron density with respect to solar radiation as observed by tweek analysis also. These low latitude ionospheric perturbations on the eclipse day are discussed and compared with other normal days.     

Height and time variability of planetary wave activity

The height–season and year-to-year regularities of parameters of first and second spatial harmonics determine the structure of the stratosphere and mesosphere circulation and its variability. In the period 1992–2002 at heights 0–55 km, the amplitudes and phases of the first and second spatial harmonics in the field of temperature, geopotential height, zonal and meridional wind were calculated by the method of harmonic decomposition. Dispersion (standard or mean square deviation) of their day-to-day and year-to-year variations was calculated by their wavelength constants. Height and season patterns of variability have been estimated. The difference in height–longitude variability for wave numbers m = 1 and 2 has been discovered. At the same time, the intensity of wave disturbances for m = 1 is less than for m = 2 excluding the polar areas, where a significant variability appears at the heights 0–55 km. There is also a tendency for the intensity of year-to-year variations to decrease in comparison with day-to-day variations. In cold and warm periods the amplitude of perturbation waves with m = 2 both for day-to-day and year-to-year variations is greater than for waves with m = 1. Transient height areas in the interval of 20–30 km are more distinct for day-to-day variations of polar area.     

One-to-one relationship between low latitude whistlers and conjugate source lightning discharges and their propagation characteristics

One-to-one relation with its causative lightning discharges and propagation features of night-time whistlers recorded at low-latitude station, Allahabad (geomag. lat. 16.05°N, L = 1.08), India, from continuous observations made during 1–7 April, 2011 have been studied. The whistler observations were made using the Automatic Whistler Detector (AWD) system and AWESOME VLF receiver. The causative lightning strikes of whistlers were checked in data provided by World-Wide Lightning Location Network (WWLLN). A total of 32 whistlers were observed out of which 23 were correlated with their causative lightnings in and around the conjugate location (geom. lat. 9.87°S) of Allahabad. A multi-flash whistler is also observed on 1 April with dispersions 15.3, 17.5 and 13.6 s^1/2. About 70% (23 out of 32) whistlers were correlated with the WWLLN detected causative lightnings in the conjugate region which supports the ducted mode of propagation at low latitude. The multi-flash and short whistlers also propagated most likely in the ducted mode to this station.     

3D Imaging of the OH mesospheric emissive layer

A new and original stereo imaging method is introduced to measure the altitude of the OH nightglow layer and provide a 3D perspective map of the altitude of the layer centroid. Near-IR photographs of the OH layer are taken at two sites separated by a 645 km distance. Each photograph is processed in order to provide a satellite view of the layer. When superposed, the two views present a common diamond-shaped area. Pairs of matched points that correspond to a physical emissive point in the common area are identified in calculating a normalized cross-correlation coefficient (NCC). This method is suitable for obtaining 3D representations in the case of low-contrast objects. An observational campaign was conducted in July 2006 in Peru. The images were taken simultaneously at Cerro Cosmos (12°09′08.2″ S, 75°33′49.3″ W, altitude 4630 m) close to Huancayo and Cerro Verde Tellolo (16°33′17.6″ S, 71°39′59.4″ W, altitude 2272 m) close to Arequipa. 3D maps of the layer surface were retrieved and compared with pseudo-relief intensity maps of the same region. The mean altitude of the emission barycenter is located at 86.3 km on July 26. Comparable relief wavy features appear in the 3D and intensity maps. It is shown that the vertical amplitude of the wave system varies as exp (Δz/2H) within the altitude range Δz = 83.5–88.0 km, H being the scale height. The oscillatory kinetic energy at the altitude of the OH layer is comprised between 3 × 10⁻⁴ and 5.4 × 10⁻⁴ J/m³, which is 2–3 times smaller than the values derived from partial radio wave at 52°N latitude.     

Variation of ionospheric total electron content in Indian low latitude region of the equatorial anomaly during May 2007–April 2008

The ionospheric total electron content (TEC), derived by analyzing dual frequency signals from the Global Positioning System (GPS) recorded near the Indian equatorial anomaly region, Varanasi (geomagnetic latitude 14°, 55′N, geomagnetic longitude 154°E) is studied. Specifically, we studied monthly, seasonal and annual variations as well as solar and geomagnetic effects on the equatorial ionospheric anomaly (EIA) during the solar minimum period from May 2007 to April 2008. It is found that the daily maximum TEC near equatorial anomaly crest yield their maximum values during the equinox months and their minimum values during the summer. Using monthly averaged peak magnitude of TEC, a clear semiannual variation is seen with two maxima occurring in both spring and autumn. Statistical studies indicate that the variation of EIA crest in TEC is poorly correlated with Dst-index (r = −0.03) but correlated well with Kp-index (r = 0.82). The EIA crest in TEC is found to be more developed around 12:30 LT.     

The heliospheric modulation of cosmic rays: Effects of a latitude dependent solar wind termination shock

Observations made with the two Voyager spacecraft confirmed that the solar wind decelerates to form the heliospheric termination shock and that it has begun its merger with the local interstellar medium. The compression ratio of this shock affects galactic cosmic rays when they enter the heliosphere. Hydrodynamic (HD) models show that the compression ratio can have a significant latitude dependence; with the largest value in the nose direction of the heliosphere, becoming significantly less towards the polar regions. The modulation effects of such large latitude dependence are studied, using a well-established numerical drift and shock modulation model. We focus on computing the modulated spectra for galactic protons with emphasis on the radial and polar gradients in the equatorial plane and at a polar angle of θ = 55°, corresponding to the heliolatitude of Voyager 1. Two sets of solutions are computed and compared each time; with and without a latitude dependence for the compression ratio. All computations are done for the two magnetic field polarity cycles assuming solar minimum conditions. Including the termination shock in the model allows the study of the re-acceleration of galactic protons in the outer heliosphere. We find that for the A < 0 polarity cycle the intensity between ∼200 MeV and ∼1 GeV in the vicinity of the shock in the heliospheric equatorial plane may exceed the local interstellar value specified at the heliopause. Unfortunately, at θ = 55°, the effect is reduced. This seems not possible during an A > 0 cycle because significant modulation is then predicted between the heliopause and the termination shock, depending on how strong global gradient and curvature drifts are in the heliosheath. The overall effect of the shock on galactic protons in the equatorial plane is to reduce the total modulation as a function of radial distance with respect to the interstellar spectrum. Making the compression ratio latitude dependent enhances these effects at energies E < 200 MeV in the equatorial plane. At larger heliolatitudes these effects are even more significant. The differences in the modulation between the two drift cycles are compelling when the compression ratio is made latitude dependent but at Earth this effect is insignificant. A general result is that the computed radial gradient changes for galactic protons at and close to the TS and that these changes are polarity dependent. In line with previous work, large polarity dependent effects are predicted for the inner heliosphere and also close to the shock’s position in the equatorial plane. In contrast, at θ = 55°, the largest polarity effect occurs in the middle heliosphere (50 AU), enhanced by the latitude dependence of the compression ratio. At this latitude, the amount of proton modulation between the heliopause and the termination shock is much reduced. If galactic cosmic rays were to experience some diffusive shock acceleration over the 100–1000 MeV range at the shock, the radial gradient should change its sign in the vicinity of the shock, how large, depends on the compression ratio and the amount of drifts taking place in the outer heliosphere. The effective polar gradient shows a strong polarity dependence at Earth but this dissipates at θ = 55°, especially with increasing radial distance. This tendency is enhanced by making the compression ratio latitude dependent.     

Very low frequency interferometry for the Chang’E-2 project

Very low frequency interferometry among two astronomical experiments has been proposed and accepted for further study for the second phase of China’s lunar exploration programme (the Chang’E Programme), which is envisaged to operate a lander and a rover on the surface of the moon. This experiment is an interferometer experiment in the very low frequency (VLF, f < 15 MHz) regime of radio frequencies with at least degree-level angular resolution. The goals include observing solar storm activities, Coronal Mass Ejections, Auroral Kilometric Radiation, and planetary radiation in the solar system, studying the origin of Cosmic Rays, spectral properties of pulsars, surveying ionized hydrogen in the Galaxy, and exploring coherent radio emissions.     

Temperature variability in the tropical mesosphere during the northern hemisphere winter

Temperature observations at 20–90 km height and 5°N–15°N during the period of December 1992–March 1993 from the WINDII and MLS experiments on the UARS satellite are analysed together with MF radar winds and UKMO assimilated fields of temperature and zonal and meridional winds. The correlation between the different datasets at the tropics and zonal mean wind data at mid latitudes is examined for period February–March 1993, when series of stratospheric warming events were observed at middle and high latitudes. Wavelet analysis is applied to investigate coupling between stationary and travelling planetary waves in the stratosphere and the upper mesosphere. Planetary waves m = 1 with periods of 4–7 days, 8–12 days and 13–18 days are found to dominate the period. Westward 7- and 16–18 day waves at the tropics appear enhanced by stationary planetary waves during sudden stratospheric warming events.     

F3 layer feature under low and medium solar activity observed at a Chinese low latitude station Fuke

Ionogram observations from the ionosonde at Fuke (9.5°N geomagnetic latitude), a Chinese low latitude station, in 2010–2012 are analyzed to present the features of F3 layer under low and moderate solar activity conditions. Structure of the ionogram, displaying the F3 layer, was more distinct and clear during MSA than LSA periods especially during spring to summer. Start time of occurrence of the F3 layer is about at 0830–0900 LT and is approximately the same for LSA and MSA conditions. The average duration time of the F3 layer occurrence was 181 min per day under F_10.7 = 75 condition, 263 min in F_10.7 = 99 and 358 min in F_10.7 = 125, respectively. The differences of h′F2 and h′F3 exhibited obvious semiannual variation observed at Fuke from March 2010 to June 2012 and increased with increasing solar activity. The difference of foF2 and foF3 in the months February, March, September, October and November is less evident in the middle solar activity period 2011–2012 than the low solar activity 2010 and in the other period it shows a slight increase (0.5 MHz) or keeps constant. The results show that the solar activity dependence of the F3 layer occurrence at low latitude away from the magnetic equator is different from that at near the magnetic equator.     

Ionospheric precursors of earthquakes recorded by VLF receiver at Tashkent IHY station

Tashkent International Heliophysical Year (IHY) station is a member of Atmospheric Weather Electromagnetic System for Observation, Modeling and Education (AWESOME) network being operated globally to study the ionosphere and the magnetosphere with the help of electromagnetic waves in Very Low Frequency (VLF) band. Regular monitoring of the D- and F-layers of ionosphere over Central Asia territory is being performed on the permanent basis starting year 2008. We have studied VLF amplitude anomalies related to the EQs occurred in 2008–2009 years with magnitude more than 5 on the path way from the VLF transmitters to the Tashkent station assuming that propagation of VLF ground-based transmitters signals can be perturbed by EQ preparation detectable from the ground-based measurements in the VLF bands. For analyzing narrowband data we have used the nighttime fluctuation (NF) method paying attention to the data obtained during the local nighttime (20:00 LT–04:00 LT) in Tashkent where the VLF receiver is operating. The mean nighttime amplitude (or trend) and nighttime fluctuation are found to increase significantly before the EQ occurred on the path way from the transmitters to the receiver. The obtained results have revealed an agreement with VLF amplitude anomalies observed in Tashkent VLF station during the strong EQs occurred on the path way from the transmitters to the receiver. Some results are presented to show the probing potentiality of VLF waves to predict short term EQs with high magnitude.     

An intercomparison of meteor radar measurements at the South Pole using two different processing systems

A new meteor radar system was installed at the Amundsen–Scott station South Pole in 2001 to further the understanding of the dynamics of the Antarctic region. The antenna array consists of four yagis pointed along the 0°, 90°, 180°, and 270° meridians and five folded crossed dipoles arranged in a cross configuration and operating as an interferometer to provide position measurements for the detected radio meteors. The four yagis are time division multiplexed and used for both transmitting and receiving while the five folded crossed dipoles are only used for reception. The current arrangement of data acquisition (DAQ) systems at the South Pole allows the collection of meteors in a configuration similar to the previous meteor radar system that operated at the South Pole in the mid 1990s while also using an interferometer to accurately determine the meteor positions in the sky, which enables the determination of the vertical structure of the observed waves. This has been accomplished through the use of two DAQ and post-processing systems: COBRA (Colorado Obninsk radar) connected to the yagis and MEDAC (meteor echo detection and collection) connected to the folded crossed dipoles. With two separate DAQ systems operating in parallel we have the ability to directly compare the results and understand the inherent variability in the derived scientific results based on different system architectures and processing assumptions. The impact of operating a system without an interferometer on the amplitudes and phases of the observed wave components is considered. We find that the lack of altitude resolution of the COBRA DAQ system leads to an underestimation of the amplitude of the s = 1 component of the semidiurnal tide of ∼20% during the summer months.     

Jovian periodicities (∼10 h, ∼40 min) on Ulysses’ Distant Jupiter Encounter observations around the Halloween CIR events

We analyzed data from four different instruments (HI-SCALE, URAP, SWOOPS, VHM/FGM) onboard Ulysses spacecraft (s/c) and we searched for possible evidence of Jovian emissions when the s/c approached Jupiter during the times of Halloween events (closest time approach/position to Jupiter: February 5, 2004/R = 1683 R_J,θ = ∼49°). In particular, we analyzed extensively the low energy ion measurements obtained by the HI-SCALE experiment in order to examine whether low energy ion/electron emissions show a symmetry, and whether they are observed at north high latitudes upstream from the jovian bow shock, as is known to occur in the region upstream from the south bow shock as well ( Marhavilas et al., 2001). We studied the period from October 2003 to March 2004, as Ulysses moved at distances 0.8–1.2 AU from the planet at north Jovicentric latitudes <75°, and we present here an example of characteristic Jovian periodicities in the measurements around a CIR observed by Ulysses on days ∼348–349/2003 (R = 1894 R_J,θ = 72°). We show that Ulysses observed low energy ion (∼0.055–4.7 MeV) and electron (>∼40 keV) flux and/or spectral modulation with the Jupiter rotation period (∼10 h) as well as variations with the same period in solar wind parameters, radio and magnetic field directional data. In addition, characteristic strong ∼40 min periodic variations were found superimposed on the ∼10 h ion spectral modulation. Both the ∼10 h and ∼40 min ion periodicities in HI-SCALE measurements were present in several cases during the whole period examined (October 2003 to March 2004) and were found to be more evident during some special conditions, for instance during enhanced fluxes around the start (forward shock) and the end (reverse shock) of CIRs. We infer that the Jovian magnetosphere was triggered by the impact of the CIRs, after the Halloween events, and it was (a) a principal source of forward and reverse shock-associated ion flux structures and (b) the cause of generation of ∼10 h quasi-periodic magnetic field and plasma modulation observed by Ulysses at those times.     

Nighttime thermospheric meridional neutral winds inferred from ionospheric h′F and hpF2 data

Nighttime thermospheric meridional winds aligned to the magnetic meridian have been inferred using h′F and hpF2 ionosonde data taken from two equatorial stations, Manaus (2.9°S, 60.0°W, dip latitude 6.0°N) and Palmas (10.17°S, 48.2°W, dip latitude 6.2°S), and one low-latitude station, Sao Jose dos Campos (23.21°S, 45.86°W, dip latitude 17.26°S), during geomagnetic quiet days of August and September, 2002. Using an extension of the ionospheric servo model and a simple formulation of the diffusive vertical drift velocity, the magnetic meridional component of the thermospheric neutral winds is inferred, respectively, at the peak (hpF2) and at the base (h′F) heights of the F region over Sao Jose dos Campos. An approach has been included in the models to derive the effects of the electrodynamic drift over Sao Jose dos Campos from the time derivative of hpF2 and h′F observed at the equatorial stations. The magnetic meridional winds inferred from the two methods, for the months of August and September, are compared with winds calculated using the HWM-90 model and with measurements from Fabry–Perot technique. The results show varying agreements and disagreements. Meridional winds calculated from hpF2 ionospheric data (servo model) may produce errors of about 59 m/s, whereas the method calculated from the F-region base height (h′F) ionospheric data gives errors of about 69 m/s during the occurrence of equatorial spread-F.     

Response of the mid-latitude D-region ionosphere to the total solar eclipse of 22 July 2009 studied using VLF signals in South Korean peninsula

In this paper, we analyze VLF signals received at Busan to study the the D-region changes linked with the solar eclipse event of 22 July 2009 for very short (∼390 km) transmitter–receiver great circle path (TRGCP) during local noon time 00:36–03:13 UT (09:36–12:13 KST). The eclipse crossed south of Busan with a maximum obscuration of ∼84%. Observations clearly show a reduction of ∼6.2 dB in the VLF signal strength at the time of maximum solar obscuration (84% at 01:53 UT) as compared to those observed on the control days. Estimated values of change in Wait ionospheric parameters: reflection height (h′) in km and inverse scale height parameter (β) in km⁻¹ from Long Wave Propagation Capability (LWPC) model during the maximum eclipse phase as compared to unperturbed ionosphere are 7 km and 0.055 km⁻¹, respectively. Moreover, the D-region electron density estimated from model computation shows 95% depletion in electron density at the height of ∼71 km. The reflection height is found to increase by ∼7 km in the D-region during the eclipse as compared to those on the control days, implying a depletion in the Lyman-α flux by a factor of ∼7. The present observations are discussed in the light of current understanding on the solar eclipse induced D-region dynamics.     

Forward and reverse CIR shocks at 4–5 AU: Ulysses

In this article, we study fast shocks at CIR boundaries during an extended interval of 15 consecutive major high speed solar wind streams in 1992–1993. Ulysses was 4–5 AU from the sun. The Abraham-Schrauner shock normal method and the Rankine-Hugoniot relations were used to determine fast shock directions and speeds. Out of 33 potential CIR shocks, 14 were determined to be fast forward shocks (FSs) and 14 were fast reverse shocks (RSs). Of the remaining 5 events, 2 were forward waves and 3 were reverse waves. CIR edges at latitudes below ∼30^o were, for the most part, bounded by fast magnetosonic shocks. The forward shocks were generally quasi-perpendicular (average θ_nBo = 67^o). The reverse shocks were more oblique (average θ_nBo = 52^o), but they extended to all angles. Both FSs and RSs had magnetosonic Mach numbers ranging from 1 to 5 or 6. The average Mach numbers were 2.4 and 2.6 for FSs and RSs, respectively. The shock Mach numbers were noted to generally decrease with increasing latitude. The non-shock events or waves were noted to occur preferentially at high (∼−30° to −35°) heliolatitudes where stream-stream interactions were presumably weaker. These results are consistent with expectations, indicating the general accuracy of the Abraham-Schrauner technique.     

Investigation of ionospheric irregularities and scintillation using TEC at high latitude

The occurrence of ionospheric irregularities at high latitudes, with dimensions of several kms down to decameter scale size shows strong correlation with geomagnetic disturbance, season and solar activity. Transionospheric radio waves propagating through these irregularities experience rapid random fluctuations in phase and/or amplitude of the signal at the receiver, termed scintillation, which can degrade GNSS services. Thus, investigation and prediction of this scintillation effect is very important. To investigate such scintillation effects, a GISTM (GPS Ionospheric Scintillation and TEC Monitoring) NovAtel dual frequency (L1/L2) GPS receiver has been installed at Trondheim, Norway (63.41°$63.41°$ N, 10.4°$10.4°$ E), capable of collecting scintillation indices at a 1 min rate as well as the raw data (phase and intensity) of the satellite signals at a 50 Hz sampling rate and TEC (Total Electron Content) at a 1 Hz rate. Many researchers have reported that both phase and amplitude scintillation is closely associated with TEC fluctuations or associated with a significant developing enhancement or depletion in the TEC. In this study, a novel analogous phase index is developed which provides samples at a 1 min rate. Generally the scintillation indices can help in estimating the irregularity scintillation effect at a one minute rate, but such procedures are time consuming if DFTs of the phase and/or amplitude at a 50 Hz data are required. In this study, instead, this analogous phase index is estimated from 1 Hz rate TEC values obtained from the raw signals and is then compared for weak, moderate and strong scintillation at Trondheim for one year of data collected from the installed GPS receiver. The spectral index of the irregularities (that is the inverse power law of their spatial spectrum) is determined from the resultant phase scintillation psd. The correlations of the scintillation indices and spectral indices with the analogous phase index have been investigated under different geomagnetic conditions (represented by the Kp index) and an approximate linear correlation of phase scintillation with the analogous phase index was found. Then a principal advantage of this index is that it achieves this correlation without requiring a high sampling data rate and the need for DFTs. Thus, the index seems a good candidate for developing a simple means of ionospheric scintillation prediction which could also be utilized in the development of alerts using regional mappings.     

Gravity wave propagation studies using the Indian MST radar observations

The Indian MST radar facility at Gadanki (13.5°N, 79.2°E) has been utilised to study the propagation of gravity waves from the troposphere/lower stratosphere to the mesosphere and their interaction with the radar backscattered signal variations. The main objective is to correlate vertically propagating gravity waves derived from the tropospheric velocity fields with the dynamics of mesospheric scattering centres. The tropospheric wind velocities and signal strengths over the entire height range have been subjected to power spectral and wavelet analysis to determine the predominant wave periods/amplitudes and the coupling between the lower atmosphere and mesosphere. Results show that (a) the gravity waves are clearly detectable near tropopause heights, (b) while relatively higher period gravity waves (20–50 min) interact with mesospheric scattering centres, the lower period waves (<20 min) are absorbed in the troposphere itself, (c) the mesospheric scattering layers are affected by gravity waves of complementary periods.     

Earthquake induced dynamics at the ionosphere in presence of magnetic storm

The modifications induced in the dynamics of the ionosphere by the major Japan earthquake (EQ) of March 11, 2011 (epicenter at 38.322°N, 142.369°E, M = 8.9) in presence of a magnetic storm are examined by mapping latitudinal variations of F-layer ionization density (NmF2) from 22 stations covering the epicenter zone. The changes forced into the Total Electron Content (TEC) by the major EQ in the magnetic storm ambiance are also examined from the GPS data collected at Guwahati (26° 10′ N, 91° 45’ E), situated in the major fault system of East Asia. The contributions of pre-seismic electric field as well as of magnetic storm time electric field in the observed density variations are brought into the ambit of discussion. The influence of lower atmosphere in shaping TEC features during the study case is highlighted. The effects of solar activity on density variations during such complex ambiances are also addressed.