The deterministic chaotic behavior and dynamical complexity of the space plasma dynamical system over Nigeria are analyzed in this study and characterized. The study was carried out using GPS (Global Positioning System) TEC (Total Electron Content) time series, measured in the year 2011 at three GPS receiver stations within Nigeria, which lies within the equatorial ionization anomaly region. The TEC time series for the five quietest and five most disturbed days of each month of the year were selected for the study. The nonlinear aspect of the TEC time series was obtained by detrending the data. The detrended TEC time series were subjected to various analyses for phase space reconstruction and to obtain the values of chaotic quantifiers like Lyapunov exponents, correlation dimension and also Tsallis entropy for the measurement of dynamical complexity. The observations made show positive Lyapunov exponents (LE) for both quiet and disturbed days, which indicates chaoticity, and for different days the chaoticity of the ionosphere exhibits no definite pattern for either quiet or disturbed days. However, values of LE were lower for the storm period compared with its nearest relative quiet periods for all the stations. The monthly averages of LE and entropy also show no definite pattern for the month of the year. The values of the correlation dimension computed range from 2.8 to 3.5, with the lowest values recorded at the storm period of October 2011. The surrogate data test shows a significance of difference greater than 2 for all the quantifiers. The entropy values remain relatively close, with slight changes in these values during storm periods. The values of Tsallis entropy show similar variation patterns to those of Lyapunov exponents, with a lot of agreement in their comparison, with all computed values of Lyapunov exponents correlating with values of Tsallis entropy within the range of 0.79 to 0.81. These results show that both quantifiers can be used together as indices in the study of the variation of the dynamical complexity of the ionosphere. The results also show a strong play between determinism and stochasticity. The behavior of the ionosphere during these storm and quiet periods for the seasons of the year are discussed based on the results obtained from the chaotic quantifiers.

The equatorial Congo has been recognized as the most active lightning chimney region in the Globe. Although the perturbation of tropospheric thunderstorms on the lower ionosphere has been noticed in the middle latitudes through their transient lightning electric fields or convective gravity waves, the effects on equatorial ionosphere and the horizontal extent of this perturbation remains a mystery because of the difficulties in extracting the effects due to the sporadic nature of the equatorial ionosphere. Here we present observational results showing solid evidence of deviations in ionospheric total electron content (TEC) and its direction of propagation associated with thunderstorms using the method of polynomial filtering, by utilizing the TEC measured from equatorial Global Positioning System (GPS) Receiver stations along the West African region-Congo Basin. The TEC deviations due to the thunderstorms were found to be mostly propagated in a specific direction from the point of the event, with the highest absolute peak TEC at ~±1.5 TECUs. The internal dynamics of the equatorial ionosphere have been found to be suppressed by large thunderstorm effects during the daytime, with negligible impact at night.

The quest to find an index for proper characterization and description of the dynamical response of the ionosphere to external influences and its various internal irregularities has led to the study of the day-to-day variations of the chaoticity and dynamical complexity of the ionosphere. This study was conducted using Global Positioning System (GPS) total electron content (TEC) time series, measured in the year 2011, from five GPS receiver stations in Nigeria, which lies within the equatorial ionization anomaly region. The non-linear aspects of the TEC time series were obtained by detrending the data. The detrended TEC time series were subjected to various analyses to obtain the phase space reconstruction and to compute the chaotic quantifiers, which are Lyapunov exponents LE, correlation dimension, and Tsallis entropy, for the study of dynamical complexity. Considering all the days of the year, the daily/transient variations show no definite pattern for each month, but day-to-day values of Lyapunov exponents for the entire year show a wavelike semiannual variation pattern with lower values around March, April, September and October. This can be seen from the correlation dimension with values between 2.7 and 3.2, with lower values occurring mostly during storm periods, demonstrating a phase transition from higher dimension during the quiet periods to lower dimension during storms for most of the stations. The values of Tsallis entropy show a similar variation pattern to that of the Lyapunov exponent, with both quantifiers correlating within the range of 0.79 to 0.82. These results show that both quantifiers can be further used together as indices in the study of the variations of the dynamical complexity of the ionosphere. The presence of chaos and high variations in the dynamical complexity, even in quiet periods in the ionosphere, may be due to the internal dynamics and inherent irregularities of the ionosphere which exhibit non-linear properties. However, this inherent dynamics may be complicated by external factors like geomagnetic storms. This may be the main reason for the drop in the values of the Lyapunov exponent and Tsallis entropy during storms. The dynamical behaviour of the ionosphere throughout the year, as described by these quantifiers, was discussed in this work.

Perturbations from thunderstorms can play a notable role in the dynamics of the ionosphere. In this work, ionospheric perturbation effects due to thunderstorms were extracted and studied. Thunderstorm-associated lightning activities and their locations were detected by the World-Wide Lightning Location Network (WWLLN). The mechanical components of ionospheric perturbations due to thunderstorms were extracted from the total electron content (TEC), which was measured at selected thunderstorm locations using the polynomial filtering method. Further analyses were conducted using wavelet analysis and Discrete Fourier Transform (DFT) to study the frequency modes and periodicities of TEC deviation. It was revealed that the highest magnitudes of TEC deviations could reach up to ~2.2 TECUs, with dominant modes of frequency in the range of ~0.2 mHz to ~1.2 mHz, falling within the gravity wave range and the second dominant mode in the acoustic range of >1 mHz to <7.5 mHz. Additionally, a 20–60 min time delay was observed between the sprite events, the other high-energy electrical discharges, and the time of occurrence at the highest peak of acoustic-gravity wave perturbations extracted from TEC deviations. The possible mechanism responsible for this phenomenon is further proposed and discussed.