THE CHARGE-DISTRIBUTION AND THE ELECTRIC-FIELD GRADIENT IN THE HIGH-TEMPERATURE SUPERCONDUCTOR TLBA2CAN-1CUNO2N+3 (N = 1 TO 3)

The net charges at atoms in the high-temperature superconductor TlBa2Can-1CunO2n+3 (n = 1 to 3) are calculated by means of the tight-binding approximation based on the EHMO method. The results indicate that the charge distribution in this kind of compounds possesses a specially layered arrangement. An insulating Ba-Ba layer is inserted between the Cu-O layer and the Tl-O layer. There may exist a weak coupling between the Cu-O layer and the Tl-O layer through the interaction of the same O(2) atom with both the Cu atom and the Tl atom. The existence of the Ca in the compounds can cause the valence fluctuation at the Cu atom. The calculated electric field gradients at atoms implies that the conducting electron or hole may move in the Cu-O layer, which is closest to the Tl-O layer, along the a-b plane.

Effective dielectric responses of graded composites of the particles with a general power-law gradation profile

The effective dielectric response of graded spherical composites having general power-law gradient inclusions is investigated under a uniform applied electric field, where the dielectric gradation profile of the spherical inclusions is modeled by the equation epsilon(i) (r) = c(b+r)(k). Analytical solutions of the local electrical potentials are derived in terms of hyper-geometric function and the effective dielectric response of the graded composites is predicted in the dilute limit. From our result, the local potentials of graded spherical composites having both simple power-law dielectric profile epsilon(i)(r) = cr(k) and linear dielectric profile epsilon(i) (r) = c(b+r) are derived exactly by taking the limits b --> 0 and k --> 1, respectively. In the dilute limit, our exact result is used to test the validity of differential effective dipole approximation (DEDA) for estimating the effective response of graded spherical composites, and it is shown that the DEDA is in excellent agreement with exact result. (C) 2005 Elsevier B.V. All rights reserved.

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

温度跃层是反映海洋温度场的重要物理特性指标，对水下通讯、潜艇活动及渔业养殖、捕捞等有重要影响。本文利用中国科学院海洋研究所“中国海洋科学数据库”在中国近海及西北太平洋（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%。因此说模型重构的温度廓线可以达到对我国陆架海域、深水区温跃层的较好估算。 基于对历史水文温度廓线观测资料的分析及自主温跃层统计预报模型，研制了实时可利用微机简单、快捷地进行温跃层估算及查询的可视化系统，这是迄今进行大范围海域温跃层统计与实时预报研究的较系统成果。

In vitro iron enrichment experiments in the Prydz Bay, the Southern Ocean: A test of the iron hypothesis

In "high nitrate, low chlorophyll" (HNLC) ocean regions, iron has been typically regarded as the limiting factor for phytoplankton production. This "iron hypothesis" needs to be tested in various oceanic environments to understand the role of iron in marine biological and biogeochemical processes. In this paper, three in vitro iron enrichment experiments were performed in Prydz Bay and at the Polar Front north of the Ross Sea, to study the role of iron on phytoplankton production. At the Polar Front of Ross Sea, iron addition significantly (P < 0.05, Student's t-test) stimulated phytoplankton growth. In Prydz Bay, however, both the iron treatments and the controls showed rapid phytoplankton growth, and no significant effect (P > 0.05, Student's t-test) as a consequence of iron addition was observed. These results confirmed the limiting role of iron in the Ross Sea and indicated that iron was not the primary factor limiting phytoplankton growth in Prydz Bay. Because the light environment for phytoplankton was enhanced in experimental bottles, light was assumed to be responsible for the rapid growth of phytoplankton in all treatments and to be the limiting factor controlling field phytoplankton growth in Prydz Bay. During the incubation experiments, nutrient consumption ratios also changed with the physiological status and the growth phases of phytoplankton cells. When phytoplankton growth was stimulated by iron addition, N was the first and Si was the last nutrient which absorption enhanced. The Si/N and Si/P consumption ratios of phytoplankton in the stationary and decay phases were significantly higher than those of rapidly growing phytoplankton. These findings were helpful for studies of the marine ecosystem and biogeochemistry in Prydz Bay, and were also valuable for biogeochemical studies of carbon and nutrients in various marine environments.

Application of Neomysis awatschensis as a standard marine toxicity test organism in China

The small mysid crustacean Neomysis awatschensis was collected in the west coast of Jiaozhou Bay, Qingdao, China in 1992 and acclimated and cultured in laboratory conditions since then. Standard acute toxicity tests using 4-6 d juvenile mysids of this species were conducted and the results were compared with Mysidopsis bahia, a standard toxicity test organism used in the US in terms of their sensitivities to reference toxins, as well as their taxonomy, morphology and geographic distributions. Because of its wide distribution along the Chinese coast, similar sensitivity to pollutants as M. bahia, short life history, small size and the case of handling, this study intended to use N. awatschensis as one of the standard marine organisms for toxicity testing in China. The species were applied to acute toxicity evaluations of drilling fluid and its additives I organotin TPT and toxic algae, and to chronic ( life cycle) toxicity assays of organotin TPT and a toxic dinofalgellate Alexandrium tamarense, respectively. Using N, awatschensis as a standard toxicity testing organism in marine pollution assessment in China is suggested.

Studies on the hemolytic activity of tentacle extracts of jellyfish Rhopilema esculentum Kishinouye: Application of orthogonal test

The present work is first reporting the hemolytic activity of venom from jellyfish Rhopilema esculentum Kishinouye extracted by different phosphate buffer solutions and incubated at different temperature according to the orthogonal test L6(1) x 3(6). Of the seven controllable independent variables, incubated temperature and phenylmethylsulfonyl fluoride (PMSF) had strongest effect on the hemolytic activity. (c) 2006 Elsevier B.V. All rights reserved.

Changes in plant biomass and species composition of alpine Kobresia meadows along altitudinal gradient on the Qinghai-Tibetan Plateau

Alpine Kobresia meadows are major vegetation types on the Qinghai-Tibetan Plateau. There is growing concern over their relationships among biodiversity, productivity and environments. Despite the importance of species composition, species richness, the type of different growth forms, and plant biomass structure for Kobresia meadow ecosystems, few studies have been focused on the relationship between biomass and environmental gradient in the Kobresia meadow plant communities, particularly in relation to soil moisture and edaphic gradients. We measured the plant species composition, herbaceous litter, aboveground and belowground biomass in three Kobresia meadow plant communities in Haibei Alpine Meadow Ecosystem Research Station from 2001 to 2004. Community differences in plant species composition were reflected in biomass distribution. The total biomass showed a decrease from 13196.96 +/- 719.69 g/m(2) in the sedge-dominated K. tibetica swamp to 2869.58 +/- 147.52 g/m(2) in the forb and sedge dominated K. pygmaea meadow, and to 2153.08 +/- 141.95 g/m(2) in the forbs and grasses dominated K. humilis along with the increase of altitude. The vertical distribution of belowground biomass is distinct in the three meadow communities, and the belowground biomass at the depth of 0-10 cm in K. tibetica swamp meadow was significantly higher than that in K. humilis and K. pygmaea meadows (P < 0.01). The herbaceous litter in K. tibetica swamp was significantly higher than those in K. pygnaeca and K. humilis meadows. The effects of plant litter are enhanced when ground water and soil moisture levels are raised. The relative importance of litter and vegetation may vary with soil water availability. In the K. tibetica swamp, total biomass was negatively correlated to species richness (P < 0.05); aboveground biomass was positively correlated to soil organic matter, soil moisture, and plant cover (P < 0.05); belowground biomass was positively correlated with soil moisture (P < 0.05). However, in the K. pygnaeca and K. humilis meadow communities, aboveground biomass was positively correlated to soil organic matter and soil total nitrogen (P < 0.05). This suggests that the distribution of biomass coincided with soil moisture and edaphic gradient in alpine meadows.

VTI各向异性叠前深度偏移及速度分析

Theoretical research, laboratory test and field observation show that most of sediment rock has anisotropic features. It will produce some notable errors when applying isotropic methods such as prestack depth migration and velocity analysis to dada acquired under anisotropic condition; it also has a bad effect on geologic interpretation. Generally speaking, the vertical transverse isotropic media is a good approximation to geologic structure, thus it has an important realistic meaning for anisotropic prestack depth migration theory researching and precise complex geologic imaging if considering anisotropic effect of seismic wave propagation. There are two indispensable parts in prestack depth migration of realistic records, one is proper prestack depth migration algorithm, and the other is velocity analysis using prestack seismic data. The paper consists of the two aspects. Based on implicit finite difference research proposed by Dietrich Ristow et al （1997） about VTI media prestack depth migration, the paper proposed split-step Fourier prestack depth migration algorithm （VTISSF） and Fourier finite difference algorithm （VTIFFD） based on wave equation for VTI media, program are designed and the depth migration method are tested using synthetic model. The result shows that VTISSF is a stable algorithm, it generally gets a good result if the reflector dip is not very steep, while undermigration phenomena appeared in steep dips case; the VTIFFD algorithm bring us better result in steep dips with lower efficiency and frequency dispersion. For anisotropic prestack depth migration velocity analysis of VTI media, The paper discussed the basic hypothesis of VTI model in velocity analysis algorithm, basis of anisotropic prestack depth migration velocity analysis and travel time table calculation of VTI media in integral prestack depth migration. Then , analyzed the P-wave common imaging gather in the case of homogeneous velocity and vertically variable velocity . studied the residual correction in common imaging gather produced by media parameter error, analyzed the condition of flat event and correct depth in common imaging gather . In this case, the anisotropic model parameter vector is , is vertical velocity of a point at top surface, is vertical velocity gradient, and are anisotropic parameter. We can get vertical velocity gradient from seismic data; then the P-wave common imaging gather of VTI media whose velocity varies in vertical and horizontal direction, the relationship between media parameter and event residual time shift of common image gather are studied. We got the condition of flattening common imaging gather with correct depth. In this case the anisotropic model parameter vector is , is velocity gradient in horizontal direction. As a result, the vertical velocity grads can be decided uniquely, but horizontal velocity grads and anisotropic parameter can’t be distinguished if no priori information available, our method is to supply parameter by velocity scanning; then, as soon as is supplied we can get another four parameters of VTI media from seismic data. Based on above analysis, the paper discussed the feasibility of migration velocity analysis in vertically and horizontally varied VTI media, synthetic record of three models are used to test the velocity analysis method . Firstly, anisotropic velocity analysis test is done using a simple model with one block, then we used a model with multiple blocks, thirdly, we analyzed the anisotropic velocity using a part of Marmousi model. The model results show that this velocity analysis method is feasible and correct.

地下水流动数值模拟的并行计算方法研究

Numerical modeling of groundwater is very important for understanding groundwater flow and solving hydrogeological problem. Today, groundwater studies require massive model cells and high calculation accuracy, which are beyond single-CPU computer’s capabilities. With the development of high performance parallel computing technologies, application of parallel computing method on numerical modeling of groundwater flow becomes necessary and important. Using parallel computing can improve the ability to resolve various hydro-geological and environmental problems. In this study, parallel computing method on two main types of modern parallel computer architecture, shared memory parallel systems and distributed shared memory parallel systems, are discussed. OpenMP and MPI (PETSc) are both used to parallelize the most widely used groundwater simulator, MODFLOW. Two parallel solvers, P-PCG and P-MODFLOW, were developed for MODFLOW. The parallelized MODFLOW was used to simulate regional groundwater flow in Beishan, Gansu Province, which is a potential high-level radioactive waste geological disposal area in China. 1. The OpenMP programming paradigm was used to parallelize the PCG (preconditioned conjugate-gradient method) solver, which is one of the main solver for MODFLOW. The parallel PCG solver, P-PCG, is verified using an 8-processor computer. Both the impact of compilers and different model domain sizes were considered in the numerical experiments. The largest test model has 1000 columns, 1000 rows and 1000 layers. Based on the timing results, execution times using the P-PCG solver are typically about 1.40 to 5.31 times faster than those using the serial one. In addition, the simulation results are the exact same as the original PCG solver, because the majority of serial codes were not changed. It is worth noting that this parallelizing approach reduces cost in terms of software maintenance because only a single source PCG solver code needs to be maintained in the MODFLOW source tree. 2. P-MODFLOW, a domain decomposition–based model implemented in a parallel computing environment is developed, which allows efficient simulation of a regional-scale groundwater flow. The basic approach partitions a large model domain into any number of sub-domains. Parallel processors are used to solve the model equations within each sub-domain. The use of domain decomposition method to achieve the MODFLOW program distributed shared memory parallel computing system will process the application of MODFLOW be extended to the fleet of the most popular systems, so that a large-scale simulation could take full advantage of hundreds or even thousands parallel processors. P-MODFLOW has a good parallel performance, with the maximum speedup of 18.32 (14 processors). Super linear speedups have been achieved in the parallel tests, indicating the efficiency and scalability of the code. Parallel program design, load balancing and full use of the PETSc were considered to achieve a highly efficient parallel program. 3. The characterization of regional ground water flow system is very important for high-level radioactive waste geological disposal. The Beishan area, located in northwestern Gansu Province, China, is selected as a potential site for disposal repository. The area includes about 80000 km2 and has complicated hydrogeological conditions, which greatly increase the computational effort of regional ground water flow models. In order to reduce computing time, parallel computing scheme was applied to regional ground water flow modeling. Models with over 10 million cells were used to simulate how the faults and different recharge conditions impact regional ground water flow pattern. The results of this study provide regional ground water flow information for the site characterization of the potential high-level radioactive waste disposal.

地震P-S波时差耦合作用的斜坡崩滑效应研究

China locates between the circum-Pacific and the Mediterranean-Himalayan seismic belt. The seismic activities in our country are very frequent and so are the collapses and slides of slope triggered by earthquakes. Many collapses and slides of slope take place mainly in the west of China with many earthquakes and mountains, especially in Sichuan and Yunnan Provinces. When a strong earthquake happening, the damage especially in mountains area caused by geological hazards it triggered such as rock collapses, landslides and debris flows is heavier than that it caused directly. A conclusion which the number of lives lost caused by geological hazards triggered by a strong earthquake in mountains area often accounts for a half even more of the total one induced by the strong earthquake can be made by consulting the statistical loss of several representative earthquakes. As a result, geological hazards such as collapses and slides of slope triggered by strong earthquakes attract wide attention for their great costs. Based on field geological investigation, engineering geological exploration and material data analysis, chief conclusions have been drawn after systematic research on formation mechanism, key inducing factors, dynamic characteristics of geological hazards such as collapses and slides of slope triggered by strong earthquakes by means of engineering geomechanics comprehensive analysis, finite difference numerical simulation test, in-lab dynamic triaxial shear test of rock, discrete element numerical simulation. Based on research on a great number of collapses and landslides triggered by Wenchuan and Xiaonanhai Earthquake, two-set methods, i.e. the method for original topography recovering based on factors such as lithology and elevation comparing and the method for reconstructing collapsing and sliding process of slope based on characteristics of seism tectonic zone, structural fissure, diameter spatial distribution of slope debris mass, propagation direction and mechanical property of seismic wave, have been gotten. What is more, types, formation mechanism and dynamic characteristics of collapses and slides of slope induced by strong earthquakes are discussed comprehensively. Firstly, collapsed and slided accumulative mass is in a state of heavily even more broken. Secondly, dynamic process of slope collapsing and sliding consists of almost four stages, i.e. broken, thrown, crushed and river blocked. Thirdly, classified according to failure forms, there are usually four types which are made up of collapsing, land sliding, land sliding-debris flowing and vibrating liquefaction. Finally, as for key inducing factors in slope collapsing and sliding, they often include characteristics of seism tectonic belts, structure and construction of rock mass, terrain and physiognomy, weathering degree of rock mass and mechanical functions of seismic waves. Based on microscopic study on initial fracturing of slope caused by seismic effect, combined with two change trends which include ratio of vertical vs. horizontal peak ground acceleration corresponding to epicentral distance and enlarging effect of peak ground acceleration along slope, key inducing factor of initial slope fracturing in various area with different epicentral distance is obtained. In near-field area, i.e. epicentral distance being less than 30 km, tensile strength of rock mass is a key intrinsic factor inducing initial fracturing of slope undergoing seismic effect whereas shear strength of rock mass is the one when epicentral distance is more than 30 km. In the latter circumstance, research by means of finite difference numerical simulation test and in-lab dynamic triaxial shear test of rock shows that initial fracture begins always in the place of slope shoulder. The fact that fracture strain and shear strength which are proportional to buried depth of rock mass in the place of slope shoulder are less than other place and peak ground acceleration is enlarged in the place causes prior failure at slope shoulder. Key extrinsic factors inducing dynamic fracture of slope at different distances to epicenter have been obtained through discrete element numerical simulation on the total process of collapsing and sliding of slope triggered by Wenchuan Earthquake. Research shows that combined action of P and S seismic waves is the key factor inducing collapsing and sliding of slope at a distance less than 64 km to initial epicenter along earthquake-triggering structure. What is more, vertical tensile action of P seismic wave plays a leading role near epicenter, whereas vertical shear action of S seismic wave plays a leading role gradually with epicentral distance increasing in this range. On the other hand, single action of P seismic wave becomes the key factor inducing collapsing and sliding of slope at a distance between 64 km and 216 km to initial epicenter. Horizontal tensile action of P seismic wave becomes the key factor gradually from combined action between vertical and horizontal tensile action of P seismic wave with epicentral distance increasing in this distance range. In addition, initial failure triggered by strong earthquakes begins almost in the place of slope shoulder. However, initial failure beginning from toe of slope relates probably with gradient and rock occurrence. Finally, starting time of initial failure in slope increases usually with epicentral distance. It is perhaps that the starting time increasing is a result of attenuating of seismic wave from epicenter along earthquake-triggering structure. It is of great theoretical and practical significance for us to construct towns and infrastructure in fragile geological environment along seism tectonic belts and conduct risk management on earthquake-triggered geological hazards by referring to above conclusions.