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.     

Controls on permafrost thaw in a coupled groundwater-flow and heat-transport system: Iqaluit Airport,Nunavut, Canada

Numerical simulations of groundwater flow and heat transport are used to provide insight into the interaction between shallow groundwater flow and thermal dynamics related to permafrost thaw and thaw settlement at the Iqaluit Airport taxiway, Nunavut, Canada. A conceptual model is first developed for the site and a corresponding two-dimensional numerical model is calibrated to the observed ground temperatures. Future climate-warming impacts on the thermal regime and flow system are then simulated based on climate scenarios proposed by the Intergovernmental Panel on Climate Change (IPCC). Under climate warming, surface snow cover is identified as the leading factor affecting permafrost degradation, including its role in increasing the sensitivity of permafrost degradation to changes in various hydrogeological factors. In this case, advective heat transport plays a relatively minor, but non-negligible, role compared to conductive heat transport, due to the significant extent of low-permeability soil close to surface. Conductive heat transport, which is strongly affected by the surface snow layer, controls the release of unfrozen water and the depth of the active layer as well as the magnitude of thaw settlement and frost heave. Under the warmest climate-warming scenario with an average annual temperature increase of 3.23 °C for the period of 2011–2100, the simulations suggest that the maximum depth of the active layer will increase from 2 m in 2012 to 8.8 m in 2100 and, over the same time period, thaw settlement along the airport taxiway will increase from 0.11 m to at least 0.17 m.

     

Meadow fragmentation and reproductive output of the SE Asian seagrass Enhalus acoroides

Flower and fruit production of the abundant, tall, long-lived, dioecious, surface-pollinating seagrass species Enhalus acoroides (L.) Royle were estimated at seven sites in the reef flats off Bolinao (NW Luzon, The Philippines) featuring different fragmentation of the seagrass meadows. Fragmentation of the seagrass meadow was quantified as cover of E. acoroides and all seagrass species present in 20×20 m plots. E. acoroides and overall seagrass cover were correlated positively. The proportion of female flowers of E. acoroides that developed a fruit increased sharply as overall seagrass cover was around 50%. Apparent sex ratio bore no relationship with overall seagrass cover. This threshold-type of relationship suggests that fragmentation of seagrass meadows can have a major effect on the reproductive output of this species. A possible mechanism underlying these results would be a non-linear increase of the efficiency of trapping the surface-dispersed pollen with increasing seagrass canopy density. This provides the first evidence based on real data that fragmentation can affect the population dynamics of seagrass species.     

Climate and climate change in western canada as simulated by the Canadian regional climate model

Abstract

A þrst climate simulation performed with the novel Canadian Regional Climate Model (CRCM) is presented. The CRCM is based on fully elastic non‐hydrostatic þeld equations, which are solved with an efþcient semi‐implicit semi‐Lagrangian (SISL) marching algorithm, and on the parametrization package of subgrid‐scale physical effects of the second‐generation Canadian Global Climate Model (GCMII). Two 5‐year integrations of the CRCM nested with GCMII simulated data as lateral boundary conditions are made for conditions corresponding to current and doubled CO₂ scenarios. For these simulations the CRCM used a grid size of 45 km on a polar‐stereographic projection, 20 scaled‐height levels and a time step of 15 min; the nesting GCMII has a spectral truncation of T32, 10 hybrid‐pressure levels and a time step of 20 min. These simulations serve to document: (1) the suitability of the SISL numerical scheme for regional climate modelling, (2) the use of GCMII physics at much higher resolution than in the nesting model, (3) the ability of the CRCM to add realistic regional‐scale climate information to global model simulations, and (4) the climate of the CRCM compared to that of GCMII under two greenhouse gases (GHG) scenarios.     

The Batu Hijau porphyry copper-gold deposit, Sumbawa Island, Indonesia

The Batu Hijau porphyry Cu---Au deposit lies in southwest Sumbawa Island, Indonesia. It is a world-class porphyry Cu deposit in an island are setting, and is typical of this deposit type in most features, including igneous association, morphology, hydrothermal alteration and mineralisation style.The region was not previously recognised as a porphyry Cu province; disseminated Cu sulphides were first recognised in float samples in southwest Sumbawa in 1987. Associated stream sediment sampling identified a broad area of anomalous Au and Cu in an area of greater than 5 km² around Batu Hijau, including 169 ppb Au in BLEG samples and 580 ppm Cu in stream silts 1 km from the deposit. Mineralisation in bedrock at surface contains > 0.1 wt % Cu and > 0.1 ppm Au over an area of 0.6 km × 1.2 km, including a zone 300 m × 900 m containing > 0.3 wt % Cu. Areas with elevated Mo (> 30 ppm) form a distinctive annulus around this Cu-rich zone.Batu Hijau mineralisation is hosted in a tonalite intrusive complex, and diorite and metavolcanic wallrocks. There are no post-mineralisation igneous intrusions or breccia pipes within the deposit. The main tonalite intrusion forms a stock in the centre of the deposit, where it generally displays intensely pervasive potassic (biotite with magnetite-quartz) alteration and hosts most of the higher grade mineralisation. Younger tonalite dykes intruding the centre of this stock are generally less altered and mineralised than the older tonalite.The core zone of potassic alteration grades outward into extensive propylitic alteration (chlorite-epidote), with both variably overprinted by widespread fracture controlled intermediate argillic alteration (sericite-chlorite), and minor phyllic (sericite-pyrite) and sodic (albite) alteration. Argillic (sericite-kaolinite) and advanced argillic (kaolinite-alunite-pyrophyllite) assemblages occur near surface.Copper and Au grades within the orebody show a positive correlation with quartz stockwork intensity, although disseminated Cu sulphides are also common. Chalcopyrite and bornite are the principle hypogenal minerals, with minor chalcocite. Oxidation extends to a depth of 5 m to 85 m below surface across the deposit, and is underlain by weak supergene mineralisation. Drill testing of the deposit down to 650 m below surface reveals a single cylindrical to conical orebody of 334 million tonnes grading 0.8 wt % Cu and 0.69 gm per tonne Au; the depth extent of mineralisation is unknown.     

On the flow structure and turbulence during sweep and ejection events in a wind-tunnel model canopy

Particle image velocimetry (PIV) data obtained in a wind-tunnel model of a canopy boundary layer is used to examine the characteristics of mean flow and turbulence. The vector spacing varies between 1.7 and 2.5 times the Kolmogorov scales. Conditional sampling based on quadrants, i.e. based on the signs of velocity fluctuations, reveals fundamental differences in flow structure, especially between sweep and ejection events, which dominate the flow. During sweeps, the downward flow generates a narrow, highly turbulent, shear layer containing multiple small-scale vortices just below canopy height. During ejections, the upward flow expands this shear layer and the associated small-scale flow structures to a broad region located above the canopy. Consequently, during sweeps the turbulent kinetic energy (TKE), Reynolds stresses, as well as production and dissipation rates, have distinct narrow peaks just below canopy height, whereas during ejections these variables have broad maxima well above the canopy. Three methods to estimate the dissipation rate are compared, including spectral fits, measured subgrid-scale (SGS) energy fluxes at different scales, and direct measurements of slightly underresolved instantaneous velocity gradients. The SGS energy flux is 40–60% of the gradient-based (direct) estimates for filter sizes inside the inertial range, while decreasing with scale, as expected, within the dissipation range. The spectral fits are within 5–30% of the direct estimates. The spectral fits exceed the direct estimates near canopy height, but are lower well above and below canopy height. The dissipation rate below canopy height increases with velocity magnitude, i.e. it has the highest values during sweep and quadrant 1 events, and is significantly lower during ejection and quadrant 3 events. Well above the canopy, ejections are the most dissipative. Turbulent transport during sweep events acts as a source below the narrow shear layer within the canopy and as a sink above it. Transport during ejection events is a source only well above the canopy. The residual term in the TKE transport equation, representing mostly the effect of pressure–velocity correlations, is substantial only within the canopy, and is dominated by sweeps.     

Surface Geophysical Methods used to Verify the Karst Geological Structure in the Built-up Area: A Case Study of Specific Engineering-Geological Conditions

This article presents a research study of complex limestone karst engineering-geological conditions in the municipality Valaská near Banská Bystrica in Slovakia. The aim of the study is to demonstrate the impossibility of spatial identification of cave spaces using surface geophysical methods due to the specific engineering-geological conditions of a thick surface layer of anthropogenic fill containing highly heterogeneous anthropogenic material. Its maximum thickness is 3 m. Another specificific condition of the study area is its location in the built-up area, due to which the applicability of geophysical methods was limited. The article contains methodological recommendations to be used in analogous geological conditions with karst structures topped with anthropogenic fill, which complicates the identification of cave spaces. The recommended solution herein is the identification of the cave system using underground mapping of the karst and its projection onto the surface for which surface geophysical methods have been combined.     

Yucatán karst features and the size of Chicxulub crater

The buried Chicxulub impact structure is marked by a dramatic ring of sinkholes (called cenotes if containing water), and adjacent less prominent partial rings, which have been shown to coincide with maxima in horizontal gravity gradients and a topographic depression. These observations, along with the discreteness and spacing of the features, suggest a formation mechanism involving faulting in the outer slump zone of the crater, which would thus have a diameter of approximately 180 km.
An opposing view, based primarily on the interpretation of gravity data, is that (he crater is much larger than the cenote ring implies. Given the association of the known cenote ring with faults, we here examine northern Yucatan for similar rings in gravity, surface features and elevation, which we might expect to be associated with outer concentric faults in the case of a larger, possibly multiring, structure.
No such outer rings have been found, although definite patterns are seen in the distribution of karst features outside the crater rim. We explain these patterns as resulting mainly from deformation related to the block fault zone that parallels tbe shelf edge of eastern Yucatan.     

shinyNORRRM: A Cross-Platform Software to Calculate the CIPW Norm

In this paper, a novelty-free software to assess an efficient CIPW Norm (± 0.006 wt.% in differences between input and output data) is presented. The package is available in the official repository for user-contributed R packages (CRAN: Comprehensive R Archive Network). The software is able to handle big data sets and considers minor and trace element compositions. The algorithm can calculate odd minerals in igneous rocks, such as cancrinite and calcite, adjust the Fe⁺³/Fe⁺² ratio in different standard approaches, and recalculate the compositions of the rocks in an anhydrous basis (100 ± 0.003 wt.% volatile-free adjusted). Furthermore, the package calculates several petrological parameters, and the graphical outputs are displayed following IUGS scheme standards. The prime aspect of shinyNORRRM is the symbiosis of native R functions with the R package’s shiny (Web Application Framework for R) to run the norm in a user-friendly interface. shinyNORRRM can be executed in any operating system and requires no previous programming knowledge, thus promising to be the universal computational program in this matter. The output data are printed in the standard comma-separated values (*.csv) format, which is highly compatible with general spreadsheet editors. In this work, the algorithm of our program is validated using already compiled whole-rock geochemical databases.

     

Geochemistry of suspended matter from the Baltic Sea 2. Results of bulk trace metal analysis by AAS

In 1984, on a transect covering the whole Baltic Sea and parts of the adjacent North Sea, 160 water samples were taken and analysed for their concentrations of particulate and dissolved metals. In addition, the suspended materials were investigated for their elemental bulk composition.The particulate fractions represented from about 5% (Cd, Cu and Ni) to 50% (Fe and Pb) of the total (particulate plus dissolved) concentrations. For some elements (Ba, Cd, Cu, Pb and Zn), the particulate matter from the surface microlayer was enriched with respect to those suspended materials taken from 0.2 m depth. This could reflect the atmospheric input of metal-rich aerosols. In anoxic deep waters, maximum contents of Zn (6400 μg g⁻¹), Cu (1330 μg g⁻¹) and Cd (12 μg g⁻¹) were observed in the particulate matter, indicating sulphidic forms. On the other hand, under oxic conditions the distribution coefficients (K_d) decreased with the water depth (Cd, Fe and Pb).Relative to global background levels, the particulate matter contained metal “excesses” amounting to more than 90% of the total contents (Cd, Mn, Pb and Zn). Automated electron probe X-ray microanalysis (EPXMA) revealed that the elemental composition of sediments is mainly governed by post-depositional processes of early diagenesis and is only weakly related to the composition of suspended matter in the overlying water body. For instance, in relation to surface mud sediments of the central Baltic net-sedimentation basins, Zn, Cd, Cu and Mn had 30–100% higher levels in the suspended materials. The general pattern of metal contents of particulate matter taken from 10 m depth on a transect between the Bothnian Bay and the North Sea were—possibly as a result of anthropogenic inputs—rather similar for Pb, Zn and Cu. For Fe and Mn, the distribution patterns along the transect were probably governed by the natural loading characteristics and by the biogeochemistry of those elements.