Indicator mineralogy of kimberlite boulders from eskers in the Kirkland Lake and Lake Timiskaming areas, Ontario, Canada

Sixteen kimberlite boulders were collected from three sites on the Munro and Misema River Eskers in the Kirkland Lake kimberlite field and one site on the Sharp Lake esker in the Lake Timiskaming kimberlite field. The boulders were processed for heavy-mineral concentrates from which grains of Mg-ilmenite, chromite, garnet, clinopyroxene and olivine were picked, counted and analyzed by electron microprobe. Based on relative abundances and composition of these mineral phases, the boulders could be assigned to six mineralogically different groups, five for the Kirkland Lake area and one for the Lake Timiskaming area. Their indicator mineral composition and abundances are compared to existing data for known kimberlites in both the Kirkland Lake and Lake Timiskaming areas. Six boulders from the Munro Esker form a compositionally homogeneous group (I) in which the Mg-ilmenite population is very similar to that of the A1 kimberlite, located 7–12 km N (up-ice), directly adjacent to the Munro esker in the Kirkland Lake kimberlite field. U–Pb perovskite ages of three of the group I boulders overlap with that of the A1 kimberlite. Three other boulders recovered from the same localities in the Munro Esker also show some broad similarities in Mg-ilmenite composition and age to the A1 kimberlite. However, they are sufficiently different in mineral abundances and composition from each other and from the A1 kimberlite to assign them to different groups (II–IV). Their sources could be different phases of the same kimberlite or—more likely—three different, hitherto unknown kimberlites up-ice of the sample localities along the Munro Esker in the Kirkland Lake kimberlite field. A single boulder from the Misema River esker, Kirkland Lake, has mineral compositions that do not match any of the known kimberlites from the Kirkland Lake field. This suggests another unknown kimberlite exists in the area up-ice of the Larder Lake pit along the Misema River esker. Six boulders from the Sharp Lake esker, within the Lake Timiskaming field, form a homogeneous group with distinct mineral compositions unmatched by any of the known kimberlites in the Lake Timiskaming field. U–Pb perovskite age determinations on two of these boulders support this notion. These boulders are likely derived from an unknown kimberlite source up-ice from the Seed kimberlite, 4 km NW of the Sharp Lake pit, since indicator minerals with identical compositions to those of the Sharp Lake boulders have been found in till samples collected down-ice from Seed. Based on abundance and composition of indicator minerals, most importantly Mg-ilmenite, and supported by U–Pb age dating of perovskite, we conclude that the sources of 10 of the 16 boulders must be several hitherto unknown kimberlite bodies in the Kirkland Lake and Lake Timiskaming kimberlite fields.     

Diet and habitat of toxodont megaherbivores (Mammalia, Notoungulata) from the late Quaternary of South and Central America

The toxodont megaherbivores Toxodon and Mixotoxodon were endemic to South and Central America during the late Quaternary. Isotopic signatures of 47 toxodont teeth were analyzed to reconstruct diet and ancient habitat. Tooth enamel carbon isotope data from six regions of South and Central America indicate significant differences in toxodont diet and local vegetation during the late Quaternary. Toxodonts ranged ecologically from C₃ forest browsers in the Amazon (mean δ¹³C = −13.4‰), to mixed C₃ grazers and/or browsers living either in C₃ grasslands, or mixed C₃ forested and grassland habitats in Honduras (mean δ¹³C = −9.3‰), Buenos Aires province, Argentina (δ¹³C = −8.7‰), and Bahia, Brazil (mean δ¹³C = −8.6‰), to predominantly C₄ grazers in northern Argentina (δ¹³C = −4.4‰), to specialized C₄ grazers in the Chaco of Bolivia (δ¹³C = −0.1‰). Although these toxodonts had very high-crowned teeth classically interpreted for grazing, the isotopic data indicate that these megaherbivores had the evolutionary capacity to feed on a variety of dominant local vegetation. In the ancient Amazon region, carbon isotope data for the toxodonts indicate a C₃-based tropical rainforest habitat with no evidence for grasslands as would be predicted from the Neotropical forest refugia hypothesis.     

Quantifying flows along hydrological pathways by applying a new filtering algorithm in conjunction with master recession curve analysis

Quantifying the proportion of the river hydrograph derived from the different hydrological pathways is essential for understanding the behaviour of a catchment. This paper describes a new approach using the output from master recession curve analysis to inform a new algorithm based on the Lyne and Hollick ‘one‐parameter’ signal analysis filtering algorithm. This approach was applied to six catchments (including two subcatchments of these) in Ireland. The conceptual model for each catchment consists of four main flow pathways: overland flow, interflow, shallow groundwater and deep groundwater. The results were compared with those of the master recession curve analysis, a recharge coefficient approach developed in Ireland and the semi‐distributed, lumped and deterministic hydrological model Nedbør‐Afstrømings‐Model. The new algorithm removes the ‘free variable’ aspect that is typically associated with filtering algorithms and provides a means of estimating the contribution of each pathway that is consistent with the results of hydrograph separation in catchments that are dominated by quick response pathways. These types of catchments are underlain by poorly productive aquifers that are not capable of providing large baseflows in the river. Such aquifers underlie over 73% of Ireland, ensuring that this new algorithm is applicable in the majority of catchments in Ireland and potentially in those catchments internationally that are strongly influenced by the quick‐responding hydrological pathways. Copyright © 2013 John Wiley & Sons, Ltd.     

Characterizing spatial and temporal variation in 18O and 2H content of New Zealand river water for better understanding of hydrologic processes

Time series of hydrogen and oxygen stable isotope ratios (δ²H and δ¹⁸O) in rivers can be used to quantify groundwater contributions to streamflow, and timescales of catchment storage. However, these isotope hydrology techniques rely on distinct spatial or temporal patterns of δ²H and δ¹⁸O within the hydrologic cycle. In New Zealand, lack of understanding of spatial and temporal patterns of δ²H and δ¹⁸O of river water hinders development of regional and national-scale hydrological models. We measured δ²H and δ¹⁸O monthly, together with river flow rates at 58 locations across New Zealand over a two-year period. Results show: (a) general patterns of decreasing δ²H and δ¹⁸O with increasing latitude were altered by New Zealand's major mountain ranges; δ²H and δ¹⁸O were distinctly lower in rivers fed from higher elevation catchments, and in eastern rain-shadow areas of both islands; (b) river water δ²H and δ¹⁸O values were partly controlled by local catchment characteristics (catchment slope, PET, catchment elevation, and upstream lake area) that influence evaporation processes; (c) regional differences in evaporation caused the slope of the river water line (i.e., the relationship between δ²H and δ¹⁸O in river water) for the (warmer) North Island to be lower than that of the (cooler, mountain-dominated) South Island; (d) δ²H seasonal offsets (i.e., the difference between seasonal peak and mean values) for individual sites ranged from 0.50‰ to 5.07‰. Peak values of δ¹⁸O and δ²H were in late summer, but values peaked 1 month later at the South Island sites, likely due to greater snow-melt contributions to streamflow. Strong spatial differences in river water δ²H and δ¹⁸O caused by orographic rainfall effects and evaporation may inform studies of water mixing across landscapes. Generally distinct seasonal isotope cycles, despite the large catchment sizes of rivers studied, are encouraging for transit time analysis applications.     

Identification of Uranyl Minerals Using Oxygen K‐Edge X‐Ray Absorption Spectroscopy

Although most of the world's uranium exists as pitchblende or uraninite, this mineral can be weathered to a great variety of secondary uranium minerals, most containing the uranyl cation. Anthropogenic uranium compounds can also react in the environment, leading to spatial–chemical alterations that could be useful for nuclear forensics analyses. Soft X‐ray absorption spectroscopy (XAS) has the advantages of being non‐destructive, element‐specific and sensitive to electronic and physical structure. The soft X‐ray probe can also be focused to a spot size on the order of tens of nanometres, providing chemical information with high spatial resolution. However, before XAS can be applied at high spatial resolution, it is necessary to find spectroscopic signatures for a variety of uranium compounds in the soft X‐ray spectral region. To that end, we collected the near edge X‐ray absorption fine structure (NEXAFS) spectra of a variety of common uranyl‐bearing minerals, including uranyl carbonates, oxyhydroxides, phosphates and silicates. We find that uranyl compounds can be distinguished by class (carbonate, oxyhydroxide, phosphate or silicate) based on their oxygen K‐edge absorption spectra. This work establishes a database of reference spectra for future spatially resolved analyses. We proceed to show scanning X‐ray transmission microscopy (STXM) data from a schoepite particle in the presence of an unknown contaminant.     

Stochastic groundwater modeling of a sedimentary aquifer: evaluation of the impacts of abstraction scenarios under conditions of reduced recharge

A transient finite difference groundwater flow model has been calibrated for the Nasia sub-catchment of the White Volta Basin. This model has been validated through a stochastic parameter randomization process and used to evaluate the impacts of groundwater abstraction scenarios on resource sustainability in the basin. A total of 1500 equally likely model realizations of the same terrain based on 1500 equally likely combinations of the data of the key aquifer input parameters were calibrated and used for the scenario analysis. This was done to evaluate model non-uniqueness arising from uncertainties in the key aquifer parameters especially hydraulic conductivity and recharge by comparing the realizations and statistically determining the degree to which they differ from each other. Parameter standard deviations, computed from the calibrated data of the key parameters of hydraulic conductivity and recharge, were used as a yardstick for evaluating model non-uniqueness. All model realizations suggest horizontal hydraulic conductivity estimates in the range of 0.03–78.4 m/day, although over 70 % of the area has values in the range of 0.03–14 m/day. Low standard deviations of the horizontal hydraulic conductivity estimates from the 1500 solutions suggest that this range adequately reflects the properties of the material in the terrain. Lateral groundwater inflows and outflows appear to constitute significant components of the groundwater budgets in the terrain, although estimated direct vertical recharge from precipitation amounts to about 7 % of annual precipitation. High potential for groundwater development has been suggested in the simulations, corroborating earlier estimates of groundwater recharge. Simulation of groundwater abstraction scenarios suggests that the domain can sustain abstraction rates of up to 200 % of the current estimated abstraction rates of 12,960 m³/day under the current recharge rates. Decreasing groundwater recharge by 10 % over a 20-year period will not significantly alter the results of this abstraction scenario. However, increasing abstraction rates by 300 % over the period with decreasing recharge by 10 % will lead to drastic drawdowns in the hydraulic head over the entire terrain by up to 6 m and could cause reversals of flow in most parts of the terrain.