Study of daytime vertical E × B drift velocities inferred from ground-based magnetometer observations of ΔH,at low latitudes under geomagnetically disturbed conditions

In this study, 30 storm sudden commencement (SSC) events during the period 2001–2007 for which daytime vertical E × B drift velocities from JULIA radar, Jicamarca (geographic latitude 11.91°S, geographic longitude 283.11°E, 0.81°N dip latitude), Peru and ΔH component of geomagnetic field measured as the difference between the magnitudes of the horizontal (H) components between two magnetometers deployed at two different locations Jicamarca (geographic latitude 11.91°S, geographic longitude 283.11°E, 0.81°N dip latitude) and Piura (geographic latitude 5.21°S, geographic longitude 279.41°E, 6.81°N dip latitude), in Peru, were considered. It is observed that a positive correlation exists between peak value of daytime vertical E × B drift velocity and peak value of ΔH for the three consecutive days of SSC. A qualitative analysis made after selecting the peak values of daytime vertical E × B drift velocity and ΔH showed that 57% of the events have daytime vertical E × B drift velocity peak in the magnitude range 20–30 m/s and 63% of the events have ΔH peak in the range 80–100 nT. The maximum probable (45%) range of time of occurrence of peak value for both vertical E × B drift velocity and ΔH during the daytime hours were found to be the same, i.e., 10:00–12:00 LT. A strong positive correlation was also found to exist between the daytime vertical E × B drift velocity and ΔH for all the three consecutive days of SSC, for all the events considered. To establish a quantitative relationship between day time vertical E × B drift velocity and ΔH, linear and polynomial (order 2 and 3) regression analysis (Least Square Method (LSM)) were carried out, considering the fully disturbed day after the commencement of the storm as ‘disturbed period’ for the SSC events selected for analysis. The formulae indicating the relationship between daytime vertical E × B drift velocity and ΔH, for the ‘disturbed periods’, obtained through the regression analysis were verified using the JULIA radar observed E × B drift velocity for 3 selected events. Root Mean Square (RMS) error analysis carried out for each case suggest that polynomial regression (order 3) analysis provides a better agreement with the observations from among the linear, polynomial (order 2 and 3) analysis.     

Rare observation of daytime whistlers at very low latitude (L = 1.08)

The source region and propagation mechanism of low latitude whistlers (Geomag. lat. <30°) have puzzled scientific community for last many decades. In view of recent reports, there is consensus on the source region of low latitude whistlers in the vicinity of the conjugate point. But the plausible conditions of ionospheric medium through which they travel are still uncertain. In addition to that, the whistlers in daytime are never observed at geomagnetic latitudes less than 20°. Here, for the first time, we present a rare observations of whistlers during sunlit hours from a very low-latitude station Allahabad (Geomag. Lat: 16.79°N, L = 1.08) in India on 04 February 2011. More than 90 whistlers are recorded during 1200–1300 UT during which the whole propagation path from lightning source region to whistler observation site is under sunlit. The favorable factors that facilitated the whistlers prior to the sunset are investigated in terms of source lightning characteristics, geomagnetic and background ionospheric medium conditions. The whistler activity period was found to be geomagnetically quiet. However, a significant suppression in ionospheric total electron content (TEC) compared to its quiet day average is found. This shows that background ionospheric conditions may play a key role in low latitude whistler propagation. This study reveals that whistlers can occur under sunlit hours at latitudes as low as L = 1.08 when the source lightning and ionospheric medium characteristics are optimally favorable.     

Multi precursors analysis associated with the powerful Ecuador (MW = 7.8) earthquake of 16 April 2016 using Swarm satellites data in conjunction with other multi-platform satellite and ground data

After DEMETER satellite mission (2004–2010), the launch of the Swarm satellites (Alpha (A), Bravo (B) and Charlie (C)) has created a new opportunity in the study of earthquake ionospheric precursors. Nowadays, there is no doubt that multi precursors analysis is a necessary phase to better understand the LAIC (Lithosphere Atmosphere Ionosphere Coupling) mechanism before large earthquakes. In this study, using absolute scalar magnetometer, vector field magnetometer and electric field instrument on board Swarm satellites, GPS (Global Positioning System) measurements, MODIS-Aqua satellite and ECMWF (European Centre for Medium-Range Weather Forecasts) data, the variations of the electron density and temperature, magnetic field, TEC (Total Electron Content), LST (Land Surface Temperature), AOD (Aerosol Optical Depth) and SKT (SKin Temperature) have been surveyed to find the potential seismic anomalies around the strong Ecuador (M_w = 7.8) earthquake of 16 April 2016. The four solar and geomagnetic indices: F10.7, D_st, K_p and a_p were investigated to distinguish whether the preliminary detected anomalies might be associated with the solar-geomagnetic activities instead of the seismo-ionospheric anomalies. The Swarm satellites (A, B and C) data analysis indicate the anomalies in time series of electron density variations on 7, 11 and 12 days before the event; the unusual variations in time series of electron temperature on 8 days preceding the earthquake; the analysis of the magnetic field scalar and vectors data show the considerable anomalies 52, 48, 23, 16, 11, 9 and 7 days before the main shock. A striking anomaly is detected in TEC variations on 1 day before earthquake at 9:00 UTC. The analysis of MODIS-Aqua night-time images shows that LST increase unusually on 11 days prior to main shock. In addition, the AOD variations obtained from MODIS measurements reach the maximum value on 10 days before the earthquake. The SKT around epicentral region presents anomalous higher value about 40 days before the earthquake. It should be noted that the different lead times of the observed anomalies could be acknowledged based on a reasonable LAIC earthquake mechanism. Our results emphasize that the Swarm satellites measurements play an undeniable role in progress the studies of the ionospheric precursors.     

Comparison of GPS-TEC measurements with NeQuick2 and IRI model predictions in the low latitude East African region during varying solar activity period (1998 and 2008–2015)

This paper examines the performances of NeQuick2, the latest available IRI-2016, IRI-2012 and IRI-2007 models in describing the monthly and seasonal mean total electron content (TEC) over the East African region. This is to gain insight into the success of the various model types and versions at characterizing the ionosphere within the equatorial ionization anomaly. TEC derived from five Global Positioning System (GPS) receivers installed at Addis Ababa (ADD, 5.33°N, 111.99°E Geog.), Asab (ASAB, 8.67°N, 116.44°E Geog.), Ambo (ABOO, 5.43°N, 111.05°E Geog.), Nairobi (RCMN, ?4.48°N, 108.46°E Geog.) and Nazret (NAZR, 4.78°N, 112.43°E Geog.), are compared with the corresponding values computed using those models during varying solar activity period (1998 and 2008–2015). We found that different models describe the equatorial and anomaly region ionosphere best depending on solar cycle, season and geomagnetic activity levels. Our results show that IRI-2016 is the best model (compared to others in terms of discrepancy range) in estimating the monthly mean GPS-TEC at NAZR, ADD and RCMN stations except at ADD during 2008 and 2012. It is also found that IRI-2012 is the best model in estimating the monthly mean TEC at ABOO station in 2014. IRI show better agreement with observations during June solstice for all the years studied at ADD except in 2012 where NeQuick2 better performs. At NAZR, NeQuick2 better performs in estimating seasonal mean GPS-TEC during 2011, while IRI models are best during 2008–2009. Both NeQuick2 and IRI models underestimate measured TEC for all the seasons at ADD in 2010 but overestimate at NAZR in 2009 and RCMN in 2008. The periodic variations of experimental and modeled TEC have been compared with solar and geomagnetic indices at ABOO and ASAB in 2014 and results indicate that the F10.7 and sunspot number as indices of solar activity seriously affects the TEC variations with periods of 16–32?days followed by the geomagnetic activity on shorter timescales (roughly periods of less than 16?days). In this case, NeQuick2 derived TEC shows better agreement with a long term period variations of GPS-TEC, while IRI-2016 and IRI-2007 show better agreement with observations during short term periodic variations. This indicates that the dependence of NeQuick2 derived TEC on F10.7 is seasonal. Hence, we suggest that representation of geomagnetic activity indices is required for better performance over the low latitude region.