汶川MemSsub8.0地震前地磁短临异常与强震预测探索l

强震前地磁垂直分量日变化异常主要表现在幅度与相位上的变化．利用中国地磁台网多年观测资料,研究了2008年5月12日四川汶川8.0级地震前,地磁垂直分量日变极小值出现的时间（低点时间）在空间分布上的异常特征,即低点位移现象;并用此方法预测了汶川8.0级地震6级余震平静2个月后的一组6级强余震．这种少见的异常现象与强震发生有较好的相关性．发震日期为异常出现后的第27天或第41天的前后4天,强震发生地区在低点位移突变分界线附近．此外,还研究了汶川8.0级地震临震前,震中周围出现的地磁日变幅异常现象．      

地磁低点位移异常界线交汇现象与地震关系研究

利用国家地磁台网中心产出的观测资料,研究中国大陆地区低点位移现象与强震之间的关系,认为地磁低点位移异常与强震之间具有较好的关联性。多次5级以上地震发生前两月内会出现一次或多次低点位移异常,地震往往发生在低点位移异常分界线300 km内。2022年1月8日青海门源6.9级地震、2016年1月21日青海门源6.4级地震和2019年9月16日甘肃肃南5.0级地震前两月内均出现了3次明显的地磁低点位移异常。将多次低点位移异常分界线叠加到同一幅图上,形成了1个或多个交汇区域,发现地震发生在低点位移分界线的交汇区域附近。综合活动构造断裂分布、历史大地震发生情况、中长期预测的强震潜在危险区及年度全国地震危险区进一步判断可能发震的区域,此项研究明显缩小了用地磁低点位移法预测地震发生地点的范围。     

汶川M_S8.0地震前地磁短临异常与强震预测探索

强震前地磁垂直分量日变化异常主要表现在幅度与相位上的变化.利用中国地磁台网多年观测资料,研究了2008年5月12日四川汶川8.0级地震前,地磁垂直分量日变极小值出现的时间(低点时间)在空间分布上的异常特征,即低点位移现象;并用此方法预测了汶川8.0级地震6级余震平静2个月后的一组6级强余震.这种少见的异常现象与强震发生有较好的相关性.发震日期为异常出现后的第27天或第41天的前后4天,强震发生地区在低点位移突变分界线附近.此外,还研究了汶川8.0级地震临震前,震中周围出现的地磁日变幅异常现象.     

地磁“低点位移”现象与昆仑山口西8.1级地震

利用中国地磁台网多年的观测资料,研究了地磁垂直分量日变极小值出现的时间(低点时间)在空间分布上的异常特征即低点位移现象．该异常现象与强震发生有较好的相关性,发震日期为异常日后的第27天或41天前后4天,地点在低点位移突变分界线附近．2001年11月14日昆仑山口西8．1级地震前的10月8日也出现了少见的地磁低点位移强震异常．讨论了地磁低点位移的“感应-触发”机制。     

地磁日变化异常与强震的关系

强震前地磁垂直分量日变化异常主要表现在幅度与相位上的变化 .本文利用中国地磁台网多年观测资料 ,研究了 2 0 0 1年 1 1月 1 4日昆仑山口西 8.1级地震、2 0 0 3年 2月 2 4日巴楚 伽师 6.8级地震和 1 989年 9月 2 2日小金 6.6级等 3 0多次 6.6级以上强震前 ,地磁垂直分量日变极小值出现的时间 (低点时间 )在空间分布上的异常特征 ,即低点位移现象 .这种少见的异常现象与强震发生有较好的相关性 .发震日期为异常日后的第 2 7天或 41天的前后 4天 ,强震发生地区在低点位移突变分界线附近 .此外 ,还研究了昆仑山口西 8.1级地震和巴楚 伽师 6.8级地震临震前 ,震中周围出现的地磁日变幅异常现象 .     

地磁日变化异常与强震的关系

强震前地磁垂直分量日变化异常主要表现在幅度与相位上的变化.本文利用中国地磁台网多年观测资料，研究了2001年11月14日昆仑山口西8.1级地震、2003年2月24日巴楚-伽师6.8级地震和1989年9月22日小金6.6级等30多次6.6级以上强震前，地磁垂直分量日变极小值出现的时间（低点时间）在空间分布上的异常特征，即低点位移现象. 这种少见的异常现象与强震发生有较好的相关性. 发震日期为异常日后的第27天或41天的前后4天，强震发生地区在低点位移突变分界线附近. 此外，还研究了昆仑山口西8.1级地震和巴楚-伽师6.8级地震临震前, 震中周围出现的地磁日变幅异常现象.      

北部湾6.1,6.2级地震前地磁,气象异常初析

探讨了北部湾地震地磁和气象的异常现象,初步结果表明,震前几个月至数天,地磁垂直分量的相关系数ｒ值,特征线斜率Ｋ值以及Ｚ日变形态双低点位移先后出现短临异常,气象因子分析结果表明,年雨量一阶差分大于４００ｍｍ可作为５级以上地震的年度预报参考指标。震前数月,－３．２ｍ的地温年变形态出现畸变。     

强震前的地磁异常现象

震前地磁Z分量的日变化异常主要表现在幅度和相位上的变化, 文中利用我国50多个地磁台多年的资料, 重点研究了1997年11月8日玛尼7.5级地震和2005年11月26日江西九江5.7级地震。 发现玛尼地震和九江地震以及大多数中强地震前, 全国地磁台站地磁Z分量日变化幅度的极小值出现时间(低点时间)在空间分布上异常特征, 即低点位移现象。 这种少见的异常现象与强震发生有较好的相关性。 发震日期为异常日后的第27天或第41天的前后4天, 强震发生地区在低点位移突变分界线附近。 结果显示, 在玛尼7.5级地震和江西九江5.7级临震前, 震中周围出现区域性Z分量日变化幅度异常现象。     

地磁低点位移现象及其与地震的关系

本文用中国地磁台网多年的资料,研究了地磁垂直分量日变极小值出现的时间(低点时间)在空间分布的异常特征,即低点位移现象。这种少见的地磁异常现象与地震的发生有较好的相关性,其发震日期为异常日以后的第27天或41天前后4天,发震地区在低点位移突变分界线附近。低点位移是一种复杂的地球物理现象,受多种因素影响,其机制仍处在探索阶段。      

新疆伽师MS6.4地震前地磁日变化异常分析

以新疆伽师M S6.4地震为案例,利用地磁日变化异常分析方法,对中国地磁台网所记录的地磁观测资料进行分析,结果显示:1、离震中最近的喀什台震前59 d和39 d出现一组逐日比高值异常,排除干扰,且地磁日变形态也出现明显的畸变,异常集中于震中附近台站,震中位于异常阈值线附近;2、震前,新疆及周边地区出现地磁垂直分量日变化幅度的极小值时间在空间上出现地磁低点位移现象,低点时间分界线横穿南天山西段,呈现近EW走向分布,分界线两侧低点时间相差超过2 h,发震日期为异常发生后的38(41±4)d,震中位于低点位移分界线附近;3、地磁逐日比和低点位移异常的发生时间和空间分布可能与伽师M S6.4地震的孕育和发生有非常高的关联性。该研究进一步推广和检验了地磁日变化异常分析方法在地震预报方面的应用。     

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.