A Bayesian trans-dimensional approach for the fusion of multiple geophysical datasets

We propose a Bayesian fusion approach to integrate multiple geophysical datasets with different coverage and sensitivity. The fusion strategy is based on the capability of various geophysical methods to provide enough resolution to identify either subsurface material parameters or subsurface structure, or both. We focus on electrical resistivity as the target material parameter and electrical resistivity tomography (ERT), electromagnetic induction (EMI), and ground penetrating radar (GPR) as the set of geophysical methods. However, extending the approach to different sets of geophysical parameters and methods is straightforward. Different geophysical datasets are entered into a trans-dimensional Markov chain Monte Carlo (McMC) search-based joint inversion algorithm. The trans-dimensional property of the McMC algorithm allows dynamic parameterisation of the model space, which in turn helps to avoid bias of the post-inversion results towards a particular model. Given that we are attempting to develop an approach that has practical potential, we discretize the subsurface into an array of one-dimensional earth-models. Accordingly, the ERT data that are collected by using two-dimensional acquisition geometry are re-casted to a set of equivalent vertical electric soundings. Different data are inverted either individually or jointly to estimate one-dimensional subsurface models at discrete locations. We use Shannon's information measure to quantify the information obtained from the inversion of different combinations of geophysical datasets. Information from multiple methods is brought together via introducing joint likelihood function and/or constraining the prior information. A Bayesian maximum entropy approach is used for spatial fusion of spatially dispersed estimated one-dimensional models and mapping of the target parameter. We illustrate the approach with a synthetic dataset and then apply it to a field dataset. We show that the proposed fusion strategy is successful not only in enhancing the subsurface information but also as a survey design tool to identify the appropriate combination of the geophysical tools and show whether application of an individual method for further investigation of a specific site is beneficial.     

How Does Subsurface Characterization Affect Simulations of Hyporheic Exchange?

We investigated the role of increasingly well‐constrained geologic structures in the subsurface (i.e., subsurface architecture) in predicting streambed flux and hyporheic residence time distribution (RTD) for a headwater stream. Five subsurface realizations with increasingly resolved lithological boundaries were simulated in which model geometries were based on increasing information about flow and transport using soil and geologic maps, surface observations, probing to depth to refusal, seismic refraction, electrical resistivity (ER) imaging of subsurface architecture, and time‐lapse ER imaging during a solute tracer study. Particle tracking was used to generate RTDs for each model run. We demonstrate how improved characterization of complex lithological boundaries and calibration of porosity and hydraulic conductivity affect model prediction of hyporheic flow and transport. Models using hydraulic conductivity calibrated using transient ER data yield estimates of streambed flux that are three orders of magnitude larger than uncalibrated models using estimated values for hydraulic conductivity based on values published for nearby hillslopes (10^?4 vs. 10^?7 m²/s, respectively). Median residence times for uncalibrated and calibrated models are 10³ and 10⁰ h, respectively. Increasingly well‐resolved subsurface architectures yield wider hyporheic RTDs, indicative of more complex hyporheic flowpath networks and potentially important to biogeochemical cycling. The use of ER imaging to monitor solute tracers informs subsurface structure not apparent from other techniques, and helps to define transport properties of the subsurface (i.e., hydraulic conductivity). Results of this study demonstrate the value of geophysical measurements to more realistically simulate flow and transport along hyporheic flowpaths.     

Integration of ground-penetrating radar and microgravimetric methods to map shallow caves

Ground-penetrating radar (GPR) and microgravimetric surveys have been conducted in the southern Jura mountains of western Switzerland in order to map subsurface karstic features. The study site, La Grande Rolaz cave, is an extensive system in which many portions have been mapped. By using small station spacing and careful processing for the geophysical data, and by modeling these data with topographic information from within the cave, accurate interpretations have been achieved.The constraints on the interpreted geologic models are better when combining the geophysical methods than when using only one of the methods, despite the general limitations of two-dimensional (2D) profiling. For example, microgravimetry can complement GPR methods for accurately delineating a shallow cave section approximately 10×10 m² in size. Conversely, GPR methods can be complementary in determining cavity depths and in verifying the presence of off-line features and numerous areas of small cavities and fractures, which may be difficult to resolve in microgravimetric data.     

Near-surface study of a hot spring site in Fiji

A geophysical investigation of a hot spring system located in Rabulu, Fiji, was carried out from October 2014 to March 2015. The investigation covered a survey area of 6075 m². Self-potential (SP), ground temperature and soil carbon dioxide (CO₂) concentrations were measured and investigated for their distribution characteristics and inter-linkages. Results indicated obvious anomalous zone at the hot spring discharge site. The SP profile analysis highlighted thermal water upwelling zones and elevation-driven subsurface groundwater pathways. Measurement of subsurface temperatures up to 1 m depths revealed increasing temperatures, indicating potentially high thermal gradients in the area. Surface soil CO₂ distributions also agreed with SP and ground temperature results. The overall result of the study demonstrated that synchronised measurements of SP, ground temperature and soil CO₂ can be instrumental in identifying anomalous zones near the hot spring sites. Other parameters such as spring water temperature, discharge rate and energy flux estimates from the spring were calculated and analysed. The high-dense multi-parameter data coverage allowed interpretation of geothermal features at a scale never conducted in Fiji before. The near-surface investigations reported in this study corroborate previously suggested steady geothermal activity in the region, deserving further detailed investigation.     

Investigations regarding Alpine talus slopes using ground‐penetrating radar (GPR) in the Bavarian Alps,Germany

The applicability of ground‐penetrating radar (GPR) for the investigation of loose debris was tested at two sites (Viererkar and Zugspitzplatt). A pulseEKKO 100 GPR system equipped with 25 MHz antennae was utilized. The aim of the investigation was to record the base of the debris layer, and thereby acquire an estimation of the backweathering rates of the adjacent rockwalls. The study areas are situated in the Northern Alps near the German–Austrian border. The sites are characterized by steep limestone rockwalls and extensive talus accumulations. A total of six profiles was surveyed. The method is suitable and effective for a quick survey in this dry, high‐ohmic substrate. The GPR system was able to deliver information about the subsurface stratigraphy to c. 70 m depth. The boundary line to the bedrock was discovered – depending upon the profile surveyed ?5 to 25 m below the surface. The base of the debris material sometimes shows no distinct reflection. Buried features (V‐shaped furrows, zones overdeepened by ice action, geological structures) could be detected. Arched structures well below the talus–bedrock interface can be interpreted as drainage systems in the karstic bedrock. A thick scree layer of Late Glacial age was separated from a thinner layer on the talus surface, which was related to the Holocene. The backweathering rates were fixed by a calculation of talus volume to c. 100 mm/10³ a during the Holocene (Viererkar) and 150–300 mm/10³ a (Zugspitzplatt). The detrital formation in north‐exposed sites is twice as intense as in south‐exposed sites. These results match the rates of recent rockfall in the same area of investigation. The calculated backweathering for the late glacial period is 150–730 mm/10³ a. The magnitude of the calculated rockwall retreat lies well within the range of previous measurements. The discrepancy between some weathering rates highlights the fact that recent and past relief formation must be differentiated. Otherwise recent removal rates may be overestimated. Copyright © 2001 John Wiley & Sons, Ltd.     

New Ways of Using Well Construction Reports for Hydrostratigraphic Analyses

Hydrostratigraphic models require subsurface data. Those data can be gleaned from Well Construction Reports (WCRs) if the reports are properly reviewed and processed. In the United States, WCRs or similar documents are generally required to be submitted to a state agency. The quality of those reports varies considerably and depends on the water well contractor's level of geologic knowledge. Nevertheless, these types of reports are usually the largest available subsurface data set. We present two examples from Wisconsin, each of which had approximately 2000 WCRs, to demonstrate a variety of new applications of WCRs including obtaining information about the depositional environment; selecting sites to perform surface geophysical measurements; and creating three‐dimensional hydrostratigraphic models for groundwater modeling.     

ON THE THEORY OF EVOLUTION OF RIVER NETS

Abstract

An investigation is made of the growth of river basins based on the two current models of the structure of such basins: the cyclic model and the random-graph model. Based on the assumption that the networks grow allometrically, the structure of the network that evolves has been calculated based upon the above two models. Nature seems to favor the random-graph model.     

Characterizing hydrological processes in a semiarid rangeland watershed: A hydrogeophysical approach

The complex ecohydrological processes of rangelands can be studied through the framework of ecological sites (ESs) or hillslope‐scale soil–vegetation complexes. High‐quality hydrologic field investigations are needed to quantitatively link ES characteristics to hydrologic function. Geophysical tools are useful in this context because they provide valuable information about the subsurface at appropriate spatial scales. We conducted 20 field experiments in which we deployed time‐lapse electrical resistivity tomography (ERT), variable intensity rainfall simulation, ground‐penetrating radar (GPR), and seismic refraction, on hillslope plots at five different ESs within the Upper Crow Creek Watershed in south‐east Wyoming. Surface runoff was measured using a precalibrated flume. Infiltration data from the rainfall simulations, coupled with site‐specific resistivity–water content relationships and ERT datasets, were used to spatially and temporally track the progression of the wetting front. First‐order constraints on subsurface structure were made at each ES using the geophysical methods. Sites ranged from infiltrating 100% of applied rainfall to infiltrating less than 60%. Analysis of covariance results indicated significant differences in the rate of wetting front progression, ranging from 0.346 m min^?1/2 for sites with a subsurface dominated by saprolitic material to 0.156 m min^?1/2 for sites with a well‐developed soil profile. There was broad agreement in subsurface structure between the geophysical methods with GPR typically providing the most detail. Joint interpretation of the geophysics showed that subsurface features such as soil layer thickness and the location of subsurface obstructions such as granite corestones and material boundaries had a large effect on the rate of infiltration and subsurface flow processes. These features identified through the geophysics varied significantly by ES. By linking surface hydrologic information from the rainfall simulations with subsurface information provided by the geophysics, we can characterize the ES‐specific hydrologic response. Both surface and subsurface flow processes differed among sites and are directly linked to measured characteristics.     

Site Response from Ambient Vibrations in the Towns of Lod and Ramle (Israel) and Earthquake Hazard Assessment

Throughout their long history, the towns of Lod and Ramle have been severely affected by strong earthquakes. The last destructive earthquake occurred on July 11, 1927 and caused the destruction of large parts of these cities, reaching a seismic intensity of VIII–IX on the MSK scale. Such a high intensity from a relatively distant earthquake (about 70 km) of magnitude 6.2 is likely to be the result of local site effects of the sedimentary layers that may have significantly enhanced earthquake ground motions. This study is focused on estimating the seismic hazard to Lod and Ramle by implementing a three-step process (1) detailed mapping of the characteristics of the H/V spectral ratios from ambient noise, (2) incorporating geological information and well data to construct subsurface models for different sites within the investigated area and (3) estimating the seismic hazard in terms of uniform hazard site-specific accelerations. The horizontal-to-vertical spectral ratios of ambient noise were used to approximate the fundamental resonance frequencies of the subsurface and their associated amplitudes. About 360 sites in Lod and Ramle were instrumented for varying periods. The soil sites exhibits H/V peak amplitudes ranging from 4 to 6 in the frequency range 0.5–2.5 Hz. These data were used to constrain 1-D subsurface models that were developed using geological data and borehole information. H/V spectral ratio observations were checked against theoretical subsurface transfer functions at locations where borehole information is available farther constraint the range of possible V_s velocities of the different layers and thus, by means of trial an error it was possible to conclude a systematic spatial distribution of the V_s velocity and thickness in the substrata that are also consistent with the spatial distribution of the fundamental resonance frequencies of the soft sediments obtained by means of the H/V spectral ratios, and other geological and geophysical information available at different locations in the study area. The evaluated subsurface models are introduced using the SEEH procedure of Shapira and van~Eck [(1993) Natural Hazards 8, 201–205] to assess Uniform Hazard Site-Specific Acceleration Spectra for different zones within the towns of Lod and Ramle. These evaluations are very important for realistic assessment of the vulnerabilities of all types of existing and newly designed structures and for urban and land use planning.     

Resolution of MRS applied to the characterization of hard-rock aquifers

The performance of the Magnetic Resonance Sounding (MRS) method applied to the investigation of heterogeneous hard-rock aquifers was studied. It was shown using both numerical modeling and field measurements that MRS could be applied to the investigation of the weathered part of hard-rock aquifers when the product of the free water content multiplied by the thickness of the aquifer is >0.2 (for example, 10-m-thick layer with a 2% water content). Using a currently available one-dimensional MRS system, the method allows the characterization of two-dimensional subsurface structures with acceptable accuracy when the size of the subsurface anomaly is equal to or greater than the MRS loop. However, the fractured part of hard-rock aquifers characterized by low effective porosity (<0.5%) cannot be resolved using currently available MRS equipment. It was found that shallow water in the weathered part of the aquifer may screen MRS signals from deeper water-saturated layers, thus further reducing the possibility of investigating deeper fractured aquifers. A field study using the NUMIS(plus) MRS system developed by IRIS Instruments was carried out on an experimental watershed in southern India. A heterogeneous unconfined aquifer in a gneissic formation was successfully localized, and MRS results were confirmed by drilling shortly after the geophysical study. The top of the aquifer revealed by MRS was found to be in a good agreement with observed static water level measurements in boreholes.     

Groundwater potential zones using a combination of geospatial technology and geophysical approach: case study in Dehradun,India

ABSTRACT

A combination of geospatial, geophysical and statistical models using satellite data, the weighted index overlay (WIO) method and two-dimensional electrical resistivity tomography (2D-ERT) is applied to generate the highest potential groundwater area and to further explore the groundwater in Dehradun, India. The results show that of 19.7 km² total basin area, 0.26% falls under the “poor” category as a prospect zone for groundwater, 4.3% is “moderate”, 10.10% “moderately good”, 4.9% “good” and 0.17% “very good”. In addition, the demonstration of the geophysical survey is presented, in which Purkal Youth Society Division (PYSD) site is categorized as a shallow aquifer zone and the Guru Nanak Fifth Centenary School (GNFCS) site is a deeper aquifer zone. Our study emphasizes remote sensing and geographic information system integrated with a geophysical survey to support prospecting the most probable area and confirm the existence of groundwater.     

Mantle rheology: radial and lateral viscosity variations inferred from microphysical creep laws

Inferences on the rheology of the mantle based on theoretical and experimental rate equations for steady state creep are discussed and compared with results from geophysical models. The radial increase of viscosity by one to three orders of magnitude across the mantle, required by inversion of postglacial rebound and geodynamic data, is confirmed by microphysical models based on the estimation of continuous and discontinuous changes of creep parameters with depth. The upper mantle (viscosity 10²⁰–10²¹ Pa s) is likely to show non-Newtonian rheology (power-law creep) for average grain sizes larger than 0.1 mm as an order of magnitude. Given the variability of both grain size and stress conditions, local regions of linear rheology can be present. The rheology of transition zone and lower mantle (viscosity 10²²–10²⁴ Pa s) cannot be definitely resolved at present. Estimation of creep parameters leads to possible nonlinear or mixed rheology, if grain sizes are not lower than 0.1 mm and flow conditions can be approximated by a constant strain rate of about 10⁻¹⁵ s⁻¹. This conclusion can be modified by different flow conditions (e.g. a decrease in strain rate or constant viscous dissipation). Furthermore, experiments on fine-grained garnetites and perovskite analogues have shown that diffusion creep is predominant at laboratory conditions. However, the pressure dependence of creep in these phases is unknown, and therefore direct extrapolation to lower mantle conditions is necessarily speculative. Lateral variations of viscosity, largest in the upper and lowermost mantle (up to 2–4 orders of magnitude) are predicted by models based on lateral temperature anomalies derived from seismic tomographic models.     

大地电磁资料精细处理和二维反演解释技术研究(六)——交错模型的大地电磁二维反演

交错模型指的是由几组走向不同的线性构造（近似二维）在空间（纵向、横向）上组合而成的模型.对于电性结构,交错模型表现为纵向或横向上的构造存在电性主轴方位的变化.从构造维性的角度来看,交错模型是一种由二维模型组合而成的特殊的三维模型,直接对其进行大地电磁二维反演,不易获得可靠的反演结果.本文针对交错模型的特点,提出分频段-分区段反演方案.该方案首先需要借助于阻抗张量成像技术,在频率域确定组成交错模型的各线性构造的电性主轴,然后,针对不同的频段、区段,选择对应电性主轴的数据进行反演,通过初始模型的构建将不同电性主轴方位的反演结果对接起来.本文通过三维理论模型的研究,系统展示了分频段-分区段反演的全过程,归纳得到：在分频段反演时先做低频段反演,在分区段反演时先做二维性更不显著一侧的反演.最后,本文将这一技术用于郯庐断裂带中南段一条实测剖面的反演中,其结果与常规二维反演结果相比较,深部的信息更为丰富,且与其他已有地质、地球物理结果的可对比性更好,表明在构造复杂地区,大地电磁分频段-分区段二维反演具有较高的模型分辨率和可靠性.     

Calibrating a Salt Water Intrusion Model with Time‐Domain Electromagnetic Data

Salt water intrusion models are commonly used to support groundwater resource management in coastal aquifers. Concentration data used for model calibration are often sparse and limited in spatial extent. With airborne and ground‐based electromagnetic surveys, electrical resistivity models can be obtained to provide high‐resolution three‐dimensional models of subsurface resistivity variations that can be related to geology and salt concentrations on a regional scale. Several previous studies have calibrated salt water intrusion models with geophysical data, but are typically limited to the use of the inverted electrical resistivity models without considering the measured geophysical data directly. This induces a number of errors related to inconsistent scales between the geophysical and hydrologic models and the applied regularization constraints in the geophysical inversion. To overcome these errors, we perform a coupled hydrogeophysical inversion (CHI) in which we use a salt water intrusion model to interpret the geophysical data and guide the geophysical inversion. We refer to this methodology as a Coupled Hydrogeophysical Inversion‐State (CHI‐S), in which simulated salt concentrations are transformed to an electrical resistivity model, after which a geophysical forward response is calculated and compared with the measured geophysical data. This approach was applied for a field site in Santa Cruz County, California, where a time‐domain electromagnetic (TDEM) dataset was collected. For this location, a simple two‐dimensional cross‐sectional salt water intrusion model was developed, for which we estimated five uniform aquifer properties, incorporating the porosity that was also part of the employed petrophysical relationship. In addition, one geophysical parameter was estimated. The six parameters could be resolved well by fitting more than 300 apparent resistivities that were comprised by the TDEM dataset. Except for three sounding locations, all the TDEM data could be fitted close to a root‐mean‐square error of 1. Possible explanations for the poor fit of these soundings are the assumption of spatial uniformity, fixed boundary conditions and the neglecting of 3D effects in the groundwater model and the TDEM forward responses.     

FIRST EXPERIENCE AND RESULTS WITH A NEW SEISMIC INTERPRETATION METHOD*

At the EAEG meeting in Venice (May 1969) Linsser presented a new method for the interpretation of seismic measurements (“Transformation of Seismic Data into Geological Models by Digital Template Analysis”). First results were presented during the S.E.G. convention in Calgary (September 1969). Meanwhile some experience concerning the possibility of the application of this method is available. Starting from the fact that the computer is able to store several physical parameters (amplitude, wave-length, direction of the seismic event, coincidence with a template and the beginning and the end of each reflection), the importance of these data for the investigation of special geophysical questions will be discussed. This method has been tested with several examples on the following problems:

• 1 Proving and exact localisation of faults.
• 2 Investigation of velocities of interesting layers.
• 3 Additional information for the geological interpretation of seismic sections.
• 4 Picking of weak reflections in the “Pre-Zechsteinian” in Northern Germany.

Finally the possibilities of further applications and developments are discussed.     

SIMPLE DATA PROCESSING OF TRIPOTENTIAL APPARENT RESISTIVITY MEASUREMENTS AS AN AID TO THE INTERPRETATION OF SUBSURFACE STRUCTURE*

The advantages of the Wenner tripotential method (Carpenter 1955) for apparent resistivity profiling are described and two new data processing techniques introduced as an aid to the interpretation of apparent resistivity sections (pseudo-sections). These techniques were developed from model data computed using a two-dimensional finite difference method. Oscillatory components present in anomalies on tripotential profiles and related to electrode spacing are shown to be effectively removed by linear filtering that also simplifies their form and aids recognition. Furthermore, it is shown that the ratio of the beta- and gamma-apparent resistivities is a good indicator of resistivity variation, and is particularly sensitive to lateral change. Model data indicates that, over a wide range of conditions, enough subsurface information can be obtained by inspection of tripotential resistivity and ratio profiles, and from space sections to make possible a useful—and sometimes semi-quantitative—interpretation. A rationale for the general interpretation of tripotential data is developed. Field data are described from an area of weathered granite basement in Nigeria. A model of the subsurface is developed using parameters derived from the processed observations. The observed and calculated apparent resistivity space sections are very similar.     

Development and application of a distributed catchment-scale hydrological model for the River Ythan,NE Scotland

Mathematical models are being used to develop a decision support system for integrated management of the Ythan catchment in NE Scotland. One component of this has involved the development of a distributed catchment-scale hydrological model. The model is based on subsurface flow routing and calculates the contribution to stream flow from each 50 m×50 m cell in the 548 km² catchment. It uses two topographic parameters, slope and distance to stream following the main line of flow, and five physical parameters. The topographic analysis and distributed flow accumulation are performed by linking the single cell model with a geographic information system. Preliminary results from a three-year simulation of daily flows indicate that the model successfully predicts the main characteristics of the catchment flow. © 1998 John Wiley & Sons, Ltd.     

Soil loss rates due to piping erosion

Compared with surface soil erosion by water, subsurface erosion (piping) is generally less studied and harder to quantify. However, wherever piping occurs, it is often a significant or even the main sediment source. In this study, the significance of soil loss due to piping is demonstrated through an estimation of soil volume lost from pipes and pipe collapses (n = 560) in 137 parcels under pasture on loess‐derived soils in a temperate humid climate (Belgium). Assuming a period of 5 to 10 years for pipe collapse to occur, mean soil loss rates of 2.3 and 4.6 t ha^?1 yr^?1 are obtained, which are at least one order of magnitude higher than surface erosion rates (0.01–0.29 t ha^?1 yr^?1) by sheet and rill erosion under a similar land use. The results obtained for the study area in the Flemish Ardennes correspond well to other measurements in temperate environments; they are, however, considerably smaller than soil loss rates due to subsurface erosion in semi‐arid environments. Although local slope gradient and drainage area largely control the location of collapsed pipes in the study area, these topographic parameters do not explain differences in eroded volumes by piping. Hence, incorporation of subsurface erosion in erosion models is not straightforward. Copyright © 2011 John Wiley & Sons, Ltd.     

Imaging DOAS for volcanological applications

The simultaneous quantitative determination of two-dimensional bromine monoxide (BrO) and sulphur dioxide (SO₂) distributions in volcanic gas plumes is described. Measurements at the fumarolic field on the island Vulcano (autumn 2004) and in the plume of Mt. Etna volcano (spring 2005) were carried out with an Imaging DOAS instrument. The SO₂ fluxes of several fumaroles were estimated from two-dimensional distributions of SO₂. Additionally, the first two-dimensional distributions of BrO within a volcanic plume were successfully retrieved. Slant column densities of up to 2.6 × 10¹⁴ molecules per square centimetre were detected in the plume of Mt. Etna. The investigation of the BrO/SO₂ ratio, calculated from the two-dimensional distributions of SO₂ and BrO, shows an increase from the centre to the edge of the volcanic plume. These results have significance for the involvement of ozone during BrO formation processes in volcanic emissions.     

Groundwater connections under a barrier beach: A case study in the Venice Lagoon

The flow of groundwater beneath barrier islands has been cited as a possible pathway for salt water and chemical exchange between a protected embayment and the open sea. Evidence is presented that identifies an exchange of groundwater through a highly permeable paleoinlet along the barrier beach of Cavallino, which separates the northern Venice Lagoon from the Adriatic Sea. We utilized both point measurements of submarine groundwater discharge (SGD) and a geophysical investigation of the subsurface resistivity to analyze the movement of saline groundwater. Discharge of groundwater and associated nutrients, was higher at the site of a former inlet than at a similar site along the barrier and modulated by the difference in tidal water level between the lagoon and Adriatic Sea. If the measured conditions are typical, storm surge barriers could potentially result in a saline groundwater flow of up to 1.5×10⁶ m³ d⁻¹ into the lagoon.