Criteria for interpreting kimberlite as coherent: insights from the Muskox and Jericho kimberlites (Nunavut,Canada)

The Jurassic Muskox and Jericho kimberlites (Northern Slave Province, Nunavut, Canada) contain a variety of facies exhibiting different geometries, contact relationships, internal organisation, country rock abundance and olivine shapes, although many have similar matrix/groundmass mineralogies and textures. Five facies are examined that either have characteristics consistent with coherent rocks in general (i.e. intrusive and extrusive non-fragmental rocks) or are mineralogically and texturally similar to kimberlite described as coherent (or apparent coherent). Three facies are interpreted as coherent on the basis of: (1) geological setting, (2) apparent-porphyritic texture, (3) sharp contacts with fragmental kimberlite, (4) relative abundance of elongate and unbroken olivine crystals and (5) paucity of country rock xenoliths, while the remaining two facies are interpreted as fragmental on the basis of: (1) the gradational contacts with demonstrably fragmental kimberlite, (2) relative abundance and range of sizes of country rock lithic clasts and (3) numerous broken olivine crystals. Comparisons are made with coherent and apparent-coherent kimberlite from the literature. Our three coherent facies are similar to literature reported coherent kimberlite dykes hosted in country rock (CK_d) in terms of internal organisation, low abundance of country rock xenoliths, and apparent-porphyritic texture. Conversely, our two fragmental facies share attributes with previously described pipe-filling coherent and apparent-coherent kimberlite (CK_pf) in terms of geometry, internal organisation and abundance of country rock xenoliths. We conclude that CK_d and most CK_pf, although similar in matrix/groundmass mineralogy and texture, can be distinguished on the basis of internal organisation, country rock lithic clast abundance, texture (e.g. apparent-porphyritic texture) and possibly olivine crystal shapes and suggest that fragmental kimberlite is more common than reported.     

A Study of Long-Term MTBE Attenuation in the Borden Aquifer, Ontario, Canada

In 1988 and 1989, a natural gradient tracer test was performed in the shallow, aerobic and aquifer at Canadian Forces Base (CFB) Borden. A mixture of ground water containing dissolved oxygenated gasoline was injected below the water table along with chloride (Cl-) as a conservative tracer. The migration of BTEX, MTBE, and Cl was monitored in detail for 16 moths. The mass of BTEX compounds in the plume diminished significantly with time due to intrinsic aerobic biodegradation, while MTBE showed only a small decrease in mass over the 16-month period. In 1995/96, a comprehensive ground water sampling program was undertaken to define the mass of MTBE still present in the aquifer. Since the plume had migrated into an unmonitored section of the Borden Aquifer, numerical modeling and geostatistical methods were applied to define an optimal sampling grid and to improve the level of confidence in the results. A drive point profiling system was used to obtain ground water samples. Numerical modeling with no consideration of degradation pedicted maximum concentrations in excess of 3000 μg/L; field sampling found maximum concentrations of less than 200 μg/L. A mass balance for the remaining MTBE mass in the aquifer eight years after injection showed that only 3% of the original mass remained. Sorption, volatilization, a biotic degradation, and plant uptake are not considered significant attenuation processes for the field conditions. Therefore, we suggest that biodegradation may have played a major role in the attenuation of MTBE within the Borden Aquifer.     

Geology of the BK9 kimberlite (Damtshaa,Botswana): implications for the formation of dark volcaniclastic kimberlite

The BK9 kimberlite consists of three overlapping pipes. It contains two dark varieties of massive volcaniclastic kimberlite, informally termed dark volcaniclastic kimberlite (DVK). DVK(ns) is present in the north and south pipes and is interbedded with lenses of basalt breccia at the margins of the pipes. DVK(c) is present within the central pipe where it is overlain by a sequence of basalt breccias with interbedded volcanogenic sediments. The features observed within the DVK units of the BK9 kimberlite provide strong evidence for gas fluidisation of the accumulating pyroclastic material. These include the massive interior of the pipes, marginal epiclastic units, well-dispersed country-rock xenoliths and small-scale heterogeneities in lithic clast abundance. The upper portions of the central pipe provide a record of the transition from pyroclastic eruption and infill to passive epiclastic infilling of the crater, after the eruption has ceased. The wall-rock of the BK9 kimberlite dips inwards and is interpreted as post pipe-fill subsidence of the adjacent country rock. The two DVK units contain interstitial, silt-sized pyroclasts. The DVK(ns) has a higher fraction of former melt and displays evidence of incipient welding, as a result of differences in eruption dynamics. These units demonstrate that whilst DVK is comparable in many respects to MVK and forms part of a spectrum of volcaniclastic rocks formed by fluidisation, it differs in frequently containing silt-sized particles and including agglutinated and welded varieties with a high melt fraction. The DVK varieties, studied here, also have a distinctive hydrothermal assemblage, resulting from the abundance of low-silica accidental lithic clasts. Both the hydrothermal alteration and the abundance of silt-sized particles contribute to the DVKs distinctive dark colour.     

Quantifying groundwater discharge to a small perennial stream in southern Ontario, Canada

A study of the interaction between groundwater and surface water was undertaken within a small agricultural watershed in southern Ontario, Canada. Groundwater contributions to streamflow were measured along a section of stream during baseflow conditions and during rainfall events. Four techniques were used to estimate the contribution of groundwater to the stream along a 450 m reach (three during baseflow and one during stormflow conditions). Under baseflow conditions, streamflow measurements using the velocity–area technique indicated that the net groundwater flux to the stream during the summer months was 10 ml s⁻¹ m⁻¹. Hydrometric measurements (i.e. hydraulic gradient and hydraulic conductivity) taken using mini-piezometers installed in the sediments beneath the stream resulted in net groundwater flux estimates that were four to five times lower. Seepage meters failed to provide any measurements of water flux into or out of the stream. Therefore, based on these results, the velocity–area technique gives the best estimate of groundwater discharge. Hydrograph separations were conducted using isotopic ratios and electrical conductivity on two large rainfall events with different antecedent moisture conditions in the catchment. Both events showed that pre-event water (generally considered groundwater) dominated streamflow and tile drain flow with 64%–80% of the total discharge contributed by pre-event water. High water table conditions within the catchment resulted in greater stream discharge and a greater contribution of event water in the streamflow than that observed under low water table conditions for similar intensity storm events. The results also showed that differences in riparian zone width, vegetation and surface saturation conditions between the upper and lower catchment can influence the relative magnitude of streamflow response from the two catchment areas.     

Reconstruction of a kimberlite eruption,using an integrated volcanological,geochemical and numerical approach: A case study of the Fox Kimberlite,NWT, Canada

An integrated approach involving volcanology, geochemistry and numerical modelling has enabled the reconstruction of the volcanic history of the Fox kimberlite pipe. The observed deposits within the vent include a basal massive, poorly sorted, matrix supported, lithic fragment rich, eruption column collapse lapilli tuff. Extensive vent widening during the climactic magmatic phase of the eruption led to overloading of the eruption column with cold dense country rock lithic fragments, dense juvenile pyroclasts and olivine crystals, triggering column collapse. > 40% dilution of the kimberlite by granodiorite country rock lithic fragments is observed both in the physical componentry of the rocks and in the geochemical signature, where enrichment in Al₂O₃ and Na₂O compared to average values for coherent kimberlite is seen. The wide, deep, open vent provided a trap for a significant proportion of the collapsing column material, preventing large scale run-away in the form of pyroclastic flow onto the ground surface, although minor flows probably also occurred. A massive to diffusely bedded, poorly sorted, matrix supported, accretionary-lapilli bearing, lithic fragment rich, lapilli tuff overlies the column collapse deposit providing evidence for a late phreatomagmatic eruption stage, caused by the explosive interaction of external water with residual magma. Correlation of pipe morphology and internal stratigraphy indicate that widening of the pipe occurred during this latter stage and a thick granodiorite cobble-boulder breccia was deposited. Ash- and accretionary lapilli-rich tephra, deposited on the crater rim during the late phreatomagmatic stage, was subsequently resedimented into the vent. Incompatible elements such as Nb are used as indicators of the proportion of the melt fraction, or kimberlite ash, retained or removed by eruptive processes. When compared to average coherent kimberlite the ash-rich deposits exhibit ~ 30% loss of fines whereas the column collapse deposit exhibits ~ 50% loss. This shows that despite the poorly sorted nature of the column collapse deposit significant elutriation has occurred during the eruption, indicating the existence of a high sustained eruption column. The deposits within Fox record a complex eruption sequence showing a transition from a probable violent sub-plinian style eruption, driven by instantaneous exsolution of magmatic volatiles, to a late phreatomagmatic eruption phase. Mass eruption rate and duration of the sub-plinian phase of the eruption have been determined based on the dimensions of milled country-rock boulders found within the intra-vent deposits. Calculations show a short lived eruption of one to eleven days for the sub-plinian magmatic phase, which is similar in duration to small volume basaltic eruptions. This is in general agreement with durations of kimberlite eruptions calculated using entirely different approaches and parameters, such as predictions of magma ascent rates in kimberlite dykes.     

Identification of changes in heavy rainfall events in Ontario,Canada

To assess whether changes in the frequency of heavy rainfall events are occurring over time, annual maximum records from 21 rainfall gauges in Ontario are examined using frequency analysis methods. Relative RMSE and related boxplots are used to characterize assessment for selecting distributions; the Gumbel distribution is verified as one of the most suitable distributions to provide accurate quantile estimates. Records were divided into two time periods, and tested using the Mann-Kendall test and lag-1 autocorrelations to ensure that data in each period are identically distributed. The confidence intervals of design rainfalls for each return period (2, 5, 10, and 25-year) are derived by using resampling method, and compared at 90 % confidence levels. The changes in heavy rainfall intensities are tested at gauges across the Province of Ontario. Several significant decreases in heavy rainfall intensities are identified in central and southern Ontario. Increases in heavy rainfall intensities are identified in gauges at Sioux Lookout and Belleville. The sensitivity analysis of changes identified with respect to the year of splitting indicates changes are occurring during the 1980s and 1990s.     

Mercury and methylmercury biogeochemistry in a thawing permafrost wetland complex,Northwest Territories,Canada

In arctic and sub‐arctic environments, mercury (Hg), more specifically toxic methylmercury (MeHg), is of growing concern to local communities because of its accumulation in fish. In these regions, there is particular interest in the potential mobilization of atmospherically deposited Hg sequestered in permafrost that is thawing at unprecedented rates. Permafrost thaw and the resulting ground surface subsidence transforms forested peat plateaus into treeless and permafrost‐free thermokarst wetlands where inorganic Hg released from the thawed permafrost and draining from the surrounding peat plateaus may be transformed to MeHg. This study begins to characterize the spatial distribution of MeHg in a peat plateau–thermokarst wetland complex, a feature that prevails throughout the wetland‐dominated southern margin of thawing discontinuous permafrost in Canada's Northwest Territories. We measured pore water total Hg, MeHg, dissolved organic matter characteristics and general water chemistry parameters to evaluate the role of permafrost thaw on the pattern of water chemistry. A gradient in vegetation composition, water chemistry and dissolved organic matter characteristics followed a toposequence from the ombrotrophic bogs near the crest of the complex to poor fens at its downslope margins. We found that pore waters in poor fens contained elevated levels of MeHg, and the water draining from these features had dissolved MeHg concentrations 4.5 to 14.5 times higher than the water draining from the bogs. It was determined through analysis of historical aerial images that the poor fens in the toposequence had formed relatively recently (early 1970s) as a result of permafrost thaw. Differences between the fens and bogs are likely to be a result of their differences in groundwater function, and this suggests that permafrost thaw in this landscape can result in hotspots for Hg methylation that are hydrologically connected to downstream ecosystems. Copyright © 2016 John Wiley & Sons, Ltd.     

Lacustrine Arcellinida (testate lobose amoebae) as bioindicators of arsenic concentration within the Yellowknife City Gold Project,Northwest Territories,Canada

Arcellinida (testate lobose amoebae) were examined in near-surface sediment samples from 30 lakes located within the Gold Terra Resource Corporation (GTRC) Yellowknife City Gold Project, Northwest Territories, Canada, to assess the applicability of using the group as a robust tool to assess Arsenic (As) contamination. Lake contamination by As is of concern in the Yellowknife area because it is derived from geogenic (bedrock) and anthropogenic sources (legacy pollution from former mining operations, particularly the Giant Mine [1948–2004]). Statistical multivariate analyses (cluster analysis, non-metric multidimensional scaling [NMDS], and redundancy analysis [RDA]) were performed on geochemical (inductively coupled plasma mass spectrometry [ICP-MS]) and organic (loss-on-ignition) data to identify geochemical and hydrogeological controls on the arcellinidan distribution in the lakes. Three distinct faunal assemblages were identified using Cluster analysis and NMDS: 1) Elevated Arsenic Assemblage (EAA; Approximately Unbiased (AU) p-value = 97 %; median As =272.4 mg/kg, range = 42–1353 mg/kg ; n = 9); 2) Centropyxid-Dominated Assemblage (CDA; AU p-value = 87 %; median As =169.7 mg/kg, range 60–233 mg/kg ; n = 9); and, 3) Difflugiid-Dominated Assemblage (DDA; AU p-value = 95 %; median As =61.3 mg/kg, range = 16–316 mg/kg ; n = 12). Six statistically significant controls on assemblage structure were identified using RDA. These were As, calcium, iron, strontium, water depth, and total organic carbon (Arcellinida variance explained = 36.2 %) — As is the most significant control (12.2 %; p-value < 0.002). Stress-tolerant centropyxids dominate in lakes characterized by elevated sedimentary As concentrations (64 % centropyxids in EAA; 40 % centropyxids in CDA), while stress-sensitive taxa thrived in samples associated with lower As concentrations (DDA). Our results corroborate the findings of previous Arcellinida studies in subarctic Canada and confirms the reliability of using Arcellinida as a reconnaissance tool for tracking As contamination in impacted lakes.     

Turkey Lakes Watershed,Ontario, Canada: 40 years of interdisciplinary whole-ecosystem research

The Turkey Lakes Watershed (TLW) study is a federal, interdepartmental study established in 1979 to investigate the effects of acid rain on terrestrial and aquatic ecosystems. The 10.5 km² watershed, located in the Eastern Temperate Mixed Forest on the Canadian Shield, has been the site of multidisciplinary studies on biogeochemical and ecological processes conducted across plot to catchment scales. The whole-ecosystem investigative approach was adopted from the outset and has allowed research to evolve from its original (and continuing) acidification focus to include investigations on the effects of climate change, forest harvesting and other forest ecosystem perturbations. The extensive scientific and support infrastructure allows for collection of a comprehensive data record essential for understanding long-term environmental trends. Data include atmospheric deposition, meteorology, stream hydrology and chemistry, soil, pore and ground water properties, understory and overstory vegetation, lake and outflow physical and chemical properties, and aquatic macroinvertebrate and fish community composition and abundance. These data have contributed to over 400 published research papers and graduate theses. The watershed has also figured prominently in many continent-wide comparisons advancing fundamental watershed theory. The knowledge gained at TLW has influenced pollutant emission and natural resource management policies provincially, nationally and internationally.     

Ammonium occurrence in a salinized lowland coastal aquifer (Ferrara,Italy)

To understand the spatial distribution and the possible sources of high ammonium concentration in the unconfined coastal aquifer of the Po River Delta lowland (Italy), a detailed characterization of inorganic nitrogen species was completed. A total of 59 boreholes were drilled to collect core samples and groundwater samples of the aquifer. Among them, four piezometers, located along a representative flow line, were chosen to monitor the inorganic nitrogen concentrations with depth. Ammonium is the prevalent nitrogen inorganic species in groundwater, and its concentration increases with depth and salinity. Very high ammonium concentrations are found in coincidence with peaty sediments in the salinized anoxic aquifer and in the low‐lying aquitard. In particular, the elevated ammonium concentration derives from mineralization of organic matter present in fine sediments deposited in paleo‐marsh environments. Copyright © 2012 John Wiley & Sons, Ltd.     

The freeze-thaw weathering regime at a section of the Niagara escarpment on the Bruce Peninsula,Southern Ontario,Canada

The bedrock freeze-thaw and moisture regimes at an actively eroding site on the Niagara Escarpment, Bruce Peninsula, southern Ontario, were monitored between December 1983 and April 1984, and the results compared with amounts of debris collected in an adjacent rock trap. Frost wedging in pre-existing rock fissures is the primary mechanism responsible for the observed rockfall events; hydration effects are negligible. Debris production was more closely related to the duration of the freezing leg of the freeze-thaw cycle than to intensity or to cycle frequency. Release also coincided with periods of high pore saturation (> 60 per cent) and the seepage of water from cracks and fissures. Fewer freeze-thaw cycles were recorded in the air than at 1 and 3.5 cm in the bedrock. The number of cycles that could be declared geomorphologically effective according to established temperature criteria was normally less than half the total number of freeze-thaw cycles recorded in both air and bedrock. Under the current temperature regime at the field site, few effective cycles are capable of penetrating more than 5 cm into the free face.     

Rainfall partitioning through a mixed cedar swamp and associated C and N fluxes in Southern Ontario,Canada

Partitioning of rainfall through a forest canopy into throughfall, stemflow, and canopy interception is a critical process in the water cycle, and the contact of precipitation with vegetated surfaces leads to increased delivery of solutes to the forest floor. This study investigates the rainfall partitioning over a growing season through a temperate, riparian, mixed coniferous‐deciduous cedar swamp, an ecosystem not well studied with respect to this process. Seasonal throughfall, stemflow, and interception were 69.2%, 1.5%, and 29.3% of recorded above‐canopy precipitation, respectively. Event throughfall ranged from a low of 31.5 ± 6.8% for a small 0.8‐mm event to a high of 82.9 ± 2.4% for a large 42.7‐mm event. Rain fluxes of at least 8 mm were needed to generate stemflow from all instrumented trees. Most trees had funnelling ratios <1.0, with an exponential decrease in funnelling ratio with increasing tree size. Despite this, stand‐scale funnelling ratios averaged 2.81 ± 1.73, indicating equivalent depth of water delivered across the swamp floor by stemflow was greater than incident precipitation. Throughfall dissolved organic carbon (DOC) and total dissolved nitrogen (TDN) averaged 26.60 ± 2.96 and 2.02 ± 0.16 mg L^?1, respectively, which were ~11 and three times above‐canopy rain levels. Stemflow DOC averaged 73.33 ± 7.43 mg L^?1, 35 times higher than precipitation, and TDN was 4.45 ± 0.56 mg L^?1, 7.5 times higher than rain. Stemflow DOC concentration was highest from Populus balsamifera and TDN greatest from Thuja occidentalis trees. Although total below‐canopy flux of TDN increased with increasing event size, DOC flux was greatest for events 20–30 mm, suggesting a canopy storage threshold of DOC was readily diluted. In addition to documenting rainfall partitioning in a novel ecosystem, this study demonstrates the excess carbon and nitrogen delivered to riparian swamps, suggesting the assimilative capacity of these zones may be underestimated.     

Transport of three herbicides in ground water at Twin Lake test site,Chalk River,Ontario, Canada

A field tracer test performed under natural flow conditions at the Twin Lake test site, Chalk River Laboratories of the Atomic Energy of Canada Ltd. in Chalk River, Ontario, Canada, using tritium and three herbicides (Chlortoluron, Terbuthylazine, and Pendimethalin) was interpreted using the dispersion equation with a combined reaction model. The reaction model couples an instantaneous equilibrium reaction governed by a linear adsorption isotherm with a reversible or irreversible kinetic reaction of the first order, and decay. An improved interpretation method consists of a simultaneous fitting of theoretical concentration and mass-recovery curves to the experimental data, which leads to a more reliable determining of reaction models and improves the accuracy of fitting. Tritium served as the reference tracer to determine the flow velocity, dispersivity, and the recovery of the herbicides. Chlortoluron was slightly delayed by equilibrium exchange with strongly reduced concentration due to an irreversible kinetic reaction and/or decay. Terbuthilazine was slightly delayed by equilibrium exchange, with strongly reduced concentration due to a reversible kinetic reaction with some influence of decay. A strong equilibrium reaction and a strong reversible kinetic reaction without degradation governed the transport of Pendimethalin, reducing considerably its concentration. The results obtained show that simulations based only on Kd and decay constant, especially if these parameters are found in the laboratory, may considerably differ from those performed with reaction parameters determined in properly performed field tests. The dominant reaction types, and the values of parameters found in the study, supply useful information on the transport of the investigated herbicides in sandy aquifers under natural flow conditions.     

Antecedent soil moisture conditions of different soil types in South‐western Ontario,Canada

Soil moisture conditions prior to input design storms need to be known in the planning and design of urban stormwater control facilities using the design storm approach. Limited information is available on these soil moisture conditions which are commonly referred to as the antecedent soil moisture conditions. In this study, a deterministic continuous simulation model was used to simulate antecedent soil moisture conditions under south‐western Ontario, Canada, climate conditions. A wide range of different soil types were investigated and various statistical analyses on the simulated antecedent soil moisture results were performed. Frequency analyses illustrated typical distributions of antecedent soil moisture conditions and the influence of finer and coarser textured soil particles. Empirical equations were developed for the estimation of average antecedent soil moisture conditions based on commonly known soil characteristics. Satisfactory performance of the empirical equations was demonstrated by comparing between field average antecedent soil moisture data and empirically estimated average antecedent soil moisture values. These equations are therefore recommended for use in urban stormwater studies incorporating the design storm approach. Copyright © 2010 John Wiley & Sons, Ltd.     

A Method for Predicting Chloride Concentrations in Leachate at Natural Attenuation Landfills in the Precambrian Shield Regions of Ontario, Canada

Natural attenuation landfill sites continue to be the preferred method of domestic waste disposal in the Precambrian Shield regions of Ontario due to economic factors. The main challenge in siting these landfills is ensuring that there will be no adverse impact on off-site water resources. Impact risk assessments are generally based on estimated volumes and strengths of chloride in the leachate. While volumes can be estimated using simple water balances, peak chloride concentration predictions are based on judgment and are quite variable. Since design chloride strengths dictate the size of the required attenuation zone, overestimating concentrations will typically make it impossible to find a suitable site, while underestimating concentrations increases the potential for adverse off-site impacts occurring.
Hydrogeological data from active and closed landfills in the Precambrian Shield region were collected to help develop a reliable method of predicting peak chloride concentrations in leachate. This study focused on 21 sites located on relatively permeable sandy soils since landfills underlain by low permeability clayey soils retain leachate similar to lined facilities.
Linear regression analyses were conducted to determine if source chloride concentrations at the "sand" sites are significantly influenced by waste thickness, fill area, waste volume, waste deposition rate, hydraulic conductivity, upgradient flow length, depth to the water table, and moisture surplus.
A strong relationship (R = 0.957) was found to exist between source chloride concentrations and waste volume. This empirical volume versus chloride regression equation can be used as the basis for establishing design chloride concentrations at new natural attenuation landfills developed over sandy soils in the Precambrian Shield regions of Ontario. An alternative risk assessment approach is required for sites developed over clay soils.     

Geology of a complex kimberlite pipe (K2 pipe,Venetia Mine,South Africa): insights into conduit processes during explosive ultrabasic eruptions

K2 is a steep-sided kimberlite pipe with a complex internal geology. Geological mapping, logging of drillcore and petrographic studies indicate that it comprises layered breccias and pyroclastic rocks of various grain sizes, lithic contents and internal structures. The pipe comprises two geologically distinct parts: K2 West is a layered sequence of juvenile- and lithic-rich breccias, which dip 20–45° inwards, and K2 East consists of a steep-sided pipe-like body filled with massive volcaniclastic kimberlite nested within the K2 pipe. The layered sequence in K2 West is present to > 900 m below present surface and is interpreted as a sequence of pyroclastic rocks generated by explosive eruptions and mass-wasting breccias generated by rock fall and sector collapse of the pipe walls: both processes occurred in tandem during the infill of the pipe. Several breccia lobes extend across the pipe and are truncated by the steep contact with K2 East. Dense pyroclastic rocks within the layered sequence are interpreted as welded deposits. K2 East represents a conduit that was blasted through the layered breccia sequence at a late stage in the eruption. This phase may have involved fluidisation of trapped pyroclasts, with loss of fine particles and comminution of coarse clasts. We conclude that the K2 kimberlite pipe was emplaced in several distinct stages that consisted of an initial explosive enlargement, followed by alternating phases of accumulation and ejection.