中低纬电离层暴的统计分析与经验建模

As a key issue of ionospheric weather study, systemic studies on ionospheric storms can not only further improve our understanding of the response of the ionosphere to solar and geomagnetic disturbances, but also help us to reveal the chemical, dynamic and electro-dynamic mechanisms during storms. Empirical modelling for regional ionospheric storm is also very useful, because it can provide us with tools and references for the forecasting and further practical application of ionospheric activity. In this thesis, we focus on describing and forecasting of ionospheric storms at middle and low latitudes. The main points of my investigations are listed as follows. (1) By using magnetic storms during the period over 50 years, the dependence of the type, onset time and time delay of the ionospheric storms on magnetic latitude, season and local time at middle and low latitudes in the East-Asian sector are studied. The results show that the occurrences of the types of ionospheric disturbances differ in latitude and season. The onset of the ionospheric storms depends on local time. At middle latitudes, most negative phase onsets are within the local time interval from night to early morning, and they rarely occurred in the local noon and afternoon sectors. At low latitudes, positive phases commence most frequently in the daytime sector as well as pre-midnight sector. The average time delays for both the positive and negative ionospheric storms increase with descending latitudes. The time delay has significant dependence on the local time of main phase onset (MPO). The time delay of positive response is shorter for daytime MPO and longer for night-time MPO, whereas the opposite applies for negative response. (2) Based on some previous researches, a primary empirical model for mid-latitude ionospheric disturbance is set up. By fitting to the observed data, we get a high accuracy with a mean RMSE of only 12-14% in summer and equinox. The model output has been compared with the output of STORM model, and the results show that, our model is much better than STORM in summer and a little better for some mid-latitude stations at equinox. Especially, for the type of two-step geomagnetic storm, our model can present twice descending of foF2 very well. In addition, our model can forecast positive ionospheric storms.

一种TEC地图构建方法与全球电子总含量的气候学研究

Global positioning system (GPS) can not only provide precise service for navigation and timing, but also be used to investigate the ionospheric variation. From the GPS observations, we can obtain total electron content (TEC), so-called GPS TEC, which is used to characterize the ionospheric structure. This thesis mainly concerns about GPS TEC data processing and ionospheric climatological analysis as follows. Firstly, develop an algorithm for high-resolution global ionospheric TEC mapping. According to current algorithms in global TEC mapping, we propose a practical way to calibrate the original GPS TEC with the existing GIM results. We also finish global/local TEC mapping by model fitting with the processed GPS TEC data; in practice, we apply it into the local TEC mapping in Southeast of China and obtain some initial results. Next, suggest a new method to calculate equivalent ionospheric global electron content (GEC). We calculate such an equivalent GEC with the TEC data along the geographic longitude 120°E. With the climatological analysis, we can see that GEC climatological variation is mainly composed of three factors: solar cycle, annual and semiannual variations. Solar cycle variation is dominant among them, which indicates the most prominent influence; both annual and semiannual variations play a secondary role and are modulated by solar activity. We construct an empirical GEC model driven by solar activity and seasonal factors on the basis of partial correlation analysis. Generally speaking, our researches not only show that GPS is advantageous in now-casting ionospheric TEC as an important observation, but also show that GEC may become a new index to describe the solar influence on the global ionosphere since the great correlation between GEC and solar activity factor indicates the close relationship between the ionosphere and solar activity.

太阳活动变化对中低纬电离层的影响

With the variations of solar activity, solar EUV and X-ray radiations change over different timescales (e.g., from solar cycle variation to solar flare burst). Since solar EUV and X-ray radiations are the primary energy sources for the ionosphere, theirs variations undoubtedly produce significant and complicated effects on the ionosphere. So the variations of solar activity significantly affect the ionosphere. It is essential for both ionospheric theory and applications to study solar activity effects on the ionosphere. The study about solar activity variations of the ionosphere is an important part of the ionospheric climatology. It can enhance the understanding for the basic processes in the ionosphere, ionospheric structure and its change, ionosphere/thermosphere coupling, and so on. As for applications, people need sufficient knowledges about solar activity variations of the ionosphere in order to improve ionospheric models so that more accurate forecast for the ionospheric environments can be made. Presently, the whole image about the modalities of ionospheric solar activity variations is still unknown, and related mechanisms still cannot be well understood. This paper is about the effects of the 11-year change in solar activity to the low- and mid-latitude ionosphere. We use multi-type ionospheric observations and model to investigate solar activity effects on the electron density and ionospheric spatial structure, and we focus on discussing some related mechanisms. The main works are as follows: Firstly, solar activity variations of ionospheric peak electron density (NmF2) around 1400 LT were investigated using ionosonde observations in the 120°E sector. The result shows that the variation trend of NmF2 with F107 depends on latitudes and seasons. There is obvious saturation trend in low latitudes in all seasons; while in middle latitudes, NmF2 increases linearly with F107 in winter but saturates with F107 at higher solar activity levels in the other seasons. We calculated the photochemical equilibrium electron density to discuss the effects induced by the changes of neutral atmosphere and dynamics processes on the solar activity variations of NmF2. We found that: (1) Seasonal variation of neutral atmosphere plays an important role in the seasonal difference of the solar activity variations of NmF2 in middle latitudes. (2) Less [O]/[N2] and higher neutral temperature are important for the saturation effect in summer, and the increase of vibrational excited N2 is also important for the saturation effect. (3) Dynamics processes can significantly weaken the increase of NmF2 when solar activity enhances, which is also a necessary factor for the saturation effect. Secondly, solar activity variations of nighttime NmF2 were investigated using ionosonde observations in the 120°E sector. The result shows that the variation trends of NmF2 with F107 in nighttime are different from that in daytime in some cases, and the nighttime variation trends depend on seasons. There is linear increase trend in equinox nighttime, and saturation trend in summer nighttime, while the increase rate of NmF2 with F107 increases when solar activity enhances in winter nighttime (we term it with “amplification trend”). We discussed the possible mechanisms which affect the solar activity variations of nighttime NmF2. The primary conclusions are as follows: (1) In the equatorial ionization anomaly (EIA) crest region, the plasma influx induced by the pre-reversal enhancement (PRE) results in the change of the variation trend between NmF2 and F107 from “saturation” to “linear” after sunset in equinoxes and winter; while the recombination process at the F2-peak is the primary factor that affects the variation trend of NmF2 with F107 in middle latitudes. (2) The recombination coefficient at the F2-peak height reaches its maximum at moderate solar activity level in winter nighttime, which induces NmF2 attenuates more quickly at moderate solar activity level. This is the main reason for the amplification trend. (3) The change of the recombination process at the F2-peak with solar activity depends on the increases of neutral parameters (temperature, density et al.) and the F2-peak height (hmF2). The seasonal differences in the changes of neutral atmosphere and hmF2 with solar activity are the primary reasons for the seasonal difference in the variation trend of nighttime NmF2 with F107. Finally, we investigated the solar activity dependence of the topside ionosphere in low latitudes using ROCSAT-1 satellite (at 600 km altitude) observations. The primary results and conclusions are as follows: (1) Latitudinal distribution of the plasma density is local time, seasonal, and solar activity dependent. In daytime, there is a plasma density peak at the dip equator. The peak is obviously enhanced at high solar activity level, and the strength of the peak strongly depends on seasons. While at sunset, two profound plasma density peaks (double-peak structure) are found in solar maximum equinox months. (2) Local time dependence of the latitudinal distribution is due to the local time variation of the equatorial dynamics processes. Double-peak structure is attributed to the fountain effect induced by strong PRE. Daytime peak enhances with solar activity since the plasma density increases with solar activity more strongly at the dip equator due to the equatorial vertical drift, and its seasonal dependence is mainly due to the seasonal variations of neutral density and the equatorial vertical drift. In the sunset sector, seasonal and solar activity dependences of the latitudinal distribution are related to the seasonal and solar activity variations of PRE. (3) The variation trend of the plasma density with solar activity shows local time, seasonal, and latitudinal differences. That is different from the changeless amplification trend at the DMSP altitude (840 km). Profound saturation effect is found in the dip equator region at equinox sunset. This saturation effect in the topside ionosphere is realated to the increase of PRE with solar activity. Solar activity variation trend of the topside plasma density was discussed quantitatively by Chapman-α function. The result shows that the effect induced by the change of the scale height is dominant at high altitudes; while the variation trend of ROCSAT-1 plasma density with solar activity is suggested to be related to the changes of the peak height, the scale height, and the peak electron density with solar activity.

太阳耀斑和日食期间电离层行为的统计和模拟研究

The ionosphere is the ionized component of the Earth's upper atmosphere. Solar EUV radiation is the source of ionospheric ionization. Thus the ionosphere is affected strongly by the variations in solar radiation. Solar flares and solar eclipses can induce remarkable short time changes in solar radiation: the solar radiation would increase suddenly during solar flares and decrease significantly during solar eclipses. Solar flare and eclipse events not only affect directly the photochemical processes, but also affect the dynamic processes, and even affect the neutral atmosphere, which is strongly coupled with the ionosphere. The study on the ionospheric response to solar flares and eclipses can advance our knowledge on the ionosphere and its photochemical and dynamic processes and help us to evaluate the ionospheric parameters (such as ion loss coefficients). In addition, the study on the ionospheric responses to solar flares and eclipses is an important part of the ionospheric space weather, which can provide guides for space weather monitoring. This thesis devotes to the study on the ionospheric responses to solar flares and solar eclipses. I have developed two models to simulate the variations of solar EUV radiation during solar flares and solar eclipses, and involved in developing a 2D mid- and low-latitude ionospheric model. On the basis of some observed data and the ionospheric model, I study the temporal and spatial variations of the ionosphere during solar flares and eclipses, and investigate the influences of solar activity, solar zenith angle, neutral gas density, and magnetic dip angle on the ionospheric responses to solar flares and solar eclipses. The main points of my works and results are summarized as follows. 1. The ionospheric response to the X17.2 solar flare on October 28, 2003 was modeled via using a one-dimension theoretical ionospheric model. The simulated variation of TEC is in accordance with the observations, though there are some differences in the amplitude of the variation. Then I carried out a series of simulations to explore the local time and seasonal dependences of the ionospheric responses to solar flares. These calculations show that the ionospheric responses are largely related with the solar zenith angle (SZA). During the daytime (small SZA), most of the increases in electron density occur at altitudes below 300 km with a peak at around 115 km; whereas around sunrise and sunset (SZA>90°), the strongest ionospheric responses occur at much higher altitudes. The TEC increases slower at sunrise than at sunset, which is caused by the difference in the evolution of SZA at sunrise and sunset: SZA decreases with time at sunrise and increase with time at sunset. The ionospheric response is largest in summer and smallest in winter, which is also related to the seasonal difference of SZA. 2. Based on the observations from the ionosondes in Europe and the ionospheric model, I investigated the differences of the ionosphere responses to solar eclipses between the E-layer and F1-layer. Both the observation and simulation show that the decrease in foF1 due to the solar eclipses is larger than that in foE. This effect is due to that the F1 region locates at the transition height between the atomic ion layer and the molecular ion layer. With the revised model of solar radiation during solar flares, our model calculates the radiations from both the inside and outside of photosphere. Large discrepancy can be found between the observations and the calculations with an unrevised model, while the calculations with the revised model consist with the observations. 3. I also explore the effects of the F2-layer height, local time, solar cycle, and magnetic dip angle on the ionospheric responses to solar eclipses via using an ionospheric model and study on the solar zenith angle and the dip dependences by analyzing the data derived from 23 ionosonde stations during seven eclipse events. Both the measured and simulated results show that these factors have significant effect on the ionospheric response. The larger F2-layer height causes the smaller decrease in foF2, which is because that the electron density response decreases with height. The larger dip results in the smaller eclipse effect on the F2 layer, because the larger dip would cause the more diffusion from the top ionosphere which can make up for the plasma loss. The foF2 response is largest at midday and decreases with the increasing SZA. The foF2 response is larger at high solar activity than at low solar activity. The simulated results show that the local time and solar activity discrepancy of the eclipse effect mainly attribute to the difference of the background neutral gas density. 4. I carried out a statistical study on the latitudinal dependence of the ionospheric response to solar eclipses and modeled this latitudinal dependence by the ionospheric model. Both the observations and simulations show that the foF2 and TEC responses have the same latitudinal dependence: the eclipse effects on foF2 and TEC are smaller at low latitudes than at middle latitudes; at the middle latitudes (>40°), the eclipse effect decreases with increasing latitude. In addition, the simulated results show the change in electron temperature at the heights of above 300 km of low latitudes is much smaller than that at the same heights of middle latitudes. This is due to the smaller decrease in photoelectron production rate at its conjugate low heights. 5. By analyzing the observed data during the October 3, 2005 solar eclipse, I find some significant disturbances in the conjugate region of the eclipse region, including a decrease in Te, an increase in foF2 and TEC, and an uprising in hmF2. I also simulated the ionosphere behavior during this eclipse using a mid-low latitude ionospheric model. The simulations reproduce the measured ionospheric disturbances mentioned above in the conjugated hemisphere. The simulations show that the great loss of arriving photoelectron heat from the eclipse region is the principal driving source for the disturbances in the conjugate hemisphere.

电离层胶片频高图数字化转换与分析研究

The ionogram acquired with the ionospheric vertical sounding method is the oldest data in the history of ionospheric research. Using of modern microelectronics and computer technology to digitalize, analyse and preserve the huge amount of historical film ionogram has become more and more important and urgent. This paper introduced the progress of the film ionogram digitalization by using digital image processing technologies to correct and repair film ionogram and convert them in an exchangeable format. An analysis and conversion software, basing on this method, has been developed for the film ionogram analysis, and then it introduces the application of this software by combining the SAO Explorer program for Wuhan film ionogram and pseudo-color ionogram in Yamagawa in Japan. It shows that our method is reliable，and the developed software is used friendly and provides a positive solution in digitalization and analysis of huge amount of historical film ionogram. Firstly, we briefly introduce the film ionogram and the process of its digitalization. By observing a amount of film ionogram, we obtain some common characteristics of the digitalized film ionogram following as: (1) the image rotation are caused by scanning; (2) the vertical axis of a large number of film ionogram exist more or less tilt and bending ; (3) coordinates of the film ionogram appear the non-uniformity phenomena result from the instability of driving motor rotation and the error of altitudinal cursor orientation. Moreover, based on the characteristics of the film ionogram and the SAO Explorer software which is widely used for the digital ionogram analysis in the world, a new method has been developed for film ionogram procession. The method contains the image geometric correction and film ionogram format conversion. The image geometric correction includes such as image rotation correction, vertical correction and coordinates scale correction. After geometric correction, the BMP file format images will be converted to the SBF file format images. Then, we also discuss the data format converting methods, which include two methods of the image data mapping basing on the normalization and logarithm, and the method of the preprocessing of the noise filtering and the threshold setting. Combining with SAO Explorer software, we successfully obtain ionospheric parameters and electron profile from the converted SBF file format digital ionograms. Based on the above method, we developed the software for the film ionogram to realize its correction analysis and conversion of the image format, and then give a introduction for its function and operation. Subsequently, the software are applied into the Wuhan film ionogram which separately observed in the high solar activity year and the low in 1980s last century. The results reveal the converted SBF digital ionogram almost preserve the all echo information of the film ionogram. Furthermore, we expressly discuss the application to the Wuhan film ionogram in 1958 in order to validate the applicability and credibility of the software. And it is showed that the important information of the film ionogram are maintained into the SBF digital ionogram. It is represented that there is credibility for conversion of the software when it applied in the older film ionogram. In sum, this software could apply to the digitalization and analysis of huge amount of historical film ionogram. Last, we extended the function of the software by bring some new conversion method and used it to apply to the pseudo-color ionogram of yamagawa in Japan. The results show that the converted ionogram information basically maintain the importantly ionogram information and the error of scaling of converted SBF file format image is almost acceptable, though there is no preprocessing for the original ionogram. Hence, we could extend the applicable range of the software and apply it to all kinds of simulative ionogram imaging by improving the method and software.

AE指数的可靠性研究

Auroral electrojet index (AE) are usually used to quantitatively describe the activity of the geomagnetic field in the polar region. AE is a means to identify the level of a substorm as well. The auroral electrojet indices (AU, AL, and AE) have served well for more than four decades as measures of magnetospheric substorm activity. However, as substorm studies have progressed considerably during the past several years, the accuracy of the present electrojet indices have become an important issue. Thus it is fortunate for us to reexamine and evaluate the accuracy of the present electrojet indices and improve them if necessary. For a better use of the present indices and for future improvement we examine the limitations of the auroral electrojet indices as an accurate quantitative measure of the auroral electrojets and magnetospheric substorms. Some of the limitations arise from the data availability and also from the present simplified scheme in deriving them, but some of them originate in the definition themselves. In the present paper, we analyze the characteristics of aurora, ionospheric current system and the AE index. It is noted that for the pseudo-breakup events the variations of AE ( > 500 nT) clearly show the sudden increase and slow decrease phases. However, aurora does not expand to poleward or equatorward, and the ionospheric currents presents the features of the magnetic convection. We mainly focus on investigating what time the AE could be used to identify the subtorm.

电离层垂直TEC映射函数的实验观测与统计特性

When used in the determining the total electron content (TEC), which may be the most important ionospheric parameter, the worldwide GPS observation brings a revolutionary change in the ionospheric science. There are three steps in the data processing to retrieve GPS TEC: (1) to estimate slant TEC from the measurements of GPS signals; (2) to map the slant TEC into vertical; and (3) to interpolate the vertical TEC into grid points. In this scientific dissertation we focus our attention on the second step, the mapping theory and method to convert slant TEC into vertical. This is conventionally done by multiplying on the slant TEC a mapping function which is usually determined by certain models of electron density profile. Study of the vertical TEC mapping function is of significance in GPS TEC measurement. This paper first reviews briefly the three steps in GPS TEC mapping process. Then we compare the vertical TEC mapping function which were respectively calculated from the electron density profiles of the ionospheric model and retrieved from the observation of worldwide GPS TEC. We also perform the statistical analysis on the observational mapping functions. The main works and results are as follows: 1. We calculated the vertical TEC mapping functions for both SLM and Chapman models, and discussed the modulation of the ionosphere height to the mapping functions. We use two simple models, single layer model (SLM) and Chapman models, of the ionospheric electron density profiles to calculate the vertical TEC mapping function. In the case of the SLM, we discuss the control of the ionospheric altitude, i.e., the layer height hipp, to the mapping function. We find that the mapping function decreases rapidly as hipp increases. For the Chapman model we study also the control mapping function by both ionospheric altitude indicated by the peak electron density height hmF2, and the scale height, H, which present the thickness of the ionosphere. It is also found that the mapping function decreases rapidly as hmF2 increases. and it also decreases as H increases. 2. Then we estimate the mapping functions from the GPS observations and compare them with those calculated from the electron density models. We first, proposed a new method to estimate the mapping functions from GPS TEC data. This method is then used to retrieve the observational mapping function from both the slant TEC (TECS) provided by International GPS Service (IGS）and vertical TEC provide by JPL Global Ionospheric Maps (GIMs). Then we compare the observational mapping function with those calculated from the electron density models, SLM and Chapman. We find that the values of the observational mapping functions are much smaller than that from the model mapping functions, when the zenith angle is large enough. We attribute this to the effect of the plasmasphere which is above about 1000 km. 3. We statistically analyze the observational mapping functions and reveal their climatological changes. Observational mapping functions during 1999-2007 are used in our statistics. The main results are as follows. (1) The observational mapping functions decrease obviously with the decrement of the solar activity which is represented by the F10.7 index; (2) In annual variations of the observational mapping functions, the semiannual component is found at low-latitudes, and the remarkable seasonal variations at mid- and high-latitudes. (3) The diurnal variation of the observational mapping functions is that they are large in daytime and small at night, they become extremely small in the early morning before sunrise. (4) The observational mapping functions change with latitudes that they are smaller at lower latitudes and larger at higher. All of the above variations of the observational mapping functions are explained by the existence of the plasmasphere, which changes more slowly with time and more rapidly with latitude than the ionosphere does . In summary, our study on the vertical TEC mapping function imply that the ionosphere height has a modulative effect on the mapping function. We first propose the concept of the 'observational mapping functions' , and provide a new method to calculate them. This is important in improving the TEC mapping. It may also possible to retrieving the plasmaspheric information from GPS observations.

Millstone Hill电离层F层等效风场的气候学研究

Neutral winds and electric fields in the ionospheric F layer play important roles in the variations of the ionosphere, and also affect the thermospheric circulation via the close coupling between the ionosphere and the thermosphere. By now, the neutral winds and electric drifts are generally observed with ground-based Fabry-Perot interferometers (FPI) and incoherent scatter radars (ISR), rockets, and satellite-borne instrument. Based on the servo theory, the ionospheric equivalent winds, which include the information of both the neutral winds and electric fields, can be derived from these characteristic parameters observed by ionosondes. This indirect derivation has potential values in climatological researches and space weather forecast. With the data set of the incoherent scatter radar observations at Millstone Hill, USA, from 1976 to 2006, we statistically analyzed the climatological variations of the vertical component of the equivalent winds (VEWs) over Millstone Hill, which are derived from the ionospheric key parameters (the peak electron number density and peak height of the F2 layer, NmF2 and hmF2) on the basis of the servo theory, Liu's method, and measurements from the ion line-of-sight velocity as well. The main results of this analysis are summarized as follows: (1) The values of VEWs over Millstone Hill during nighttime are stronger than in the daytime, and the upward drift dominates most of the day. In 1993, Hagan found that the component of the neutral winds in the magnetic meridion in daytime is weaker than during nighttime under both solar maximum and minimum conditions; he also found that the equatorward winds dominate most of the day. Both results suggest that the thermosphere in Millstone Hill is modulated by the aurorally driven high-latitude circulation cell; that is, during geomagnetic quiet periods, the average auroral activity is strong enough to drive thermospheric circulation equatorward for most of the day at Millstone Hill. Moreover, since ion drag is the strongest during daytime when F region densities are enhanced by photoionization, the wind speeds are smaller during the daytime than in the nighttime. (2) There is equinoctial symmetry in VEWs at Millstone Hill. The amplitudes and phases of VEWs in spring are quite similar to those in autumn. In contrast, the nighttime upward drift in winter is weaker than in summer and the difference becomes more significant with increasing solar activity. This solstice asymmetry indicates that, the aurorally driven circulation in the northern hemisphere at Millstone Hill latitude is weaker in winter due to arctic darkness, because the subsolar point is in the southern hemisphere. (3) The comparison of the VEWs derived from three methods, i.e., the servo theory, Liu's method, and the ISR ion line-of-sight velocity measurements, indicates that the amplitudes and main phase tendencies of these VEWs accord well with each other during nighttime hours. However, the case in the daytime is relatively worse. This daytime discrepancy can be explained in terms of the effects of photochemical processes and the choices of the servo constants. A larger servo constant gives a stronger plasma drift in daytime. Therefore, this study tells how important to choose a suitable constant for deriving VEWs at Millstone Hill.

磁层－电离层电流的暴时变化以及主磁场对空间环境的影响

Two problems are studied in this thesis, the relationship of the magneto-spheric - ionospheric current systems during storms, and the effects of the main field to the space environment. The thesis includes three parts. 1. Magnetic disturbances caused by magnetospheric - ionospheric current systems Transient variations of the geomagnetic field at middle-low latitudes are mainly caused by the ionospheric dynamo current (IDC), the symmetric ring current (SRC), the partial ring current-region II field-aligned current-ionospheric current system (PRFI), and the region I field-aligned current-ionospheric current system (FACI). The storm on May 1 ～ 6, 1998 is analyzed. Firstly, the S_q-field caused by IDC current is removed by using the modified Hibberd's method in which the effect of SRC is considered. The neglect of SRC-field can give as much as 40% error in S_q-field evaluation. Secondly, the disturbance fields at the middle and low latitudes are separated according to their origins. As a result, the disturbance caused by FACI-current is an important part of the asymmetrical depression of H-component in middle and low latitudes during storms. The results show that the relative intensity of the Sq-field increases in the main phase of the storm and decreases in the recovery phase. The latitudinal gradient of the Sq-field is positive during the whole storm. The storm of May 1 ～ 6, 1998 contains two events. In the first event on May 2, the SRC-field is similar to Dst index. But in the second event on May 4 ～ 5, the SRC-field delays to Dst index, and the SRC-field depresses while the PRFI- and FACI-fields recovery. 2. Analysis of S_q~p variation in CGM coordinates In order to study the conjugation of geomagnetic variations between northern and southern hemispheres, we use the corrected geomagnetic coordinates (CGM) instead of the geomagnetic coordinates (GM) to analyze the S_q~P equivalent current system. The CGM coordinates are built up by International Geomagnetic Reference Field (IGRF) model. The S_q~p variations and equivalent current systems in the northern and southern polar regions are more symmetrical in CGM coordinates than in GM co-ordinates. This fact implies that the current distributions in polar regions are governed by the configuration of the geomagnetic field lines. As the elaborate structure of S_q~p current system in quiet time is obtained, we summarize the seasonal variation of the electrojet in quiet time. 3. The magnetospheric configuration of non-parallel-dipole model The magnetospheric configurations are calculated for two possible geomag-netic field models during the geomagnetic field reversals. These models are the dipole field with the axis to the sun and the quadrupole field model. We use the finite element method to solve the magnetic equation, and use the surface evolution method to solve the equilibrium equation. The results show that the main field greatly affects the space environment.