用时空转换系数判断震后趋势

本文定义地震的时空距ＤＳＴ格时空转换系数Ｃ为：Ｄ＾２ＳＴ＝ｄ＾２＋Ｃ＾２（△ｔ）＾２，ｄ为序列中连续发生的两次地震的空间距，△ｔ为时间间隔。采用６０年代来发生在我国５级以上有详细目录的地震序列２３个，分别计算了主震后１天，２－３天，４－７天内Ｃ的值Ｃ１，Ｃ３，Ｃ７，计算中取ＤＳＴ为相应时间段内地震空间距的最大值。结果表明：（１）若Ｃ１＞Ｃ３＞Ｃ７，则无强余震发生；（２）若Ｃ１≌Ｃ３，或Ｃ１＜Ｃ３＜     

大华北10次中强震前后广义时空分维的变化及其在中短期预报中的应用

本文应用高安秀树提出的广义分维方法分析计算了１９７０年以来大华北地区１０次Ｍｓ≥５．８级地震前后震中周围小震活动广义时（Ｄｔ）空（Ｄｓ）分维的变化，结果表明中强震前普遍检测到Ｄｔ，Ｄｓ的低值异常。具有实用意义的是至少有一条曲线可以统一地取Ｄｔ，Ｄｓ低于震前数年一倍均方差线，而另外一条曲线较同步地低于均值线达半年之久作为异常判别标准，中强震一般于异常结束后３－９个月发生，其中个别地震于异常之中发生。     

前震活动特征及其识别的研究（Ⅰ）

分析了近２０年来中国大陆强震（Ｍ≥７．０）和华北地区中强地震前广义前震序列时间和空间分布的多分形特征，计算其广义分维谱Ｄｑ和标度指数谱Ｆ（α），并探讨了它们随时间的变化．结果表明：（１）对广义前震活动，在所研究地区的时、空分布，均呈现多分形结构；（２）大地震发生前半年至一年，地震活动时、空分布的高次分维Ｄ∞呈下降趋势．这种降维特征，亦在实验室含水大样本岩石破坏前，声发射序列Ｄ∞的变化中观测到；（３）与中强地震相比，大地震孕育区地震活动的空间高次分维偏低，表明地壳应力分布的高度不均匀性，是板内大地震孕育环境中的本质因素之一．     

地震序列自仿射分形及自相似分形的计算与比较

本用标准偏差法和标度变换法，计算了云南省五个地震区带上１９６５－１９９５年间地震序列的自仿射分形的分维值Ｄ△ｔ（地震时间间隔序列自仿射分形）ＤＭ－ｔ（地震震级一时间序列自仿射分形）及自相似分形的分维值Ｄ０（地震时间分布自相似分形）得到了相应的时间变化曲线，ＤＭ－ｔ～ｔ，Ｄ△ｔ～ｔ和Ｄ０～ｔ。通过相互比较发现，其中ＤＭ－ｔ的随时间变化在中强地震前的反应最为敏感。     

前震活动特征及其识别的研究（Ⅰ）

分析了近２０年来中国大陆强震（Ｍ≥７．０）和华北地区中强地震前广义前震序列时间和空间分布的多分形特征，计算其广义分维谱Ｄ_ｑ和标度指数谱Ｆ（α），并探讨了它们随时间的变化．结果表明：（１）对广义前震活动，在所研究地区的时、空分布，均呈现多分形结构；（２）大地震发生前半年至一年，地震活动时、空分布的高次分维Ｄ_∞呈下降趋势．这种降维特征，亦在实验室含水大样本岩石破坏前，声发射序列Ｄ_∞的变化中观测到；（３）与中强地震相比，大地震孕育区地震活动的空间高次分维偏低，表明地壳应力分布的高度不均匀性，是板内大地震孕育环境中的本质因素之一．     

华北地区地震时间结构的多分形特征

利用固定半径法和直接计算法系统研究了华北地区（Ｎ３０°～４３°，λＥ１０６°～１２５°）１３个Ｋ表性子区地震时间结构的广义维数Ｄｑ—ｑ及标度指数谱ｆ（ａ）—ａ的特征，并与单—Ｄｑ—ｔ曲线进行了对比。结果表明；（１）华北地区地震时间结构为多分形结构，单一分形不能完整地描述其分形特征，应该用Ｄｑ—ｑ或ｆ（ａ）—ａ来加以分析；（２）华北强震前的０．５～２．０年．周围地区地震活动的时间多分形结构出现异常，主要表现为ｑ为负时Ｄｑ增大、ｑ为正时Ｄｑ减小，或ｆ（ａ）两端点下降，ａ跨度增大及顶点右偏等现象，这是对地震密集—平静现象的另一种描述；（３）华北地区地震时间多分形方法回顾性预报评分Ｒ＝０．４７，说明其有一定预报效能，是有前景的预报方法。     

1986年德都中强震群地震活动性再研究

本文运用“八五”以来地震短临预报攻关成果,对黑龙江省１９８６年德都中强地震前后地震活动性资料进行了再研究,结果表明,震前ｂ值无异常,η值,∑（ｔ）等地震活动性参数在一至一年半时段内有不同程度的异常显示,地震强度因子Ｍｆ值未出现高值异常,地震演化指数ＹＨ值和地震算法复杂性Ａｃ值的异常发展与“八五”短临地震预报判据相反,广义时间维Ｄ１和广义空间维Ｄｓ参数的异常发展过程符合“八五”地震短临预报判定指标。     

广东及邻区地震学参量b、η、Dt值与地震的关系

应用ｂ,η和广义时间分维Ｄｔ值的时间扫描计算方法,对广东省及邻区１９７０年以来发生的１５次Ｍｓ≥４．６级地震中周围的小震活动进行单项和综合分析研究,结果表明,震前２个月至１年半内,这三种参量大都存在明显的低值异常过程,与地震的对应率在６９％,表明该组参量可作为本区地震预测的有效判据,利用其综合预测效能,可提取本区地震的中短期异常标志。     

中国大陆的特征地震活动及其中—长期预测研究（一）

中国大陆不同地点２４个特征地震轮回的震级－时间序列分析表明，在一个完整轮回中，相继特征地震事件ｉ和ｊ之间时段的前１／２极少有破坏性地震发生；而在轮回的后１／２，大多数情况下至少有一次中－强地震发生，其中，中－强地震发生在轮回最后１／３时段的约占５４％。称这种地震为“填隙式中强震”（事件ｆ），其发生可能标志着当地自从上一次特征地震后，应变已重新积累 到达距下一次和下地震极限应变不远的背景，故对于发生     

大震前小震活动的时空分形特征

应用分形理论与方法，计算了唐山，海域，澜沧－耿马地震前小震活动能量分布的信息维及其在空间上的差异和时间上的变化，结果表明：（１）局域空间发维可以反映大震前小震活动能量分布的空间变化特征，大震发生在信息维中等值区；（２）局域时间分数维显示大震前小震活动存在的降维过程，且大震发生在降维之后的升维过程中或发生在降维过程结束之时。     

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.