Refine
Has Fulltext
- no (49)
Year of publication
Document Type
- Article (49) (remove)
Language
- English (49)
Is part of the Bibliography
- yes (49)
Keywords
- Ground-penetrating radar (4)
- Electromagnetics (2)
- Global inversion (2)
- Ground penetrating radar (2)
- 3D modeling (1)
- Applied geophysics (1)
- Aquifer (1)
- Austrian Alps (1)
- Brilliant blue (1)
- Civil engineering (1)
Institute
According to Dooge (1986) intermediate-scale catchments are systems of organized complexity, being too organized and yet too small to be characterized on a statistical/conceptual basis, but too large and too heterogeneous to be characterized in a deterministic manner. A key requirement for building structurally adequate models precisely for this intermediate scale is a better understanding of how different forms of spatial organization affect storage and release of water and energy. Here, we propose that a combination of the concept of hydrological response units (HRUs) and thermodynamics offers several helpful and partly novel perspectives for gaining this improved understanding. Our key idea is to define functional similarity based on similarity of the terrestrial controls of gradients and resistance terms controlling the land surface energy balance, rainfall runoff transformation, and groundwater storage and release. This might imply that functional similarity with respect to these specific forms of water release emerges at different scales, namely the small field scale, the hillslope, and the catchment scale. We thus propose three different types of "functional units" - specialized HRUs, so to speak - which behave similarly with respect to one specific form of water release and with a characteristic extent equal to one of those three scale levels. We furthermore discuss an experimental strategy based on exemplary learning and replicate experiments to identify and delineate these functional units, and as a promising strategy for characterizing the interplay and organization of water and energy fluxes across scales. We believe the thermodynamic perspective to be well suited to unmask equifinality as inherent in the equations governing water, momentum, and energy fluxes: this is because several combinations of gradients and resistance terms yield the same mass or energy flux and the terrestrial controls of gradients and resistance terms are largely independent. We propose that structurally adequate models at this scale should consequently disentangle driving gradients and resistance terms, because this optionally allow sequifinality to be partly reduced by including available observations, e. g., on driving gradients. Most importantly, the thermodynamic perspective yields an energy-centered perspective on rainfall-runoff transformation and evapotranspiration, including fundamental limits for energy fluxes associated with these processes. This might additionally reduce equifinality and opens up opportunities for testing thermodynamic optimality principles within independent predictions of rainfall-runoff or land surface energy exchange. This is pivotal to finding out whether or not spatial organization in catchments is in accordance with a fundamental organizing principle.
In an attempt to map the shallow geometry of the Maleme Fault Zone (North Island, New Zealand) and estimate vertical displacements of selected fault strands, we have collected 2D and 3D georadar data using 100 MHz antennae. The 2D data consisted of three parallel georadar lines recorded perpendicular to the axis of the well-defined graben of the Maleme Fault Zone. These similar to 160 in long lines, which were 7.5 m apart, crossed several fault strands on either side of the graben axis. The processed georadar sections revealed two prominent parallel reflections that originated from the boundaries of Late Pleistocene lacustrine and tephra deposits. Distinct vertical offsets of these reflections allowed us to estimate displacernents at individual fault strands across the entire inner graben. The total displacements represented by these offsets was similar to 10-20% greater than that inferred from geomorphological studies, thus demonstrating the limitations of surface observations for determining cumulative fault movements. The 3D georadar data set, recorded across an area of similar to 70x similar to 20 in to one side of the graben axis, provided key details on individual fault strands. For the 3D visualization of fault-related structures, various spatial attribute analyses based on the cosine of the instantaneous phase proved to be useful
Modern mobile devices (i.e. smartphones and tablet computers) are widespread, everyday tools, which are equipped with a variety of sensors including three-axis magnetometers. Here, we investigate the feasibility and the potential of using such mobile devices to mimic geophysical experiments in the classroom in a table-top setup. We focus on magnetic surveying and present a basic setup of a table-top experiment for collecting three-component magnetic data across well-defined source bodies and structures. Our results demonstrate that the quality of the recorded data is sufficient to address a number of important basic concepts in the magnetic method. The shown examples cover the analysis of magnetic data recorded across different kinds of dipole sources, thus illustrating the complexity of magnetic anomalies. In addition, we analyze the horizontal resolution capabilities using a pair of dipole sources placed at different horizontal distances to each other. Furthermore, we demonstrate that magnetic data recorded with a mobile device can even be used to introduce filtering, transformation, and inversion approaches as they are typically used when processing magnetic data sets recorded for real-world field applications. Thus, we conclude that such table-top experiments represent an easy-to-implement experimental procedure (as student exercise or classroom demonstration) and can provide first hands-on experience in the basic principles of magnetic surveying including the fundamentals of data acquisition, analysis and processing, as well as data evaluation and interpretation.
Crosshole traveltime tomography using particle swarm optimization a near-surface field example
(2012)
Particle swarm optimization (PSO) is a relatively new global optimization approach inspired by the social behavior of bird flocking and fish schooling. Although this approach has proven to provide excellent convergence rates in different optimization problems, it has seldom been applied to inverse geophysical problems. Until today, published geophysical applications mainly focus on finding an optimum solution for simple, 1D inverse problems. We have applied PSO-based optimization strategies to reconstruct 2D P-wave velocity fields from crosshole traveltime data sets. Our inversion strategy also includes generating and analyzing a representative ensemble of acceptable models, which allows us to appraise uncertainty and nonuniqueness issues. The potential of our strategy was tested on field data collected at a well-constrained test site in Horstwalde, Germany. At this field site, the shallow subsurface mainly consists of sand- and gravel-dominated glaciofluvial sediments, which, as known from several boreholes and other geophysical experiments, exhibit some well-defined layering at the scale of our crosshole seismic data. Thus, we have implemented a flexible, layer-based model parameterization, which, compared with standard cell-based parameterizations, allows for significantly reducing the number of unknown model parameters and for efficiently implementing a priori model constraints. Comparing the 2D velocity fields resulting from our PSO strategy to independent borehole and direct-push data illustrated the benefits of choosing an efficient global optimization approach. These include a straightforward and understandable appraisal of nonuniqueness issues as well as the possibility of an improved and also more objective interpretation.
Near-surface geophysical techniques are extensively used in a variety of engineering, environmental, geologic, and hydrologic applications. While many of these applications ask for detailed, quantitative models of selected material properties, geophysical data are increasingly used to estimate such target properties. Typically, this estimation procedure relies on a two-step workflow including (1) the inversion of geophysical data and (2) the petrophysical translation of the inverted parameter models into the target properties. Standard deterministic implementations of such a quantitative interpretation result in a single best-estimate model, often without considering and propagating the uncertainties related to the two steps. We address this problem by using a rather novel, particle-swarm-based global joint strategy for data inversion and by implementing Monte Carlo procedures for petrophysical property estimation. We apply our proposed workflow to crosshole ground-penetrating radar, P-, and S-wave data sets collected at a well-constrained test site for a detailed geotechnical characterization of unconsolidated sands. For joint traveltime inversion, the chosen global approach results in ensembles of acceptable velocity models, which are analyzed to appraise inversion-related uncertainties. Subsequently, the entire ensembles of inverted velocity models are considered to estimate selected petrophysical properties including porosity, bulk density, and elastic moduli via well-established petrophysical relations implemented in a Monte Carlo framework. Our results illustrate the potential benefit of such an advanced interpretation strategy; i.e., the proposed workflow allows to study how uncertainties propagate into the finally estimated property models, while concurrently we are able to study the impact of uncertainties in the used petrophysical relations (e.g., the influence of uncertain, user-specified parameters). We conclude that such statistical approaches for the quantitative interpretation of geophysical data can be easily extended and adapted to other applications and geophysical methods and might be an important step toward increasing the popularity and acceptance of geophysical tools in engineering practice.
Vertical radar profiling (VRP) is a single-borehole geophysical technique, in which the receiver antenna is located within a borehole and the transmitter antenna is placed at one or various offsets from the borehole. Today, VRP surveying is primarily used to derive 1D velocity models by inverting the arrival times of direct waves. Using field data collected at a well-constrained test site in Germany, we evaluated a VRP workflow relying on the analysis of direct-arrival traveltimes and amplitudes as well as on imaging reflection events. To invert our VRP traveltime data, we used a global inversion strategy resulting in an ensemble of acceptable velocity models, and thus, it allowed us to appraise uncertainty issues in the estimated velocities as well as in porosity models derived via petrophysical translations. In addition to traveltime inversion, the analysis of direct-wave amplitudes and reflection events provided further valuable information regarding subsurface properties and architecture. The used VRP amplitude preprocessing and inversion procedures were adapted from raybased crosshole ground-penetrating radar (GPR) attenuation tomography and resulted in an attenuation model, which can be used to estimate variations in electrical resistivity. Our VRP reflection imaging approach relied on corridor stacking, which is a well-established processing sequence in vertical seismic profiling. The resulting reflection image outlines bounding layers and can be directly compared to surface-based GPR reflection profiling. Our results of the combined analysis of VRP, traveltimes, amplitudes, and reflections were consistent with independent core and borehole logs as well as GPR reflection profiles, which enabled us to derive a detailed hydro-stratigraphic model as needed, for example, to understand and model groundwater flow and transport.
Ground-based magnetic surveying is a common geophysical method to explore near-surface environments in a non-destructive manner. In many typical applications (such as archaeological prospection), the resulting anomaly maps are often characterized by low signal-to-noise ratios and, thus, the suppression of noise is a key step in data processing. Here, we propose the steering kernel regression (SKR) method to denoise magnetic data sets. SKR has been recently developed to suppress random noise in images and video sequences. The core of the method is the steering kernel function which represents a robust estimate of local image structure. Using such a kernel within an iterative regression based denoising framework, helps to minimize image blurring and to preserve the underlying structures such as edges and corners. Because such filter characteristics are desirable for random noise attenuation in potential field data sets, we apply the SKR method for processing high-resolution ground-based magnetic data as they are typically collected in archaeological applications. We test and evaluate the SKR method using synthetic and field data examples and also compare it to more commonly employed denoising strategies relying, for example, on fixed filter masks (e.g., Gaussian filters). Our results show that the SKR method is successful in removing random and acquisition related noise present in our data. Concurrently, it preserves the local image structure including the amplitudes of anomalies. As demonstrated by derivative based transformations, the mentioned filter characteristics significantly impact subsequent processing steps and, therefore, result in an improved analysis and interpretation of magnetic data. Thus, the method can be considered as a promising and novel approach for denoising ground-based magnetic data.
In near- surface geophysics, ground-based mapping surveys are routinely used in a variety of applications including those from archaeology, civil engineering, hydrology, and soil science. The resulting geophysical anomaly maps of, for example, magnetic or electrical parameters are usually interpreted to laterally delineate subsurface structures such as those related to the remains of past human activities, subsurface utilities and other installations, hydrological properties, or different soil types. To ease the interpretation of such data sets, we have developed a multiscale processing, analysis, and visualization strategy. Our approach relies on a discrete redundant wavelet transform (RWT) implemented using cubic-spline filters and the a trous algorithm, which allows to efficiently compute a multiscale decomposition of 2D data using a series of 1D convolutions. The basic idea of the approach is presented using a synthetic test image, whereas our archaeogeophysical case study from northeast Germany demonstrates its potential to analyze and process rather typical geophysical anomaly maps including magnetic and topographic data. Our vertical-gradient magnetic data show amplitude variations over several orders of magnitude, complex anomaly patterns at various spatial scales, and typical noise patterns, whereas our topographic data show a distinct hill structure superimposed by a microtopographic stripe pattern and random noise. Our results demonstrate that the RWT approach is capable to successfully separate these components and that selected wavelet planes can be scaled and combined so that the reconstructed images allow for a detailed, multiscale structural interpretation also using integrated visualizations of magnetic and topographic data. Because our analysis approach is straightforward to implement without laborious parameter testing and tuning, computationally efficient, and easily adaptable to other geophysical data sets, we believe that it can help to rapidly analyze and interpret different geophysical mapping data collected to address a variety of near-surface applications from engineering practice and research.
Ground-penetrating radar (GPR) is an established geophysical method to explore near-surface sedimentary environments. Interpreting GPR images is largely based on manual procedures following concepts known as GPR facies analysis. We have developed a novel strategy to distinguish GPR facies in a largely automated and more objective manner. First, we calculate 13 textural attributes to quantify GPR reflection characteristics. Then, this database is reduced using principal component analysis. Finally, we image the dominating principal components using composite imaging and classify them using standard clustering methods. The potential of this work-flow is evaluated using a 2D GPR field example collected across stratified glaciofluvial deposits. Our results demonstrate that the derived facies images are well correlated with the composition and the porosity of the sediments as known from independent borehole logs. Our analysis strategy eases and improves the interpretability of GPR data and will help in a variety of geologic and hydrological problems.
Mapping hydrological parameter distributions in high resolution is essential to understand and simulate groundwater flow and contaminant transport. Of particular interest is surface-based ground-penetrating radar (GPR) reflection imaging in electrically resistive sediments because of the expected close link between the subsurface water content and the dielectric permittivity, which controls GPR wave velocity and reflectivity. Conventional tools like common midpoint (CMP) velocity analysis provide physical parameter models of limited resolution only. We present a novel reflection amplitude inversion workflow for surface-based GPR data capable of resolving the subsurface dielectric permittivity and related water content distribution with markedly improved resolution. Our scheme is an adaptation of a seismic reflection impedance inversion scheme to surface-based GPR data. Key is relative-amplitude-preserving data preconditioning including GPR deconvolution, which results in traces with the source-wavelet distortions and propagation effects largely removed. The subsequent inversion for the underlying dielectric permittivity and water content structure is constrained by in situ dielectric permittivity data obtained by direct-push logging. After demonstrating the potential of our novel scheme on a realistic synthetic data set, we apply it to two 2-D 100 MHz GPR profiles acquired over a shallow sedimentary aquifer resulting in water content images of the shallow (3-7 m depth) saturated zone having decimeter resolution.
Three-dimensional hydrostratigraphic models from ground-penetrating radar and direct-push data
(2011)
Three-dimensional models of hydraulic conductivity and porosity are essential to understand and simulate groundwater flow in heterogeneous geological environments. However, considering the inherent limitations of traditional hydrogeological field methods in terms of resolution, alternative field approaches are needed to establish such 3-D models with sufficient accuracy. In this study, we developed a workflow combining 3-D structural information extracted from ground penetrating radar (GPR) images with 1-D in situ physical-property estimates from direct-push (DP) logging to construct a 3-D hydrostratigraphic model. To illustrate this workflow, we collected an similar to 70 m x 90 m 100 MHz 3-D GPR data set over a shallow sedimentary aquifer system resolving six different GPR facies down to similar to 15 m depth. DP logs of the relative dielectric permittivity, the relative hydraulic conductivity, the cone resistance, the sleeve friction and the pore pressure provided crucial data (1) to establish a GPR velocity model for 3-D depth migration and to check the time-to-depth conversion of the GPR data, and (2) to construct a 3-D hydrostratigraphic model. This model was built by assigning porosity values, which were computed from the DP relative dielectric permittivity logs, and DP relative hydraulic conductivity estimates to the identified GPR facies. We conclude that the integration of 3-D GPR structural images and 1-D DP logs of target physical parameters provides an efficient way for detailed 3-D subsurface characterization as needed, for example, for groundwater flow simulations.
Enhancing the resolution and accuracy of surface ground-penetrating radar (GPR) reflection data by inverse filtering to recover a zero-phased band-limited reflectivity image requires a deconvolution technique that takes the mixed-phase character of the embedded wavelet into account. In contrast, standard stochastic deconvolution techniques assume that the wavelet is minimum phase and, hence, often meet with limited success when applied to GPR data. We present a new general-purpose blind deconvolution algorithm for mixed-phase wavelet estimation and deconvolution that (1) uses the parametrization of a mixed-phase wavelet as the convolution of the wavelet's minimum-phase equivalent with a dispersive all-pass filter, (2) includes prior information about the wavelet to be estimated in a Bayesian framework, and (3) relies on the assumption of a sparse reflectivity. Solving the normal equations using the data autocorrelation function provides an inverse filter that optimally removes the minimum-phase equivalent of the wavelet from the data, which leaves traces with a balanced amplitude spectrum but distorted phase. To compensate for the remaining phase errors, we invert in the frequency domain for an all-pass filter thereby taking advantage of the fact that the action of the all-pass filter is exclusively contained in its phase spectrum. A key element of our algorithm and a novelty in blind deconvolution is the inclusion of prior information that allows resolving ambiguities in polarity and timing that cannot be resolved using the sparseness measure alone. We employ a global inversion approach for non-linear optimization to find the all-pass filter phase values for each signal frequency. We tested the robustness and reliability of our algorithm on synthetic data with different wavelets, 1-D reflectivity models of different complexity, varying levels of added noise, and different types of prior information. When applied to realistic synthetic 2-D data and 2-D field data, we obtain images with increased temporal resolution compared to the results of standard processing.
Ice complex deposits are characteristic, ice-rich formations in northern East Siberia and represent an important part in the arctic carbon pool. Recently, these late Quaternary deposits are the objective of numerous investigations typically relying on outcrop and borehole data. Many of these studies can benefit from a 3D structural model of the subsurface for upscaling their observations or for constraining estimations of inventories, such as the local carbon stock. We have addressed this problem of structural imaging by 3D ground-penetrating radar (GPR), which, in permafrost studies, has been primarily used for 2D profiling. We have used a 3D kinematic GPR surveying strategy at a field site located in the New Siberian Archipelago on top of an ice complex. After applying a 3D GPR processing sequence, we were able to trace two horizons at depths below 20 m. Taking available borehole and outcrop data into account, we have interpreted these two features as interfaces of major lithologic units and derived a 3D cryostratigraphic model of the subsurface. Our data example demonstrated that a 3D surveying and processing strategy was crucial at our field site and showed the potential of 3D GPR to image geologic structures in complex ice-rich permafrost landscapes.
For a detailed characterization of near-surface environments, geophysical techniques are increasingly used to support more conventional point-based techniques such as borehole and direct-push logging. Because the underlying parameter relations are often complex, site-specific, or even poorly understood, a remaining challenging task is to link the geophysical parameter models to the actual geotechnical target parameters measured only at selected points. We propose a workflow based on nonparametric regression to establish functional relationships between jointly inverted geophysical parameters and selected geotechnical parameters as measured, for example, by different borehole and direct-push tools. To illustrate our workflow, we present field data collected to characterize a near-surface sedimentary environment Our field data base includes crosshole ground penetrating radar (GPR), seismic P-, and S-wave data sets collected between 25 m deep boreholes penetrating sand- and gravel dominated sediments. Furthermore, different typical borehole and direct-push logs are available. We perform a global joint inversion of traveltimes extracted from the crosshole geophysical data using a recently proposed approach based on particle swarm optimization. Our inversion strategy allows for generating consistent models of GPR, P-wave, and S-wave velocities including an appraisal of uncertainties. We analyze the observed complex relationships between geophysical velocities and target parameter logs using the alternating conditional expectation (ACE) algorithm. This nonparametric statistical tool allows us to perform multivariate regression analysis without assuming a specific functional relation between the variables. We are able to explain selected target parameters such as characteristic grain size values or natural gamma activity by our inverted geophysical data and to extrapolate these parameters to the inter-borehole plane covered by our crosshole experiments. We conclude that the ACE algorithm is a powerful tool to analyze a multivariate petrophysical data base and to develop an understanding of how a multi-parameter geophysical model can be linked and translated to selected geotechnical parameters.
Assessing uncertainty in refraction seismic traveltime inversion using a global inversion strategy
(2015)
To analyse and invert refraction seismic travel time data, different approaches and techniques have been proposed. One common approach is to invert first-break travel times employing local optimization approaches. However, these approaches result in a single velocity model, and it is difficult to assess the quality and to quantify uncertainties and non-uniqueness of the found solution. To address these problems, we propose an inversion strategy relying on a global optimization approach known as particle swarm optimization. With this approach we generate an ensemble of acceptable velocity models, i.e., models explaining our data equally well. We test and evaluate our approach using synthetic seismic travel times and field data collected across a creeping hillslope in the Austrian Alps. Our synthetic study mimics a layered near-surface environment, including a sharp velocity increase with depth and complex refractor topography. Analysing the generated ensemble of acceptable solutions using different statistical measures demonstrates that our inversion strategy is able to reconstruct the input velocity model, including reasonable, quantitative estimates of uncertainty. Our field data set is inverted, employing the same strategy, and we further compare our results with the velocity model obtained by a standard local optimization approach and the information from a nearby borehole. This comparison shows that both inversion strategies result in geologically reasonable models (in agreement with the borehole information). However, analysing the model variability of the ensemble generated using our global approach indicates that the result of the local optimization approach is part of this model ensemble. Our results show the benefit of employing a global inversion strategy to generate near-surface velocity models from refraction seismic data sets, especially in cases where no detailed a priori information regarding subsurface structures and velocity variations is available.
Assessing the human and economic threat introduced by sliding or creeping masses is of major importance in landslide hazard assessment and mitigation. Especially, in the densely populated alpine region unstable hillslopes represent a major hazard to men and infrastructure. Detailed knowledge, especially, of the dominant site-specific controlling factors such as subsurface architecture and geology is thereby key in assessing slope vulnerability. In order to quantify the geological variations at a creeping hillslope in the Austrian Alps, we have collected six 2D refraction seismic profiles. We propose using a layer-based inversion strategy to reconstruct P-wave velocity models from first arrival times. Considering the geological complexity at such sites, the selected inversion approach eases the interpretability of geological structures given intrinsic optimization for only a discrete, user-defined, number of layers. As the applied layer-based inversion approach fits our travel time data equally well as traditional smooth inversion approaches, it represents a feasible mean to summarize the structural complexity often present at such sites. Analysis of the inversion results illustrates that bedrock topography clearly deviates from a previously assumed planar surface and exhibits distinct variations across the slope extension. Bedrock topography additionally impacts the intermediate geological units and, thus, this information is critical for further analyses such as geomechanical modeling. (C) 2012 Elsevier B.V. All rights reserved.
Inversions of an individual geophysical data set can be highly nonunique, and it is generally difficult to determine petrophysical parameters from geophysical data. We show that both issues can be addressed by adopting a statistical multiparameter approach that requires the acquisition, processing, and separate inversion of two or more types of geophysical data. To combine information contained in the physical-property models that result from inverting the individual data sets and to estimate the spatial distribution of petrophysical parameters in regions where they are known at only a few locations. we demonstrate the potential of the fuzzy c-means (FCM) clustering technique. After testing this new approach on synthetic data, we apply it to limited crosshole georadar, crosshole seismic, gamma-log, and slug-test data acquired within a shallow alluvial aquifer. The derived multiparameter model effectively outlines the major sedimentary units observed in numerous boreholes and provides plausible estimates for the spatial distributions of gamma-ray emitters and hydraulic conductivity
In many near-surface geophysical studies it is now common practice to collect co-located disparate geophysical data sets to explore subsurface structures. Reconstruction of physical parameter distributions underlying the available geophysical data sets usually requires the use of tomographic reconstruction techniques. To improve the quality of the obtained models, the information content of all data sets should be considered during the model generation process, e.g., by employing joint or cooperative inversion approaches. Here, we extend the zonal cooperative inversion methodology based on fuzzy c-means cluster analysis and conventional single-input data set inversion algorithms for the cooperative inversion of data sets with partially co-located model areas. This is done by considering recent developments in fuzzy c-means cluster analysis. Additionally, we show how supplementary a priori information can be incorporated in an automated fashion into the zonal cooperative inversion approach to further constrain the inversion. The only requirement is that this a priori information can be expressed numerically; e.g., by physical parameters or indicator variables. We demonstrate the applicability of the modified zonal cooperative inversion approach using synthetic and field data examples. In these examples, we cooperatively invert S- and P-wave traveltime data sets with partially co-located model areas using water saturation information expressed by indicator variables as additional a priori information. The approach results in a zoned multi-parameter model, which is consistent with all available information given to the zonal cooperative inversion and outlines the major subsurface units. In our field example, we further compare the obtained zonal model to sparsely available borehole and direct-push logs. This comparison provides further confidence in our zonal cooperative inversion model because the borehole and direct-push logs indicate a similar zonation.
Geophysical techniques offer the potential to tomographically image physical parameter variations in the ground in two or three dimensions. Due to the limited number and accuracy of the recorded data, geophysical model generation by inversion suffers ambiguity. Linking the model generation process of disparate data by jointly inverting two or more data sets allows for improved model reconstruction. Fully nonlinear inversion using optimization techniques searching the solution space of the inverse problem globally enables quantitative assessment of the ambiguity inherent to the model reconstruction. We used two different multiobjective particle swarm optimization approaches to jointly invert synthetic crosshole tomographic data sets comprising radar and P-wave traveltimes, respectively. Beginning with a nonlinear joint inversion founded on the principle of Pareto optimality and game theoretic concepts, we obtained a set of Pareto-optimal solutions comprising commonly structured radar and P-wave velocity models for low computational costs. However, the efficiency of the approach goes along with some risk of achieving a final model ensemble not adequately illustrating the ambiguity inherent to the model reconstruction process. Taking advantage of the results of the first approach, we inverted the database using a different nonlinear joint-inversion approach reducing the multiobjective optimization problem to a single-objective one. Computational costs were significantly higher, but the final models were obtained mutually independently allowing for objective appraisal of model parameter determination. Despite the high computational effort, the approach was found to be an efficient nonlinear joint-inversion formulation compared to what could be extracted from individual nonlinear inversions of both data sets.
Ground-penetrating radar (GPR) is a standard geophysical technique used to image near-surface structures in sedimentary environments. In such environments, GPR data acquisition and processing are increasingly following 3D strategies. However, the processed GPR data volumes are typically still interpreted using selected 2D slices and manual concepts such as GPR facies analyses. In seismic volume interpretation, the application of (semi-)automated and reproducible approaches such as 3D attribute analyses as well as the production of attribute-based facies models are common practices today. In contrast, the field of 3D GPR attribute analyses and corresponding facies models is largely untapped. We have developed and applied a workflow to produce 3D attribute-based GPR facies models comprising the dominant sedimentary reflection patterns in a GPR volume, which images complex sandy structures on the dune island of Spiekeroog (Northern Germany). After presenting our field site and details regarding our data acquisition and processing, we calculate and filter 3D texture attributes to generate a database comprising the dominant texture features of our GPR data. Then, we perform a dimensionality reduction of this database to obtain meta texture attributes, which we analyze and integrate using composite imaging and (also considering additional geometric information) fuzzy c-means cluster analysis resulting in a classified GPR facies model. Considering our facies model and a corresponding GPR facies chart, we interpret our GPR data set in terms of near-surface sedimentary units, the corresponding depositional environments, and the recent formation history at our field site. Thus, we demonstrate the potential of our workflow, which represents a novel and clear strategy to perform a more objective and consistent interpretation of 3D GPR data collected across different sedimentary environments.