Is there an abnormal enhancement of atmospheric aerosol before the 2008 Wenchuan earthquake?

There are extensive reports of ionospheric disturbances before the great 2008 Wenchuan earthquake, which are possibly explained by seismogenic electric field hypotheses linked with the aerosols injected in atmosphere. This paper attempts to investigate the possible change of atmospheric aerosol optical depth (AOD) associated with this earthquake by using MODIS data from both Terra and Aqua satellites. The result shows a clear enhancement of AOD along the Longmenshan faults 7 days before the quake, which is 1 day and 4 days earlier than the reported negative and positive ionospheric disturbances, respectively, and is 1 day earlier than or quasi-synchronism with other reported atmospheric anomalies including air temperature, outgoing longwave radiation and relative humidity. Particularly, the spatial distribution of AOD enhancement is very local and it is correlated well with the active faults and surface ruptures. We suggest that this unique enhancement could be associated with the Lithosphere–Atmosphere–Ionosphere coupling process during the preparation of the Wenchuan earthquake.     

The propagation features of LF radio waves at topside ionosphere and their variations possibly related to Wenchuan earthquake in 2008

The artificial low frequency (LF) signals recorded at topside ionosphere from DEMETER satellite were analyzed in this paper, and the typical diurnal and seasonal variations were illustrated around 162?kHz in electric field spectra. The larger power spectrum density (PSD) values in electric field at local nighttime and in winter season all demonstrate the correlation feature of lower plasma content to higher penetration of LF waves into ionosphere. Around Wenchuan earthquake, the comparison of signal-noise-ratio (SNR) values in electric field with each half month during January to May in 2008 and the same half month in May from 2005 to 2007 revealed their lowest values and small covering area around the preparation region of Wenchuan earthquake in 2008. Combined with other researches in VLF radio waves and geochemical observations from satellite, the interaction of ion accumulation and upward movement from gas-water release at surface might be a key factor to disturb the ionospheric plasma density, and then possibly leading to the decrease of low energy penetration of LF radio waves from the artificial transmitted source at ground.     

Anomalous enhancement of electric field derived from ionosonde data before the great Wenchuan earthquake

The seismo-ionospheric coupling is a hot topic of modern scientific research. One of the most reasonable mechanisms is the hypothesis of quasistatic electric field. With a preliminary analytical method, electric field in F2-region before Wenchuan earthquake was derived from the five lower latitude ground-based ionosondes: Lhasa (29.63°N, 91.17°E), Chongqing (29.50°N, 106.40°E), Kunming (25.00°N, 102.70°E), Guangzhou (23.15°N, 113.35°E) and Haikou (23.15°N, 110.35°E). The results show that there was an anomalous enhancement of electric field in F2-region close to the epicenter with maximal amplitude ∼2 mV/m (about 10 times of the background) at 15:00 LT on 9 May, which is in line with the amplitudes of anomalous electric field related to some previous earthquakes observed by satellites.     

Formation mechanism of great positive TEC disturbances prior to Wenchuan earthquake on May 12, 2008

The problems of physical explanation and possible mechanisms of the seismo-ionospheric effects formation are under discussion now. There are proposed different mechanisms of such effects, for example, large- and small-scale internal gravity waves (IGWs), atmospheric electric field, electromagnetic fields and emissions. However, the appearance of local large-scale seismo-ionospheric anomalies in Total Electron Content (TEC) is possible to explain only by two mechanisms: an atmospheric electric field and/or small-scale IGWs. In this paper, the simulation results for reproduction of the observed seismo-ionospheric great positive effects in TEC prior to strong Wenchuan earthquake are presented. The obtained results confirm the proposed mechanism of seismo-ionospheric effects formation by the penetration of the seismogenic electric field from the atmosphere into the ionosphere. It is suggested that so great TEC enhancement observed 3 days prior to Wenchuan earthquake could be explained by combined action of seismogenic vertical electric field and IGWs generated by the solar terminator.     

Satellite and ground detection of very dense smoke clouds produced on the islands of the Paraná river delta that affected a large region in Central Argentina

Intense fires were produced on the Paraná river delta islands, Argentina, during most part of 2008, by a combination of an exceptionally dry period and the farmers’ use of a fire land-cleaning technique. In April 2008, those fires significantly affected the nearby regions and their inhabitants, from Rosario city to Buenos Aires mega-city. In this work we present satellite as well as ground Aerosol Optical Depth (AOD) at 550 nm data obtained during the propagation of pollution clouds to the central zone of Argentina. The highest value (1.18) was registered at Buenos Aires by atmospheric remote sensing, using the satellite instrument MODIS/Terra on April 18th 2008 at 10:35 local time (= UT − 3 h). On the same day, ground air quality detectors also measured in this city the highest Total Suspended Particle (TSP) value of the month, 2.02 mg/m³. The AOD(550) daily variation at Rosario Astronomical Observatory, which is located near the Paraná riverside, was derived by combining solar ultraviolet erythemal irradiance data (measured with a YES biometre) with model calculations. On April 25th 2008, from 12:00 to 15:30 local time, a rather high and constant AOD(550) value was registered, with a mean value of (0.90 ± 0.21). Cities located on the side of the Rosario–Buenos Aires highway (San Nicolás, Baradero and San Pedro) were also affected, showing a mean AOD(550) between the Rosario and Buenos Aires values. The particulate matter was collected with gridded samplers placed on the Paraná river islands as well as at the Rosario Observatory. They were analysed with a Scanning Electron Microscope (SEM) and mainly showed a biological origin. Even if normally large particles travel small distances from the source, organic aerosol in the range of 40–100 μm and complex asymmetric structures were registered several kilometres away from the aerosol sources on the islands. Another event of intense UV index attenuation (98.6%) occurred on September 18th 2008, due to very dense smoke clouds that extended over the Rosario area for several hours. The clouds were driven away from the fires by East–northeast and East–southeast winds. The minimum value of this index measured around noon allows to derive a maximum AOD(550)_max = (3.65 ± 0.90) at 12:45 local time. Soot clouds extended over the Paraná river, transporting Burned Biomass Debris (BBD) that deposited on Rosario. In particular, burned leaves and small branches with dimensions of 1–20 cm were collected. The mean (BBD) particles deposited on the ground from 7:00 to 19:00 local time were (0.92 ± 0.20) BBD/(m² h).     

Use of principal component analysis in the identification of the spatial pattern of an ionospheric total electron content anomalies after China’s May 12, 2008, M = 7.9 Wenchuan earthquake

This paper uses principal component analysis (PCA) to determine the spatial pattern of total electron content (TEC) anomalies in the ionosphere post the China’s Wenchuan Earthquake of 12 May, 2008 (UTC) (M_w = 7.9). PCA and image processing are applied to the global ionospheric map (GIM) with transforms conducted for the time period 08:00–10:00 UT on 12 May 2008. Results show that at a height of approximately 200 km the anomaly is widespread and less intense; however, it becomes more localized with height reaching maximum intensity and localization at an altitude of 300 km. The spatial distribution is remarkably similar to that reported for a TEC anomaly previously identified as a precursor anomaly on May 9, 2008 for the same time period. Potential causes of the results are discussed with emphasis given to vertical acoustic gravity waves based on the spatial pattern identified.     

Electric field and ion density anomalies in the mid latitude ionosphere: Possible connection with earthquakes?

The problem of earthquake prediction has stimulated the search for a correlation between seismic activity and ionospherical anomalies. We found observational evidence of possible earthquake effects in the near-equatorial and low latitude ionosphere; these ionospheric anomalies have been proposed by Gousheva et al. [Gousheva, M., Glavcheva, R., Danov, D., Angelov P., Hristov, P., Influence of earthquakes on the electric field disturbances in the ionosphere on board of the Intercosmos-Bulgaria-1300 satellite. Compt. Rend. Acad. Bulg. Sci. 58 (8) 911–916, 2005a; Gousheva, M., Glavcheva, R., Danov, D., Angelov, P., Hristov, P., Kirov, B., Georgieva, K., Observation from the Intercosmos-Bulgaria-1300 satellite of anomalies associated with seismic activity. In: Poster Proceeding of 2nd International Conference on Recent Advances in Space Technologies: Space in the Service of Society, RAST ‘2005, June 9–11, Istanbul, Turkey, pp. 119–123, 2005b; Gousheva, M., Glavcheva, R., Danov, D., Angelov, P., Hristov, P., Kirov, B., Georgieva, K., Satellite monitoring of anomalous effects in the ionosphere probably related to strong earthquakes. Adv. Space Res. 37 (4), 660–665, 2006]. This paper presents new results from observations of the quasi-static electric field and ion density on board INTERCOSMOS-BULGARIA-1300 satellite in the mid latitude ionosphere above sources of moderate earthquakes. Data from INTERCOSMOS-BULGARIA-1300 satellite and seismic data (World Data Center, Denver, Colorado, USA) for magnetically quiet and medium quiet days are juxtaposed in time-space domain. For satellite’s orbits in the time period 15.09–01.10.1981 an increase in the horizontal and vertical components of the quasi-static electric field and fluctuations of the ion density are observed over zones of forthcoming seismic events. Some similar post effects are observed too. The emphasis of this paper is put on the anomalies which specify the mid latitude ionosphere. The obtained results contain important information because they confirm our previous results for near-equatorial and low latitude regions.     

The asymmetrical features in electron density during extreme solar minimum

The variations of plasma density in topside ionosphere during 23rd/24th solar cycle minimum attract more attentions in recently years. In this analysis, we use the data of electron density (Ne) from DEMETER (Detection of Electromagnetic Emissions Transmitted from Earthquake Regions) satellite at the altitude of 660–710 km to investigate the solstitial and equinoctial asymmetry under geomagnetic coordinate system at LT (local time) 1030 and 2230 during 2005–2010, especially in solar minimum years of 2008–2009. The results reveal that ΔNe (December–June) is always positive over Southern Hemisphere and negative over northern part whatever at LT 1030 or 2230, only at 0–10°N the winter anomaly occurs with ΔNe (December–June) > 0, and its amplitude becomes smaller with the declining of solar flux from 2005 to 2009. The ΔNe between September and March is completely negative during 2005–2008, but in 2009, it turns to be positive at latitudes of 20°S–40°N at LT 1030 and 10°S–20°N at LT 2230. Furthermore, the solstitial and equinoctial asymmetry index (AI) are calculated and studied respectively, which all depends on local time, latitude and longitude. The notable differences occur at higher latitudes in solar minimum year of 2009 with those in 2005–2008. The equinoctial AI at LT 2230 is quite consistent with the variational trend of solar flux with the lowest absolute AI occurring in 2009, the extreme solar minimum, but the solstitial AI exhibits abnormal enhancement during 2008 and 2009 with bigger AI than those in 2005–2007. Compared with the neutral compositions at 500 km altitude, it illustrates that [O/N₂] and [O] play some roles in daytime and nighttime asymmetry of Ne at topside ionosphere.     

A nonlinear background removal method for seismo-ionospheric anomaly analysis under a complex solar activity scenario: A case study of the M9.0 Tohoku earthquake

A precise determination of ionospheric total electron content (TEC) anomaly variations that are likely associated with large earthquakes as observed by global positioning system (GPS) requires the elimination of the ionospheric effect from irregular solar electromagnetic radiation. In particular, revealing the seismo-ionospheric anomalies when earthquakes occurred during periods of high solar activity is of utmost importance. To overcome this constraint, a multiresolution time series processing technique based on wavelet transform applicable to global ionosphere map (GIM) TEC data was used to remove the nonlinear effect from solar radiation for the earthquake that struck Tohoku, Japan, on 11 March, 2011. As a result, it was found that the extracted TEC have a good correlation with the measured solar extreme ultraviolet flux in 26–34 nm (EUV_26–34) and the 10.7 cm solar radio flux (F_10.7). After removing the influence of solar radiation origin in GIM TEC, the analysis results show that the TEC around the forthcoming epicenter and its conjugate were significantly enhanced in the afternoon period of 8 March 2011, 3 days before the earthquake. The spatial distributions of the TEC anomalous and extreme enhancements indicate that the earthquake preparation process had brought with a TEC anomaly area of size approximately 1650 and 5700 km in the latitudinal and longitudinal directions, respectively.     

Seismo-ionospheric anomalies in total electron content of the GIM and electron density of DEMETER before the 27 February 2010 M8.8 Chile earthquake

In this paper we examine the pre-earthquake ionospheric anomalies by the total electron content (TEC) extracted from GIM (global ionospheric map) and the electron density (Ne) observed by the DEMETER (Detection of Electro-Magnetic Emissions Transmitted from Earthquake Regions) satellite during the 2010 M8.8 Chile earthquake. Temporal variations show the nighttime TEC and Ne simultaneously increase 9–19 days before the earthquake. A cross-comparison of data recorded during the period of 1 February to 3 March in 2006–2010 confirms the above temporal anomalies specifically appear in 2010. The spatial analyses show that the anomalies tend to appear over the epicenter.     

Surface latent heat flux anomalies quasi-synchronous with ionospheric disturbances before the 2007 Pu’er earthquake in China

One of various mechanisms of pre-earthquake lithosphere–atmosphere–ionosphere coupling as possible explanation of the seismo-ionospheric effects is connected with the release of latent heat. Abnormal variations of ionospheric electromagnetic parameters possibly related to the 2007 Ms 6.4 Pu’er earthquake in China were reported. This paper attempts to examine whether there were abnormal changes of surface latent heat flux (SLHF) linked with this pre-earthquake ionospheric disturbances. The spatio-temporal statistical analyzes of multi-years SLHF data from USA NCEP/NCAR Reanalysis Project reveal that local SLHF enhancements appeared 11, 10 and 7 days before the Pu’er earthquake, respectively. As contrasted to the formerly reported local ionospheric Ne enhancement 9 days before the shock observed by DEMETER satellite, it is discovered that the SLHF and Ne anomalies are quasi-synchronous and have good spatial correspondence with the epicentre and the local active faults. This is valuable for understanding the seismogenic coupling processes and for recognizing earthquake anomaly with multiple parameters from integrated Earth observation system.     

Magnetospheric electric field variations caused by storm-time shock fronts

On January 20, 2005 there was an X 7.1 solar flare at 0636 UT with an accompanied halo coronal mass ejection (CME). The resultant interplanetary shock impacted earth ∼36 h later. Near earth, the Advanced Composition Explorer (ACE) spacecraft observed two impulses with a staircase structure in density and pressure. The estimated earth-arrival times of these impulses were 1713 UT and 1845 UT on January 21, 2005. Three MINIature Spectrometer (MINIS) balloons were aloft on January 21st; one in the northern polar stratosphere and two in the southern polar stratosphere. MeV relativistic electron precipitation (REP) observed by all three balloons is coincident (<3 min) with the impulse arrivals and magnetospheric compression observed by both GOES 10 and 12. Balloon electric field data from the southern hemisphere show no signs of the impulse electric field directly reaching the ionosphere. Enhancement of the balloon-observed convection electric field by as much as 40 mV/m in less than 20 min during this time period is consistent with typical substorm growth. Precipitation-induced ionospheric conductivity enhancements are suggested to be (a) the result of both shock arrival and substorm activity and (b) the cause of rapid (<6 min) decreases in the observed electric field (by as much as 40 mV/m). There is poor agreement between peak cross polar cap potential in the northern hemisphere calculated from Super Dual Auroral Radar Network (SuperDARN) echoes and horizontal electric field at the MINIS balloon locations in the southern hemisphere. Possible reasons for this poor agreement include (a) a true lack of north–south conjugacy between measurement sites, (b) an invalid comparison between global (SuperDARN radar) and local (MINIS balloon) measurements and/or (c) radar absorption resulting from precipitation-induced D-region ionosphere density enhancements.     

Study of the 2013 Lushan M7.0 earthquake coseismic ionospheric disturbances

On April 20, 2013, an earthquake of M7.0 occurred in Lushan, Sichuan province, China. This paper investigates the coseismic ionospheric anomalies using GPS (Global Positioning System) data from 23 reference stations in Sichuan province that are a part of the Crustal Movement Observation Network of China (CMONOC). The recorded results show that a clear ionospheric anomaly occurred within 15 min after the earthquake near the epicenter, and the occurrence time of the anomalies recorded by various stations is related to the distance from the epicenter. The maximum anomaly is 0.25 TECu, with a 2 min duration and the distance of the recording station to the epicenter is 83 km. Acoustic waves generated by the crustal vertical movement during the earthquake propagate up to the height of the ionosphere lead to the ionospheric anomaly, and the propagation speed of the acoustic wave is calculated as 0.72 ± 0.04 km/s based on the propagation time and propagation distance, consistent with the average speed of sound waves within a 0–450 km atmospheric height.     

Variations of ionospheric total electron content before three strong earthquakes in the Qinghai-Tibet region

We investigate the ionospheric total electron content (TEC) anomalies occurred in the Qinghai-Tibet region before three large earthquakes (M > 7.0). The temporal and spatial TEC variations were used to detect the ionospheric possible precursors of these earthquakes. We identified two TEC enhancements in the afternoon local time 9 days and 2–3 days before each earthquake, between which a TEC decrement occurred 3–6 days before earthquakes. These anomalies happened in the area of about 30° in latitude and the maximum is localized equatorward from the epicenters. These TEC anomalies can be found in all three earthquakes regardless the geomagnetic conditions. The features of these anomalies have something in common and may have differences from those caused by geomagnetic storms. Our results suggest that these ionospheric TEC perturbations may be precursors of the large earthquakes.     

Time Transfer by Laser Link: Data analysis and validation to the ps level

The Time Transfer by Laser Link (T2L2) is a very high resolution time transfer technique based on the recording of arrival times of laser pulses at the satellite. T2L2 was designed to achieve time stability in the range of 1 ps over 1000 s and an accuracy better than 100 ps. The project is in operation onboard the Jason-2 satellite since June 2008. The principle is based on the Satellite Laser Ranging (SLR) technology; it uses the input of 20–25 SLR stations of the international laser network which participate in the tracking. This paper focuses on the data reduction process which was developed specifically to transform the raw information given by both space instrument and ground network: first to identify the triplets (ground and onboard epochs and time of flight of the laser pulse), second to estimate a usable product in terms of ground-to-space time transfer (including instrumental corrections), and thirdly to produce synchronization between any pair of remote ground clocks. In describing the validation of time synchronizations, the paper opens a way for monitoring the time difference between ultra-stable clocks thanks to a laser link at a few ps level for Common View passes. It highlights however that without accurately characterizing the onboard oscillator of Jason-2 and knowing the unavailability of time calibrations of SLR stations generally, time transfer over intercontinental distances remain difficult to be accurately estimated.     

Dynamics of the convection in the inner magnetosphere by observations of the diffuse aurora,SAR arc and ionospheric drift

The first results of the comparison of subauroral luminosity dynamics in 557,7 and 630,0 nm emission with simultaneous measurements of the ionospheric drift in the F2 region with a digisonde DPS-4 at the Yakutsk meridian (CGMC: 55–60N, 200°E) at Kp = 2–6 are presented. It is shown from the analysis of individual events that during the magnetospheric convection intensification after the turn of the IMF B_z – component to the south the equatorward extension of diffuse aurora takes place. At the same time the westward ionospheric drift velocity increases both in the diffuse aurora region and much equatorward of it due to the occurrence of the northward polarization electric field. We suppose that the generation of polarization field can be associated with the development of the region 2 FAC during the intensification of magnetospheric convection. The comparison of ground-based observations with measurements of the plasma drift aboard the DMSP-F15 satellite has been carried out.     

Ionospheric precursors of earthquakes recorded by VLF receiver at Tashkent IHY station

Tashkent International Heliophysical Year (IHY) station is a member of Atmospheric Weather Electromagnetic System for Observation, Modeling and Education (AWESOME) network being operated globally to study the ionosphere and the magnetosphere with the help of electromagnetic waves in Very Low Frequency (VLF) band. Regular monitoring of the D- and F-layers of ionosphere over Central Asia territory is being performed on the permanent basis starting year 2008. We have studied VLF amplitude anomalies related to the EQs occurred in 2008–2009 years with magnitude more than 5 on the path way from the VLF transmitters to the Tashkent station assuming that propagation of VLF ground-based transmitters signals can be perturbed by EQ preparation detectable from the ground-based measurements in the VLF bands. For analyzing narrowband data we have used the nighttime fluctuation (NF) method paying attention to the data obtained during the local nighttime (20:00 LT–04:00 LT) in Tashkent where the VLF receiver is operating. The mean nighttime amplitude (or trend) and nighttime fluctuation are found to increase significantly before the EQ occurred on the path way from the transmitters to the receiver. The obtained results have revealed an agreement with VLF amplitude anomalies observed in Tashkent VLF station during the strong EQs occurred on the path way from the transmitters to the receiver. Some results are presented to show the probing potentiality of VLF waves to predict short term EQs with high magnitude.     

DEMETER卫星探测到的强震前O^＋浓度变化

利用DEMETER卫星观测数据,采用level-2图像目视解译的方法和重访轨道对比方法,对2006年5月3日发生在汤加群岛的7．9级地震进行分析．从震前6天、震前2天、震前1天和震后2天的level-2图像可看到,电场和等离子体参量在震中上空偏北的位置出现扰动信号,震前的扰动信号幅度随地震时间临近有所增强,扰动持续时间也相应增加;震后2天扰动信号依然存在,其幅度并未减弱．通过对比震前1天的轨道在2006年至2008年内的重访轨道,发现O^＋浓度存在三类较规律的变化特征,即在夏季时,北半球的值高于南半球;冬季时则相反;二分季节时,南北半球含量相当．本文给出了地震和磁暴对O^＋浓度的不同影响,地震影响最靠近震中位置轨道上的O^＋浓度,而磁暴影响全球轨道上O^＋浓度分布．同时获取了DEMETER卫星上朗缪尔探针在上述重访轨道同步探测的电子浓度数据,发现电子浓度与O^＋浓度存在同步增加．等离子体浓度增加的现象可能与地震产生异常电场进而影响电离层结构有关．     

GSFC DORIS contribution to ITRF2008

The NASA GSFC DORIS analysis center has provided weekly DORIS solutions from November 1992 to January 2009 (839 SINEX files) of station positions and Earth Orientation Parameters for inclusion in the DORIS contribution to ITRF2008. The NASA GSFC GEODYN orbit determination software was used to process the orbits and produce the normal equations. The weekly SINEX gscwd10 submissions included DORIS data from Envisat, TOPEX/Poseidon, SPOT-2, SPOT-3, SPOT-4, SPOT-5. The orbits were mostly seven days in length (except for weeks with data gaps or maneuvers). The processing used the GRACE-derived EIGEN-GL04S1 gravity model, updated modeling for time-variable gravity, the GOT4.7 ocean tide model and tuned satellite-specific macromodels for SPOT-2, SPOT-3, SPOT-4, SPOT-5 and TOPEX/Poseidon. The University College London (UCL) radiation pressure model for Envisat improves nonconservative force modeling for this satellite, reducing the median residual empirical daily along-track accelerations from 3.75 × 10⁻⁹ m/s² with the a priori macromodel to 0.99 × 10⁻⁹ m/s² with the UCL model. For the SPOT and Envisat DORIS satellite orbits from 2003 to 2008, we obtain average RMS overlaps of 0.8–0.9 cm in the radial direction, 2.1–3.4 cm cross-track, and 1.7–2.3 cm along-track. The RMS orbit differences between Envisat DORIS-only and SLR & DORIS orbits are 1.1 cm radially, 6.4 cm along-track and 3.7 cm cross-track and are characterized by systematic along-track mean offsets due to the Envisat DORIS system time bias of ±5–10 μs. We obtain a good agreement between the geometrically-determined geocenter parameters and geocenter parameters determined dynamically from analysis of the degree one terms of the geopotential. The intrinsic RMS weekly position repeatability with respect to the IDS-3 combination ranges from 2.5 to 3.0 cm in 1993–1994 to 1.5 cm in 2007–2008.     

Statistical study of electrostatic solitary waves associated with reconnection: Geotail observations

The role of waves in the dynamics of the magnetotail has long been a topic of interest in magnetospheric physics. The characteristics of Electrostatic Solitary Waves (ESWs) associated with reconnection have been studied statistically in the magnetotail by surveying the large amounts data obtained from Waveform Capture (WFC) which is an important component of Plasma Wave Instrument (PWI) on the Geotail spacecraft. About 150 reconnection events with WFC data available are selected, and approximately 10 thousands of ESW waveforms are picked up by hands for statistical study. The ESWs are observed near diffusion region and near the plasma sheet boundary layer (PSBL). Two kinds of waveforms of ESWs are observed: bi-polar and tri-polar pulses. It is found that the pulse width of the ESWs is in the order of 1–5 ms and the peak-to-peak amplitude is in the order of 0.1–5 mV/m. The amplitudes of ESWs are larger in the near-earth tail region than that in deep tail region. ESWs have been observed with or without guide magnetic field 〈By〉. The characteristics of ESWs in different reconnection region and under different strength of guild magnetic field, their possible generation mechanism will be discussed.