Lead, zinc, and antimony contamination of the Rio Chilco-Rio Tupiza drainage system, Southern Bolivia

An intense, but localized rainfall event in February 2003, led to the severe erosion and failure of a tailings disposal impoundment at the Abarόa Antimony Mine in southern Bolivia. The failure released approximately 5,500 m³ of contaminated tailings into the Rio Chilco-Rio Tupiza drainage system. The impacts of the event on sediment quality are examined and compared to contamination resulting from historic mining operations in the headwaters of the basin. Of primary concern are contaminated floodplain soils located along downstream reaches of the Rio Tupiza which were found to contain lead (Pb), zinc (Zn), and antimony (Sb) concentrations that locally exceed Canadian, German, and Dutch guidelines for agricultural use. Spatial patterns in sediment-borne trace metal concentrations, combined with Pb isotopic data, indicate that Pb, Zn, and Sb are derived from three tributary basins draining the Abarόa, Chilcobija, and Tatasi-Portugalete mining districts. Downstream of each tributary, geographical patterns in trace metal concentrations reflect local geomorphic changes throughout the drainage system. Trace metal concentrations within the Rio Chilco decrease rapidly downstream as a result of dilution by uncontaminated sediments and storage of metal enriched particles (e.g., sulfide minerals) in the channel bed as a result of ongoing aggradation. Storage in the floodplains is limited. These processes significantly reduced the dispersal and, thus, the relative environmental affects of tailings eroded from the Abarόa Mine during the 2003 flood. In contrast, storage of Pb, Zn, and Sb in floodplains along the Rio Tupiza is significant, the majority of which is derived from historic mining operations, particularly mining within the Tatasi-Portugalete district.     

Fast computation of Hankel Transform using orthonormal exponential approximation of complex kernel function

The computation of electromagnetic (EM) fields, for 1-D layered earth model, requires evaluation of Hankel Transform (HT) of the EM kernel function. The digital filtering is the most widely used technique to evaluate HT integrals. However, it has some obvious shortcomings. We present an alternative scheme, based on an orthonormal exponential approximation of the kernel function, for evaluating HT integrals. This approximation of the kernel function was chosen because the analytical solution of HT of an exponential function is readily available in literature. This expansion reduces the integral to a simple algebraic sum. The implementation of such a scheme requires that the weights and the exponents of the exponential function be estimated. The exponents were estimated through a guided search algorithm while the weights were obtained using Marquardt matrix inversion method. The algorithm was tested on analytical HT pairs available in literature. The results are compared with those obtained using the digital filtering technique with Anderson filters. The field curves for four types (A-, K-, H-and Q-type) of 3-layer earth models are generated using the present scheme and compared with the corresponding curves obtained using the Anderson sc heme. It is concluded that the present scheme is more accurate than the Anderson scheme     

Centrifuge modelling-based verification of response spectra design methods for some vulnerable structures

Residential RC framed structures suffered heavily during the 2001 Bhuj earthquake in Gujarat, India. These types of structures also saw severe damage in other earthquakes such as the 1999 Kocaeli earthquake in Turkey and 921 Ji-Ji earthquake in Taiwan. In this paper the seismic response of residential structures was investigated using physical modelling. Idealised soft storey and top heavy, two degrees of freedom (2DOF) portal frame structures were developed and tested on saturated and dry sand models at 25 g using the Schofield Centre 10-m Beam Centrifuge. It was possible to recreate observed field behaviour using these models. As observed in many of the recent earthquakes, soft storey structures were found to be particularly vulnerable to seismic loads. Elastic response spectra methods are often used in the design of simple portal frame structures. The seismic risk of these structures can be significantly increased due to modifications such as removal of a column or addition of heavy water tanks on the roof. The experimental data from the dynamic centrifuge tests on such soft storey or top-heavy models was used to evaluate the predictions obtained from the response spectra. Response spectra were able to predict seismic response during small to moderate intensity earthquakes, but became inaccurate during strong earthquakes and when soil structure interaction effects became important. Re-evaluation of seismic risk of such modified structures is required and time domain analyses suggested by building codes such as IBC, UBC or NEHRP may be more appropriate.     

Early Estimation of Seismic Hazard for Strong Earthquakes in Taiwan

A shakemap system providing rapid estimates of strong ground shaking could be useful for emergency response providers in a damaging earthquake. A hybrid procedure, which combines site-dependent ground motion prediction models and the limited observations of the Real-Time Digital stream output system (RTD system operated by Central Weather Bureau, CWB), was set up to provide a high-resolution shakemap in a near-real-time manner after damaging earthquakes in Taiwan. One of the main factors that affect the result of ground motion prediction analysis is the existence of site effects. The purpose of this paper is to investigate the local site effects and their influence in the ground shaking and then establish an early estimation procedure of potential hazard for damaging earthquakes. Based on the attenuation law, the site effects of each TSMIP station are discussed in terms of a bias function that is site and intensity-level dependent function. The standard deviation of the site-dependent ground motion prediction model can be significantly reduced. The nonlinear behavior of ground soil is automatically taken into account in the intensity-level dependent bias function. Both the PGA and the spectral acceleration are studied in this study. Based on the RTD data, event correctors are calculated and applied to precisely estimate the shakemap of damaging earthquakes for emergency response.     

二维电子地图和三维虚拟场景的互响应

本文提出了基于二维电子地图与三维虚拟场景之间相互响应的思想来克服它们各自的局限性.简要的描述了相互响应的原理、表现形式、实现流程及实现方法.最后提出了一些结论和展望.     

基于预测的边缘检测方法

提出了一种基于预测的、有一定自适应性的边缘检测方法。介绍了其原理和算法实现过程。通过计算两个相邻像素的灰度值的均值、均方差和梯度来预测下一个像素的灰度值，比较预测值和真实值来判断下一个像素是否是边界点。通过与现有同类算法的比较，证明了该算法的优越性。     

地震属性分析中水平切片的应用

利用地震数据体中含有的丰富地质信息,可以有效地进行储层预测。层拉平水平切片地层学解释方法是先将地震剖面目的层段做层拉平处理,然后对切片顺地层层面追踪,勾绘出各种地层现象,并以地层学的观点加以解释。该方法在吉林探区扶新南部地区的应用中,对于识别河道的展布趋势取得了较好的效果,而且发现了多期叠置的扇型沉积体。另外,还可指导反演数据体的解释,做到对岩性体从定性到定量的描述。     

一种层状大地瞬变电磁响应正演计算的改进方法

频谱法是瞬变电磁法正演计算中常用的计算方法之一,其关键在于内层含贝塞尔函数的积分计算。对于重叠回线装置的瞬变电磁法正演的内层积分,以往采用在积分区间内寻找贝塞尔函数的零点分布,依次在相邻零点之间采用一般的数值积分法在求得各自的积分值后再叠加的方法进行计算。这种方法精度较高,但效率低。利用贝塞尔函数的大宗量渐近特性,提出一种新的计算方法。计算结果表明,该方法计算效率高、方法简单、精度可靠。     

现阶段地表油气化探技术难题与发展对策

虽然地表油气化探方法技术已取得显著的应用效果, 但是由于成功率计算标准不同及方法技术本身的缺陷, 导致人们对方法的认同度较低。当前, 地表油气化探技术的发展, 一方面要以可靠的实验依据、严谨的数理推导及理论模型应对来自方法技术本身的挑战, 另一方面也要利用其快速、廉价和直接的优势, 加强与地质、遥感、非震物探等方法综合应用。     

低频泥石流特征及防治 --以四川汶川县茶园沟为例

低频泥石流指30a以上发生1次的泥石流。尽管发生频率很低，但由于自身的特征及人们认识上的原因，其造l戒的人员伤亡和财产损失却居各类泥石流之首。对低频泥石流的研究在预防和减轻泥石流灾害方面有着重要意义。文章以2003年8月9日发生在四川汶川县茶园沟的泥石流为蒯，采用形态调查法和称重法等方法确定了流量、流速及重度。分析了其物质来源及形成的自然和人为原因，探讨了低频泥石流的特征。针对低频泥石流的低频性、破坏性和隐蔽性特点，同时考虑到低频泥石流沟的防治受保护对象的重要程度及资金投入等因素的限制，提出了进行泥石流危险性分区、搞好小流域水土深持、普及泥石流防灾减灾知识及加强监测和预警预报等措施，这对减轻低频泥石流区的灾害有一定的借鉴意义。