Mössbauer, X-ray and Magnetic Studies of Black Sand from the Italian Mediterranean Sea

The study of natural magnetic sands is instrumental to investigate the geological aspects of their formation and of the origin of their territory. In particular, Mössbauer spectroscopy provides unique information on their iron content and on the oxidation state of iron in their mineral composition. The Italian coast on the Mediterranean Sea near Rome is known for the presence of highly magnetic black sands of volcanic origin. A study of the room temperature Mössbauer spec- trum, powder X-ray diffraction, energy dispersive X-ray spectroscopy, and magnetic measurements of a sample of black magnetic sand collected on the seashore of the town of Ladispoli is performed. This study reveals magnetite as main constituent with iron in both tetrahedral and octahedral sites. Minor constituents are the iron minerals hematite and ilmenite, the iron containing minerals diopsite, gossular, and allanite, as well as ubiquitous sanidine, quartz, and calcite.

Magmatic–hydrothermal molybdenum isotope fractionation and its relevance to the igneous crustal signature

We analysed the Mo isotope composition of a comprehensive series of molybdenite samples from the porphyry- type Questa deposit (NM, USA), as well as one rhyolite and one granite sample, directly associated with the Mo mineralization. The δ98Mo of the molybdenites ranges between −0.48‰ and +0.40‰, with a median at −0.05‰. The median Mo isotope composition increases from early magmatic (−0.29‰) to hydrothermal (−0.05‰) breccia mineralization (median bulk breccia = −0.17‰) to late stockwork veining (+0.22‰). Moreover, variations of up to 0.34‰ are found between different molybdenite crystals within an individual hand specimen. The rhyolite sample with 0.12 μg g−1 Mo has δ98Mo = −0.57‰ and is lighter than all molybde- nites from the Questa deposit, interpreted to represent the igneous leftover after aqueous ore fluid exsolution. We recognize three Mo isotope fractionation processes that occur between about 700 and 350 °C, affecting the Mo iso- tope composition of magmatic–hydrothermal molybdenites. Δ1Mo: Minerals preferentially incorporate light Mo isotopes during progressive fractional crystallization in subvolcanic magma reservoirs, leaving behind a melt enriched in heavy Mo isotopes. Δ2Mo: Magmatic–hydrothermal fluids preferentially incorporate heavy Mo iso- topes upon fluid exsolution. Δ3Mo: Light Mo isotopes get preferentially incorporated in molybdenite during crys- tallization from an aqueous fluid, leaving behind a hydrothermal fluid that gets heavier with progressive molybdenite crystallization. The sum of all three fractionation processes produces molybdenites that record heavier δ98Mo compositions than their source magmas. This implies that the mean δ98Mo of molybdenites published so far (~0.4‰) likely represents a maximum value for the Mo isotope composition of Phanerozoic igneous upper crust.

Mo isotopic composition of the mid-Neoproterozoic ocean: an iron formation perspective

The Neoproterozoic was a major turning point in Earth's surficial history, recording several widespread glaciations, the first appearance of complex metazoan life, and a major increase in atmospheric oxygen. Marine redox proxies have resulted in many different estimates of both the timing and magnitude of the increase in free oxygen, although the consensus has been that it occurred following the Marinoan glaciation, the second globally recorded “snowball Earth” event. A critically understudied rock type of the Neoproterozoic is iron formation associated with the Sturtian (first) glaciation. Samples from the <716 Ma Rapitan iron formation were analysed for their Re concentrations and Mo isotopic composition to refine the redox history of its depositional basin. Rhenium concentrations and Re/Mo ratios are consistently low throughout the bottom and middle of the iron formation, reflecting ferruginous to oxic basinal conditions, but samples from the uppermost jasper layers of the iron formation show significantly higher Re concentrations and Re/Mo ratios, indicating that iron formation deposition was terminated by a shift towards a sulfidic water column. Similarly, the δ98Mo values are close to 0.0‰ throughout most of the iron formation, but rise to ~+0.7‰ near the top of the section. The δ98Mo from samples of ferruginous to oxic basinal conditions are the product of adsorption to hematite, indicating that the Neoproterozoic open ocean may have had a δ98Mo of ~1.8‰. Together with the now well-established lack of a positive Eu anomaly in Neoproterozoic iron formations, these results suggest that the ocean was predominantly oxygenated at 700 Ma.

Prediction of water–rock interaction and porosity evolution in a granitoid-hosted enhanced geothermal system, using constraints from the 5 km Basel-1 well

Numerical simulations based on plans for a deep geothermal system in Basel, Switzerland are used here to understand chemical processes that occur in an initially dry granitoid reservoir during hydraulic stimulation and long-term water circulation to extract heat. An important question regarding the sustainability of such enhanced geothermal systems (EGS), is whether water–rock reactions will eventually lead to clogging of flow paths in the reservoir and thereby reduce or even completely block fluid throughput. A reactive transport model allows the main chemical reactions to be predicted and the resulting evolution of porosity to be tracked over the expected 30-year operational lifetime of the system. The simulations show that injection of surface water to stimulate fracture permeability in the monzogranite reservoir at 190 °C and 5000 m depth induces redox reactions between the oxidised surface water and the reduced wall rock. Although new calcite, chlorite, hematite and other minerals precipitate near the injection well, their volumes are low and more than compensated by those of the dissolving wall-rock minerals. Thus, during stimulation, reduction of injectivity by mineral precipitation is unlikely. During the simulated long-term operation of the system, the main mineral reactions are the hydration and albitization of plagioclase, the alteration of hornblende to an assemblage of smectites and chlorites and of primary K-feldspar to muscovite and microcline. Within a closed-system doublet, the composition of the circulated fluid changes only slightly during its repeated passage through the reservoir, as the wall rock essentially undergoes isochemical recrystallization. Even after 30 years of circulation, the calculations show that porosity is reduced by only ∼0.2%, well below the expected fracture porosity induced by stimulation. This result suggests that permeability reduction owing to water–rock interaction is unlikely to jeopardize the long-term operation of deep, granitoid-hosted EGS systems. A peculiarity at Basel is the presence of anhydrite as fracture coatings at ∼5000 m depth. Simulated exposure of the circulating fluid to anhydrite induces a stronger redox disequilibrium in the reservoir, driving dissolution of ferrous minerals and precipitation of ferric smectites, hematite and pyrite. However, even in this scenario the porosity reduction is at most 0.5%, a value which is unproblematic for sustainable fluid circulation through the reservoir.

Reactive mass transport modelling of a three-dimensional vertical fault zone with a finger-like convective flow regime

For a three-dimensional vertically-oriented fault zone, we consider the coupled effects of fluid flow, heat transfer and reactive mass transport, to investigate the patterns of fluid flow, temperature distribution, mineral alteration and chemically induced porosity changes. We show, analytically and numerically, that finger-like convection patterns can arise in a vertically-oriented fault zone. The onset and patterns of convective fluid flow are controlled by the Rayleigh number which is a function of the thermal properties of the fluid and the rock, the vertical temperature gradient, and the height and the permeability of the fault zone. Vigorous fluid flow causes low temperature gradients over a large region of the fault zone. In such a case, flow across lithological interfaces becomes the most important mechanism for the formation of sharp chemical reaction fronts. The degree of rock buffering, the extent and intensity of alteration, the alteration mineralogy and in some cases the formation of ore deposits are controlled by the magnitude of the flow velocity across these compositional interfaces in the rock. This indicates that alteration patterns along compositional boundaries in the rock may provide some insights into the convection pattern. The advective mass and heat exchanges between the fault zone and the wallrock depend on the permeability contrast between the fault zone and the wallrock. A high permeability contrast promotes focussed convective flow within the fault zone and diffusive exchange of heat and chemical reactants between the fault zone and the wallrock. However, a more gradual permeability change may lead to a regional-scale convective flow system where the flow pattern in the fault affects large-scale fluid flow, mass transport and chemical alteration in the wallrocks

The Yanque prospect (Peru): from polymetallic Zn-Pb mineralization to a nonsulfide deposit

The Yanque nonsulfide Pb-Zn deposit (inferred resources 12.5 Mt @ 3.7% Pb and @ 3.5% Zn) is located in the Andahuaylas-Yauri ore province (Cuzco, southern Peru). The deposit occurs within a base metal mineralized district, centered on the medium-sized Dolores porphyry copper. A thorough geological, mineralogical and geochemical study has carried out in order to define: the relationships between the Dolores Cu-porphyry ore and the Yanque Zn-Pb polymetallic mineralization, and the characteristics of the economic nonsulfide concentrations. Both sedimentary and igneous rocks constitute the backbone of the Yanque-Dolores area. The sedimentary lithologies belong to the Soraya, Mara and Ferrobamba Fms. (upper Jurassic-middle Cretaceous). The Yanque orebody is hosted by the Mara Fm., which prevailingly consists of a siliciclastic sedimentary breccia. The original sulfide mineralization consisted of galena, pyrite and sphalerite. The host rock has been affected by a strong hydrothermal alteration, characterized by prevailing sericite/illite, as in the typical porphyry-related phyllic-argillic alteration stage, and by minor kaolinite, dolomite and quartz. Minor element geochemistry, characterized by Sb, As, Mn, Ag and locally also by Cu, points to magmatic-hydrothermal related mineralizing fluids. The Pb isotopic compositions from Dolores and Yanque sulfides are similar, and are typical of the Tertiary magmatically-derived ores in this part of Peru. The hydrothermally altered rocks at Yanque have the same Pb isotopic compositions as the sulfides, thus confirming the hypothesis that the Yanque primary Zn-Pb mineralization may have been produced by hydrothermal circulation related to the emplacement of the Dolores Cu-porphyry, as it is the case of other porphyry Cu systems associated with polymetallic mineralization elsewhere. However, no simple genetic model for the mineralization involving just one fluid circulation episode is able to explain the data. The Yanque economic nonsulfide ore association consists of sauconite, hemimorphite, smithsonite and cerussite, which result from the weathering and alteration of the original sulfide mineralization. Zinc is allocated mainly in sauconite (Zn-smectite), rather than in carbonates: a factor strictly related to the prevailing siliciclastic character of the host rock. Distinctive features of the Yanque orebody are the comparable ore grades for both Pb and Zn (3.5% Zn and 3.7% Pb), and the inverse supergene chemical zoning. In fact, contrary to other supergene ores of this type, zinc prevails in the top zone of the Yanque deposit, whereas lead content increases with depth. Considering the different mobility of the two metals in solution, it may be assumed that most of the primary zinc that was the source for the Yanque nonsulfides was originally located far from the position occupied by the galena mineralization, whose remnants have been observed on site. Zinc sulfides may have been originally contained in the eroded rock volumes that surrounded the actual deposit: the zinc-rich solutions have possibly migrated through the siliciclastic Mara Fm. and precipitated the nonsulfide minerals by porosity filling and replacement processes. In this sense, the Yanque secondary Zn-Pb deposit could be considered as a special type of “Exotic” mineralization.

Elemental and Lead Isotopic Data of Copper Finds from the Singen Cemetery, Germany – a Methodological Approach to Investigate Early Bronze Age Trade Networks

We report a trace element - Pb isotope analytical (LIA) database on the "Singen Copper", a peculiar type of copper found in the North Alpine realm, from its type locality, the Early Bronze Age Singen Cemetery (Germany). What distinguishes “Singen Copper” from other coeval copper types? (i) is it a discrete metal lot with a uniform provenance (if so, can its provenance be constrained)? (ii) was it manufactured by a special, unique metallurgical process that can be discriminated from others? Trace element concentrations can give clues on the ore types that were mined, but they can be modified (more or less intentionally) by metallurgical operations. A more robust indicator are the ratios of chemically similar elements (e.g. Co/Ni, Bi/Sb, etc.), since they should remain nearly constant during metallurgical operations, and are expected to behave homogeneously in each mineral of a given mining area, but their partition amongst the different mineral species is known to cause strong inter-element fractionations. We tested the trace element ratio pattern predicted by geochemical arguments on the Brixlegg mining area. Brixlegg itself is not compatible with the Singen Copper objects, and we only report it because it is a rare instance of a mining area for which sufficient trace element analyses are available in the literature. We observe that As/Sb in fahlerz varies by a factor 1.8 above/below median; As/Sb in enargite varies by a factor of 2.5 with a 10 times higher median. Most of the 102 analyzed metal objects from Singen are Sb-Ni-rich, corresponding to “antimony-nickel copper” of the literature. Other trace element concentrations vary by > 100 times, ratios by factors > 50. Pb isotopic compositions are all significantly different from each other. They do not form a single linear array and require > 3 ore batches that certainly do not derive from one single mining area. Our data suggest a heterogeneous provenance of “Singen copper”. Archaeological information limits the scope to Central European sources. LIA requires a diverse supply network from many mining localities, including possibly Brittany. Trace element ratios show more heterogeneity than LIA; this can be explained either by deliberate selection of one particular ore mineral (from very many sources) or by processing of assorted ore minerals from a smaller number of sources, with the unintentional effect that the quality of the copper would not be constant, as the metallurgical properties of alloys would vary with trace element concentrations.

Predictive modelling of mineral scaling, corrosion and the performance of solute geothermometers in a granitoid-hosted, enhanced geothermal system

Experience is lacking with mineral scaling and corrosion in enhanced geothermal systems (EGS) in which surface water is circulated through hydraulically stimulated crystalline rocks. As an aid in designing EGS projects we have conducted multicomponent reactive-transport simulations to predict the likely characteristics of scales and corrosion that may form when exploiting heat from granitoid reservoir rocks at ∼200 °C and 5 km depth. The specifications of an EGS project at Basel, Switzerland, are used to constrain the model. The main water–rock reactions in the reservoir during hydraulic stimulation and the subsequent doublet operation were identified in a separate paper (Alt-Epping et al., 2013b). Here we use the computed composition of the reservoir fluid to (1) predict mineral scaling in the injection and production wells, (2) evaluate methods of chemical geothermometry and (3) identify geochemical indicators of incipient corrosion. The envisaged heat extraction scheme ensures that even if the reservoir fluid is in equilibrium with quartz, cooling of the fluid will not induce saturation with respect to amorphous silica, thus eliminating the risk of silica scaling. However, the ascending fluid attains saturation with respect to crystalline aluminosilicates such as albite, microcline and chlorite, and possibly with respect to amorphous aluminosilicates. If no silica-bearing minerals precipitate upon ascent, reservoir temperatures can be predicted by classical formulations of silica geothermometry. In contrast, Na/K concentration ratios in the production fluid reflect steady-state conditions in the reservoir rather than albite–microcline equilibrium. Thus, even though igneous orthoclase is abundant in the reservoir and albite precipitates as a secondary phase, Na/K geothermometers fail to yield accurate temperatures. Anhydrite, which is present in fractures in the Basel reservoir, is predicted to dissolve during operation. This may lead to precipitation of pyrite and, at high exposure of anhydrite to the circulating fluid, of hematite scaling in the geothermal installation. In general, incipient corrosion of the casing can be detected at the production wellhead through an increase in H2(aq) and the enhanced precipitation of Fe-bearing aluminosilicates. The appearance of magnetite in scales indicates high corrosion rates.