Adakite-like signature of Late Miocene intrusions at the Los Pelambres giant porphyry copper deposit in the Andes of central Chile: metallogenic implications

Adakite-like features are recognized in the Late Miocene (~10 Ma) porphyritic intrusions of the Los Pelambres giant porphyry copper deposit, central Chile (32°S). Located within the southern portion of the flat-slab segment (28–33°S) of the Chilean Andes, the Al- and Na-rich porphyries of Los Pelambres display distinctly higher Sr/Y (~100–300) and La_N/Yb_N (~25–60) ratios than contemporaneous and barren magmatic units (e.g., La Gloria pluton, Cerro Aconcagua volcanic rocks) of the same Andean magmatic belt. Strong fractionation of heavy rare earth elements (HREE), absence of Eu anomalies, high Sr/Y and Zr/Sm and low Nb/Ta ratios suggest melt extraction from a garnet-amphibolite source. The Late-Miocene adakite-like porphyritic intrusions at Los Pelambres formed closely related in time and space to the subduction of the Juan Fernández Ridge (JFR) hotspot chain along the Chilean margin. Current tectonic reconstructions reveal that, at the time of formation of the Los Pelambres rocks, a W-E segment of the JFR started to subduct beneath them, producing a slow-down of a previously rapid southward migration of a NE-ridge—trench collision. These particular tectonic conditions are favorable for the origin of the Los Pelambres porphyry suite by melting of subducting young hotspot rocks under flat-slab conditions. The incorporation of crustal components into the oceanic lithopheric magma source by subduction erosion is evidenced by the Sr-Nd isotope composition of the Los Pelambres rocks different from the MORB signatures of true adakites. A close relationship apparently exists between the origin of this adakite-like magmatism and the source of the mineralization in the Los Pelambres porphyry copper deposit.Editorial handling: R.J. Goldfarb     

Stable isotope (C,O, S) systematics of the mercury mineralization at Idrija,Slovenia: constraints on fluid source and alteration processes

The world-class Idrija mercury deposit (western Slovenia) is hosted by highly deformed Permocarboniferous to Middle Triassic sedimentary rocks within a complex tectonic structure at the transition between the External Dinarides and the Southern Alps. Concordant and discordant mineralization formed concomitant with Middle Triassic bimodal volcanism in an aborted rift. A multiple isotopic (C, O, S) investigation of host rocks and ore minerals was performed to put constraints on the source and composition of the fluid, and the hydrothermal alteration. The distributions of the ¹³C and ¹⁸O values of host and gangue carbonates are indicative of a fracture-controlled hydrothermal system, with locally high fluid-rock ratios. Quantitative modeling of the ¹³C and ¹⁸O covariation for host carbonates during temperature dependent fluid-rock interaction, and concomitant precipitation of void-filling dolomites points to a slightly acidic hydrothermal fluid (¹³C–4 and ¹⁸O+10), which most likely evolved during isotopic exchange with carbonates under low fluid/rock ratios. The ³⁴S values of hydrothermal and sedimentary sulfur minerals were used to re-evaluate the previously proposed magmatic and evaporitic sulfur sources for the mineralization, and to assess the importance of other possible sulfur sources such as the contemporaneous seawater sulfate, sedimentary pyrite, and organic sulfur compounds. The ³⁴S values of the sulfides show a large variation at deposit down to hand-specimen scale. They range for cinnabar and pyrite from –19.1 to +22.8, and from –22.4 to +59.6, respectively, suggesting mixing of sulfur from different sources. The peak of ³⁴S values of cinnabar and pyrite close to 0 is compatible with ore sulfur derived dominantly from a magmatic fluid and/or from hydrothermal leaching of basement rocks. The similar stratigraphic trends of the ³⁴S values of both cinnabar and pyrite suggest a minor contribution of sedimentary sulfur (pyrite and organic sulfur) to the ore formation. Some of the positive ³⁴S values are probably derived from thermochemical reduction of evaporitic and contemporaneous seawater sulfates.Editorial handling: P. Lattanzi     

Composite synvolcanic intrusions associated with Precambrian VMS-related hydrothermal systems

Large subvolcanic intrusions are recognized within most Precambrian VMS camps. Of these, 80% are quartz diorite–tonalite–trondhjemite composite intrusions. The VMS camps spatially associated with composite intrusions account for >90% of the aggregate sulfide tonnage of all the Precambrian, intrusion-related VMS camps. These low-alumina, low-K, and high-Na composite intrusions contain early phases of quartz diorite and tonalite, followed by more voluminous trondhjemite. They have a high proportion of high silica (>74% SiO ₂) trondhjemite which is compositionally similar to the VMS-hosting rhyolites within the volcanic host-rock successions. The quartz-diorite and possibly tonalite phases follow tholeiitic fractionation trends whereas the trondhjemites fall within the composition field for primitive crustal melts. These transitional M-I-type primitive intrusive suites are associated with extensional regimes within oceanic-arc environments. Subvolcanic composite intrusions related to the Archean Sturgeon Lake and Noranda, and Paleoproterozoic Snow Lake VMS camps range in volume from 300 to 1,000 km ³. Three have a sill morphology with strike lengths between 15 and 22 km and an average thickness between 1,500 and 2,000 m. The fourth has a gross stock-like shape. The VMS deposits are principally restricted to the volcanic strata above the strike length of the intrusions, as are areally extensive, thin exhalite units. The composite intrusions contain numerous internal phases which are commonly clustered within certain parts of the composite intrusion. These clusters underlie eruptive centers surrounded by areas of hydrothermal alteration and which contain most of the VMS deposits. Early quartz-diorite and tonalite phases appear to have intruded in rapid succession. Evidence includes gradational contacts, magma mixing and disequilibrium textures. They appear to have been emplaced as sill-dike swarms. These early phases are present as pendants and xenoliths within later trondhjemite phases. The trondhjemite phases contain numerous internal contacts indicating emplacement as composite sills. Common structural features of the composite intrusions include early xenolith phases, abundant small comagmatic dikes, fractures and veins and, in places, columnar jointing. Internal phases may differ greatly in texture from fine- to coarse-grained, aphyric and granophyric through seriate to porphyritic. Mineralogical and isotopic evidence indicates that early phases of each composite intrusion are affected by pervasive to fracture-controlled high-temperature (350–450 °C) alteration reflecting seawater-rock interaction. Trondhjemite phases contain hydrothermal-magmatic alteration assemblages within miarolitic cavities, hydrothermal breccias and veins. This hydrothermal-magmatic alteration may, in part, be inherited from previously altered wall rocks. Two of the four intrusions are host to Cu-Mo-rich intrusive breccias and porphyry-type mineralization which formed as much as 14 Ma after the main subvolcanic magmatic activity. The recognition of these Precambrian, subvolcanic composite intrusions is important for greenfields VMS exploration, as they define the location of thermal corridors within extensional oceanic-arc regimes which have the greatest potential for significant VMS mineralization. The VMS mineralization may occur for 2,000 m above the intrusions. In some cases, VMS mineralization has been truncated or enveloped by late trondhjemite phases of the composite intrusions. Evidence that much of the trondhjemitic magmatism postdates the principal VMS activity is a critical factor when developing heat and fluid flow models for these subseafloor magmatic-hydrothermal systems.     

Regional-scale hydrothermal alteration in the Central Blake River Group,western Abitibi subprovince,Canada: implications for VMS prospectivity

The Late Archean Blake River Group is a thick succession of predominantly mafic volcanic rocks within the southern zone of the Abitibi greenstone belt. It contains a number of silicic volcanic centers of different size, including the large Noranda volcanic complex, which is host to 17 past-producing volcanogenic massive sulfide deposits. The Noranda complex consists of a 7- to 9-km-thick succession of bimodal mafic and felsic volcanic rocks erupted during five major cycles of volcanism. Massive sulfide formation coincided with a period of intense magmatic activity (cycle III) and the formation of the Noranda cauldron. Hydrothermal alteration in these rocks is interpreted to reflect large-scale hydrothermal fluid flow associated with rapid crustal extension and rifting of the volcanic complex. The alteration includes abundant albite, chlorite, epidote and quartz (silicification), which exhibit broad stratigraphic and structural control and correlate with previously mapped whole-rock oxygen isotope zonation. The Mine Sequence volcanic rocks are characterized by abundant iron-rich chlorite (Fe/Fe+Mg >0.5), hydrothermal amphibole (ferroactinolite) and coarse-grained epidote of clinozoisite composition (<10 wt% Fe ₂O ₃). Volcanic rocks of the pre-cauldron sequences, which contain only subeconomic stringer mineralization, are characterized by less abundant chlorite and mainly fine-grained epidote (>10 wt% Fe ₂O ₃) lacking the clinozoisite solid solution. Alteration in the Mine Sequence volcanic rocks persists along strike well beyond the limits of the main ore deposits (as far as several tens of kilometers) and can be readily distinguished from greenschist facies metamorphic assemblages at a regional scale. The lack of similar alteration in the pre-cauldron sequences is consistent with limited ¹⁸O-depletion and suggests that the early history of the volcanic complex did not support large-scale, high-temperature fluid flow in these rocks. Comparisons with a much smaller, barren volcanic complex in nearby Ben Nevis township reveal important differences in the alteration mineralogy between volcanoes of different size, with implications for area selection during regional-scale mineral exploration. The Ben Nevis Complex consists of a 3- to 4-km-thick succession of mafic, intermediate and felsic volcanic rocks centered on a small subvolcanic intrusion. Alteration of the volcanic rocks comprises mainly low-temperature assemblages of prehnite, pumpellyite, magnesium-rich chlorite (Fe/Fe+Mg <0.5), iron-rich epidote (>10 wt% Fe ₂O ₃) and calcite. Actinolite ± magnetite alteration occurs proximal to the intrusive core of the complex, but the limited extent of this alteration indicates only local high-temperature fluid circulation adjacent to the intrusion. A distal zone of carbonate alteration is located 4–6 km from the center of the volcano. Although iron-bearing carbonates are present locally within this zone, the absence of siderite argues against a high-temperature origin for this alteration. These observations do not offer positive encouragement for the existence of a fossil geothermal system of sufficient size or intensity to have produced a large massive sulfide deposit.     

Geochemistry and mineralogy of Platinum-group elements in the Ransko gabbro–peridotite massif,Bohemian Massif (Czech Republic)

The Ransko gabbro-peridotite massif in Eastern Bohemia is a strongly differentiated intrusive complex of Lower Cambrian age. The complex hosts low grade Ni-Cu ores mainly developed close to the contact of olivine-rich rocks with gabbros, in troctolites and, to a much lesser extent, in both pyroxene and olivine gabbros and plagioclase-rich peridotites. The ore zone is characterized by strong serpentinization and uralitization. The total Ni + Cu locally reaches up to 4 wt%. Anomalous concentrations of platinum-group elements (PGE's) (maximum 532 ppb Pd, 182 ppb Pt, 53 ppb Rh, 15 ppb Ru, 41 ppb Ir) were detected in samples of Cu-Ni and Ni-Cu ores (maximum 2.63 wt% Ni and 2.31 wt% Cu) from the Jezírka orebody. The main ore paragenesis includes pyrrhotite, pentlandite, chalcopyrite, cubanite, pyrite, magnetite, mackinawite, valleriite, ilmenite and sphalerite. During this work, michenerite, froodite, sperrylite, gold, native bismuth, altaite, tsumoite, hessite, an unnamed Bi-Ni telluride, cobaltite-gersdorffite and galena were newly identified. The host rocks originated through partial melting of a slightly depleted mantle source with noble metals scavenged from this primitive magma prior to the development of these rocks.     

Application of matrix analysis in delineating sinkhole risk areas along highway (I-70 near Frederick,Maryland)

Sinkhole collapse in the area of Maryland Interstate 70 (I-70) and nearby roadways south of Frederick, Maryland, has been posing a threat to the safety of the highway operation as well as other structures. The occurrence of sinkholes is associated with intensive land development. However, the geological conditions that have been developing over the past 200 million years in the Frederick Valley control the locations of the sinkholes. Within an area of approximately 8 km², 138 sinkholes are recorded and their spatial distribution is irregular, but clustered. The clustering indicates the existence of an interaction between the sinkholes. The point pattern of sinkholes is considered to be a sample of a Gibbsian point process from which the hard-core Strauss Model is developed. The radius of influence is calculated for the recorded sinkholes which are most likely to occur within 30 m of an existing sinkhole. The stochastic analysis of the existing sinkholes is biased toward the areas with intensive land use. This bias is adjusted by considering (1) topography, (2) proximity to topographic depressions, (3) interpreted rock formation, (4) soil type, (5) geophysical anomalies, (6) proximity to geologic structures, and (7) thickness of overburden. Based on the properties of each factor, a scoring system is developed and the average relative risk score for individual 30-m segments of the study area is calculated. The areas designated by higher risk levels would have greater risk of a sinkhole collapse than the areas designated by lower risk levels. This risk assessment approach can be updated as more information becomes available.     

Groundwater and volcanoes: examples from the Azores archipelago

Groundwater at the Azores archipelago is a strategic resource for the freshwater supply. Freshwater, mineral and thermal water discharges occur in the archipelago, and especially at the Fogo and Furnas volcanoes (São Miguel). These discharges provide data for case studies of groundwater chemistry from volcanic monitoring due to the stable composition of the sampled waters. The mineral and thermal discharges are mainly of sodium bicarbonate types and present a large range of temperatures, from cold springs to waters at about 90 °C. Some boiling discharges have a sulfate-dominated composition, suggesting a steam-heating mechanism. Geochemical studies on these mineral and thermal waters began in the 19th century. Data gathered since these earlier studies provide a baseline for pH, temperature, CO₂ and major-element composition. Weekly measurements of pH and temperature also denote a rather stable behavior.     

Finite element simulation of seawater intrusion into a quarry-site coastal aquifer,Kocaeli-Darıca,Turkey

A numerical simulation has been carried out to investigate the effects of below-sea-level (b.s.l.) excavation on the raw material quality of a cement quarry in Turkey. The model simulates variations in the hydrodynamic and hydrogeochemical mechanisms in the coastal aquifer upon BSL excavation. In this context, behavior of the seawater intrusion zone, changes in water levels, and salt concentrations have been simulated. In the development of the model, previous geological and hydrogeological reports of the quarry site and the near vicinity have been considered. Eleven new wells (BH-1,...BH-11) have been drilled to reveal hydrogeological features of the area and also for periodical observations of the water levels and hydrogeochemical monitoring. These wells were utilized to develop and calibrate the model to the field conditions. Physical and hydrogeochemical parameters used in the model have been evaluated using available hydrogeological data, the field test results and the related literature. The model has been verified using the field observations. It is based on the virgin conditions of the aquifer as well as on the data for years 1990 and 2001. An average raw material production rate for the cement factory was considered during development of the model, and for making future predictions. Two alternative production scenarios have been considered and probable effects of above-sea-level (a.s.l.) and b.s.l. excavations on seawater intrusion into the aquifer have been studied. Future prediction studies are based on these two production scenarios that assume 43 years of total production (30 years of a.s.l. and 13 years of b.s.l. production) in the quarry. The first scenario, Scenario I, assumes that starting from 2001, the next 30 years would be devoted only to ASL and then the remaining 13 years would be used for b.s.l. production. Scenario II, on the other hand, assumes simultaneous operations both at a.s.l. and b.s.l. levels for the next 43 years after 2001. Effects of b.s.l. production in the quarry site have been simulated accordingly, and seawater intrusion into the aquifer as well as water discharge rates have been predicted for –10, –20 and –30 m production levels.     

Dissolved oxygen and nitrate of groundwater in Choshui Fan-Delta,western Taiwan

Dissolved oxygen (DO) and nitrate concentrations of groundwater from a total of 191 new monitoring wells, with a maximum depth up to 300 m, in the Choshui fan delta of Taiwan were measured to delineate the regional patterns of the two components. Unconfined aquifers of the proximal part of the fan delta contain 1~7 mg/l of DO, and a 250-m deep well still has an unusual value of 1.83 mg/l. DO concentrations decrease downgradient along flow paths to below detection limit (0.01 mg/l) in shallow confined aquifers during the last 40 years, possibly due to oxidation of organic matters in sediments. A plume of high nitrate-N concentrations (0.5–17 mg/l) originating from the proximal part of the fan delta has been developing, possibly since 1956, due to agricultural fertilizers. From 1997 to 2001, the concentration has increased at a rate of ca. 0.28 mg/l per year.     

Improving groundwater-flow modeling using optimal zoning methods

Hydraulic conductivity sometimes exhibits complicated spatial variation over a site. A thorough understanding of the spatial distributions of hydraulic conductivity helps to make deterministic models of groundwater more accurate. This study presents a novel procedure that combines simulated annealing algorithms (SA) and the shortest distance method (SD) with the modular three-dimensional groundwater flow model (MODFLOW). The procedure is applied to a hypothetical site with groundwater-monitoring wells to minimize the difference between simulated and observed hydraulic head for optimal zoning of the spatial distribution of hydraulic conductivity. The results of this optimal zoning method indicate that this new procedure not only improves the efficiency of optimization, but also increases the probability of finding the global optimum, minimizing the errors of the hydraulic head simulated by MODFLOW in two scenarios, one with known and the other with unknown hydraulic conductivity. The results also illustrated that the procedure can effectively determine and delineate hydrogeological zones.