Mineralogical and chemical evolution of the Ernest Henry Fe oxide–Cu–Au ore system, Cloncurry district, northwest Queensland, Australia

The Ernest Henry Cu–Au deposit was formed within a zoned, post-peak metamorphic hydrothermal system that overprinted metamorphosed dacite, andesite and diorite (ca 1740–1660 Ma). The Ernest Henry hydrothermal system was formed by two cycles of sodic and potassic alteration where biotite–magnetite alteration produced in the first cycle formed ca 1514±24 Ma, whereas paragenetically later Na–Ca veining formed ca 1529 +11/−8 Ma. These new U–Pb_titanite age dates support textural evidence for incursion of hydrothermal fluids after the metamorphic peak, and overlap with earlier estimates for the timing of Cu–Au mineralization (ca 1540–1500 Ma). A distal to proximal potassic alteration zone correlates with a large (up to 1.5 km) K–Fe–Mn–Ba enriched alteration zone that overprints earlier sodic alteration. Mass balance analysis indicates that K–Fe–Mn–Ba alteration—largely produced during pre-ore biotite- and magnetite-rich alteration—is associated with K–Rb–Cl–Ba–Fe–Mn and As enrichment and Na, Ca and Sr depletion. The aforementioned chemical exchange almost precisely counterbalances the mass changes associated with regional Na–Ca alteration. This initial transition from sodic to potassic alteration may have been formed during the evolution of a single fluid that evolved via alkali exchange during progressive fluid-rock interaction. Cu–Au ore, dominated by co-precipitated magnetite, minor specular hematite, and chalcopyrite as breccia matrix, forms a pipe-like body at the core of a proximal alteration zone dominated by K-feldspar alteration. Both the core and K-feldspar alteration overprint Na–Ca alteration and biotite–magnetite (K–Fe) alteration. Ore was associated with the concentration of a diverse range of elements (e.g. Cu, Au, Fe, Mo, U, Sb, W, Sn, Bi, Ag, F, REE, K, S, As, Co, Ba and Ca). Mineralization also involved the deposition of significant barite, K(–Ba)–feldspar, calcite, fluorite and complexly zoned pyrite. The complexly zoned pyrite and variable K–(Ba)–feldspar versus barite associations are interpreted to indicate fluctuating sulphur and/or barium supply. Together with the alteration zonation geochemistry and overprinting criteria, these data are interpreted to indicate that Cu–Au mineralization occurred as a result of fluid mixing during dilation and brecciation, in the location of the most intense initial potassic alteration. A link between early alteration (Na–Ca and K–Fe) and the later K-feldspathization and the Cu–Au ore is possible. However, the ore-related enrichments in particular elements (especially Ba, Mn, As, Mo, Ag, U, Sb and Bi) are so extreme compared with earlier alteration that another fluid, possibly magmatic in origin, contributed the diverse element suite geochemically independently of the earlier stages. Structural focussing of successive stages produced the distinctive alteration zoning, providing a basis both for exploration for similar deposits, and for an understanding of ore genesis.     

超级计算机上矩阵乘的并行计算与实现

数值预报系统中经常要用到矩阵乘运算。在分布式超级计算机 (如IBM-SP) 上, 矩阵乘的并行计算需要较多的数据移动, 有效的数据传输对矩阵乘的实现至关重要。该文讨论了两种矩阵乘的并行算法, 一种是基于矩阵的列-行划分方式, 一种是基于矩阵的网格划分方式。在IBM-SP计算机上的实验结果表明, 网格划分的矩阵乘并行算法通讯开销更小, 并行效率更高, 其并行加速比较列-行并行算法改善约10 %。     

Near Wall Flow over Urban-like Roughness

In this study, comprehensive measurements over a number of urban-type surfaces with the same area density of 25% have been performed in a wind tunnel. The experiments were conducted at a free stream velocity of 10 m s^-1 and the main instrumentation was 120 ° x-wire anemometry, but measurement accuracy was checked using laser Doppler anemometry.The results haveconfirmed the strong three-dimensionalityof the turbulent flow inthe roughness sublayer and the depths of the inertial sublayer (log-law region) and roughness sublayer for each surface have been determined. Spatial averaging has been used to remove the variability of the flow in the roughness sublayer due to individual obstacles and it is shown that the spatially averaged mean velocity in the inertial sublayer and roughness sublayer can,together, be described by a single log-law with a mean zero-plane displacement and roughness length for the surface, provided that the proper surface stress is known. The spatially averaged shear stresses in the inertial sublayer and roughness sublayer are compared with the surface stress deduced from form drag measurements on the roughness elements themselves.The dispersive stress arising from the spatial inhomogeneity in the mean flow profiles was deduced from the data and is shown to be negligible compared with the usual Reynolds stresses in the roughness sublayer. Comparisons have been made between a homogeneous (regular element array) surface and one consisting of random height elements of the same total volume. Although the upper limits of the inertial sublayer for both surfaces were almost identical at equivalent fetch, the roughness sublayer was much thicker for the random surface than for the uniform surface, the friction velocity and the roughness length were significantly larger and the `roughness efficiency' was greater. It is argued that the inertial sublayer may not exist at all in some of the more extreme rough urban areas. These results will provide fundamental information for modelling urban air quality and forecasting urban wind climates.     

A comparison of tectonic ambient shear stress value in China with that in western USA

Introduction The stress status of the earths crust is closely related to global structure, the direction of plate movement, the drive mechanism of plate, earthquake cause, etc. Therefore the researches of the stress status of the earths crust and the process of dynamic evolution are paid a great attention to by many people around the world (Zoback, Zoback, 1980; Zoback, et al, 1989; Zoback, 1989; Spence, 1997; Fuchs, 1997; Plenefish, Bonjer, 1997; Muller, et al, 1997; XU, et al, 1989; WAN…     

中国与美国西部平均构造剪应力值的比较

提出了一种估计一个地区平均构造剪应力值的方法．用这种方法，选取哈佛大学公布的1977～1999年共15 993次地震的矩张量数据，估计了美国西部19个地区、中国和邻区43个地区(每个地区为1010范围)的平均环境剪应力值．结果表明，美国西部南加州西部海域和南加州应力最高，达13.7和12.0 MPa，然后向北、向南和向东逐渐递减，但最小也达8.7 MPa，是最高值的63%．中国新疆西北部地区和西藏察隅地区应力水平最高，达17.2和12.9 MPa，比美国高．中国的华北、云南、四川、台湾和美国南加州的应力水平差不多．中国南北地震带的应力水平为13 MPa左右，比南加州略高．两个重要地区的平均剪应力值分布图，提供了地学的基本数据．这些结果可为研究地震活动的大背景提供依据，对研究强地面运动参数(如峰值加速度及反应谱等)的衰减关系也是有用的．      

东营市森林火险天气等级预报系统

利用9210传输的气象资料，选用国家气象中心T106和HLAFS数值预报产品中的要素作为火险天气因子，根据森林火险原因，建立了东营市森林火险天气等级预报系统。     

帕隆藏布河流纵剖面演化的最小功模式

河流地貌最小功原理表明河流纵误时面形态将向力示使流速增大的方向演化，据此推导出矩形流域全程流速的平均值（u^-）与河流纵误时面形态指数（N）的关系式：u^∝-f(N)={1/2-1/[(N 1)(N 2)]}^1/2，f(N)被称为流速函数。这一最小功模式表明河流纵剖面演化方向是N由小变大。以西藏帕隆藏布中上游流域系统、干流全程及Ⅰ级阶地为例，由最小功模式计算的全程流速均值与实测值相吻合，从而检验了上述最小功原理和数学模式。     

Optimum multiple tuned mass dampers for structures under the ground acceleration based on the uniform distribution of system parameters

The five MTMD models, with natural frequencies being uniformly distributed around their mean frequency, have been recently presented by the first author. They are shown to have the near‐zero optimum average damping ratio (more precisely, for a given mass ratio there is an upper limit on the total number, beyond which the near‐zero optimum average damping ratio occurs). In this paper, the eight new MTMD models (i.e. the UM‐MTMD1～UM‐MTMD3, US‐MTMD1～US‐MTMD3, UD‐MTMD1 and UD‐MTMD2), with the system parameters (mass, stiffness and damping coefficient) being, respectively, uniformly distributed around their average values, have been, for the first time here, proposed to seek for the MTMD models without the near‐zero optimum average damping ratio. The structure is represented by the mode‐generalized system corresponding to the specific vibration mode that needs to be controlled. Through minimization of the minimum values of the maximum dynamic magnification factors (DMF) of the structure with the eight MTMD models (i.e. through the implementation of Min.Min.Max.DMF), the optimum parameters and values of Min.Min.Max.DMF for these eight MTMD models are investigated to evaluate and compare their control performance. The optimum parameters include the optimum mass spacing, stiffness spacing, damping coefficient spacing, frequency spacing, average damping ratio and tuning frequency ratio. The six MTMD models without the near‐zero optimum average damping ratio (i.e. the UM‐MTMD1～UM‐MTMD3, US‐MTMD1, US‐MTMD2 and UD‐MTMD2) are found through extensive numerical analyses. Likewise, the optimum UM‐MTMD3 offers the higher effectiveness and robustness and requires the smaller damping with respect to the rest of the MTMD models in reducing the responses of structures subjected to earthquakes. Additionally, it is interesting to note, by comparing the optimum UM‐MTMD3 with the optimum MTMD‐1 recently investigated by the first author, that the effectiveness and robustness for the optimum UM‐MTMD3 is almost identical to that for the optimum MTMD‐1 (without inclusion of the optimum MTMD‐1 with the near‐zero optimum average damping ratio). Recognizing these performance benefits, it is preferable to employ the optimum UM‐MTMD3 or the optimum MTMD‐1 without the near‐zero optimum average damping ratio, when installing the MTMD for the suppression of undesirable oscillations of structures under earthquakes. Copyright © 2003 John Wiley & Sons, Ltd.     

System identification of instrumented bridge systems

Several recorded motions for seven bridge systems in California during recent earthquakes were analysed using parametric and non‐parametric system identification (SI) methods. The bridges were selected considering the availability of an adequate array of accelerometers and accounting for different structural systems, materials, geometry and soil types. The results of the application of SI methods included identification of modal frequencies and damping ratios. Excellent fits of the recorded motion in the time domain were obtained using parametric methods. The multi‐input/single‐output SI method was a suitable approach considering the instrumentation layout for these bridges. Use of the constructed linear filters for prediction purposes was also demonstrated for three bridge systems. Reasonable prediction results were obtained considering the various limitations of the procedure. Finally, the study was concluded by identifying the change of the modal frequencies and damping of a particular bridge system in time using recursive filters. Copyright © 2003 John Wiley & Sons, Ltd.     

Average recurrence intervals of strong earthquakes and potential risky segments along the Taiyuan-Linfen portion of the Shanxi graben system

Introduction In AD 1303, the great Hongtong, Shanxi, earthquake of magnitude 8 caused a very serious disaster, which killed over one hundred thousands people at least (Department of Earthquake Dis- *aster Prevention, State Seismological Bureau, 1995). On the occasion of commemorating this ca-tastrophe having occurred for 700 years, we have important problems that need to be answered: How long the average recurrence interval of the grea…