TY - JOUR A1 - Klose, Tim A1 - Guillemoteau, Julien A1 - Vignoli, Giulio A1 - Walter, Judith A1 - Herrmann, Andreas A1 - Tronicke, Jens T1 - Structurally constrained inversion by means of a Minimum Gradient Support regularizer BT - examples of FD-EMI data inversion constrained by GPR reflection data JF - Geophysical journal international N2 - Many geophysical inverse problems are known to be ill-posed and, thus, requiring some kind of regularization in order to provide a unique and stable solution. A possible approach to overcome the inversion ill-posedness consists in constraining the position of the model interfaces. For a grid-based parameterization, such a structurally constrained inversion can be implemented by adopting the usual smooth regularization scheme in which the local weight of the regularization is reduced where an interface is expected. By doing so, sharp contrasts are promoted at interface locations while standard smoothness constraints keep affecting the other regions of the model. In this work, we present a structurally constrained approach and test it on the inversion of frequency-domain electromagnetic induction (FD-EMI) data using a regularization approach based on the Minimum Gradient Support stabilizer, which is capable to promote sharp transitions everywhere in the model, i.e., also in areas where no structural a prioriinformation is available. Using 1D and 2D synthetic data examples, we compare the proposed approach to a structurally constrained smooth inversion as well as to more standard (i.e., not structurally constrained) smooth and sharp inversions. Our results demonstrate that the proposed approach helps in finding a better and more reliable reconstruction of the subsurface electrical conductivity distribution, including its structural characteristics. Furthermore, we demonstrate that it allows to promote sharp parameter variations in areas where no structural information are available. Lastly, we apply our structurally constrained scheme to FD-EMI field data collected at a field site in Eastern Germany to image the thickness of peat deposits along two selected profiles. In this field example, we use collocated constant offset ground-penetrating radar (GPR) data to derive structural a priori information to constrain the inversion of the FD-EMI data. The results of this case study demonstrate the effectiveness and flexibility of the proposed approach. KW - Controlled source electromagnetics (CSEM) KW - Inverse theory KW - Electrical properties KW - Ground penetrating radar KW - Frequency Domain Electromagnetics KW - Inversion Y1 - 2023 U6 - https://doi.org/10.1093/gji/ggad041 SN - 0956-540X SN - 1365-246X VL - 233 IS - 3 SP - 1938 EP - 1949 PB - Oxford Univ. Press CY - Oxford ER - TY - CHAP A1 - Simon, François-Xavier A1 - Papadopoulos, Nikos A1 - Guillemoteau, Julien A1 - Oikonomou, Dimitris A1 - Simirdanis, Kleanthis T1 - Multi-frequency loop electromagnetic system measurement on shallow offshore archaeological site of Oulos T2 - ArcheoSciences : revue d'archéométrie / Groupe des Méthodes Pluridisciplinaires Contribuant à l'Archéologie (GMPCA) KW - hallow offshore KW - multi-frequency KW - electromagnetic KW - modelling KW - case study Y1 - 2021 U6 - https://doi.org/10.4000/archeosciences.9690 SN - 1960-1360 SN - 2104-3728 VL - 45 IS - 1 SP - 215 EP - 218 PB - Presses Universitaires de Rennes CY - Rennes ER - TY - JOUR A1 - Klose, Tim A1 - Guillemoteau, Julien A1 - Vignoli, Giulio A1 - Tronicke, Jens T1 - Laterally constrained inversion (LCI) of multi-configuration EMI data with tunable sharpness JF - Journal of applied geophysics N2 - Frequency-domain electromagnetic (FDEM) data are commonly inverted to characterize subsurface geoelectrical properties using smoothness constraints in 1D inversion schemes assuming a layered medium. Smoothness constraints are suitable for imaging gradual transitions of subsurface geoelectrical properties caused, for example, by varying sand, clay, or fluid content. However, such inversion approaches are limited in characterizing sharp interfaces. Alternative regularizations based on the minimum gradient support (MGS) stabilizers can, instead, be used to promote results with different levels of smoothness/sharpness selected by simply acting on the so-called focusing parameter. The MGS regularization has been implemented for different kinds of geophysical data inversion strategies. However, concerning FDEM data, the MGS regularization has only been implemented for vertically constrained inversion (VCI) approaches but not for laterally constrained inversion (LCI) approaches. We present a novel LCI approach for FDEM data using the MGS regularization for the vertical and lateral direction. Using synthetic and field data examples, we demonstrate that our approach can efficiently and automatically provide a set of model solutions characterized by different levels of sharpness and variable lateral consistencies. In terms of data misfit, the obtained set of solutions contains equivalent models allowing us also to investigate the non-uniqueness of FDEM data inversion. KW - frequency-domain electromagnetics KW - laterally constrained inversion KW - minimum gradient support regularization KW - peat characterization Y1 - 2022 U6 - https://doi.org/10.1016/j.jappgeo.2021.104519 SN - 0926-9851 VL - 196 PB - Elsevier CY - Amsterdam ER - TY - JOUR A1 - Guillemoteau, Julien A1 - Vignoli, Giulio A1 - Barreto, Jeniffer A1 - Sauvin, Guillaume T1 - Sparse laterally constrained inversion of surface-wave dispersion curves via minimum gradient support regularization JF - Geophysics N2 - We have developed a 1D laterally constrained inversion of surface-wave dispersion curves based on the minimum gradient support regularization, which allows solutions with tunable sharpness in the vertical and horizontal directions. The forward modeling consists of a finite-elements approach incorporated in a flexible nonparametric gradient-based inversion scheme, which has already demonstrated good stability and convergence capabilities when tested on other kinds of data. Our deterministic inversion procedure is performed in the shear-wave velocity log space as we noticed that the associated Jacobian indicates a reduced model dependency, and this, in turn, decreases the risks of local nonconvexity. We show several synthetics and one field example to demonstrate the effectiveness and the applicability of the proposed approach. KW - surface wave, inversion, near surface Y1 - 2022 U6 - https://doi.org/10.1190/GEO2021-0247.1 SN - 0016-8033 SN - 1942-2156 VL - 87 IS - 3 SP - R281 EP - R289 PB - Society of Exploration Geophysicists CY - Tulsa, Okla. ER - TY - JOUR A1 - Tronicke, Jens A1 - Allroggen, Niklas A1 - Biermann, Felix A1 - Fanselow, Florian A1 - Guillemoteau, Julien A1 - Krauskopf, Christof A1 - Lück, Erika T1 - Rapid multiscale analysis of near-surface geophysical anomaly maps BT - application to an archaeogeophysical data set JF - Geophysics N2 - 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. KW - archaeology KW - case history KW - near surface KW - magnetics KW - decomposition Y1 - 2020 U6 - https://doi.org/10.1190/GEO2019-0564.1 SN - 0016-8033 SN - 1942-2156 VL - 85 IS - 4 SP - B109 EP - B118 PB - Society of Exploration Geophysicists CY - Tulsa, Okla. ER - TY - JOUR A1 - Arboleda-Zapata, Mauricio A1 - Guillemoteau, Julien A1 - Tronicke, Jens T1 - A comprehensive workflow to analyze ensembles of globally inverted 2D electrical resistivity models JF - Journal of applied geophysics N2 - Electrical resistivity tomography (ERT) aims at imaging the subsurface resistivity distribution and provides valuable information for different geological, engineering, and hydrological applications. To obtain a subsurface resistivity model from measured apparent resistivities, stochastic or deterministic inversion procedures may be employed. Typically, the inversion of ERT data results in non-unique solutions; i.e., an ensemble of different models explains the measured data equally well. In this study, we perform inference analysis of model ensembles generated using a well-established global inversion approach to assess uncertainties related to the nonuniqueness of the inverse problem. Our interpretation strategy starts by establishing model selection criteria based on different statistical descriptors calculated from the data residuals. Then, we perform cluster analysis considering the inverted resistivity models and the corresponding data residuals. Finally, we evaluate model uncertainties and residual distributions for each cluster. To illustrate the potential of our approach, we use a particle swarm optimization (PSO) algorithm to obtain an ensemble of 2D layer-based resistivity models from a synthetic data example and a field data set collected in Loon-Plage, France. Our strategy performs well for both synthetic and field data and allows us to extract different plausible model scenarios with their associated uncertainties and data residual distributions. Although we demonstrate our workflow using 2D ERT data and a PSObased inversion approach, the proposed strategy is general and can be adapted to analyze model ensembles generated from other kinds of geophysical data and using different global inversion approaches. KW - Near-surface geophysics KW - Electrical resistivity tomography KW - Non-uniqueness KW - Global inversion KW - Particle swarm optimization KW - Ensemble KW - analysis Y1 - 2021 U6 - https://doi.org/10.1016/j.jappgeo.2021.104512 SN - 0926-9851 SN - 1879-1859 VL - 196 PB - Elsevier CY - Amsterdam ER - TY - JOUR A1 - Garcin, Yannick A1 - Schildgen, Taylor F. A1 - Acosta, Veronica Torres A1 - Melnick, Daniel A1 - Guillemoteau, Julien A1 - Willenbring, Jane A1 - Strecker, Manfred T1 - Short-lived increase in erosion during the African Humid Period BT - evidence from the northern Kenya Rift JF - Earth & planetary science letters N2 - The African Humid Period (AHP) between similar to 15 and 5.5 cal. kyr BP caused major environmental change in East Africa, including filling of the Suguta Valley in the northern Kenya Rift with an extensive (similar to 2150 km(2)), deep (similar to 300 m) lake. Interfingering fluvio-lacustrine deposits of the Baragoi paleo-delta provide insights into the lake-level history and how erosion rates changed during this time, as revealed by delta-volume estimates and the concentration of cosmogenic Be-10 in fluvial sand. Erosion rates derived from delta-volume estimates range from 0.019 to 0.03 mm yr(-1). Be-10-derived paleo-erosion rates at similar to 11.8 cal. kyr BP ranged from 0.035 to 0.086 mm yr(-1), and were 2.7 to 6.6 times faster than at present. In contrast, at similar to 8.7 cal. kyr BP, erosion rates were only 1.8 times faster than at present. Because Be-10-derived erosion rates integrate over several millennia; we modeled the erosion-rate history that best explains the 10Be data using established non-linear equations that describe in situ cosmogenic isotope production and decay. Two models with different temporal constraints (15-6.7 and 12-6.7 kyr) suggest erosion rates that were 25 to 300 times higher than the initial erosion rate (pre-delta formation). That pulse of high erosion rates was short (similar to 4 kyr or less) and must have been followed by a rapid decrease in rates while climate remained humid to reach the modern Be-10-based erosion rate of,similar to 0.013 mm yr(-1). Our simulations also flag the two highest Be-10-derived erosion rates at 11.8 kyr BP related to nonuniform catchment erosion. These changes in erosion rates and processes during the AHP may reflect a strong increase in precipitation, runoff, and erosivity at the arid-to-humid transition either at 15 or similar to 12 cal. kyr BP, before the landscape stabilized again, possibly due to increased soil production and denser vegetation. KW - northern Kenya Rift KW - Baragoi KW - paleo-delta KW - African Humid Period KW - erosion KW - Be-10 Y1 - 2017 U6 - https://doi.org/10.1016/j.epsl.2016.11.017 SN - 0012-821X SN - 1385-013X VL - 459 SP - 58 EP - 69 PB - Elsevier CY - Amsterdam ER - TY - JOUR A1 - Klose, Tim A1 - Guillemoteau, Julien A1 - Simon, Francois-Xavier A1 - Tronicke, Jens T1 - Toward subsurface magnetic permeability imaging with electromagnetic induction sensors BT - Sensitivity computation and reconstruction of measured data JF - Geophysics N2 - In near-surface geophysics, small portable loop-loop electro-magnetic induction (EMI) sensors using harmonic sources with a constant and rather small frequency are increasingly used to investigate the electrical properties of the subsurface. For such sensors, the influence of electrical conductivity and magnetic permeability on the EMI response is well-understood. Typically, data analysis focuses on reconstructing an electrical conductivity model by inverting the out-of-phase response. However, in a variety of near-surface applications, magnetic permeability (or susceptibility) models derived from the in-phase (IP) response may provide important additional information. In view of developing a fast 3D inversion procedure of the IP response for a dense grid of measurement points, we first analyze the 3D sensitivity functions associated with a homogeneous permeable half-space. Then, we compare synthetic data computed using a linear forward-modeling method based on these sensitivity functions with synthetic data computed using full nonlinear forward-modeling methods. The results indicate the correctness and applicability of our linear forward-modeling approach. Furthermore, we determine the advantages of converting IP data into apparent permeability, which, for example, allows us to extend the applicability of the linear forward-modeling method to high-magnetic environments. Finally, we compute synthetic data with the linear theory for a model consisting of a controlled magnetic target and compare the results with field data collected with a four-configuration loop-loop EMI sensor. With this field-scale experiment, we determine that our linear forward-modeling approach can reproduce measured data with sufficiently small error, and, thus, it represents the basis for developing efficient inversion approaches. KW - Electromagnetics KW - Imaging KW - Magnetic+Susceptibility KW - Near+Surface KW - Modeling Y1 - 2018 U6 - https://doi.org/10.1190/GEO2017-0827.1 SN - 0016-8033 SN - 1942-2156 VL - 83 IS - 5 SP - E335 EP - E345 PB - Society of Exploration Geophysicists CY - Tulsa ER - TY - GEN A1 - Arnous, Ahmad A1 - Zeckra, Martin A1 - Venerdini, Agostina A1 - Alvarado, Patricia A1 - Arrowsmith, Ramón A1 - Guillemoteau, Julien A1 - Landgraf, Angela A1 - Gutiérrez, Adolfo Antonio A1 - Strecker, Manfred T1 - Neotectonic Activity in the Low-Strain Broken Foreland (Santa Bárbara System) of the North-Western Argentinean Andes (26°S) T2 - Postprints der Universität Potsdam : Mathematisch Naturwissenschaftliche Reihe N2 - Uplift in the broken Andean foreland of the Argentine Santa Bárbara System (SBS) is associated with the contractional reactivation of basement anisotropies, similar to those reported from the thick-skinned Cretaceous-Eocene Laramide province of North America. Fault scarps, deformed Quaternary deposits and landforms, disrupted drainage patterns, and medium-sized earthquakes within the SBS suggest that movement along these structures may be a recurring phenomenon, with yet to be defined repeat intervals and rupture lengths. In contrast to the Subandes thrust belt farther north, where eastward-migrating deformation has generated a well-defined thrust front, the SBS records spatiotemporally disparate deformation along structures that are only known to the first order. We present herein the results of geomorphic desktop analyses, structural field observations, and 2D electrical resistivity tomography and seismic-refraction tomography surveys and an interpretation of seismic reflection profiles across suspected fault scarps in the sedimentary basins adjacent to the Candelaria Range (CR) basement uplift, in the south-central part of the SBS. Our analysis in the CR piedmont areas reveals consistency between the results of near-surface electrical resistivity and seismic-refraction tomography surveys, the locations of prominent fault scarps, and structural geometries at greater depth imaged by seismic reflection data. We suggest that this deformation is driven by deep-seated blind thrusting beneath the CR and associated regional warping, while shortening involving Mesozoic and Cenozoic sedimentary strata in the adjacent basins was accommodated by layer-parallel folding and flexural-slip faults that cut through Quaternary landforms and deposits at the surface. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 1008 Y1 - 2020 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-480183 SN - 1866-8372 IS - 1008 SP - 1 EP - 25 ER - TY - JOUR A1 - Arnous, Ahmad A1 - Zeckra, Martin A1 - Venerdini, Agostina A1 - Alvarado, Patricia A1 - Arrowsmith, Ramón A1 - Guillemoteau, Julien A1 - Landgraf, Angela A1 - Gutiérrez, Adolfo Antonio A1 - Strecker, Manfred T1 - Neotectonic Activity in the Low-Strain Broken Foreland (Santa Bárbara System) of the North-Western Argentinean Andes (26°S) JF - Lithosphere N2 - Uplift in the broken Andean foreland of the Argentine Santa Bárbara System (SBS) is associated with the contractional reactivation of basement anisotropies, similar to those reported from the thick-skinned Cretaceous-Eocene Laramide province of North America. Fault scarps, deformed Quaternary deposits and landforms, disrupted drainage patterns, and medium-sized earthquakes within the SBS suggest that movement along these structures may be a recurring phenomenon, with yet to be defined repeat intervals and rupture lengths. In contrast to the Subandes thrust belt farther north, where eastward-migrating deformation has generated a well-defined thrust front, the SBS records spatiotemporally disparate deformation along structures that are only known to the first order. We present herein the results of geomorphic desktop analyses, structural field observations, and 2D electrical resistivity tomography and seismic-refraction tomography surveys and an interpretation of seismic reflection profiles across suspected fault scarps in the sedimentary basins adjacent to the Candelaria Range (CR) basement uplift, in the south-central part of the SBS. Our analysis in the CR piedmont areas reveals consistency between the results of near-surface electrical resistivity and seismic-refraction tomography surveys, the locations of prominent fault scarps, and structural geometries at greater depth imaged by seismic reflection data. We suggest that this deformation is driven by deep-seated blind thrusting beneath the CR and associated regional warping, while shortening involving Mesozoic and Cenozoic sedimentary strata in the adjacent basins was accommodated by layer-parallel folding and flexural-slip faults that cut through Quaternary landforms and deposits at the surface. Y1 - 2019 U6 - https://doi.org/10.2113/2020/8888588 SN - 1947-4253 SN - 1941-8264 VL - 2020 IS - 1 SP - 1 EP - 25 PB - GSA CY - Boulder, Colo. ER -