Generalized-positioning for mixed-frequency of mixed-GNSS and its preliminary applications

Modernized GPS and GLONASS, together with new GNSS systems, BeiDou and Galileo, offer code and phase ranging signals in three or more carriers. Traditionally, dual-frequency code and/or phase GPS measurements are linearly combined to eliminate effects of ionosphere delays in various positioning and analysis. This typical treatment method has imitations in processing signals at three or more frequencies from more than one system and can be hardly adapted itself to cope with the booming of various receivers with a broad variety of singles. In this contribution, a generalized-positioning model that the navigation system independent and the carrier number unrelated is promoted, which is suitable for both single- and multi-sites data processing. For the synchronization of different signals, uncalibrated signal delays (USD) are more generally defined to compensate the signal specific offsets in code and phase signals respectively. In addition, the ionospheric delays are included in the parameterization with an elaborate consideration. Based on the analysis of the algebraic structures, this generalized-positioning model is further refined with a set of proper constrains to regularize the datum deficiency of the observation equation system. With this new model, uncalibrated signal delays (USD) and ionospheric delays are derived for both GPS and BeiDou with a large dada set. Numerical results demonstrate that, with a limited number of stations, the uncalibrated code delays (UCD) are determinate to a precision of about 0.1 ns for GPS and 0.4 ns for BeiDou signals, while the uncalibrated phase delays (UPD) for L1 and L2 are generated with 37 stations evenly distributed in China for GPS with a consistency of about 0.3 cycle. Extra experiments concerning the performance of this novel model in point positioning with mixed-frequencies of mixed-constellations is analyzed, in which the USD parameters are fixed with our generated values. The results are evaluated in terms of both positioning accuracy and convergence time.

A reference station-based GNSS computing mode to support unified precise point positioning and real-time kinematic services

Currently, the GNSS computing modes are of two classes: network-based data processing and user receiver-based processing. A GNSS reference receiver station essentially contributes raw measurement data in either the RINEX file format or as real-time data streams in the RTCM format. Very little computation is carried out by the reference station. The existing network-based processing modes, regardless of whether they are executed in real-time or post-processed modes, are centralised or sequential. This paper describes a distributed GNSS computing framework that incorporates three GNSS modes: reference station-based, user receiver-based and network-based data processing. Raw data streams from each GNSS reference receiver station are processed in a distributed manner, i.e., either at the station itself or at a hosting data server/processor, to generate station-based solutions, or reference receiver-specific parameters. These may include precise receiver clock, zenith tropospheric delay, differential code biases, ambiguity parameters, ionospheric delays, as well as line-of-sight information such as azimuth and elevation angles. Covariance information for estimated parameters may also be optionally provided. In such a mode the nearby precise point positioning (PPP) or real-time kinematic (RTK) users can directly use the corrections from all or some of the stations for real-time precise positioning via a data server. At the user receiver, PPP and RTK techniques are unified under the same observation models, and the distinction is how the user receiver software deals with corrections from the reference station solutions and the ambiguity estimation in the observation equations. Numerical tests demonstrate good convergence behaviour for differential code bias and ambiguity estimates derived individually with single reference stations. With station-based solutions from three reference stations within distances of 22–103 km the user receiver positioning results, with various schemes, show an accuracy improvement of the proposed station-augmented PPP and ambiguity-fixed PPP solutions with respect to the standard float PPP solutions without station augmentation and ambiguity resolutions. Overall, the proposed reference station-based GNSS computing mode can support PPP and RTK positioning services as a simpler alternative to the existing network-based RTK or regionally augmented PPP systems.

GNSS ambiguity resolution with controllable failure rate for long baseline network RTK

Many large-scale GNSS CORS networks have been deployed around the world to support various commercial and scientific applications. To make use of these networks for real-time kinematic positioning services, one of the major challenges is the ambiguity resolution (AR) over long inter-station baselines in the presence of considerable atmosphere biases. Usually, the widelane ambiguities are fixed first, followed by the procedure of determination of the narrowlane ambiguity integers based on the ionosphere-free model in which the widelane integers are introduced as known quantities. This paper seeks to improve the AR performance over long baseline through efficient procedures for improved float solutions and ambiguity fixing. The contribution is threefold: (1) instead of using the ionosphere-free measurements, the absolute and/or relative ionospheric constraints are introduced in the ionosphere-constrained model to enhance the model strength, thus resulting in the better float solutions; (2) the realistic widelane ambiguity precision is estimated by capturing the multipath effects due to the observation complexity, leading to improvement of reliability of widelane AR; (3) for the narrowlane AR, the partial AR for a subset of ambiguities selected according to the successively increased elevation is applied. For fixing the scalar ambiguity, an error probability controllable rounding method is proposed. The established ionosphere-constrained model can be efficiently solved based on the sequential Kalman filter. It can be either reduced to some special models simply by adjusting the variances of ionospheric constraints, or extended with more parameters and constraints. The presented methodology is tested over seven baselines of around 100 km from USA CORS network. The results show that the new widelane AR scheme can obtain the 99.4 % successful fixing rate with 0.6 % failure rate; while the new rounding method of narrowlane AR can obtain the fix rate of 89 % with failure rate of 0.8 %. In summary, the AR reliability can be efficiently improved with rigorous controllable probability of incorrectly fixed ambiguities.

Multi-GNSS precise point positioning with raw single-frequency and dual-frequency measurement models

The emergence of multiple satellite navigation systems, including BDS, Galileo, modernized GPS, and GLONASS, brings great opportunities and challenges for precise point positioning (PPP). We study the contributions of various GNSS combinations to PPP performance based on undifferenced or raw observations, in which the signal delays and ionospheric delays must be considered. A priori ionospheric knowledge, such as regional or global corrections, strengthens the estimation of ionospheric delay parameters. The undifferenced models are generally more suitable for single-, dual-, or multi-frequency data processing for single or combined GNSS constellations. Another advantage over ionospheric-free PPP models is that undifferenced models avoid noise amplification by linear combinations. Extensive performance evaluations are conducted with multi-GNSS data sets collected from 105 MGEX stations in July 2014. Dual-frequency PPP results from each single constellation show that the convergence time of undifferenced PPP solution is usually shorter than that of ionospheric-free PPP solutions, while the positioning accuracy of undifferenced PPP shows more improvement for the GLONASS system. In addition, the GLONASS undifferenced PPP results demonstrate performance advantages in high latitude areas, while this impact is less obvious in the GPS/GLONASS combined configuration. The results have also indicated that the BDS GEO satellites have negative impacts on the undifferenced PPP performance given the current “poor” orbit and clock knowledge of GEO satellites. More generally, the multi-GNSS undifferenced PPP results have shown improvements in the convergence time by more than 60 % in both the single- and dual-frequency PPP results, while the positioning accuracy after convergence indicates no significant improvements for the dual-frequency PPP solutions, but an improvement of about 25 % on average for the single-frequency PPP solutions.

Wide area kinematic positioning with BeiDou triple frequency signals

A key challenge of wide area kinematic positioning is to overcome the effects of the varying hardware biases in code signals of the BeiDou system. Based on three geometryfree/ionosphere-free combinations, the elevation-dependent code biases are modelled for all BeiDou satellites. Results from the data sets of 30-day for 5 baselines of 533 to 2545 km demonstrate that the wide-lane (WL) integer-fixing success rates of 98% to 100% can be achieved within 25 min. Under the condition of HDOP of less than 2, the overall RMS statistics show that ionospheric-free WL single-epoch solutions achieve 24 to 50 cm in the horizontal direction. Smoothing processing over the moving window of 20 min reduces the RMS values by a factor of about 2. Considering distance-independent nature, the above results show the potential that reliable and high precision positioning services could be provided in a wide area based on a sparsely distributed ground network.

Observations of the solar wind-magnetosphere-ionosphere coupling

In this thesis, the solar wind-magnetosphere-ionosphere coupling is studied observationally, with the main focus on the ionospheric currents in the auroral region. The thesis consists of five research articles and an introductory part that summarises the most important results reached in the articles and places them in a wider context within the field of space physics. Ionospheric measurements are provided by the International Monitor for Auroral Geomagnetic Effects (IMAGE) magnetometer network, by the low-orbit CHAllenging Minisatellite Payload (CHAMP) satellite, by the European Incoherent SCATter (EISCAT) radar, and by the Imager for Magnetopause-to-Aurora Global Exploration (IMAGE) satellite. Magnetospheric observations, on the other hand, are acquired from the four spacecraft of the Cluster mission, and solar wind observations from the Advanced Composition Explorer (ACE) and Wind spacecraft. Within the framework of this study, a new method for determining the ionospheric currents from low-orbit satellite-based magnetic field data is developed. In contrast to previous techniques, all three current density components can be determined on a matching spatial scale, and the validity of the necessary one-dimensionality approximation, and thus, the quality of the results, can be estimated directly from the data. The new method is applied to derive an empirical model for estimating the Hall-to-Pedersen conductance ratio from ground-based magnetic field data, and to investigate the statistical dependence of the large-scale ionospheric currents on solar wind and geomagnetic parameters. Equations describing the amount of field-aligned current in the auroral region, as well as the location of the auroral electrojets, as a function of these parameters are derived. Moreover, the mesoscale (10-1000 km) ionospheric equivalent currents related to two magnetotail plasma sheet phenomena, bursty bulk flows and flux ropes, are studied. Based on the analysis of 22 events, the typical equivalent current pattern related to bursty bulk flows is established. For the flux ropes, on the other hand, only two conjugate events are found. As the equivalent current patterns during these two events are not similar, it is suggested that the ionospheric signatures of a flux rope depend on the orientation and the length of the structure, but analysis of additional events is required to determine the possible ionospheric connection of flux ropes.

Ionosphere-atmosphere interaction during solar proton events

Among the most striking natural phenomena affecting ozone are solar proton events (SPE), during which high-energy protons precipitate into the middle atmosphere in the polar regions. Ionisation caused by the protons results in changes in the lower ionosphere, and in production of neutral odd nitrogen and odd hydrogen species which then destroy ozone in well-known catalytic chemical reaction chains. Large SPEs are able to decrease the ozone concentration of upper stratosphere and mesosphere, but are not expected to significantly affect the ozone layer at 15--30~km altitude. In this work we have used the Sodankylä Ion and Neutral Chemistry Model (SIC) in studies of the short-term effects caused by SPEs. The model results were found to be in a good agreement with ionospheric observations from incoherent scatter radars, riometers, and VLF radio receivers as well as with measurements from the GOMOS/Envisat satellite instrument. For the first time, GOMOS was able to observe the SPE effects on odd nitrogen and ozone in the winter polar region. Ozone observations from GOMOS were validated against those from MIPAS/Envisat instrument, and a good agreement was found throughout the middle atmosphere. For the case of the SPE of October/November 2003, long-term ozone depletion was observed in the upper stratosphere. The depletion was further enhanced by the descent of odd nitrogen from the mesosphere inside the polar vortex, until the recovery occurred in late December. During the event, substantial diurnal variation of ozone depletion was seen in the mesosphere, caused mainly by the the strong diurnal cycle of the odd hydrogen species. In the lower ionosphere, SPEs increase the electron density which is very low in normal conditions. Therefore, SPEs make radar observations easier. In the case of the SPE of October, 1989, we studied the sunset transition of negative charge from electrons to ions, a long-standing problem. The observed phenomenon, which is controlled by the amount of solar radiation, was successfully explained by considering twilight changes in both the rate of photodetachment of negative ions and concentrations of minor neutral species. Changes in the magnetic field of the Earth control the extent of SPE-affected area. For the SPE of November 2001, the results indicated that for low and middle levels of geomagnetic disturbance the estimated cosmic radio noise absorption levels based on a magnetic field model are in a good agreement with ionospheric observations. For high levels of disturbance, the model overestimates the stretching of the geomagnetic field and the geographical extent of SPE-affected area. This work shows the importance of ionosphere-atmosphere interaction for SPE studies. By using both ionospheric and atmospheric observations, we have been able to cover for the most part the whole chain of SPE-triggered processes, from proton-induced ionisation to depletion of ozone.

Low-frequency variability and the galactic distribution of electron-density fluctuation

A sample of 96 compact flat-spectrum extragalactic sources, spread evenly over all galactic latitudes, has been studied at 327 MHz for variability over a time interval of about 15 yr. The variability shows a dependence on galactic latitude being less both at low and high latitudes and peaking around absolute value of b approximately 15-degrees. The latitude dependence is surprisingly similar in both the galactic centre and anticentre directions. Assuming various single and multi-component distributions for the ionized, irregular interstellar plasma, we have tried to generate the observed dependence using a semi-qualitative treatment of refractive interstellar scintillations. We find that it is difficult to fit our data with any single or double component cylindrical distribution. Our data suggests that the observed variability could be influenced by the spiral structure of our Galaxy.

Some observational evidence of Alfven surface waves induced magnetic reconnection

The surface wave induced magnetic reconnection (SWIMR) model based on Alfven Resonance theory will be discussed briefly both for collisional and collisionless plasmas. It is shown that the spatial scales and time delays associated with Flux Transfer Events and Pulsed Ionospheric Flows, as observed by satellites and SuperDARN radars and the magnetic bubbles, observed at the high latitude boundary of the magnetopause, can be explained by the SWIMR model.

On the linkage of mesospheric planetary waves with those of the lower atmosphere and ionosphere: A case study from Indian low latitudes

In this paper we study the planetary-scale wave features using concurrent observations of mesospheric wind and temperature, ionospheric h'F, and tropospheric wind from Tirunelveli, Gadanki, and Kolhapur, all located in the Indian low latitudes, made during February 2009. Our investigations reveal that 3 to 5 day periodicity, characterized as ultrafast Kelvin (UFK) waves, was persistent throughout the atmosphere during this period. These waves show clear signatures of upward wave propagation from troposphere to the upper mesosphere, linking the ionosphere through a clear correlation between mesospheric winds and h'F variations. We also note that the amplitude of this wave decreased as we moved away from the equator. These results are the first of their kind from Indian sector, portraying the vertical as well as latitudinal characteristics of the 3 to 5 day UFK waves simultaneously from the troposphere to the ionosphere.

中低纬度变化磁场的暴时特征研究

Magnetic storm is a kind of severe disturbances in the whole solar-earth electromagnetic space. It has significant effects on communication, electric power, oil transport pipe and human activities in space. Therefore, magnetic storms are worth for applications systems, not only being a favorable issue for scientists. In this paper, the spatial and temporal distributions of the magnetic fields produced by the magnetosphere-ionosphere current systems during storms are studied. Four parts are included in this paper decomposion of different disturbances with different origins, topological structure of the ring current, the asymmetric characteristics of the ring current, and the statistic peculiarities of the day-to-day variability (DTD) of Sq. 1 The decomposition of magnetic disturbances at mid-low latitudes and its evolutions during storms Transient variations in the geomagnetic field recorded at mid-low latitudes mainly include the storm-time variation (Dst), solar quiet daily variation (Sq) and disturbance daily variation (SD). With the data of the geomagnetic meridian chain observatories in China, 25 storms during the period of 1997 to 1999 have been analyzed. According to the features of different variations, a method of “three-steps decomposition” is developed by using the method of Natural Orthogonal Components (NOC), Correlation Analysis and Fourier Analysis to separate those three components in turn. The results show that, the first eigenmode by the MNOC clearly describing the special distribution and temporal evolution of storm-time variation, in addition, Correlation Analysis and Fourier Analysis offer a useful method to extract the Sq and SD variations. The latitudinal shift of the Sq current focus seems to be the principal reason of the day-to-day variaitons in the daily range of Sq. The magnitude of SD reaches a maximum during the main phase, and then gradually decreases. 2 The topology structure of the ring current during storms Both the mechanism of the ring current and the geomagnetic data suggest that the central plane of the ring current is declining to the geomagnetic equator plane with a tilt angle δ. Using the H and Z component data at two stations in a meridian chain, we deduce a new parameter describing the invariable peculiarity of different storms. Then the δ angle is calculated by using the data from a meridian chain and tested with the ERC model. Finally the deduced tilt angles are used to modify Dst index. 3 The asymmetric characteristics of the ring current during storms The variations of the geomagnetic field at mid-low latitudes show a significant dawn-dusk asymmetry, resulting from the superposition of the fields from the symmetric ring current and the partial ring current. On the basis of the data from the 20°E, 30°E meridian chains and 30°N latitudinal chain, the dawn-dusk asymmetry is investigated by using three methods, namely, statistic analysis, ring current model calculation and typical event analysis. This characteristic implies the asymmetry of the spatial distribution of the ring current. In addition, during the main phase after the sudden commencement (SC), H field increases and reaches maximum around noontime, implying the effect of the Chapman-Ferraro current. 4 The statistic characteristics of the day-to-day variability and its mechanism The day-to-day variability of the geomagnetic Sq field is studied by using the magnetic data from a meridian chain of magnetometers along 120° E longitude. The method of NOC is applied to separate the Sq variation from complicated disturbances. The first eigenmode with the largest eigenvalue represents fairly well the Sq variation with a conspicuous day-to-day variability in the daily range. For the stations on the same north- or south-side of the Sq current system focus, the day-to-day variations show a positive correlation. In contrast, for the stations on the different sides of the Sq focus, they show a negative correlation, suggesting an important role of latitudinal shift of the Sq current system focus to the day-to-day variability of the Sq daily range. The Sq daily range is correlated with the magnetic indices Ap and Dst in a peculiar way: on some severe disturbed days, noticeably enhancements of the Sq are observed, implying increases of the ionospheric conductivities and/or tidal wind velocities; on other severe disturbed days, however, dramatically reduced Sq variations occur, suggesting dominant effects of the ‘disturbance dynamo’ process.

电离层行进式扰动的GPS台网监测研究

In this dissertation, we investigated two types of traveling ionospheric disturbances (TIDs)/gravity waves (GWs) triggered separately by auroral energy input during super geomagnetic storms and solar terminator (ST) under quiet geomagnetic conditions (kp<3+) using TEC measurements from the global network of GPS receivers. Research into the generation and propagation of TIDs/GWs during storms greatly enhance our understandings on the evolution processes of energy transportation from the high-latitude’s magnetosphere to the low-latitude ionosphere and the conjugated effect of TIDs propagation between the northern and southern hemispheres. Our results revealed that the conjugacy of propagation direction between the northern and southern hemispheres was subject to the influence of Coriolis force. We also figure out the evolution processes of ionospheric disturbances at the global scale. These are important topics that had not been well addressed previously. In addition, we also obtained thee wave structures of medium scale TIDs excited by the solar terminator (ST) moving over the northern America and physical mechanisms involved. Our observations confirm that the ST is a stable and repetitive source of ionospheric wave disturbances and the evidence of solar terminator generated disturbances has been demonstrated experimentally via the GPS TEC measurement. The main researches and results of this dissertation are as follows. First, the global traveling ionospheric disturbances (TIDs) during the drastic magnetic storms of October 29–31, 2003 were analyzed using the Global Position System (GPS) total electron content (TEC) data observed in the Asian-Australian, European and North American sectors. We collected the most comprehensive set of the TEC data from more than 900 GPS stations on the International GNSS Services (IGS) website and introduce here a strategy that combines polynomial fitting and multi-channel maximum entropy spectral analysis to obtain TID parameters. Moreover, in collaboration with my thesis advisor, I have developed an imaging technique of 2-dimensional map of TIDs structures to obtain spatial and temporal maps of large scale traveling ionospheric disturbances (LSTIDs). The clear structures of TEC perturbations map during the passage of TIDs were displayed. The results of our study are summarized as follows: (1) Large-scale TIDs (LSTIDs) and medium-scale TIDs (MSTIDs) were detected in all three sectors after the sudden commencement (SC) of the magnetic storm, and their features showed longitudinal and latitudinal dependences. The duration of TIDs was longer at higher latitudes than at middle latitudes, with a maximum of about 16 h. The TEC variation amplitude of LSTIDs was larger in the North American sector than in the two other sectors. At the lower latitudes, the ionospheric perturbations were more complicated, and their duration and amplitude were relatively longer and larger. (2) The periods and phase speeds of TIDs were different in these three sectors. In Europe, the TIDs propagated southward; in North America and Asia, the TIDs propagated southwestward; in the near-equator region, the disturbances propagated with the azimuth (the angle of the propagation direction of the LSTIDs measured clockwise from due north with 0°) of 210° showing the influence of Coriolis force; in the Southern Hemisphere, the LSTIDs propagated conjugatedly northwestward. Both the southwestward and northeastward propagating LSTIDs are found in the equatorial region. These results mean that the Coriolis effect cannot be ignored for the wave propagation of LSTIDs and that the propagation direction is correlated with the polar magnetic activity. (3) The day (day of year: 301) before the SC (sudden commencement) of magnetic storm, we observed a sudden TEC skip disturbances (±10 TECU). It should be a response for the high flux of proton during the solar flare event, but not the magnetic storms. Next, the most comprehensive and dense GPS network’s data from North-America region were used in this paper to analyze the medium scale traveling ionospheric disturbances (MSTIDs) which were generated by the moving solar terminator during the quiet days in 2005. We applied the multi-channel maximum entropy spectral analysis to calculated TID parameters, and found that the occurrence of ST-MSTIDs depends on the seasonal variations. The results of our study are summarized as follows: (1) MSTIDs stimulated by the moving ST (ST-MSTIDs) are detected at mid-latitudes after the passage of the solar terminator with the life time of 2～3 hours and the variation amplitude of 0.2～0.8 TECU. Spectral analysis indicated that the horizontal wavelength, average period, horizontal phase velocity of the MSTIDs are around 300±150 km，150±80 m/s and 25±15 min, respectively. In addition, ST-MSTIDs have wave fronts elongating the moving ST direction and almost parallel to ST. (2) The statistical results demonstrate that the dusk MSTIDs stimulated by ST is more obvious than the dawn MSTIDs in summer. On the contrary, the more-pronounced dawn MSTIDs occurs in winter. (3) Further analysis indicates that the seasonal variations of ST-MSTIDs occurrence frequency are most probably related to the seasonal differences of the variations of EUV flux in the ionosphere region and recombination process during sunrise and sunset period at mid-latitudes. Statistical study of occurrence characteristics of TIDs using the GPS network in North-American and European during solar maximum, In conclusion, statistical studies of the propagation characteristics of TIDs, which excited by the two common origins including geomagnetic storms and moving solar terminator, were involved with global GPS TEC databasein this thesis. We employed the multichannel maximum entropy spectral analysis method to diagnose the characteristics of propagation and evolvement of ionospheric disturbances, also, the characteristics of their regional distribution and climatological variations were revealed by the statistic analysis. The results of these studies can improve our knowledge about the energy transfer in the solar-terrestrial system and the coupling process between upper and lower atmosphere (thermosphere-ionosphere-mesosphere). On the other hand, our results of the investigation on TIDs generated by particular linear origin such as ST are important for developing ionospheric irregularity physics and modeling the transionosphere radio wave propagation. Besides, the GPS TEC representation of the ST-generated ionospheric structure suggests a better possibility for investigating this phenomenon. Subsequently, there are scientific meaning of the result of this dissertation to deeply discuss the energy transfer and coupling in the ionosphere, as well as realistic value to space weather forecast in the ionosphere region.

太阳辐射对热层和电离层变化性的影响

Solar ultraviolet (UV) radiation at wavelengths less than 400 nm is an important source of energy for aeronomic processes throughout the solar system. Solar UV photons are absorbed in planetary atmospheres, as well as throughout the heliosphere, via photodissociation of molecules, photoionization of molecules and atoms, and photoexcitation toexcitation including resonance scattering. In this paper, the solar irradiances data measured by TIMED SEE, as well as the solar proxies such as F10.7 and Mg II, thermosphere neutral density of CHAMP measurements and topside ionospheric plasmas densities from DMSP, are used to analyze solar irradiance effects on the variabilities of the thermosphere and the ionosphere. First, thermosphere densities near 410 km altitude are analyzed for solar irradiance variability effects during the period 2002-2004. Correlations between the densities and the solar irradiances for different spectral lines and wavelength ranges reveal significantly different characteristics. The density correlates remarkably well with all the selected solar irradiances except the lower chromospheric O I (130.4 nm) emission. Among the chosen solar proxies, the Mg II core-to-wing ratio index, EUV (30-120 nm) and F10.7 show the highest correlations with the density for short-term (< ~27 days) variations. For both long- (> ~27 days) and short-term variations, linear correlation coefficients exhibit a decreasing trend from low latitudes towards high latitudes. The density variability can be effectively modeled (capturing 71% of the variance) using multiple solar irradiance indices, including F10.7, SEUV (the EUV 30-120 nm index), and SFUV (the FUV 120-193 nm index), in which a lag time of 1 day was used for both F10.7 and SEUV, and 5 days for SFUV. In our regression formulation SEUV has the largest contribution to the density variation (40%), with the F10.7 having the next largest contribution (32%) and SFUV accounting for the rest (28%). Furthermore, a pronounced period of about 27.2 days (mean period of the Sun's rotation) is present in both density and solar irradiance data of 2003 and 2004, and a pronounced period of about 54.4 days (doubled period of the solar rotation) is also revealed in 2004. However, soft X-ray and FUV irradiances did not present a pronounced 54.4 day period in 2004, in spite of their high correlation with the densities. The Ap index also shows 54-day periodicities in 2004, and magnetic activity, together with solar irradiance, affects the 54-day variation in density significantly. In addition, NRLMSISE00, DTM-2000 and JB2006 model predictions are compared with density measurements from CHAMP to assess their accuracy, and the results show that these models underestimate the response of the thermosphere to variations induced by solar rotation. Next, the equatorial topside ionospheric plasmas densities Ni are analyzed for solar irradiance variability effects during the period 2002-2005. Linear correlations between Ni and the solar irradiances for different wavelength ranges reveal significantly different characteristics. XUV (0-35 nm) and EUV (115-130 nm) show higher correlation with Ni for the long-term variations, whereas EUV (35-115 nm) show higher correlation for the short-term variations. Moreover, partial correlation analysis shows that the long-term variations of Ni are affected by both XUV (0-35 nm) and EUV (35-115 nm), whereas XUV (0-35 nm) play a more important role; the short-term variations of Ni are mostly affected by EUV (35-115 nm). Furthermore, a pronounced period of about 27 days is present in both Ni and solar irradiance data of 2003 and 2004, and a pronounced period of about 54 days is also revealed in 2004. Finally, prompted by previous studies that have suggested solar EUV radiation as a means of driving the semiannual variation, we investigate the intra-annual variation in thermosphere neutral density near 400 km during 2002-2005. The intra-annual variation, commonly referred to as the ‘semiannual variation’, is characterized by significant latitude structure, hemispheric asymmetries, and inter-annual variability. The magnitude of the maximum yearly difference, from the yearly minimum to the yearly maximum, varies by as much as 60% from year to year, and the phases of the minima and maxima also change by 20-40 days from year to year. Each annual harmonic of the intra-annual variation, namely, annual, semiannual, ter-annual and quatra-annual, exhibits a decreasing trend from 2002 through 2005 that is correlated with the decline in solar activity. In addition, some variations in these harmonics are correlated with geomagnetic activity, as represented by the daily mean value of Kp. Recent empirical models of the thermosphere are found to be deficient in capturing most of the latitude dependencies discovered in our data. In addition, the solar flux and geomagnetic activity proxies that we have employed do not capture some latitude and inter-annual variations detected in our data. It is possible that these variations are partly due to other effects, such as seasonal-latitudinal variations in turbopause altitude (and hence O/N2 composition) and ionosphere coupling processes that remain to be discovered in the context of influencing the intra-annual variations depicted here. Our results provide a new dataset to challenge and validate thermosphere-ionosphere general circulation models that seek to delineate the thermosphere intra-annual variation and to understand the various competing mechanisms that may contribute to its existence and variability. We furthermore suggest that the term “intra-annual” variation be adopted to describe the variability in thermosphere and ionosphere parameters that is well-captured through a superposition of annual, semiannual, ter-annual, and quatra-annual harmonic terms, and that “semiannual’ be used strictly in reference to a pure 6-monthly sinusoidal variation. Moreover, we propose the term “intra-seasonal” to refer to those shorter-term variations that arise as residuals from the above Fourier representation.

电离层参量M(3000)F2和hmF2的经验正交函数方法建模研究

The ionospheric parameter M(3000)F2 (the so-called transmission factor or the propagation factor) is important not only in practical applications such as frequency planning for radio-communication but also in ionospheric modeling. This parameter is strongly anti-correlated with the ionospheric F2-layer peak height hmF2，a parameter often used as a key anchor point in some widely used empirical models of the ionospheric electron density profile (e.g., in IRI and NeQuick models). Since hmF2 is not easy to obtain from measurements and M(3000)F2 can be routinely scaled from ionograms recorded by ionosonde/digisonde stations distributed globally and its data has been accumulated for a long history, usually the value of hmF2 is calculated from M(3000)F2 using the empirical formula connecting them. In practice, CCIR M(3000)F2 model is widely used to obtain M(3000)F2 value. However, recently some authors found that the CCIR M(3000)F2 model has remarkable discrepancies with the measured M(3000)F2, especially in low-latitude and equatorial regions. For this reason, the International Reference Ionosphere (IRI) research community proposes to improve or update the currently used CCIR M(3000)F2 model. Any efforts toward the improvement and updating of the current M(3000)F2 model or newly development of a global hmF2 model are encouraged. In this dissertation, an effort is made to construct the empirical models of M(3000)F2 and hmF2 based on the empirical orthogonal function (EOF) analysis combined with regression analysis method. The main results are as follows: 1. A single station model is constructed using monthly median hourly values of M(3000)F2 data observed at Wuhan Ionospheric Observatory during the years of 1957–1991 and compared with the IRI model. The result shows that EOF method is possible to use only a few orders of EOF components to represent most of the variance of the original data set. It is a powerful method for ionospheric modeling. 2. Using the values of M(3000)F2 observed by ionosondes distributed globally, data at grids uniformly distributed globally were obtained by using the Kriging interpolation method. Then the gridded data were decomposed into EOF components using two different coordinates: (1) geographical longitude and latitude; (2) modified dip (Modip) and local time. Based on the EOF decompositions of the gridded data under these two coordinates systems, two types of the global M(3000)F2 model are constructed. Statistical analysis showed that the two types of the constructed M(3000)F2 model have better agreement with the observational M(3000)F2 than the M(3000)F2 model currently used by IRI. The constructed models can represent the global variations of M(3000)F2 better. 3. The hmF2 data used to construct the hmF2 model were converted from the observed M(3000)F2 based on the empirical formula connecting them. We also constructed two types of the global hmF2 model using the similar method of modeling M(3000)F2. Statistical analysis showed that the prediction of our models is more accurate than the model of IRI. This demonstrated that using EOF analysis method to construct global model of hmF2 directly is feasible. The results in this thesis indicate that the modeling technique based on EOF expansion combined with regression analysis is very promising when used to construct the global models of M(3000)F2 and hmF2. It is worthwhile to investigate further and has the potential to be used to the global modeling of other ionospheric parameters.

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

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.