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.     

The evolution of groundwater in the Tyrrell catchment, south-central Murray Basin, Victoria, Australia

The Tyrell catchment lies on the western margin of the Riverine Province in the south-central Murray Basin, one of Australia’s most important groundwater resources. Groundwater from the shallow, unconfined Pliocene Sands aquifer and the underlying Renmark Group aquifer is saline (total dissolved solids up to 150,000 mg/L) and is Na-Cl-Mg type. There is no systematic change in salinity along hydraulic gradients implying that the aquifers are hydraulically connected and mixing during vertical flow is important. Stable isotopes (¹⁸O+²H) and Cl/Br ratios indicate that groundwater is entirely of meteoric origin and salts in this system have largely been derived by evapotranspiration of rainfall with only minor halite dissolution, rock weathering (mainly feldspar dissolution), and ion exchange between Na and Mg on clays. Similarity in chemistry of all groundwater in the catchment implies relative consistency in processes over time, independent of any climatic variation. Groundwater in both the Pliocene Sands and Renmark Group aquifers yield ages of up to 25 ka. The Tyrrell Catchment is arid to semi-arid and has low topography. This has resulted in relatively low recharge rates and hydraulic gradients that have resulted in long groundwater residence times.     

Quantification of karst aquifer discharge components during storm events through end-member mixing analysis using natural chemistry and stable isotopes as tracers

Karst aquifer components that contribute to the discharge of a water supply well in the Classical Karst (Kras) region (Italy/Slovenia) were quantitatively estimated during storm events. Results show that water released from storage within the epikarst may comprise as much as two-thirds of conduit flow in a karst aquifer following rainfall. Principal components analysis (PCA) and end-member mixing analysis (EMMA) were performed using major ion chemistry and the stable isotopes of water (δ¹⁸O, δ²H) and of dissolved inorganic carbon (δ¹³C_DIC) to estimate mixing proportions among three sources: (1) allogenic river recharge, (2) autogenic recharge, and (3) an anthropogenic component stored within the epikarst. The sinking river most influences the chemical composition of the water-supply well under low-flow conditions; however, this proportion changes rapidly during recharge events. Autogenic recharge water, released from shallow storage in the epikarst, displaces the river water and is observed at the well within hours after the onset of precipitation. The autogenic recharge end member is the second largest component of the well chemistry, and its contribution increases with higher flow. An anthropogenic component derived from epikarstic storage also impacts the well under conditions of elevated hydraulic head, accounting for the majority of the chemical response at the well during the wettest conditions.     

Hydrochemistry and isotope geochemistry as tools for groundwater hydrodynamic investigation in multilayer aquifers: a case study from Lomellina, Po plain, South-Western Lombardy, Italy

A multicriteria approach in studying hydrodynamics of a multilayer aquifer system has been used in the Lomellina region (Northern Italy). It involves the reconstruction of the hydrogeological framework coupled to the definition of the hydrochemical and isotopic features of the aquifers. A shallow phreatic aquifer, reaching depths of about 60–80 m from the surface, and deeper aquifers containing confined groundwater, were distinguished. Groundwater generally shows mineralisation decreasing with depth; dissolved ions depict calcium-bicarbonate hydrochemical facies and stable isotopes define the recharge mechanisms, the origin of groundwater, and the hydraulic confinement of deep aquifers. The phreatic aquifer is fed by local infiltration and by streams and irrigation channels. Tritium and Carbon-14 groundwater dating indicate long residence times (on the order of thousands of years) for confined aquifers. The confined aquifers show essentially passive hydrodynamic conditions and maintain a higher piezometric level than the phreatic aquifer. This inhibits the possibility of recent water penetrating far below the surface. The hydrogeological setting of the Lomellina region displays features which are common to other sectors of the Po plain. As a consequence, the results of this study, although conducted on a restricted area, are highly illustrative of groundwater hydrodynamics in large sedimentary aquifers.     

Observations Of Nocturnal Drainage Flow In A Shallow Gully

The formation of cold air drainage flows in a shallow gully is studied during CASES-99 (Cooperative Atmosphere-Surface Exchange Study). Fast and slow response wind and temperature measurements were obtained on an instrumented 10-m tower located in the gully and from a network of thermistors and two-dimensional sonic anemometers, situated across the gully. Gully flow formed on clear nights even with significant synoptic flow. Large variations in surface temperature developed within an hour after sunset and in situ cooling was the dominant factor in wind sheltered locations. The depth of the drainage flow and the height of the down-gully wind speed maximum were found to be largest when the external wind speed above the gully flow is less than 2 m s^-1. The shallow drainage current is restricted to a depth of a few metres, and is deepest when the stratification is stronger and the external flow is weaker. During the night the drainage flow breaks down, sometimes on several occasions, due to intermittent turbulence and downward fluxes of heat and momentum. The near surface temperature may increase by 6 ° C in less than 30 min due to the vertical convergence of downward heat flux. The mixing events are related to acceleration of the flow above the gully flow and decreased Richardson number. These warming events also lead to warming of the near surface soil and reduction of the upward soil heat flux. To examine the relative importance of different physical mechanisms that could contribute to the rapid warming, and to characterize the turbulence generated during the intermittent turbulent periods, the sensible heat budget is analyzed and the behaviour of different turbulent parameters is discussed.     

Stable isotope values of Costa Rican surface waters

Stable isotope data of surface waters from the humid tropics in general, and Costa Rica in particular, are scarce. To improve our understanding of the spatial distribution of stable isotopes in surface waters, we measured δ¹⁸O and δD in river and lake (n=63) and precipitation (n=3) samples from Costa Rica. We also present data from the IAEA/WMO isotopes in precipitation network as context for our study. Surface water isotope values do not strongly correlate with elevation, stream head elevation, stream length, distance from Caribbean Sea, or estimated mean annual precipitation for the country as a whole. However, the data show distinct regional trends. The δ¹⁸O and δD values downwind of mountain ranges are inversely related to the altitude of the ranges the air masses traverse. In the lee of the high Talamanca Range, δ¹⁸O values are 6–8‰ lower, while in the lee of the lower Tilarán Range δ¹⁸O values are 2–3‰ lower than upwind sites along the Caribbean Slope. An altitude effect of −1.4‰ δ¹⁸O/km is present on the Pacific slope of southern Costa Rica, equivalent to a temperature effect of −0.3‰/°C. The Nicoya and Osa Peninsulas have higher values than upwind sites, suggesting input of Pacific-sourced moisture, evaporative enrichment, or decreased condensation temperatures. Elevated and increasing d-excess values inland along the Nicaragua Trough suggest a recycled component may be an important contributor to the water budget. These data provide preliminary stable isotope information for Costa Rica, and will benefit paleoclimatic research in the region. More detailed studies would be beneficial to our understanding of the controls on stable isotope composition of tropical waters.     

Extracting flow structures from tracer data

A method of visualizing structures in closed chaotic flows out of homogenous particle distributions is presented in the example of models of a meandering jet. To this end, the system will be leaked or opened up by defining a region of the flow, so that a particle is considered to be escaped if it leaves this region. By applying this method to an ensemble of nonescaped tracers, we are able to characterize mixing processes by visualizing the converging and stretching filamentations (stable and unstable manifolds) in the flow without using additional mathematical tools. The possibility of applying the algorithm to analyze buoy data, and a comparison with the finite time manifolds are discussed.     

Diagenetic history and porosity evolution of Upper Carboniferous sandstones from the Spring Valley #1 well, Maritimes Basin, Canada–implications for reservoir development

Eighty-two core samples were collected from the Spring Valley #1 well which penetrates the Upper Carboniferous strata in the Late Devonian–Early Permian Maritimes Basin. The strata consist of alternating sandstones and mudstones deposited in a continental environment. The objective of this study is to characterize the relationship of sandstone porosity with depth, and to investigate the diagenetic processes related to the porosity evolution. Porosity values estimated from point counting range from 0% to 27.8%, but are mostly between 5% and 20%. Except samples that are significantly cemented by calcite, porosity values clearly decrease with depth. Two phases of calcite cement were distinguished based on Cathodoluminescence, with the early phase being largely dissolved and preserved as minor relicts in the later phase. Feldspar dissolution was extensive and contributed significantly to the development of secondary porosity. Quartz cementation was widespread and increased with depth. Fluid inclusions recorded in calcite and quartz cements indicate that interstitial fluids in the upper part of the stratigraphic column were dominated by waters with salinity lower than that of seawater, the middle part was first dominated by low-salinity waters, then invaded by brines, and the lower part was dominated by brines. Homogenization temperatures of fluid inclusions generally increase with depth and suggest a paleogeothermal gradient of 25 °C/km, which is broadly consistent with that indicated by vitrinite reflectance data. An erosion of 1.1–2.4 (mean 1.75) km of strata is inferred to have taken place above the stratigraphic column. δ¹⁸O values of calcite cements (mainly from the late phase) decrease with depth, implying increasing temperatures of formation, as also suggested by fluid-inclusion data. δ¹³C values of calcite cements range from −13.4‰ to −5.7‰, suggesting that organic matter was an important carbon source for calcite cements. A comparison of the porosity data with a theoretical compaction curve indicates that the upper and middle parts of the stratigraphic column show higher-than-normal porosity values, which are related to significant calcite and feldspar dissolution. Meteoric incursion and carboxylic acids generated from organic maturation were probably responsible for the abundant dissolution events.     

Hydrochemical and isotopic characteristics of groundwater in the Souss Upstream Basin, southwestern Morocco

Heterogeneous shallow Plio-Quaternary formations of the Souss Plain represent the most important aquifer in southern High Atlas Mountains in Morocco. The present work was conducted in the Souss Upstream Basin to identify the chemical characteristics and the origin of groundwater in an aquifer under semi-arid climate. Isotopic and hydrochemical compositions combined with geological and hydrogeological data were used for this purpose. The total dissolved solids vary from 239 to 997 mg l⁻¹, and the following groundwater types are recognized: Ca²⁺–Mg²⁺–HCO₃⁻, Ca²⁺–Mg²⁺–SO₄²⁻ and Ca²⁺–Mg²⁺–Cl⁻. The groundwater is saturated and slightly supersaturated with respect to carbonate minerals and undersaturated with respect to evaporite minerals, which means that the groundwater composition is largely controlled by the dissolution of carbonate rocks known in the basin. The isotopic contents of groundwaters ranged from −8‰ to −5.2‰ for δ¹⁸O, from −52‰ to −34‰ for δD, and from 0 to 5.5 TU for tritium. The hydrogen (δD) and oxygen (δ¹⁸O) isotope signatures reveal a significant infiltration before evaporation takes place, indicating a major recharge directly from fractures in the crystalline and limestone formations of Atlas Mountains (above 800 m a.s.l.) and infiltration of surface water in the alluvial cones at the border of the Atlas basins. The very low tritium values suggest that the groundwater recharge follows a long flow path and a mixing between old and modern water is shown. However, a slight evaporation effect is noted in the southern part of the basin close to the Anti-Atlas Mountains.     

A biogenic-silica δO record of climatic change during the last glacial–interglacial transition in southwestern Alaska

Despite growing evidence for environmental oscillations during the last glacial–interglacial transition from high latitude, terrestrial sites of the North Pacific rim, oxygen-isotopic records of these oscillations remain sparse. The lack of data is due partially to the paucity of lakes that contain carbonate sediment suitable for oxygen-isotopic analysis. We report here the first record of oxygen-isotopic composition in diatom silica (δ¹⁸O_Si) from a lake in that region. δ¹⁸O_Si increases gradually from 19.0 to 23.5‰ between 12,340 and 11,000 ¹⁴C yr B.P., reflecting marked climatic warming at the end of the last glaciation. Around 11,000 ¹⁴C yr B.P., δ¹⁸O_Si decreases by 1.7‰, suggesting a temperature decrease of 3.5–8.9 °C at the onset of the Younger Dryas (YD) in southwestern Alaska. Climatic recovery began ca. 10,740 ¹⁴C yr B.P., as inferred from the increase of δ¹⁸O_Si to a maximum of 23.9‰ near the end of the YD. Our data reveal that a YD climatic reversal in southwestern coastal areas of Alaska occurred, but the YD climate did not return to full-glacial conditions.