3D joint inversion of magnetotelluric and airborne tipper data: a case study from the Morrison porphyry Cu–Au–Mo deposit,British Columbia,Canada

Z‐axis tipper electromagnetic and broadband magnetotelluric data were used to determine three‐dimensional electrical resistivity models of the Morrison porphyry Cu–Au–Mo deposit in British Columbia. Z‐axis tipper electromagnetic data are collected with a helicopter, thus allowing rapid surveys with uniform spatial sampling. Ground‐based magnetotelluric surveys can achieve a greater exploration depth than Z‐axis tipper electromagnetic surveys, but data collection is slower and can be limited by difficult terrain. The airborne Z‐axis tipper electromagnetic tipper data and the ground magnetotelluric tipper data show good agreement at the Morrison deposit despite differences in the data collection method, spatial sampling, and collection date. Resistivity models derived from individual inversions of the Z‐axis tipper electromagnetic tipper data and magnetotelluric impedance data contain some similar features, but the Z‐axis tipper electromagnetic model appears to lack resolution below a depth of 1 km, and the magnetotelluric model suffers from non‐uniform and relatively sparse spatial sampling. The joint Z‐axis tipper electromagnetic inversion solves these issues by combining the dense spatial sampling of the airborne Z‐axis tipper electromagnetic technique and the deeper penetration of the lower frequency magnetotelluric data. The resulting joint resistivity model correlates well with the known geology and distribution of alteration at the Morrison deposit. Higher resistivity is associated with the potassic alteration zone and volcanic country rocks, whereas areas of lower resistivity agree with known faults and sedimentary units. The pyrite halo and ≥0.3% Cu zone have the moderate resistivity that is expected of disseminated sulphides. The joint Z‐axis tipper electromagnetic inversion provides an improved resistivity model by enhancing the lateral and depth resolution of resistivity features compared with the individual Z‐axis tipper electromagnetic and magnetotelluric inversions. This case study shows that a joint Z‐axis tipper electromagnetic–magnetotelluric approach effectively images the interpreted mineralised zone at the Morrison deposit and could be beneficial in exploration for disseminated sulphides at other porphyry deposits.     

Overstrength of dowelled CLT connections under monotonic and cyclic loading

This paper presents an evaluation of overstrength based on an experimental study on dowelled connections in Cross Laminated Timber (CLT). Connection overstrength needs to be well understood in order to ensure that ductile system behaviour and energy dissipation can be achieved under seismic loading. Overstrength is defined as the difference between the code-based strength, using characteristic material strengths, and the 95th percentile of the true strength distribution. Many aspects contribute to total connection overstrength, which makes its definition challenging. In this study, half-hole embedment tests were performed on CLT to establish embedment strength properties and three point bending tests were performed to determine the fastener yield moment. Different connection layouts, making use of mild steel dowels and an internal steel plate, were tested under monotonic and cyclic loading to evaluate theoretically determined overstrength values and study the influence of cyclic loading on overstrength. Experimental results were compared with strength predictions from code provisions and analytical models for ductile response under monotonic loading. It was found that cyclic loading does not significantly influence overstrength for connections that respond in a mixed-mode ductile way indicating that in future more expedient monotonic test campaigns could be used. This work also provides further experimental data and theoretical considerations necessary for the estimation of a generally applicable overstrength factor for dowelled CLT connections.     

Simulating a century of soil erosion for agricultural catchment management

Agricultural land management requires strategies to reduce impacts on soil and water resources while maintaining food production. Models that capture the effects of agricultural and conservation practices on soil erosion and sediment delivery can help to address this challenge. Historic records of climatic variability and agricultural change over the last century also offer valuable information for establishing extended baselines against which to evaluate management scenarios. Here, we present an approach that combines centennial‐scale reconstructions of climate and agricultural land cover with modelling across four lake catchments in the UK where radiometric dating provides a record of lake sedimentation. We compare simulations using MMF‐TWI, a catchment‐scale model developed for humid agricultural landscapes that incorporates representation of seasonal variability in vegetation cover, soil water balance, runoff and sediment contributing areas. MMF‐TWI produced mean annual sediment exports within 9–20% of sediment core‐based records without calibration and using guide parameter values to represent vegetation cover. Simulations of land management scenarios compare upland afforestation and lowland field‐scale conservation measures to reconstructed historic baselines. Oak woodland versus conifer afforestation showed similar reductions in mean annual surface runoff (8–16%) compared to current moorland vegetation but a larger reduction in sediment exports (26–46 versus 4–30%). Riparian woodland buffers reduced upland sediment yields by 15–41%, depending on understorey cover levels, but had only minor effect on surface runoff. Planting of winter cover crops in the lowland arable catchment halved historic sediment exports. Permanent grass margins applied to sets of arable fields across 15% or more of the catchment led to further significant reduction in exports. Our findings show the potential for reducing sediment delivery at the catchment scale with land management interventions. We also demonstrate how MMF‐TWI can support hydrologically‐informed decision making to better target conservation measures in humid agricultural environments. Copyright © 2018 John Wiley & Sons, Ltd.     

Incorporating spatially heterogeneous infiltration capacity into hydrologic models with applications for simulating post‐wildfire debris flow initiation

Soils in post‐wildfire environments are often characterized by a low infiltration capacity with a high degree of spatial heterogeneity relative to unburned areas. Debris flows are frequently initiated by run‐off in recently burned steeplands, making it critical to develop and test methods for incorporating spatial variability in infiltration capacity into hydrologic models. We use Monte Carlo simulations of run‐off generation over a soil with a spatially heterogenous saturated hydraulic conductivity (K_s) to derive an expression for an aerially averaged saturated hydraulic conductivity ( ) that depends on the rainfall rate, the statistical properties of K_s, and the spatial correlation length scale associated with K_s. The proposed method for determining is tested by simulating run‐off on synthetic topography over a wide range of spatial scales. Results provide a simplified expression for an effective saturated hydraulic conductivity that can be used to relate a distribution of small‐scale K_s measurements to infiltration and run‐off generation over larger spatial scales. Finally, we use a hydrologic model based on to simulate run‐off and debris flow initiation at a recently burned catchment in the Santa Ana Mountains, CA, USA, and compare results to those obtained using an infiltration model based on the Soil Conservation Service Curve Number.     

High-resolution modelling of a coastal harbour in the presence of strong tides and significant river runoff

Ocean Dynamics - In the context of Canada’s Ocean Protection Plan (OPP), improved coastal and near-shore modelling is needed to enhance marine safety and emergency response capacity in the...     

Structural and Hydrogeological Controls on Hydrocarbon and Brine Migration into Drinking Water Aquifers in Southern New York

Environmental concerns regarding the potential for drinking water contamination in shallow aquifers have accompanied unconventional energy development in the northern Appalachian Basin. These activities have also raised several critical questions about the hydrogeological parameters that control the naturally occurring presence and migration of hydrocarbon gases in shallow aquifers within petroliferous basins. To interrogate these factors, we analyzed the noble gas, dissolved ion, and hydrocarbon gas (molecular and isotopic composition) geochemistry of 98 groundwater samples from south‐central New York. All samples were collected ?1km from unconventional drilling activities and sample locations were intentionally targeted based on their proximity to various types of documented fault systems. In agreement with studies from other petroliferous basins, our results show significant correlations between elevated levels of radiogenic [⁴He], thermogenic [CH₄], and dissolved ions (e.g., Cl, Br, Sr, Ba). In combination, our data suggest that faults have facilitated the transport of exogenous hydrocarbon‐rich brines from Devonian source rocks into overlying Upper Devonian aquifer lithologies over geologic time. These data conflict with previous reports, which conclude that hydrodynamic focusing regulates the occurrence of methane and salt in shallow aquifers and leads to elevated levels of these species in restricted flow zones within valley bottoms. Instead, our data suggest that faults in Paleozoic rocks play a fundamental role in gas and brine transport from depth, regulate the distribution of their occurrence in shallow aquifers, and influence the geochemistry of shallow groundwater in this petroliferous basin.     

Linking physiography and evaporation using the isotopic composition of river water in 16 Canadian boreal catchments

Highly seasonal boreal catchments are hydrologically complex and generally data poor and, hence, are ripe for investigation using tracer‐aided hydrologic models. The influence of physiography on isotopic metrics was assessed to identify the catchment characteristics dominating evaporative enrichment. A multiyear stable isotope of water dataset was collected at the outlets of 16 boreal catchments in central Canada ranging in area from 12 to 15,282 km². Physiographic characteristics were obtained through raster analysis of freely available land cover images, stream networks, and digital elevation models. Correlation analysis indicated that as the percentage coverage of open water increased, so too did the evaporative effects observed at the catchment outlet. Correlation to wetland metrics indicated that increasing the percentage coverage of wetlands can reduce or increase evaporative effects observed, depending on the isotopic metric used and the corresponding drainage density, catchment slope, and presence of headwater lakes. The slopes of river evaporative‐mixing lines appear to reflect multifaceted relationships, strongest between catchment slope, headwater lakes, and connected wetlands, whereas mean line‐conditioned excess is more directly linked to physiographic variables. Hence, the slopes of river evaporative‐mixing lines and mean line‐conditioned excess are not interchangeable metrics of evaporative enrichment in a catchment. Relationships identified appear to be independent of catchment scale. These results suggest that adequate inclusion of the distribution of open water throughout a catchment, adequate representation of wetland processes, catchment slope, and drainage density are critical characteristics to include in tracer‐aided hydrologic models in boreal environments in order to minimize structural uncertainty.     

Seismology software: state of the practice

We analyzed the state of practice for software development in the seismology domain by comparing 30 software packages on four aspects: product, implementation, design, and process. We found room for improvement in most seismology software packages. The principal areas of concern include a lack of adequate requirements and design specification documents, a lack of test data to assess reliability, a lack of examples to get new users started, and a lack of technological tools to assist with managing the development process. To assist going forward, we provide recommendations for a document-driven development process that includes a problem statement, development plan, requirement specification, verification and validation (V&V) plan, design specification, code, V&V report, and a user manual. We also provide advice on tool use, including issue tracking, version control, code documentation, and testing tools.