南华北地区地震活动特征分析

南华北盆地与秦岭大别造山带在构造上有耦合关系,地震活动水平相当。将南华北盆地和秦岭大别造山带组成的南华北地区作为研究对象,考察了其地震活动特征和青藏高原强震对其影响作用。发现南华北地区地震活动水平显著低于北华北地区,略低于下扬子地区。并且,通过对比分析青藏块体7.0级以上地震和南华北地区4.0级以上地震的对应关系,发现两个地区之间存在明显的成组对应现象,并表现2年左右的时间差,表明青藏块体的地震活动对南华北地区有明显的构造应力场传递作用。该结果为南华北地区中长期地震预测提供依据。2017年四川阿坝州九寨沟县发生7.0级地震,根据时间延迟的存在,对南华北地区未来中等地震的发生具有一定的指示意义。     

大华北地区应力场在汶川地震时的变化研究

利用大华北地区在汶川8级地震时有水位阶变的31口并孔的观测资料,通过计算各井水位系统的观测格值,进而估算了震时引起大华北地区应力场的调整变化.研究发现汶川地震后大华北中部地区的压应力得到增强,大华北东部和南部地区则表现为应力降低,对未来大华北地区的地震危险性预测应具有一定参考价值.     

华北地震区强震与月相的关系

计算了1587年以来的华北地区200个5级以上地震发生日的月相,其结果表明二者的相关性不明显,但分布在华北境内的中国东部重力梯级带上的101个较强地震,却与朔望日有着较密切的关系。     

大同—阳高2次地震前后震源应力场的变化过程

利用1989、1991年大同－阳高两次地震前后几年的地震记录资料，详细分析了其震源应力场的变化，结果表明，6.1级地震前震源应力场中两个主应力轴的方向曾发生逆转，主压应力轴逆转时间约10个月，主张应力轴逆转时间为2个月。6.1级主震前4个月，两个主应力轴的方向、仰角都恢复到正常状态，且与华北区域应力场的特征基本一致，而1991年5.8级地震前却没有出现上述现象，说明大同－阳高2次地震是在同一震源应力场作用下孕育、发展的，亦说明5.8级地震是6.1级地城的晚期强余震，指出山西地震带6级以上地震的孕育发生似是在华北区域应力场的作用之下完成的。     

从中等地震活动分析华北应力水平变化过程

分析汶川地震以来,华北地区M≥4.0中等地震的时间序列、空间分布特征,以及华北应力水平在汶川M8.0地震、东日本MW9.0地震前后的变化过程。结果显示,汶川地震对华北的影响基本止于华北西部的汾渭地震带,影响最大的时间段为2010年;东日本MW9.0地震对河北平原地震带及以东地区的地震活动影响大,特别是距离板块边界越近,影响越大,影响最大的时间段为2013年;华北地区的中等以上地震活动与周边较大尺度的构造运动具有一定关系;华北周边发生大地震,但并没有引起强震发生,可能与华北处于第四地震活动期的平静幕、应力水平低有重要关系。     

1998年张北地震序列构造应力张量变化特征的研究

根据张北地震序列的地震活动特征,利用张北震区1998年1月～1999年3月间145个3 级左右地震的震源机制解资料,应用Gephart(1990)的应力张量反演方法,研究张北地震序列构造应力张量的总体变化特征和时序变化特征;研究发现强余震发生前,应力张量的波动幅度不是很大,震源及附近地区的构造应力作用较强,强余震发生后,应力因发生破裂而重新分布。而且震源区的构造应力场方向与华北地区的构造应力场方向基本一致,说明张北主震和强余震都是在华北统一构造应力场的作用下发生的。     

华北地区5级地震平静现象研究

分析了华北地区1484年以来M≥5级地震平静现象：（1）华北地区1484-1730年及1815年以来2个地震活动期所有M≥5级地震平均平静时间为17．8个月,1731-1814年地震平静期的所有5级地震的平均平静时间约90．1个月,二者有明显差异;（2）华北地区2个地震活动期中相邻活动幕的幕间M≥5地震平均平静时间约86．4个月,与地震平静期大体一致,这样,90个月可以为华北地区M≥5级地震长期平静的标志性时间;（3）当前华北地区M≥5级地震平静时间达77个月,已接近长期平静的标志性时间,发生5级地震以上地震的可能已经很大;（4）上述工作成文后发生了2013年1月23日辽宁省辽阳市与沈阳市交界处5．1级地震,应注意华北进入新的地震活动幕。     

大同─阳高2次地震前后震源应力场的变化过程

利用１９８９、１９９１年大同－阳高两次地震前后几年的地震记录资料，详细分析了其震源应力场的变化，结果表明，６．１级地震前震源应力场中两个主应力轴的方向曾发生逆转，主压应力袖逆转时间约１０个月，主张应力轴逆转时间为２个月。６．１级主震前４个月，两个主应力轴的方向、仰角都恢复到正常状态，且与华北区域应力场的特征某本一致，而１９９１年５．８级地震前却没有出现上述现象。说明大同－阳高２次地震是在同一震源应力场作用之下孕育、发展的，亦说明５．８级地震是６．１级地震的晚期强余震。指出山西地震带６级以上地震的孕育发生似是在华北区域应力场的作用之下完成的。     

华北地区库仑应力演化与强震活动关系

采用更加符合长期变形和震后短期变形的Burgers流变模型,模拟了华北地区1303年以来6.5级以上地震引起的同震和震后库仑应力演化.结合华北地区历史地震期幕活动特征,分析了不同活跃期之间、同一活跃期内不同地震之间的库仑应力加卸载效应.结果表明:华北地区不同活跃期的强震主体活动区受控于历史强震的库仑应力加载作用;每个活跃期内强震活动主体区在空间上的迁移与该活跃期内首个7.5级以上地震的触发作用有关;华北地区每个活跃期内强震活动在时间上表现出的"平静期—活跃警告期—高潮期—剩余释放期"应该是区域动力加载过程的一种表现.岩石圈流变松弛效应引起的库仑应力变化对华北地区强震活动时空演化有显著的促进作用.本研究讨论了库仑应力变化在华北地区历史地震活动时空演化过程中可能扮演的角色,为探索华北地区强震活动空间迁移和韵律特征蕴含的构造动力学过程,以及与之相关的地震危险性判定提供参考依据.     

华北地区库仑应力演化与强震活动关系

采用更加符合长期变形和震后短期变形的Burgers流变模型,模拟了华北地区1303年以来6.5级以上地震引起的同震和震后库仑应力演化.结合华北地区历史地震期幕活动特征,分析了不同活跃期之间、同一活跃期内不同地震之间的库仑应力加卸载效应.结果表明：华北地区不同活跃期的强震主体活动区受控于历史强震的库仑应力加载作用;每个活跃期内强震活动主体区在空间上的迁移与该活跃期内首个7.5级以上地震的触发作用有关;华北地区每个活跃期内强震活动在时间上表现出的"平静期—活跃警告期—高潮期—剩余释放期"应该是区域动力加载过程的一种表现.岩石圈流变松弛效应引起的库仑应力变化对华北地区强震活动时空演化有显著的促进作用.本研究讨论了库仑应力变化在华北地区历史地震活动时空演化过程中可能扮演的角色,为探索华北地区强震活动空间迁移和韵律特征蕴含的构造动力学过程,以及与之相关的地震危险性判定提供参考依据.     

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.