The perimeter-area fractal model and its application to geology

Perimeters and areas of similarly shaped fractal geometries in two-dimensional space are related to one another by power-law relationships. The exponents obtained from these power laws are associated with, but do not necessarily provide, unbiased estimates of the fractal dimensions of the perimeters and areas. The exponent (D_AL) obtained from perimeter-area analysis can be used only as a reliable estimate of the dimension of the perimeter (D_L) if the dimension of the measured area is D_A=2. If D_A<2, then the exponent D_AL=2D_L/D_A>D_L. Similar relations hold true for area and volumes of three-dimensional fractal geometries. The newly derived results are used for characterizing Au associated alteration zones in porphyry systems in the Mitchell-Sulphurets mineral district, northwestern British Columbia.     

Fractal analysis of the fracture strength of lava dome material based on the grain size distribution of block-and-ash flow deposits at Unzen Volcano, Japan

A fractal theory of rock fragmentation is applied to block-and-ash flow deposits from the Fugendake dome, Unzen Volcano, Kyushu, Japan, in order to analyze the material strength and the energy required for size reduction of the source lava dome. Two fractal dimensions h and D_s, which are mutually interchangeable, represent the relative strength and energy for particles reduced to a given size. They can be theoretically estimated from the power relations of a reference grain size to the cumulative mass and number of fragments smaller than that size. The Unzen–Fugendake block-and-ash flow deposits have been further modified by size sorting and secondary fragmentation that occurred during flowage, so that the h value decreases (or D_s value increases) with increasing distance from the source. Coarse, reversely graded deposits are, however, found to retain the original size population relatively well. The D_s values estimated from deposits of this type are compatible with those previously reported from decompression–fragmentation experiments conducted on the same dome material. The employed fractal approach could thus give insights into the potential mode of dome collapse that generates block-and-ash flows.     

The use of the fractal dimension to quantify the morphology of irregular-shaped particles

For several decades, sedimentologists have had difficulty in obtaining an efficient index of particle form that can be used to specify adequately irregular morphology of sedimentary particles. Mandelbrot has suggested the use of the fractal dimension as a single value estimate of form, in order to characterize morphologically closed loops of an irregular nature. The concept of fractal dimension derives from Richardson's unpublished suggestion that a stable linear relationship appears when the logarithm of the perimeter estimate of an irregular outline is plotted against the logarithm of the unit of measurement (step length). Decreases in step length result in an increase in perimeter by a constant weight (b) for particles whose morphological variations are the same at all measurement scales (self-similarity). The fractal dimension (D) equals 1.0-(b), where b is the slope coefficient of the best-fitting linear regression of the plot. The value of D lies between 1.0 and 2.0, with increasing values of D correlating with increasing irregularity of the outline. In practice, particle outline morphology is not always self-similar, such that two or possibly more fractal elements can occur for many outlines. Two fractal elements reflect the morphological difference between micro-scale edge textural effects (D₁) and macro-scale particle structural effects (D₂) generated by the presence of crenellate-edge morphology (re-entrants). Fractal calibration on a range of regular/irregular particle outline morphologies, plus examination of carbonate beach, pyroclastic and weathered quartz particles indicates that this type of analysis is best suited for morphological characterization of irregular and crenellate particles. In this respect, fractal analysis appears as the complementary analytical technique to harmonic form analysis in order to achieve an adequate specification of all types of particles on a continuum of irregular to regular morphology.     

Application of Fractal Models to Distinguish between Different Mineral Phases

Fractal modeling is demonstrated to be an effective and rapid tool to distinguish between mineral phases in rock samples. It supplements work that previously could be performed only by observing the interpenetrational or metasomatic phenomena between different minerals with the aid of mineralographic microscope. The Gejiu tin district in southwestern China was chosen as a study area for the recognition and characterization of the spatial distribution of two phases (Types I and II) of cassiterite. Vector patterns used for this study were extracted from digital photomicrographs and analyzed with the aid of MapGIS. Perimeter–area fractal dimension, cumulative number–area exponent, and shape index were determined in order to quantify geometrical irregularities and spatial cassiterite phase distribution characteristics. The results show that fractal dimensions based on area and perimeter are larger for crystals of Type I than for those of Type II. The mean shape index (SI) increases from 0.54 (Type I) to 0.64 (Type II), indicating an increase in regularity. The number–area exponent also increases from 0.88 to 1.15, indicating the smaller crystals of Type II. The cumulative number–shape index log–log plot shows two separate straight-line segments. One of these probably represents a background shape realized during the original process of natural crystallization, whereas the other likely represents anomalous shapes because of weathering or other superimposed processes. Two parallel lines can be constructed on the perimeter–area log–log plots. The upper line, with a larger intercept, represents crystals with lower SI. The lower line represents crystals with higher SI, indicating that the intercept provides a measure of the irregularity of grains. By combining the perimeter–area model with cumulative number–area plot and shape index, the two phases of cassiterite can be distinguished and characterized. One phase has fewer crystals of large size, and the other has smaller crystals. This difference can be explained by assuming that under higher-temperature conditions, the large cassiterite crystals formed earlier than the smaller crystals. Consequently, the large cassiterites underwent longer, high-intensive weathering than the small crystals so that their shapes became more irregular. The younger, more abundant small cassiterites retained their original regular shapes.     

Fractal dimension of pore space in carbonate samples from Tushka area (Egypt)

To investigate inhomogeneous and porous structures in nature, the concept of fractal dimension was established. This paper briefly introduces the definition and measurement methods of fractal dimension. Three different methods including mercury injection capillary pressure (MICP), nuclear magnetic resonance (NMR), and nitrogen adsorption (BET) were applied to determine the fractal dimensions of the pore space of eight carbonate rock samples taken from West Tushka area, Egypt. In the case of fractal behavior, the capillary pressure P _c and cumulative fraction V _c resulting from MICP are linearly related with a slope of D-3 in a double logarithmic plot with D being the value of fractal dimension. For NMR, the cumulative intensity fraction V _c and relaxation time T ₂ show a linear relation with a slope of 3-D in a double logarithmic plot. Fractal dimension can also be determined by the specific surface area S _por derived from nitrogen adsorption measurements and the effective hydraulic radius. The fractal dimension D shows a linear relation with the logarithm of S _por . The fractal dimension is also used in models of permeability prediction. To consider a more comprehensive data set, another 34 carbonate samples taken from the same study area were integrated in the discussion on BET method and permeability prediction. Most of the 42 rock samples show a good agreement between measured permeability and predicted permeability if the mean surface fractal dimension for each facies is used.     

Morphometry of Quartz Aggregates in Granites: Fractal Images Referring to Nucleation and Growth Processes

Knowledge of mineral aggregate morphologies is of importance to analyze characteristic differences in rock-forming features. For quantifying these differences, the fractal geometry of quartz aggregate cuts digitized from polish sections of different types of granites has been studied. As an approach to measure fractal dimension (D), a power-law dependence of square of aggregate cuts on their linear size has been used. The D values thus calculated mainly increase from 1.48–1.62 for amazonite granites to 1.63–1.70 for alaskite granites and 1.75–1.81 for standard granites. To account for the data of morphometry, the model of nucleation and growth as applied to silicate melt freezing has been considered. For comparison between the nature and model textures, the fractal properties of cluster cuts in the system of overlapping spheres randomly distributed with random radii have been investigated through computer simulation. It has been demonstrated that the distributions of quartz aggregates in granites may be simulated by homogeneous or heterogeneous Poisson models, and both order of crystallization and metamorphic recrystallization should be taken for explaining textural variability. The results of the simulation have enabled the granitic texture to be discussed with respect to the random configuration of the spatial percolation cluster.     

Seismic hazards assessment of Kumaun Himalaya and adjacent region

The ongoing continent?Ccontinent collision between Indian and Eurasian plates houses a seismic gap in the geologically complex and tectonically active central Himalaya. The seismic gap is characterized by unevenly distributed seismicity. The highly complex geology with equally intricate structural elements of Himalaya offers an almost insurmountable challenge to estimating seismogenic hazard using conventional methods of Physics. Here, we apply integrated unconventional hazard mapping approach of the fractal analysis for the past earthquakes and the box counting fractal dimension of structural elements in order to understand the seismogenesis of the region properly. The study area extends from latitude 28°N?C33°N and longitude 76°E?C81°E has been divided into twenty-five blocks, and the capacity fractal dimension (D ₀) of each block has been calculated using the fractal box counting technique. The study of entire blocks reveal that four blocks are having very low value of D ₀ (0.536, 0.550, 0.619 and 0.678). Among these four blocks two are characterized by intense clustering of earthquakes indicated by low value of correlation fractal dimension (D _c ) (0.245, 0.836 and 0.946). Further, these two blocks are categorized as highly stressed zones and the remaining two are characterized by intense clustering of structural elements in the study area. Based on the above observations, integrated analysis of the D _c of earthquakes and D ₀ of structural elements has led to the identification of diagnostic seismic hazard pattern for the four blocks.     

Collision leading to multiple-stage large-scale extrusion in the Qinling orogen: Insights from the Mianlue suture

The geologic framework of the Phanerozoic Qinling–Dabie orogen was built up through two major suturing events of three blocks. From north to south these include the North China craton (including the north Qinling block), the Qinling–Dabie microblock, and the South China craton (including the Bikou block), separated by the Shangdan and Mianlue sutures. The Mianlue suture zone contains evidence for Mesozoic extrusion tectonics in the form of major strike–slip border faults surrounding basement blocks, a Late Paleozoic ophiolite and a ca. 240–200 Ma thrust belt that reformed by 200–150 Ma thrusts during A-type (intracontinental) subduction. The regional map pattern shows that the blocks are surrounded by complexly deformed Devonian to Early Triassic metasandstones and metapelites, forming a regional-scale block-in-matrix mélange fabric. Five distinct tectonic units have been recognized in the belt: (1) basement blocks including two types of Precambrian basement, crystalline and transitional; (2) continental margin slices including Early Paleozoic strata, and Late Paleozoic fluviodeltaic sedimentary rocks, proximal and distal fan clastics, reflecting the development of a north-facing rift margin on the edge of the South China plate; (3) out of sequence oceanic crustal slices including strongly deformed postrift, deep-water sedimentary rocks, sheeted dikes, basalts, and mafic–ultramafic cumulates of a Late Paleozoic ophiolite suite, developing independent of the rift margin in a separate basin; (4) out-of-sequence island-arc slices; (5) accretionary wedge slices. All the tectonic units were deformed during three geometrically distinct deformation episodes (D₁, D₂ and D₃ during 240–200 Ma). Units 2–4 involved southward thrusting and vertical then southward extrusion of about 20 km of horizontal displacement above the autochthonous basement during the D₁ episode. Thrust slices 20 km south of the Mianlue suture are related to this vertical extrusion due to the same rock assemblages, ages and kinematics. The D₂ and D₃ episodes folded all the units in a thick-skinned style about east–west (D₂) and west–northwest (D₃) axes in the Mianlue suture zone. An early foreland propagating sequence of accretion of Late Paleozoic rocks deposited above the Yangtze craton is not involved in D₁ deformation but is temporally equivalent to the D₂ and D₃ deformation in the Mianlue suture. Two stages of strike–slip faulting mainly occurred at the end of D₂ and D₃, respectively. During D₂ deformation, the Bikou block was obliquely indented to the ESE into the Mianlue suture, rather than being thrust over the Mianlue suture from the north as a part of the Qinling–Dabie microblock. During D₃ deformation, however, the Bikou block was bounded by the south boundary fault of the Mianlue suture, and the Yangpingguan fault on the south. These faults are coeval strike–slip faults, but of opposite senses, and accommodated minor southwestward extrusion of the Bikou block into Songpan–Ganze orogen. The other basement blocks north of the Mianlue suture were extruded eastward by about 20 km of lateral displacement, based on the offset of the Wudang dome, during the D₃ episode due to the northeastward indentation of the Hannan complex of the South China craton. Post-D₃ emplacement of granite, cutting across the strike–slip faults such as the Mianlue suture, provides a minimum age of 200 Ma for D₃ deformation. Therefore, based on insights from the evolution of the Mianlue suture, the D₂ and D₃ episodes in the Mianlue suture and its neighbors are not responsible for and associated with the two-stage extrusion of the Dabie UHP-HP terranes from the Foping dome to the present erosional surface (more than 350 km).     

The determination of trace elements in Fe–Mn oxide coatings on pebbles using LA-ICP-MS

Iron–manganese oxide coatings form on a wide range of geologic samples where they have the ability to adsorb elements and potentially act as a mineral exploration/environmental monitoring tool. In this study, Fe–Mn oxide coatings on stream pebbles were collected from streams in four study areas located across the province of Newfoundland and Labrador, Canada. The study locations were in areas of former copper mines (Tilt Cove and Betts Cove), carbonate geology (Robinsons River), and a metropolitan area (Rennies River). Collected pebbles underwent a simple sample preparation procedure and were then analyzed for a wide range of elements by LA-ICP-MS after optimization of the operating conditions. Water samples accompanied the pebbles, and these were analyzed for pH, dissolved oxygen, conductivity, and a large selection of elements by ICP-MS. Multivariate statistics, in the form of Principal Component Factor Analysis (PCFA) was performed on both data sets. Graphs of the factor scores from the PCFA produced groupings of the samples that were related to geologic/environmental inputs. The loading of variables in each factor was related to the adsorption of the element either to the MnO₂ or Fe₂O₃ phase with most elements except Cr and Cu displaying preferential adsorption to MnO₂. Elemental Fe–Mn oxide coating concentrations were a result of the element's affinity (chalcophile, lithophile, or siderophile), pH of the environment, stream water concentration, and amount of each oxide phase present. Even with these complications, LA-ICP-MS analysis of Fe–Mn oxides was able to identify areas of heavy metal pollution and locate geologic inputs.     

Cordierite gneisses and associated lithologies of the Guri area,Northwest Guayana Shield,Venezuela

Gneisses in the Guri area of the Venezuelan Guayana Shield contain mineral assemblages with cordierite, garnet, sillimanite, hypersthene, biotite and Fe-Ti oxide intergrowths.Analysis of mineral assemblages and compositional relationships in the light of experimental data indicate metamorphic conditions of 725–800° C, 5–6 kb P _T , <P _T for the highest grade rocks and 650–700° C, 5–7 kb P _T , approximating P _T for the lowest grade rocks. Oxygen fugacities in different lithologies ranged between those of the MH and QFM buffers.The distribution coefficient K _D (Mg-Fe) (gar-bio), decreases by 0.006 per atom percent increase in (Mn/Mn+Mg+Fe)gar, falls in the range of K _D typical of the sillimanite-K feldspar zone and granulite facies, and is systematically lower in lower grade rocks-all in accord with observation in other localities. K _D (Mg-Fe) (cord-bio) ranges from 3.0 in the highest grade rocks to 10.0 in the lowest grade rocks, appears independent of FeO/MgO of cordierite or biotite, and varies systematically with grade. In contrast with conclusions based on observation in other localities, data from the Guri area suggest -K_D(cord-bio) may be a sensitive index of grade.A number of mineralogic and geologic observations are difficultly reconciled with existing experimental data.     

Timing of deformation in the Norseman-Wiluna Belt, Yilgarn Craton, Western Australia

Establishing relative and absolute time frameworks for the sedimentary, magmatic, tectonic and gold mineralisation events in the Norseman-Wiluna Belt of the Archean Yilgarn Craton of Western Australia, has long been the main aim of research efforts. Recently published constraints on the timing of sedimentation and absolute granite ages have emphasized the shortcomings of the established rationale used for interpreting the timing of deformation events. In this paper the assumptions underlying this rationale are scrutinized, and it is shown that they are the source of significant misinterpretations. A revised time chart for the deformation events of the belt is established. The first shortening phase to affect the belt, D₁, was preceded by an extensional event D_1e and accompanied by a change from volcanic-dominated to plutonic-dominated magmatism at approximately 2685–2675 Ma. Later extension (D_2e) controlled deposition of the ca 2655 Ma Kurrawang Sequence and was followed by D₂, a major shortening event, which folded this sequence. D₂ must therefore have started after 2655 Ma—at least 20 Ma later than previously thought and after the voluminous 2670–2655 Ma high-Ca granite intrusion. Younger transcurrent deformation, D₃–D₄, waned at around 2630 Ma, suggesting that the crustal shortening deformation cycle D₂–D₄ lasted approximately 20–30 Ma, contemporaneous with low-volume 2650–2630 Ma low-Ca granites and alkaline intrusions. Time constraints on gold deposits suggest a late mineralisation event between 2640–2630 Ma. Thus, D₂–D₄ deformation cycle and late felsic magmatism define a 20–30 Ma long tectonothermal event, which culminated with gold mineralisation. The finding that D₂ folding took place after voluminous high-Ca granite intrusion led to research into the role of competent bodies during folding by means of numerical models. Results suggest that buoyancy-driven doming of pre-tectonic competent bodies trigger growth of antiforms, whereas non-buoyant, competent granite bodies trigger growth of synforms. The conspicuous presence of pre-folding granites in the cores of anticlines may be a result from active buoyancy doming during folding.     

Mineralogy and geochemistry of boehmite-rich coals: New insights from the Haerwusu Surface Mine, Jungar Coalfield, Inner Mongolia, China

Boehmite-rich coal of Pennsylvanian age was discovered earlier at the Heidaigou Surface Mine, Jungar Coalfield, Inner Mongolia, China. This paper reports new results on 29 bench samples of the no. 6 coal from a drill core from the adjacent Haerwusu Surface Mine, and provides new insights into the origin of the minerals and elements present. The results show that the proportion of inertinite in the no. 6 coal is higher than in other Late Paleozoic coals in northern China. Based on mineral proportions (boehmite to kaolinite ratio) and major element concentrations in the coal benches of the drill core, the no. 6 coal may be divided into five sections (I to V). Major minerals in Sections I and V are kaolinite. Sections II and IV are mainly kaolinite with a trace of boehmite, and Section III is high in boehmite. The boehmite is derived from bauxite in the weathered surface (Benxi Formation) in the sediment-source region. The no. 6 coal is rich in Al₂O₃ (8.89%), TiO₂ (0.47%), Li (116 μg/g), F (286 μg/g), Ga (18 μg/g), Se (6.1 μg/g), Sr (350 μg/g), Zr (268 μg/g), REEs (172 μg/g), Pb (30 μg/g), and Th (17 μg/g). The elements are classified into five associations by cluster analysis, i.e. Groups A, B, C, D, and E. Group A (ash–SiO₂–Al₂O₃–Na₂O–Li) and Group B (REE–Sc–In–Y–K₂O–Rb–Zr–Hf–Cs–U–P₂O₅–Sr–Ba–Ge) are strongly correlated with ash yield and mainly have an inorganic affinity. The elements that are negatively or less strongly correlated with ash yield (with exceptions of Fe₂O₃, Be, V, and Ni) are grouped in the remaining three associations: Group C, Se–Pb–Hg–Th–TiO₂–Bi–Nb–Ta–Cd–Sn; Group D, Co–Mo–Tl–Be–Ni–Sb–MgO–Re–Ga–W–Zn–V–Cr–F–Cu; and Group E, S–As–CaO–MnO–Fe₂O₃. Aluminum is mainly distributed in boehmite, followed by kaolinite. The high correlation coefficients of the Li–ash, Li–Al₂O₃, and Li–SiO₂ pairs indicate that Li is related to the aluminosilicates in the coal. The boehmite-rich coal is high in gallium and F, which occur in boehmite and the organic matter. Selenium and Pb are mainly in epigenetic clausthalite fillings in fractures. The abundant rare earth elements in the coal benches were supplied from two sources: the bauxite on the weathered surface of the Benxi Formation and from adjacent partings by groundwater leaching during diagenesis. The light rare earth elements (LREEs) are more easily leached from the partings and incorporated into the organic matter than the heavy REEs, leading to a higher ratio of LREEs to HREEs in the coal benches than in the overlying partings.     

Partitioning of Cu,Sn, Mo,W, U,and Th between melt and aqueous fluid in the systems haplogranite-H2O−HCl and haplogranite-H2O−HF

The partition coefficients K_D=c_fluid/c_melt of Cu, Sn, Mo, W, U, and Th between aqueous fluid and melt were measured in the systems haplogranite-H₂O–HCl and haplogranite-H₂O–HF at 2kbars, 750°C, and Ni–NiO buffer conditions using rapid-quench cold seal bombs, with many reversed runs. Concentrations of trace elements (1–1000 ppm) in the quenched aqueous fluid and in the glass were determined by plasma emission spectrometry (DCP). K_D of F is close to 1 in the system studied. K_D of Cu and Sn strongly increases with increasing Cl concentration due to the formation of chloride complexes in the aqueous fluid, while HF has no effect. However, in 2M HCl, K_D of Cu approaches 100, while K_D of Sn is below 0.1 under the same conditions. The partition coefficients of Mo and W are high if water is the only volatile present (Mo: 5.5, W: 3.5), but strongly decrease with increasing HCl and HF, due to the destabilization of hydroxy complexes. K_D of U and Th is very low in the absence of complexing agents, but strongly increases with increasing HF concentration. K_D of U also increases with increasing HCl concentration and with increasing CO₂ concentration in the system haplogranite-H₂O–CO₂, indicating the stability of chloride and carbonate complexes of U at magmatic temperatures. The data suggest a stoichiometric ratio of Cl: U=3:1 and of F:U=2:1 in these complexes. Cl-rich fluids are responsible for the formation of porphyry Cu deposits, but are much less effective in the transport of Sn. F appears not to be essential for the concentration of Mo and W in fluids evolving from a granitic magma. The different complexing behavior of U and Th in aqueous fluids may account for their fractionation during magma genesis.     

Geology and geochemistry of bauxite deposits in Lushoto District, Usambara Mountains, Tanzania

Bauxite deposits in the Usambara Mountains of north eastern Tanzania occur as remnants of residual deposits on two geomorphologically related plateaus of Mabughai-Mlomboza and Kidundai at Magamba in Lushoto, Usambara Mountains. The parent rocks for the deposits are mainly granulites and feldspathic gneisses of Neoproterozoic Mozambique belt. The plateaus represent a preserved Late Cretaceous–Lower Tertiary old land surface (African surface). Other parts of the Usambara Mountains and the neighbouring Pare Mountains are covered mostly by red–brown lateritic soils and impure reddish-brown kaolinitic clays. The bauxite deposits contain mainly Al₂O₃ (40–69 wt.%), Fe₂O₃ (3–10 wt.%), SiO₂ (0.16–7 wt.%) and other elements occur in quantities not substantial to affect the quality or processing of the bauxite, and are attributed to the presence of relic minerals. Gibbsite makes up to 98 vol.% of the bauxite ore in special cases. Gibbsite is accompanied by goethite in the ore. Boehmite occurs in small amounts and is usually accompanied by hematite. Impurities include goethite, hematite, kaolinite, and minor relic quartz and microcline. Kaolinite is the sole clay mineral encountered in the bauxite ore, suggesting mature soil profiles and a development of the bauxite deposits on a well-drained peneplanation. Ore reserve estimates from the drilling data and surface geological mapping of the deposits yielded bauxite reserves of about 37 million tonnes.     

Magmatic evolution of the differentiated ultramafic, alkaline and carbonatite intrusion of Vuoriyarvi (Kola Peninsula, Russia). A LA-ICP-MS study of apatite

The nature of the petrogenetic links between carbonatites and associated silicate rocks is still under discussion (i.e., [Gittins J., Harmer R.E., 2003. Myth and reality of the carbonatite–silicate rock “association”. Period di Mineral. 72, 19–26.]). In the Paleozoic Kola alkaline province (NW Russia), the carbonatites are spatially and temporally associated to ultramafic cumulates (clinopyroxenite, wehrlite and dunite) and alkaline silicate rocks of the ijolite–melteigite series [(Kogarko, 1987), (Kogarko et al., 1995), (Verhulst et al., 2000), (Dunworth and Bell, 2001) and (Woolley, 2003)]. In the small (≈ 20 km²) Vuoriyarvi massif, apatite is typically a liquidus phase during the magmatic evolution and so it can be used to test genetic relationships. Trace elements contents have been obtained for both whole rocks and apatite (by LA-ICP-MS). The apatites define a single continuous chemical evolution marked by an increase in REE and Na (belovite-type of substitution, i.e., 2Ca²⁺ = Na⁺ + REE³⁺). This evolution possibly reflects a fractional crystallisation process of a single batch of isotopically homogeneous, mantle-derived magma.The distribution of REE between apatite and their host carbonatite have been estimated from the apatite composition of a carbonatite vein, belonging to the Neskevara conical-ring-like vein system. This carbonatite vein is tentatively interpreted as a melt. So, the calculated distribution coefficients are close to partition coefficients. Rare earth elements are compatible in apatite (D > 1) with a higher compatibility for the middle REE (D_Sm : 6.1) than for the light (D_La : 4.1) and the heavy (D_Yb : 1) REE.     

Genetic significance of fluid inclusions in the CSA Cu–Pb–Zn deposit, Cobar, Australia

CSA mine exploits a ‘Cobar-type’ Cu–Pb–Zn±Au±Ag deposit within a cleaved and metamorphosed portion of the Cobar Supergroup, central New South Wales. The deposit comprises systems of ‘lenses’ that encompass veins, disseminations and semi-massive to massive Cu–Pb–Zn ores. The systems and contained lenses truncate bedding, are approximately coplanar with regional cleavage and similarly oriented shear zones and plunge parallel to the elongation lineation. Systems have extreme vertical continuity (>1000 m), short strike length (400 m) and narrow width (100 m), exhibit vertical and lateral ore-type variation and have alteration haloes. Models of ore formation include classical hydrothermalism, structurally controlled remobilisation and polymodal concepts; syntectonic emplacement now holds sway.Fluid inclusions were examined from quartz±sulphide veins adjacent to now-extracted ore, from coexisting quartz–sulphide within ore, and from vughs in barren quartz veins. Lack of early primary inclusions precluded direct determination of fluids associated with D₂–D₃ ore and vein emplacement. Similarly, decrepitation (by near-isobaric heating) of the two oldest secondary populations precluded direct determination of fluid phases immediately following D₂–D₃ ore and vein emplacement. Post-decrepitation outflow (late D₃ to early post-D₃) is recorded by monophase CH₄ inclusions. Entrained outflow of deeply circulated meteoric fluid modified the CH₄ system; modification is recorded by H₂O+CH₄ and H₂O+(trace CH₄) secondary populations and by an H₂O+(trace CH₄) primary population. The contractional tectonics (D₂–D₃) of ore emplacement was superseded by relaxational tectonics (D_4P) that facilitated meteoric water penetration and return flow.Under D₂ prograde metamorphism, entrapment temperatures (Tt) and pressures (Pt) for pre-decrepitation secondary inclusions are estimated as Tt300–330 °C and Pt1.5–2 kbar≈Plith (the lithostatic pressure). Decrepitation accompanied peak metamorphism (T350–380 °C) in mid- to late-D₃, while in late-D₃ to early post-D₃, essentially monophase CH₄ secondary inclusions were entrapped at Tt350 °C and Pt=1.5–2 kbar≈Plith. Subsequently, abundant CH₄ and entrained meteoric water were entrapped as H₂O+CH₄ secondaries under slowly decreasing temperature (Tt330–350 °C) and constant pressure (Pt1.5–2 kbar). Finally, with increasingly dominant meteoric outflow, H₂O+(trace CH₄) populations record decreasing temperatures (Tt>300 to <350 down to 275–300 °C) at pressures of Phydrostatic<Pt (1 kbar) <Plith (1.5 kbar).The populations of inclusions provide insight into fluid types, flow regimes and P–T conditions during parts of the deposit's evolution. They indirectly support the role of basin-derived CH₄ fluids in ore formation, but provide no insight into a basement-sourced ore-forming fluid. They fully support post-ore involvement of meteoric water. The poorly constrained entrapment history is believed to span 10 Ma from 395 to 385 Ma.     

Gold mineralisation throughout about 45 Ma of Archaean orogenesis: protracted flux of gold in the Golden Mile, Yilgarn craton, Western Australia

The Golden Mile deposit was discovered in 1893 and represents today the largest Archaean orogenic lode gold system in the world (50 M oz produced gold). The Golden Mile deposit comprises three major styles of gold mineralisation: Fimiston, Oroya and Charlotte styles. Fimiston-style lodes formed at 250 to 350 °C and 100 to 200 MPa and are controlled by brittle–ductile fault zones, their subsidiary fault zone and vein networks including breccias and open-cavity-infill textures and hydrothermally altered wall rock. Fimiston lodes were formed late D₁, prior to D₂ regional upright folding. Hydrothermal alteration haloes comprise a progression toward the lode of diminishing chlorite, an increase in sericite and in Fe content of carbonates. Lodes contain siderite, pyrite, native gold, 17 different telluride minerals (Au–Ag tellurides contain ~25% of total gold), tourmaline, haematite, sericite and V-rich muscovite. Oroya-style lodes formed at similar P–T conditions as the Fimiston lodes and are controlled by brittle–ductile shear zones, associated dilational jogs that are particularly well developed at the contact between Paringa Basalt and black shale interflow sedimentary rocks and altered wall rock. The orebodies are characterised by micro-breccias and zones of intense shear zone foliation, very high gold grades (up to 100,000 g/t Au) and the common association of tellurides and vanadian mica (green leader). Oroya lodes crosscut Fimiston lodes and are interpreted to have formed slightly later than Fimiston lodes as part of one evolving hydrothermal system spanning D₁ and D₂ deformation (ca. 2,675–2,660 Ma). Charlotte-style lodes, exemplified by the Mt Charlotte deposit, are controlled by a sheeted vein (stockwork) complex of north-dipping quartz veins and hydrothermally altered wall rock. The Mt Charlotte orebody formed at 120 to 440 °C and 150 to 250 MPa during movement along closely spaced D₄ (2,625 Ma) and reactivated D₂ faults with the quartz granophyre in the Golden Mile Dolerite exerting a strong lithological control on gold mineralisation. Veins consist of quartz–carbonate–minor scheelite, and wall-rock alteration comprises chlorite destruction and growth of ferroan carbonate–sericite–pyrite–native gold. Pyrite–pyrrhotite is zoned on the scale of vein haloes and of the entire mine, giving a vertical temperature gradient of 50–100 °C over 1,000 vertical metres. The structural–hydrothermal model proposed consists of four major stages: (1) D₁ thrusting and formation of Fimiston-style lodes, (2) D₂ reverse faulting and formation of Oroya-style lodes, (3) D₃ faulting and dissecting of Fimiston- and Oroya-style lodes, and (4) D₄ faulting and formation of Mt Charlotte-style sheeted quartz vein system. The giant accumulation of gold in the Golden Mile deposit was formed due to protracted gold mineralisation throughout episodes of an Archaean orogeny that spanned about 45 Ma. Fluid conduits formed early in the tectonic history and persisted throughout orogenesis with the plumbing system showing a rare high degree of focussing, efficiency and duration. In addition to the long-lasting fluid plumbing system, the wide variety of transient structural and geochemical traps, multiple fluid sources and precipitation mechanism contributed towards the richest golden mile in the world.Editorial handling: B. Lehmann     

Tectonomagmatism in continental arcs: evidence from the Sark arc complex

The island of Sark (Channel Islands, UK) exposes syntectonic plutons and country rock gneisses within a Precambrian (Cadomian) continental arc. This Sark arc complex records sequential pulses of magmatism over a period of 7 Ma (ca. 616–609 Ma). The earliest intrusion (ca. 616 Ma) was a composite sill that shows an ultramafic base overlain by a magma-mingled net vein complex subsequently deformed at near-solidus temperatures into the amphibolitic and tonalitic Tintageu banded gneisses. The deformation was synchronous with D₂ deformation of the paragneissic envelope, with both intrusion and country rock showing flat, top-to-the-south LS fabrics. Later plutonism injected three homogeneous quartz diorite–granodiorite sheets: the Creux–Moulin pluton (150–250 m; ca. 614 Ma), the Little Sark pluton (>700 m; 611 Ma), and the Northern pluton (>500 m; 609 Ma). Similar but thinner sheets in the south (Derrible–Hogsback–Dixcart) and west (Port es Saies–Brecqhou) are interpreted as offshoots from the Creux–Moulin pluton and Little Sark pluton, respectively. All these plutons show the same LS fabric seen in the older gneisses, with rare magmatic fabrics and common solid state fabrics recording syntectonic crystallisation and cooling. The cooling rate increased rapidly with decreasing crystallisation age: >9 Ma for the oldest intrusion to cool to lower amphibolite conditions, 7–8 Ma for the Creux Moulin pluton, 5–6 Ma for the Little Sark pluton, and <3 Ma for the Northern pluton. This cooling pattern is interpreted as recording extensional exhumation during D₂. The initiation of the D₂ event is suggested to have been a response to the intrusion of the Tintageu magma which promoted a rapid increase in strain rate (>10⁻¹⁴ s⁻¹) that focussed extensional deformation into the Sark area. The increased rates of extension allowed ingress of the subsequent quartz diorite–granodiorite sheets, although strain rate slowly declined as the whole complex cooled during exhumation. The regional architecture of syntectonic Cadomian arc complexes includes flat-lying “Sark-type” and steep “Guernsey-type” domains produced synchronously in shear zone networks induced by oblique subduction: a pattern seen in other continental arcs such as that running from Alaska to California.     

Structural control and hydrothermal alteration at the BIF-hosted Raposos lode-gold deposit, Quadrilátero Ferrífero, Brazil

In the Raposos orogenic gold deposit, hosted by banded iron-formation (BIF) of the Archean Rio das Velhas greenstone belt, the hanging wall rocks to BIF are hydrothermally-altered ultramafic schists, whereas metamafic rocks and their hydrothermal schistose products represent the footwall. Planar and linear structures at the Raposos deposit define three ductile to brittle deformational events (D₁, D₂ and D₃). A fourth group of structures involve spaced cleavages that are considered to be a brittle phase of D₃. The orebodies constitute sulfide-bearing D₁-related shear zones of BIF in association with quartz veins, and result from the sulfidation of magnetite and/or siderite. Pyrrhotite is the main sulfide mineral, followed by lesser arsenopyrite and pyrite. At level 28, the hydrothermal alteration of the mafic and ultramafic wall rocks enveloping BIF define a gross zonal pattern surrounding the ore zones. Metabasalt comprises albite, epidote, actinolite and lesser Mg/Fe–chlorite, calcite and quartz. The incipient stage includes the chlorite and chlorite-muscovite alteration zone. The least-altered ultramafic schist contains Cr-bearing Mg-chlorite, actinolite and talc, with subordinate calcite. The incipient alteration stage is subdivided into the talc–chlorite and chlorite–carbonate zone. For both mafic and ultramafic wall rocks, the carbonate–albite and carbonate–muscovite zones represent the advanced alteration stage.Rare earth and trace element analyses of metabasalt and its alteration products suggest a tholeiitic protolith for this wall rock. In the case of the ultramafic schists, the precursor may have been peridotitic komatiite. The Eu anomaly of the Raposos BIF suggests that it was formed proximal to an exhalative hydrothermal source on the ocean floor. The ore fluid composition is inferred by hydrothermal alteration reactions, indicating it to having been H₂O-rich containing CO₂ + Na⁺ and S. Since the distal alteration halos are dominated by hydrated silicate phases (mainly chlorite), with minor carbonates, fixation of H₂O is indicated. The CO₂ is consumed to form carbonates in the intermediate alteration stage, in halos around the chlorite-dominated zones. These characteristics suggest variations in the H₂O to CO₂-ratio of the sulfur-bearing, aqueous-carbonic ore fluid, which interacted at varying fluid to rock ratios with progression of the hydrothermal alteration.     

The effects of surface area, grain size and mineralogy on organic matter sedimentation and preservation across the modern Squamish Delta, British Columbia: the potential role of sediment surface area in the formation of petroleum source rocks

Surface sediment samples were collected from the Squamish River Delta, British Columbia, in order to determine the role of sediment surface area in the preservation of organic matter (OM) in a paralic sedimentary environment. The Squamish Delta is an actively prograding delta, located at the head of Howe Sound.Bulk total organic carbon (TOC) values across the Squamish Delta are low, ranging from 0.1 to 1.0 wt.%. The carbon/total nitrogen ratio (C_org/N) ranges from 6 to 17, which is attributed to changes in OM type and facies variations. The <25-μm fraction has TOC concentrations up to 2.0 wt.%, and a C_org/N ratio that ranges from 14 to 16. The 53–106-μm fraction has higher TOC concentrations and C_org/N ratios relative to the 25–53-μm fraction. The C_org/N ratio ranges from 9 to 18 in the 53–106-μm fraction and 5.5–10.5 in the 25–53-μm fraction. Surface area values for bulk sediments are low (0.5–3.0 m²/g) due to the large proportion of silt size material. Good correlation between surface area and TOC in bulk samples suggests that OM is adsorbed to mineral surfaces. Similar relationships between surface area and TOC were observed in size-fractionated samples. Mineralogy and elemental composition did not correlate with TOC concentration.The relationships between surface area, TOC and total nitrogen (TN) can be linked to the hydrodynamic and sedimentological conditions of the Squamish Delta. As a result, the Squamish Delta is a useful modern analogue for the formation of petroleum source rocks in ancient deltaic environments, where TOC concentrations are often significantly lower than those in source rocks formed in other geological settings.