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.     

Triangulation hand‐held laser‐scanning (TriHaLaS) for micro‐ and meso‐habitat surveys in streams

A novel hand‐held laser‐based stream bed survey system is presented. The system facilitates the capture of detailed 3D mapping of shallow (< 0.7 m) riverbed topography in sections approximately 4 m by 2 m. The system includes a stationary reference system, which projects three laser sheets (two at offset angles), within which a hand‐held monitoring pole is moved. The unique configuration of the light sheets intercepts with the pole as it moves within the survey area providing an exact horizontal location. Pole tilt is compensated for by an inertial measurement unit on the pole, and the height above the bed of the pole and submerged scanning laser are monitored relative to the horizontal laser sheet. Verification and application measurements demonstrate high resolution and accuracy in the horizontal (~5 mm) and vertical (~1 mm) direction. The system can be applied at sites where a free view is blocked and other optical through‐water methods fail. It is appropriate for studies on riverbed statistics and dynamics, which necessitate non‐invasive in‐situ surveys. Copyright © 2017 John Wiley & Sons, Ltd.     

Revisiting “An Exercise in Groundwater Model Calibration and Prediction” After 30 Years: Insights and New Directions

In 1988, an important publication moved model calibration and forecasting beyond case studies and theoretical analysis. It reported on a somewhat idyllic graduate student modeling exercise where many of the system properties were known; the primary forecasts of interest were heads in pumping wells after a river was modified. The model was calibrated using manual trial-and-error approaches where a model's forecast quality was not related to how well it was calibrated. Here, we investigate whether tools widely available today obviate the shortcomings identified 30 years ago. A reconstructed version of the 1988 true model was tested using increasing parameter estimation sophistication. The parameter estimation demonstrated the inverse problem was non-unique because only head data were available for calibration. When a flux observation was included, current parameter estimation approaches were able to overcome all calibration and forecast issues noted in 1988. The best forecasts were obtained from a highly parameterized model that used pilot points for hydraulic conductivity and was constrained with soft knowledge. Like the 1988 results, however, the best calibrated model did not produce the best forecasts due to parameter overfitting. Finally, a computationally frugal linear uncertainty analysis demonstrated that the single-zone model was oversimplified, with only half of the forecasts falling within the calculated uncertainty bounds. Uncertainties from the highly parameterized models had all six forecasts within the calculated uncertainty. The current results outperformed those of the 1988 effort, demonstrating the value of quantitative parameter estimation and uncertainty analysis methods.     

Groundwater Modeling with Nonlinear Uncertainty Analyses to Enhance Remediation Design Confidence

Soil and groundwater contamination are often managed by establishing on‐site cleanup targets within the context of risk assessment or risk management measures. Decision‐makers rely on modeling tools to provide insight; however, it is recognized that all models are subject to uncertainty. This case study compares suggested remediation requirements using a site‐specific numerical model and a standardized analytical tool to evaluate risk to a downgradient wetland receptor posed by on‐site chloride impacts. The base case model, calibrated to observed non‐pumping and pumping conditions, predicts a peak concentration well above regulatory criteria. Remediation scenarios are iteratively evaluated to determine a remediation design that adheres to practical site constraints, while minimizing the potential for risk to the downgradient receptor. A nonlinear uncertainty analysis is applied to each remediation scenario to stochastically evaluate the risk and find a solution that meets the site‐owner risk tolerance, which in this case required a risk‐averse solution. This approach, which couples nonlinear uncertainty analysis with a site‐specific numerical model provides an enhanced level of knowledge to foster informed decision‐making (i.e., risk‐of‐success) and also increases stakeholder confidence in the remediation design.     

Using the size and position of drumlins to understand how they grow,interact and evolve

Drumlins are subglacial bedforms streamlined in the direction of ice flow. Common in deglaciated landscapes, they have been widely studied providing rich information on their internal geology, size, shape, and spacing. In contrast with bedform investigations elsewhere in geomorphology (aeolian and fluvial dunes and ripples for example) most drumlin studies derive observations from relict, and thus static features. This has made it difficult to gain information and insights about their evolution over time, which likely hampers our understanding of the process(es) of drumlin formation. Here we take a morphological approach, studying drumlin size and spacing metrics. Unlike previous studies which have focussed on databases derived from entire ice sheet beds, we adopt a space‐for‐time substitution approach using individual drumlin flow‐sets distributed in space as proxies for different development times/periods. Framed and assisted by insights from aeolian and fluvial geomorphology, we use our metric data to explore possible scenarios of drumlin growth, evolution and interaction. We study the metrics of the size and spacing of 36 222 drumlins, distributed amongst 71 flow‐sets, left behind by the former British‐Irish Ice Sheet, and ask whether behaviour common to other bedform phenomena can be derived through statistical analysis. Through characterizing and analysing the shape of the probability distribution functions of size and spacing metrics for each flow‐set we argue that drumlins grow, and potentially migrate, as they evolve leading to pattern coarsening. Furthermore, our findings add support to the notion that no upper limit to drumlin size exists, and to the idea that perpetual coarsening could occur if given sufficient time. We propose that the framework of process and patterning commonly applied to non‐glacial bedforms is potentially powerful for understanding drumlin formation and for deciphering glacial landscapes. © 2017 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Experiments on the morphological controls of velocity inversions in bedrock canyons

A better understanding of bedrock incision mechanisms and processes is essential to the study of long‐term landscape evolution. Yet, little is known about flow dynamics in bedrock rivers, limiting our ability to make realistic predictions of local bedrock incision rates. A recent investigation of flow through bedrock canyons of the Fraser River revealed that plunging flows, defined by the downward‐directed movement of near surface flow toward the channel bed, occur in channels that have low width‐to‐depth ratios. Plunging flows occur into deep scour pools, which are often coincident with lateral constrictions and channel spanning submerged ridges (sills). A phenomenological investigation was undertaken to reproduce the flow fields observed in the Fraser canyons and to explore morphological controls on the occurrence and relative strength of plunging flow in bedrock canyons. Our observations show that the plunging flow structure can be produced along a scour pool entrance slope by accelerating the flow at the canyon entrance either over submerged sills or through lateral constrictions. Plunging flow appears to be a function of convective deceleration into a scour pool which can be enhanced by sill height, the amount of the channel width that is constricted, pool entrance slope, discharge, and a reduction in channel width‐to‐depth ratio. Plunging flow greatly enhances the potential for incision to occur along the channel bed and is an extreme departure from the assumptions of steady, uniform flow in bedrock incision models, highlighting the need for improved formulations that account for fluid flow. Copyright © 2017 John Wiley & Sons, Ltd.     

Sensitivity of Simulated Hyporheic Exchange to River Bathymetry: The Steinlach River Test Site

This study determines the aspects of river bathymetry that have the greatest influence on the predictive biases when simulating hyporheic exchange. To investigate this, we build a highly parameterized HydroGeoSphere model of the Steinlach River Test Site in southwest Germany as a reference. This model is then modified with simpler bathymetries, evaluating the changes to hyporheic exchange fluxes and transit time distributions. Results indicate that simulating hyporheic exchange with a high-resolution detailed bathymetry using a three-dimensional fully coupled model leads to nested multiscale hyporheic exchange systems. A poorly resolved bathymetry will underestimate the small-scale hyporheic exchange, biasing the simulated hyporheic exchange towards larger scales, thus leading to overestimates of hyporheic exchange residence times. This can lead to gross biases in the estimation of a catchment's capacity to attenuate pollutants when extrapolated to account for all meanders along an entire river within a watershed. The detailed river slope alone is not enough to accurately simulate the locations and magnitudes of losing and gaining river reaches. Thus, local bedforms in terms of bathymetric highs and lows within the river are required. Bathymetry surveying campaigns can be more effective by prioritizing bathymetry measurements along the thalweg and gegenweg of a meandering channel. We define the gegenweg as the line that connects the shallowest points in successive cross-sections along a river opposite to the thalweg under average flow conditions. Incorporating local bedforms will likely capture the nested nature of hyporheic exchange, leading to more physically meaningful simulations of hyporheic exchange fluxes and transit times.     

An examination of the rheology of flocculated clay suspensions

A dense cohesive sediment suspension, sometimes referred to as fluid mud, is a thixotropic fluid with a true yield stress. Current rheological formulations struggle to reconcile the structural dynamics of cohesive sediment suspensions with the equilibrium behaviour of these suspensions across the range of concentrations and shear. This paper is concerned with establishing a rheological framework for the range of sediment concentrations from the yield point to Newtonian flow. The shear stress equation is based on floc fractal theory, put forward by Mills and Snabre (1988). This results in a Casson-like rheology equation. Additional structural dynamics is then added, using a theory on the self-similarity of clay suspensions proposed by Coussot (1995), giving an equation which has the ability to match the equilibrium and time-dependent viscous rheology of a wide range of suspensions of different concentration and mineralogy.     

Using UAV acquired photography and structure from motion techniques for studying glacier landforms: application to the glacial flutes at Isfallsglaciären

Glacier and ice sheet retreat exposes freshly deglaciated terrain which often contains small‐scale fragile geomorphological features which could provide insight into subglacial or submarginal processes. Subaerial exposure results in potentially rapid landscape modification or even disappearance of the minor‐relief landforms as wind, weather, water and vegetation impact on the newly exposed surface. Ongoing retreat of many ice masses means there is a growing opportunity to obtain high resolution geospatial data from glacier forelands to aid in the understanding of recent subglacial and submarginal processes. Here we used an unmanned aerial vehicle to capture close‐range aerial photography of the foreland of Isfallsglaciären, a small polythermal glacier situated in Swedish Lapland. An orthophoto and a digital elevation model with ~2 cm horizontal resolution were created from this photography using structure from motion software. These geospatial data was used to create a geomorphological map of the foreland, documenting moraines, fans, channels and flutes. The unprecedented resolution of the data enabled us to derive morphological metrics (length, width and relief) of the smallest flutes, which is not possible with other data products normally used for glacial landform metrics mapping. The map and flute metrics compare well with previous studies, highlighting the potential of this technique for rapidly documenting glacier foreland geomorphology at an unprecedented scale and resolution. The vast majority of flutes were found to have an associated stoss‐side boulder, with the remainder having a likely explanation for boulder absence (burial or erosion). Furthermore, the size of this boulder was found to strongly correlate with the width and relief of the lee‐side flute. This is consistent with the lee‐side cavity infill model of flute formation. Whether this model is applicable to all flutes, or multiple mechanisms are required, awaits further study. © 2016 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

The dynamics of mountain erosion: Cirque growth slows as landscapes age

Glacial cirques are widely used palaeoenvironmental indicators, and are key to understanding the role of glaciers in shaping mountain topography. However, notable uncertainty persists regarding the rate and timing of cirque erosion. In order to address this uncertainty, we analyse the dimensions of 2208 cirques in Britain and Ireland and model ice accumulation to investigate the degree of coupling between glacier occupation times and cirque growth. Results indicate that during the last ~120 ka, cirques were glacier-free for an average of 52.0 ± 21.2 ka (43 ± 18%); occupied by small (largely cirque-confined) glaciers for 16.2 ± 9.9 ka (14 ± 8%); and occupied by large glaciers, including ice sheets, for 51.8 ± 18.6 ka (43 ± 16%). Over the entire Quaternary (i.e. 2.6 Ma), we estimate that cirques were glacier-free for 1.1 ± 0.5 Ma; occupied by small glaciers for 0.3 ± 0.2 Ma; and occupied by large glaciers for 1.1 ± 0.4 Ma. Comparing occupation times to cirque depths, and calculating required erosion rates, reveals that continuous cirque growth during glacier occupation is unlikely. Instead, we propose that cirques attained much of their size during the first occupation of a non-glacially sculpted landscape (perhaps during the timeframe of a single glacial cycle). During subsequent glacier occupations, cirque growth may have slowed considerably, with the highest rates of subglacial erosion focused during periods of marginal (small glacier) glaciation. We propose comparatively slow rates of growth following initial cirque development because a ‘least resistance’ shape is formed, and as cirques deepen, sediment becomes trapped subglacially, partly protecting the bedrock from subsequent erosion. In support of the idea of rapid cirque growth, we present evidence from northern British Columbia, where cirques of comparable size to those in Britain and Ireland developed in less than 140 ka. © 2019 John Wiley & Sons, Ltd. © 2019 John Wiley & Sons, Ltd.