磁层亚暴期间高纬地区三维电流体系

确定亚暴期间高纬地区的三维电流体系是磁层和电离层物理的基本问题之一,本文简述了近年来发展的根据地面地磁活动的记录及演三维电流体系的两种较新的方法:KRM（Kamide-Richmond-Matsushita）方法和K（Kisabeth）方法.在KRM方法和许多其它类似的计算中,都假定磁力线是垂直于电离层沿地球径向的直线.本文介绍了一种递推方法,可计算场向电流沿弯曲的偶极子场线流动的情况。同时,还介绍了高纬三维电流体系对亚暴期间中低纬度地磁扰动的贡献。最后介绍了在计算电流体系时所需的电离层电导率模式。     

无线电掩星反演大气参数误差分析及应用进展

本文详细分析了影响空间GPS无线电掩星反演大气参数方法的主要误差源，其中包括信噪比误差，精密定轨误差，电离层模制误差，多路径效应和时钟改正等误差，并且介绍了国内外在无线电掩星领域的一些应用研究进展。     

Three-dimensional finite-element modelling of magnetotelluric data with a divergence correction

A finite-element method for computing the electric field in a 3-D conductivity model of the Earth for plane wave sources, thus enabling magnetotelluric responses to be calculated, is presented. The method incorporates in the iterative solution of the electric-field system of equations the divergence correction technique introduced for finite-difference solutions by Smith (1996). The correction technique accelerates the development of the discontinuity of the normal component of the approximate electric field across conductivity discontinuities. The convergence rate of the iterative solution is improved significantly, especially for low frequencies. The correction technique involves computing the divergence of the current density for the approximate electric field, computing the static potential whose source is this divergence of the current density, and ‘correcting’ the approximate electric field by subtracting from it the gradient of the potential. This is repeated at regular intervals during the iterative solution of the electric-field system of equations. For the method presented here, the Earth model is discretised using a rectilinear mesh comprising uniform cells. Edge-element basis functions are used to approximate the electric field and nodal basis functions are used to approximate the correction potential. The Galerkin method is used to derive the systems of equations for the approximate electric field and correction potential from the respective differential equations. A bi-conjugate gradient solver was found to be adequate for the system of equations for the correction potential; a generalised minimum residual solver was found to be better for the electric-field system of equations. The method is illustrated using the COMMEMI 3D-1A and 3D-2A models.     

Effect of horizontal gradients on ionospherically reflected or transionospheric paths using a precise homing-in method

A homing-in method is presented for determining ionospheric reflected or transionospheric paths between fixed transmitter and receiver locations in the presence of ionospheric gradients or ripples. Both initial elevation and azimuth are automatically adjusted to find the path that arrives exactly at the receiver. The method can be used for any 3D ionospheric model to find precise ray paths and phase and group delays for both magneto-ionic modes. The method takes full account of path location, geomagnetic field orientation and the bending of the ray path resulting from horizontal as well as vertical gradients of electron density. It can also find multiple paths e.g. low and high angle, 1- and 2-hops for both ordinary and extraordinary modes. Examples of its use are given for both terrestrial HF links and Earth to Satellite paths. For paths reflected from the ionosphere, the effect of gradients of both critical frequency and height of maximum electron density are determined and the comparative effect of gradients on high and low angle and 1- and 2-hops paths for both magneto-ionic modes investigated. Path variation with frequency for a fixed link is also studied and the bandwidth of the ionospheric background channel (dispersive bandwidth) and its reciprocal (the pulse rise time), important for wideband digital HF broadcasting or spread spectrum HF communications, is estimated for a range of frequencies, for high- and low-angle rays and 1- and 2-hop paths. For Earth–satellite paths, the effect of the ionosphere and horizontal ionospheric gradients is determined for a range of frequencies and elevation angles. It is shown that the method can also enable the determination of second-order errors in satellite navigation methods, such as GPS, due to ionospheric gradients and the effect of the geomagnetic field.     

Influence of seasonal variation of neutral winds on the coupling between North and South Hemispheres in a realistic Earth main field

The neutral wind dynamo contributes significantly to the ionospheric electrodynamics, and the variation of the neutral winds thus affects the ionosphere. Here we study the effects of the seasonal variation of the winds in a realistic Earth main field. The two-dimensional ionospheric dynamo equation is expressed in the framework of a revision of the International Geomagnetic Reference Field (IGRF) under the assumptions of equipotential field lines and conservation of current. A revision of IGRF and typical u...     

Geomagnetic field effect on the ionospheric contribution to the error of navigation satellite systems

When using global navigation satellite systems (GNSSs) for high-precision measurements, one should consider high-order errors. The ionospheric second-order error caused by the geomagnetic field is approximately proportional to the total electron content. Therefore, this error can be taken into account by modifying the coefficients in an “ionosphere-free” combination of GNSS measurements at two frequencies. This study checks the approximations underlying this modification. We reveal that these approximations are valid and the results depend weakly on the accuracy of ionospheric parameters used a priori for calculating the coefficients of the modified two-frequency formula. In addition, we investigate how the choice of a model of the Earth’s magnetic field affects the second-order ionospheric error.     

A Method for the Ionospheric Delay Estimation and Interpolation in a Local GPS Network

To estimate ionospheric delays from the Global Positioning System (GPS) measurements, satellite and receiver equipment biases have to be modeled. This paper presents a procedure based on the least squares (LS) approach, which implicitly takes into account these equipment biases in the estimation of the ionospheric effect. The second part of this work deals with the interpolation of the ionospheric correction from a permanent GPS network to a single frequency GPS user. The results obtained show that for 10-cm position accuracy the ionospheric delay can be successfully interpolated when the GPS user is within 40 km of the GPS permanent network.     

Auroral ionospheric and thermospheric measurements using the incoherent scatter technique

The incoherent scatter technique has been applied since 1965 to study the ionosphere and thermosphere in different regions of the Earth. The analysis of the received signal gives access to several ionospheric parameters as a function of height: electron density, electron and ion temperatures and ion velocity. The derivation of these parameters is usually a complicated mathematical procedure that requires a non-linear regression program. A lot of research has been done in the ionospheric and atmospheric science using this technique.In this paper we describe how one derives the ion-neutral collision frequency and the ion composition parameters. It is usually difficult to retrieve these parameters with the incoherent scatter technique; as a result, in the standard data analysis procedure, an ionospheric model is used instead. However the numerical values chosen in the model have an influence on the other derived parameters. For instance the choice of a wrong ion composition leads to erroneous plasma temperatures. It is therefore important to assess by how much the standard procedure deviates from reality. For this reason we compare the ion composition and collision frequency retrieved from a sophisticated analysis scheme with the values that are derived from models under similar geophysical conditions.It also possible to derive from the observed ionospheric parameters the neutral concentrations, temperatures and winds, by using the energy and momentum equations for the ions and the neutrals. In this paper the different methods and the corresponding assumptions involved in the data analysis are discussed. We describe the influence of the frictional heating, of the vertical neutral wind and of the ionospheric perturbations on the derivation of the neutral atmospheric parameters. Our discussion of the processes involved are drawn from results obtained by Chatanika, Sondrestrom and EISCAT radars.     

Ionospheric specification and forecasting based on observations from European ionosondes participating in DIAS project

There are two main objectives of the DIAS (European Digital Upper Atmosphere Server) project. First, it establishes a pan-European repository of raw and derived digital data describing the state of ionospheric part of the upper atmosphere, which is capable of ingesting real-time information and maintaining historical data collections provided by most operating ionospheric stations in Europe. Second, the DIAS system produces and distributes, based on the raw data collection, several operational products required by various user groups for ionospheric nowcasting, prediction and forecasting purposes. The project completed on May 2006 and the DIAS server operates since then continuously. The basic products that are delivered are real-time and historical ionograms from all DIAS-affiliated ionospheric stations, frequency plots and maps of the ionosphere over Europe based on the foF2, M(3000)F2, MUF and electron density parameters, as well as long term prediction and short term forecasting up to 24 hour ahead. The paper describes use of the ionospheric measurements in the DIAS modelling techniques for specification, predict-tion and forecasting of the ionosphere over the European region, and details the final products available to the DIAS user community.     

电离层暴时的foF2预报技术研究

电离层是一个非常复杂的非线性系统.本文利用BP网络非线性输入-输出映射的特点预报电离层F2层临界频率.由于foF2在太阳活动高低年、不同经纬度、不同季节有着不同变化特点，本文仅以海口站、长春站的数据为例分别对1994年和2001年的foF2进行预报，预报结果表明发生电离层暴时用神经网络法预报foF2比国际参考电离层更有优势.此外，由于不同用户对预报精度的需求不同，本文以海口站2000、2001年的预报结果为例构建了残差模型.     

Assessment of NeQuick ionospheric model for Galileo single-frequency users

The ionosphere is the main error source in GNSS measurements and in extreme cases can degrade the positioning significantly, with errors exceeding 100 m; therefore, modelling and predicting of this type of error is crucial and critical. The ionospheric effect can be reduced using different techniques, such as dual-frequency receiver or suitable augmentation system (DGPS, SBAS); the aforesaid approaches involve the use of expensive devices and/or complex architectures. Single frequency stand-alone receivers are the cheapest and most widespread GNSS devices; they can estimate and partially correct the error due to the ionosphere, through adequate algorithms, which use parameters broadcasted by the navigation message. The aim of this paper is performance assessment of the ionospheric model NeQuick, adopted by the European GNSS Galileo for single frequency receivers. The analysis is performed in measurements domain and the data are collected in different geographical locations and in various geomagnetic conditions.     

Ionospheric quiet-day currents and their effects on the interpretation of upper mantle conductivity – A review

A review of ionospheric currents,the external part of geomagnetic solar quiet (Sq)variation and the internal part of the induced currentwithin the Earth, has been carried out.The theoretical background has been reviewed,and the method of analysis, specificallythe spherical harmonic analysis (SHA),is presented. Various works are reviewed andfully discussed. It has been noted that anew model of the Sq current system could beused to obtain a mantle electrical conductivityprofile for hemispheres. It is suggested thatfurther research work be carriedout in other hemispheres using similarmodels for more robust interpretations.     

Unbiased total electron content (UTEC), their fluctuations,and correlation with seismic activity over Japan

Earthquakes are among the most dangerous events that occur on earth and many scientists have been investigating the underlying processes that take place before earthquakes occur. These investigations are fueling efforts towards developing both single and multiple parameter earthquake forecasting methods based on earthquake precursors. One potential earthquake precursor parameter that has received significant attention within the last few years is the ionospheric total electron content (TEC). Despite its growing popularity as an earthquake precursor, TEC has been under great scrutiny because of the underlying biases associated with the process of acquiring and processing TEC data. Future work in the field will need to demonstrate our ability to acquire TEC data with the least amount of biases possible thereby preserving the integrity of the data. This paper describes a process for removing biases using raw TEC data from the standard Rinex files obtained from any global positioning satellites system. The process is based on developing an unbiased TEC (UTEC) data and model that can be more adaptable to serving as a precursor signal for earthquake forecasting. The model was used during the days and hours leading to the earthquake off the coast of Tohoku, Japan on March 11, 2011 with interesting results. The model takes advantage of the large amount of data available from the GPS Earth Observation Network of Japan to display near real-time UTEC data as the earthquake approaches and for a period of time after the earthquake occurred.     

利用精密单点定位求解电离层延迟

近年来,高时空分辨率的全球导航卫星系统(GNSS)观测信号已成为电离层研究的重要资源.利用GNSS研究电离层,需首先将观测资料转换成包含电离层信息的可观测量(Ionospheric 0bservables,称之为"电离层观测值").目前,最常用的电离层观测值一般采用联合无几何影响组合的码和相位观测,利用相位平滑伪距方法...     

基于陆态网络GPS数据的电离层空间天气监测与研究

中国大陆构造环境监测网络(简称陆态网络)是以全球卫星导航定位系统(GNSS)为主,辅以多种空间观测技术,实时动态监测大陆构造环境变化,探求其对资源、环境和灾害的影响的地球科学综合观测网络.基于陆态网络约200个基准站的GPS观测数据,本文探讨了其在电离层空间天气监测与研究方面的应用.包括磁暴期间电离层暴扰动形态,大尺度电离层行进式扰动,太阳耀斑引起的电离层骚扰和低纬电离层不规则体结构等.研究结果表明:陆态网络布局合理,观测数据质量良好,完全可用于中国及周边地区电离层空间天气监测与研究,为进一步开展我国电离层空间天气预警和预报奠定了观测基础.     

Interdependent perturbations of the Earth’s surface, atmosphere, and ionosphere

The results of original experiments performed with a ground-based geophysical laser interferometer and a GPS-based satellite ionospheric profilometer are given. Synchronous growth was recorded for deformations of the Earth’s surface and variations in the atmospheric pressure and in the level of spatiotemporal modifications of the electron content within the ionospheric F2 layer with characteristic space scales of 10²–10³ km and periods of 10²–10³ s. The relationship between the revealed phenomena and the Earth’s seismic activity is analyzed.     

Magnetic Pulsations: Their Sources and Relation to Solar Wind and Geomagnetic Activity

Ultra low frequency (ULF) waves incident on the Earth are produced by processes in the magnetosphere and solar wind. These processes produce a wide variety of ULF hydromagnetic wave types that are classified on the ground as either Pi or Pc pulsations (irregular or continuous). Waves of different frequencies and polarizations originate in different regions of the magnetosphere. The location of the projections of these regions onto the Earth depends on the solar wind dynamic pressure and magnetic field. The occurrence of various waves also depends on conditions in the solar wind and in the magnetosphere. Changes in orientation of the interplanetary magnetic field or an increase in solar wind velocity can have dramatic effects on the type of waves seen at a particular location on the Earth. Similarly, the occurrence of a magnetospheric substorm or magnetic storm will affect which waves are seen. The magnetosphere is a resonant cavity and waveguide for waves that either originate within or propagate through the system. These cavities respond to broadband sources by resonating at discrete frequencies. These cavity modes couple to field line resonances that drive currents in the ionosphere. These currents reradiate the energy as electromagnetic waves that propagate to the ground. Because these ionospheric currents are localized in latitude there are very rapid variations in wave phase at the Earth’s surface. Thus it is almost never correct to assume that plane ULF waves are incident on the Earth from outer space. The properties of ULF waves seen at the ground contain information about the processes that generate them and the regions through which they have propagated. The properties also depend on the conductivity of the Earth underneath the observer. Information about the state of the solar wind and the magnetosphere distributed by the NOAA Space Disturbance Forecast Center can be used to help predict when certain types and frequencies of waves will be observed. The study of ULF waves is a very active field of space research and much has yet to be learned about the processes that generate these waves.     

GNSS测站坐标非线性变化影响机制的研究

对影响地球参考框架稳定度与精确度的地球动力学因素进行深入研究,能更好地促进毫米级地球参考架的建立.而研究GNSS测站坐标的非线性变化及其相关机制,是实现毫米级地球参考框架的前提,也为提升GNSS坐标精度、研究地球物理因素等提供参考.首先研究了GNSS测站坐标非线性变化的影响机制与理论改正模型,给出了季节性温度效应、地表物质负荷效应、GNSS周年性系统误差和高阶电离层延迟这四种因素的研究进展,并提出了在机制和模型研究方面一些改进的方向.然后通过以高阶电离层延迟对GNSS测站坐标时间序列的影响为例,进行分组实验,研究了高阶电离层延迟对GNSS测站坐标非线性变化的影响,并进一步将影响进行了量化.研究结果有助于对测站坐标非线性变化的机制做出更精确合理的解释,同时也有助于削弱与限制高阶电离层延迟对建立高精度地球参考框架的影响;为我国北斗地球参考框架的建立与维持提供了可靠的理论基础,具有一定的参考价值.     

The influence of the color of visibility targets on the visibility

Reduced visibilities have been simulated in the laboratory by means of a hydrosol. It was illuminated homogeneously by iodine quartz lamps producing an illuminance similar to a bright day. In the hydrosol several visibility targets were suspended at different distances and the visibility observation was performed by counting the number of visible objects. By this method a very accurate and fast determination of the visibility was possible. The extinction coefficient of the hydrosol was determined with a long path photometer: the intrinsic brightness of the objects in the illuminated hydrosol was determined for several wavelengths by means of a fiber optic. The correlation between the measured visibility and the extinction coefficient shows that the frequently used visual range underestimates the visibility, since the highest perception for an atmospheric aerosol is not at 550 nm but at 580 nm wavelength. Therefore it is suggested to use the atmospheric extinction coefficient at 580 nm for the calculation of the visual range. All colored objects could be seen less far, since they had a smaller contrast at the wavelength of maximum perception than the black object. Dark colored objects can be seen further than bright ones. Especially objects that have a high reflection in the yellow and organe have a small visibility. If the intrinsic brightness of the objects and the extinction coefficient (both as a function of wavelength) are known, the wavelength of maximum perception can be calculated and from this the visibility. The visibility thus calculated agrees with the observed visibility. If visibility observations in the atmosphere have to be performed with non-black objects, methods for corrections to the visibilty of a black object are given.     

A methodology for using optimal MSIS parameters retrieved from SSULI data to compute satellite drag on LEO objects

A key application to be derived from Space Weather research will be to forecast atmospheric drag on low Earth orbit (LEO) satellites with significantly better accuracy than is attainable today. The recently launched STP P91-1 ARGOS mission will serve as a testbed for the use of future operational ultraviolet remote-sensing data to achieve such an improvement. This paper describes the associated methodology, which uses discrete inverse theory in conjunction with the data to derive correction factors in near real time for the MSISE-90 empirical thermospheric model. To simulate the application of this technique to orbit prediction, we use the Jacchia-71 operational model to generate an evolving “ground truth” upper atmospheric state over a 48 h time period. This permits a state-of-the-art Satellite Tool Kit orbit propagator to synthesize a corresponding “ground truth” orbit on a standard LEO test object at 350 km altitude. Our tests show that, for orbit prediction, the “data-enhanced” MSIS density specification can provide significant improvement over the uncorrected MSIS specification. However, for orbit prediction, the results are sensitive to the strategy selected for applying the correction factors. We contrast our results for orbit prediction with those of Marcos et al. (1998. Astrodynamics. Vol. 97(1). AAS, San Diego, pp. 501–513) for precision orbit determination. An important result in the context of Space Weather is that the Jacchia and MSIS models can show significant point-to-point disagreement, which has major implications for operational specification of thermospheric drag.