玉树M_S7.1地震前电离层异常扰动分析

本文基于美国JPL实验室提供的全球电离层GPS TEC(GIM)数据和COSMIC掩星探测电子密度剖面,利用滑动四分位异常检测方法,综合分析了2010年4月13日青海玉树MS7.1地震前后15天内电离层异常的时空变化情况。研究结果表明,震前11天电离层TEC开始出现扰动异常,4月13日出现了明显的TEC增强现象,且峰值电子密度也存在相似的变化趋势,综合考虑空间活动水平等因素,认为这些异常很可能与玉树地震有关。     

2016年缅甸M7.2地震前中国区域电离层TEC扰动的时空分布

基于GPS TEC数据, 以2016年4月13日发生在缅甸的M7.2地震为例, 通过限定上、 下阈值的扰动提取方法分析了地震电离层效应在中国区域的时空分布。 结果显示, 该地震前2~4天, 中国区域的TEC扰动为地震电离层前兆; 其中4月11日(震前2天)电离层扰动受地震的影响最大, 扰动范围可达2816 km; 同时还发现4月11日异常最明显的区域在震中偏北方向(中国南部), 单独提取该区域经度链上的JPL TEC数据, 发现异常区域内北驼峰北移的现象, 并用静电场理论对异常区域向北偏移的物理机制做了分析解释。     

2020年8月18日印度尼西亚MS7.0地震电离层扰动研究

基于GPS TEC单站数据、 GIM TEC数据、垂测数据以及中国电磁监测试验卫星观测的原位电子密度数据,利用时序和空间的扰动提取方法,对2020年8月18日(世界时)印度尼西亚M_S7.0地震所引起的电离层扰动开展了研究。通过分析,我们发现多数据源的异常主要集中出现在8月12日—8月15日。震中附近8月12日GPS TEC单站数据和GIM TEC数据有明显的正异常现象,8月13日卫星原位电子密度和GIM TEC数据提取到正异常扰动,8月14日观测到卫星原位电子密度和GPS单站TEC数值有所增加,8月15日卫星原位电子密度、 GIM TEC和GPS单站TEC均出现正异常现象。此次地震引起的电离层异常不仅出现在震中上空,同时也反映在震中磁共轭区域。对应于8月12日震中附近的TEC异常,磁共轭区的单站TEC出现同步扰动,8月13日卫星数据、 GIM TEC和单站TEC数据在磁共轭区同步响应,8月15日磁共轭区的卫星原位电子密度、 GIM TEC、 GPS单站TEC和垂测数据均提取到了正异常扰动,并且与震中附近扰动出现的时间和幅度都有一定的共性。     

四川及周边地区TEC异常与印度尼西亚地震关联性研究

使用中国地壳运动观测网络和美国喷气动力实验室的GPS TEC数据,通过滑动四分位法,对2019年8月2日印度尼西亚苏门答腊岛南部海域6.8级地震进行分析,研究印度尼西亚地震与四川及周边地区TEC异常的关联性.基于对四川及周边地区以及全球TEC的监测,我们发现在2019年7月25日(震前一周),印度尼西亚地震震源区及其北半球的磁共轭区域曾发生不同程度的TEC异常扰动现象.结合2011-2019年四川及周边地区GPS TEC数据与印尼6.0级以上地震目录统计分析发现,大部分印尼6.0级以上地震会在四川及周边地区出现TEC异常扰动现象.在时间上,异常集中在震前10 d至震后4 d的范围内,主要发生在地方时6:00-11:00,在7:00(LT)和10:00(LT)两个时刻形成两个峰值,而在夜间出现相关的异常最少;在空间上,在四川地区存在TEC异常的印尼强震集中在印尼西侧,印尼西侧与四川地区TEC异常表现出良好的关联性,且TEC异常扰动现象面积约为32×104 km2,异常面积巨大,易于监测.研究结果表明,四川及周边地区TEC异常与印尼6.0级以上强震具有很好的时空关联性,结合全球TEC异常分布可以对印尼地区地震监测发挥重要作用.     

日本M_W 9.0级特大地震前电离层扰动初步分析

利用JPL提供的TEC地图数据和日本NICT提供的测高仪数据分析了2011年3月11日日本MW9.0级特大地震(Tohoku-oki Earthquake)前上空的电离层变化.GIM TEC地图数据分析发现,3月1～4日GPS TEC在赤道异常区域出现了不同程度的负扰动,其中3月1日和2日的扰动幅度和范围较大,以上很可能是由于3月1日爆发的9小时磁暴所致;此外,震前3天即3月8日04-14UT期间赤道异常区电离层TEC出现了明显的增强现象,异常最大幅度达30TECU,并伴随有南半球磁共轭区增强现象.日本境内的4个电离层测高仪数据分析显示,3月8日4个台站同时出现了foF2增强现象,其中Okinawa和Yamagawa两个站增幅明显.分析发现震前3天地震上空区域电离层参量异常增强,认为这可能与大气层和电离层垂直电场的变化有关,这种电场的变化有可能是地震产生的电磁信号传播所致.     

基于地基GPS-TEC的中国区域地震电离层效应研究

基于中国陆态网络地基GPS-TEC观测,针对2008～2019年发生在中国区域的7个M_S≥6.0地震,采用滑动四分位距法分析了地震前后的电离层扰动时空分布特征。结果显示,5个地震的震前2～6天,GPS-TEC值出现负异常扰动,地震发生期间及震后电离层TEC出现正异常扰动,主要集中在震后2～7天;GPS观测站距离震中越近,垂直上空的TEC扰动越明显,扰动空间最大范围可达2000km。随着震级的增加,震前电离层TEC异常扰动的发生率有所增加,且异常覆盖的范围也有所扩大。因此,认为震前一周内的电离层TEC变化可能提供揭示电离层扰动与地震活动之间关系的线索。     

2011年3月24日缅甸地震期间电离层TEC异常分析

基于IGS提供的2011年3月全球TEC数据,利用滑动时窗进行异常识别的方法,分析电离层赤道异常区TEC在缅甸地震期间的异常情况。结果表明：在震前10 d内,地震上空电离层TEC出现了明显的异常扰动,电离层异常北驼峰在16日至18日期间出现了明显的较大区域负异常现象,且在18日时北驼峰明显向磁赤道漂移,负异常在19日时基本消失;而在地震发生的当天和震后一天出现了小范围的正异常,随后正异常消失。在地震期间地磁活动水平总体都相当平静,且太阳处于一个相对较稳定的中低活动水平,电离层出现的异常很可能与此次缅甸地震有关。     

广州地区低纬电离层TEC扰动与两次强地震的关系研究

利用广州地区GPS-TEC监测网获得的TEC数据,分析了两次7.0级以上地震与低纬电离层TEC扰动存在的可能联系。对TEC的时序分析表明:两次地震发生前TEC的日变幅度均出现下降。但两次地震异常之间存在一定差异。台湾南部邻海7.2级地震前10 d内出现了负异常和正异常,但汶川8.0级地震前8 d和后6d内只出现负异常。此外,还讨论了低层电离层TEC扰动与地磁活动及热层动力学之间可能存在的联系。     

基于GPS和COSMIC数据分析汶川地震TEC和NmF2扰动

使用相关方法和welch法计算功率谱密度研究汶川地震前后GPS TEC在时域和频域上的变化,并进一步分析N_mF2空间分布变化.得出主要结论, 震前5月9日开始TEC平常具有良好的相关性被打破,24 h全日波、12 h半日波成分功率谱密度幅度值变化明显.5月12日至6月11日N_mF2变化幅度明显增大.研究结果显示,TEC和N_mF2时空分布变化存在电离层震前扰动和同震响应, 并且TEC震后的异常变化相比震前更加显著.表明使用的方法可以有效检测到地震发生前后电离层异常扰动信号.     

2020年新疆于田6.4级地震前电离层扰动现象分析

基于中国区域GPS反演TEC、 JPL TEC mapping、 张衡一号卫星探测等离子体参量数据, 着重分析了2020年6月26日新疆于田6.4级地震前的电离层异常现象, 结合之前于田发生的两次7级以上地震, 研究认为于田地震前异常集中出现在震前一周内, 以上升异常为主, 异常有明显的局地效应, 部分异常在磁共轭区有同步效应。 多参量综合分析增强了异常的判识能力, 并提高了异常的可靠性。 太阳及空间磁扰活动对电离层地震判识有较大影响, 会极大增加全球异常频次, 但日食现象引起的TEC扰动与地震电离层异常有明显差异, 较易区分, 弱磁扰活动下的地震电离层异常判识能力有待加强。     

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.