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.     

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.     

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.     

Magnetic Resonance Imaging of Nonaqueous Phase Liquid Distribution in Unconsolidated Soils Using Discriminatory Freezing

Organic contaminants present as nonaqueous phase liquids (NAPLs) in the subsurface often pose a long-term risk to human health and the environment. Investigating the distribution of NAPLs in porous media remains a major challenge in risk assessment and management of contaminated sites. Conventional soil coring and monitoring wells have been widely used over past decades as the primary means of subsurface investigation to determine NAPL extent. Known limitations of conventional approaches have led us to explore an alternative or a complementary technique to provide high-quality information of NAPL source zone architecture. This work advances an imaging tool for a variety of organic NAPL contaminants in unconsolidated soils through magnetic resonance imaging (MRI) of frozen cores. Using trichloroethylene (TCE) and o-xylene as model species, we illustrate that discriminatory freezing of water, while keeping the NAPL in a liquid state, enables high-resolution qualitative delineation of NAPL distribution within porous media. This novel approach may help improve site conceptual models and consequentially lead to highly tailored, more efficient remedial measures.     

Geomorphological evolution of a debris-covered glacier surface

There exists a need to advance our understanding of debris-covered glacier surfaces over relatively short timescales due to rapid, climatically induced areal expansion of debris cover at the global scale, and the impact debris has on mass balance. We applied unpiloted aerial vehicle structure-from-motion (UAV-SfM) and digital elevation model (DEM) differencing with debris thickness and debris stability modelling to unravel the evolution of a 0.15 km² region of the debris-covered Miage Glacier, Italy, between June 2015 and July 2018. DEM differencing revealed widespread surface lowering (mean 4.1 ± 1.0 m a^-1; maximum 13.3 m a^-1). We combined elevation change data with local meteorological data and a sub-debris melt model, and used these relationships to produce high resolution, spatially distributed maps of debris thickness. These maps were differenced to explore patterns and mechanisms of debris redistribution. Median debris thicknesses ranged from 0.12 to 0.17 m and were spatially variable. We observed localized debris thinning across ice cliff faces, except those which were decaying, where debris thickened. We observed pervasive debris thinning across larger, backwasting slopes, including those bordered by supraglacial streams, as well as ingestion of debris by a newly exposed englacial conduit. Debris stability mapping showed that 18.2–26.4% of the survey area was theoretically subject to debris remobilization. By linking changes in stability to changes in debris thickness, we observed that slopes that remain stable, stabilize, or remain unstable between periods almost exclusively show net debris thickening (mean 0.07 m a^-1) whilst those which become newly unstable exhibit both debris thinning and thickening. We observe a systematic downslope increase in the rate at which debris cover thickens which can be described as a function of the topographic position index and slope gradient. Our data provide quantifiable insights into mechanisms of debris remobilization on glacier surfaces over sub-decadal timescales, and open avenues for future research to explore glacier-scale spatiotemporal patterns of debris remobilization. © 2020 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd     

Constraining sediment provenance for tsunami deposits using distributions of grain size and foraminifera from the Kujukuri coastline and shelf,Japan

Tsunami deposits preserved in the geological record provide a more comprehensive understanding of their patterns of frequency and intensity over longer timescales; but recognizing tsunami deposits can prove challenging due to post-depositional changes, lack of contrast between the deposits and surrounding sedimentary layers, and differentiating between tsunami and storm deposition. Modern baseline studies address these challenges by providing insight into modern spatial distributions that can be compared with palaeotsunami deposits. This study documents the spatial fingerprint of grain size and foraminifera from Hasunuma Beach and the Kujukuri shelf to provide a basis from which tsunami deposits can be interpreted. At Hasunuma Beach, approximately 50 km east of Tokyo, the spatial distribution of three common proxies (foraminiferal taxonomy, foraminiferal taphonomy and sediment grain size) for tsunami identification were mapped and clustered using Partitioning Around Medoids cluster analysis. Partitioning Around Medoids cluster analysis objectively discriminated two coastal zones corresponding to onshore and offshore sample locations. Results show that onshore samples are characterized by coarser grain sizes (medium to coarse sand) and higher abundances of Pararotalia nipponica (27 to 63%) than offshore samples, which are characterized by finer grain sizes (fine to medium sand), lower abundances of Pararotalia nipponica (2 to 19%) and Ammonia parkinsoniana (0 to 10%), higher abundances of planktonics (15 to 58%) and species with fragile tests including Uvigerinella glabra. When compared to grain-size and foraminiferal taxonomy, foraminiferal taphonomy; i.e. surface condition of foraminifera, a proxy not commonly used to identify tsunami deposits, was most effective in discriminating modern coastal zones (identified supratidal, intertidal and offshore environments) and determining sediment provenance for tsunami deposits at Kujukuri. This modern baseline study assists the interpretation of tsunami deposits in the geological record because it provides a basis for sediment provenance to be determined.