质子地磁垂直截止刚度在2004年11月7-8日两个时刻处的数值模拟

利用单粒子轨道理论, 在T04和IGRF2000磁场模式建立的磁层模型基础上, 应用四阶龙格库塔方法, 模拟计算宇宙线带电质子在地球磁层中的运动以及沿天顶方向入射到达地球磁场内的某一特定位置, 得到了距地球表面450km高度处全球质子垂直截止刚度在2004年11月7-8日中两个时刻的计算值. 根据计算得到的684个不同位置处的截止刚度值, 分析了同一时刻地磁垂直截止刚度随磁纬和磁经的变化. 与此同时, 模拟计算了相同时刻下磁纬为30°, 磁经为0°, 45°, 90°, 135°, 180°, 225°, 270°, 315°处, 地磁有效垂直截止刚度随高度的变化情况. 结果表明, 在磁纬30°处, 其质子有效垂直截止刚度随距离地心高度的变化沿向阳处向磁尾处方向, 变化越来越缓慢.      

Latitudinal features of Total Electron Content over the African and European longitude sector following the St. Patrick’s day storm of 2015

GNSS TEC values have been obtained from 18 stations distributed from the magnetic equator to nearly 80°N magnetic dip in the African and west-European longitude sector corresponding to the March 17–18, 2015 geomagnetic storm. Significantly depleted ionosphere have been observed at stations north of 50°N geographic on March 18, 2015 following the above storm over a longitude swath 11.9°–21°E covering the Eastern Africa and Western European longitude sector. High ROTI values were noted on March 17^th at locations around 80°N magnetic dip. Two prominent peaks in PCN were noted around 09:00 UT and 14:00 UT on March 17, 2015 and around 15:00 UT on March 18, 2015. Daytime thermospheric (O/N₂) ratio was markedly less on March 18th at latitudes above 60°N geographic which is suggested to be the major driver behind depleted high latitude ionosphere during the recovery phase of the storm on March 18, 2015.     

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.     

The local time dependence of the anisotropic solar cosmic ray flux.

The distribution of the solar cosmic radiation flux over the earth is not uniform, but the result of complex phenomena involving the interplanetary magnetic field, the geomagnetic field and latitude and longitude of locations on the earth. The latitude effect relates to the geomagnetic shield; the longitude effect relates to local time. For anisotropic solar cosmic ray events the maximum particle flux is always along the interplanetary magnetic field direction, sometimes called the Archimedean spiral path from the sun to the earth. During anisotropic solar cosmic ray event, the locations on the earth viewing "sunward" into the interplanetary magnetic field direction will observe the largest flux (when adjustments are made for the magnetic latitude effect). To relate this phenomena to aircraft routes, for anisotropic solar cosmic ray events that occur during "normal quiescent" conditions, the maximum solar cosmic ray flux (and corresponding solar particle radiation dose) will be observed in the dawn quadrant, ideally at about 06 hours local time.     

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).     

Responses of the African equatorial ionization anomaly (EIA) to some selected intense geomagnetic storms during the maximum phase of solar cycle 24

This study investigates the morphology of the GPS TEC responses in the African Equatorial Ionization Anomaly (EIA) region to intense geomagnetic storms during the ascending and maximum phases of solar cycle 24 (2012–2014). Specifically, eight intense geomagnetic storms with Dst ≤ ?100 nT were considered in this investigation using TEC data obtained from 13 GNSS receivers in the East African region within 36–42°E geographic longitude; 29°N–10°S geographic latitude; ± 20°N magnetic latitude. The storm-time behavior of TEC shows clear positive and negative phases relative to the non-storm (median) behavior, with amplitudes being dependent on the time of sudden commencement of the storm and location. When a storm starts in the morning period, total electron content increases for all stations while a decrease in total electron content is manifested for a storm that had its sudden commencement in the afternoon period. The TEC and the EIA crest during the main phase of the storm is significantly impacted by the geomagnetic storm, which experiences an increase in the intensity of TEC while the location and spread of the crest usually manifest a poleward expansion.     

Effects of CME and CIR induced geomagnetic storms on low-latitude ionization over Indian longitudes in terms of neutral dynamics

This paper presents the response of the ionosphere during the intense geomagnetic storms of October 12–20, 2016 and May 26–31, 2017 which occurred during the declining phase of the solar cycle 24. Total Electron Content (TEC) from GPS measured at Indore, Calcutta and Siliguri having geomagnetic dips varying from 32.23°N, 32°N and 39.49°N respectively and at the International GNSS Service (IGS) stations at Lucknow (beyond anomaly crest), Hyderabad (between geomagnetic equator and northern crest of EIA) and Bangalore (near magnetic equator) in the Indian longitude zone have been used for the storms. Prominent peaks in diurnal maximum in excess of 20–45 TECU over the quiet time values were observed during the October 2016 storm at Lucknow, Indore, Hyderabad, Bangalore and 10–20 TECU for the May 2017 storm at Siliguri, Indore, Calcutta and Hyderabad. The GUVI images onboard TIMED spacecraft that measures the thermospheric O/N₂ ratio, showed high values (O/N₂ ratio of about 0.7) on October 16 when positive storm effects were observed compared to the other days during the storm period. The observed features have been explained in terms of the O/N₂ ratio increase in the equatorial thermosphere, CIR-induced High Speed Solar Wind (HSSW) event for the October 2016 storm. The TEC enhancement has also been explained in terms of the Auroral Electrojet (AE), neutral wind values obtained from the Horizontal Wind Model (HWM14) and equatorial electrojet strength from magnetometer data for both October 2016 and May 2017 storms. These results are one of the first to be reported from the Indian longitude sector on influence of CME- and CIR-driven geomagnetic storms on TEC during the declining phase of solar cycle 24.     

Statistical average estimates of high latitude field-aligned currents from the STARE and SABRE coherent VHF radar systems

Two bistatic VHF radar systems, STARE and SABRE, have been employed to estimate ionospheric electric fields in the geomagnetic latitude range 61.1 – 69.3° (geographic latitude range 63.8 – 72.6°) over northern Scandinavia. 173 days of good backscatter from all four radars have been analysed during the period 1982 to 1986, from which the average ionospheric divergence electric field versus latitude and time is calculated. The average magnetic field-aligned currents are computed using an AE-dependent empirical model of the ionospheric conductance. Statistical Birkeland current estimates are presented for high and low values of the Kp and AE indices as well as positive and negative orientations of the IMF B_z component. The results compare very favourably to other ground-based and satellite measurements.     

大气中子实时计算模式

为实时评估0~100km高度范围内的大气中子全球分布,对宇宙线在地磁场和大气中的传输过程进行了分析.利用蒙特卡罗方法工具包Geant4,预先计算不同能量的粒子在大气层中产生的次级粒子能谱分布,形成大气次级粒子数据库,并与相关模型进行对比,验证了该数据库的有效性和可靠性.以实测或预报的空间环境参数作为输入,计算同步轨道银河宇宙线和太阳质子事件能谱以及100km高度上的地磁垂直截止刚度,最终得到大气层顶上的粒子能谱.通过对大气次级粒子数据库的线性插值,实现1h分辨率的大气中子能谱和辐射剂量全球分布的实时计算.      

A method for operative calculation of charged particle penetration into the LEO

We develop a simple method for calculating the effective vertical cutoff rigidity of charged particles, taking into account the Kp-index and the local time, on the basis of generalization of the results of extensive trajectory calculations for trial particles moving in the geomagnetic field. The vertical cutoff rigidities, calculated by the Tsyganenko-89 model, are presented as an International Geomagnetic Reference Field (IGRF) model calculated and thereafter corrected in accordance with the geomagnetic disturbance and local time conditions. The fits from the proposed method agree with the results of cutoff rigidity measurements carried out by satellites. The method is intended for applications using cutoff calculations, such as evaluating particle penetration of spatial boundaries, calculating magnetospheric transmissions for low-orbital spacecrafts flights and interpreting the results of orbital experiments.     

一种基于GNSS数据的全球电离层不规则结构活动表征方法

针对如何利用GNSS(Global Navigation Satellite System)数据进行电离层扰动监测的问题,提出了一种基于GNSS数据表征全球电离层扰动的方法.利用大约400个GNSS地面站点的观测数据,计算总电子含量(Total Electron Content,TEC)变化率的标准差——ROTI(Ra...     

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.     

The regional peculiarities of the total ozone content variations caused by solar/geomagnetic phenomena

The ozone variations possibly caused by solar electromagnetic radiation, geomagnetic storms and solar particle events depend on the latitude and longitude. The results of the statistical analysis on the base of TOMS total ozone content (TOC) measurements are compared for the regions with the same geographical or geomagnetic latitude but with different stratospheric and/or tropospheric dynamics. The atmospheric circulation could be the intermediate link of a chain of solar/geomagnetic influence on the TOC.     

GPS derived ionospheric TEC response to geomagnetic storm on 24 August 2005 at Indian low latitude stations

Results pertaining to the response of the low latitude ionosphere to a major geomagnetic storm that occurred on 24 August 2005 are presented. The dual frequency GPS data have been analyzed to retrieve vertical total electron content at two Indian low latitude stations (IGS stations) Hyderabad (Geographic latitude 17°20′N, Geographic longitude 78°30′E, Geomagnetic latitude 8.65°N) and Bangalore (Geographic latitude 12°58′N, Geographic longitude 77°33′E, Geomagnetic latitude 4.58°N). These results show variation of GPS derived total electron content (TEC) due to geomagnetic storm effect, local low latitude electrodynamics response to penetration of high latitude convection electric field and effect of modified fountain effect on GPS–TEC in low latitude zone.     

Calculations and observations of solar particle enhancements to the radiation environment at aircraft altitudes.

Solar particle events can give greatly enhanced radiation at aircraft altitudes, but are both difficult to predict and to calculate retrospectively. This enhanced radiation can give significant dose to aircrew and greatly increase the rate of single event effects in avionics. Validation of calculations is required but only very few events have been measured in flight. The CREAM detector on Concorde detected the event of 29 September 1989 and also four periods of enhancement during the events of 19-24 October 1989. Instantaneous rates were enhanced by up to a factor ten compared with quiet-time cosmic rays, while flight-averages were enhanced by up to a factor six. Calculations are described for increases in radiation at aircraft altitudes using solar particle spectra in conjunction with Monte Carlo radiation transport codes. In order to obtain solar particle spectra with sufficient accuracy over the required energy range it is necessary to combine space data with measurements from a wide range of geomagnetically dispersed, ground-level neutron monitors. Such spectra have been obtained for 29 September 1989 and 24 October 1989 and these are used to calculate enhancements that are compared with the data from CREAM on Concorde. The effect of cut-off rigidity suppression by geomagnetic activity is shown to be significant. For the largest event on record on 23 February 1956, there are no space data but there are data from a number of ground-level cosmic-ray detectors. Predictions for all events show very steep dependencies on both latitude and altitude. At high latitude and altitude (17 km) calculated increases with respect to cosmic rays are a factor 70 and 500 respectively for 29 September 1989 and 23 February 1956. The levels of radiation for high latitude, subsonic routes are calculated, using London to Los Angeles as an example, and can exceed 1 mSv, which is significantly higher than for Concorde routes from Europe to New York. The sensitivity of the calculations to spectral fitting, geomagnetic activity and other assumptions demonstrates the requirement for widespread carriage of radiation monitors on aircraft.     

HXMT卫星空间环境监测器及初步观测结果

硬X射线调制望远镜（HXMT）卫星是中国首个专门进行天文探测的空间科学实验卫星,运行于高度约550km、倾角约43°的低地球轨道.星载空间环境监测器为星上科学任务开展提供背景辐射实测资料.该监测器采用固体探测器望远镜系统和扇形阵列全新组合设计,可获取轨道空间高能质子和高能电子能谱、方向综合动态结果,给出更为全面的粒子辐射分布图像.初步探测结果显示,卫星运行轨道遭遇的带电粒子辐射集中分布在经度80°W-20°E,纬度0°-40°S的南大西洋异常区,粒子辐射在该区域表现出不同程度的方向差异分布,高能电子方向差异分布显著强于高能质子.2017年9月空间环境扰动期间,爆发的太阳质子事件并未对该轨道粒子辐射产生影响,而地磁活动导致该轨道穿越经度120°W-60°E,纬度40°-43°N的北美上空和经度60°-120°E,纬度43°-40°S的澳大利亚西南区域时遭遇增强粒子辐射影响,增强的粒子辐射表现出极强的方向分布.      

Observation of the north-south asymmetry in the equatorial spread F intensity at the Brazilian longitude

This study presents an analysis of the observed north-south asymmetry of the range spread F (RSF) intensity at the low latitude region during an equinoctial month of different solar epochs (2002, 2015 and 2017). The ionospheric parameters were obtained during geomagnetic quiet days from four digisonde stations located along the Brazilian longitude, which include a dip equator station (Sao Luiz (SL: 2.33 S, 44.2 W)), conjugate stations (Campo Grande (CG: 20.5°S, 55°W) and Boa Vista (BV: 2.8°N, 60.7°W)) and another low latitude station (Cachoeira Paulista (CP: 22.7°S, 45°W)). The results highlight the competing effect of the post-sunset electric field strength and the trans-equatorial wind on the latitudinal distribution of the irregularity intensity at both hemispheres under varying background ionospheric condition. The RSF intensity was seen to reduce as the solar flux index decreased and the latitudinal peak shifted closer to the dip equator. This was dependent on the variation of the field line mapped irregularity spectrum and the density gradient. Likewise, the north-south asymmetry in the irregularity occurrence was seen to become more significant as a denser ionosphere was observed at the hemisphere with the equatorward meridional wind. This has further proven that the non-linear cascading of the plasma irregularity across the low latitude region is strongly influenced by the local electric field.     

Modelling bottom and topside electron density and TEC with profile data from topside ionograms

The paper describes the technique that has been implemented to model the electron density distribution above and below the F2 peak making use of only the profiles obtained from the INTERCOSMOS-19 topside ionograms. Each single profile from the satellite height to the ionosphere peak has been fitted by a semi-Epstein layer function of the type used in the DGR model with shape factor variable with altitude. The topside above the satellite height has been extrapolated to match given values of plasmaspheric electron densities to obtain the full topside profile. The bottomside electron density has been calculated by using the maximum electron density and its altitude estimated from the topside ionogram as input for a modified version of the DGR derived profiler that uses model values for the foF1 and foE layers of the ionosphere. Total electron content has also been calculated. Longitudinal cross sections of vertical profiles from latitudes 50° N to 50° S latitude are shown for low and high geomagnetic activity. These cross sections indicate the equatorial anomaly effect and the changes of the shape of low latitude topside ionosphere during geomagnetic active periods. These results and the potentiality of the technique are discussed.     

A new dynamical atmospheric ionizing radiation (AIR) model for epidemiological studies.

A new Atmospheric Ionizing Radiation (AIR) model is currently being developed for use in radiation dose evaluation in epidemiological studies targeted to atmospheric flight personnel such as civilian airlines crewmembers. The model will allow computing values for biologically relevant parameters, e.g. dose equivalent and effective dose, for individual flights from 1945. Each flight is described by its actual three dimensional flight profile, i.e. geographic coordinates and altitudes varying with time. Solar modulated primary particles are filtered with a new analytical fully angular dependent geomagnetic cut off rigidity model, as a function of latitude, longitude, arrival direction, altitude and time. The particle transport results have been obtained with a technique based on the three-dimensional Monte Carlo transport code FLUKA, with a special procedure to deal with HZE particles. Particle fluxes are transformed into dose-related quantities and then integrated all along the flight path to obtain the overall flight dose. Preliminary validations of the particle transport technique using data from the AIR Project ER-2 flight campaign of measurements are encouraging. Future efforts will deal with modeling of the effects of the aircraft structure as well as inclusion of solar particle events.