非构造形变对GPS连续站位置时间序列的影响和修正

GPS观测得到的地壳形变场通常包含有构造形变与非构造形变二类信息, 去除其中的非构造形变信息对于有效运用GPS数据研究构造形变场至关重要. 本文运用国际卫星对地观测资料及各类地球物理模型, 定量计算海潮、大气、积雪和土壤水、海洋非潮汐4项负荷效应造成的地壳非构造形变, 并以此研究和修正这些非构造形变对中国地壳运动观测网络GPS基准站位置时间序列的影响. 研究发现此4项负荷效应, 特别是大气、积雪和土壤水, 对于测站垂向位置的影响显著. 通过模型改正可以使测站垂向位置的RMS降低～1 mm, 占其总量的～11％. 对于垂向时间序列的周年项部分, 这一改正可降低其振幅的37％. 研究还表明经过地球物理模型改正和周年、半周年谐波拟合改正的时间序列比起仅经过周年、半周年谐波拟合改正的时间序列更为平滑, 表明地球物理模型改正对于消除非构造形变场的作用不是周年、半周年谐波拟合改正所能替代的.     

Spatial and temporal variation of picoplankton distribution in the Yellow Sea, China

Seven surveys were carried out in April, September, October, December 2006 and March, May, August 2007 in the Yellow Sea, China. Variations in the spatial and temporal distribution of Synechococcus, picoeukaryotes and heterotrophic bacteria are quantified using flow cytometry. Synechococcus and heterotrophic bacteria are most abundant from late spring to autumn, while picoeukaryotes concentration is high in spring. Synechococcus and heterotrophic bacteria concentrated high in the northwest part of the Yellow Sea in spring and autumn, while picoeukaryotes distributed evenly over the whole study area except for a small frontal zone in the coastal area on the west (in spring) and central Yellow Sea (in autumn). Under mixing conditions, the vertical distribution of the three picoplankton groups showed a well-mixed pattern. Upon a well-established stratification, the maximum abundance of picoplankton occurred above the mixed layer depth (~30 m). Cell sizes of Synechococcus and picoeukaryotes were estimated by converting forward scatter signals (FSC) from cytometry analysis to cell diameter, showing the results of 0.65–0.82 μm for Synechococcus and 0.85–1.08 μm for picoeukaryotes. The average integrated carbon biomasses ranged 15.26–312.62 mgC/m2 for Synechococcus, 18.54–61.57 mgC/m2 for picoeukaryotes, and 402.63–818.46 mgC/m2 for heterotrophic bacteria. The distribution of Synechococcus and heterotrophic bacteria was temperature dependent, and picoplankton presence was poor in the Yellow Sea Cold Water Mass.     

A comparison of different mass elements for use in gravity gradiometry

Topographic and isostatic mass anomalies affect the external gravity field of the Earth. Therefore, these effects also exist in the gravity gradients observed, e.g., by the satellite gravity gradiometry mission GOCE (Gravity and Steady-State Ocean Circulation Experiment). The downward continuation of the gravitational signals is rather difficult because of the high-frequency behaviour of the combined topographic and isostatic effects. Thus, it is preferable to smooth the gravity field by some topographic-isostatic reduction. In this paper the focus is on the modelling of masses in the space domain, which can be subdivided into different mass elements and evaluated with analytical, semi-analytical and numerical methods. Five alternative mass elements are reviewed and discussed: the tesseroid, the point mass, the prism, the mass layer and the mass line. The formulae for the potential, the attraction components and the Marussi tensor of second-order potential derivatives are provided. The formulae for different mass elements and computation methods are checked by assuming a synthetic topography of constant height over a spherical cap and the position of the computation point on the polar axis. For this special situation an exact analytical solution for the tesseroid exists and a comparison between the analytical solution of a spherical cap and the modelling of different mass elements is possible. A comparison of the computation times shows that modelling by tesseroids with different methods produces the most accurate results in an acceptable computation time. As a numerical example, the Marussi tensor of the topographic effect is computed globally using tesseroids calculated by Gauss–Legendre cubature (3D) on the basis of a digital height model. The order of magnitude in the radial-radial component is about  ± 8 E.U. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Activation energy for diffusion of helium in water-saturated, compacted Na-montmorillonite

It is important to clarify the migration behavior of hydrogen gas dissolved in water-saturated, compacted bentonite, which is a promising material for geologic disposal of high-level waste and TRU waste disposal. The diffusion coefficients of helium, which can be detected under extremely low background conditions, in water-saturated, compacted Na-montmorillonite were determined as a function of temperature by a transient diffusion method. The activation energies for diffusion of helium were then obtained. The activation energies were from 6.9 ± 4.8 to 19 ± 2.8 kJ mol^− 1 and were regarded to be independent of dry density. The activation energies of helium in water-saturated Na-montmorillonite were roughly equal to those in bulk water, 14.9 kJ mol^− 1, and in ice, from 11 to 13 kJ mol^− 1. It is possible that helium diffuses not only in pore water but also in interlayer water.     

基于BIF条带状铁建造矿化特征的国际铁矿石盈亏平衡运营成本定量研究

国内外铁矿石价格对标基准多采用离岸价或到岸价,而非盈亏平衡运营成本,难以揭示我国铁矿石所面对的真实市场承压价格。为了厘清国际一线生产商的铁矿石盈亏平衡运营成本价格,本文对世界上最重要的条带状铁建造(BIF)矿产地——西澳哈默斯利盆地高品位赤铁矿矿床的矿化特征及代表性铁矿石产品展开系统研究,同时引入巴西铁四角地区的铁英岩型赤铁矿矿石作为对照,分析全球典型高品位赤铁矿矿石经济指标。结合前人研究成果,将西澳哈默斯利盆地与BIF相关的高品位赤铁矿的富集矿化类型划分为假象赤铁矿-针铁矿、微板状赤铁矿与河道沉积型赤铁矿,巴西铁四角主要为铁英岩型赤铁矿。上述各矿化类型对应的铁矿石产品的铁元素含量均高于56%;在杂质元素含量上,假象赤铁矿-针铁矿的磷含量高,微板状赤铁矿的磷、硫含量较高,河道沉积型赤铁矿的磷、硫含量较低,铁英岩型赤铁矿含锰。经定量估算,西澳力拓、必和必拓、FMG和巴西淡水河谷的铁矿石盈亏平衡运营成本价格分别为34.66、36.76、47.35、38.07美元/干吨,可为中国海外权益铁矿项目开发提供运营成本的参考。     

Behaviour of Fe-oxides relevant to contaminant uptake in the environment

The behaviour of Fe-oxides was investigated during precipitation and co-precipitation, phase transformation and dissolution, while their ability to adsorb and incorporate trace components was examined. Some samples were synthesised and studied under controlled laboratory conditions and other samples were taken from experiments designed to test the effectiveness of waste treatment strategies using iron. Surface-sensitive and high-resolution techniques were used to complement information gathered from classical, macroscopic methods.

Adsorption isotherms for Ni²⁺ uptake on synthetic ferrihydrite (Fe₅HO₈·4H₂O, often written simply Fe(OH)₃), goethite (-FeOOH), hematite (-Fe₂O₃) and magnetite (Fe₃O₄) were all similar, increasing as expected at higher pH. Desorption behaviour was also similar, but one third or more of the Ni²⁺ failed to return to solution. In the past, “irreversible sorption” has been blamed on uptake into micro-fractures or pores, but during examination (using Atomic force microscopy, AFM) of hundreds of Fe-oxide particles, no evidence for such features could be found, leading to the conclusion that Ni²⁺ must become incorporated onto or into the solids. When solutions of Fe(II) are oxidised in controlled laboratory conditions or during treatment of ash from municipal waste incinerators, two-line ferrihydrite forms rapidly and on never-dried samples, AFM shows abundant individual particles with diameter ranging from 0.5 to several tens of nanometers. Aging in solution at 70°C promotes growth of the particles into hematite and goethite and their identification (by X-ray powder diffraction, XRPD, with Rietveld refinement) becomes possible at the same aging stage as mineral morphology becomes recognisable by AFM. In other experiments that were designed to mimic natural attack by organic acids, colloidal lepidocrocite (γ-FeOOH) was observed in situ by AFM, while reductive dissolution removed material on specific crystal faces. Lath ends are eroded fastest while basal planes are more stable.

In order to help elucidate mechanisms of contaminant immobilisation by Fe-oxides, we examined samples from a reactive barrier made with 90% quartz sand, 5% bentonite and 5% zero-valent iron filings that had reacted with a solution typical of leachate from coal-burning fly ash using time-of-flight secondary ion mass spectroscopy (TOF-SIMS). Fe(0) oxidised to Fe(III), while soluble and toxic Cr(VI) was reduced to insoluble Cr(III). Chemical maps show Fe-oxide coatings on bentonite; Cr is associated with Fe-oxides to some extent but its association with Ca in a previously undescribed phase is much stronger. Other samples taken from municipal waste incinerator ash that had been treated by aeration in Fe(II) solutions were examined with transmission electron microscopy (TEM), selected area electron diffraction (SAED) and energy dispersive X-ray spectroscopy (EDS). Pb and some Zn are seen to be dispersed throughout two-line ferrihydrite aggregates, whereas Sn and some Zn are incorporated simply as a result of entrainment of individual ZnSn-oxide crystallites.

Geochemical speciation models that fail to account for contaminant uptake in solid solutions within major phases or as thin coatings or entrained crystals of uncommon phases such as those described here risk to underestimate contaminant retardation or immobilisation.     

Seismic risk assessment of buildings in urban areas: a case study for Denizli,Turkey

This study aims to carry out a seismic risk assessment for a typical mid-size city based on building inventory from a field study. Contributions were made to existing loss estimation methods for buildings. In particular, a procedure was introduced to estimate the seismic quality of buildings using a scoring scheme for the effective parameters in seismic behavior. Denizli, a typical mid-size city in Turkey, was used as a case study. The building inventory was conducted by trained observers in a selected region of Denizli that had the potential to be damaged from expected future earthquakes according to geological and geotechnical studies. Parameters that are known to have some effect on the seismic performance of the buildings during past earthquakes were collected during the inventory studies. The inventory includes data of about 3,466 buildings on 4,226 parcels. The evaluation of inventory data provided information about the distribution of building stock according to structural system, construction year, and vertical and plan irregularities. The inventory data and the proposed procedure were used to assess the building damage, and to determine casualty and shelter needs during the M6.3 and 7.0 scenario earthquakes, representing the most probable and maximum earthquakes in Denizli, respectively. The damage assessment and loss studies showed that significant casualties and economic losses can be expected in future earthquakes. Seismic risk assessment of reinforced concrete buildings also revealed the priorities among building groups. The vulnerability in decreasing order is: (1) buildings with 6 or more stories, (2) pre-1975 constructed buildings, and (3) buildings with 3–5 stories. The future studies for evaluating and reducing seismic risk for buildings should follow this priority order. All data of inventory, damage, and loss estimates were assembled in a Geographical Information System (GIS) database.     

二元混合模型在广东塘蓬花岗岩成因和物质来源研究中的应用

广东塘蓬岩体中发育两类不同成因及时代的花岗岩,采用Ｒｂ－Ｓｒ等时线年龄测定,确定塘蓬主岩体的成岩时代为海西期（２０４Ｍａ）;应用Ｓｍ－Ｎｄ同位素,采用二元混合模型确定塘蓬岩体黑云母二长花岗岩为改造型花岗岩,黑云母花岗闪长岩为同熔型花岗岩。     

Dynamic process analysis for the initiation and movement of the Donghekou landslide-debris flow triggered by the Wenchuan earthquake

The Donghekou landslide-debris flow was a remarkable geological disaster triggered by the Wenchuan earthquake in 2008. The dynamic process of a rapid landslide-debris flow is very complicated and can be divided into two aspects: the slope dynamic response of the earthquake and the mass movement and accumulation process. A numerical method combined with a finite difference method (FDM) and discrete element method (DEM) for simulation of landslide-debris flow under seismic loading is presented. The FDM and DEM are coupled through the critical sliding surface, initiation time and velocity. The dynamic response of the slope is simulated by the finite difference method, and critical sliding surface is determined using the earthquake response spectrum method. The landslide initiation time and the velocity are determined by time–history analysis. The mass movement and accumulation process is simulated using the discrete element method. Simulation results demonstrate that the maximum amplification coefficient of dynamic acceleration for the Donghekou slope is approximately 3.909, the initiation time of landslide is approximately 6.0 s, and the average initial velocity of the sliding mass is approximately 0.85 m/s. The failure of the slope is the result of elevation-orientated amplification effect and the sliding mass triggered with a small initial velocity. The numerical simulated result of the maximum sliding velocity is approximately 66.35 m/s, and the mass is disintegrated rapidly because of collision and free fall. The landslide velocity decreases when the flowing mass reaches a lower slope angle and gradually comes to a stop, and the total travel distance is approximately 2400 m.     

Mass movements in tropical volcanic terrains: the case of Teziutlán (México)

During the last decade, soil degradation coupled with global climate changes has increased hydrogeological hazards in Mexico. In tropical volcanic terrains, alteration processes have enhanced the formation of clay minerals that promote water retention and result in soil/rock weakness. Intense seasonal rainfall can trigger the liquefaction and remobilization of these low-resistance terrains. During the first week of October 1999, heavy rains affected eastern Mexico, including Puebla State. As a consequence, approximately 3000 mass movements, consisting of rock and soil slides and slips, debris flows and avalanches were generated in this area. In the town of Teziutlán (Puebla), which is located on volcanic deposits, a single mass-movement event caused approximately 150 deaths. In the present work we identified two types of mass movements in the Teziutlán area—Type 1: superficial erosion of an unwelded ignimbritic sequence forming small detrital fans, and Type 2: thin soil slide/debris flow from the remobilization of a volcanic sequence composed of clay-rich paleosols interbedded with ashfall horizons. The clay-rich volcanic paleosols favored the formation of perched water tables on a hydraulic aquiclude. Positive pore-water pressures triggered the failure. Based on these results, the principal human settlement in the Teziutlán area may be threatened by future debris flows, which could cause serious harm to the dense population and severe damage to its infrastructure. It is necessary to prevent future deaths and damage by installation of mitigative measures based on detailed studies. Without any further study, it will not be possible to prevent and mitigate a natural disaster with the same magnitude as the 1999 catastrophic hydrogeological phenomena.