Unusually strong magnetic fields in Titan’s ionosphere: T42 case study

Observations of unusually large magnetic fields in the ionosphere indicate periods of maximum stress on Titan’s ionosphere and potentially of the strongest loss rates of ionospheric plasma. During Titan flyby T42, the observed magnetic field attained a maximum value of 37 nT between an altitude of 1200 and 1600 km, about 20 nT stronger than on any other Titan pass and close to five times greater in magnetic pressure. The strong fields occurred near the corotation-flow terminator rather than at the sub-flow point, suggesting that the flow which magnetized the ionosphere was from a direction far from corotation and possibly towards Saturn. Extrapolation of solar wind plasma conditions from Earth to Saturn using the University of Michigan MHD code predicts an enhanced solar wind dynamic pressure at Saturn close to this time. Cassini’s earlier exits from Saturn’s magnetosphere support this prediction because the Cassini Plasma Spectrometer instrument saw a magnetopause crossing three hours before the strong field observation. Thus it appears that Titan’s ionosphere was magnetized when the enhanced solar wind dynamic pressure compressed the Saturnian magnetosphere, and perhaps the magnetosheath magnetic field, against Titan. The solar wind pressure then decreased, leaving a strong fossil field in the ionosphere. When observed, this strong magnetic flux tube had begun to twist, further enhancing its strength.     

Two-dimensional ionospheric total electron content map (TEC) seismo-ionospheric anomalies through image processing using principal component analysis

This paper examines China’s Wenchuan Earthquake of 12 May 2008 (UTC) (M_w = 7.9) using principal component analysis and image processing of the global ionospheric map (GIM) for the region. Transforms are conducted for 4, 8, and 9 May 2008. The GIMs are subdivided into 100 (36° in Long. and 18° in Lat.) smaller maps. The smaller maps (71 × 71 pixels) form the transform matrices of corresponding dimensions (2 × 1) through image processing. The transform allows for principle eigenvalues to be assigned to TEC anomalies for May 8 and 9. These may represent the seismo-ionospheric signature described by Zhao et al. (2008). The May 4 result shows no evidence of TEC anomalies. These results are in keeping with the findings of Liu et al. (2009). It is evident in this research that PCA could have the capacity to detect both the seismo-ionospheric signature and determine the approximate location of an earthquake’s epicenter prior to nucleation.     

Shape and motion estimate of LEO debris using light curves

We have derived a tri-axial ellipsoidal model of an LEO object, a Cosmos 2082 rocket body, including its rotational axis direction, rotation period, precession, and a compositional parameter, using only light curve data from an optical telescope. The brightness of the object was monitored for two days and least-squares fitting was used to determine these values. The derived axial ratios are 100:18:18, the coordinates of the rotational axis direction on the celestial sphere are R.A. = 305.8° and Dec. = 2.6°, and its observed average rotation period is 41 s. When precession is considered, its amplitude and precession period are 30.5° and 29.4 min. These results show that optical light curve data are sufficient to determine the shapes and the motions of LEO objects.     

Determination of station positions and velocities from laser ranging observations to Ajisai,Starlette and Stella satellites

The positions and velocities of the four Satellite Laser Ranging (SLR) stations: Yarragadee (7090), Greenbelt (7105), Graz (7839) and Herstmonceux (7840) from 5-year (2001–2005) SLR data of low orbiting satellites (LEO): Ajisai, Starlette and Stella were determined. The orbits of these satellites were computed from the data provided by 20 SLR stations. All orbital computations were performed by means of NASA Goddard’s GEODYN-II program. The geocentric coordinates were transformed to the topocentric North–South, East–West and Vertical components in reference to ITRF2005. The influence of the number of normal points per orbital arc and the empirical acceleration coefficients on the quality of station coordinates was studied. To get standard deviation of the coordinates determination lower than 1 cm, the number of the normal points per site had to be greater than 50. The computed positions and velocities were compared to those derived from LAGEOS-1/LAGEOS-2 data. Three parameters were used for this comparison: station coordinates stability, differences from ITRF2005 positions and velocities. The stability of coordinates of LEO satellites is significantly worse (17.8 mm) than those of LAGEOS (7.6 mm), the better results are for Ajisai (15.4 mm) than for Starlette/Stella (20.4 mm). The difference in positions between the computed values and ITRF2005 were little bit worse for Starlette/Stella (6.6 mm) than for LAGEOS (4.6 mm), the results for Ajisai were five times worse (29.7 mm) probably due to center of mass correction of this satellite. The station velocities with some exceptions were on the same level (≈1 mm/year) for all satellites. The results presented in this work show that results from Starlette/Stella are better than those from Ajisai for station coordinates determination. We can applied the data from LEO satellites, especially Starlette and Stella for determination of the SLR station coordinates but with two times lower accuracy than when using LAGEOS data.     

Investigation of GPS-TEC measurements using ANN method indicating seismo-ionospheric anomalies around the time of the Chile (Mw = 8.2) earthquake of 01 April 2014

In 94 km NW of Iquique in Chile (19.610°S, 70.776°W) a powerful earthquake of M_w = 8.2 took place at 23:46:47 UTC (20:46:47 LT) on April 01, 2014. Using GPS-TEC (Total Electron Content) measurements, potential unusual variations around the time and location of the Chile earthquake have been detected based on the median and Artificial Neural Network (ANN) methods. The indices D_st, K_p, A_p and F10.7 were used to distinguish pre-earthquake anomalies from the other anomalies related to the solar-geomagnetic activities. Using the median method, striking anomalies in time series of TEC data are observed 4 days before the earthquake at 14:00 and 16:00 UTC. The ANN method detected a number of anomalies, 4 (02:00 and 16:00 UTC) and 13 (24:00 UTC) days preceding the earthquake. The results indicate that the ANN method due to its capability of non linear learning is quite promising and deserves serious attention as a robust predictor tool for seismo-ionospheric anomalies detection.     

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.     

An icy mineralogy package (IMP) for in-situ studies of Titan’s surface

We describe the scientific case for and preliminary design of an instrument whose primary goal is to determine the chemistry (element abundance) and mineralogy (compound identity and abundance) of Titan’s surface using a combination of energy dispersive X-ray fluorescence spectroscopy (EDXRF) and X-ray diffraction (XRD). XRD is capable of identifying any crystalline substance present on Titan’s surface at relative abundances greater than ∼1 wt%, allowing unambiguous identification of, for example, structure I and II clathrates (even in the presence of ice), and various organic solids, which may include C₂H₂, C₂H₄, C₄H₂, HCN, CH₃CN, HC₃N, and C₄N₂). The XRF component of the instrument will obtain elemental abundances for 16 < Z < 60 with minimum detection limits better than 10 ppm (including detection of atmospheric noble gas isotopes), and may achieve detection limits of 0.01–1% for lighter elements down to Z = 6 (carbon). The instrument is well suited to integration with other analytical tools as part of a light-weight surface chemistry and mineralogy package. Although considerably less sensitive to elemental abundance than GC–MS (10⁻² vs. 10⁻⁸) it is likely to be significantly lighter (<0.5 kg vs. 10 kg).     

Mars Sample Return: Do Australians trust NASA?

Mars Sample Return (MSR) represents an important scientific goal in space exploration. Any sample return mission will be extremely challenging from a scientific, economic and technical standpoint. But equally testing, will be communicating with a public that may have a very different perception of the mission. A MSR mission will generate international publicity and it is vital that NASA acknowledge the nature and extent of public concern about the mission risks and, perhaps equally importantly, the public’s confidence in NASA’s ability to prepare for and manage these risks. This study investigated the level of trust in NASA in an Australian population sample, and whether this trust was dependent on demographic variables. Participants completed an online survey that explored their attitudes towards NASA and a MSR mission. The results suggested that people believe NASA will complete the mission successfully but have doubts as to whether NASA will be honest when communicating with the public. The most significant finding to emerge from this study was that confidence in NASA was significantly (p < 0.05) related to the respondent’s level of knowledge regarding the risks and benefits of MSR. These results have important implications for risk management and communication.     

Monthly mean foF2 model for an African low-latitude station and comparison with IRI

We have employed the hourly values of the ionospheric F-region critical frequency (foF2) obtained from Ouagadougou ionosonde, Burkina Faso (geographic coordinates 12° N, 1.8° W) during the interval of 1985–1995 (solar cycle 22) and solar radio flux of 10 cm wavelength (F10.7) to develop a local model (LM) for the African low-latitude station. The model was developed from regression analysis method, using the two-segmented regression analysis. We validated LM with foF2 data from Korhogo observatory, Cote d’Ivorie (geographical coordinates 9.3° N, 5.4° W). LM as well as the International Reference Ionosphere (IRI) agrees well with observations. LM gave some improvement on the IRI-predicted foF2 values at the sunrise (06 LT) at all solar flux levels and in all seasons except June solstice. The performance of the models at the representing the salient features of the equatorial foF2 was presented. Considering daytime and nighttime performances, LM and IRI are comparable in low solar activity (LSA), LM performed better than IRI in moderate solar activity (MSA), while IRI performed better than LM in high solar activity (HSA). CCIR has a root mean square error (r.m.s.e), which is only 0.10 MHz lower than that of LM while LM has r.m.s.e, which is about 0.05 MHz lower than that of URSI. In general, our result shows that performance of IRI, especially the CCIR option of the IRI, is quite comparable with the LM. The improved performance of IRI is a reflection of the numerous contributions of ionospheric physicists in the African region, larger volume of data for the IRI and the diversity of data sources, as well as the successes of the IRI task force activities.     

Comparison of mesopause region meteor radar winds,medium frequency radar winds and low frequency drifts over Germany

During 2004 and 2005 measurements of mesospheric/lower thermospheric (80–100 km) winds have been carried out in Germany using three different ground-based systems, namely a meteor radar (36.2 MHz) at the Collm Observatory (51.3°N, 13°E), a MF radar (3.18 MHz) at Juliusruh (54.6°N, 13.4°E) and the LF D1 measurements using a transmitter (177 kHz) at Zehlendorf near Berlin and receivers at Collm with the reflection point at 52.1°N, 13.2°E. This provides the possibility of comparing the results of different radar systems in nearly the same measuring volume. Meteor radar winds are generally stronger than the winds observed by MF and especially by LF radars. This difference is small near 80 km but increases with height. The difference between meteor radar and medium frequency radar winds is larger during winter than during summer, which might indicate an indirect influence of gravity waves on spaced antenna measurements.     

Radio beacons/IMU integrated navigation for Mars entry

High precision entry navigation capability is essential for future Mars pinpoint landing missions, together with the entry guidance and aerodynamic lift control. This paper addresses the issue of Mars entry navigation using inertial measurement unit (IMU) and orbiting or surface radiometric beacons. The range and Doppler information sensed from orbiting or surface radio beacons and the entry vehicle state information derived from IMU are integrated in Unscented Kalman filter to correct the inertial constant bias and suppress the navigation measurement noise. Computer simulations show that the integrated navigation algorithm proposed in this paper can achieve 50 m position error and 2 m/s velocity error, which satisfies the need of future pinpoint Mars landing missions.     

Effective non-vertical and apparent cutoff rigidities for a cosmic ray latitude survey from Antarctica to Italy in minimum of solar activity

In this paper we will report the results of the computation of cutoff rigidities of vertical and non-vertical incident cosmic ray particles. Non-vertical effective cutoff rigidities have been computed by tracing particle trajectories through the “real” geomagnetic magnetic field comprising the International Geomagnetic Reference Field model (IGRF95, IAGA Division 5 Working Group 8, 1996: Sabaka, T.J., Langel, R.A., Baldwin, R.T., Conrad, J.A. The geomagnetic field, 1900–1995, including the large scale fields from magnetospheric sources and NASA candidate models for the 1995 IGRF revision. J. Geomag. Geoelect. 49, 157–206, 1997.) and the Tsyganenko [Tsyganenko, N.A. A magnetospheric magnetic field model with a warped tail current sheet. Planet. Space Sci. 37, 5–20, 1989.] magnetosphere model. The computation have been done for the backward route (from Antarctica to Italy) of the Italian Antarctic ship survey 1996–1997, for geographic points corresponding to the daily average coordinates of the ship; for zenith angles 15°, 30°, 45° and 60°, and azimuth angles from 0° to 360° in steps of 45°. By means of the obtained non-vertical cutoffs the apparent cutoff rigidities have been calculated. The information on integral multiplicities of secondary neutrons detected by the neutron monitor in dependence of the zenith angle of incoming primary cosmic ray particles have also been used. This information is based on the theoretical calculations of meson-nuclear cascades of primary protons with different rigidities arriving to the Earth’s atmosphere at the zenith angles of 0°, 15°, 30°, 45°, 60° and 75°. The difference between the computed apparent and vertical cutoff rigidities reaches ∼1 GV at rigidities >7–8 GV. At rigidities of 10–16 GV, the difference between the apparent and vertical cutoff rigidities is larger than that obtained earlier by Clem et al. [Clem, J.M., Bieber, J.W., Duldig, M., Evenson, P., Hall, D., Humble, J.E. Contribution of obliquely incident particles to neutron monitor counting rate. J. Geophys. Res. 102, 26919–26926, 1997.] and Dorman et al. [Dorman, L.I., Villoresi, G., Iucci, N., Parisi, M., Tyasto, M.I., Danilova, O.A., Ptitsyna, N.G. Cosmic ray survey to Antarctica and coupling functions for neutron component near solar minimum (1996–1997), 3. Geomagnetic effects and coupling functions. J. Geophys. Res. 105, 21047–21056, 2000.].     

Plasma and electromagnetic wave simulations of meteors

Every day billions of meteoroids impact and disintegrate in the Earth’s atmosphere. Current estimates for this global meteor flux vary from 2000 to 200,000 tons per year, and estimates for the average velocity range between 10 km/s and 70 km/s. The basic properties of this global meteor flux, such as the average mass, velocity, and chemical composition remain poorly constrained. We believe much of the mystery surrounding the basic parameters of the interplanetary meteor flux exists for the following reason, the unknown sampling characteristics of different radar meteor observation techniques, which are used to derive or constrain most models. We believe this arises due to poorly understood radio scattering characteristics of the meteor plasma, especially in light of recent work showing that plasma turbulence and instability greatly influences meteor trail properties at every stage of evolution. We present our results on meteor plasmas simulations of head echoes using particle in cell (PIC) ions, which show that electric fields strongly influence early stage meteor plasma evolution, by accelerating ions away from the meteoroid body. We also present the results of finite difference time domain electromagnetic simulations (FDTD), which can calculate the radar cross section of the simulated meteor plasmas. These simulations have shown that the radar cross section depends in a complex manner on a number of parameters. These include the angle between radar and meteor entry, a large dependence on radar frequency, which shows that for a given meteor plasma size and density, the reflectivity as a function of probing radar frequency varies, but typically peaks below 100 MHz.     

MCAV/IMU integrated navigation for the powered descent phase of Mars EDL

Pin-point landing is considered as a key technology for future manned Mars landing and Mars base missions. The traditional inertial navigation system (INS) based guidance, navigation and control (GNC) mode used in the Mars entry, descent and landing (EDL) phase has no ability to achieve the precise and safe Mars landing, so novel EDL GNC methodologies should be investigated to meet this goal. This paper proposes the MCAV/IMU integrated navigation scheme for the powered descent phase of Mars EDL. The Miniature Coherent Altimeter and Velocimeter (MCAV) is adopted to correct the inertial bias and drift and improve the performance of integrated navigation. Altitude and velocity information derived from MCAV and the lander’s state information sensed by inertial measurement unit (IMU) are integrated in extended Kalman filter algorithm. The validity of the proposed navigation scheme is confirmed by computer simulation.     

Spin rate and spin axis orientation of LARES spectrally determined from Satellite Laser Ranging data

Satellite Laser Ranging (SLR) is a powerful technique able to measure spin rate and spin axis orientation of the fully passive, geodetic satellites. This work presents results of the spin determination of LARES – a new satellite for testing General Relativity. 529 SLR passes measured between February 17 and June 9, 2012, were spectrally analyzed. Our results indicate that the initial spin frequency of LARES is f₀ = 86.906 mHz (RMS = 0.539 mHz). A new method for spin axis determination, developed for this analysis, gives orientation of the axis at RA = 12^h22^m48^s (RMS = 49^m), Dec = −70.4° (RMS = 5.2°) (J2000.0 celestial reference frame), and the clockwise (CW) spin direction. The half-life period of the satellite’s spin is 214.924 days and indicates fast slowing down of the spacecraft.     

Use of principal component analysis in the identification of the spatial pattern of an ionospheric total electron content anomalies after China’s May 12, 2008, M = 7.9 Wenchuan earthquake

This paper uses principal component analysis (PCA) to determine the spatial pattern of total electron content (TEC) anomalies in the ionosphere post the China’s Wenchuan Earthquake of 12 May, 2008 (UTC) (M_w = 7.9). PCA and image processing are applied to the global ionospheric map (GIM) with transforms conducted for the time period 08:00–10:00 UT on 12 May 2008. Results show that at a height of approximately 200 km the anomaly is widespread and less intense; however, it becomes more localized with height reaching maximum intensity and localization at an altitude of 300 km. The spatial distribution is remarkably similar to that reported for a TEC anomaly previously identified as a precursor anomaly on May 9, 2008 for the same time period. Potential causes of the results are discussed with emphasis given to vertical acoustic gravity waves based on the spatial pattern identified.     

Landforms of Europa and selection of landing sites

Three major features make Europa a unique scientific target for a lander-oriented interplanetary mission: (1) the knowledge of the composition of the surface of Europa is limited to interpretations of the spectral data, (2) a lander could provide unique new information about outer parts of the solar system, and (3) Europa may have a subsurface ocean that potentially may harbor life, the traces of which may occur on the surface and could be sampled directly by a lander. These characteristics of Europa bring the requirement of safe landing to the highest priority level because any successful landing on the surface of this moon will yield scientific results of fundamental importance. The safety requirements include four major components. (1) A landing site should preferentially be on the anti-Jovian hemisphere of Europa in order to facilitate the orbital maneuvers of the spacecraft. (2) A landing site should be on the leading hemisphere of Europa in order to extend the lifetime of a lander and sample pristine material of the planet. (3) Images with the highest possible resolution must be available for the selection of landing sites. (4) The terrain for landing must have morphology (relief) that minimizes the risk of landing and represents a target that is important from a scientific point of view. These components severely restrict the selection of regions for landing on the surface of Europa. After the photogeologic analysis of all Galileo images with a resolution of better than about 70 m/pixel taken for the leading hemisphere of Europa, we propose one primary and two secondary (backup) landing sites. The primary site (51.8°S, 177.2°W) is within a pull-apart zone affected by a small chaos. The first backup site (68.1°S, 196.7°W) is also inside of a pull-apart zone and is covered by images of the lower resolution (51.4 m/pixel). The second backup site (2.4°N, 181.1°W) is imaged by relatively low-resolution images (∼70 m/pixel) and corresponds to a cluster of small patches of dark and probably smooth plains that may represent landing targets of the highest scientific priority from the scientific point of view. The lack of the high-resolution images for this region prevents, however, its selection as the primary landing target.     

Debris flux comparisons from the Goldstone Radar,Haystack Radar,and Hax Radar prior,during, and after the last solar maximum

The continual monitoring of the low Earth orbit (LEO) debris environment using highly sensitive radars is essential for an accurate characterization of these dynamic populations. Debris populations are continually evolving since there are new debris sources, previously unrecognized debris sources, and debris loss mechanisms that are dependent on the dynamic space environment. Such radar data are used to supplement, update, and validate existing orbital debris models. NASA has been utilizing radar observations of the debris environment for over a decade from three complementary radars: the NASA JPL Goldstone radar, the MIT Lincoln Laboratory (MIT/LL) Long Range Imaging Radar (known as the Haystack radar), and the MIT/LL Haystack Auxiliary radar (HAX). All of these systems are highly sensitive radars that operate in a fixed staring mode to statistically sample orbital debris in the LEO environment. Each of these radars is ideally suited to measure debris within a specific size region. The Goldstone radar generally observes objects with sizes from 2 mm to 1 cm. The Haystack radar generally measures from 5 mm to several meters. The HAX radar generally measures from 2 cm to several meters. These overlapping size regions allow a continuous measurement of cumulative debris flux versus diameter from 2 mm to several meters for a given altitude window. This is demonstrated for all three radars by comparing the debris flux versus diameter over 200 km altitude windows for 3 nonconsecutive years from 1998 to 2003. These years correspond to periods before, during, and after the peak of the last solar cycle. Comparing the year to year flux from Haystack for each of these altitude regions indicate statistically significant changes in subsets of the debris populations. Potential causes of these changes are discussed. These analysis results include error bars that represent statistical sampling errors.     

Azimuthal propagation of Pc5 geomagnetic field pulsations in the southern polar cap

A statistical analysis of low frequency geomagnetic fluctuations at the two Antarctic stations Mario Zucchelli Station (geographic coordinates: 74.7°S, 164.1°E; corrected geomagnetic coordinates: 80.0°S, 306.8°E) and Dumont D’Urville (geographic coordinates: 66.7°S, 140.0°E; corrected geomagnetic coordinates: 80.4°S, 236.0°E) is shown. The analysis focuses on power spectra, coherence and phase difference between the stations, which are both located in the polar cap, with a 5-h magnetic local time displacement along a geomagnetic parallel; in this situation, the phase difference between geomagnetic fluctuations indicates the direction of their azimuthal propagation. Coherent fluctuations have been found to occur preferably when both stations are on the same side (dawnward or duskward) with respect to the polar cusp; moreover, around local magnetic midnight, they occur essentially during open magnetospheric conditions. The phase difference for coherent fluctuations indicates a propagation direction away from local geomagnetic noon and midnight. Also the analysis of three individual pulsation events, occurring at different times during the day, is shown; they are characterized at the two stations by simultaneous, coherent fluctuations, whose phase difference finds correspondence with the statistical behaviour.     

Study of mid-latitude ionospheric convection during quiet and disturbed periods using the SuperDARN Hokkaido radar

Westward ionospheric convective flows around midnight are frequently observed at mid-latitudes. They can be generated by so-called disturbance dynamo mechanisms working mainly in the mid-latitudes. To understand the influence of disturbance dynamo effects in the mid-latitudes, we studied the latitudinal distribution of westward flows in association with several kinds of geomagnetic disturbances using the SuperDARN Hokkaido radar. This radar creates high temporal resolution (1 s to 2 min), two-dimensional observations measuring the line-of-sight velocities of ionospheric plasma irregularities, which can be regarded as line-of-sight velocities of ionospheric convection in the mid-latitude region from 40° to 50°. This region could not be monitored using preexisting SuperDARN radars. In this study, we used ionospheric echo data obtained by the SuperDARN Hokkaido radar over 5 years (December 2006 to November 2011). We identified westward flows around midnight at about 40° to 55° geomagnetic latitude. Additionally, the data showed that the westward flow around midnight intensified under high geomagnetic activity (high Kp). This suggests that the disturbance dynamo could affect the mid-latitude ionospheric convection. We performed Superposed Epoch Analysis (SEA) to study the influences from the geomagnetic disturbances on mid-latitude ionospheric convection. We found no obvious influence during major storms (minimum Dst below −60 nT). SEA was also used to study the temporal and latitudinal dependence on the influences from substorms. From analysis of 36 events of AL-defined substorms, we saw that the influence of substorms lasted from 5 to 20 h after the onset between 44° and 53° geomagnetic latitude. The westward flow at mid-latitude grew to a maximum at 12 h after the geomagnetic substorm onset. This is consistent with the results of past numerical simulation studies of the disturbance dynamo effects.