Evaluation of the dynamic cutoff rigidity model using dosimetry data from the STS-28 flight.

We have a developed a dynamic cutoff rigidity model based on computed world grids of vertical cutoff rigidities derived from employing the Tsyganenko magnetospheric model. The dynamic range of this model covers all magnetic activity levels specified by integer values of the Kp magnetic index. We present comparisons of the measured dose observed on the space shuttle during the August 1989 solar proton event with the dose computed from solar particles predicted to be allowed through the magnetosphere to the space shuttle position. We find a one-to-one correspondence between the portion of the orbit predicted to be subjected to solar protons and the portion of the orbit where solar particle dose measurements were obtained.     

The space-developed dynamic vertical cutoff rigidity model and its applicability to aircraft radiation dose.

We have developed a dynamic geomagnetic vertical cutoff rigidity model that predicts the energetic charged particle transmission through the magnetosphere. Initially developed for space applications, we demonstrate the applicability of this library of cutoff rigidity models for computing aircraft radiation dose. The world grids of vertical cutoff rigidities were obtained by particle trajectory tracing in a magnetospheric model. This reference set of world grids of vertical cutoff rigidities calculated for satellite altitudes covers all magnetic activity levels from super quiet to extremely disturbed (i.e., Kp indices ranging from 0 to 9+) for every three hours in universal time. We utilize the McIlwain "L" parameter as the basis of the interpolation technique to reduce these initial satellite altitude vertical cutoff rigidities to cutoff rigidity values at aircraft altitudes.     

Variations in cosmic ray cutoff rigidities during the great geomagnetic storm of November 2004

A very strong interplanetary and magnetospheric disturbance observed on 7–13 November 2004 can be regarded as one of the strongest events during the entire period of space observations. In this paper we report on the studies of cosmic ray cutoff rigidity variations during 7–13 November 2004 showing how storm conditions can affect the direct cosmic ray access to the inner magnetosphere. Effective cutoff rigidities have been calculated for selected points on the ground by tracing trajectories of cosmic ray particles through the magnetospheric magnetic field of the “storm-oriented” Tsyganenko 2003 model. Cutoff rigidity variations have also been determined by the spectrographic global survey method on the basis of experimental data of the neutron monitor network. Relations between the calculated and experimental cutoff rigidities and the geomagnetic Dst-index and interplanetary parameters have been investigated. Correlation coefficients between the cutoff rigidities obtained by the trajectory tracing method and the spectrographic global survey method have been found to be in the limits 0.76–0.89 for all stations except the low-latitude station Tokyo (0.35). The most pronounced correlation has been revealed between the cutoff rigidities that exhibited a very large variation of ∼1–1.5 GV during the magnetic storm and the Dst index.     

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

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

On the possibility to check the magnetosphere’s model by CR: The strong geomagnetic storm in November 2003

We report on a study of cosmic ray cutoff rigidity variations during the strong geomagnetic storm of 18–24 November 2003. We employed the most recent Tsyganenko magnetospheric model to represent the very strong disturbed Magnetosphere. We used this magnetic field for the cosmic ray trajectory calculations to determine the geomagnetic cutoff rigidity throughout this period of severe geomagnetic disturbance. We determine the cutoff rigidity changes during this period by two methods, by trajectory calculations and by the spectrographic global survey method. The values of the change of cutoff rigidities obtained by two different methods are correlated with the D_st and interplanetary magnetic field and plasma parameters and result in correlation coefficients in the range 0.63–0.84 for the various cosmic ray stations. The result of this study indicates that the most significant contributions to the cutoff rigidity changes are due to D_st variation although the influence of solar wind density and B_z and B_y components of IMF variations is significant.     

磁扰期间的三维漂移壳分离

基于磁层粒子动力学理论,首先对比了计算漂移壳分离的引导中心法和磁力线追踪法,计算表明两种方法的计算结果一致.然后分别采用T89c和T96磁层磁场模式,用磁力线追踪法数值计算了不同初始位置(≤9Re)、不同初始投掷角、不同Kp指数和不同太阳风压力下,带电粒子的漂移壳分离.计算结果揭示了漂移壳分离随初始位置、投掷角、Kp指数和太阳风压力的变化.其具体特征如下. (1)随着径向距离的增大,漂移壳分离效应愈加显著,由正午出发的粒子将被稳定捕获,而午夜出发的径向距离≥7Re的部分大投掷角粒子将沿磁层顶逃逸. (2)正午出发的粒子,漂移到午夜时其漂移壳随投掷角减小向外排列;午夜出发的粒子,漂移到正午时其漂移壳随投掷角增大排列; 90°投掷角粒子在磁赤道面的漂移壳沿着磁场等值线排列. (3)漂移壳分离随Kp指数和太阳风压力增大变得显著,且随这两种扰动参数的变化特征和趋势是基本相似的.      

Geomagnetic cutoffs: a review for space dosimetry applications.

The earth's magnetic field acts as a shield against charged particle radiation from interplanetary space, technically described as the geomagnetic cutoff. The cutoff rigidity problem (except for the dipole special case) has "no solution in closed form". The dipole case yields the Stormer equation which has been repeatedly applied to the earth in hopes of providing useful approximations of cutoff rigidities. Unfortunately the earth's magnetic field has significant deviations from dipole geometry, and the Stormer cutoffs are not adequate for most applications. By application of massive digital computer power it is possible to determine realistic geomagnetic cutoffs derived from high order simulation of the geomagnetic field. Using this technique, "world-grids" of directional cutoffs for the earth's surface and for a limited number of satellite altitudes have been derived. However, this approach is so expensive and time consuming it is impractical for most spacecraft orbits, and approximations must be used. The world grids of cutoff rigidities are extensively used as lookup tables, normalization points and interpolation aids to estimate the effective geomagnetic cutoff rigidity of a specific location in space. We review the various options for estimating the cutoff rigidity for earth-orbiting satellites.     

基于机器学习相似度算法的Kp指数预报

基于机器学习中的相似度算法,建立了在历史太阳风数据中寻找与当前太阳风特征相近事例的推荐模型,用来预报地磁Kp指数.使用1998-2019年间随机选择的120个太阳风事例作为测试数据集,该模型能够推荐得到历史上与输入太阳风造成相似地磁影响的太阳风事例,最优事例的Kp指数与实际值的均方根误差为0.79,相关系数为0.93....     

Simultaneous observations of the plasma density on the same field line by the CPMN ground magnetometers and the Cluster satellites

By applying the cross-phase method and the amplitude-ratio method to magnetic field data obtained from two ground stations located close to each other, we can determine the frequency of the field line resonance (FLR), or the field line eigenfrequency, for the field line running through the midpoint of the two stations. From thus identified FLR frequency we can estimate the equatorial plasma mass density (ρ)$(\rho )$ by using the T05s magnetospheric field model [Tsyganenko, N.A., Sitnov, M.I. Modeling the dynamics of the inner magnetosphere during strong geomagnetic storms, J. Geophys. Res. 110, A03208, 2005] and the equation of Singer et al. [Singer, H.J., Southwood, D.J., Walker, R.J., Kivelson, M.G. Alfven wave resonances in a realistic magnetospheric magnetic field geometry, J. Geophys. Res. 86 (A6) 4589–4596, 1981].     

Modeling solar proton access to geostationary spacecraft with geomagnetic cutoffs

Solar energetic particle (SEP) cutoffs at geosynchronous orbit are sensitive to moderate geomagnetic activity and undergo daily variations due to the day–night asymmetry of the magnetosphere. At geosynchronous orbit, cutoff rigidity also has a large directional dependence, with the highest cutoff rigidity corresponding to ions arriving from magnetic east and lowest cutoff rigidity corresponding to ions incident from the west. Consequently, during geomagnetically quiet periods, the SEP flux observed by an eastward facing particle detector is significantly lower than observed by a westward facing particle detector. During geomagnetically disturbed periods the cutoff is suppressed allowing SEPs access well inside of geosynchronous, so that the east–west SEP flux ratio approaches unity. Variations in the east–west SEP flux ratio observed by GOES Energetic Particle Sensors (EPS) have recently been reported by Rodriguez et al. (2010). In NOAA’s operational processing of EPS count rates into differential fluxes, the differential flux is treated as isotropic and flat over the energy width of the channel. To compare modeled SEP flux with GOES EPS observations, the anisotropy of the flux over the EPS energy range and field of view must be taken into account. A technique for making direct comparisons between GOES EPS observations and SEP flux modeled using numerically computed geomagnetic cutoffs is presented. Initial results from a comparison between modeled and observed flux during the 6–11 December 2006 SEP event are also presented. The modeled cutoffs reproduce the observed flux variations well but are in general too high.     

Case study of September 24–26, 1998 magnetic storm

We compute global magnetospheric parameters based upon solar wind data obtained from the WIND spacecraft upstream. Using the paraboloid magnetospheric model, calculations of the dynamic global magnetospheric current systems have been made. The solar wind dynamic pressure, the interplanetary magnetic field, the strength of the tail current, and the ring current control the polar cap and auroral oval size and location during the magnetic storm. The model calculations demonstrate that the polar cap and the auroral oval areas are mainly controlled by the tail current. The substorm onset at 0630 UT on September 25, 1998 happened near the minimum in the main phase field depression. The substorm expansion onset time is also marked by a sudden enhancement in the solar wind dynamic pressure and an enhancement in the tail current. The magnetic signatures of these two effects cancel each other, which explains why the D_st profile shows no strong time variation during the substorm. Evidence for the substorm expansion includes not only the signature in the AL index but also the strong asymmetry of the low latitude magnetic disturbances (substorm positive bay signature). Model calculations were checked by comparison with the GOES 8 and 10 magnetic field measurements.     

The Earth’s magnetosphere response to interplanetary medium conditions on January 21–22, 2005 and on December 14–15, 2006

The Earth’s magnetosphere response to interplanetary medium conditions on January 21–22, 2005 and on December 14–15, 2006 has been studied. The analysis of solar wind parameters measured by ACE spacecraft, of geomagnetic indices variations, of geomagnetic field measured by GOES 11, 12 satellites, and of energetic particle fluxes measured by POES 15, 16, 17 satellites was performed together with magnetospheric modeling based in terms of A2000 paraboloid model. We found the similar dynamics of three particle populations (trapped, quasi-trapped, and precipitating) during storms of different intensities developed under different external conditions: the maximal values of particle fluxes and the latitudinal positions of the isotropic boundaries were approximately the same. The main sources caused RC build-up have been determined for both magnetic storms. Global magnetospheric convection controlled by IMF and substorm activity driven magnetic storm on December 14–15, 2006. Extreme solar wind pressure pulse was mainly responsible for RC particle injection and unusual January 21, 2005 magnetic storm development under northward IMF during the main phase.     

Geomagnetic transmission disturbances and heavy-ion fluences observed in low Earth orbit during the solar energetic particle events of October 1989.

The large solar energetic particle (SEP) events and simultaneous large geomagnetic disturbances observed during October 1989 posed a significant, rapidly evolving space radiation hazard. Using data from the GOES-7, NOAA-10, IMP-8 and LDEF satellites, we determined the geomagnetic transmission, heavy ion fluences, mean Fe ionic charge state, and effective radiation hazard observed in low Earth orbit (LEO) for these SEPs. We modeled the geomagnetic transmission by tracing particles through the combination of the internal International Geomagnetic Reference Field (IGRF) and the Tsyganenko (1989) magnetospheric field models, extending the modeling to large geomagnetic disturbances. We used our results to assess the radiation hazard such very large SEP events would pose in the anticipated 52 degrees inclination space station orbit.     

Neutron monitor asymptotic directions of viewing during the event of 13 December 2006

During the recent ground level enhancement of 13 December 2006, also known as GLE70, solar cosmic ray particles of energy bigger that ∼500 MeV/nucleon propagated inside the Earth’s magnetosphere and finally accessed low-altitude satellites and ground level neutron monitors. The magnitude and the characteristics of this event registered at different neutron monitor stations of the worldwide network can be interpreted adequately on the basis of an estimation of the solar particle trajectories in the near Earth interplanetary space. In this work, an extended representation of the Earth’s magnetic field was realized applying the Tsyganenko 1989 model. Using a numerical back-tracing technique the solar proton trajectories inside the magnetospheric field of the Earth were calculated for a variety of particles, initializing their travel at different locations, covering a wide range of energies. In this way, the asymptotic directions of viewing were calculated for a significant number of neutron monitor stations, providing crucial information on the Earth’s “magnetospheric optics” for primary solar cosmic rays, on the top of the atmosphere, during the big solar event of December 2006. The neutron monitor network has been treated, therefore, as a multidimensional tool that gives insights into the arrival directions of solar cosmic ray particles as well as their spatial and energy distributions during extreme solar events.     

Prediction of Magnetospheric Disturbances Caused by a Quasi-Stationary Solar Wind

When predicting parameters of quasi-stationary Solar Wind (SW) streams at 1 AU, it is customary to use, as the indicator of solar sources, the Bases of Open Magnetic Tubes (BOMT) on the solar surface obtained via a calculation relying on a new Bd-technique of harmonic expansion of the magnetic field from daily full-disk magnetograms developed by Rudenko[4]. By considering an example of 17 events, it is shown that the correspondence between fast SW streams at the Earth's orbit and the BOMT, calculated with ≤ 24 h time resolution, makes up about 94%, while the correspondence of SW stereams with the CH in the light of the 10830 A line is about 29%. With this technique, the predictability of maxima of the Kp index of magnetospheric disturbance caused by a fast quasi-stationary SW, is over 90%, and the prediction accuracy of the maximun velocity vm of the stream is ±15%.      

Evaluation of geomagnetic field models using magnetometer measurements for satellite attitude determination system at low earth orbits: Case studies

In this study, different geomagnetic field models are compared in order to study the errors resulting from the representation of magnetic fields that affect the satellite attitude system. For this purpose, we used magnetometer data from two Low Earth Orbit (LEO) spacecraft and the geomagnetic models IGRF-12 (Thébault et al., 2015) and T89 (Tsyganenko, 1989) models to study the differences between the magnetic field components, strength and the angle between the predicted and observed vector magnetic fields. The comparisons were made during geomagnetically active and quiet days to see the effects of the geomagnetic storms and sub-storms on the predicted and observed magnetic fields and angles. The angles, in turn, are used to estimate the spacecraft attitude and hence, the differences between model and observations as well as between two models become important to determine and reduce the errors associated with the models under different space environment conditions. We show that the models differ from the observations even during the geomagnetically quiet times but the associated errors during the geomagnetically active times increase. We find that the T89 model gives closer predictions to the observations, especially during active times and the errors are smaller compared to the IGRF-12 model. The magnitude of the error in the angle under both environmental conditions was found to be less than 1°. For the first time, the geomagnetic models were used to address the effects of the near Earth space environment on the satellite attitude.     

Variations of geomagnetic cutoff rigidity in the southern hemisphere close to 70°W (South-Atlantic Anomaly and Antarctic zones) in the period 1975–2010

We report the existence of rapid variations in (effective) geomagnetic cutoff rigidity (Rc) between the equatorial and Antarctic zones adjacent to the Andes Mountains, revealed by the variation rate of geomagnetic cutoff rigidity (VRc) in the period 1975–2010. Our analysis is based on empirical records and theoretical models of the variations in cosmic rays and on the structure of geomagnetic fields. These have given us a different view of variations in Rc in time and space along the 70°W meridian, where secular variations in the geomagnetic field are strongly influenced by the proximity of the South Atlantic Magnetic Anomaly (SAMA), one of the most important characteristics of the terrestrial magnetic field that affects our planet, close from the equator to the 50°S parallel and from South America to South Africa. The VRc presents rapid changes in mid-latitudes where SAMA exerts its influence despite the existence of smooth changes in the geomagnetic field. This shows that these changes occur mainly in the spatial configuration, rather than in the temporal evolution of Rc. The analysis was performed using measurements from the Chilean Network of Cosmic Rays and Geomagnetism Observatories, equipped with BF-3 and latest generation He-3 neutron monitors, Fluxgate magnetometers, geomagnetic reference field (IGRF) and Tsyganenko 2001 model (just for completeness).     

Substorm activity during the main phase of magnetic storms induced by the CIR and ICME events

In this work, the relation of high-latitude indices of geomagnetic activity (AE, Kp) with the rate of storm development and a solar wind electric field during the main phase of magnetic storm induced by the CIR and ICME events is investigated. 72 magnetic storms induced by CIR and ICME events have been selected. It is shown that for the CIR and ICME events the increase of average value of the Kp index (Kp_aver) is observed with the growth of rate of storm development. The value of Kp_aver index correlates with the magnitude of minimum value of Dst index (|Dst_min|) only for the ICME events. The analysis of average values of AE and Kp indices during the main phase of magnetic storm depending on the SW electric field has shown that for the CIR events, unlike the ICME events, the value of AE_aver increases with the growth of average value of the electric field (Esw_aver). The value of Kp_aver correlates with the Esw_aver only for the ICME events. The relation between geomagnetic indices and the maximum value of SW electric field (Esw_max) is weak. However, for the ICME events Kp_aver correlates with Esw_max.     

Interconnection of high-latitude and low-latitude boundary layers when IMF BY is dominant

The antiparallel merging model places the location of the reconnection region for a dominant interplanetary magnetic field (IMF) B_Y at high latitudes at the dayside magnetopause and predicts that the low-latitude boundary layer (LLBL) is located on open field lines of the magnetospheric flanks. Interball-1 data obtained in the wide local time range near the low-latitude magnetopause makes it possible to analyze the LLBL plasma population and to find a link between possible reconnection at high latitudes and LLBL occurrence. We found that no boundary layer was observed in the regions which have no topological connection with the merging site. All cases of LLBL observations are located downstream from a specific boundary. This boundary coincides with the first magnetospheric field line touching the reconnection region and can be located in a wide local time region depending on the instant IMF direction. Even the LLBL on closed field lines shows the tendency to be concentrated in the vicinity of this boundary. Thus we show that all types of observed LLBLs are linked to reconnection sites predicted by the antiparallel merging model.