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.     

Evaluating suspended sediment and turbidity reduction projects in a glacially conditioned catchment through dynamic regression and fluvial process-based modelling

Elevated turbidity (T_n) and suspended sediment concentrations (SSC) during and following flood events can degrade water supply quality and aquatic ecosystem integrity. Streams draining glacially conditioned mountainous terrain, such as those in the Catskill Mountains of New York State, are particularly susceptible to high levels of T_n and SSC sourced from erosional contact with glacial-related sediment. This study forwards a novel approach to evaluate the effectiveness of stream restoration best management practices (BMPs) meant to reduce stream T_n and SSC, and demonstrates the approach within the Stony Clove sub-basin of the Catskills, a water supply source for New York City. The proposed approach is designed to isolate BMP effects from natural trends in T_n and SSC caused by trends in discharge and shifts in average T_n or SSC per unit discharge (Q) following large flood events. We develop Dynamic Linear Models (DLMs) to quantify how T_n-Q and SSC-Q relationships change over time at monitoring stations upstream and downstream of BMPs within the Stony Clove and in three other sub-basins without BMPs, providing observational evidence of BMP effectiveness. A process-based model, the River Erosion Model, is then developed to simulate natural, hydrology-driven SSC-Q dynamics in the Stony Clove sub-basin (absent of BMP effects). We use DLMs to compare the modelled and observed SSC-Q dynamics and isolate the influence of the BMPs. Results suggest that observed reductions in SSC and T_n in the Stony Clove sub-basin have been driven by a combination of declining streamflow and the installed BMPs, confirming the utility of the BMPs for the monitored hydrologic conditions.     

Assessing the institutional capacity to adapt to climate change: a case study in the Cambodian health and water sectors

Institutional capacity is an important element for climate change adaptation (CCA) and the development of such capacity is a great challenge in a Least Developed Country like Cambodia where resources are limited. An important first step to increasing capacity is via an understanding of the level of existing capacity; future priorities can then be subsequently identified. This study aimed to assess the capacity of organizations to implement climate change activities in Cambodia in order to provide such a basis for building capacity. Four elements of capacity were investigated in this research: (1) financial resources, (2) cooperation and coordination of stakeholders, (3) availability and quality of information on vulnerability and adaptation to climate change, and (4) the level of understanding of climate change vulnerability and adaptation. The data were collected through semistructured interviews with a wide range of government and non-government informants across a number of sectors. Results of the study showed that informants perceived capacity for CCA to be very constrained, especially in terms of financial resources and cooperation, and addressing these factors was ranked as the highest climate change capacity priority. Institutional capacity constraints were considered to relate more generally to weak governance of CCA. In light of our research findings, the absence of local higher education institutions in CCA activities should be addressed. The support of such institutions would provide an important mechanism to progress both capacity development as well as partnerships and coordination between different types of organizations and relevant sectors.

Policy relevance

Capacity for CCA within Cambodian health and water sectors was perceived to be very constrained across a range of interdependent factors. Increasing funding was ranked as the highest priority for building capacity for CCA; however, governance factors such as ‘improved cooperation’ were also ranked highly. Improving stakeholders' awareness of the availability of adaptation funds and resources, and their responsiveness to funding criteria, is an important implication of our research, as is improving the mobilization of local resources and the private sector. To address the issue of weak cooperation among stakeholders, improving the coordination function of the National Climate Change Committee (NCCC) regarding stakeholder engagement and capacity building is crucial. Ensuring that CCA activities are based on sound information and knowledge from across different disciplines and, importantly, include the perspectives of vulnerable people themselves, ultimately underpins and supports the realization of the above priorities.     

The influence of the radiation pressure force on possible critical surfaces in binary systems

Using a spherically symmetric approximation for the radiation pressure force to compute a possible critical surface for binary systems, previous authors found that the surface opens up at the far side of the companion. It is shown that this effect may be unreal, and could be a consequence of the simple approximation for the radiation pressure force. Due to the influence of the radiation force, mass will be lost over the whole surface of the star. In that way much mass could leave the system in massive binary systems. On the basis of evolutionary models, including mass loss by stellar wind, our results were applied on the X-ray binaries 3U 1700-37 and HD 77581.     

Stellar dynamics in young clusters: the formation of massive runaways and very massive runaway mergers

In the present paper we combine an N-body code that simulates the dynamics of young dense stellar systems with a massive star evolution handler that accounts in a realistic way for the effects of stellar wind mass loss. We discuss two topics.

1. The formation and the evolution of very massive stars (with masses >120 M_⊙) is followed in detail. These very massive stars are formed in the cluster core as a consequence of the successive (physical) collisions of the 10–20 most massive stars in the cluster (this process is known as ‘runaway merging’). The further evolution is governed by stellar wind mass loss during core hydrogen and core helium burning (the WR phase of very massive stars). Our simulations reveal that, as a consequence of runaway merging in clusters with solar and supersolar values, massive black holes can be formed, but with a maximum mass ≈70 M_⊙. In low-metallicity clusters, however, it cannot be excluded that the runaway-merging process is responsible for pair-instability supernovae or for the formation of intermediate-mass black holes with a mass of several 100 M_⊙.
2. Massive runaways can be formed via the supernova explosion of one of the components in a binary system (the Blaauw scenario), or via dynamical interaction of a single star and a binary or between two binaries in a star cluster. We explore the possibility that the most massive runaways (e.g. ζ Pup, λ Cep, BD+43°3654) are the product of the collision and merger of two or three massive stars.

     

Common-reflection-point time migration

Since the early days of seismic processing, time migration has proven to be a valuable tool for a number of imaging purposes. Main motivations for its widespread use include robustness with respect to velocity errors, as well as fast turnaround and low computation costs. In areas of complex geology, in which it has well-known limitations, time migration can still be of value by providing first images and also attributes, which can be of much help in further, more comprehensive depth migration. Time migration is a very close process to common-midpoint (CMP) stacking and, more recently, to zero-offset commonreflection- surface (CRS) stacking. In fact, Kirchhoff time migration operators can be readily formulated in terms of CRS parameters. In the nineties, several studies have shown advantages in the use of common-reflection-point (CRP) traveltimes to replace conventional CMP traveltimes for a number of stacking and migration purposes. In this paper, we follow that trend and introduce a Kirchhoff-type prestack time migration and velocity analysis algorithm, referred to as CRP time migration. The algorithm is based on a CRP operator together with optimal apertures, both computed with the help of CRS parameters. A field-data example indicates the potential of the proposed technique.     

Implementation of Tidal Constituent Interpolation Method for the Israeli Coastline

Many ship-borne geodetic surveys at sea, such as Global Navigation Satellite System (GNSS)-based sea surface height (SSH) observation, acoustic profiling of the bottom, and others, deal with a dynamic topography which undergoes several changes during the survey campaign (e.g., changes in tide, salinity and currents). Those changes affect the measurements and may causes for some variations in the results. There are several methods for tidal variations correction, being the most dominant phenomena, such as tidal zoning, tidal constituent interpolation or ocean tidal models. In this study, we have implemented the tidal constituent interpolation method for the Israeli coastline in order to assess its quality and determine whether it is suitable for use in this particular region. This paper depicts the interpolation method, discusses some difficulties in the implementation for the Israeli coast and presents results from exemplary processing. In addition, we compare the results to those obtained using global and regional tidal models.     

A log‐normal spectral analysis of inorganic grain‐size distributions from a Canadian boreal lake core: Towards refining depositional process proxy data from high latitude lakes

Better methods for interpreting grain‐size spectra will enhance current understanding of past transport–depositional processes. A high‐resolution inorganic grain‐size dataset has been measured from a freeze core extracted from ‘Alberta Lake E’ a boreal fresh water lake 40 km east of the Athabasca Oil Sands in north‐eastern Alberta, Canada. The grain‐size spectra are remarkably consistent throughout the core, exhibiting a structure comprising six persistent grain‐size distributions below ca 250 μm, plus a rare medium‐sand distribution. Automated deconvolution of the grain‐size spectra produced poor results. Constraining the modes of two of the distributions produced deconvolution solutions that were statistically excellent and consistent with the structure of each spectrum. Statistical analysis of the ‘constrained’ solutions indicates that deconvolution successfully extracted independent grain‐size populations. Conversely, the multimodal spectra generate traditional measures (for example, mean grain size) that are inconsistent combinations of different individual populations and thus are poor proxies of transport–depositional processes. Alberta Lake E is situated in a boreal wetland landscape where sediment delivery is dominated by overland flow transport during spring melt. This context means that the Alberta Lake E grain‐size spectra can be interpreted to reflect: (i) a bedload component transported during short‐duration high discharge events that reflect the intensity of the melt; and (ii) a finer suspended load component representing material whose magnitude is controlled by the volume of the spring melt. Stratigraphically, bedload and suspended load populations demonstrate different short‐wavelength and long‐wavelength cyclicity, suggesting that spring melt is likely to be driven by cyclic external forcing factors. The links between the grain‐size spectra and spring melt have potential for generating proxy records that better capture the external controls over spring melt in boreal systems and the risks associated with these energetic hydrodynamics. This is exemplified by the coarsest Alberta Lake E distributions, which indicate that more intense spring‐melt dynamics occurred in pre‐historical times.     

A Comparative Study of Five Reference Materials and the Lombard Meteorite for the Determination of the Platinum‐Group Elements and Gold by LA‐ICP‐MS

A range of independently characterised reference materials (RMs) for LA‐ICP‐MS, used for the determination of the platinum‐group elements (PGE) and Au in a sulfide matrix, were analysed and compared: 8b, PGE‐A, NiS‐3, Po727‐T1, Po724‐T and the Lombard meteorite. The newly developed RM NiS‐3 was used as the RM for the calibration of all LA‐ICP‐MS analyses and the measured concentrations of the other RMs compared against their published concentrations. This data were also used to assess the consistency of concentrations calibrated against the different RMs. It was found that Po727‐T1 and 8b produced results that were comparable, within uncertainty, for all elements. Po727‐T1 also produced consistent results with NiS‐3 for all elements. All other RMs showed differences for some elements, especially Ru in Po724‐T, and Os, Ir and Au in PGE‐A. The homogeneity of the PGE and Au in each RM was assessed, by comparing the precision of multiple LA‐ICP‐MS spot analyses with the average uncertainty of the signal. Po724‐T, Po727‐T1 and the Lombard meteorite were found to be homogeneous for all elements, but 8b, PGE‐A and NiS‐3 were heterogeneous for some elements. This is the first direct comparison between a range of independently characterised PGE and Au LA‐ICP‐MS RMs.     

Multivariate statistical analysis of geochemical data as indicative of the hydrogeochemical evolution of groundwater in a sedimentary rock aquifer system

The study of groundwater hydrogeochemistry of the Paleozoic Basses-Laurentides sedimentary rock aquifer system in Québec produced a large geochemical dataset. Groundwater samples were collected at 153 sites over a 1500 km² study area and analyzed for major and minor ions. The large number of data can lead to difficulties in the integration, interpretation and representation of the results. Two multivariate statistical methods, hierarchical cluster analysis (HCA) and principal components analysis (PCA), were applied to a subgroup of the dataset to evaluate their usefulness to classify the groundwater samples, and to identify geochemical processes controlling groundwater geochemistry. This subgroup consisted of 144 samples and 14 parameters (Ca²⁺, Mg²⁺, Na⁺, K⁺, , Cl⁻, , Fe²⁺, Mn²⁺, Br⁻, Sr²⁺, F⁻, Ba²⁺, HS⁻). Seven geochemically distinct clusters, C1–C7, resulted from the HCA. Samples from clusters C3, C4, C6 and C7 are mostly located in preferential recharge areas. The majority of these samples have Ca–Mg–HCO₃ recharge groundwater (C3, C6, C7) and Na–HCO₃ evolved groundwater (C4). Samples from the other three clusters (C1, C2, C5) are characteristic of an aquifer system under confined conditions. The majority of these samples have Na–HCO₃ evolved groundwater (C1, C5) and Na–Cl ancient groundwater that exhibits elevated concentrations in Br⁻ (C2). In addition to recognizing the importance of hydrogeological conditions on groundwater geochemistry, the distribution of clusters also showed the importance of the geological formations on minor and trace elements, such as Fe²⁺, Mn²⁺, Sr²⁺, F⁻ and Ba²⁺. The first five components of the PCA account for 78.3% of the total variance in the dataset. Component 1 is defined by highly positive loadings in Na⁺, Cl⁻ and Br⁻ and is related to groundwater mixing with Champlain Sea water and solute diffusion from the marine clay aquitard. The high positive loadings in Ca²⁺ and Mg²⁺ of component 2 suggest the importance of dissolution of carbonate rocks in this aquifer system. From their characteristic loadings, the first two components are defined as the “salinity” and “hardness” components, respectively. Components 3–5 are related to more local and geological effects. The integration of the HCA and the PCA, with conventional classification of groundwater types, as well as with the hydrogeological and geological contexts, allowed the division of the region into four main geochemical areas, providing an improved regional picture of the aquifer system dynamics and hydrogeochemical evolution of groundwater. The following factors were recognized as influencing the evolution of groundwater identified in every geochemical area: (1) geological characteristics including sedimentary rock type and till mineralogy; (2) hydrogeological characteristics represented by the level of confinement and the hydraulic gradient; and (3) the geological history including the latest glaciation and the Champlain Sea invasion. With its integrated approach, this hydrogeochemical study contributes to the characterization and understanding of complex groundwater flow systems, and provides an example of the long-term geochemical evolution of hydrogeological systems after a major perturbation, in this case seawater invasion.