Joint parameter estimation from magnetic resonance and vertical electric soundings using a multi‐objective genetic algorithm

Magnetic resonance sounding (MRS) has increasingly become an important method in hydrogeophysics because it allows for estimations of essential hydraulic properties such as porosity and hydraulic conductivity. A resistivity model is required for magnetic resonance sounding modelling and inversion. Therefore, joint interpretation or inversion is favourable to reduce the ambiguities that arise in separate magnetic resonance sounding and vertical electrical sounding (VES) inversions. A new method is suggested for the joint inversion of magnetic resonance sounding and vertical electrical sounding data. A one‐dimensional blocky model with varying layer thicknesses is used for the subsurface discretization. Instead of conventional derivative‐based inversion schemes that are strongly dependent on initial models, a global multi‐objective optimization scheme (a genetic algorithm [GA] in this case) is preferred to examine a set of possible solutions in a predefined search space. Multi‐objective joint optimization avoids the domination of one objective over the other without applying a weighting scheme. The outcome is a group of non‐dominated optimal solutions referred to as the Pareto‐optimal set. Tests conducted using synthetic data show that the multi‐objective joint optimization approximates the joint model parameters within the experimental error level and illustrates the range of trade‐off solutions, which is useful for understanding the consistency and conflicts between two models and objectives. Overall, the Levenberg‐Marquardt inversion of field data measured during a survey on a North Sea island presents similar solutions. However, the multi‐objective genetic algorithm method presents an efficient method for exploring the search space by producing a set of non‐dominated solutions. Borehole data were used to provide a verification of the inversion outcomes and indicate that the suggested genetic algorithm method is complementary for derivative‐based inversions.     

Appraisal of a new 1D weighted joint inversion of ground based and helicopter‐borne electromagnetic data

In order to couple spatial data from frequency‐domain helicopter‐borne electromagnetics with electromagnetic measurements from ground geophysics (transient electromagnetics and radiomagnetotellurics), a common 1D weighted joint inversion algorithm for helicopter‐borne electromagnetics, transient electromagnetics and radiomagnetotellurics data has been developed. The depth of investigation of helicopter‐borne electromagnetics data is rather limited compared to time‐domain electromagnetics sounding methods on the ground. In order to improve the accuracy of model parameters of shallow depth as well as of greater depth, the helicopter‐borne electromagnetics, transient electromagnetics, and radiomagnetotellurics measurements can be combined by using a joint inversion methodology. The 1D joint inversion algorithm is tested for synthetic data of helicopter‐borne electromagnetics, transient electromagnetics and radiomagnetotellurics. The proposed concept of the joint inversion takes advantage of each method, thus providing the capability to resolve near surface (radiomagnetotellurics) and deeper electrical conductivity structures (transient electromagnetics) in combination with valuable spatial information (helicopter‐borne electromagnetics). Furthermore, the joint inversion has been applied on the field data (helicopter‐borne electromagnetics and transient electromagnetics) measured in the Cuxhaven area, Germany. In order to avoid the lessening of the resolution capacities of one data type, and thus balancing the use of inherent and ideally complementary information content, a parameter reweighting scheme that is based on the exploration depth ranges of the specific methods is proposed. A comparison of the conventional joint inversion algorithm, proposed by Jupp and Vozoff ( 1975 ), and of the newly developed algorithm is presented. The new algorithm employs the weighting on different model parameters differently. It is inferred from the synthetic and field data examples that the weighted joint inversion is more successful in explaining the subsurface than the classical joint inversion approach. In addition to this, the data fittings in weighted joint inversion are also improved.     

In‐Situ Sonication for Enhanced Recovery of Aquifer Microbial Communities

Sampling methods for characterization of microbial communities in aquifers should target both suspended and attached microorganisms (biofilms). We investigated the effectiveness and reproducibility of low‐frequency (200 Hz) sonication pulses on improving extraction efficiency and quality of microorganisms from a petroleum‐contaminated aquifer in Studen (Switzerland). Sonication pulses at different power levels (0.65, 0.9, and 1.1 kW) were applied to three different groundwater monitoring wells. Groundwater samples extracted after each pulse were compared with background groundwater samples for cell and adenosine tri‐phosphate concentration. Turbidity values were obtained to assess the release of sediment fines and associated microorganisms. The bacterial community in extracted groundwater samples was analyzed by terminal‐restriction‐fragment‐length polymorphism and compared with communities obtained from background groundwater samples and from sediment cores. Sonication enhanced the extraction efficiency up to 13‐fold, with most of the biomass being associated with the sediment fines extracted with groundwater. Consecutive pulses at constant power were decreasingly effective, while pulses with higher power yielded the best results both in terms of extraction efficiency and quality. Our results indicate that low‐frequency sonication may be a viable and cost‐effective tool to improve the extraction of microorganisms from aquifers, taking advantage of existing groundwater monitoring wells.     

Can GPS-Derived Surface Loading Bridge a GRACE Mission Gap?

We investigated two ‘gap-filler’ methods based on GPS-derived low-degree surface loading variations (GPS-I and GPS-C) and a more simple method (REF-S) which extends a seasonal harmonic variation into the expected Gravity Recovery and Climate Experiment (GRACE) mission gap. We simulated two mission gaps in a reference solution (REF), which is derived from a joint inversion of GRACE (RL05) data, GPS-derived surface loading and simulated ocean bottom pressure. The GPS-I and GPS-C methods both have a new type of constraint applied to mitigate the lack of GPS station network coverage over the ocean. To obtain the GPS-C solution, the GPS-I method is adjusted such that it fits the reference solution better in a 1.5 year overlapping period outside of the gap. As can be expected, the GPS-I and GPS-C solutions contain larger errors compared to the reference solution, which is heavily constrained by GRACE. Within the simulated gaps, the GPS-C solution generally fits the reference solution better compared to the GPS-I method, both in terms of spherical harmonic loading coefficients and in terms of selected basin-averaged hydrological mass variations. Depending on the basin, the RMS-error of the water storage variations (scaled for leakage effects) ranges between 1.6 cm (Yukon) and 15.3 cm (Orinoco). In terms of noise level, the seasonal gap-filler method (REF-S) even outperforms the GPS-I and GPS-C methods, which are still affected by spatial aliasing problems. However, it must be noted that the REF-S method cannot be used beyond the study of simple harmonic seasonal variations.     

Frictional Properties and Microstructure of Calcite-Rich Fault Gouges Sheared at Sub-Seismic Sliding Velocities

We report an experimental and microstructural study of the frictional properties of simulated fault gouges prepared from natural limestone (96 % CaCO₃) and pure calcite. Our experiments consisted of direct shear tests performed, under dry and wet conditions, at an effective normal stress of 50 MPa, at 18–150 °C and sliding velocities of 0.1–10 μm/s. Wet experiments used a pore water pressure of 10 MPa. Wet gouges typically showed a lower steady-state frictional strength (μ = 0.6) than dry gouges (μ = 0.7–0.8), particularly in the case of the pure calcite samples. All runs showed a transition from stable velocity strengthening to (potentially) unstable velocity weakening slip above 80–100 °C. All recovered samples showed patchy, mirror-like surfaces marking boundary shear planes. Optical study of sections cut normal to the shear plane and parallel to the shear direction showed both boundary and inclined shear bands, characterized by extreme grain comminution and a crystallographic preferred orientation. Cross-sections of boundary shears, cut normal to the shear direction using focused ion beam—SEM, from pure calcite gouges sheared at 18 and 150 °C, revealed dense arrays of rounded, ~0.3 μm-sized particles in the shear band core. Transmission electron microscopy showed that these particles consist of 5–20 nm sized calcite nanocrystals. All samples showed evidence for cataclasis and crystal plasticity. Comparing our results with previous models for gouge friction, we suggest that frictional behaviour was controlled by competition between crystal plastic and granular flow processes active in the shear bands, with water facilitating pressure solution, subcritical cracking and intergranular lubrication. Our data have important implications for the depth of the seismogenic zone in tectonically active limestone terrains. Contrary to recent claims, our data also demonstrate that nanocrystalline mirror-like slip surfaces in calcite(-rich) faults are not necessarily indicative of seismic slip rates.     

Regional gravity field modeling using free-positioned point masses

A two-step free-positioned point mass method is used for regional gravity field modeling together with the remove-compute-restore (RCR) technique. The Quasi-Newton algorithm (L-BFGS-B) is implemented to solve the nonlinear problem with bound constraints in the first step, while in the second step the magnitudes of the point masses are re-adjusted with known positions in the least-squares sense. In order to reach a good representation of the gravity field, a number of parameter sets have to be defined carefully before the computations. The effects of four important parameter sets (depth limits, number of point masses, original/reduced basis functions and optimization directions) are investigated for regional gravity field modeling based on two numerical test cases with synthetic and real data. The results show that the selection of the initial depth and depth limits is of most importance. The number of point masses for obtaining a good fit is affected by the data distribution, while a dependency on the data variability (signal variation) is negligible. Long-wavelength errors in the predicted height anomalies can be reduced significantly by using reduced basis functions, and the radial-direction optimization proves to be stable and reliable for regular and irregular data scenarios. If the parameter sets are defined properly, the solutions are similar to the ones computed by least-squares collocation (LSC), but require fewer unknowns than LSC.     

Impact of channel bar form and grain size variability on estimated stranding risk of juvenile brown trout during hydropeaking

Hydropeaking leads to artificial fluctuations in discharge and corresponding water levels with pronounced dewatering areas between base and peak flow along gravel bars and channel banks. In the present study, 16 hydropeaking reaches in Austria were investigated to assess possible differences in the estimated stranding risk for young of the year brown trout according to different gravel bar types and differences in microtopography roughness. Based on hydrodynamic‐numerical modelling, a predictive habitat modelling approach was implemented in the study design. Accompanied by grain size sampling along the various channel bars, a conceptual stranding risk model (SRM) was developed. The results showed that a high variability in estimated stranding risk exists for the tested sites considering discharge ratios of 1:3, 1:5 and 1:10. With respect to the discussion of establishing legal thresholds for ramping ratios in discharge, it was documented that, exemplarily, a discharge ratio base flow/peak flow of 1:5 (winter base flow conditions) could cause minor differences in the spatial extent of dewatering areas and the related estimated stranding risk for juvenile brown trout compared to a ratio of 1:2 for summer base flow conditions. Microtopographic roughness was addressed due to sampling and analysis of grain size distributions. Statistical testing of grain size distributions revealed significant differences between the surface material compositions of the investigated gravel bars. Those differences are evident, particularly for the coarser fraction (d₉₀), which is important as cover for young of the year brown trout. These aspects of grain size in habitat use and hydraulics have been addressed in the conceptual SRM. The results showed that point bar morphology, in particular, was less sensitive to the risk of stranding compared to, for example, alternating gravel bars. Considering the multiple pressures for alpine rivers, the improvement of structural features due to bar formation and related self‐forming processes is discussed as a possible alternative for future mitigation measures to reduce the negative impacts of hydropeaking. Copyright © 2014 John Wiley & Sons, Ltd.     

The use of mycorrhiza for eco‐engineering measures in steep alpine environments: effects on soil aggregate formation and fine‐root development

Steep erosion‐prone and vegetation‐free slopes are widespread in alpine areas and are often discussed since they have a high socio‐economic damage potential. We present an eco‐engineering approach to test whether a mycorrhizal inoculum improves the establishment of hedge brush layers and in turn soil structural stability on a steep, coarse‐grained vegetation‐free slope in the eastern Swiss Alps. We established (i) mycorrhizal and (ii) non‐mycorrhizal treated eco‐engineered research plots on a field experimental scale, covering a total area of approximately 1000 m² on an east‐northeast (ENE) exposed slope, where many environmental parameters can be regarded as homogeneous. After a full vegetation period, we quantified soil aggregate stability, the formation of water stable aggregates and the fine‐root development. Our results illustrate that the establishment of brush layers without mycorrhizal inoculum increased aggregate stability significantly. Against our expectation and glasshouse experiments, the addition of mycorrhizal inoculum did not have a statistically significant effect after one vegetation period although it tended to increase aggregate stability. Analogously, root length density (RLD) tended to be higher at the non‐mycorrhizal treated site. Aggregate stability was significantly correlated with RLD. Studies on a bigger field experimental scale are inevitable, complement glasshouse studies and lead to a better understanding for a successful application of sustainable eco‐engineering measures in alpine environments. Based on our results and considering the fact that the response time in natural ecosystems may be slower than in laboratory approaches, we conclude that long‐term field studies are necessary to validate results gained through laboratory experiments. Copyright © 2014 John Wiley & Sons, Ltd.     

Influence of river channel morphology and bank characteristics on water surface boundary delineation using high‐resolution passive remote sensing and template matching

Accurate mapping of water surface boundaries in rivers is an important step for monitoring water stages, estimating discharge, flood extent, and geomorphic response to changing hydrologic conditions, and assessing riverine habitat. Nonetheless, it is a challenging task in spatially and spectrally heterogeneous river environments, commonly characterized by high spatiotemporal variations in morphology, bed material, and bank cover. In this study, we investigate the influence of channel morphology and bank characteristics on the delineation of water surface boundaries in rivers using high spatial resolution passive remote sensing and a template‐matching (object‐based) algorithm, and compare its efficacy with that of Support Vector Machine (SVM) (pixel‐based) algorithm. We perform a detailed quantitative evaluation of boundary‐delineation accuracy using spatially explicit error maps in tandem with the spatial maps of geomorphic and bank classes. Results show that template matching is more successful than SVM in delineating water surface boundaries in river sections with spatially challenging geomorphic landforms (e.g. sediment bar structures, partially submerged sediment deposits) and shallow water conditions. However, overall delineation accuracy by SVM is higher than that of template matching (without iterative hierarchical learning). Vegetation and water indices, especially when combined with texture information, improve the accuracy of template matching, for example, in river sections with overhanging trees and shadows – the two most problematic conditions in water surface boundary delineation. By identifying the influence of channel morphology and bank characteristics on water surface boundary mapping, this study helps determine river sections with higher uncertainty in delineation. In turn, the most suitable methods and data sets can be selectively utilized to improve geomorphic/hydraulic characterization. The methodology developed here can also be applied to similar studies on other geomorphic landforms including floodplains, wetlands, lakes, and coastlines. Copyright © 2014 John Wiley & Sons, Ltd.