@article{AllroggenBeiterTronicke2020, author = {Allroggen, Niklas and Beiter, Daniel and Tronicke, Jens}, title = {Ground-penetrating radar monitoring of fast subsurface processes}, series = {Geophysics}, volume = {85}, journal = {Geophysics}, number = {3}, publisher = {Society of Exploration Geophysicists}, address = {Tulsa}, issn = {0016-8033}, doi = {10.1190/GEO2019-0737.1}, pages = {A19 -- A23}, year = {2020}, abstract = {Earth and environmental sciences rely on detailed information about subsurface processes. Whereas geophysical techniques typically provide highly resolved spatial images, monitoring subsurface processes is often associated with enormous effort and, therefore, is usually limited to point information in time or space. Thus, the development of spatial and temporal continuous field monitoring methods is a major challenge for the understanding of subsurface processes. We have developed a novel method for ground-penetrating-radar (GPR) reflection monitoring of subsurface flow processes under unsaturated conditions and applied it to a hydrological infiltration experiment performed across a periglacial slope deposit in northwest Luxembourg. Our approach relies on a spatial and temporal quasicontinuous data recording and processing, followed by an attribute analysis based on analyzing differences between individual time steps. The results demonstrate the ability of time-lapse GPR monitoring to visualize the spatial and temporal dynamics of preferential flow processes with a spatial resolution in the order of a few decimeters and temporal resolution in the order of a few minutes. We observe excellent agreement with water table information originating from different boreholes. This demonstrates the potential of surface-based GPR reflection monitoring to observe the spatiotemporal dynamics of water movements in the subsurface. It provides valuable, and so far not accessible, information for example in the field of hydrology and pedology that allows studying the actual subsurface processes rather than deducing them from point information.}, language = {en} } @article{AllroggenBoothBakeretal.2019, author = {Allroggen, Niklas Robin and Booth, Adam D. and Baker, Sandra E. and Ellwood, Stephen A. and Tronicke, Jens}, title = {High-resolution imaging and monitoring of animal tunnels using 3D ground-penetrating radar}, series = {Near surface geophysics}, volume = {17}, journal = {Near surface geophysics}, number = {3}, publisher = {Wiley}, address = {Hoboken}, issn = {1569-4445}, doi = {10.1002/nsg.12039}, pages = {291 -- 298}, year = {2019}, abstract = {Ground-penetrating radar is widely used to provide highly resolved images of subsurface sedimentary structures, with implications for processes active in the vadose zone. Frequently overlooked among these structures are tunnels excavated by fossorial animals (e.g., moles). We present two repeated ground-penetrating radar surveys performed a year apart in 2016 and 2017. Careful three-dimensional data processing reveals, in each data set, a pattern of elongated structures that are interpreted as a subsurface mole tunnel network. Our data demonstrate the ability of three-dimensional ground-penetrating radar imaging to non-invasively delineate the small animal tunnels (similar to 5 cm diameter) at a higher spatial and geolocation resolution than has previously been achieved. In turn, this makes repeated surveys and, therefore, long-term monitoring possible. Our results offer valuable insight into the understanding of the near-surface and showcase a potential new application for a geophysical method as well as a non-invasive method of ecological surveying.}, language = {en} } @article{AllroggenTronicke2016, author = {Allroggen, Niklas and Tronicke, Jens}, title = {Attribute-based analysis of time-lapse ground-penetrating radar data}, series = {Geophysics}, volume = {81}, journal = {Geophysics}, publisher = {Society of Exploration Geophysicists}, address = {Tulsa}, issn = {0016-8033}, doi = {10.1190/GEO2015-0171.1}, pages = {H1 -- H8}, year = {2016}, abstract = {Analysis of time-lapse ground-penetrating radar (GPR) data can provide information regarding subsurface hydrological processes, such as preferential flow. However, the analysis of time-lapse data is often limited by data quality; for example, for noisy input data, the interpretation of difference images is often difficult. Motivated by modern image-processing tools, we have developed two robust GPR attributes, which allow us to distinguish amplitude (contrast similarity) and time-shift (structural similarity) variations related to differences between individual time-lapse GPR data sets. We tested and evaluated our attributes using synthetic data of different complexity. Afterward, we applied them to a field data example, in which subsurface flow was induced by an artificial rainfall event. For all examples, we identified our structural similarity attribute to be a robust measure for highlighting time-lapse changes also in data with low signal-to-noise ratios. We determined that our new attribute-based workflow is a promising tool to analyze time-lapse GPR data, especially for imaging subsurface hydrological processes.}, language = {en} } @article{AllroggenTronickeDelocketal.2015, author = {Allroggen, Niklas and Tronicke, Jens and Delock, Marcel and B{\"o}niger, Urs}, title = {Topographic migration of 2D and 3D ground-penetrating radar data considering variable velocities}, series = {Near surface geophysics}, volume = {13}, journal = {Near surface geophysics}, number = {3}, publisher = {European Association of Geoscientists \& Engineers}, address = {Houten}, issn = {1569-4445}, doi = {10.3997/1873-0604.2014037}, pages = {4}, year = {2015}, abstract = {We present a 2D/3D topographic migration scheme for ground-penetrating radar (GPR) data which is able to account for variable velocities by using the root mean square (rms) velocity approximation. We test our migration scheme using a synthetic 2D example and compare our migrated image to the results obtained using common GPR migration approaches. Furthermore, we apply it to 2D and 3D field data. These examples are recorded across common subsurface settings including surface topography and variations in the GPR subsurface velocity field caused by a shallow ground water table. In such field settings, our migration strategy provides well focused images of commonoffset GPR data without the need for a detailed interval velocity model. The synthetic and field examples demonstrate that our topographic migration scheme allows for accurate GPR imaging in the presence of variations in surface topography and subsurface velocity.}, language = {en} } @article{AllroggenvanSchaikTronicke2015, author = {Allroggen, Niklas and van Schaik, N. Loes M. B. and Tronicke, Jens}, title = {4D ground-penetrating radar during a plot scale dye tracer experiment}, series = {Journal of applied geophysics}, volume = {118}, journal = {Journal of applied geophysics}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0926-9851}, doi = {10.1016/j.jappgeo.2015.04.016}, pages = {139 -- 144}, year = {2015}, abstract = {Flow phenomena in the unsaturated zone are highly variable in time and space. Thus, it is challenging to measure and monitor such processes under field conditions. Here, we present a new setup and interpretation approach for combining a dye tracer experiment with a 4D ground-penetrating radar (GPR) survey. Therefore, we designed a rainfall experiment during which we measured three surface-based 3D GPR surveys using a pair of 500 MHz antennas. Such a survey setup requires accurate acquisition and processing techniquesto extract time-lapse information supporting the interpretation of selected cross-sections photographed after excavating the site. Our results reveal patterns of traveltime changes in the measured GPR data, which are associated with soil moisture changes. As distinct horizons are present at our site, such changes can be quantified and transferred into changes in total soil moisture content. Our soil moisture estimates are similar to the amount of infiltrated water, which confirms our experimental approach and makes us confident for further developing this strategy, especially, with respect to improving the temporal and spatial resolution. (C) 2015 Elsevier B.V. All rights reserved.}, language = {en} } @article{AngelopoulosOverduinWestermannetal.2020, author = {Angelopoulos, Michael and Overduin, Pier Paul and Westermann, Sebastian and Tronicke, Jens and Strauss, Jens and Schirrmeister, Lutz and Biskaborn, Boris and Liebner, Susanne and Maksimov, Georgii and Grigoriev, Mikhail N. and Grosse, Guido}, title = {Thermokarst lake to lagoon transitions in Eastern Siberia}, series = {Journal of geophysical research : Earth surface}, volume = {125}, journal = {Journal of geophysical research : Earth surface}, number = {10}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9003}, doi = {10.1029/2019JF005424}, pages = {21}, year = {2020}, abstract = {As the Arctic coast erodes, it drains thermokarst lakes, transforming them into lagoons, and, eventually, integrates them into subsea permafrost. Lagoons represent the first stage of a thermokarst lake transition to a marine setting and possibly more saline and colder upper boundary conditions. In this research, borehole data, electrical resistivity surveying, and modeling of heat and salt diffusion were carried out at Polar Fox Lagoon on the Bykovsky Peninsula, Siberia. Polar Fox Lagoon is a seasonally isolated water body connected to Tiksi Bay through a channel, leading to hypersaline waters under the ice cover. The boreholes in the center of the lagoon revealed floating ice and a saline cryotic bed underlain by a saline cryotic talik, a thin ice-bearing permafrost layer, and unfrozen ground. The bathymetry showed that most of the lagoon had bedfast ice in spring. In bedfast ice areas, the electrical resistivity profiles suggested that an unfrozen saline layer was underlain by a thick layer of refrozen talik. The modeling showed that thermokarst lake taliks can refreeze when submerged in saltwater with mean annual bottom water temperatures below or slightly above 0 degrees C. This occurs, because the top-down chemical degradation of newly formed ice-bearing permafrost is slower than the refreezing of the talik. Hence, lagoons may precondition taliks with a layer of ice-bearing permafrost before encroachment by the sea, and this frozen layer may act as a cap on gas migration out of the underlying talik.}, language = {en} } @misc{AngermannJackischAllroggenetal.2017, author = {Angermann, Lisa and Jackisch, Conrad and Allroggen, Niklas and Sprenger, Matthias and Zehe, Erwin and Tronicke, Jens and Weiler, Markus and Blume, Theresa}, title = {Form and function in hillslope hydrology}, series = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {658}, issn = {1866-8372}, doi = {10.25932/publishup-41916}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-419161}, pages = {22}, year = {2017}, abstract = {The phrase form and function was established in architecture and biology and refers to the idea that form and functionality are closely correlated, influence each other, and co-evolve. We suggest transferring this idea to hydrological systems to separate and analyze their two main characteristics: their form, which is equivalent to the spatial structure and static properties, and their function, equivalent to internal responses and hydrological behavior. While this approach is not particularly new to hydrological field research, we want to employ this concept to explicitly pursue the question of what information is most advantageous to understand a hydrological system. We applied this concept to subsurface flow within a hillslope, with a methodological focus on function: we conducted observations during a natural storm event and followed this with a hillslope-scale irrigation experiment. The results are used to infer hydrological processes of the monitored system. Based on these findings, the explanatory power and conclusiveness of the data are discussed. The measurements included basic hydrological monitoring methods, like piezometers, soil moisture, and discharge measurements. These were accompanied by isotope sampling and a novel application of 2-D time-lapse GPR (ground-penetrating radar). The main finding regarding the processes in the hillslope was that preferential flow paths were established quickly, despite unsaturated conditions. These flow paths also caused a detectable signal in the catchment response following a natural rainfall event, showing that these processes are relevant also at the catchment scale. Thus, we conclude that response observations (dynamics and patterns, i.e., indicators of function) were well suited to describing processes at the observational scale. Especially the use of 2-D time-lapse GPR measurements, providing detailed subsurface response patterns, as well as the combination of stream-centered and hillslope-centered approaches, allowed us to link processes and put them in a larger context. Transfer to other scales beyond observational scale and generalizations, however, rely on the knowledge of structures (form) and remain speculative. The complementary approach with a methodological focus on form (i.e., structure exploration) is presented and discussed in the companion paper by Jackisch et al. (2017).}, language = {en} } @article{AngermannJackischAllroggenetal.2017, author = {Angermann, Lisa and Jackisch, Conrad and Allroggen, Niklas and Sprenger, Matthias and Zehe, Erwin and Tronicke, Jens and Weiler, Markus and Blume, Theresa}, title = {Form and function in hillslope hydrology: characterization of subsurface flow based on response observations}, series = {Hydrology and earth system sciences : HESS}, volume = {21}, journal = {Hydrology and earth system sciences : HESS}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1027-5606}, doi = {10.5194/hess-21-3727-2017}, pages = {3727 -- 3748}, year = {2017}, abstract = {The phrase form and function was established in architecture and biology and refers to the idea that form and functionality are closely correlated, influence each other, and co-evolve. We suggest transferring this idea to hydrological systems to separate and analyze their two main characteristics: their form, which is equivalent to the spatial structure and static properties, and their function, equivalent to internal responses and hydrological behavior. While this approach is not particularly new to hydrological field research, we want to employ this concept to explicitly pursue the question of what information is most advantageous to understand a hydrological system. We applied this concept to subsurface flow within a hillslope, with a methodological focus on function: we conducted observations during a natural storm event and followed this with a hillslope-scale irrigation experiment. The results are used to infer hydrological processes of the monitored system. Based on these findings, the explanatory power and conclusiveness of the data are discussed. The measurements included basic hydrological monitoring methods, like piezometers, soil moisture, and discharge measurements. These were accompanied by isotope sampling and a novel application of 2-D time-lapse GPR (ground-penetrating radar). The main finding regarding the processes in the hillslope was that preferential flow paths were established quickly, despite unsaturated conditions. These flow paths also caused a detectable signal in the catchment response following a natural rainfall event, showing that these processes are relevant also at the catchment scale. Thus, we conclude that response observations (dynamics and patterns, i.e., indicators of function) were well suited to describing processes at the observational scale. Especially the use of 2-D time-lapse GPR measurements, providing detailed subsurface response patterns, as well as the combination of stream-centered and hillslope-centered approaches, allowed us to link processes and put them in a larger context. Transfer to other scales beyond observational scale and generalizations, however, rely on the knowledge of structures (form) and remain speculative. The complementary approach with a methodological focus on form (i.e., structure exploration) is presented and discussed in the companion paper by Jackisch et al. (2017).}, language = {en} } @article{ArboledaZapataAngelopoulosOverduinetal.2022, author = {Arboleda-Zapata, Mauricio and Angelopoulos, Michael and Overduin, Pier Paul and Grosse, Guido and Jones, Benjamin M. and Tronicke, Jens}, title = {Exploring the capabilities of electrical resistivity tomography to study subsea permafrost}, series = {The Cryosphere}, volume = {16}, journal = {The Cryosphere}, publisher = {Copernicus}, address = {Katlenburg-Lindau}, issn = {1994-0424}, doi = {10.5194/tc-16-4423-2022}, pages = {4423 -- 4445}, year = {2022}, abstract = {Sea level rise and coastal erosion have inundated large areas of Arctic permafrost. Submergence by warm and saline waters increases the rate of inundated permafrost thaw compared to sub-aerial thawing on land. Studying the contact between the unfrozen and frozen sediments below the seabed, also known as the ice-bearing permafrost table (IBPT), provides valuable information to understand the evolution of sub-aquatic permafrost, which is key to improving and understanding coastal erosion prediction models and potential greenhouse gas emissions. In this study, we use data from 2D electrical resistivity tomography (ERT) collected in the nearshore coastal zone of two Arctic regions that differ in their environmental conditions (e.g., seawater depth and resistivity) to image and study the subsea permafrost. The inversion of 2D ERT data sets is commonly performed using deterministic approaches that favor smoothed solutions, which are typically interpreted using a user-specified resistivity threshold to identify the IBPT position. In contrast, to target the IBPT position directly during inversion, we use a layer-based model parameterization and a global optimization approach to invert our ERT data. This approach results in ensembles of layered 2D model solutions, which we use to identify the IBPT and estimate the resistivity of the unfrozen and frozen sediments, including estimates of uncertainties. Additionally, we globally invert 1D synthetic resistivity data and perform sensitivity analyses to study, in a simpler way, the correlations and influences of our model parameters. The set of methods provided in this study may help to further exploit ERT data collected in such permafrost environments as well as for the design of future field experiments.}, language = {en} } @misc{ArboledaZapataAngelopoulosOverduinetal.2022, author = {Arboleda-Zapata, Mauricio and Angelopoulos, Michael and Overduin, Pier Paul and Grosse, Guido and Jones, Benjamin M. and Tronicke, Jens}, title = {Exploring the capabilities of electrical resistivity tomography to study subsea permafrost}, series = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {1285}, issn = {1866-8372}, doi = {10.25932/publishup-57123}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-571234}, pages = {4423 -- 4445}, year = {2022}, abstract = {Sea level rise and coastal erosion have inundated large areas of Arctic permafrost. Submergence by warm and saline waters increases the rate of inundated permafrost thaw compared to sub-aerial thawing on land. Studying the contact between the unfrozen and frozen sediments below the seabed, also known as the ice-bearing permafrost table (IBPT), provides valuable information to understand the evolution of sub-aquatic permafrost, which is key to improving and understanding coastal erosion prediction models and potential greenhouse gas emissions. In this study, we use data from 2D electrical resistivity tomography (ERT) collected in the nearshore coastal zone of two Arctic regions that differ in their environmental conditions (e.g., seawater depth and resistivity) to image and study the subsea permafrost. The inversion of 2D ERT data sets is commonly performed using deterministic approaches that favor smoothed solutions, which are typically interpreted using a user-specified resistivity threshold to identify the IBPT position. In contrast, to target the IBPT position directly during inversion, we use a layer-based model parameterization and a global optimization approach to invert our ERT data. This approach results in ensembles of layered 2D model solutions, which we use to identify the IBPT and estimate the resistivity of the unfrozen and frozen sediments, including estimates of uncertainties. Additionally, we globally invert 1D synthetic resistivity data and perform sensitivity analyses to study, in a simpler way, the correlations and influences of our model parameters. The set of methods provided in this study may help to further exploit ERT data collected in such permafrost environments as well as for the design of future field experiments.}, language = {en} }