空间电流体系地磁反演及理论研究

The space currents definitely take effects on electromagnetic environment and also are scientific highlight in the space research. Space currents as a momentum and energy provider to Geospace Storm, disturb the varied part of geomagnetic field, distort magnetospheric configuration and furthermore take control of the coupling between magnetosphere and ionosphere. Due to both academic and commercial objectives above, we carry on geomagnetic inverse and theoretical studies about the space currents by using geomagnetic data from INTERMAGNET. At first, we apply a method of Natural Orthogonal Components (NOC) to decomposition the solar daily variation, especially for (solar quiet variation). NOC is just one of eign mode analysis, the most advantage of this method is that the basic functions (BFs) were not previously designated, but naturally came from the original data so that there are several BFs usually corresponding to the process really happened and have more physical meaning than the traditional spectrum analysis with the fixed BFs like Fourier trigonometric functions. The first two eign modes are corresponding to the and daily variation and their amplitudes both have the seasonal and day-to-day trend, that will be useful for evaluating geomagnetic activity indices. Because of the too strict constraints of orthogonality, we try to extend orthogonal contraints to the non-orthogonal ones in order to give more suitable and appropriate decomposition of the real processes when the most components did not satisfy orthogonality. We introduce a mapping matrix which can transform the real physical space to a new mathematical space, after that process, the modified components which associated with the physical processes have satisfied the orthogonality in the new mathematical space, furthermore, we can continue to use the NOC decomposition in the new mathematical space, and then all the components inversely transform back to original physical space, so that we would have finished the non-orthogonal decomposition which more generally in the real world. Secondly, geomagnetic inverse of the ring current’s topology is conducted. Configurational changes of the ring current in the magnetosphere lead to different patterns of disturbed ground field, so that the global configuration of ring current can be inferred from its geomagnetic perturbations. We took advantages of worldwide geomagnetic observatories network to investigate the disturbed geomagnetic field which produced by ring current. It was found that the ring current was not always centered at geomagnetic equator, and significantly deviated off the equator during several intense magnetic storms. The deviation owing to the tilting and latitudinal shifting of the ring current with respect to the earth’s dipole can be estimated from global geomagnetic survey. Furthermore those two configurational factors which gave a quantitative description of the ring current configuration, will be helpful to improve the Dst calibration and understand the dependence of ring current’s configuration on the plasma sheet location relative to the equator when magnetotail field warped. Thirdly, the energization and physical acceleration process of ring current during magnetic storm has been proposed. When IMF Bz component increase, the enhanced convection electric field drive the plasma injection into the inner magnetosphere. During the transport process, a dynamic heating is happened which make the particles more ‘hot’ when the injection is more deeply inward. The energy gradient along the injection path is equivalent to a kind of force, which resist the plasma more earthward injection, as a diamagnetic effect of the magnetosphere anti and repellent action to the exotically injected plasma. The acceleration efficiency has a power law form. We use analytical way to quantitatively describe the dynamical process by introducing a physical parameter: energization index, which will be useful to understand how the particle is heated. At the end, we give a scheme of how to get the from storm time geomagnetic data. During intense magnetic storms, the lognormal trend of geomagnetic Dst decreases depend on the heating dynamic of magnetosphere controlling ring current. The descending pattern of main phase is governed by the magnetospheric configuration, which can be describled by the energization index. The amplitude of Dst correlated with convection electric field or south component of the solar wind. Finally, the Dst index is predicted by upstream solar wind parameter. As we known space weather have posed many chanllenges and impacts on techinal system, the geomagnetic index for evaluating the activity space weather. We review the most popular Dst prediction method and repeat the Dst forecasting model works. A concise and convnient Key Points model of the polar region is also introduced to space weather. In summary, this paper contains some new quantitative and physical description of the space currents with special focus on the ring current. Whatever we do is just to gain a better understanding of the natural world, particularly the space environment around Earth through analytical deduction, algorithm designing and physical analysis, to quantitative interpretation. Applications of theoretical physics in conjunction with data analysis help us to understand the basic physical process govering the universe.

磁层-电离层耦合现象的分析研究

The magnetosphere-ionosphere coupling is mainly manifested by the trans- porting processes of energy into the ionosphere , the energy is carried by solar wind and firstly accumulate at the magnetosphere, and the coupling processes also significantly include the interaction between the magnetosphere and ionosphere for mass and energy. At the quiet condition, energy is delivered by the large-scale convection of the geomagnetic field; the huge energy from solar wind bulk will be injected into and consumed at the near magnetosphere and ionosphere by the geomagnetic storm and substorm activities. Aurorae and FACs (Field-aligned currents) are the important phenomena in the coupling processes. In the present work, firstly, we analyze the activity characteristics of auroral precipitating particle, secondly, we study the distribution characters of large-scale field aligned currents (LS FACs) at storm-time using the observations from different satellites at different altitudes. Finally, we investigate the evolution of the geomagnetic field configuration at the nightside sector on the onset of the expansion phase in a substorm event, the substorm event happened at 0430UT to 0630UT on 8th Nov. 2004. The main results as follows： At the first, the data of the estimated power input (EPI) of auroral particles from NOAA/POES (Polar orbiting environmental satellite) for some 30 years have been analyzed. The variation tendencies of the EPI generally coincide with aa, AE and Dst indices. The annual variation of EPI shows equinox peaks and an asymmetric-activity with a higher peak in the winter-hemisphere than in the summer-hemisphere. The diurnal UT variations are different from north and south hemisphere: for north hemisphere, the peak appears at 1200UT, and the relative deviation is 22% to the daily average of the north hemisphere. For south hemisphere, the maximal deviation is 22% at 2000UT. So the diurnal variation of EPI is more dominant than the annual variation which maximal deviation is 3% to 12% for different seasons. Studies on correlations of the hourly average of EPI, Pa, with AE and Dst indices show a correlation coefficient r=0.74 of Pa and AE, and r=-0.55 of Pa and Dst. The hourly EPIs for north and south polar regions, NPa and SPa, show a north-south asymmetry with a higher correlation of SPa and AE (or Dst). Time delays of EPI with respect to magnetic indices are examined, the maximum correlation coefficient of Pa with AE (r=0.78) occurs when the time delay =0, suggesting a synchronous activity of auroral electrojet and auroral precipitating particles, while =1-2h, the correlation coefficient of Pa with Dst is maximum (r=0.57), suggesting that the activity of auroral particle precipitating may influence the ring current on some extent. Sencondly, we use the high-resolution magnetic field vector data of the CHAMP satellite to investigate the distribution of large-scale FACs during the great magnetic storm on 7th to 8th Nov. 2004. The results show that, whether in the northern or southern hemisphere, the number and density of large-scale FACs during the main-phase are more and bigger than these during the recover-phase, and the number of large-scale FACs in morning sector obviously is more than that in afternoon sector. In terms of the magnetic indices, we find that large-scale FACs in morning sector significantly affected by the substorm activities, while in afternoon sector the large-scale FACs mainly indicate the fluctuations of the ring-current in storm time. Accordingly to the former studies, similarly, we find that in the morning sector, the scale of the large-scale FACs move to the high-latitude region, and in the afternoon sector, large-scale FACs distinctly expand to the low-latitude region. During the time periods that the NOAA/POES auroral precipitating particle power data temporally correspond to the large-scale FACs, the more the power of auroral particle is, the more and bigger the number and density of FACs are. At the same time, we use the magnetic field vector data of POLAR obtain a good form of region 1, region 2, and three pieces of cusp FACs during a single transit at 1930UT-2006UT on 07th. And the characteristics of simultaneous electric field and energy particles observations on Polar are coincide with the five FACs pieces. Finally, by means of the observation of Cluster 4 and Goes 10、 Goes 12, we analyze the evolution process of the change of the magnetic field configuration at night sector at the expansion phase of a substorm event which happened during 0430UT to 0630UT on 8th Nov. 2004, we find that the times of the beginning of the polarizations of magnetic field are observed from Goes 10 to Goes 12 then to Cluster 4. So, at the synchronous orbit ( 6.6 RE) to 10RE distance scale of the neutral sheet, the current disruption spread tailward. Simultaneously, the strengthen of the FACs deduced from these satellites’ magnetic field observations are almost consistent with the times of polarizations, as well as the high energy particles injection and the electric field dominant variation. The onset times determined by the magnetic field polarizations from these satellites are all ahead of the onset time that confirmed from the auroral electrojet indices. So, these characters of different observations can be used as the criterions to determine the onset time for the substorms of such type as we studied.

基于广义Lippmann-Schwinger方程的地震波场模拟算法与应用

In the last several decades, due to the fast development of computer, numerical simulation has been an indispensable tool in scientific research. Numerical simulation methods which based on partial difference operators such as Finite Difference Method (FDM) and Finite Element Method (FEM) have been widely used. However, in the realm of seismology and seismic prospecting, one usually meets with geological models which have piece-wise heterogeneous structures as well as volume heterogeneities between layers, the continuity of displacement and stress across the irregular layers and seismic wave scattering induced by the perturbation of the volume usually bring in error when using conventional methods based on difference operators. The method discussed in this paper is based on elastic theory and integral theory. Seismic wave equation in the frequency domain is transformed into a generalized Lippmann-Schwinger equation, in which the seismic wavefield contributed by the background is expressed by the boundary integral equation and the scattering by the volume heterogeneities is considered. Boundary element-volume integral method based on this equation has advantages of Boundary Element Method (BEM), such as reducing one dimension of the model, explicit use the displacement and stress continuity across irregular interfaces, high precision, satisfying the boundary at infinite, etc. Also, this method could accurately simulate the seismic scattering by the volume heterogeneities. In this paper, the concrete Lippmann-Schwinger equation is specifically given according to the real geological models. Also, the complete coefficients of the non-smooth point for the integral equation are introduced. Because Boundary Element-Volume integral equation method uses fundamental solutions which are singular when the source point and the field are very close，both in the two dimensional and the three dimensional case, the treatment of the singular kernel affects the precision of this method. The method based on integral transform and integration by parts could treat the points on the boundary and inside the domain. It could transform the singular integral into an analytical one both in two dimensional and in three dimensional cases and thus it could eliminate the singularity. In order to analyze the elastic seismic wave scattering due to regional irregular topographies, the analytical solution for problems of this type is discussed and the analytical solution of P waves by multiple canyons is given. For the boundary reflection, the method used here is infinite boundary element absorbing boundary developed by a pervious researcher. The comparison between the analytical solutions and concrete numerical examples validate the efficiency of this method. We thoroughly discussed the sampling frequency in elastic wave simulation and find that, for a general case, three elements per wavelength is sufficient, however, when the problem is too complex, more elements per wavelength are necessary. Also, the seismic response in the frequency domain of the canyons with different types of random heterogeneities is illustrated. We analyzed the model of the random media, the horizontal and vertical correlation length, the standard deviation, and the dimensionless frequency how to affect the seismic wave amplification on the ground, and thus provide a basis for the choice of the parameter of random media during numerical simulation.

中国大陆及邻区地震活动背景、环境触发与地震预测研究

This dissertation that includes most of the P. PH.D research work during 2001~2002 covers the large-scale distribution of continental earthquakes in mainland China, the mechanism and statistic features of grouped strong earthquakes related to the tidal triggering, some results in earthquake prediction with correlativity analysis methods, and the flushes from the two strong continental earthquakes in South Asia in 2001. Mainland China is the only continental sub-plate that is compressed by collision boundaries at its two sides, within which earthquakes are dispersive and distributed as seismic belts with different widths. The control capability of the continental block boundaries on the strong earthquakes and seismic hazards is calculated and analyzed in this dissertation. By mapping the distribution of the 31282 ML:3s2,0 earthquakes, I found that the depth of continental earthquakes depend on the tectonic zonings. The events on the boundaries of relatively integrated blocks are deep and those on the new-developed ruptures are shallow. The average depth of earthquakes in the West of China is about 5km deeper than that in the east. The western and southwestern brim of Tarim Basin generated the deepest earthquakes in mainland China. The statistic results from correlation between the grouped M7 earthquakes and the tidal stress show that the strong events were modulated by tidal stress in active periods. Taking Taiwan area as an example, the dependence of moderate events on the moon phase angles (￡>) is analyzed, which shows that the number of the earthquakes in Taiwan when D is 50° ,50° +90° and 50° +180° is more than 2 times of standard deviation over the average frequency at each degree, corresponding to the 4th, 12th and 19th solar day after the new moon. The probability of earthquake attack to the densely populated Taiwan island on the 4th solar day is about 4 times of that on other solar days. On the practice of earthquake prediction, I calculated and analyzed the temporal correlation of the earthquakes in Xinjinag area, Qinghai-Tibet area, west Yunnan area, North China area and those in their adjacent areas, and predicted at the end of 2000 that it would be a special time interval from 2001 to 2003, within which moderate to strong earthquakes would be more active in the west of China. What happened in 2001 partly validated the prediction. Within 10 months, there were 2 great continental earthquakes in south Asia, i.e., the M7.8 event in India on Jan 26 and M8.1 event in China on Nov. 14, 2001, which are the largest earthquake in the past 50 years both for India and China. No records for two great earthquakes in Asia within so short time interval. We should speculate the following aspects from the two incidences: The influence of the fallacious deployment of seismic stations on the fine location and focal mechanism determination of strong earthquakes must be affronted. It is very important to introduce comparative seismology research to seismic hazard analysis and earthquake prediction research. The improvement or changes in real-time prediction of strong earthquakes with precursors is urged. Methods need to be refreshed to protect environment and historical relics in earthquake-prone areas.

东南极格罗夫山地区的土壤特征及其环境研究

The Grove Mountains, including 64 nunataks, is situated on an area about 3200km2 in the inland ice cap of east Antarctica in Princess Elizabeth land (72o20'-73°101S, 73°50'-75o40'E), between Zhongshan station and Dome A, about 450km away from Zhongshan station (69°22'S, 76°22'E). Many workers thought there was no pedogenesis in the areas because of the less precipitation and extreme lower temperature. However, during the austral summer in 1999-2000, the Chinaer 16 Antarctic expedition teams entered the inland East Antarctica and found three soil spots in the Southern Mount Harding, Grove Mountains, East Antarctica. It is the first case that soils are discovered in the inland in East Antarctica. Interestingly, the soils in this area show clay fraction migration, which is different from other cold desert soils. In addition, several moraine banks are discovered around the Mount Harding. The soil properties are discussed as below. Desert pavement commonly occurs on the three soil site surfaces, which is composed of pebbles and fragments formed slowly in typical desert zone. Many pebbles are subround and variegated. These pebbles are formed by abrasion caused by not only wind and wind selective transportation, but also salt weathering and thaw-freezing action on rocks. The wind blows the boulders and bedrocks with snow grains and small sands. This results in rock disintegration, paved on the soil surface, forming desert pavement, which protects the underground soil from wind-blow. The desert pavement is the typical feature in ice free zone in Antarctica. There developed desert varnish and ventifacts in this area. Rubification is a dominant process in cold desert Antarctic soils. In cold desert soils, rubification results in relatively high concentrations of Fed in soil profile. Stained depth increases progressively with time. The content of Fed is increasing up to surface in each profile. The reddish thin film is observed around the margin of mafic minerals such as biotite, hornblende, and magnetite in parent materials with the microscope analyzing on some soil profiles. So the Fed originates from the weathering of mafic minerals in soils. Accumulations of water-soluble salts, either as discrete horizons or dispersed within the soil, occur in the soil profiles, and the salt encrustations accumulate just beneath surface stones in this area. The results of X-ray diffraction analyses show that the crystalline salts consist of pentahydrite (MgSO4-5H2O), hexahydrite (MgSO4-6H2O), hurlbutite (CaBe2(PO4)2), bloedite (Na2Mg(S04)2-4H2O), et al., being mainly sulfate. The dominant cations in 1:5 soil-water extracts are Mg2+ and Na+, as well as Ca2+ and K+, while the dominant anion is SO42-, then NO3-, Cl- and HCO3-. There are white and yellowish sponge materials covered the stone underside surface, of which the main compounds are quartz (SiO2, 40.75%), rozenite (FeSOKkO, 37.39%), guyanaite (Cr2O3-1.5H2O, 9.30%), and starkeyite (MgSO4-4H2O, 12.56%). 4) The distribution of the clay fraction is related to the maximum content of moisture and salts. Clay fraction migration occurs in the soils, which is different from that of other cold desert soils. X-ray diffraction analyses show that the main clay minerals are illite, smectite, then illite-smectite, little kaolinite and veirniculite. Mica was changed to illite, even to vermiculite by hydration. Illite formed in the initial stage of weathering. The appearance of smectite suggests that it enriched in magnesium, but no strong eluviation, which belongs to cold and arid acid environment. 5) Three soil sites have different moisture. The effect moisture is in the form of little ice in site 1. There is no ice in site 2, and ice-cement horizon is 12 cm below the soil surface in site 3. Salt horizon is 5-10 cm up to the surface in Site 1 and Site 2, while about 26cm in site 3. The differentiation of the active layer and the permafrost are not distinct because of arid climate. The depth of active layer is about 10 cm in this area. Soils and Environment: On the basis of the characteristics of surface rocks, soil colors, horizon differentiation, salt in soils and soil depth, the soils age of the Grove Mountains is 0.5-3.5Ma. No remnants of glaciations are found on the soil sites of Mount Harding, which suggests that the Antarctic glaciations have not reached the soil sites since at least 0.5Ma, and the ice cap was not much higher than present, even during the Last Glacial Maximum. The average altitude of the contact line of level of blue ice and outcrop is 2050m, and the altitude of soil area is 2160m. The relative height deviation is about 110m, so the soils have developed and preserved until today. The parental material of the soils originated from alluvial sedimentary of baserocks nearby. Sporepollen were extracted from the soils, arbor pollen grains are dominant by Pinus and Betula, as well as a small amount Quercus, Juglans, Tilia and Artemisia etc. Judging from the shape and colour, the sporepollen group is likely attributed to Neogene or Pliocene in age. This indicates that there had been a warm period during the Neogene in the Grove Mountains, East Antarctica.

A new method of predicting of gas chromatographic retention indices for polychlorinated dibenzofurans (PCDFs)

A new method has been developed to describe the quantitative relationship between molecular structures of PCDFs and their gas chromatographic retention indices on a 30-m fused silica column coated with DB-5 stationary phase. The regression equation is derived with a multiple correlation coefficient greater than 0.9995. The highest residual is 20 index units. The standard deviation is less than 7 index units. Using this regression equation, the retention indices of PCDFs for which data is not available have also been predicted. (C) 2000 Elsevier Science Ltd. All rights reserved.

Separation of peptides by strong cation-exchange capillary electrochromatography

Separation of small peptides on ion-exchange capillary electrochromatography (IE-CEC) with strong cation-exchange packing (SCX) as stationary phase was investigated. It was observed that the number of theoretical plates for small peptides varied from 240 000 to 460 000/m, and the relative standard deviation for t(0) and the migration time of peptides were less than 0.57% and 0.27%, respectively for ten consecutive runs. Unusually high column efficiency has been explained by the capillary electrophoretic stacking and chromatofocusing phenomena during the injection and separation of positively charged peptides. The sample buffer concentration had a marked effect on the column efficiency and peak area of the retained peptides. The influences of the buffer concentration and pH value as well as the applied voltage on the separation were investigated. It has been shown that the electrostatic interaction between the positively charged peptides and the SCX stationary phase played a very important role in IE-CEC, which provided the different separation selectivity from those in the capillary electrophoresis and reversed-phase liquid chromatography. A fast separation of ten peptides in less than 3.5 min on IE-CEC by adoption of the highly applied voltage was demonstrated. (C) 2000 Elsevier Science B.V. All rights reserved.

Estimation of soil organic partition coefficients: from retention factors measured by soil column chromatography with water as eluent

The retention factors (k) of 104 hydrophobic organic chemicals (HOCs) were measured in soil column chromatography (SCC) over columns filled with three naturally occurring reference soils and eluted with Milli-Q water. A novel method for the estimation of soil organic partition coefficient (K-oc) was developed based on correlations with k in soil/water systems. Strong log K-oc versus log k correlations (r>0.96) were found. The estimated K-oc values were in accordance with the literature values with a maximum deviation of less than 0.4 log units. All estimated K-oc values from three soils were consistent with each other. The SCC approach is promising for fast screening of a large number of chemicals in their environmental applications. (C) 2002 Elsevier Science B.V. All rights reserved.

Algebraic approach to stretching vibrational spectrum of H2S

The stretching vibrational spectra of H2S have been calculated by using the algebraic model, and the local mode characteristics have been analyzed. Within the vibrational quanta v=5, the standard deviation is 1.71 cm(-1), which is better than that from the local model HCAO model calculation.

不同干扰条件下社会认知偏离者两半球协同功能研究

Previous studies have witnessed some psychological or behavioral deviation (such as aggressive behavior) might have an association with cerebral hemisphere cooperative dysfunction, however, it is still unclear whether there is an association between individuals with social cognitive bias and their hemispheric cooperative functions especially while the interhemisphere cooperative processing is under the conditions of emotional interferences. The purpose of this study is to explore the differences between the social cognitive bias group and the normal group’s interhemispheric cooperative functional activity under the conditions of with or without interferences. Methods: According to Dodge’s (1993) model of “social-cognitive mechanisms in the development of conduct disorder and depression”, a 51 items of “social cognitive bias scale” was created and was used to screen the high score group. 20 male subjects was composed of high score group and other 23 matched the control group. Stimulus tachistoscopically presented to the bilateral visual field and compared with the central. Both group’s interhemispheric cooperative functional activity were observed and compared under the conditions of without interference- i.e. base level and with the emotional interferences of white noise level and negative evaluative feedback speech level while finishing: experiment one: Chinese word-figure Stroop analogue task; experiment two: two single Chinese Characters combination task. Heart rate and respiratory rate were simultaneously recorded as index of emotional changes. Results: ① The high score group showed a decrease in processing accuracy compared with the normal group under the condition of white noise interference level in experiment one. ② Still under the condition of white noise interference level, there were more reaction time and more errors were observed in high score group than normal in experiment two. ③ Both groups showed speed up effect and the strategic processing tendency of speed-accuracy trade-off effect under the condition of white noise interference level in both experiments. ④ Between group differences of interhemipheric cooperative function were not observed under the conditions of base level and the negative evaluative feedback speech level within both experiments. Conclusion: The results suggested that interhemispheric cooperative functional differences exists between the two groups, characterized as ① differences existed in interhemispheric cooperative processing strategy between the two groups, with the high score group presented “hierarchic” deficiency strategy. ② the appearance of the differences between the two groups were condition specified , and in this research it was only under the white noise interference condition. ③ the features of the differences between the two groups were the differences on multidimensional performances and with a deficit orientation in high score group. ④ the varieties of the differences were changing with cooperative tasks, as in this research the high score group performed worse in complementary cooperative task. In addition, both group adjusted the processing strategy respectively under the condition of white noise evoked emotional interference implied that the interaction between the interhemisphere cooperative processing and emotion might exist.