中国近海和西北太平洋温跃层时空变化分析、模拟及预报

温度跃层是反映海洋温度场的重要物理特性指标，对水下通讯、潜艇活动及渔业养殖、捕捞等有重要影响。本文利用中国科学院海洋研究所“中国海洋科学数据库”在中国近海及西北太平洋（110ºE-140ºE,10ºN-40ºN）的多年历史资料（1930-2002年，510143站次），基于一种改进的温跃层判定方法，分析了该海域温跃层特征量的时空分布状况。同时利用Princeton Ocean Model(POM)，对中国近海，特别是东南沿海的水文结构进行了模拟，研究了海洋水文环境对逆温跃层的影响。最后根据历史海温观测资料，利用EOF分解统计技术，提出了一种适于我国近海及毗邻海域，基于现场有限层实测海温数据，快速重构海洋水温垂直结构的统计预报方法，以达到对现场温跃层的快速估计。 历史资料分析结果表明，受太阳辐射和风应力的影响，20°N以北研究海域，温跃层季节变化明显，夏季温跃层最浅、最强，冬季相反，温跃层厚度的相位明显滞后于其他变量，其在春季最薄、秋季最厚。12月份到翌年3月份，渤、黄及东海西岸，呈无跃层结构，西北太平洋部分海域从1月到3月份，也基本无跃层结构。在黄海西和东岸以及台湾海峡附近的浅滩海域，由于风力搅拌和潮混合作用，温跃层出现概率常年较低。夏季，海水层化现象在近海陆架海域得到了加强，陆架海域温跃层强度季节性变化幅度(0.31°C/m)明显大于深水区(约0.05°C/m)，而前者温跃层深度和厚度的季节性变化幅度小于后者。20°N以南研究海域，温跃层季节变化不明显。逆温跃层主要出现在冬、春季节（10月-翌年5月）。受长江冲淡水和台湾暖流的影响，东南沿海区域逆温跃层持续时间最长，出现概率最大，而在山东半岛北及东沿岸、朝鲜半岛西及北岸，逆温跃层消长过程似乎和黄海暖流有关。多温跃层结构常年出现于北赤道流及对马暖流区。在黑潮入侵黄、东、南海的区域，多温跃层呈现明显不同的季节变化。在黄海中部，春季多温跃层发生概率高于夏季和秋季，在东海西部，多跃层主要出现在夏季，在南海北部，冬季和春季多温跃层发生概率大于夏季和秋季。这些变化可能主要受海表面温度变化和风力驱动的表层流的影响。 利用Princeton Ocean Model(POM)，对中国东南沿海逆温跃层结构进行了模拟，模拟结果显示，长江冲淡水的季节性变化以及夏季转向与实际结果符合较好，基本再现了渤、黄、东海海域主要的环流、温盐场以及逆温跃层的分布特征和季节变化。通过数值实验发现，若无长江、黄河淡水输入，则在整个研究海域基本无逆温跃层出现，因此陆源淡水可能是河口附近逆温跃层出现的基本因素之一。长江以及暖流（黑潮和台湾暖流）流量的增加，均可在不同程度上使逆温跃层出现概率及强度、深度和厚度增加，且暖流的影响更加明显。长江对东南沿海逆温跃层的出现，特别是秋季到冬季初期，有明显的影响，使长江口海域逆温跃层位置偏向东南。暖流对于中国东南沿海的逆温跃层结构，特别是初春时期，有较大影响，使长江口海域的逆温跃层位置向东北偏移。 通过对温跃层长期变化分析得出，黄海冷水团区域，夏季温跃层强度存在3.8年左右的年际变化及18.9年左右的年代际变化，此变化可能主要表现为对当年夏季和前冬东亚地区大气气温的热力响应。东海冷涡区域，夏季温跃层强度存在3.7年的年际变化，在El Nino年为正的强度异常，其可能主要受局地气旋式大气环流变异所影响。谱分析同时表明，该海域夏季温跃层强度还存在33.2年的年代际变化，上世纪70年代中期，温跃层强度由弱转强，而此变化可能与黑潮流量的年代际变化有关。 海洋水温垂直结构的统计预报结果显示，EOF分解的前四个主分量即能够解释原空间点温度距平总方差的95%以上，以海洋表层附近观测资料求解的特征系数推断温度垂直结构分布的结果最稳定。利用东海陆架区、南海深水区和台湾周边海域三个不同区域的实测CTD样本廓线资料，对重构模型的检验结果表明，重构与实测廓线的相关程度超过95%的置信水平。三个区重构与实测温度廓线值的平均误差分别为0.69℃，0.52℃，1.18℃，平均重构廓线误差小于平均气候偏差，统计模式可以很好的估算温度廓线垂直结构。东海陆架海区温度垂直重构廓线与CTD观测廓线获得的温跃层结果对比表明，重构温跃层上界、下界深度和强度的平均绝对误差分别为1.51m、1.36m和0.17℃/m，它们的平均相对误差分别为24.7%、8.9%和22.6%，虽然温跃层深度和强度的平均相对误差较大，但其绝对误差量值较小。而在南海海区，模型重构温跃层上界、下界和强度的平均绝对预报误差分别为4.1m、27.7m和0.007℃/m，它们的平均相对误差分别为16.1%、16.8%和9.5%，重构温跃层各特征值的平均相对误差都在20%以内。虽然南海区温跃层下界深度平均绝对预报误差较大，但相对于温跃层下界深度的空间尺度变化而言（平均温跃层下界深度为168m），平均相对误差仅为16.8%。因此说模型重构的温度廓线可以达到对我国陆架海域、深水区温跃层的较好估算。 基于对历史水文温度廓线观测资料的分析及自主温跃层统计预报模型，研制了实时可利用微机简单、快捷地进行温跃层估算及查询的可视化系统，这是迄今进行大范围海域温跃层统计与实时预报研究的较系统成果。

Numerical study on the velocity structure around tidal fronts in the Yellow Sea

The velocity components across tidal fronts are examined using the Blumberg and Mellor 3-D nonlinear numerical coastal circulation model incorporated with the Mellor and Yamada level 2.5 turbulent closure model based on the reasonable model output of the M-2 tide and density residual currents. In the numerical experiments, upwelling motion appears around all the fronts with different velocity structures, accounting for surface cold water around the fronts. The experiments also suggest that the location and formation of fronts are closely related to topography and tidal mixing, as is the velocity structure around the front.

Preliminary study of the dynamic origin of the distribution pattern of bottom sediments on the continental shelves of the Bohai Sea, Yellow Sea and East China Sea

The bottom sediment types in the Bohai Sea, Yellow Sea and East China Sea (BYECS) are diversified, and their distribution pattern is very complicated. However, the bottom sediment types can be simplified to be sandy sediment, clayey sediment and mixed sediment, which comprise the complicated distribution pattern of bottom sediment in the BYECS. The continental shelves of the BYECS are broad, with shallow water depths and tidal currents which are permanent and dominate the marine dynamics in the BYECS. Based on numerical simulation of tidal elevations and currents in the BYECS, the rates of suspended load transport and bed load transport during a single tidal cycle for sediments of eight different grain size ranges are calculated. The results show that any sediment, whose threshold velocity is less than that of tidal current, has the same transport trend. Suspended load transport rare, bed load transport rate, and the ratio of the former to the latter decrease with grain size becoming coarser and coarser. The erosion/accretion patterns of sediments with different grain sizes are determined by the sediment transport rate divergences, and the results show that the patterns are the same for sediments with different grain sizes. Three main bottom sediment types, i.e. sandy sediment mainly composed of fine sand, clayey sediment mainly composed of silty clay, and mixed sediment mainly composed of fine sand, silt, and clay, are obtained by computation. The three bottom sediment types and their distribution pattern are consistent not only with sediment transport field and the sea bed erosion/accretion pattern obtained by simulation, but also with field data of bottom sediment types and divisions. In the BYECS, sand ridges form mainly in the areas with strong rectilinear tidal currents, sand sheets form mainly in the areas dominated by strong rotatory tidal currents, and clayey sediments, i.e. mud patches, form mainly in the areas with weak tidal currents. Hence, not only the sandy sediments but also the clayey sediments in the BYECS are formed under the control of the whole tidal current field of the BYECS. The three main bottom sediment types are not isolated respectively-in fact, they constitute a whole tidal depositional system. Under the condition with no cyclonic cold eddy, the clayey sediments in the BYECS can form in weak tidal current environments. Therefore, a cold eddy is not necessary for the deposition of clayey sediments in the BYECS. (C) 2000 Academic Press.

Restructuring of a zooplankton community by perturbation from a wind-forced coastal jet

The impact of transient wind events on an established zooplankton community was observed during a, field survey in a, coastal region off northern Norway in May 2002. A transient wind event induced a coastal jet/filament intrusion of warm, saline water into our survey area where a semi-permanent eddy was present. There was an abrupt change in zooplankton community structure within 4-7 days of the wind event, with a change in the size structure, an increase in lower size classes less than 1 mm in equivalent spherical diameter (ESD) and a decrease in larger size classes greater than 1.5 mm in ESD. The slope of zooplankton biovolume spectra changed from -0.6 to -0.8, consistent with the size shifting towards smaller size classes. This study shows that even well established zooplankton communities are susceptible to restructuring during transient wind events, and in particular when wind forcing induces horizontal currents or filaments.

A three-dimensional coupled physical-biological model study in the spring of 1993 in the Bohai Sea of China

A three-dimensional (3-D) coupled physical and biological model was used to investigate the physical processes and their influence on the ecosystem dynamics of the Bohai Sea of China. The physical processes include M-2 tide, time - varying wind forcing and river discharge. Wind records from I to 31 May in 1993 were selected to force the model. The biological model is based on a simple, nitrate and phosphate limited, lower trophic food web system. The simulated results showed that variation of residual currents forced by M, tide, river discharge and time-varying wind had great impact on the distribution of phytoplankton biomass in the Laizhou Bay. High phytoplankton biomass appeared in the upwelling region. Numerical experiments based on the barotropic model and baroclinic model with no wind and water discharge were also conducted. Differences in the results by the baroclinic model and the barotropic model were significant: more patches appeared in the baroclinic model comparing with the barotropic model. And in the baroclinic model, the subsurface maximum phytoplankton biomass patches formed in the stratified water.

细粒沉积物多组分粒度特征及其成因机制研究

By using high-resolution laser grain size instrument Mastersizer 2000, the grain size distribution of windblown depositions (loess and sandy dunes), aqueous sediments (lake, river, riverside and foreshore sand), weathering crust, sloping materials and other fine-grain sediments are systemically measured. The multimodal characteristics of grain size distribution of these sediments are carefully studied. The standard patterns and their grain size characteristics of various sediments are systemically summarized. The discrepancies of multimodal distribution among windblown depositions, aqueous sediments and other sediments are concluded and the physical mechanisms of grain size multimodal distribution of various sediments are also discussed in this paper. The major conclusions are followed: 1. The multimodal characteristic of grain size distribution is a common feature in all sediments and results from properties of transportation medium, dynamic intensity, transportation manner and other factors. 2. The windblown depositions are controlled by aerodynamic forcing, resulting in that the median size of the predominant mode gradually decreases form sandy dunes to loess. Similarly, the aqueous sediments are impacted by dynamic forces of water currents and the median grain size of the predominant mode decreases gradually from river to lake sediments. Because the kinetic viscidity of air is lower than of water, the grain size of modes of windblown depositions is usually finer than that of corresponding modes of aqueous sediments. Typical characteristics of sediments grain size distribution of various sediments have been summarized in the paper: (1) Suspended particles which diameters are less than 75μm are dominant in loess and dust. There are three modes in loess’ grain size distribution: fine, median and coarse (the median size is <1μm、1-10μm、10-75μm, respectively). The coarse mode which percentage is larger than that of others is controlled by source distance and aerodynamic intensity of dust source areas. Some samples also have a saltation mode which median size is about 300-500μm. Our analysis demonstrates that the interaction of wind, atmospheric turbulence, and dust grain gravity along the dust transportation path results in a multimodal grain size distribution for suspended dust. Changes in the median sizes of the coarse and medium modes are related to variation in aerodynamic forcing (lift force related to vertical wind and turbulence) during dust entrainment in the source area and turbulence intensity in the depositional area. (2) There is a predominant coarse saltation mode in grain size distribution of sandy dunes, which median size is about 100-300μm and the content is larger than that of other modes. The grain size distribution curve is near axis symmetric as a standard logarithm normal function. There are some suspended particles in some samples of sandy dunes, which distribution of the fine part is similar to that of loess. Comparing with sandy samples of river sediments, the sorting property of sandy dunes is better than of river samples although both they are the saltation mode. Thus, the sorting property is a criterion to distinguish dune sands and river sands. (3) There are 5~6 modes (median size are <1μm, 1-10μm, 10-70μm, 70-150μm, 150-400μm, >400μm respectively) in grain size distribution of lacustrine sediments. The former 4 modes are suspensive and others are saltated. Lacustrine sediments can be divided into three types: lake shore facies, transitional facies and central lake facies. The grain size distributions of the three facies are distinctly different and, at the same time, the transition among three modes is also clear. In all these modes, the third mode is a criteria to identify the windblown deposition in the watershed. In lake shore sediments, suspended particles are dominant, a saltation mode sometimes occurs and the fourth mode is the most important mode. In the transitional facies, the percentage of the fourth mode decreases and that of the second mode increases from lake shore to central lake. In the central lake facies, the second mode is dominant. A higher content of the second mode indicates its position more close to the central lake. (4) The grain size distribution of river sediments is the most complex. It consist of suspension, saltation and rolling modes. In most situations, the percentage of the saltation mode is larger than that of other modes. The percentage of suspension modes of river sediments is more than of sandy dunes. The grain size distribution of river sediments indicates dynamic strength of river currents. If the fourth mode is dominant, the dynamic forcing of river is weaker, such as in river floodplain. If the five or sixth mode is dominant, the water dynamic forcing of rivers is strong. (5) Sediments can be changed by later forcing in different degree to form some complicated deposition types. In the paper, the grain size distribution of aqueous sediments of windblown deposition, windblown sediments of aqueous deposition, weathering crust and slope materials are discussed and analyzed. 3. The grain size distribution characteristics of different sediments are concluded: (1) Modal difference: Usually there are suspended and saltation modes in the windblown deposition. The third mode is dominant in loess dust and the fifth mode is predominant in sandy dunes. There are suspended, saltation and rolling particles in aqueous sediments. In lacustrine sediments, the second and fourth mode are predominant for central lake facies or lake shore facies, respectively. In river sediments, the fourth, or fifth, or sixth mode is predominant. Suspended modes: the grain size of suspended particles of windblown depositions usually is less than 75μm. The content of suspended particles is lower or none in sandy dunes. However, suspended particles of aqueous sediments may reach 150μm. Difference in grain size of suspended modes represents difference between transitional mediums and the strength of dynamic forcing. Saltation modes: the median size of saltation mode of sandy dunes fluctuates less than that of river sediments. (2) Loess dust and lacustrine sediment: Their suspended particles are clearly different. There is an obvious pit between the second and the third modes in grain size distribution of lacustrine sediments. The phenomenon doesn’t occur in loess dust. In lacustrine sediments, the second mode can be a dominant mode, such as central lake facies, and contents of the second and the third modes change reversely. However, the percentage of the third mode is always the highest in loess dust. (3) Dune Sand and fluvial sand: In these two depositions, the saltation particles are dominant and the median sizes of their saltation modes overlay in distribution range. The fifth mode of dune sand fluctuates is sorted better than that of fluvial sand. (4) Lacustrine and fluvial sediments: In lacustrine sediments, there are 5-6 modes and suspended particles can be predominant. The second mode is dominant in central lake facies and the third mode is dominant in lake shore facies. Saltation or roll modes occurred in central lake facies may indicate strong precipitation events. In fluvial sediments, saltation particles (or rolling particles) usually dominant. 4. A estimation model of lake depth is firstly established by using contents of the second, the third and the fourth modes. 5. The paleoenvironmental history of the eastern part of SongLiao basin is also discussed by analyzing the grain size distribution of Yushu loess-like sediments in Jilin. It was found that there is a tectonic movement before 40ka B.P. in SongLiao basin. After the movement, loess dust deposited in Yushu area as keerqin desert developed. In recent 2000 years, the climate became drier and more deserts activated in the eastern part of Song-Liao basin.

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

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.

AE指数的可靠性研究

Auroral electrojet index (AE) are usually used to quantitatively describe the activity of the geomagnetic field in the polar region. AE is a means to identify the level of a substorm as well. The auroral electrojet indices (AU, AL, and AE) have served well for more than four decades as measures of magnetospheric substorm activity. However, as substorm studies have progressed considerably during the past several years, the accuracy of the present electrojet indices have become an important issue. Thus it is fortunate for us to reexamine and evaluate the accuracy of the present electrojet indices and improve them if necessary. For a better use of the present indices and for future improvement we examine the limitations of the auroral electrojet indices as an accurate quantitative measure of the auroral electrojets and magnetospheric substorms. Some of the limitations arise from the data availability and also from the present simplified scheme in deriving them, but some of them originate in the definition themselves. In the present paper, we analyze the characteristics of aurora, ionospheric current system and the AE index. It is noted that for the pseudo-breakup events the variations of AE ( > 500 nT) clearly show the sudden increase and slow decrease phases. However, aurora does not expand to poleward or equatorward, and the ionospheric currents presents the features of the magnetic convection. We mainly focus on investigating what time the AE could be used to identify the subtorm.

低感应数条件下大地电导率测量方法研究

A frequency domain electromagnetic (conductivity) method for near surface soundings at low frequencies is discussed in this thesis. Its elementary principle is to detect the conductivity of the earth by the secondary magnetic fields induced by a current dipole on the earth. According to the EM induction theory, a coil with alternating current on the earth will generate a magnetic field in whole space which is referred to as the primary field Hp. The primary field would induce secondary currents in the earth which go down to depth like a batch of smoking rings. These currents further produce secondary magnetic field Hs .The primary and secondary magnetic fields are collected together by a receiver coil. Generally speaking，the secondary magnetic field is a complicated function of coil spacing, transmitting frequency and earth conductivity. But at low induction numbers, the secondary field is deduced to as a simple function of frequency, spacing and conductivity. Especially the ratio of secondary to primary field shares a linear proportion to the apparent conductivity. The earth conductivity can be interpreted by proper inversions with the apparent conductivity. The method is discussed at three steps: (1)Derivation of primary and secondary magnetic fields arising from vertical and horizontal magnetic dipoles on the earth based on the basic EM induction theory. (2)Field techniques and equipment developed for the method. (3)An interpretation technique was introduced using a cumulative and relative response function. Finally a test example is presented for examining the effectiveness of the method.

华南震旦纪-寒武纪两期成磷事件及其地球动力学意义

Phosphorus is an important biological and ecological element that to a certain degree constrains ecological environment and nutrient (including carbon) cycling. Marine sedimentary phosphorites are the principal phosphorus supply of the mankind. In the eastern to southern margins of the Yangtze Craton, South China, there are two phosphogenetic events at the Doushantuo stage of the Late Sinian and the Meishucun stage of the Early Cambrian respectively, corresponding two explosion events of life across the Precambrian\Cambrian boundary. Phosphorus ores from the Sinian and Cambrian phosphate in South China can be classified roughly into two categories, namely, grained and non-grained phosphorites. Grained phosphorites, hosted in dolostone type of phosphogenetic sequences and with larger industrial values, occur mainly in margins of the Upper Yangtze Platform, formed in shallow-water environments with high hydraulic energy and influenced by frequent sea-level change. Non-grained phosphorites, hosted principally in black-shale type of phosphogenetic sequences and with smaller industrial values, are distributed mainly in the Jiangnan region where deeper-water sub-basins with low hydraulic energy were prevailing at the time of phosphogenesis. Secular change ofδ~(13)C, δ~(18) O, ~(86)Sr/~(87)Sr values of carbonates from Sinian and Cambrian sequences were determined. A negative abnormal ofδ~(13)C, δ~(18)O values and positive abnormal of 86Sr/87Sr values from the fossiliferous section of the Lowest Cambrian Meishucun Formation implies life depopulation and following explosion of life across the PrecambriamCambrian boundary. Based on a lot of observations, this paper put forward a six-stage genetic model describing the whole formational process of industrial phosphorites: 1) Phosphorus was transported from continental weathering products and stored in the ocean; 2) dissolved phosphates in the seawater were enriched in specific deep seawater layer; 3) coastal upwelling currents took this phosphorus-rich seawater to a specific coastal area where phosphorus was captured by oceanic microbes; 4) clastic sediments in this upwelling area were enriched in phosphorus because of abundant phosphorus-rich organic matters and because of phosphorus absorption on grain surfaces; 5) during early diagenesis, the phosphorus enriched in the clastic sediments was released into interstitial water by decomposition and desorption, and then transported to the oxidation-reduction interface where authigenic phosphates were deposited and enriched; 6) such authigenic phosphate-rich layers were scoured, broken up, and winnowed in shallow-water environments resulting in phosphate enrichment. The Sinian-Cambrian phosphorites in South China are in many aspects comparable with coastal-upwelling phosphorites of younger geological ages, especially with phosphorites from modern coastal upwelling areas. That implies the similarities between the Sinian-Cambrian ocean and the modern ocean. Although Sinian-Cambrian oceanic life was much simpler than modern one, but similar oceanic planktons prevail, because oceanic planktons (particularly phytoplanktons) are crucial for phosphate enrichment related to coastal upwelling. It implies also a similar seawater-layering pattern between the Sinian-Cambrian ocean and the modern ocean. The two global phosphate-forming events and corresponding life-explosion events at the Sinian and Cambrian time probably resulted from dissolved-phosphate accumulation in seawater over a critical concentration during the Earth's evolution. Such an oceanic system with seawater phosphorus supersaturation is evidently unstable, and trends to return to normal state through phosphate deposition. Accordingly, this paper put forward a new conception of "normal state <=> phosphorus-supersaturation state" cycling of oceanic system. Such "normal state <=> phosphorus-supersaturation state" cycling was not only important for the three well-known global phosphate-forming events, also related to the critical moments of life evolution on the Earth. It might be of special significance. The favorable paleo-oceanic orientation in regard to coastal-upwelling phosphorite formation suggests a different orientation of the Yangtze Craton between the Sinian time and the present time (with a 135° clockwise difference), and a 25° anti-clockwise rotation of the Yangtze Craton from late Sinian to early Cambrian. During the Sinian-Cambrian time, the Yangtze Craton might be separated from the Cathaysia Block, but might be still associated with the North China Craton.