东大别山构造带六安—英山地磁剖面观测与地震研究

本文以东大别构造带内六安—英山地磁剖面磁异常数据为基础,结合研究区域地质学和地震学研究成果,反演给出研究区域居里面埋深及地下介质磁性结构模型.通过对磁异常变化特征曲线、介质磁性结构模型分析发现:东大别构造带北部磁异常较强且波动较大、地下磁性结构特征相对复杂,北淮阳地区磁异常较弱、地下介质结构成叠瓦状形态分布,合肥盆地磁异常曲线变化平缓、埋深较深,为连续均匀的磁性基底,结果较好地反应了合肥盆地、北淮阳、东大别构造带北部之间的差异;"霍山地震窗"地壳下方6～12km处存在显著地磁性差异,结合东大别构造带内70年代以来地震和历史强震,对"霍山地震窗"地震活动性与地下构造、场源关系进行了研究,初步探讨了其动力学成因.     

郯庐断裂带江苏段及邻区航磁异常与深部构造特征

收集了郯庐断裂带江苏段及周边地区较大比例尺的航磁异常数据,通过水平及垂向求导、欧拉反褶积处理、居里面反演及剖面建模,系统分析了区域地壳磁性结构、断层空间展布及深部构造特征,并揭示了其与地热、地震以及岩浆活动的相互关系.研究表明:(1)区域航磁异常场具有明显的分区性,线性特征明显的航磁异常高值区基本由断裂主控,为岩浆沿断层侵入、灌入导致,而线性特征不明显的航磁异常高值或者低值区主要与磁性基底的埋深有关,其中郯庐断裂带形成NNE向的串珠状高磁性异常条带.(2)航磁场及其方向导数图显示区域发育17条主要断裂,其中郯庐断裂带与周边的10条断裂处岩浆活动剧烈;欧拉反褶积结果显示研究区磁性基底埋深在郯庐断裂带处及其两侧不同构造单元处存在差异;居里面反演结果表明郯庐断裂带为居里面坳陷区,该处为热流值及地温梯度低值区.剖面建模结果揭示郯庐断裂带底部可能存在大范围的基性及超基性岩浆岩,且推测断层剖面延伸形态为张扭性的"负花状构造".(3)以无锡—宿迁断裂为界,郯庐断裂带南北两侧在原始航磁场、航磁梯度场以及欧拉反褶积结果中均体现出了明显的分段性,表现为北侧航磁异常值更高、航磁梯度场波动更加剧烈且磁性基底埋深更浅.     

利用福建地区磁测资料计算分析居里面深度

采用2014—2015年福建地区流动地磁观测资料分析福建区域岩石圈磁场变化情况,利用向上延拓处理方法对福建地区岩石圈磁场变化信息进行了局部与区域异常场分离,结果表明在漳州南部地区磁异常变化比较明显。基于区域异常场分离结果,反演居里等温面深度,初步得到福建地区居里面埋深大约为16~28 km。     

麻城1932年6.0级地震区地球物理场与深部环境研究

研究了麻城地区重、磁异常场特征、经反演计算获得了深部结构特征。对重力布格异常利用Ｐａｒｄｅｒ－Ｏｌｄｅｎｂｕｒｇ快速位场反演方法计算出莫霍面的展布，对航磁异常利用对率谱法和矩谱法计算出居里面深度。在此基础上，探讨该震区深部蕴震条件。主要结论是：①麻城－团凤断裂带为该区切割较深的主要控震构造；②１９３２年麻城６．０级地震发生在区域重力异常等值线梯度最大处，磁力异常区局部正异常间的低值带，莫霍面及居面     

基于融合重力异常的大别造山带东段莫霍面反演研究

利用实测布格异常和EGM2008重力异常融合结果,采用频率域位场反演方法计算大别造山带东段莫霍面三维空间分布,结合区域地质构造和地震活动等资料讨论大别造山带东段莫霍面分布特征及其构造含义等。研究结果显示:①莫霍面空间分布体现了块体构造差异,大别造山带莫霍面最深,最大深度达42km,显示东大别造山带存在明显山根,扬子地块深度次之,华北地块最浅;郯庐断裂带及其东侧区域存在莫霍面上隆带;②大别造山带北侧和南侧莫霍面陡变带分别位于青山—晓天断裂附近和襄樊—广济断裂以北,体现华北地块和扬子地块向大别造山带之下俯冲的构造特征,指示了深部构造缝合带位置;莫霍面深度特征表明郯庐断裂带区域构造边界带属性明显,其切割深度至少达到壳幔边界,大别造山带商城—麻城断裂两侧经历了差异隆升;③研究区域绝大多数地震发生在莫霍面以上深度,莫霍面深度陡变带、上隆带及不同莫霍面深度特征区的转换带为区域地震活动提供了深部构造条件。     

大陆深俯冲带的地壳速度结构——东大别造山带深地震宽角反射/折射研究

在安徽大别山(东大别)进行的深地震宽角反射/折射探测获得6条二维地壳速度结构剖面. 结果显示，东大别造山带地壳为一高速穹隆构造，在其核部中、下地壳变质岩出露于地表，波速高达5.0km/s；在其翼部，上、中地壳发育速度约6.1km/s的壳内低速层(体). 莫霍面的起伏变化较大，中心部位深达41km左右，周边地区则抬升到32～34km. 在晓天—磨子潭断裂一线下方莫霍面垂向错断，断距约4km. 东大别造山带具有大陆深俯冲-碰撞造山带地壳结构的典型式样. 莫霍面错断与扬子陆块深俯冲有关，错断处表征扬子与华北陆块碰撞缝合的深部位置. 高速穹隆构造可能是两陆块碰撞挤压的产物，穹隆翼部上、中地壳发育的低速滑脱带(面)可能在碰撞期之后的地壳伸展、超高压变质岩从中地壳抬升出露于地表过程中起到重要作用.     

地球物理综合观测揭示秦岭—桐柏—大别复合造山带地壳及上地幔结构

秦岭—桐柏—大别复合造山带(以下称为秦岭大别造山带)属于中国中央造山带的一部分,由华北克拉通与扬子克拉通汇聚形成.对于秦岭大别造山带及其周缘地区的研究,可以为这一大陆碰撞造山带的形成与演化过程提供重要信息.本文整合研究区域的接收函数与背景噪声数据,采用H-κ叠加分析、接收函数与背景噪声联合反演、克希霍夫偏移成像等方法,得到了沿秦岭东西方向具有高分辨率的地壳及上地幔结构.研究结果显示:(1)莫霍面深度由西向东逐步抬升,由剖面西侧最深约55 km上升至剖面东侧最浅约30 km;莫霍面于东西秦岭之间起伏明显;桐柏以及东大别下方莫霍面局部加深.(2)西秦岭中下地壳观测到的高速异常阻隔了青藏高原东北缘地壳低速异常的向东扩张,反映了青藏高原东北缘的中下地壳流没有通过西秦岭继续向东流动.(3)西秦岭岩石圈地幔顶部高速异常延伸至100 km深度(剖面底部),桐柏—西大别岩石圈地幔顶部高速延伸至70 km深度,东大别、东秦岭岩石圈地幔顶部未见较大深度范围的高速异常.     

基于重磁资料研究贵州省深部构造特征

贵州省在大地构造位置上属于羌塘—扬子—华南板块下的扬子陆块,受不同构造旋回期的发展演化控制,早期地质构造被掩盖,而研究深部构造有利于构造单元划分、成矿作用等研究.本文利用区域重磁资料,结合其他地质资料,利用优选向上延拓分离不同尺度重磁异常,并对不同尺度重磁异常进行归一化总水平导数垂向导数处理,提取线性构造信息,分析区内主要断裂构造特征,进行定性分析;根据分离出的大尺度重磁异常反演研究区莫霍面、居里面起伏特征,莫霍面起伏介于34～46 km之间,总体上具有西高东低的特征;居里面深度介于22～28.5 km之间,西部以坳陷区为主,东部以隆起区为主,以区内两条深大断裂为界,深部构造具有南北分带、东西分块的特征,反映出贵州省上扬子地块与江南造山带两个三级构造分区的深部构造特征差异.     

基于磁异常数据的居里面反演方法研究

本文总结了基于磁异常频谱特征反演居里面的基本理论和方法,分析了磁性体场源满足随机分布及分形分布两种假设下的磁异常频谱特征及其对数径向能谱与磁性体场源埋深之间的关系,并给出居里面埋深估计的方法.基于以上方法对华北地区磁异常数据进行分析并计算了相应的磁性体场源埋深,实际计算表明,基于磁性体场源随机分布假设的居里面深度计算结果与该地区地热结构有良好的对应关系.实际数据处理的结果验证了这些方法的可靠性和实用性.估算居里面深度的两种不同方法适用于不同的地质模型,实际应用应结合区域磁异常源分布的特点,以获取合理的居里面深度估计值.     

安徽霍山地震区深部电性结构和发震构造特征

霍山地震区位于大别造山带北缘华北板块与扬子板块接触带上,是大别造山带及周边地震活动最频繁、最集中的地区.83个大地电磁测点组成的大地电磁三维阵列覆盖了整个霍山地震区.用多重网格法、印模迭代重构法和非线性共轭梯度法对阵列数据进行三维带地形反演,获得了地震区深部三维电性结构.电性结构显示,北大别、北淮阳区的中上地壳为电阻率1000Ωm以上的高阻区,中下地壳为电阻率数十欧姆米的相对低阻区;六安盆地电阻率整体较低,中地壳存在显著的电阻率为几欧姆米的壳内高导层.北西向的晓天—磨子潭断裂分隔了北大别高阻层和北淮阳高阻层,在浅部向NE倾,深部向SW倾;北东向的落儿岭—土地岭断裂切穿北大别上地壳高阻层.小震双差定位结果表明,地震主要发生在NE向延伸的落儿岭—土地岭断裂附近的北大别、北淮阳中上地壳的高阻区,并集中于NW向的晓天—磨子潭断裂运动所造成的构造薄弱带中;2014年M S4.3霍山地震震源深度较深,位于北大别高阻区内部的电性梯度较大的区域.综合上述结果我们认为,霍山地震区的主要发震断裂为落儿岭—土地岭断裂,断裂的运动变形充分利用了晓天—磨子潭断裂早先活动所形成的构造薄弱带,断裂下方壳源高导体中的流体沿断层传播使断层强度弱化,使得这些薄弱带区易于发生小地震.由于北大别、北淮阳构造区显著高阻层的存在,我们认为霍山地震区存在发生6级以上中强震的深部孕震环境.     

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.