1980．0年东亚地区地磁场的勒让德多项式模型

根据中国、原苏联等地区的地磁复测点和地磁台资料，利用勒让德多项式方法，计算出１９８０．０年东亚地区的地磁场模型，绘制出总强度（Ｆ）、磁偏角（Ｄ）和磁倾角（Ｉ）地磁图．地磁场模型的均方偏差分别为：总强度１４４．６ｎＴ，磁偏角１９．１′，磁倾角８．５′．将勒让德多项式模型与ＤＧＲＦ１９８０进行了比较，并绘制出总强度的差值分布图．本文计算出来的地磁场勒让德多项式模型能较好地表示东亚地区地磁场的分布．

东亚地磁场模型的计算与分析

根据中国、前苏联和蒙古等地区的地磁复测点和地磁台站资料 ,使用泰勒多项式方法和冠谐分析方法 ,计算东亚地磁场 (X ,Y ,Z)的泰勒多项式模型和冠谐模型 ,以及东亚剩余磁场 (ΔX ,ΔY ,ΔZ)的冠谐模型和泰勒多项式模型 ,并绘制了相应的理论地磁图 .泰勒多项式模型的展开原点位于 45°N和 1 0 0°E ,截断阶数为 7.地磁场泰勒多项式模型的均方偏差为 :X分量是1 3 3 0nT ,Y分量是 1 0 7 4nT ,Z分量是 1 48 0nT .球冠极点位于 45°N和 1 0 0°E ,球冠半角为 42°,冠谐模型的截断阶数为 1 0 .剩余磁场冠谐模型的均方偏差分别为 1 3 1 2nT(ΔX) ,1 1 2 6nT(ΔY)和1 3 8 7nT(ΔZ) .对比分析了上述两种模型 .提出了确定区域模型截断阶数的判据 .

Online computation of International Reference Ionosphere Extended to Plasmasphere(IRI-Plas) model for space weather

Ionosphere is the most challenging part of Space Weather with its spatio-temporal variability and dispersive nature. Ionospheric models are very important in reducing positioning error in GNSS system.International Reference Ionosphere(IRI) is an empirical, deterministic and climatic model of ionosphere up to 2000 km in height. Recently, IRI Extended to Plasmasphere(IRI-Plas) model has been developed to extend the interest region of IRI to the GPS orbital height of 20,000 km. Both IRI and IRI-Plas provide ionospheric parameters such as electron density, electron and ion temperatures according to their height profiles. In order to update the model to current ionospheric conditions, IRI-Plas can input F2 layer critical frequency(foF2), maximum ionization height(hmF2), and also Total Electron Content(TEC).Online IRI-Plas is developed for the ionospheric community to run multiple tasks at various locations,dates and times with optional foF2, hmF2 and TEC inputs in a user-friendly manner. In this paper, we are going to present the capabilities of the Online IRI-Plas service and provide some comparisons between IRI-Plas outputs and ionosonde measurements. The comparison between online IRI-Plas foF2 outputs and ionosonde foF2 measurements indicates that the model with TEC input can significantly improve the representation of the current ionospheric state, which is very successful especially in the geomagnetically disturbed days.

Distribution of Helical Properties of Solar Magnetic Fields

We summarize studies of helical properties of solar magnetic fields such as current helicity and twist of magnetic fields in solar active regions (ARs), that are observational tracers of the alpha-effect in the solar convective zone (SCZ). Information on their spatial distribution is obtained by analysis of systematic magnetographic observations of active regions taken at Huairou Solar Observing Station of National Astronomical Observatories of Chinese Academy of Sciences. The main property is that the tracers of the alpha-effect are antisymmetric about the solar equator. Identifying longitudinal migration of active regions with their individual rotation rates and taking into account the internal differential rotation law within the SCZ known from helioseismology, we deduce the distribution of the effect over depth. We have found evidence that the alpha-effect changes its value and sign near the bottom of the SCZ, and this is in accord with the theoretical studies and numerical simulations. We discuss other regularities which can be revealed by further analysis such as possible dependence on longitude, time, and magnetic field strength, etc.

Local sunspot oscillations and umbral dots

Data analysis of sunspot oscillations based on a 6-hr SDO run of an observation showed that low frequency （0.2 〈ω 〈 1 mHz） oscillations are locally similar to three and five minute oscillations. The oscillations in the sunspot are concentrated in cells of a few arcsec, each of which has its own oscillation spectrum. The analysis of two scenarios for sunspot oscillations leads to a conclusion that local sunspot oscillations occur due to a subphotospheric resonator for slow MHD waves. Empirical models of a sunspot atmosphere and the theory of slow waves in thin magnetic flux tubes are applied to modeling the subphotospheric resonator. The spectrum of local oscillations consists of a great number of lines. This kind of spectrum can occur only if the subphotospheric resonator is a magnetic tube with a rather weak magnetic field. Magnetic tubes of this sort are umbral dots that appear due to the convective tongues in monolithic sunspots. The interrelation of local oscillations with umbral dots and wavefronts of traveling waves in sunspots is discussed.

Recent results of zebra patterns in solar radio bursts

This review covers the most recent experimental results and theoretical research on zebra patterns （ZPs） in solar radio bursts. The basic attention is given to events with new peculiar elements of zebra patterns received over the last few years. All new properties are considered in light of both what was known earlier and new theoretical models. Large-scale ZPs consisting of small-scale fiber bursts could be explained by simultaneous inclusion of two mechanisms when whistler waves ＂high- light＂ the levels of double plasma resonance （DPR）. A unique fine structure was observed in the event on 2006 December 13： spikes in absorption formed dark ZP stripes against the absorptive type Ⅲ-like bursts. The spikes in absorption can appear in accordance with well known mechanisms of absorptive bursts. The additional injection of fast particles filled the loss-cone （breaking the loss-cone distribution）, and the generation of the continuum was quenched at these moments. The maximum absorptive effect occurs at the DPR levels. The parameters of millisecond spikes are determined by small dimensions of the particle beams and local scale heights in the radio source. Thus, the DPR model helps to understand severai aspects of unusual elements of ZPs. However, the simultaneous existence of several tens of the DPR levels in the corona is impossible for any realistic profile of the plasma density and magnetic field. Three new theories of ZPs are examined. The formation of eigenmodes of transparency and opac- ity during the propagation of radio waves through regular coronal inhomogeneities is the most natural and promising mechanism. Two other models （nonlinear periodic space - charge waves and scattering of fast protons on ion-sound harmonics） could happen in large radio bursts.

The Role of Nuclei-Nuclei Interactions in the Production of Gamma-ray Lines in Solar Flares

Dramatic extensions of experimental possibilities (spacecraft RHESSI, CORONAS-F and others) in solar gamma-ray astronomy call for urgent, detailed theoretical consideration of a set of physical problems of solar activity and solar-terrestrial relationships that earlier may have only been outlined. Here we undertake a theoretical analysis of issues related to the production of gamma-radiation in the processes of interactions of energetic (accelerated) heavy and middle nuclei with the nuclei of the solar atmosphere (the so-called i-j interactions). We also make an estimate of the contribution of these interactions to the formation of nuclear and isotopic abundances of the solar atmosphere in the range of light and rare elements. The analysis is carried out for solar flares in the wide range of their intensities. We compare our theoretical estimates with RHESSI observations for the flare of 2002 July 23. It was shown that the 24Mg gamma-ray emission in this event was produced by the newly generated Mg nuclei. With a high probability, the gamma-ray line emission of 28Si nuclei from this flare was generated by the same processes.

Changes in Environmental Parameters and Their Impact on Forest Growth in Northern Eurasia

We performed an empirical investigation of forest growth for two types of forests in northern Eurasia (larches and spruces) in order to show that the sensitivity of trees to the variable climate and geomagnetic field can be seen even under the large-scale average. The main purpose of this research was to model a forest growth rate V for each forest type on the basis of several environmental parameters influencing the tree growth in a high degree and to find the differences and similarities of the larches and spruces’ response to changing environment. We showed that V, which is related to the annual tree-ring width, could be derived from the Normalized Difference Vegetation Index (NDVI) data. Averaged yearly by species for 1982-2006, it displayed a long-term decrease (most likely related to the global climate change) as well as short-term variations with periods of 2.2, 4 and 8 years. A composite function method was used for modeling. We selected several tree growth drivers (the temperature, precipitation, insolation and the geomagnetic field intensity) that were highly correlated with V, and a function was modeled that described the behavior of V. The correlation coefficients between the derived function and experimental time series were 0.8 for larches and 0.9 for spruces. Compared with spruces, larches demonstrated higher thermo-sensitivity. A loss of tree sensitivity to temperature changes is puzzling for dendroclimatology, as a similar process might have occurred during previous periods of sharp global climate changes (as observed currently). Sensitivity of trees to geomagnetic field changes is confirmed both at long- and short-timescales. Spruces are found to be more magnetosensitive than larches.

Performance of GPS slant total electron content and IRI-Plas-STEC for days with ionospheric disturbance

Total Electron Content （TEC） is an important observable parameter of the ionosphere which forms the main source of error for space based navigation and positioning systems. Since the deployment of Global Navigation Satellite Systems （GNSS）, cost-effective esti- mation of TEC between the earth based receiver and Global Positioning System （GPS） sat- ellites became the major means of investigation of local and regional disturbance for earthquake precursor and augmentation system studies. International Reference Iono- sphere （IRI） extended to plasmasphere （IRI-Plas） is the most developed ionospheric and plasmaspheric climatic model that provides hourly, monthly median of electron density distribution globally. Recently, IONOLAB group {www.ionolab.org） has presented a new online space weather service that can compute slant TEC （STEC） on a desired ray path for a given date and time using IRI-Plas model （IRI-Plas-STEC）. In this study, the performance of the model based STEC is compared with GPS-STEC computed according to the estimation method developed by the IONOLAB group and includes the receiver bias as IONOLAB-BIAS （IONOLAB-STEC）. Using Symmetric Kullback-Leibler Distance （SKLD）, Cross Correlation （CC） coefficient and the metric norm （L2N） to compare IRI-Plas-STEC and IONOLAB-STEC for the month of October 2011 over the Turkish National Permanent GPS Network （TNPGN- Active）, it has been observed that SKLD provides a good indicator of disturbance for both earthquakes and geomagnetic storms.

Envelope Soliton in Solar Radio Emission

Several envelope soliton fine structures have been observed in solar radio metric-wave emission.present a model of longitudinal modulational instability to explain these fine structures.It is found that this instability can only occur in the condition of sound velocity being larger than Alfvin velocity in corona.Therefore,the envelope soliton fine structures should display in the coronal region with high temperature and low magnetic 6eld,which corresponds to the solar radio emission in the region of meter and decameter wavelength.

On the Positronium Continuum and 0.511 MeV Line in Solar Flares

We have studied the influence of the density of the annihilation region on the positronium continuum. A relation between the ratio 3γ/2γ and the density is explicitly given, with which one can derive directly from the observed 3γ/2γ the density where the annihilation occurs. A unique solution may be found from the observed width of the 0.511 MeV line. We applied the method to three flares observed by GRS/SMM. It is shown that due to the measuring uncertainties in the 0.511 MeV line width, we cannot distinguish a chromospheric source from a coronal source, though both accurately localized. To improve the measuring accuracy of the 0.511 MeV line and the ratio 3γ/2γ will be an important step for a better understanding of the annihilation process in solar flares.

Flare evolution and polarization changes in fine structures of solar radio emission in the 2013 April 11 event

The measurement of positions and sizes of radio sources in observations is important for un- derstanding of the flare evolution. For the first time, solar radio spectral fine structures in an M6.5 flare that occurred on 2013 April 11 were observed simultaneously by several radio instruments at four different observatories： Chinese Solar Broadband Radio Spectrometer at Huairou （SBRS/Huairou）, Ondrejov Radio Spectrograph in the Czech Republic （ORSC/Ondrejov）, Badary Broadband Microwave Spectropolarimeter （BMS/Irkutsk）, and spectrograph/IZMIRAN （Moscow, Troitsk）. The fine structures included microwave zebra patterns （ZPs）, fast pulsations and fiber bursts. They were observed during the flare brightening lo- cated at the tops of a loop arcade as shown in images taken by the extreme ultraviolet （EUV） telescope onboard NASA＇s satellite Solar Dynamics Observatory （SDO）. The flare occurred at 06：58-07：26 UT in solar active region NOAA 11719 located close to the solar disk center. ZPs appeared near high frequency boundaries of the pulsations, and their spectra observed in Huairou and Ondrejov agreed with each other in terms of details. At the beginning of the flare＇s impulsive phase, a strong narrowband ZP burst occurred with a moderate left-handed circular polarization. Then a series of pulsations and ZPs were observed in almost unpolarized emission. After 07：00 UT a ZP appeared with a moderate right-handed polarization. In the flare decay phase （at about 07：25 UT）, ZPs and fiber bursts become strongly right-hand polarized. BMS/Irkutsk spectral observations indicated that the background emission showed a left-handed circular polarization （similar to SBRS/Huairou spectra around 3 GHz）. However, the fine structure appeared in the right-handed polarization. The dynamics of the polarization was associated with the motion of the flare ex- citer, which was observed in EUV images at 171 A and 131 A by the SDO Atmospheric Imaging Assembly （AIA）. Combining magnetograms observed by the SDO Helioseismic and Magnetic Imager （HMI） with the homologous assumption of EUV flare brightenings and ZP bursts, we deduced that the observed ZPs correspond to the ordinary radio emission mode. However, future analysis needs to verify the assumption that zebra radio sources are really related to a closed magnetic loop, and are located at lower heights in the solar atmosphere than the source of pulsations.

The importance of source positions during radio fine structure observations

The measurement of positions and sizes of radio sources in the observations of the fine structure of solar radio bursts is a determining factor for the selection of the radio emission mechanism. The identical parameters describing the radio sources for zebra structures（ZSs） and fiber bursts confirm there is a common mechanism for both structures. It is very important to measure the size of the source in the corona to determine if it is distributed along the height or if it is point-like. In both models of ZSs（the double plasma resonance（DPR） and the whistler model） the source must be distributed along the height, but by contrast to the stationary source in the DPR model, in the whistler model the source should be moving. Moreover, the direction of the space drift of the radio source must correlate with the frequency drift of stripes in the dynamic spectrum. Some models of ZSs require a local source, for example,the models based on the Bernstein modes, or on explosive instability. The selection of the radio emission mechanism for fast broadband pulsations with millisecond duration also depends on the parameters of their radio sources.

Concerning spikes in emission and absorption in the microwave range

In some events, weak fast solar bursts （near the level of the quiet Sun） were observed in the background of numerous spikes in emission and absorption. In such a case, the background contains the noise signals of the receiver. In events on 2005 September 16 and 2002 April 14, the solar origin of fast bursts was confirmed by simultaneous recording of the bursts at several remote observatories. The noisy background pixels in emission and absorption can be excluded by subtracting a higher level of continuum when constructing the spectra. The wavelet spectrum, noisy pro- files in different polarization channels and a spectrum with continuum level greater than zero demonstrates the noisy character of pixels with the lowest levels of emission and absorption. Thus, in each case, in order to judge the solar origin of all spikes, it is necessary to determine the level of continuum against the background of which the solar bursts are observed. Several models of microwave spikes are discussed. The electron cyclotron maser emission mechanism runs into serious problems with the in- terpretation of microwave millisecond spikes： the main obstacles are too high values of the magnetic field strength in the source （Pe 〈 uB）. The probable mechanism is the interaction of plasma Langmuir waves with ion-sound waves （l ＋ s → t） in a source related to shock fronts in the reconnection region.

Proton activity of the Sun in current solar cycle 24

We present a study of seven large solar proton events in the current solar cycle 24（from 2009 January up to the current date）. They were recorded by the GOES spacecraft with the highest proton fluxes being over 200 pfu for energies 〉10 Me V. In situ particle measurements show that：（1） The profiles of the proton fluxes are highly dependent on the locations of their solar sources, namely flares or coronal mass ejections（CMEs）, which confirms the ＂heliolongitude rules＂ associated with solar energetic particle fluxes;（2） The solar particle release（SPR） times fall in the decay phase of the flare emission, and are in accordance with the times when the CMEs travel to an average height of 7.9 solar radii; and（3） The time differences between the SPR and the flare peak are also dependent on the locations of the solar active regions. The results tend to support the scenario of proton acceleration by the CME-driven shock,even though there exists a possibility of particle acceleration at the flare site, with subsequent perpendicular diffusion of accelerated particles in the interplanetary magnetic field. We derive the integral time-of-maximum spectra of solar protons in two forms： a single power-law distribution and a power law roll-over with an exponential tail. It is found that the unique ground level enhancement that occurred in the event on 2012 May 17 displays the hardest spectrum and the largest roll-over energy which may explain why this event could extend to relativistic energies.

Distribution of apparent magnetization for Asia

Magsat total field anomalies over Asia were used to construct an equivalent magnetization model, which represents the apparent magnetization distribution within an equivalent layer 40 km thick and correlates well with large-scale tectonics, for example, the Kazakhstan, Tarim, Yangtze, India, Sino-Korea and Indochina blocks. The basin, plain, sea basin, and islands are delineated by magnetization lows whereas the plateau and marine ridge correspond to magnetization highs. The boundary between Tibetan Plateau and India marked by a strong gradient along its length coincides with the Yarlung Zangbo River fault roughly. The Tanlu fault belt is the boundary between positive and negative anomalies. This boundary stretches in southwest direction and joins Sanjiang fault belt. The boundary between the Southeast China block and the Yangtze block is also clearly delineated by the magnetization anomalies. Generally, the magnetization boundaries are consistent with the collisional suture of blocks.