青海南部熔积岩的发现:对寻找VMS型矿床的重要启示

熔积岩指的是侵入、混合到未固结或弱固结的湿沉积物中的熔浆分解、原位形成的一类特殊岩石。正确地认识该类岩石,有利于增进人们对岩浆-水(沉积物)相互作用过程的理解,恢复古环境。在青海南部沱沱河地区发现了一套角砾为撕片状、锯齿状及浑圆状的安山岩,胶结物为铁硅质组合的特殊熔积岩。研究表明,该熔积岩的角砾为岩浆遇水后快速淬火、裂解的产物,铁硅质组合为海底喷气沉积形成的含铁建造;且安山岩与含铁建造发生混合时,含铁建造尚未固结。该套熔积岩的发现,改变了长期以来对开心岭铁矿为火山热液交代安山岩而形成的认识,对于在矿区寻找VMS型矿床、区域内寻找海底热水喷流沉积型矿床具有重要的启示意义。     

The vibrational frequencies of forsterite Mg2SiO4: an all-electron ab initio study with the CRYSTAL code

The vibrational spectrum of Mg₂SiO₄ olivine was calculated at the Γ point by using the periodic ab initio CRYSTAL program. An all electron localized Gaussian-type basis set and the B3LYP Hamiltonian were employed. The full set of frequencies (35 IR active, 36 Raman active, 10 “silent” modes) was computed and compared to experimental data from different sources (four for IR and four for Raman). A generally good agreement is observed with experiment (the mean absolute difference ranging from 7 to 10 cm⁻¹ for the various sets), when some of the experimental frequencies, whose attribution is uncertain or appears to be affected by large errors, are not taken into account. A small number of observed peaks are not consistent with calculated frequencies, and a few theoretical peaks do not correspond to measured values. The implications are discussed in detail. The full set of modes are characterized using different tools, namely isotopic substitution, direct inspection of the eigenvectors and graphical representation, so as to obtain a consistent mode assignment.     

Paleoliquefaction study of the Clarendon–Linden fault system, western New York State

The lack of earthquake-induced liquefaction features in Late Wisconsin and Holocene sediments in Genesee, Wyoming, and Allegany Counties suggests that the Clarendon–Linden fault system (CLF) did not generate large, moment magnitude, M≥6 earthquakes during the past 12,000 years. Given that it was the likely source of the 1929 M 4.9 Attica earthquake, however, the Clarenden–Linden fault system probably is capable of producing future M5 events. During this study, we reviewed newspaper accounts of the 1929 Attica earthquake, searched for earthquake-induced liquefaction features in sand and gravel pits and along tens of kilometers of river cutbanks, evaluated numerous soft-sediment deformation structures, compiled geotechnical data and performed liquefaction potential analysis of saturated sandy sediments. We found that the 1929 M 4.9 Attica earthquake probably did not induce liquefaction in its epicentral area and may have been generated by the western branch of the Clarendon–Linden fault system. Most soft-sediment deformation structures found during reconnaissance did not resemble earthquake-induced liquefaction features, and even the few that did could be attributed to non-seismic processes. Our analysis suggests that the magnitude threshold for liquefaction is between M 5.2 and 6, that a large (M≥6) earthquake would liquefy sediments at many sites in the area, and that a moderate earthquake (M 5–5.9) would liquefy sediments at some sites but perhaps not at enough sites to have been found during reconnaissance. We conclude that the Clarendon–Linden fault system could have produced small and moderate earthquakes, but probably not large events, during the Late Wisconsin and Holocene.     

长江中下游成矿带陆内斑岩型矿床的成岩成矿作用

陆内环境斑岩型矿床的发现对斑岩成矿理论的完善具有重要意义。长江中下游成矿带作为中国东部重要的陆内成矿带之一,成矿带内发育多个重要的斑岩型矿床,如铜山口Cu-Mo矿床、鸡冠嘴Cu-Au矿床、白云山Cu矿、城门山Cu-Mo矿床、武山Cu-Mo矿床、丰山洞Cu-Au矿床、丁家山Cu矿、洋鸡山Au矿、沙溪Cu-Au矿床、冬瓜山Cu-Au矿床、舒家店Cu矿床和安基山Cu矿床等。本文选取成矿带内典型的、具有代表性的斑岩型矿床,对其地质特征(地层、构造、含矿斑岩、脉体特征和围岩蚀变)、成岩成矿年代、成矿岩体的岩石化学和成岩成矿地球化学等方面的研究资料和成果进行了系统总结,讨论和试图阐明长江中下游成矿带陆内斑岩型矿床的成岩成矿作用与成矿模式。研究显示,长江中下游成矿带形成于燕山期陆内造山过程,成矿斑岩岩浆活动和成矿作用主要发生于149~105Ma之间,进一步可以分为早、中、晚三阶段:149~135Ma、133~125Ma和123~105Ma,三阶段岩浆活动和成矿作用主要发生于成矿带中的断隆区,早阶段(149~135Ma)和晚阶段(123~105Ma)多为斑岩-矽卡岩型矿化,中阶段(133~125Ma)矿化为典型的斑岩型矿化。长江中下游成矿带内斑岩型矿床的含矿斑岩为高钾钙碱性-钙碱性系列岩石,大部分具有埃达克岩的地球化学特征,可能为源自富集地幔的岩浆和加厚下地壳部分熔融的岩浆混合的产物,源自富集地幔的基性岩浆对成矿具有至关重要的作用,它的混入使得混合岩浆富水、硫和金属(Cu、Au)等。进一步通过与岩浆弧环境的斑岩型矿床对比研究发现,长江中下游成矿带斑岩型矿床一般不发育高级泥化岩帽(advanced argillic liithocaps)以及浅部的高-中硫矿化蚀变系统,含矿岩浆源区性质和成矿物质来源等与岩浆弧环境的斑岩型矿床明显不同。     

The importance of iron speciation (Fe+2/Fe+3) in determining mineral assemblages: an example from the high‐grade aluminous metapelites of southeastern Madagascar

Metapelitic granulites from the Anosyen domain of southeastern Madagascar are exposed in three intercalated formations: the Amparihy, Bakika and Ihosy formations. Although mineralogically distinct from each other, the rocks from these formations show very similar bulk‐rock compositions when measured on a Fe_T basis. The preserved mineral assemblages thus do not reflect differences in the ratios of the main rock‐forming oxides (i.e. Al₂O₃:Fe_T:MgO), but instead reflect variations in the pre‐metamorphic oxidation state of the protolith rocks. These differences in oxidation state are manifested via differences in iron speciation – either Fe⁺² or Fe⁺³. The relatively reduced rocks of the Amparihy Formation preserve the assemblage bi–sp–sill–g–cd, which contrasts markedly with the mostly garnet and spinel‐absent bi–cd–sill–mt assemblages preserved in the strongly oxidized rocks of the Ihosy Formation. Compositionally intermediate rocks of the Bakika Formation are garnet bearing, but sillimanite‐absent, and contain the assemblage sp–g–cd–mag. Modelling of these rocks in the Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O–TiO₂–O system suggests that they evolved along a heating and cooling P–T path with only limited decompression accompanying cooling on the retrograde path. Peak temperatures and pressures of ～880–920 °C and 6–6.5 kbar are inferred for the majority of the Anosyen domain, with slightly lower peak temperatures (～840 °C) estimated in the extreme northwest of the area. The high‐temperature and relatively low‐pressure nature of metamorphism suggests high geothermal gradients existed during orogenesis, which in southern Madagascar is related to the amalgamation of Gondwana (580–520 Ma). Although metamorphic temperatures may have been augmented via thermal advection from the emplacement of the syn‐ to post‐tectonic Ambalavao suite, the high geothermal gradients nevertheless suggest thin and consequently hot lithosphere existed prior to orogenesis.     

Petrogenetic modelling of strongly residual metapelitic xenoliths within the southern Platreef,Bushveld Complex,South Africa

Xenoliths of quartz‐absent Fe‐rich aluminous metapelite are common within the platinum group element‐rich mafic/ultramafic magmatic rocks of the Platreef. Relative to well‐characterized protoliths, the xenoliths are strongly depleted in K₂O and H₂O, and have lost a substantial amount of melt (>50 vol.%). Mineral equilibria calculations in the NCKFMASHTO system yield results that are consistent with observations in natural samples. Lower‐grade rocks that lack staurolite constrain peak pressures to ～2.5 kbar in the southern Platreef. Smaller xenoliths and the margins of larger xenoliths comprise micro‐diatexite rich in coarse acicular corundum and spinel, which record evidence for the metastable persistence of lower‐grade hydrous phases and rapid melting consequent on a temperature overstep of several hundred degrees following their incorporation in the mafic/ultramafic magmas. In the cores of larger xenoliths, temperatures increased more slowly enabling progressive metamorphism by continuous prograde equilibration and the loss of H₂O by subsolidus dehydration; the H₂O migrated to xenolith margins where it may have promoted increased melting. According to variations in the original compositional layering, layers became aluminosilicate‐ and/or cordierite‐rich, commonly with spinel but only rarely with corundum. The differing mineralogical and microstructural evolution of the xenoliths depends on heating rates (governed by their size and, therefore, proximity to the Platreef magmas) and the pre‐intrusive metamorphic grade of the protoliths. The presence or absence of certain phases, particularly corundum, is strongly influenced by the degree of metastable retention of lower‐grade hydrates in otherwise identical protolith bulk compositions. The preservation of fine‐scale compositional layering that is inferred to be relict bedding in xenolith cores implies that melt loss by compaction was extremely efficient.     

PRINCIPLES AND APPLICATION OF MAXIMUM KURTOSIS PHASE ESTIMATION1

Methods of minimum entropy deconvolution (MED) try to take advantage of the non-Gaussian distribution of primary reflectivities in the design of deconvolution operators. Of these, Wiggins’(1978) original method performs as well as any in practice. However, we present examples to show that it does not provide a reliable means of deconvolving seismic data: its operators are not stable and, instead of whitening the data, they often band-pass filter it severely. The method could be more appropriately called maximum kurtosis deconvolution since the varimax norm it employs is really an estimate of kurtosis. Its poor performance is explained in terms of the relation between the kurtosis of a noisy band-limited seismic trace and the kurtosis of the underlying reflectivity sequence, and between the estimation errors in a maximum kurtosis operator and the data and design parameters. The scheme put forward by Fourmann in 1984, whereby the data are corrected by the phase rotation that maximizes their kurtosis, is a more practical method. This preserves the main attraction of MED, its potential for phase control, and leaves trace whitening and noise control to proven conventional methods. The correction can be determined without actually applying a whole series of phase shifts to the data. The application of the method is illustrated by means of practical and synthetic examples, and summarized by rules derived from theory. In particular, the signal-dominated bandwidth must exceed a threshold for the method to work at all and estimation of the phase correction requires a considerable amount of data. Kurtosis can estimate phase better than other norms that are misleadingly declared to be more efficient by theory based on full-band, noise-free data.     

Epithermal environments and styles of mineralization: Variations and their causes, and guidelines for exploration

Epithermal precious- and base-metal deposits are diverse, reflecting the different tectonic, igneous and structural settings in which they occur, the complexities of their local setting, and the many processes involved in their formation. Most epithermal deposits form at shallow crustal levels where abrupt changes in physical and chemical conditions result in metal deposition and attendant hydrothermal alteration. The principal factors that influence the conditions prevailing in the epithermal environment, and which ultimately determine the sites and character of mineralization, include: geology (structure, stratigraphy, intrusions and rock type, which affect the style and degree of permeability and the reactivity of the host); pressure and temperature (which in the epithermal environment are related on the boiling point with depth curve); hydrology (the relationship between permeability and topography which governs fluid flow, and discharge/recharge characteristics, as well as access of steam-heated waters); chemistry of the mineralizing fluid (which determines the metal-carrying capacity, as well as the associated vein and alteration assemblage); and syn-hydrothermal development of permeability and/or changes in hydraulic gradients.Many attempts have been made to classify epithermal deposits based on mineralogy and alteration, the host rocks, deposit form, genetic models, and standard deposits. All have their strengths and weaknesses. We prefer a simple approach using the fundamental fluid chemistry (high or low sulfidation, reflecting relatively oxidized or reduced conditions, respectively) as readily inferred from vein and alteration mineralogy and zoning, together with the form of the deposit, and using comparative examples to clarify the character of the deposit.Guidelines for exploration vary according to the scale at which work is conducted, and are commonly constrained by a variety of local conditions. On a regional scale the tectonic, igneous and structural settings can be used, together with assessment of the depth of erosion, to select areas for project area scale exploration. At project area scale, direct (i.e. geochemical) or indirect guidelines may be used. Indirect methods involve locating and interpreting hydrothermal alteration as a guide to ore, with the topographic and hydrologic reconstruction of the system being of high priority. These pursuits may involve mineralogic, structural, geophysical or remote sensing methods. On a prospect scale, both direct and indirect methods may be used; however, they can only be effective in the framework of a sound conceptual understanding of the processes that occur in the epithermal environment, and the signatures they leave.     

SEISMIC CHARACTER MAPPING OVER RESERVOIR INTERVALS1

This paper presents some results from an investigation into the utility of pattern recognition methods in seismic interpretation. The seismic instantaneous attributes of amplitude, phase and frequency provide a way of quantifying the character of a simple reflection. Measures of character can be developed from cross-plots and cluster analysis of these attributes. It is demonstrated that such seismic character can produce better-defined maps than a single attribute. These procedures can be extended to attributes derived from seismic trace segments, such as trace energy and centre frequency, and to multitrace attributes, but more effort is then needed to analyse the attributes and search out useful ones. An introduction is given to projection pursuit which has proved a useful exploratory tool for the anlysis of attribute relationships. It is important to stress that pattern recognition techniques simply help bring relationships and patterns in the data to the attention of the interpreter and the most persistent problem in applying these techniques is the evaluation of potentially interesting patterns. The decision on what use can be made of them is highly interpretive and their calibration is difficult. Well control is vital but it normally allows only very limited supervision of a seismic classifier. An example is presented to illustrate these problems.     

Fitting variogram models by weighted least squares

The method of weighted least squares is shown to be an appropriate way of fitting variogram models. The weighting scheme automatically gives most weight to early lags and down-weights those lags with a small number of pairs. Although weights are derived assuming the data are Gaussian (normal), they are shown to be still appropriate in the setting where data are a (smooth) transform of the Gaussian case. The method of (iterated) generalized least squares, which takes into account correlation between variogram estimators at different lags, offer more statistical efficiency at the price of more complexity. Weighted least squares for the robust estimator, based on square root differences, is less of a compromise.