地下水测年方法及其在地震监测中的应用展望

对地下水动力学特征的观测与研究,是地震前兆监测的重要方法之一。地下水动力学研究中,利用地下水测年方法,可以对地下水的补给、径流特征及更新能力作出定性和定量的分析。本文介绍了近年来国际上常用的地下水测年方法,总结了前人在该领域的主要研究进展,重点分析了地下水更新能力与地震前兆信息的关系、地下水运动规律与地震构造活动的关系以及地下水浅层补给与异常干扰排除的方法等。已有的研究成果表明,地下水年龄的测定与分析方法对于了解监测点映震能力、评价构造活动与地震活动程度以及在观测资料异常变化的现场核实等方面,可以发挥重要作用。     

华北地区地下流体中期和转折型短期异常成因分析

结合观测井孔的水文地质条件、大气降水资料 ,深入剖析了华北地区水位、水氡的中期和转折型短期异常变化的直接原因 ,认为水位异常是由于降水量的增大与减少造成的 ,其形成机制主要为降水渗入补给 ;水氡异常的形成原因主要为降水渗入 -混合机制。进而提出“深浅构造的互动作用” ,用以解释流体异常与地震的关系 ,其论点是 :降水、地下水变化可诱发 (活化 )浅层的构造变动 ,而浅层构造活动又引发深部构造活动和地震的发生 ;反之 ,深部构造的活动和地震的孕育 ,使浅层岩石的应力增强 ,在降水、地下水的诱发作用下 ,可导致浅层的构造活动。降水引起的地下水异常是一种广义的地震前兆     

水化学分析方法在地下水异常核实中的应用

地震前兆识别中的异常核实工作是提高地震分析预报能力的重要环节。"四图一法"是地下水异常核实中有效的水化学分析方法,其中Piper图和Schoeller图用于判定地下水化学类型和补给关系,Gibbs图用于判定地下水化学成分来源,Giggenbach三角图用于判定地下水的水岩平衡状态和是否有深部活动信息,以及混合比例法用于计算地下水的混合程度。同时,介绍了水化学分析方法在异常核实应用中的水样采集、测试方法和基本原理。使用水化学分析方法,有助于判定异常的构造或非构造活动影响因素,提高异常判定的科学水平。     

标准年动态曲线形态分析方法在地震地下流体氡数据处理中的应用

台站观测数据主要包含台站固有信息与前兆信息,对前兆信息的准确识别和提取,是提高运用氡观测站开展地震活动监测预测有效性的关键环节。采用标准年动态曲线形态分析方法,对甘肃省内8个地下流体台站的水(气)氡异常信息进行统计分析,给出8个地下流体台站水(气)氡的前兆异常映震能力与空间响应特征。研究表明,标准年动态曲线形态分析方法是识别地下水(气)氡前兆短临异常信息的有效方法,有望在地下流体前兆数据分析及异常落实工作中发挥重要的技术支撑作用。     

滇西实验场及邻区前兆群体异常与成组中强震的关系

本文使用前兆群体异常信息提取的方法,即应用前兆群体非均匀度ID和前兆综合信息量S数学分析,以及相关计算软件,对滇西地震预报实验场及邻区(N24～28°,E98～102°)的断层形变和地下流体各种监测手段的多台项观测资料进行了系统计算分析。研究结果表明:多台多测项前兆群体异常的映震效果优于多台单测项前兆群体异常的应震效果;前兆群体异常与中强地震的对应关系较好,且一组前兆群体异常一般会对应2次以上中强地震,显示了前兆群体异常与地震成组对应关系,而不是以往单次前兆与单个地震一一对应关系。     

金沙江水网及其地震地下水前兆监测能力的讨论

金沙江下游梯级水电站水库地震地下水动态监测网(简称金沙江水网),是中国第2个专门用于水库区地震前兆监测与研究为目的的企业地下水动态监测网.文中简要介绍了该网的布设、观测网的技术构成、观测井与观测含水层的基本特征及试运行的基本结果,从观测井-含水层特征及观测到的水位与水温多月、月、日动态特征,评估了其地震前兆监测能力.     

地壳形变与地下水动态异常特征研究进展

地壳变形与地下流体之间既相互联系又相互影响。地壳形变异常是地壳构造活动最直接的表现,是物理含义最为明确的地震前兆之一;地下流体异常能够提供多种型式的构造活动过程(包括地震)的前兆信息。两类异常以不同的方式和特点反映出地震孕育过程中伴随的前兆信息,从而成为地震预报的重要手段。本文综述了有关地壳形变与地下水动态异常特征方面的研究成果,从机理和前兆特征方面归纳了两种异常的相互关系。     

华北地区地震短期综合预测方法研究

在对华北地区114项骨干前兆观测资料进行地震短期异常前兆信息量提取与计算基础上，结合华北地区地质构造、地震活动、地壳厚度，特别是现今动力学环境因子等因素，研究了华北地区地震短期异常前兆信息量的时、空演化特征及其与地震的关系，给出了华北地区Ms≥5地震短期综合预测的指标和方法。经R值评分检验表明，该方法具有一定的地震预测效果。     

大同6.1级地震前后河北省地下水位前兆异常分析

识别地下水动态地震前兆异常,必须根据观测井地质及水文地质条件,分析地下水的多年动态特征,影响水位的干扰因素及异常变化与地震活动的关系。本文通过上述分析,认为大同6.1级地震前,河北省地下水网中万全、雄县、完县三口井存在水位异常变化。     

景东井水温异常变化与地震关系初探

系统分析景东观测井1998年以来的水温观测数据,研究水温前兆异常与地震之间的对应关系。结果表明,该井水源补给可能来自深部,水温变化不受季节性降雨渗入影响;震前异常全部为短期异常,以高温异常为主,短临异常突出,地震对应率50%,低温异常地震对应率80%;水温异常对发生在红河断裂和澜沧江断裂之间的地震指示意义最好,使用资料时考虑区域构造特点能很好地提高单个观测点的地震监视能力。     

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.