•3D estimation and visualization of waste volumes.•Understanding 3D effects in geophysical images.•Improving waste characterization by frequency-dependence of electrical properties.•First 3D ...frequency-dependence survey in landfills.•Resolving capabilities of different configurations for 3D electrical surveys.
There is growing interest in the use of spectral induced polarization (SIP) surveys to characterize the near-surface environment. Few attempts have been made to perform field SIP surveys in a 3D configuration; when done, they are typically conducted using a series of parallel 2D electrode lines with collinear measurements. However, such measurements are limited in the resolution between the lines which is critical in the case of heterogeneous subsurface conditions, such as in landfills. To overcome this, we investigate here the enhanced resolution in SIP measurements through true 3D measurements, i.e., the resolving capabilities of different electrode configurations distributed across measuring planes. First, we investigate, through a synthetic study, the difference between results from using 2D parallel collinear electrode arrays and true 3D configurations. Second, we collected SIP data (in the frequency range between 1 and 240 Hz) using 2D and 3D configurations in two landfills to evaluate the application of our results in real field conditions. Both the synthetic and the field experiments demonstrate that measurements of parallel 2D collinear arrays result in the creation of artifacts and the loss of resolution in the 3D structure, especially of polarizable features. In contrast, the 3D configurations are able to resolve polarizable anomalies in synthetic and field measurements, resulting in a better delineation of the geometry of waste units. Our results also demonstrate that 3D configurations are better suited to recover the frequency-dependence of the electrical properties; thus, permitting an improved interpretation of waste composition and the quantification of waste volume.
We propose a new approach for the inversion of anisotropic P-wave data based on Monte Carlo methods combined with a multigrid approach. Simulated annealing facilitates objective minimization of the ...functional characterizing the misfit between observed and predicted traveltimes, as controlled by the Thomsen anisotropy parameters (ɛ, δ). Cycling between finer and coarser grids enhances the computational efficiency of the inversion process, thus accelerating the convergence of the solution while acting as a regularization technique of the inverse problem. Multigrid perturbation samples the probability density function without the requirements for the user to adjust tuning parameters. This increases the probability that the preferred global, rather than a poor local, minimum is attained. Undertaking multigrid refinement and Monte Carlo search in parallel produces more robust convergence than does the initially more intuitive approach of completing them sequentially. We demonstrate the usefulness of the new multigrid Monte Carlo (MGMC) scheme by applying it to (a) synthetic, noise-contaminated data reflecting an isotropic subsurface of constant slowness, horizontally layered geologic media and discrete subsurface anomalies; and (b) a crosshole seismic data set acquired by previous authors at the Reskajeage test site in Cornwall, UK. Inverted distributions of slowness (s) and the Thomson anisotropy parameters (ɛ, δ) compare favourably with those obtained previously using a popular matrix-based method. Reconstruction of the Thomsen ɛ parameter is particularly robust compared to that of slowness and the Thomsen δ parameter, even in the face of complex subsurface anomalies. The Thomsen ɛ and δ parameters have enhanced sensitivities to bulk-fabric and fracture-based anisotropies in the TI medium at Reskajeage. Because reconstruction of slowness (s) is intimately linked to that ɛ and δ in the MGMC scheme, inverted images of phase velocity reflect the integrated effects of these two modes of anisotropy. The new MGMC technique thus promises to facilitate rapid inversion of crosshole P-wave data for seismic slownesses and the Thomsen anisotropy parameters, with minimal user input in the inversion process.
This study demonstrates the potential value of a combined unmanned aerial vehicle (UAV) Photogrammetry and ground penetrating radar (GPR) approach to map snow water equivalent (SWE) over large ...scales. SWE estimation requires two different physical parameters (snow depth and density), which are currently difficult to measure with the spatial and temporal resolution desired for basin‐wide studies. UAV photogrammetry can provide very high‐resolution spatially continuous snow depths (SD) at the basin scale, but does not measure snow densities. GPR allows nondestructive quantitative snow investigation if the radar velocity is known. Using photogrammetric snow depths and GPR two‐way travel times (TWT) of reflections at the snow‐ground interface, radar velocities in snowpack can be determined. Snow density (RSN) is then estimated from the radar propagation velocity (which is related to electrical permittivity of snow) via empirical formulas. A Phantom‐4 Pro UAV and a MALA GX450 HDR model GPR mounted on a ski mobile were used to determine snow parameters. A snow‐free digital surface model (DSM) was obtained from the photogrammetric survey conducted in September 2017. Then, another survey in synchronization with a GPR survey was conducted in February 2019 whilst the snowpack was approximately at its maximum thickness. Spatially continuous snow depths were calculated by subtracting the snow‐free DSM from the snow‐covered DSM. Radar velocities in the snowpack along GPR survey lines were computed by using UAV‐based snow depths and GPR reflections to obtain snow densities and SWEs. The root mean square error of the obtained SWEs (384 mm average) is 63 mm, indicating good agreement with independent SWE observations and the error lies within acceptable uncertainty limits.
Photogrammetric snow depths are combined with measurements of dielectric properties of the snowpack from a common offset GPR system to fully exploit the benefits of each for snow water equivalent (SWE) measurement. Radar velocities in the snowpack along GPR survey lines are computed by using UAV‐derived snow depths and GPR reflections to obtain snow densities and SWE.The root mean square error of the obtained SWEs are in good agreement with independent SWE observations demonstrating the potential value of this approach.
Forest planting is increasingly being incorporated into land management policies to mitigate diffuse pollution and localised flooding because forest soils are associated with enhanced hydraulic ...properties and lower surface runoff compared to soils under other vegetation types. Despite this, our understanding of the effects of different tree species and forest land use on soil hydraulic properties is limited. In this study we tested for the effects of two tree species, sycamore (Acer pseudoplatanus) and Scots pine (Pinus sylvestris), subject to contrasting land use systems, namely ungrazed forest and livestock grazed forest, on soil surface saturated hydraulic conductivity (Kfs) at a long term (23year) experimental site in Scotland. Additionally these forest land use systems were compared to grazed pasture. Kfs was found to be significantly higher under ungrazed Scots pine forest (1239mmhr−1) than under ungrazed sycamore forest (379mmhr−1) and under both of these forest types than under pasture (32mmhr−1). However, this measure did not differ significantly between the sycamore and Scots pine grazed forest and pasture. It was inferred, from comparison of measured Kfs values with estimated maximum rainfall intensities for various return periods at the site, that surface runoff, as infiltration excess overland flow, would be generated in pasture and grazed forest by storms with a return period of at least 1 in 2years, but that surface runoff is extremely rare in the ungrazed forests, regardless of tree species. We concluded that, although tree species with differing characteristics can create large differences in soil hydraulic properties, the influence of land use can mask the influence of trees. The choice of tree species may therefore be less important than forest land use for mitigating the effects of surface runoff.
•Forest land cover had greater soil hydraulic conductivity than pasture.•Contrasting tree species can mitigate the effects of surface runoff.•Return period indicates surface runoff would be extremely rare in ungrazed forest.•Livestock grazing under trees negates their mitigative potential against flooding.
SUMMARY Induced polarization (IP) has been widely used to non-invasively characterize electrical conduction and polarization in the subsurface resulting from an applied electric field. Earth ...materials exhibit a lossy capacitance defined by a negative intrinsic phase in frequency-domain IP (FDIP) or a positive intrinsic chargeability in time-domain IP (TDIP). However, error-free positive apparent phase or negative apparent chargeability (i.e. negative IP effects) can occur in IP measurements over heterogeneous media. While negative IP effects in TDIP data sets have been discussed, no studies have addressed this topic in detail for FDIP measurements. We describe theory and numerical modelling to explain the origin of negative IP effects in FDIP measurements. A positive apparent phase may occur when a relatively high polarizability feature falls into negative sensitivity zones of complex resistivity measurements. The polarity of the apparent phase is determined by the distribution of subsurface intrinsic phase and resistivity, with the resistivity impacting the apparent phase polarity via its control on the sensitivity distribution. A physical explanation for the occurrence of positive apparent phase data is provided by an electric circuit model representing a four-electrode measurement. We also show that the apparent phase polarity will be frequency dependent when resistivity changes significantly with frequency (i.e. in the presence of significant IP effects). Consequently, negative IP effects manifest themselves in the shape of apparent phase spectra recorded with multifrequency (spectral IP) data sets. Our results imply that positive apparent phase measurements should be anticipated and should be retained during inversion and interpretation of single frequency and spectral IP data sets.
SUMMARY
Induced polarization (IP) has been widely used to non-invasively characterize electrical conduction and polarization in the subsurface resulting from an applied electric field. Earth ...materials exhibit a lossy capacitance defined by a negative intrinsic phase in frequency-domain IP (FDIP) or a positive intrinsic chargeability in time-domain IP (TDIP). However, error-free positive apparent phase or negative apparent chargeability (i.e. negative IP effects) can occur in IP measurements over heterogeneous media. While negative IP effects in TDIP data sets have been discussed, no studies have addressed this topic in detail for FDIP measurements. We describe theory and numerical modelling to explain the origin of negative IP effects in FDIP measurements. A positive apparent phase may occur when a relatively high polarizability feature falls into negative sensitivity zones of complex resistivity measurements. The polarity of the apparent phase is determined by the distribution of subsurface intrinsic phase and resistivity, with the resistivity impacting the apparent phase polarity via its control on the sensitivity distribution. A physical explanation for the occurrence of positive apparent phase data is provided by an electric circuit model representing a four-electrode measurement. We also show that the apparent phase polarity will be frequency dependent when resistivity changes significantly with frequency (i.e. in the presence of significant IP effects). Consequently, negative IP effects manifest themselves in the shape of apparent phase spectra recorded with multifrequency (spectral IP) data sets. Our results imply that positive apparent phase measurements should be anticipated and should be retained during inversion and interpretation of single frequency and spectral IP data sets.
Groundwater‐surface water exchange within the hyporheic zone is widely recognized as a key mechanism controlling the fate of nutrients within catchments. In gaining river systems, groundwater‐surface ...water interactions are constrained by upwelling groundwater but there is increasing evidence that a rapid rise in river stage during storm events can result in a temporary reversal of vertical hydraulic gradients, leading to surface water infiltration into the subsurface and supply of surface‐borne reactive solutes to this biogeochemically active interface. At a UK study site, using logged hydraulic heads in the surface water, riverbed, and riverbanks and logged electrical conductivity at multiple depths in the riverbed we show that storm events can lead to a temporary reversal of vertical hydraulic gradient with mixing evident up to 30 cm beneath the riverbed. Cross‐channel variability is evident, with the center of the channel consistently having shorter reversals of hydraulic gradient, compared to the channel margins. The direction of shallow subsurface riverbank flow at the site is also reactive to storm events, temporarily aligning with the surface flow direction and then reverting back to preevent conditions. Such a transition of flow paths during events is also likely to lead to expansion of lateral hyporheic exchange. This study provides evidence that storm events can be a key driver of enhanced hyporheic exchange in gaining river systems, which may support nutrient reactions beyond the duration of event‐driven change. Our observations demonstrate the dynamic nature of the hyporheic zone, which should be considered when evaluating its biogeochemical function.
Key Points:
Hyporheic zone is temporarily expanded due to storm events
Short‐lived downwelling during events may fuel biogeochemical processes
Electrical conductivity is highlighted as a tracer of hyporheic travel times