3D numerical modelling of fluid flow in the Val-d’Or orogenic gold district: major crustal shear zones drain fluids from overpressured vein fields

Fluid flow patterns have been determined using oxygen isotope isopleths in the Val-d’Or orogenic gold district. 3D numerical modelling of fluid flow and oxygen isotope exchange in the vein field shows that the fluid flow patterns can be reproduced if the lower boundary of the model is permeable, which represents middle or lower crustal rocks that are infiltrated by a metamorphic fluid generated at deeper levels. This boundary condition implies that the major crustal faults so conspicuous in vein fields do not act as the only major channel for upward fluid flow. The upper model boundary is impermeable except along the trace of major crustal faults where fluids are allowed to drain out of the vein field. This upper impermeable boundary condition represents a low-permeability layer in the crust that separates the overpressured fluid from the overlying hydrostatic fluid pressure regime. We propose that the role of major crustal faults in overpressured vein fields, independent of tectonic setting, is to drain hydrothermal fluids out of the vein field along a breach across an impermeable layer higher in the crust and above the vein field. This breach is crucial to allow flow out of the vein field and accumulation of metals in the fractures, and this breach has major implications for exploration for mineral resources. We propose that tectonic events that cause episodic metamorphic dehydration create a short-lived pulse of metamorphic fluid to rise along zones of transient permeability. This results in a fluid wave that propagates upward carrying metals to the mineralized area. Earthquakes along crustal shear zones cause dilation near jogs that draw fluids and deposit metals in an interconnected network of subsidiary shear zones. Fluid flow is arrested by an impermeable barrier separating the hydrostatic and lithostatic fluid pressure regimes. Fluids flow through the evolving and interconnected network of shear zones and by advection through the rock matrix. Episodic breaches in the impermeable barrier along the crustal shear zones allow fluid flow out of the vein field.     

Alpujarride attribution of the supposed ‘Almagride Complex’ (Betic Internal Zone,Almerı́a Province,Spain)Appartenance Alpujarride du prétendu « Complexe almagride » (Zones internes bétiques,province d'Almeria,Espagne)

The differentiation of units in the Sierra de Almagro has been a source of controversy. There were defined the Almagride and Ballabona–Cucharón complexes, the former considered by several authors as part of a Subbetic metamorphosed and outcropping in a tectonic window. In this study, the units of Ballabona, Almagro and Cucharón are integrated into a single one, that of Tres Pacos, because they correspond to different parts of the same stratigraphic series. This unit is tectonically over the Nevado–Filabride Complex. The existence of the Almagride and Ballabona–Cucharón complexes is discarded and their units form part of the Alpujarride Complex. To cite this article: C. Sanz de Galdeano, F.J. Garc??a Tortosa, C. R. Geoscience 334 (2002) 355–362.     

Computational study of tetrahedral Al–Si and octahedral Al–Mg ordering in phengite

 As part of a wider study of the nature and origins of cation order–disorder in micas, a variety of computational techniques have been used to investigate the nature of tetrahedral and octahedral ordering in phengite, K₂ ^[6](Al₃Mg)^[4](Si₇Al)O₂₀(OH)₄. Values of the atomic exchange interaction parameters J _n used to model the energies of order–disorder were calculated. Both tetrahedral Al–Si and octahedral Al–Mg ordering were studied and hence three types of interaction parameter were necessary: for T–T, O–O and T–O interactions (where T denotes tetrahedral sites and O denotes octahedral sites). Values for the T–T and O–O interactions were taken from results on other systems, whilst we calculated new values for the T–O interactions. We have demonstrated that modelling the octahedral and tetrahedral sheets alone and independently produces different results from modelling a whole T–O–T layer, hence justifying the inclusion of the T–O interactions. Simulations of a whole T–O–T layer of phengite indicated the presence of short-range order, but no long-range order was observed. Received: 8 August 2002 / Accepted: 14 February 2003 Acknowledgements The authors are grateful to EPSRC (EJP) and the Royal Society (CIS) for financial support. Monte Carlo simulations were performed on the Mineral Physics Group's Beowulf cluster and the University of Cambridge's High Performance Computing Facility.     

The impact of hydrological changes on travertine deposits related to thermal springs in the Pamukkale area (SW Turkey)

Pamukkale thermal waters (35 °C), exhibiting calcium-bicarbonate-sulfate composition and high carbon dioxide concentration, are of a predominantly meteoric origin. The meteoric fluid, circulating through faults and fractures, is heated by magmatic intrusions at great depth, and ascends from deep reservoirs to the surface. Mixing with relatively cold groundwater in the near surface zone promotes different saturation conditions with respect to calcium carbonate that later precipitates at depth and/or the surface. Dissolution-deposition processes of calcium carbonate both at surface and depth environments may help to reconstruct past climate direction in the field. During wet climate conditions a high-rate of calcium carbonate accumulation would be expected to occur at the surface because thermal fluid would be under-saturated with respect to calcium carbonate at depth because of a relatively higher mixing ratio with cold groundwater. During dry climate conditions the thermal fluid would be super-saturated at depth because of the highly acidic environment. Hydrometeorological studies reveal that the annual precipitation at the Pamukkale hydrothermal field tends to decrease with time. This climatic change in the area was also detected from geological records. While humid climate conditions prevailed during the late Quaternary, the area has recently been affected by arid/semi-arid climate conditions, followed by some episodic transitions. This study has shown how the system has possibly reacted to different climate conditions since antiquity.     

Diagenetic paths in the margin of a Triassic Basin: NW zone of the Iberian Chain, Spain

Buntsandstein deposits generated in a slowly subsiding basin on the western margin of the Iberian Chain are represented by a stratigraphic succession of fluvial deposits less than 100 m thick (conglomerates, sandstones, and shales). Diagenetic processes in sandstones can be grouped as eodiagenetic, mesodiagenetic, and telodiagenetic. Eodiagenesis can be associated with Muschelkalk, Keuper, and probably early Jurassic times. Mesodiagenesis is probably related to Jurassic times. Diagenetic chemical reactions suggest a maximum burial less than 1.5 km and low temperatures (<120°C). Patterns of porosity reduction by compaction and cementation suggest four diagenetic stages: (1) Loss of primary porosity by early mechanical compaction; (2) early cementation (K-feldspar and dolomite); (3) dissolution of cements; and (4) framework collapse by re-compaction. These stages are manifested by the presence of two types of sandstone. Type I sandstones present high intergranular volume (mean, 30%). Type II sandstones are characterized by high compactional porosity loss and exhibit low values of intergranular volume (mean, 16.9%). Type II sandstones are associated with the dissolution of cement and later re-compaction of type I sandstones. An intermediate telodiagenetic phase is deduced and related to the sharp unconformity between Lower Cretaceous sediments and the underlying sediments. This suggests that a mechanically unstable framework collapsed during the Cretaceous, generating type II sandstones. The analyzed diagenetic paths have a wide applicability on similar marginal areas of rift basins.     

Emplacement age and thermal footprint of the diamondiferous Ellendale E9 lamproite pipe, Western Australia

The diamondiferous Ellendale 9 (E9) pipe is a funnel-shaped maar-diatreme volcano consisting of inward-dipping tuff sequences intruded by lamproite plugs and dykes. The host rocks for the E9 pipe are Permian sandstones. The multiple lithological contacts exposed within the mined maar volcano provide a natural laboratory in which to study the effect of volcanic processes on U–Th–Pb–He systematics. Zircon from the regional sandstone and E9 lamproite display a bimodal distribution of ages on (U–Th)/He–U/Pb plots. The zircon U/Pb ages for the E9 pipe (n?=?52) range from 440 to 2,725 Ma, while the cluster of (U–Th)/He ages for the lamproite dyke zircon indicate that dyke emplacement occurred at 20.6?±?2.8 Ma, concordant with a maximum emplacement age of about ≤22 Ma from phlogopite ⁴⁰Ar/³⁹Ar. These ages indicate a xenocrystic origin for the zircon entrained in the E9 dyke. The U/Pb ages of detrital zircon from the regional sandstone host (373–3,248 Ma; n?=?41) are indistinguishable from those of the lamproite zircon xenocrysts, whereas the detrital zircon in the host sandstone yield (U–Th)/He ages from 260 to 1,500 Ma. A thermochronology traverse across the E9 lamproite dyke reveals that the zircon (U–Th)/He ages in the host sandstone have not been significantly thermally reset during dyke emplacement, even at the contact. The capability of the zircon (U–Th)/He method to distinguish deep, mantle source lithologies from upper crustal source lithologies could be used in geochemical exploration for diamonds. Pre-screening of detrital samples using etching and helium assay methods will improve the efficiency and decrease the cost of greenfields exploration.     

Some Insights into Topographic,Elastic and Self-gravitation Interaction in Modelling Ground Deformation and Gravity Changes in Active Volcanic Areas

Surface displacements and gravity changes due to volcanic sources are influenced by medium properties. We investigate topographic, elastic and self-gravitation interaction in order to outline the major factors that are significant in data modelling. While elastic-gravitational models can provide a suitable approximation to problems of volcanic loading in areas where topographic relief is negligible, for prominent volcanoes the rough topography could affect deformation and gravity changes to a greater extent than self-gravitation. This fact requires the selection, depending on local relief, of a suitable model for use in the interpretation of surface precursors of volcanic activity. We use the three-dimensional Indirect Boundary Element Method to examine the effects of topography on deformation and gravity changes in models of magma chamber inflation/deflation. Topography has a significant effect on predicted surface deformation and gravity changes. Both the magnitude and pattern of the geodetic signals are significantly different compared to half-space solutions. Thus, failure to account for topographic effects in areas of prominent relief can bias the estimate of volcanic source parameters, since the magnitude and pattern of deformation and gravity changes depend on such effects.     

North Atlantic atmospheric circulation and surface wind in the Northeast of the Iberian Peninsula: uncertainty and long term downscaled variability

The variability and predictability of the surface wind field at the regional scale is explored over a complex terrain region in the northeastern Iberian Peninsula by means of a downscaling technique based on Canonical Correlation Analysis. More than a decade of observations (1992–2005) allows for calibrating and validating a statistical method that elicits the main associations between the large scale atmospheric circulation over the North Atlantic and Mediterranean areas and the regional wind field. In an initial step the downscaling model is designed by selecting parameter values from practise. To a large extent, the variability of the wind at monthly timescales is found to be governed by the large scale circulation modulated by the particular orographic features of the area. The sensitivity of the downscaling methodology to the selection of the model parameter values is explored, in a second step, by performing a systematic sampling of the parameters space, avoiding a heuristic selection. This provides a metric for the uncertainty associated with the various possible model configurations. The uncertainties associated with the model configuration are considerably dependent on the spatial variability of the wind. While the sampling of the parameters space in the model set up moderately impact estimations during the calibration period, the regional wind variability is very sensitive to the parameters selection at longer timescales. This fact illustrates that downscaling exercises based on a single configuration of parameters should be interpreted with extreme caution. The downscaling model is used to extend the estimations several centuries to the past using long datasets of sea level pressure, thereby illustrating the large temporal variability of the regional wind field from interannual to multicentennial timescales. The analysis does not evidence long term trends throughout the twentieth century, however anomalous episodes of high/low wind speeds are identified.