scholarly journals Applicability of time-lapse refraction seismic tomography for the detection of ground ice degradation

2010 ◽  
Vol 4 (1) ◽  
pp. 77-119 ◽  
Author(s):  
C. Hilbich
Keyword(s):  
Seismic Tomography ◽  
Time Lapse ◽  
Swiss Alps ◽  
P Wave ◽  
Permafrost Degradation ◽  
Travel Times ◽  
Ground Ice ◽  
Refraction Seismic ◽  
Velocity Changes ◽  
Ice Content

Abstract. The ice content of the subsurface is a major factor controlling the natural hazard potential of permafrost degradation in alpine terrain. Monitoring of changes in ground ice content is therefore similarly important as temperature monitoring in mountain permafrost. Although electrical resistivity tomography monitoring (ERTM) has proved to be a valuable tool for the observation of ground ice degradation, results are often ambiguous or contaminated by inversion artefacts. In theory, the P-wave velocity of seismic waves is similarly sensitive to phase changes between unfrozen water and ice. Provided that the general conditions (lithology, stratigraphy, state of weathering, pore space) remain unchanged over the observation period, temporal changes in the observed travel times of repeated seismic measurements should indicate changes in the ice and water content within the pores and fractures of the subsurface material. In this paper, the applicability of refraction seismic tomography monitoring (RSTM) as an independent and complementary method to ERTM is analysed for two test sites in the Swiss Alps. The development and validation of an appropriate RSTM approach involves a) the comparison of time-lapse seismograms and analysis of reproducibility of the seismic signal, b) the analysis of time-lapse travel time curves with respect to shifts in travel times and changes in P-wave velocities, and c) the comparison of inverted tomograms including the quantification of velocity changes. Results show a high potential of the RSTM approach concerning the detection of altered subsurface conditions caused by freezing and thawing processes. For velocity changes on the order of 3000 m/s even an unambiguous identification of significant ground ice loss is possible.

scholarly journals Time-lapse refraction seismic tomography for the detection of ground ice degradation

2010 ◽  
Vol 4 (3) ◽  
pp. 243-259 ◽  
Author(s):  
C. Hilbich
Keyword(s):  
Seismic Tomography ◽  
Time Lapse ◽  
Phase Changes ◽  
Swiss Alps ◽  
P Wave ◽  
Permafrost Degradation ◽  
Travel Times ◽  
Refraction Seismic ◽  
Velocity Changes ◽  
Ice Content

Abstract. The ice content of the subsurface is a major factor controlling the natural hazard potential of permafrost degradation in alpine terrain. Monitoring of changes in ice content is therefore similarly important as temperature monitoring in mountain permafrost. Although electrical resistivity tomography monitoring (ERTM) proved to be a valuable tool for the observation of ice degradation, results are often ambiguous or contaminated by inversion artefacts. In theory, the sensitivity of P-wave velocity of seismic waves to phase changes between unfrozen water and ice is similar to the sensitivity of electric resistivity. Provided that the general conditions (lithology, stratigraphy, state of weathering, pore space) remain unchanged over the observation period, temporal changes in the observed travel times of repeated seismic measurements should indicate changes in the ice and water content within the pores and fractures of the subsurface material. In this paper, a time-lapse refraction seismic tomography (TLST) approach is applied as an independent method to ERTM at two test sites in the Swiss Alps. The approach was tested and validated based on a) the comparison of time-lapse seismograms and analysis of reproducibility of the seismic signal, b) the analysis of time-lapse travel time curves with respect to shifts in travel times and changes in P-wave velocities, and c) the comparison of inverted tomograms including the quantification of velocity changes. Results show a high potential of the TLST approach concerning the detection of altered subsurface conditions caused by freezing and thawing processes. For velocity changes on the order of 3000 m/s even an unambiguous identification of significant ice loss is possible.

Ice or rock matrix? Improved quantitative imaging of Alpine permafrost evolution through time-lapse petrophysical joint inversion

2021 ◽  
Author(s):  
Johanna Klahold ◽  
Christian Hauck ◽  
Florian Wagner
Keyword(s):  
Liquid Water ◽  
Joint Inversion ◽  
Time Lapse ◽  
Swiss Alps ◽  
Permafrost Degradation ◽  
Rock Matrix ◽  
Validation Data ◽  
Water Ice ◽  
Borehole Temperature ◽  
Ice Content

<p>Quantitative estimation of pore fractions filled with liquid water, ice and air is one of the prerequisites in many permafrost studies and forms the basis for a process-based understanding of permafrost and the hazard potential of its degradation in the context of global warming. The volumetric ice content is however difficult to retrieve, since standard borehole temperature monitoring is unable to provide any ice content estimation. Geophysical methods offer opportunities to image distributions of permafrost constituents in a non-invasive manner. A petrophysical joint inversion was recently developed to determine volumetric water, ice, air and rock contents from seismic refraction and electrical resistivity data. This approach benefits from the complementary sensitivities of seismic and electrical data to the phase change between ice and liquid water. A remaining weak point was the unresolved petrophysical ambiguity between ice and rock matrix. Within this study, the petrophysical joint inversion approach is extended along the time axis and respective temporal constraints are introduced. If the porosity (and other time-invariant properties like pore water resistivity or Archie exponents) can be assumed invariant over the considered time period, water, ice and air contents can be estimated together with a temporally constant (but spatially variable) porosity distribution. It is hypothesized that including multiple time steps in the inverse problem increases the ratio of data and parameters and leads to a more accurate distinction between ice and rock content. Based on a synthetic example and a field data set from an Alpine permafrost site (Schilthorn, Swiss Alps) it is demonstrated that the developed time-lapse petrophysical joint inversion provides physically plausible solutions, in particular improved estimates for the volumetric fractions of ice and rock. The field application is evaluated with independent validation data including thaw depths derived from borehole temperature measurements and shows generally good agreement. As opposed to the conventional petrophysical joint inversion, its time-lapse extension succeeds in providing reasonable estimates of permafrost degradation at the Schilthorn monitoring site without <em>a priori </em>constraints on the porosity model.</p>

Upscaling of geophysical measurements: A methodology for the estimation of the total ground ice content at two study sites in the dry Andes of Chile and Argentina

2020 ◽  
Author(s):  
Tamara Mathys ◽  
Christin Hilbich ◽  
Cassandra E.M. Koenig ◽  
Lukas Arenson ◽  
Christian Hauck
Keyword(s):  
Hydrological Cycle ◽  
Spatial Scales ◽  
Geophysical Methods ◽  
Central Andes ◽  
Permafrost Degradation ◽  
Rock Glacier ◽  
Ground Ice ◽  
Refraction Seismic ◽  
Geophysical Surveys ◽  
Ice Content

<p>With climate change and the associated continuing recession of glaciers, water security, especially in regions depending on the water supply from glaciers, is threatened. In this context, the understanding of permafrost distribution and its degradation is of increasing importance as it is currently debated whether ground ice can be considered as a significant water reservoir and as an alternative resource of fresh water that could potentially moderate water scarcity during dry seasons in the future. Thus, there is a pressing need to better understand how much water is stored as ground ice in areas with extensive permafrost occurrence and how meltwater from permafrost degradation may contribute to the hydrological cycle in the region.</p><p>Although permafrost and permafrost landforms in the Central Andes are considered to be abundant and well developed, the data is scarce and understanding of the Andean cryosphere lacking, especially in areas devoid of glaciers and rock glaciers.</p><p>In the absence of boreholes and test pits, geophysical investigations are a feasible and cost-effective technique to detect ground ice occurrences within a variety of landforms and substrates. In addition to the geophysical surveys themselves, upscaling techniques are needed to estimate ground ice content, and thereby future water resources, on larger spatial scales. To contribute to reducing the data scarcity regarding ground ice content in the Central Andes, this study focuses on the permafrost distribution and the ground ice content (and its water equivalent) of two catchments in the semi-arid Andes of Chile and Argentina. Geophysical methods (Electrical Resistivity Tomography, ERT and Refraction Seismic Tomography, RST) were used to detect and quantify ground ice in the study regions in the framework of environmental impact assessments in mining areas. Where available, ERT and RST measurements were quantitatively combined to estimate the volumetric ground ice content using the Four Phase Model (Hauck et al., 2011). Furthermore, we developed one of the first methodologies for the upscaling of these geophysical-based ground ice quantifications to an entire catchment in order to estimate the total ground ice volume in the study areas.</p><p>In this contribution we will present the geophysical data, the upscaling methodology used to estimate total ground ice content (and water equivalent) of permafrost areas, and some first estimates of total ground ice content in rock glacier and rock glacier free areas and compare them to conventional estimates using remotely sensed data.</p><p> </p><p>Hauck, C., Böttcher, M., and Maurer, H. (2011). A new model for estimating subsurface ice content based on combined electrical and seismic datasets, The Cryosphere, 5: 453-468.</p>

scholarly journals Resolution capacity of geophysical monitoring regarding permafrost degradation induced by hydrological processes

2017 ◽  
Vol 11 (6) ◽  
pp. 2957-2974 ◽  
Author(s):  
Benjamin Mewes ◽  
Christin Hilbich ◽  
Reynald Delaloye ◽  
Christian Hauck
Keyword(s):  
Electrical Resistivity ◽  
Seismic Tomography ◽  
Field Data ◽  
Electrical Resistivity Tomography ◽  
Geophysical Methods ◽  
Phase Model ◽  
Permafrost Degradation ◽  
Resistivity Tomography ◽  
Refraction Seismic ◽  
Ice Content

Abstract. Geophysical methods are often used to characterize and monitor the subsurface composition of permafrost. The resolution capacity of standard methods, i.e. electrical resistivity tomography and refraction seismic tomography, depends not only on static parameters such as measurement geometry, but also on the temporal variability in the contrast of the geophysical target variables (electrical resistivity and P-wave velocity). Our study analyses the resolution capacity of electrical resistivity tomography and refraction seismic tomography for typical processes in the context of permafrost degradation using synthetic and field data sets of mountain permafrost terrain. In addition, we tested the resolution capacity of a petrophysically based quantitative combination of both methods, the so-called 4-phase model, and through this analysed the expected changes in water and ice content upon permafrost thaw. The results from the synthetic data experiments suggest a higher sensitivity regarding an increase in water content compared to a decrease in ice content. A potentially larger uncertainty originates from the individual geophysical methods than from the combined evaluation with the 4-phase model. In the latter, a loss of ground ice can be detected quite reliably, whereas artefacts occur in the case of increased horizontal or vertical water flow. Analysis of field data from a well-investigated rock glacier in the Swiss Alps successfully visualized the seasonal ice loss in summer and the complex spatially variable ice, water and air content changes in an interannual comparison.

scholarly journals Towards accurate quantification of ice content in permafrost of the Central Andes – Part II: an upscaling strategy of geophysical measurements to the catchment scale at two study sites

2021 ◽  
Author(s):  
Tamara Mathys ◽  
Christin Hilbich ◽  
Lukas U. Arenson ◽  
Pablo A. Wainstein ◽  
Christian Hauck
Keyword(s):  
Geophysical Data ◽  
Distribution Model ◽  
Permafrost Degradation ◽  
Rock Glacier ◽  
Ground Ice ◽  
Data Set ◽  
Refraction Seismic ◽  
Geophysical Surveys ◽  
Rock Glaciers ◽  
Ice Content

Abstract. With ongoing climate change, there is a pressing need to better understand how much water is stored as ground ice in areas with extensive permafrost occurrence and how the regional water balance may alter in response to the potential generation of melt water from permafrost degradation. However, field-based data on permafrost in remote and mountainous areas such as the South-American Andes is scarce and most current ground ice estimates are based on broadly generalised assumptions such as volume-area scaling and mean ground ice content estimates of rock glaciers. In addition, ground ice contents in permafrost areas outside of rock glaciers are usually not considered, resulting in a significant uncertainty regarding the volume of ground ice in the Andes, and its hydrological role. In part I of this contribution, Hilbich et al. (submitted) present an extensive geophysical data set based on Electrical Resistivity Tomography (ERT) and Refraction Seismic Tomography (RST) surveys to detect and quantify ground ice of different landforms and surface types in several study regions in the semi-arid Andes of Chile and Argentina with the aim to contribute to the reduction of this data scarcity. In part II we focus on the development of a methodology for the upscaling of geophysical-based ground ice quantification to an entire catchment to estimate the total ground ice volume (and its estimated water equivalent) in the study areas. In addition to the geophysical data, the upscaling approach is based on a permafrost distribution model and classifications of surface and landform types. Where available, ERT and RST measurements were quantitatively combined to estimate the volumetric ground ice content using petrophysical relationships within the Four Phase Model (Hauck et al., 2011). In addition to introducing our upscaling methodology, we demonstrate that the estimation of large-scale ground ice volumes can be improved by including (i) non-rock glacier permafrost occurrences, and (ii) field evidence through a large number of geophysical surveys and ground truthing information. The results of our study indicate, that (i) conventional ground ice estimates for rock-glacier dominated catchments without in-situ data may significantly overestimate ground ice contents, and (ii) substantial volumes of ground ice may also be present in catchments where rock glaciers are lacking.

scholarly journals Resolution capacity of geophysical monitoring regarding permafrost degradation induced by hydrological processes

2016 ◽  
Author(s):  
Benjamin Mewes ◽  
Christin Hilbich ◽  
Reynald Delaloye ◽  
Christian Hauck
Keyword(s):  
Electrical Resistivity ◽  
Seismic Tomography ◽  
Field Data ◽  
Electrical Resistivity Tomography ◽  
Geophysical Methods ◽  
Phase Model ◽  
Permafrost Degradation ◽  
Resistivity Tomography ◽  
Refraction Seismic ◽  
Ice Content

Abstract. Geophysical methods are often used to characterise and monitor the subsurface composition of permafrost. The resolution capacity of standard methods, i.e. Electrical Resistivity Tomography and Refraction Seismic Tomography, depends hereby not only on static parameters such as measurement geometry, but also on the temporal variability in the contrast of the geophysical variables (electrical resistivity and P-wave velocity). Our study analyses the resolution capacity of Electrical Resistivity Tomography and Refraction Seismic Tomography for typical processes in the context of permafrost degradation using synthetic and field data sets of mountain permafrost terrain. In addition, we tested especially the resolution capacity of a petrophysically-based quantitative combination of both methods, the so-called 4-phase model, and by this analysed the expected changes in water and ice content upon permafrost thaw. The results from the synthetic data experiments suggest a higher sensitivity regarding increasing water content compared to decreased ice content, and potentially larger uncertainty for the individual geophysical methods than for the combined evaluation with the 4-phase model. In the latter, ground ice loss can be detected quite reliably, whereas artefacts occur in the case of increased horizontal or vertical water flow. Analysis of field data from a well-investigated rock glacier in the Swiss Alps successfully visualised the seasonal ice loss in summer, and the complex spatially variable ice-, water- and air content changes in an interannual comparison.

P-wave travel times from shallow earthquakes recorded in India and inferred upper mantle structure

1968 ◽  
Vol 58 (6) ◽  
pp. 1879-1897
Author(s):  
K. L. Kaila ◽  
P. R. Reddy ◽  
Hari Narain
Keyword(s):  
Upper Mantle ◽  
P Wave ◽  
Mantle Structure ◽  
Travel Times ◽  
The North ◽  
Shallow Earthquakes ◽  
Upper Mantle Structure ◽  
Wave Travel ◽  
Excess Value ◽  
Velocity Changes

ABSTRACT P-wave travel times of 39 shallow earthquakes and three nuclear explosions with epicenters in the North in Himalayas, Tibet, China and USSR as recorded in Indian observatories have been analyzed statistically by the method of weighting observations. The travel times from Δ = 2° to 50° can be represented by four straight line segments indicating abrupt velocity changes around 19°, 22° and 33° respectively. The P-wave velocity at the top of the mantle has been found to be 8.31 ± 0.02 km/sec. Inferred upper mantle structure reveals three velocity discontinuities in the upper mantle at depths (below the crust) of 380 ± 20, 580 ± 50 and 1000 ± 120 km with velocities below the discontinuities as 9.47 ± 0.06, 10.15 ± 0.07 and 11.40 ± 0.08 km/sec respectively. The J-B residuals up to Δ = 19° are mostly negative varying from 1 to 10 seconds with a dependence on Δ values indicating a different upper mantle velocity in the Himalayan region as compared to that used by Jeffreys-Bullen in their tables (1940). Between 19° to 33° there is a reasonably good agreement between the J-B curve and the observation points. From Δ = 33° to 50° the J-B residuals are mostly positive with an average excess value of about 4 sec.

Post-glacial dynamics of alpine Little Ice Age glacitectonized frozen landforms (Swiss Alps)

2020 ◽  
Author(s):  
Julie Wee ◽  
Reynald Delaloye ◽  
Chloé Barboux
Keyword(s):  
Little Ice Age ◽  
Spatial Configuration ◽  
Sar Interferometry ◽  
Swiss Alps ◽  
Ice Age ◽  
Permafrost Degradation ◽  
Mass Wasting ◽  
Ground Ice ◽  
Slope Movements ◽  
Glacier Ice

<p>Glaciers and frozen debris landforms have coexisted and episodically interacted throughout the Holocene, the former having altered the development, spatial distribution and thermal regime of the latter. In the Alps, the apogee of last interaction phase occurred during the Little Ice Age (LIA). Since then, due to glacier shrinkage, interactions between glaciers and LIA pre-existing frozen debris have gradually diminished and are leaning towards being non-existent. Post-LIA glacier forefields in permafrost environments, including associated glacitectonized frozen landforms (GFL) have shifted from a thermal and mechanical glacier dominant regime towards a periglacial or even post-periglacial regime. GFL are undergoing thermal and mechanical readjustments in response to both the longer-term glacier recession and the more recent drastic climatic warming. They can be expressed by a combination of mass-wasting processes and thaw-induced subsidence.</p><p> </p><p>In various regions of the Swiss Alps, slope movements occurring in a periglacial context have been inventoried in previous works using differential SAR interferometry (DInSAR) (Barboux et al., 2014). In the scope of this study, and focusing solely on mass-wasting GFL, the former inventory allowed the identification of the latter under various spatial configurations within LIA glacier forefields. While most observed GFL are disconnected from the associated glacier, some are still connected. Additionally, ground ice occurs as interstitial or massive (buried) glacier ice. This potentially infers the ongoing of non-uniform morphodynamical readjustments.</p><p> </p><p>To understand the site-specific behaviour of GFL, the analysis of long-term time-series of permafrost monitoring and multi-temporal high-resolution Digital Elevation Models will allow the assessment of the recent evolution of the Aget and Ritord/Challand LIA glacier forefields (46°00’32’’ N, 7°14’20’’ E and 45°57’10’’ N, 7°14’52’’ E, respectively) and their associated GFL (i.e. push-moraines). Both glacier forefields present a contrasting spatial configuration, making their morphodynamical evolution to differ partly from one another. The Aget push-moraine is a back-creeping GFL, which has been disconnected from the Aget glacier since the 1940s at latest. For the last two decades, surface displacement velocities have decelerated in comparison to the accelerating regional trend (PERMOS, 2019). Additionally, a 30% decrease of the electrical resistivity of the frozen ground, combined with locally observed thaw-induced subsidence of up to 10 cm/year suggest an advanced permafrost degradation. The Ritord/Challand system presents a push-moraine disconnected from its glacier as well as several push-moraines connected to a still existing debris-covered glacier. Between 2016 and 2019, surface lowering up to 10 m attesting massive ice melt has been locally detected in the former where buried glacier ice was visually observed. Whereas in the latter, subtle surface displacements ranging from 10 to 30 cm/year occur. This confirms the heterogeneity of the morphodynamical processes occurring in GFL, expressed as a function of both their spatial configuration and ground ice properties.</p><p> </p><p>Barboux, C., Delaloye R. and Lambiel, C. (2014). Inventorying slope movements in an Alpine environment using DInSAR. Earth Surface Processes and Landforms, 39/15, 2087-2099.</p><p>PERMOS 2019. Permafrost in Switzerland 2014/2015 to 2017/2018. Noetzli, J., Pellet, C., and Staub, B. (eds.), Glaciological Report (Permafrost) No. 16-19 of the Cryospheric Commission of the Swiss Academy of Sciences, 104.</p>

Toward a full 4D seismic tomography: a case study of an active mine

2021 ◽  
Author(s):  
Nicola Piana Agostinetti ◽  
Christina Dahnér-Lindkvist ◽  
Savka Dineva
Keyword(s):  
Monte Carlo ◽  
Seismic Tomography ◽  
Water Injection ◽  
Elastic Field ◽  
Monte Carlo Sampling ◽  
P Wave ◽  
Travel Times ◽  
S Wave ◽  
Posterior Probability Distribution ◽  
Space And Time

<p>Rock elasticity in the subsurface can change in response to natural phenomena (e.g. massive precipitation, magmatic processes) and human activities (e.g. water injection in geothermal wells, ore-body exploitation). However, understanding and monitoring the evolution of physical properties of the crust is a challenging due to the limited possibility of reaching such depths and making direct measurements of the state of the rocks. Indirect measurements, like seismic tomography, can give some insights, but are generally biased by the un-even distribution (in space and time) of the information collected from seismic observations (travel-times and/or waveforms). Here we apply a Bayesian approach to overcome such limitations, so that data uncertainties and data distribution are fully accounted in the reconstruction of the posterior probability distribution of the rock elasticity  We compute a full 4D local earthquake tomography based on trans-dimensional Markov chain Monte Carlo sampling of 4D elastic models, where the resolution in space and time is fully data-driven. To test our workflow, we make use of a “controlled laboratory”: we record seismic data during one month of mining activities across a 800x700x600 m volume of Kiruna mine (LKAB, Sweden). During such period, we obtain about 260 000 P-wave and 240 000 S-wave travel-times coming from about 36000  seismic events. We operate a preliminary selection of the well-located events, using a Monte Carlo search. Arrival-times of about 19 000 best-located events (location errors less than 20m) are used as input to the tomography workflow. Preliminary results indicate that: (1) short-term (few hours) evolutions of the elastic field are mainly driven by seismic activation, i.e. the occurrence of a seismic swarm, close to the mine ore-passes. Such phenomena partially mask the effects of explosions; (2) long-term (2-3 days) evolutions of the elastic field closely match the local measurements of the stress field at a colocated stress cell. </p>

Sign in / Sign up

Export Citation Format

Share Document