Spatial and temporal variations of sound speed at the PN section

Gridded sound speed data were calculated using Del Grosso's formulation from the temperature and salinity data at the PN section in the East China Sea covering 92 cruises between February 1978 and October 2000. The vertical gradients of sound speed are mainly related to the seasonal variations, and the strong horizontal gradients are mainly related to the Kuroshio and the upwelling. The standard deviations show that great variations of sound speed exist in the upper layer and in the slope zone. Empirical orthogonal function analysis shows that contributions of surface heating and the Kuroshio to sound speed variance are almost equivalent.

胶州湾及其邻近海域的浮游植物－种类组成与数量变化

研究胶州湾浮游植物的物种组成与时空分布的特点，对于了解该湾生态系统的现状与历史变化趋势以及生态系统对自然条件变化和人类干扰的响应具有重要意义。 本文根据2004年每月一次采集的浮游植物样品，分析了胶州湾浮游植物的物种组成和优势种的时空分布情况。调查发现浮游植物共142种，分属于6门53属，其中硅藻门40属113种，占总物种数的79.6%；甲藻门10属24种，占总物种数的16.9%；其它为金藻门1属2种，裸藻1属1种，绿藻1属2种。浮游植物丰度周年波动范围为11.12－14602.39×104cells/m3，全年平均为1857.55×104cells/m3。全年丰度最高的藻为环纹劳德藻(Lauderia annulata)，出现在2月份，而全年的丰度最高值也出现在2月份。在胶州湾中，硅藻所占比例最大，平均为97.44%，最高为99.95%，最低为82.55%。浮游植物的丰度的周年变化：有两个高峰，分别出现在2月份和10月份，是典型的温带海域双高峰的分布形式。 胶州湾浮游植物的优势种分析，采用Kikvidze等（2002）提出的确定优势种数量的计算方法，再依据各物种的优势度排序最终确定优势种。其中在多样性最低的6月和10月，优势种数目均为1种，分别为丹麦细柱藻（Leptocylindrus danicus）和中肋骨条藻（Skeletonema costatum）。而在多样性指数较低的1月、2月，优势种数目分别为2和3；而在物种多样性比较高的5月、7月、11月，由于分布相对比较分散，所以优势种数量较多。这说明这个方法可以较好地完成对胶州湾浮游植物群落分析时确定优势种的目标。优势种出现频率较多的种类为中肋骨条藻（1、3、4、5、8、10、11、12月）、洛氏角毛藻(Chaetoceros lorenzianus)（7、8、9、11、12月）、尖刺拟菱形藻(Pseudonitzschia pungens)（2、3、5、7月）、加拉星杆藻(Asteronella kariana)（1、2、3、5月）、密联角毛藻(Chaetoceros densus)（5、11、12月）、扭鞘藻(Streptothece thamesis)（7、8、9月）、夜光藻(Noctiluca scintillancs)（4、5、7月）、奇异菱形藻(Nitzschia paradoxa)（5、11、12月）,对这些种的时空分布进行了分析。 对2004年胶州湾的浮游植物数据进行分层聚类分析后，发现可以将胶州湾划分为三个海区：湾南与湾外海区、湾中西部海区、湾东海区。

Spatial and temporal variations in pH and total alkalinity at the beginning of the rainy season in the Changjiang Estuary, China

The results of field observation carried out in May 2003 were used to examine pH and total alkalinity behaviors in the Changjiang Estuary. It was showed that PH and total alkalinity took on clear spatial variations in values with the minima in the low salinity region. Like salinity, transect distributions of PH and total alkalinity (TA) in a downriver direction had a sharp gradient each. These gradients appeared in such a sequence that the TA gradient was earlier than salinity and PH gradients, and the salinity gradient was earlier than the PH gradient. These distribution characteristics seemed to be strongly influenced by the mixing process of freshwater and seawater, for both PH and total alkalinity had significant linear relationships with salinity and temperature. For PH, phytoplankton activities also had a significant impact upon its spatial distribution. During a period of 48 h, PH and total alkalinity changed within wide ranges for every layer of the two anchor stations, namely, Stas 13 and 20, which were located at the mixed water mass and seawater mass, respectively. For both Stas 13 and 20, PH and TA fluctuation of every layer could be very wide during a 4 h period. As a whole, the data of the two anchor stations showed that neither variations in salinity and temperature nor phytoplankton activities were the main factors strongly influencing the total alkalinity temporal variability on a small time scale. The data of Sta. 20 implied that both salinity variation and phytoplankton activities had a significant influence on PH temporal variability, but the same conclusion could not be drawn from the data of Sta. 13.

Investigation of heavy metals in sediments and Manila clams Ruditapes philippinarum from Jiaozhou Bay, China

Manila clams (Ruditapes philippinarum) and sediments were collected bimonthly during 2007 at five locations in Jiaozhou Bay near Qingdao, China, to determine heavy metal concentrations and to assess the validation of R. philippinarum as a metal biomonitor. Concentrations of heavy metals in clam soft tissues ranged between 0.75 and 3.31, 0.89 and 15.20, 5.70 and 26.03, 52.12 and 110.33, 10.30 and 72.34, 9.64 and 28.60, and 3.15 and 52.75 mu g g (-aEuro parts per thousand 1) dry weight for Cd, Pb, Cu, Zn, Mn, Cr, and Ni, respectively. Most of the highest values occurred at the northeast bay and the lowest values occurred at the western part. Regarding seasonal variation, relatively high tissue metal concentrations were observed during October or December. A similar pattern was also found in habitat sediments. There was a strong correlation between the concentrations of Cd, Pb, Zn, Mn, Cr, and Ni in soft tissues and surrounding sediments. It is indicated that R. philippinarum could be used as a biomonitor for heavy metal contamination in Jiaozhou Bay.

Temporal and spatial variability in nitrogen uptake kinetics during harmful dinoflagellate blooms in the East China Sea

Temporal and spatial variability in the kinetic parameters of uptake of nitrate (NO3-), ammonium (NH4+), urea, and glycine was measured during dinoflagellate blooms in Changjiang River estuary and East China Sea coast, 2005. Karenia mikimotoi was the dominant species in the early stage of the blooms and was succeeded by Prorocentrum donghaiense. The uptake of nitrogen (N) was determined using N-15 tracer techniques. The results of comparison kinetic parameters with ambient nutrients confirmed that different N forms were preferentially taken up during different stages of the bloom. NO3- (V-max 0.044 h(-1); K-s 60.8 mu M-N) was an important N source before it was depleted. NH4+ (V-max 0.049 h(-1); K-s 2.15 mu M-N) was generally the preferred N. Between the 2 organic N sources, urea was more preferred when K. mikimotoi dominated the bloom (V-max 0.020 h(-1); K-s 1.35 mu M-N) and glycine, considered as a dominant amino acid, was more preferred when P. donghaiense dominated the bloom (V-max 0.025 h(-1); K-s 1.76 mu M-N). The change of N uptake preference by the bloom-forming algae was also related to the variation in ambient N concentrations. Published by Elsevier B.V.

Spatial and temporal relationships between precipitation and ANPP of four types of grasslands in northern China

Precipitation is considered to be the primary resource limiting terrestrial biological activity in water-limited regions. Its overriding effect on the production of grassland is complex. In this paper, field data of 48 sites (including temperate meadow steppe,temperate steppe, temperate desert steppe and alpine meadow) were gathered from 31 published papers and monographs to analyze the relationship between above-ground net primary productivity (ANPP) and precipitation by the method of regression analysis. The results indicated that there was a great difference between spatial pattern and temporal pattern by which precipitation influenced grassland ANPP. Mean annual precipitation (MAP) was the main factor determining spatial distribution of grassland ANPP (r~2 = 0.61,P ＜ 0.01); while temporally, no significant relationship was found between the variance of AN PP and inter-annual precipitation for the four types of grassland. However, after dividing annual precipitation into monthly value and taking time lag effect into account, the study found significant relationships between ANPP and precipitation. For the temperate meadow steppe, the key variable determining inter-annual change of ANPP was last August-May precipitation (r~2= 0.47, P = 0.01); for the temperate steppe, the key variable was July precipitation (r~2 = 0.36, P = 0.02); for the temperate desert steppe, the key variable was April-June precipitation (r~2 = 0.51, P ＜0.01); for the alpine meadow, the key variable was last September-May precipitation (r~2 = 0.29, P ＜ 0.05). In comparison with analogous research, the study demonstrated that the key factor determining inter-annual changes of grassland ANPP was the cumulative precipitation in certain periods of that year or the previous year.

Breeding behavior under temporal risk of predation in male root voles (Microtus oeconomus)

We examined breeding behavior responses of male root votes (Microtus oeconomus) to temporal risk of predation by using acute and chronic exposure to predator odor. The 2 series of exposure experiments provided 2 types of temporal patterns of risk: continuous safety with a brief period of risk and Sustained risk with a brief period of safety. Male root votes that were acutely exposed to predator odor for I h suppressed their breeding behavior, but bred immediately after exposure to control odor for I h. Those chronically exposed to predator odor for 20 days maintained behavioral suppression during the 1-h period of exposure to control odor. Acutely exposed males did not change their physiological patterns of breeding, but those chronically exposed to predator odor had reduced testosterone concentration and epididymis index. Our results indicate that breeding behavior in a given situation depends on the overall patterns of risk experienced by male root votes, and the acute and chronic stress responses that affect reproduction are responsible for different behavioral responses to the 2 types of temporal patterns of risk. We also discuss the reasons for conflicting results about breeding suppression of votes between previous studies in the laboratory and the field.

裂缝储层预测技术应用研究

In practice,many fracture reservoir was found,which has giant potential for exporation. For example,in limestone fracture reservoir,igneous rock fracture reservoir and shale fracture reservoir ,there are high yield oil wells found. The fracture reservoir has strong anisotropy and is very difficult to explore and produce.Since 1990’s,the techniques that use structure information and P-WAVE seismic attributes to detect fracture developed very rapidly,include stress and strain analysis,using amplitude,interval velocity,time-difference,azimuthal AVO analysis etc. Based on research and develop these advanced techniques of fracture detect,this paper selected two typical fracture reservoir as target area,according to the characters of research area,selected different techniques to pridect the fracture azimuth and density of target,and at last ,confirmed the favored area. This paper includes six parts:the first chapter mainly addresses the domestic and international research actuality about the fracture prediction and the evolement in ShengLi oil field,then according to the temporal exploration requirement,a research route was established; Based on the close relationship between structural fracture and the geotectonic movement and the procedure of rock distortion，the second chapter research the structural fracture predicting technique which is realized by computing the strain in every geotectonic movement ,which is by use of the forward and inversion of the growing history of structure; The third chapter discussed many kind of traditional techniques for fracture reservoir prediction,and point out their disadvantages.then research and develop the coherence volume computing technique which can distinguish from faults,the seismic wave absorbing technique,and other fracture predicting technique which is by use of seismic attributes ,such as azimuthal AVO FVO etc; The fourth chapter first establish the geological and petrophysical model by use of the existed log and drill well information, then research the variation of amplitude and seismic wave which is caused by fractures.based on it , the fracture predicting technique which is by use of variation of azimuthal impedance is researched;The fifth chapter is a case study,it selects shale fracture reservoir in LuoJia area as target,selects several kind of techniques to apply ,at last ,the fracture distribution of target reservoir and favored area were gotten;the sixth chapter is another case study,it selects limestone fracture reservoir in BoShen6 buried hill as target,selects several kind of techniques to apply,similarly favored area were gotten. Based on deeply research and development of the new techniques for fracture reservoir exploration, This paper selects two fracture reservoirs the most typical in ShengLi as targets to be applied ,good results show up a good application way ,which can be used for reference for future fracture exploration,and it can bring materially economic and social benefit.

中低纬度变化磁场的暴时特征研究

Magnetic storm is a kind of severe disturbances in the whole solar-earth electromagnetic space. It has significant effects on communication, electric power, oil transport pipe and human activities in space. Therefore, magnetic storms are worth for applications systems, not only being a favorable issue for scientists. In this paper, the spatial and temporal distributions of the magnetic fields produced by the magnetosphere-ionosphere current systems during storms are studied. Four parts are included in this paper decomposion of different disturbances with different origins, topological structure of the ring current, the asymmetric characteristics of the ring current, and the statistic peculiarities of the day-to-day variability (DTD) of Sq. 1 The decomposition of magnetic disturbances at mid-low latitudes and its evolutions during storms Transient variations in the geomagnetic field recorded at mid-low latitudes mainly include the storm-time variation (Dst), solar quiet daily variation (Sq) and disturbance daily variation (SD). With the data of the geomagnetic meridian chain observatories in China, 25 storms during the period of 1997 to 1999 have been analyzed. According to the features of different variations, a method of “three-steps decomposition” is developed by using the method of Natural Orthogonal Components (NOC), Correlation Analysis and Fourier Analysis to separate those three components in turn. The results show that, the first eigenmode by the MNOC clearly describing the special distribution and temporal evolution of storm-time variation, in addition, Correlation Analysis and Fourier Analysis offer a useful method to extract the Sq and SD variations. The latitudinal shift of the Sq current focus seems to be the principal reason of the day-to-day variaitons in the daily range of Sq. The magnitude of SD reaches a maximum during the main phase, and then gradually decreases. 2 The topology structure of the ring current during storms Both the mechanism of the ring current and the geomagnetic data suggest that the central plane of the ring current is declining to the geomagnetic equator plane with a tilt angle δ. Using the H and Z component data at two stations in a meridian chain, we deduce a new parameter describing the invariable peculiarity of different storms. Then the δ angle is calculated by using the data from a meridian chain and tested with the ERC model. Finally the deduced tilt angles are used to modify Dst index. 3 The asymmetric characteristics of the ring current during storms The variations of the geomagnetic field at mid-low latitudes show a significant dawn-dusk asymmetry, resulting from the superposition of the fields from the symmetric ring current and the partial ring current. On the basis of the data from the 20°E, 30°E meridian chains and 30°N latitudinal chain, the dawn-dusk asymmetry is investigated by using three methods, namely, statistic analysis, ring current model calculation and typical event analysis. This characteristic implies the asymmetry of the spatial distribution of the ring current. In addition, during the main phase after the sudden commencement (SC), H field increases and reaches maximum around noontime, implying the effect of the Chapman-Ferraro current. 4 The statistic characteristics of the day-to-day variability and its mechanism The day-to-day variability of the geomagnetic Sq field is studied by using the magnetic data from a meridian chain of magnetometers along 120° E longitude. The method of NOC is applied to separate the Sq variation from complicated disturbances. The first eigenmode with the largest eigenvalue represents fairly well the Sq variation with a conspicuous day-to-day variability in the daily range. For the stations on the same north- or south-side of the Sq current system focus, the day-to-day variations show a positive correlation. In contrast, for the stations on the different sides of the Sq focus, they show a negative correlation, suggesting an important role of latitudinal shift of the Sq current system focus to the day-to-day variability of the Sq daily range. The Sq daily range is correlated with the magnetic indices Ap and Dst in a peculiar way: on some severe disturbed days, noticeably enhancements of the Sq are observed, implying increases of the ionospheric conductivities and/or tidal wind velocities; on other severe disturbed days, however, dramatically reduced Sq variations occur, suggesting dominant effects of the ‘disturbance dynamo’ process.

巴丹吉林沙漠近40a风沙活动规律浅析

Based on the data of wind from 11 meteorological observation stations in recent 40 years, the paper analyzed the temporal and spatial distribution characteristics of resultant drift potential in the Badain Jaran Desert. The 84 sand samples were collected and analysed from three compound crescentic dunes developped on southeast part of The Badain Jaran Desert. Through the statistics and computation of wind data and the analysis of Grain Size over compound crescentic dunes in The Badain Jaran Desert, this paper analyzed the characteristics of the blown sand activities in in this area. The following viewpoints were obtained: 1　Most parts of The Badain Jaran Desert are dominated by northwest wind and sub－dominated by west wind, while prevailing wind in southeast part of The Badain Jaran Desert is southeast, different froms other places. 2 The Badain Jaran Desert is strongly zoned from its surounding meteorological observation stations. It has impressively high DP and RDP value. Except for Ding Xin Station in medium wind power environment, the left stations are all belongs to high wind－power districts. 3 From the seventies till now, the value of DP and RDP in The Badain Jaran Desert shows a decending trend except for the RDP value of Guai Zihu Station in the north. This trend indicates that the wind power for the moment of The Badain Jaran Desert is far lower than it in the seventies. 4 Adopting sand rose map for classification of wind environments in defferent parts of The Badain Jaran Desert, Ejinaqi is Narrow unimodal wind regimes, Ding xin, Bayinmaodao and Ayouqi districts obtuse bimodal and Guaizihu complex. 5 Resultant drift directions in most parts of the Badain Jaran desert are southeast, which indicates annual net sand in these areas drifted from the northwest to the southeast of the desert, which consists with the analysis results of grain size of predecessors. RDD in Ayouqi staion shows the net sand drifted from the southeast to the northwest in the recent 40 years, which may attribute to the reason that, the Badain Jaran desert lies in the zone influenced by both northwest and southeast wind, the southeast monsoon was strengthened and northwest wind became weaker than a former time. So, if sand mareria is sufficient, there are two different drift directions, namely two sorts of sand source, existed in Badain Jaran desert at least. 6 Based on analysis results of meteorologic data and grain size, blown sand activities is that, in most parts of Badain Jaran desert in the past 40 years, northwest wind transported net sand from the northwest part to the southeast part of The Badain Jaran Desert.And besides this, analysis of collected sand samples also shows that, in the southeast part of Badain Jaran desert, grain size do fined from southeast to northwest. Combining results of meteorologic data in the southeast part, it can infered that southeast monsoon took blown sand to the inner desert from its southeast part.

特大地震激发地球自由振荡的数值模拟研究

Large earthquakes, such as the Chile earthquake in 1960 and the Sumatra-Andaman earthquake on Dec 26, 2004 in Indonesia, have generated the Earth’s free oscillations. The eigenfrequencies of the Earth’s free oscillations are closely related to the Earth’s internal structures. The conventional methods, which mainly focus on calculating the eigenfrequecies by analytical ways, and the analysis on observations can not easily study the whole processes from earthquake occurrence to the Earth’s free oscillation inspired. Therefore, we try to use numerical method incorporated with large-scale parallel computing to study on the Earth’s free oscillations excited by giant earthquakes. We first give a review of researches and developments of the Earth’s free oscillation, and basical theories under spherical coordinate system. We then give a review of the numerical simulation of seismic wave propagation and basical theories of spectral element method to simulate global seismic wave propagation. As a first step to study the Earth’s free oscillations, we use a finite element method to simulate the propagation of elastic waves and the generation of oscillations of the chime bell of Marquis Yi of Zeng, by striking different parts of the bell, which possesses the oval crosssection. The bronze chime bells of Marquis Yi of Zeng are precious cultural relics of China. The bells have a two-tone acoustic characteristic, i.e., striking different parts of the bell generates different tones. By analysis of the vibration in the bell and the spectrum analysis, we further help the understanding of the mechanism of two-tone acoustic characteristics of the chime bell of Marquis Yi of Zeng. The preliminary calculations have clearly shown that two different modes of oscillation can be generated by striking different parts of the bell, and indicate that finite element numerical simulation of the processes of wave propagation and two-tone generation of the chime bell of Marquis Yi of Zeng is feasible. These analyses provide a new quantitative and visual way to explain the mystery of the two-tone acoustic characteristics. The method suggested by this study can be applied to simulate free oscillations excited by great earthquakes with complex Earth structure. Taking into account of such large-scale structure of the Earth, small-scale low-precision numerical simulation can not simply meet the requirement. The increasing capacity in high-performance parallel computing and progress on fully numerical solutions for seismic wave fields in realistic three-dimensional spherical models, Spectral element method and high-performance parallel computing were incorporated to simulate the seismic wave propagation processes in the Earth’s interior, without the effects of the Earth’s gravitational potential. The numerical simulation shows that, the results of the toroidal modes of our calculation agree well with the theoretical values, although the accuracy of our results is much limited, the calculated peaks are little distorted due to three-dimensional effects. There exist much great differences between our calculated values of spheroidal modes and theoretical values, because we don’t consider the effect the Earth’ gravitation in numerical model, which leads our values are smaller than the theoretical values. When , is much smaller, the effect of the Earth’s gravitation make the periods of spheroidal modes become shorter. However, we now can not consider effects of the Earth’s gravitational potential into the numerical model to simulate the spheroidal oscillations, but those results still demonstrate that, the numerical simulation of the Earth’s free oscillation is very feasible. We make the numerical simulation on processes of the Earth’s free oscillations under spherically symmetric Earth model using different special source mechanisms. The results quantitatively show that Earth’s free oscillations excited by different earthquakes are different, and oscillations at different locations are different for free oscillation excited by the same earthquake. We also explore how the Earth’s medium attenuation will take effects on the Earth’s free oscillations, and take comparisons with the observations. The medium attenuation can make influences on the Earth’s free oscillations, though the effects on lower-frequency fundamental oscillations are weak. At last, taking 2008 Wenchuan earthquake for example, we employ spectral element method incorporated with large-scale parallel computing technology to investigate the characteristics of seismic wave propagation excited by Wenchuan earthquake. We calculate synthetic seismograms with one-point source model and three-point source model respectively. Full 3-D visualization of the numerical results displays the profile of the seismic wave propagation with respect to time. The three-point source, which was proposed by the latest investigations through field observation and reverse estimation, can better demonstrate the spatial and temporal characteristics of the source rupture processes than one-point source. Primary results show that those synthetic signals calculated from three-point source agree well with the observations. This can further reveal that the source rupturing process of Wenchuan earthquake is a multi-rupture process, which is composed by at least three or more stages of rupture processes. In conclusion, the numerical simulation can not only solve some problems concluding the Earth’s ellipticity and anisotropy, which can be easily solved by conventional methods, but also finally solve the problems concluding topography model and lateral heterogeneity. We will try to find a way to fully implement self-gravitation in spectral element method in future, and do our best to continue researching the Earth’s free oscillations using the numerical simulations to see how the Earth’ lateral heterogeneous will affect the Earth’s free oscillations. These will make it possible to bring modal spectral data increasingly to bear on furthering our understanding of the Earth’s three-dimensional structure.

乡镇级人口统计数据空间化研究

Population data which collected and saved according to administrative region is a kind of statistical data. As a traditional method of spatial data expression, average distribution in every administrative region brings population data on a low spatial and temporal precision. Now, an accurate population data with high spatial resolution is becoming more and more important in regional planning, environment protection, policy making and rural-urban development. Spatial distribution of population data is becoming more important in GIS study area. In this article, the author reviewed the progress of research on spatial distribution of population. Under the support of GIS, correlative geographical theories and Grid data model, Remote Sensing data, terrain data, traffic data, river data, resident data, and social economic statistic were applied to calculate the spatial distribution of population in Fujian province, which includes following parts: (1) Simulating of boundary at township level. Based on access cost index, land use data, traffic data, river data, DEM, and correlative social economic statistic data, the access cost surface in study area was constructed. Supported by the lowest cost path query and weighted Voronoi diagram, DVT model (Demarcation of Villages and Towns) was established to simulate the boundary at township level in Fujian province. (2) Modeling of population spatial distribution. Based on the knowledge in geography, seven impact factors, such as land use, altitude, slope, residential area, railway, road, and river were chosen as the parameters in this study. Under the support of GIS, the relations of population distribution to these impact factors were analyzed quantificationally, and the coefficients of population density on pixel scale were calculated. Last, the model of population spatial distribution at township level was established through multiplicative fusion of population density coefficients and simulated boundary of towns. (3) Error test and analysis of population spatial distribution base on modeling. The author not only analyzed the numerical character of modeling error, but also its spatial distribution. The reasons of error were discussed.