An Efficient Numerical Method for Computing Synthetic Seismograms for a Layered Half-space with Sources and Receivers at Close or Same Depths

— It is difficult to compute synthetic seismograms for a layered half-space with sources and receivers at close to or the same depths using the generalized R/T coefficient method (Kennett, 1983; Luco and Apsel, 1983; Yao and Harkrider, 1983; Chen, 1993), because the wavenumber integration converges very slowly. A semi-analytic method for accelerating the convergence, in which part of the integration is implemented analytically, was adopted by some authors (Apsel and Luco, 1983; Hisada, 1994, 1995). In this study, based on the principle of the Repeated Averaging Method (Dahlquist and Björck, 1974; Chang, 1988), we propose an alternative, efficient, numerical method, the peak-trough averaging method (PTAM), to overcome the difficulty mentioned above. Compared with the semi-analytic method, PTAM is not only much simpler mathematically and easier to implement in practice, but also more efficient. Using numerical examples, we illustrate the validity, accuracy and efficiency of the new method.     

A data‐based mechanistic modelling (DBM) approach to understanding dynamic sediment transmission through Wyresdale Park Reservoir,Lancashire, UK

This paper outlines the application of a new data‐based mechanistic (DBM) modelling methodology to the characterization of the sediment transmission dynamics in a small upland reservoir, Wyresdale Park, Lancashire. The DBM modelling strategy exploits advanced statistical procedures to infer the dynamic model structure and its associated parameters directly from the instrumented data, producing a parametrically efficient, continuous time, transfer function model which relates suspended sediment load at the reservoir inflow to the outflow at the event scale. The associated differential equation model parameters have physical attributes which can be interpreted in terms of sediment transmission processes and associated reservoir trap efficiency. Sedigraph analysis suggests that wind‐induced resuspension episodically supplies an additional load to the reservoir outlet. The stochastic nature of the DBM model makes it ideal for evaluating the effects of uncertainty through Monte Carlo simulations (MCS) for discharge and sediment transmission. Copyright © 2000 John Wiley & Sons, Ltd.     

The Influence of Topology on Hydraulic Conductivity in a Sand-and-Gravel Aquifer

A field experiment consisting of geophysical logging and tracer testing was conducted in a single well that penetrated a sand-and-gravel aquifer at the U.S. Geological Survey Toxic Substances Hydrology research site on Cape Cod, Massachusetts. Geophysical logs and flowmeter/pumping measurements were obtained to estimate vertical profiles of porosity ϕ, hydraulic conductivity K, temperature, and bulk electrical conductivity under background, freshwater conditions. Saline-tracer fluid was then injected into the well for 2 h and its radial migration into the surrounding deposits was monitored by recording an electromagnetic-induction log every 10 min. The field data are analyzed and interpreted primarily through the use of Archie's (1942) law to investigate the role of topological factors such as pore geometry and connectivity, and grain size and packing configuration in regulating fluid flow through these coarse-grained materials. The logs reveal no significant correlation between K and ϕ, and imply that groundwater models that link these two properties may not be useful at this site. Rather, it is the distribution and connectivity of the fluid phase as defined by formation factor F, cementation index m, and tortuosity α that primarily control the hydraulic conductivity. Results show that F correlates well with K, thereby indicating that induction logs provide qualitative information on the distribution of hydraulic conductivity. A comparison of α, which incorporates porosity data, with K produces only a slightly better correlation and further emphasizes the weak influence of the bulk value of ϕ on K.     

多维体视方法论及其在地球科学研究中的意义

本文阐述”多维体视方法论”,主张立体透视与动态、多维地看问题;它起源于医学CT(计算机层析学)透视成像、核磁共振(MRI)等体视学方法研究疾病的生理解剖结构、病理特征、监测人体有关器官的各种变化,是行之有效的方法;已成功地应用于检测工业器械、航天器及地下结构等物体的内部结构,从而发展和演绎出来的认识论方法。在医学领域将CT、MRI所得二维信息和断面图像、对病灶区引入“体显示”和“面显示”等可视化技术,将能充分反映不同疾病的病理、药理特征;在航天、地学等领域同样应用着CT多维“可视化技术”。本文试图发扬综合预测的学术思想,将医学临床应用CT、核磁共振、体视学方法结合研究其他制约因素指标,从而区分出疾病的原发、继发分期,导致绝症的进程和前兆,概括为“多维体视预测方法”;也适用于地震孕育、强震发生与火山喷发预测的研究。     

大数据环境下数字填图数据集成服务技术

应用数字填图技术形成了大量地质填图图幅数据,这些数据空间结构化和非结构化特征并存,如何在网络环境下提供高效的数据服务是急需解决的一个难题。大数据技术的发展为数字填图、数据集成服务提供了一种新的途径。通过对数字填图数据特征的分析,在研究地质调查信息网格大数据处理框架的基础上,提出了结构化和非结构化数据相结合的有序化组织管理、发布与服务方法,并对关键技术进行了研究与试验,取得了良好的效果。     

华北地块南缘花山、五丈山岩体LA-ICP-MS锆石U-Pb年龄、地球化学特征及成因

花山复式岩体位于华北地块南缘,采用LA-ICP-MS技术测得复式岩体中的中细粒似斑状黑云二长花岗岩的锆石U-Pb年龄为128±1Ma,此年龄被解释为成岩年龄。该岩体及邻近的五丈山岩体中花岗岩具有高硅、富碱特征,为准铝-弱过铝质高钾钙碱性系列;在稀土和微量元素上,富集LREE、Rb、Ba、K、Sr等大离子亲石元素,亏损HREE、Zr、Hf、Ta、Nb、P、Ti等高场强元素,总体上Eu异常不明显,显示出I-A过渡型花岗岩的特征。锆石Lu-Hf同位素组成表明,花山花岗岩的εHf(t)主要变化于-25.8~-19.6之间,tDM2为2.81~2.43Ga,显示源区物质以古老的壳源物质为主,有年轻组分的参与。与岩体外围多金属成矿年龄的对比分析发现,外围金、钼矿床与花山复式岩体的年龄一致,比五丈山岩体晚20~25Ma,表明金、钼矿的形成与花山岩体在时空上联系密切,该岩体可能为成矿提供了热源。     

Lake Winnipeg: geological setting and sediment seismostratigraphy

Lake Winnipeg, the seventh largest lake in North America, is located at the boundary between the Interior Plains and the Canadian Shield in Manitoba, Canada. Seismic profiles were obtained in Lake Winnipeg on two geoscientific cruises in 1994 and 1996. These data indicate the morphology of the bedrock surface. In most cases, a clear distinction between low relief Paleozoic carbonate rock and high relief Precambrian rock can be made. In northern Lake Winnipeg, the eastern limit of Paleozoic rock is clearly demarcated 30 km west of the previous estimate of its position. In southern Lake Winnipeg, all or most of the Paleozoic sequence terminates at a prominent buried escarpment in the centre of the lake. This indicates that Paleozoic rock on the eastern shore, known from drilling and outcrops, is an outlier. Major moraines are apparent as abrupt, large ridges having a chaotic internal reflection pattern. These include the Pearson Reef Moraine, the George Island Moraine and the offshore extension of The Pas Moraine. Little evidence for extensive or thick till was observed. Instead, fine-grained sediments deposited in glacial Lake Agassiz rest directly on bedrock over most of the lake basin. Hence an episode of erosion to bedrock was associated with glaciation and/or deglaciation. The Agassiz Sequence sediments are well-stratified, drape underlying relief and in some areas are over 100 m thick. In places, stratification in these sediments is disrupted, perhaps by dewatering. Evidence of erosion of Agassiz Sequence sediments by recent currents was observed. The contact between the Agassiz Sequence and the overlying Winnipeg Sequence sediments is a marked angular unconformity. The Agassiz Unconformity indicates up to 10 m of erosion in places. The low-relief character of this unconformity precludes subaerial erosion and the lack of till, moraines, or extensive deformation precludes glacial erosion. Waves appear to be the most likely erosional agent, either in waning Lake Agassiz or early Lake Winnipeg time. Winnipeg Sequence sediments, in places very thin, mantle most of the lakefloor. These sediments were deposited in the present Lake Winnipeg and are faintly stratified to massive and reach about 10 m in thickness in deep water. On the surface of the Winnipeg Sequence, vigorous, episodic currents are thought to contribute to the construction of flow-transverse sand waves as much as 6 m high in a deep, narrow constriction in the lake.     

Holocone paleovegetation and paleohydrology of a prairie pothold in southern Saskatchewan, Canada

The Andrews site represents one of countless prairie potholes found in areas of hummocky moraine on the northern Great Plains. Sediments from a depth of 5.8 to 3.1 m at this 'kettle-fill' site in the Missouri Coteau upland of southern Saskatchewan, Canada, provides a record of vegetation, climate, and hydrologic changes within a small, ca 30 m diameter, closed-drainage basin from ca 10.2 to 5.8 ka. Plant macrofossil analyses of 67 samples, 6 ¹⁴C ages, and stratigraphy were used to identify 5 zones, representing the paleohydrological changes that followed deglaciation in southern Saskatchewan.Results of this study indicate that with the melting of residual stagnant ice a pond (>2 m deep) with abundant aquatic, emergent, and shoreline plants developed in the basin at ca 10.2 ka and persisted until at least ca 8.8 ka. During this time there was a shift in upland vegetation from a white spruce forest (Zone II) to a deciduous parkland at ca 10 ka (Zone III). As climate warmed, brackish and alkaline conditions developed coincident with shallowing of the pond at the end of Zone III. The perennial water phase ended at ca 8.8 ka and was followed by a low-water stand lasting ca 1100 years. Prairie fires and slopewash from unstable slopes were dominant sedimentological processes until ca 7.7 ka (Zone IV). Water levels began to rise and between ca 7.7 and 5.8 ka a semi-permanent pond was established in a grassland setting (Zone V). After ca 5.8 ka this prairie pothole wetland became ephemeral, to the point that plant macrofossils could not be preserved, a situation continuing today. Interactions between climate change, variability in local groundwater supply, and sedimentological processes likely account for the paleohydrologic events reconstructed at the Andrews site.     

基于CFD技术的河谷型城市风环境模拟——以兰州市城关区为例

本文应用CFD软件对多种流场的强大模拟功能,模拟小尺度的兰州市城关区风环境。选用SRTM资料表征兰州市城关区的地形,应用ArcGIS生成兰州市地形模型,计算相关参数,将计算得到的地形粗糙度和兰州多年气候资料转化为CFD模型的数字参数及模拟边界条件等,首次模拟得到河谷型地区兰州城关区的风环境。结果表明:（1）河谷型地区中,气流过山时风速达最大值,依山体高度增大,而在背风坡风速值明显降低,达到0～1 m·s-1;（2）气流过境时在沿河谷前进,且在谷底风速值达0～1 m·s-1;（3）在模拟城市热岛效应的对比实验中,地块加热对气流过境有停滞作用。