差分GPS在中国地磁监测网地理方位角测量中的应用

中国地磁监测网由地磁台站和流动地磁测点构成,实现对观测网内地磁场的相对记录和绝对观测。它的绝对观测采用世界上最先进的F、D、I测量仪器组合,其中地磁偏角D是地理北和地磁北的夹角。对地磁偏角D的绝对测量要求在地磁台站投入正常观测前必须完成观测墩和观测标志的地理方位角测量,流动地磁测点在每次地磁偏角D观测前后至少进行地理方位角测量各一次。据地磁监测网《观测规范》要求,地磁偏角D的观测精度必须≤6″,所以中国地磁监测网采用天文和差分GPS地理方位角两种测量方法,且两种方法观测精度相当,均≤6″。与天文测量相比,差分GPS观测不受天气等自然因素的限制,操作简单,观测快捷,数据处理计算机程序化。随着GPS技术的普及,中国地磁监测网差分GPS地理方位角测量方法已基本取代了天文方法。本文主要介绍差分GPS在中国地磁监测网地理方位角测量中的应用。     

GPS测量方位角在琼中地磁台建设中的应用

介绍了GPS的工作原理和用GPS测量方位角的原理与方法,使用GPS测量方位角方法测量了琼中台6个仪器墩的方位角,取得较好的结果,并对测量结果精度和误差来源进行了分析.测量结果证明GPS测量方位角方法不受观测时间限制,受天气影响较小,对人员要求不高,观测时间只需2～3天,具有较大的应用前景.地磁台在新建、改造观测室后都必须进行方位角的测定工作.经过若干年后,为了检验观测墩和方位标的稳定性和可靠性,需要定期或不定期地进行方位角的复测工作.     

铁磁钻具引起测斜仪方位角误差的校正方法研究

基于煤矿井下近水平钻孔姿态测量的实际需求,阐述了采用三轴加速度计与三轴磁强计构成磁测斜仪测量钻孔倾斜角、工具面向角和方位角的具体方法,分析了影响测量精度的各种误差来源.针对铁磁钻具引入的附加干扰磁场会造成方位角测量误差,提出了借助当地地磁场对方位角误差进行校正的方法.测试结果表明,相比于使用前、后无磁钻具时方位角的测量结果,去掉测斜仪前后的无磁钻具并采用该校正方法,仍可获得高精度的方位角.该校正方法只需在软件算法上做改进就可以使测斜仪钻铤抛掉配套的前、后无磁钻具,在保证测量精度的同时大幅节约了仪器成本并拓宽了测斜仪的应用.     

GPS测量方位角方法在琼中地磁台建设中的应用

介绍了GPS测量方位角的原理与方法。使用GPS测量了琼中地震台4个观测墩的方位角,取得较好的结果,并对测量结果的精度和误差来源进行了分析。结果证明,GPS测量方位角方法不受观测时间限制,受天气影响较小,对人员要求不高,观测时间只需2～3天,具有较大的应用前景。     

肇庆地磁台天文方位角测量及精度估计

介绍了肇庆地磁台天文方位角的测量方法、方案及测量结果。同时，做了测量精度的估计，测量精度完全符合《地震及前兆数字观测技术规范》的要求。     

北京地磁台天文方位角测量及误差分析

论述了北京地磁台天文方位标志的选点定位、造标埋设、方位角的测量方法。对８个观测墩位一个测点进行了方位角测量。通过计算方位角绝对精度ΔＡ＝２．４９″，平均观测均方差Ｍ＝±１．１４″，相当于三等天文观测精度，能满足地磁偏角６″的精度要求。     

浅谈新疆强震台站传感器极性检测和方位角校正

通过新疆强震台站的加速度传感器极性检测和方位角校正工作,发现新疆强震台站加速度传感器方位角最大误差为83.7°。按照国家强震动台网中心强震台站加速度传感器安装规范,加速度传感器方位角误差应该严格控制在2°以内,经过校正后,新疆区内全部强震台站的加速度传感器方位角误差都控制在1°以内,在很大程度上提高了新疆强震动台网的数据质量和可靠性。     

大连地磁台天文方位角测量

论述大连地磁台天文方位角的测量原理、方法及观测结果,在地磁观测室未封顶的条件下,能在观测墩直接观测北极星,避免比测墩转角换算等因素增加的观测误差,观测精度提高.     

基于半解析法有效和快速估计GRACE全球重力场的精度

首先基于半解析法建立了新的GRACE卫星K波段测量系统星间测速、GPS接收机轨道位置和加速度计非保守力误差联合影响累计大地水准面的误差模型;其次,基于各关键载荷精度指标的匹配关系,论证了误差模型的可靠性;最后,基于美国喷气动力实验室（JPL）公布的2006年的GRACE Level 1B实测误差数据,有效和快速地估计了120阶全球重力场的精度,在120阶处累计大地水准面的精度为18.368 cm,其结果和德国地学研究中心（GFZ）公布的EIGEN-GRACE02S全球重力场模型符合较好. 本文的研究为将来国际卫星重力测量计划（如GRACE Follow-On, 360阶）中高阶全球重力场模型精度的有效和快速估计提供了理论基础和计算保证.     

乌鲁木齐地磁台标志方位角的测定

介绍了乌鲁木齐地磁台测量标志方位角的方案、测量方法、原理。利用已有的观测成果,采用几何转换的方法精密测量了新建成观测室内的6个观洲墩相对于2个方位标志的方位角,成功的将旧观测室原标准墩方位角数据迁移到新观测室内各观测墩上,观测精度符合技术规范的要求,使新观测设施能够及时投入到观测工作中。并对地磁观测环境改造完成后的全面测量工作作出规划,在知识、技术上都做好了准备。     

A Monte Carlo study of rainfall forecasting with a stochastic model

A procedure for short-term rainfall forecasting in real-time is developed and a study of the role of sampling on forecast ability is conducted. Ground level rainfall fields are forecasted using a stochastic space-time rainfall model in state-space form. Updating of the rainfall field in real-time is accomplished using a distributed parameter Kalman filter to optimally combine measurement information and forecast model estimates. The influence of sampling density on forecast accuracy is evaluated using a series of a simulated rainfall events generated with the same stochastic rainfall model. Sampling was conducted at five different network spatial densities. The results quantify the influence of sampling network density on real-time rainfall field forecasting. Statistical analyses of the rainfall field residuals illustrate improvement in one hour lead time forecasts at higher measurement densities.     

Earthquake Engineering and Engineering Vibration Aims And Scope

正This journal is established by the Institute of Engineering Mechanics(IEM),China Earthquake Administration,to promote scientific exchange between Chinese and foreign scientists and engineers so as to improve the theory and practice of earthquake hazards mitigation,preparedness,and recovery.To accomplish this purpose,the journal aims to attract a balanced number of papers between Chinese and     

Solute concentrations in the lower niger river and the source rock contribution

Water quality analyses for the Niger River for the 1980/81 hydrological year are presented. The samples were collected from the main river at Lokoja, and from two main tributaries, the Kaduna and the Benue Rivers. Different water types were distinguished by the concentrations of major ions. The type Ca > Na > Mg > K - HCO₃ > SO₄ > Cl was represented at all stations during at least part of the year. Chloride was found to dominate the sulphate ion in the Kaduna and Niger, while the Benue maintained a higher concentration of sulphate relative to chloride all year round. Distinct patterns of seasonal variation in the ion concentrations were observed, particularly for the samples collected at Lokoja. Low ion concentrations were prominent during periods of high discharge, while low flow periods coincided with high dissolved ion concentrations. The contribution of rainwater to the total dissolved solids in the river waters was assessed indirectly using rainwater chemistry data from the Gulf of Guinea. The estimated rainwater contribution to the Lower Niger amounts to 5.15 mg 1^?1. Geochemical weathering calculations involving reactions of the four major minerals of granitic rocks - anorthite, biotite, albite, and K-feldspar - with carbon dioxide and water, can account for the average water composition of the Lower Niger. The proportion of the ionic components was also related to the occurrence of the respective element in the minerals.     

Scenario predictions for potential near-field and far-field earthquakes affecting Hong Kong

Historical earthquakes noted in the written records of the South China region, including Hong Kong, are not well delineated along identified prominent fault sources. Despite the lack of any definitive, localised trend in the spatial distribution of seismic activity in the region, there does appear to be some major disparity in the seismic activity rates (especially for large magnitude earthquakes) between the near-field and the far-field regions of Hong Kong. Despite this observation, previous studies of the regional seismic activity and seismic ground motion hazard (the latter using a probabilistic seismic hazard assessment, PSHA) have considered very broad source zone regions, in which uniform levels of seismic activity have been assumed. The present paper further scrutinises this broad source zone (BSZ) approach by adopting a novel expanding circular disc (ECD) method to determine the rates of earthquake recurrence. Such a method is intended to counter-check previously developed models by determining earthquake scenario events in terms of magnitude–distance (M–R) pairs or combinations, having defined values of average return period. Unlike the BSZ approach, the ECD method specifically accounts for the supposed variations in the seismic activity rates between events in the near-field and the far-field of Hong Kong.The form of the developed method is particularly suited to the determination of design-level earthquake ground motions for bedrock sites, since it assumes a directionally-independent attenuation model as described in the companion paper. It is found that, whilst the BSZ approach may indicate the overall average levels of hazard that are representative of the South China region as a whole, it does not capture the large disparity in seismic activity rates between near-field and far-field events. This important feature is expected to have a significant impact upon engineering assessments of the seismic safety of structures in Hong Kong and elsewhere in the South China region. For example, it is found that for events with M≥6, the seismic activity rate (normalised by time and area) in the very far-field is around 3.5 times larger than in the near-field and medium-field of Hong Kong. The resulting design M-R combinations, covering a range of return periods from 70 to 2500 years, are limited, for very long return periods and for distant events, by the maximum credible earthquake (MCE) magnitude. Intensive research to determine this seismic hazard parameter is recommended, in order to refine further the results of the ECD analysis, which presently conservatively assumes the MCE to range between M=6 in the near-field of Hong Kong to M=8 in the very far-field, at distances greater than 280 km from Hong Kong.     

Depositional ages and characteristics of Middle–Upper Jurassic and Lower Cretaceous lacustrine deposits in southeastern Mongolia

Lower Cretaceous lacustrine oil shales are widely distributed in southeastern Mongolia. Due to the high organic carbon content of oil shale, many geochemical studies and petroleum exploration have been conducted. Although most of the oil shales are considered to be Early Cretaceous in age, a recent study reveals that some were deposited in the Middle Jurassic. The present study aims at establishing depositional ages and characteristics of the Jurassic and Cretaceous lacustrine deposits in Mongolia. The Lower Cretaceous Shinekhudag Formation is about 250 m thick and composed of alternating beds of shale and dolomite. The Middle Jurassic Eedemt Formation is about 150 m thick and composed of alternating beds of shale, dolomitic marl, and siltstone. The alternations of shale and dolomite in both formations were formed by lake level changes, reflecting precipitation changes. Shales were deposited in the center of a deep lake during highstand, while dolomites were formed by primary precipitation during lowstand. Based on the radiometric age dating, the Shinekhudag Formation was deposited between 123.8 ±2.0 Ma and 118.5 ±0.9 Ma of the early Aptian. The Eedemt Formation was deposited at around 165–158 Ma of Callovian–Oxfordian. The calculated sedimentation rate of the Shinekhudag Formation is between 4.7 ±2.6 cm/ky and 10.0 ±7.6 cm/ky. Shales in the Shinekhudag Formation show micrometer‐scale lamination, consisting of algal organic matter and detrital clay mineral couplets. Given the average thickness of micro‐laminae and calculated sedimentation rate, the micro‐lamination is most likely of varve origin. Both Middle–Upper Jurassic and Lower Cretaceous lacustrine oil shales were deposited in intracontinental basins in the paleo‐Asian continent. Tectonic processes and basin evolution basically controlled the deposition of these oil shales. In addition, enhanced precipitation under humid climate during the early Aptian and the Callovian–Oxfordian was another key factor inducing the widespread oil shale deposition in Mongolia.     

THEORETICAL APPROACH TO THE SOLUTION OF THE INFILTRATION PROBLEM

ABSTRACT

The one-dimensional transient downward entry of water in unsaturated soils is investigated theoretically. The mathematical equation describing the infiltration process is derived by combining Darcy's dynamic equation of motion with the continuity and thermodynamic state equations adjusted for the unsaturated flow conditions. The resulting equation together with the corresponding initial and boundary conditions constitues a mathematical initial boundary value problem requiring the solution of a nonlinear partial differential equation of the parabolic type. The volumetric water content is taken as the dependent variable and the time and the position along the vertical direction are taken as the independent variables. The governing equation is of such nature that a solution exists for t > 0 and is uniquely determined if two relationships are defined, together with the specified state of the system, at the initial time t = 0 and at the two boundaries. The two required relations are those of pressure versus permeability and pressure versus volumetric water content.

Since the partial differential equation has strong non-linear terms, a discrete solution is obtained by approximating the derivatives with finite-differences at discrete mesh points in the solution domain and integrated for the corresponding initial and boundary conditions. The use of an implicit difference scheme is employed in order to generate a system of simultaneous non-linear equations that has to be solved for each time increment. For n mesh points the two boundary conditions provide two equations and the repetition of the recurrence formula provides n—2 equations, the total being n equations for each time increment. The solution of the system is obtained by matrix inversion and particularly with a back-substitution technique. The FORTRAN statements used for obtaining the solution with an electronic digital computer (IBM 704) are presented together with the input data.

Analysis of the errors involved in the numerical solution is made and the stability and convergence of the solution of the approximate difference equation to that of the differential equation is investigated. The method applied is that of making a Fourier series expansion of a whole line of errors and then following the progress of the general term of the series expansion and also the behavior of each constituent harmonic. The errors (forming a continuous function of points in an abstract Banach space) are represented by vectors with the Fourier coefficients constituting a second Banach space. The amplification factor of the difference equation is shown to be always less than unity which guarantees the stability of the employed implicit recurrence scheme.

Experiments conducted on a vertical column packed uniformly with very fine sand, show a satisfactory agreement between the theoretically and experimentally obtained values. Many experimental results are shown in an attempt to explain the infiltration phenomenon with emphasis on the shape and movement of the wet front, and the effects of the degree of compaction, initial water content and deaired water on the infiltration rate.     

Late Permian post‐ophiolitic trondhjemites from Central Iran: a mark of subduction role in growth of Paleozoic continental crust

Late Permian trondhjemites in the Anarak area occur as stocks and dykes, which cross cut the Anarak ophiolite and its overlying metasedimentary rocks, and are exposed along the northern Anarak east–west main faults. These leucocratic intrusive bodies have enclaves of all ophiolitic units and metamorphic rocks. They are composed of amphibole, plagioclase (oligoclase), quartz, zircon and muscovite. Secondary minerals are chlorite (pycnochlorite), epidote, albite, magnetite and calcite. Whole‐rock major‐ and trace‐element analyses reveal that they are characterized by high SiO ₂ (67.8–71.0 wt%), Al ₂ O ₃ (14.9–17.1 wt%) and Na ₂ O (5.3–8.6 wt%), low K ₂ O (0.1–1.5 wt%; average: 0.8 wt%), low Rb/Sr ratio (0.01–0.40; average: 0.09), low Y (3–6 ppm), negative Ti, Nb and Ta anomalies, slightly negative or positive Eu anomaly, LREE enrichment and fractionated HREE. These rocks present 2 to 40 times enrichment in inclined chondrite‐normalized REE patterns. Geochemical characteristics of the Anarak trondhjemites all reflect melting of a mafic protolith at more than 10 kbar. The field evidence and whole‐rock chemistry reveal that these rocks have been crystallized from magmas derived from melting of subducted Anarak oceanic crust. This study reveals that melting of garnet amphibolite was an important element of continent formation in the study area.