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.     

Spectral Simulation and Conditioning to Local Continuity Trends in Geologic Modeling

Spectral simulation has gained application in building geologic models due to the advantage of better honoring the spatial continuity of petrophysical properties, such as reservoir porosity and shale volume. Distinct from sequential simulation methods, spectral simulation is a global algorithm in the sense that a global density spectrum is calculated once and the inverse Fourier transform is performed on the Fourier coefficient also only once to generate a simulation realization. The generated realizations honor the spatial continuity structure globally over the whole field instead of only within a search neighborhood, as with sequential simulation algorithms. However, the disadvantage of global spectral simulation is that it traditionally cannot account for the local information such as the local continuity trends, which are often observed in reservoirs and hence are important to be accounted for in geologic models. This disadvantage has limited wider application of spectral simulation in building geologic models. In this paper, we present ways of conditioning geologic models to the relevant local information. To account for the local continuity trends, we first scale different frequency components of the original model with local-amplitude spectrum ratios that are specific to the local trend. The sum of these scaled frequency components renders a new model that displays the desired local continuity trend. The implementation details of this new method are discussed and examples are provided to illustrate the algorithm.     

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.     

基于高光谱卫星遥感数据的UPDM分析方法

提出了一种用于处理多高光谱卫星数据的UPDM分析方法。研究结果证明．该方法应用干Landsat／TM(ETM^-)、Terra／MODIS和ADEOS—II／GLI等高光谱卫星传感器时．光谱重构均方根误差小于0．029适用于研究高光谱卫星遥感数据。     

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

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

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

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

基于DMN的高光谱图像分割方法研究

简述高光谱遥感光谱特征体系(包括光谱曲线特征、光谱变换特征和光谱度量特征3个层面)。研究马尔科夫网的概念和方法,生成基于光谱角(SA)特征度量的DMN,并以DMN为证据对高光谱图像进行分割;研究和实验表明基于SA信息的马尔科夫网可以很好地综合高光谱数据空间特征与光谱特征间的关系,为进一步数据处理提供优化控制(其实质是概率神经网络)。最后提出未来应用和研究方向。     

Can A Cup Anemometer `Underspeed'?

An analysis of cup-anemometer dynamics has been carried out inorder to determine whether the mean-wind velocity can have anegative bias. This would be contrary to the general belief thatcup anemometers always overspeed. Compared to prior analyses, theeffect of a possible nonlinearity of the calibration function isincluded. The conclusion is that neither longitudinal nor lateralvelocity fluctuations can contribute significantly to a negativebias. However, if a cup anemometer has an angular response thatfalls below the ideal cosine response, there will, as demonstratedin the concluding discussion, be a negative contribution from thevertical velocity fluctuations to the total bias, and thiscontribution may even outbalance the positive contributions fromthe longitudinal velocity fluctuations. Concrete evidence of suchexotic cup anemometer behaviour has not been reported in theliterature.