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.     

Celestite formation, bacterial sulphate reduction and carbonate cementation of Eocene reefs and basinal sediments (Igualada, NE Spain)

Petrographic and geochemical studies of an Upper Eocene reef and associated basinal sediments from the mixed carbonate–siliciclastic fill of the south‐eastern Pyrenean foreland basin near Igualada (NE Spain) provide new insights into the evolution of subsurface hydrology during the restriction of a marine basin. The reef deposits are located on delta‐lobe sandstones and prodelta marls, which are overlain by hypersaline carbonates and Upper Eocene evaporites. Authigenic celestite (SrSO₄) is an important component in the observed diagenetic sequences. Celestite is a significant palaeohydrological indicator because its low solubility constrains transportation of Sr²⁺ and SO₄^2? in the same diagenetic fluid. Stable isotopic analyses of carbonates in the reef indicate that meteoric recharge was responsible for aragonite stabilization and calcite cementation. Sulphur and oxygen isotope geochemistry of the celestite demonstrates that it formed from residual sulphate after bacterial sulphate reduction, but also requires that there was a prior episode of sulphate recycling. Meteoric water reaching the reef and basinal areas was most probably charged with SO₄^2? from the dissolution of younger Upper Eocene marine evaporites. This sulphate, combined with organic matter present in the sediments, fuelled bacterial sulphate reduction in the meteoric palaeoaquifer. Strontium for celestite precipitation was partly derived in situ from dissolution of aragonite corals in the reef and basinal counterparts. However, ⁸⁷Sr/⁸⁶Sr data also suggest that Sr²⁺ was partly derived from dissolution of overlying evaporites. Mixing of these two fluids promoted celestite formation. The carbonate stable isotopic data suggest that the local meteoric water was enriched in ¹⁸O compared with that responsible for stabilization of other reefs along the basin margin. Furthermore, meteoric recharge at Igualada post‐dated evaporite deposition in the basin, whereas other parts of the same reef complex were stabilized before evaporite formation. This discrepancy resulted from the spatial distribution of continental siliciclastic units that acted as groundwater conduits.     

Induced Microearthquake Patterns in Hydrocarbon and Geothermal Reservoirs: Six Case Studies

—?The injection or production of fluids can induce microseismic events in hydrocarbon and geothermal reservoirs. By deploying sensors downhole, data sets have been collected that consist of a few hundred to well over 10,000 induced events. We find that most induced events cluster into well-defined geometrical patterns. In many cases, we must apply high-precision, relative location techniques to observe these patterns. At three sedimentary sites, thin horizontal strands of activity are commonly found within the location patterns. We believe this reflects fracture containment between stratigraphic layers of differing mechanical properties or states of stress. At a massive carbonate and two crystalline sites, combinations of linear and planar features indicate networks of intersecting fractures and allow us to infer positions of aseismic fractures through their influence on the location patterns. In addition, the fine-scale seismicity patterns often evolve systematically with time. At sedimentary sites, migration of seismicity toward the injection point has been observed and may result from slip-induced stress along fractures that initially have little resolved shear. In such cases, triggering events may be critical to generate high levels of seismic activity. At one crystalline site, the early occurrence of linear features that traverse planes of activity indicate permeable zones and possible flow paths within fractures. We hope the continued development of microseismic techniques and refinement of conceptual models will further increase our understanding of fluid behavior and lead to improved resource management in fractured reservoirs.     

A mineral equilibria study of the hydrothermal alteration in mafic greenschist facies rocks at Kalgoorlie, Western Australia

The influx of a H₂O–CO₂‐dominated fluid into actinolite‐bearing metabasic rocks during greenschist facies metamorphism in the Kalgoorlie area of Western Australia resulted in a zoned alteration halo around inferred fluid conduits that contain gold mineralisation. The alteration halo is divided into two outer zones, the chlorite zone and the carbonate zone, and an inner pyrite zone adjacent to the inferred fluid conduits. Reaction between the fluid and the protolith resulted in the breakdown of actinolite and the development of chlorite, dolomite, calcite and siderite. In addition, rocks in the pyrite zone developed muscovite‐bearing assemblages as a consequence of the introduction of potassium by the fluid. Mineral equilibria calculations undertaken using the computer software thermocalc in the model system Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O–CO₂ show that mineral assemblages in the outer zones of the alteration halo are consistent with equilibrium of the protoliths with a fluid of composition X_CO2 = CO₂/(CO₂ + H₂O) = 0.1–0.25 for temperatures of 315–320 °C. The inner zone of the alteration halo reflect equilibrium with a fluid of composition X_CO2≈ 0.25. Fluid‐rock buffering calculations show that the alteration halo is consistent with interaction with a single fluid composition and that the zoned structure of the halo reflects the volume of this fluid with which the rocks reacted. This fluid is likely to have also been the one responsible for the gold mineralisation at Kalgoorlie.     

Toxic Cyanobacteria in Greek Freshwaters, 1987—2000: Occurrence,Toxicity, and Impacts in the Mediterranean Region

In a survey in Greece from 1987 to 2000 hepatotoxic cyanobacterial blooms were observed in 9 out of 33 freshwaters. Microcystins (MCYSTs) were detected by HPLC in 7 of these lakes, and the total MCYST concentration per scum dry weight ranged from 50.3 to 1638 ± 464 μg g^—1. Cyanobacterial genera (Microcystis, Anabaena, Anabaenopsis, Aphanizomenon, Cylindrospermopsis) with known toxin producing taxa were present in 31 freshwaters. From our data and a review of the literature, it would appear that Mediterranean countries are more likely 1) to have toxic cyanobacterial blooms consisting of Microcystis spp. and 2) to have higher intracellular MCYST concentrations. A case study in Lake Kastoria is used to highlight seasonal patterns of cyanobacterial and MCYST‐LR occurrence and to assess cyanotoxin risk. Cyanobacterial biovolume was high (> 11 μL L^—1) throughout the year and was in excess of Guidance Level 2 (10 μL L^—1) proposed by WHO for recreational waters and Alert Level 2 for drinking water. Further, surface water samples from April to November exceeded Guidance Level 3, with the potential for acute cyanobacterial poisoning. Intracellular MCYST‐LR concentrations (max 3186 μg L^—1) exceeded the WHO guideline for drinking water (1 μg L^—1) from September to November with a high risk of adverse health effects. Preliminary evidence indicates that in 3 lakes microcystins are accumulated in some aquatic organisms. Generally, a high risk level can be deduced from the data for the Mediterranean region.     

Modeling Spatially and Temporally Complex Land‐Cover Change: The Case of Western Honduras*

This article presents an econometric analysis of land‐cover change in western Honduras. Ground‐truthed satellite image analysis indicates that between 1987 and 1996 net reforestation occurred in the 1,015‐km² study region. While some reforestation can be attributed to a 1987 ban on logging, the area of reforestation greatly exceeds that of previously clear‐cut areas. Further, new area was also deforested between 1987 and 1996. Thus, the observed land‐cover changes represent a complex mosaic of changing land‐use patterns across time and space. The analysis contributes to the literature on land‐cover change modeling in that: (1) it compares two econometric approaches to capture complex and often bidirectional changes in land cover from 1987 to 1996 as a function of agricultural suitability and transportation costs, and (2) it addresses techniques to identify and correct for spatial autocorrelation in a categorical regression framework.     

The Elusive AD 1826 Tsunami, South Westland, New Zealand

In AD 1826 sealers reported earthquake and tsunami activity in Fiordland, although, contemporary or near‐contemporary accounts of tsunami inundation at the time are elusive. A detailed analysis of recent sediments from Okarito Lagoon builds on contextual evidence provided by earlier research concerning past tsunami inundation. Sedimentological, geochemical, micropalaeontological and geochronological data are used to determine palaeoenvironments before, during and after what was most probably tsunami inundation in AD 1826. The most compelling chronological control is provided fry a young cohort of trees growing on a raised shoreline bench stranded fry a drop in the lagoon water level following tsunami inundation.     

Channel adjustments of the lower Trinity River,Texas, downstream of Livingston Dam

Channel cross‐sectional changes since construction of Livingston Dam and Lake Livingston in 1968 were studied in the lower Trinity River, Texas, to test theoretical models of channel adjustment, and to determine controls on the spatial extent of channel response. High and average flows were not significantly modified by the dam, but sediment transport is greatly reduced. The study is treated as an opportunistic experiment to examine the effects of a reduction in sediment supply when discharge regime is unchanged. Channel scour is evident for about 60 km downstream, and the general phenomena of incision, widening, coarsening of channel sediment and a decrease in channel slope are successfully predicted, in a qualitative sense, by standard models of channel response. However, there is no consistent channel response within this reach, as various qualitatively different combinations of increases, decreases or no change in width, depth, slope and roughness occur. These multiple modes of adjustment are predicted by the unstable hydraulic geometry model. Between about 60 km and the Trinity delta 175 km downstream of the dam, no morphological response to the dam is observed. Rather than a diminution of the dam's effects on fluvial processes, this is due to a fundamental change in controls of the fluvial system. The downstream end of the scour zone corresponds to the upstream extent of channel response to Holocene sea level rise. Beyond 60 km downstream, the Trinity River is characterized by extensive sediment storage and reduced conveyance capacity, so that even after dam construction sediment supply still exceeds transport capacity. The channel bed of much of this reach is near or below sea level, so that sea level rise and backwater effects from the estuary are more important controls on the fluvial system than upstream inputs. Copyright © 2005 John Wiley & Sons, Ltd.     

Entropy Analysis of Multiple Scale Causality and Qualitative Causal Shifts in Spatial Systems*

Spatial systems are typically characterized by multiple controlling factors and processes operating at different spatial and temporal scales (multiple scale causality [MSC]). An entropy decomposition‐based approach to MSC is presented here in two contexts. First, given maps or distributions of an observed phenomenon at two or more scales, the contribution at more local or global (relative to the primary scale of observation) controls to the observed entropy can be estimated. Second, a theoretical treatment of the entropy decomposition equations shows that as the range of scale is increased by broadening or narrowing resolutions or by incorporating more controls, the influence of larger or smaller‐scale influences not only changes, but may change qualitatively, e.g., in terms of having positive (entropy‐increasing) or negative (information‐increasing) effects. Such qualitative causal shifts have implications for efforts to use any single causal explanation across the molecular to planetary spatial and instantaneous to geological range of scales relevant to physical geography. The entropy decomposition method is illustrated with an application to soil landscapes in the Ouachita Mountains, Arkansas.