scholarly journals 3D analytical and numerical modelling of the regional topography influence on the surface deformation due to underground heat source

2011 ◽  
Vol 41 (3) ◽  
pp. 251-265
Author(s):  
Ladislav Brimich ◽  
María Charco ◽  
Igor Kohút ◽  
José Fernández
Keyword(s):  
Numerical Model ◽  
Numerical Modelling ◽  
Heat Source ◽  
Surface Deformation ◽  
Ground Deformation ◽  
Source Parameters ◽  
Viscoelastic Model ◽  
Thermal Source ◽  
Topographic Effects ◽  
Regional Topography

3D analytical and numerical modelling of the regional topography influence on the surface deformation due to underground heat source Thermo-elastic strains and stresses play a considerable role in the stress state of the lithosphere and its dynamics, especially at pronounced positive geothermal anomalies. Topography has a significant effect on ground deformation. In this paper we describe two methods for including the topographic effects in the thermo-viscoelastic model. First we use an approximate methodology which assumes that the main effect of the topography is due to distance from the source to the free surface and permits to have an analytical solution very attractive for solving the inverse problem. A numerical solution using Finite Element Method (FEM) is also computed. The numerical method allows to include the local shape of the topography in the modelling. In the numerical model the buried magmatic body is represented by a finite volume thermal source. The temperature distribution is computed by the higher-degree FEM. For analytical as well as numerical model solution only the forces of thermal origin are considered. The comparison of the results obtained using both analytical and numerical techniques shows the qualitative agreement of the vertical displacements. In the numerical values small differences were obtained. The results show that for the volcanic areas with an important relief the perturbation of the thermo-viscoelastic solution (deformation and total gravity anomaly) due to the topography can be quite significant. In consequence, neglecting topography could give erroneous results in the estimated source parameters.

scholarly journals Surface displacements, deformations and gravity changes due to underground heat source

2018 ◽  
Vol 48 (3) ◽  
pp. 271-279
Author(s):  
Ladislav Brimich ◽  
Igor Kohút
Keyword(s):  
Numerical Model ◽  
Ground Deformation ◽  
Source Parameters ◽  
Viscoelastic Model ◽  
Thermal Source ◽  
Topographic Effects ◽  
Surface Displacements ◽  
Magmatic Body ◽  
Thermal Origin ◽  
Plain Strain

Abstract Thermo-elastic strains and stresses play a considerable role in the stress state of the lithosphere and its dynamics, especially at pronounced positive geothermal anomalies. Topography has a significant effect on ground deformation. Two methods for including the topographic effects in the thermo-viscoelastic model are described. First we use an approximate methodology which assumes that the main effect of topography is due to the distance from the source to the free surface and permits to have an analytical solution very attractive for solving the inverse problem. A numerical solution (for 2D plain strain case) is also computed using finite element method (FEM). The numerical method allows to include the local shape of the topography in the modeling. In the numerical model the buried magmatic body is represented by a finite volume thermal source. The temperature distribution is computed by the higher-degree FEM. For analytical as well as numerical model solution only the forces of thermal origin are considered. The results show that for the volcanic areas with prominent topography, its effect on the perturbation of the thermo-viscoelastic solution (deformation and total gravity anomaly) can be quite significant. In consequence, neglecting the topography could give erroneous results in the estimated source parameters.

scholarly journals Ground Deformation and Source Geometry of the 30 October 2016 Mw 6.5 Norcia Earthquake (Central Italy) Investigated Through Seismological Data, DInSAR Measurements, and Numerical Modelling

2018 ◽  
Vol 10 (12) ◽  
pp. 1901 ◽  
Author(s):  
Emanuela Valerio ◽  
Pietro Tizzani ◽  
Eugenio Carminati ◽  
Carlo Doglioni ◽  
Susi Pepe ◽  
...  
Keyword(s):  
Numerical Modelling ◽  
Surface Deformation ◽  
Fault Plane ◽  
Ground Deformation ◽  
Central Italy ◽  
Vertical Component ◽  
Deformation Pattern ◽  
Distribution Analysis ◽  
Seismological Data ◽  
Modelling Approach

We investigate the Mw 6.5 Norcia (Central Italy) earthquake by exploiting seismological data, DInSAR measurements, and a numerical modelling approach. In particular, we first retrieve the vertical component (uplift and subsidence) of the displacements affecting the hangingwall and the footwall blocks of the seismogenic faults identified, at depth, through the hypocenters distribution analysis. To do this, we combine the DInSAR measurements obtained from coseismic SAR data pairs collected by the ALOS-2 sensor from ascending and descending orbits. The achieved vertical deformation map displays three main deformation patterns: (i) a major subsidence that reaches the maximum value of about 98 cm near the epicentral zones nearby the town of Norcia; (ii) two smaller uplift lobes that affect both the hangingwall (reaching maximum values of about 14 cm) and the footwall blocks (reaching maximum values of about 10 cm). Starting from this evidence, we compute the rock volumes affected by uplift and subsidence phenomena, highlighting that those involved by the retrieved subsidence are characterized by significantly higher deformation values than those affected by uplift (about 14 times). In order to provide a possible interpretation of this volumetric asymmetry, we extend our analysis by applying a 2D numerical modelling approach based on the finite element method, implemented in a structural-mechanic framework, and exploiting the available geological and seismological data, and the ground deformation measurements retrieved from the multi-orbit ALOS-2 DInSAR analysis. In this case, we consider two different scenarios: the first one based on a single SW-dipping fault, the latter on a main SW-dipping fault and an antithetic zone. In this context, the model characterized by the occurrence of an antithetic zone presents the retrieved best fit coseismic surface deformation pattern. This result allows us to interpret the subsidence and uplift phenomena caused by the Mw 6.5 Norcia earthquake as the result of the gravitational sliding of the hangingwall along the main fault plane and the frictional force acting in the opposite direction, consistently with the double couple fault plane mechanism.

scholarly journals Optimisation on Thermoforming of Biodegradable Poly (Lactic Acid) (PLA) by Numerical Modelling

Polymers ◽  
2021 ◽  
Vol 13 (4) ◽  
pp. 654
Author(s):  
Huidong Wei
Keyword(s):  
Lactic Acid ◽  
Numerical Modelling ◽  
Temperature Profile ◽  
Food Packaging ◽  
Strong Dependence ◽  
Viscoelastic Model ◽  
Poly Lactic Acid ◽  
Strain History ◽  
Optimal Temperature ◽  
Broad Perspective

Poly (lactic acid) (PLA) has a broad perspective for manufacturing green thermoplastic products by thermoforming for its biodegradable properties. The mechanical behaviour of PLA has been demonstrated by its strong dependence on temperature and strain rate at biaxial deformation. A nonlinear viscoelastic model by the previous study was employed in a thermoforming process used for food packaging. An optimisation approach was developed by achieving the optimal temperature profile of specimens by defining multiple heating zones based on numerical modelling with finite element analysis (FEA). The forming process of a PLA product was illustrated by modelling results on shape evolution and biaxial strain history. The optimal temperature profile was suggested in scalloped zones to achieve more even thickness distribution. The sensitivity of the optimal results was addressed by checking the robustness under perturbation.

scholarly journals Fault source investigation of the 6 December 2016 MwMw 6.5 Pidie Jaya, Indonesia, earthquake based on GPS and its implications of the geological survey result

2020 ◽  
Vol 14 (4) ◽  
pp. 405-412
Author(s):  
Endra Gunawan ◽  
Takuya Nishimura ◽  
Susilo Susilo ◽  
Sri Widiyantoro ◽  
Nanang T. Puspito ◽  
...  
Keyword(s):  
Surface Deformation ◽  
Source Parameters ◽  
Secondary Effects ◽  
Near Surface ◽  
Morphological Attributes ◽  
Ground Shaking ◽  
The North ◽  
Coulomb Failure ◽  
North And South ◽  
Fault Source

AbstractOn 6 December 2016 at 22:03 UTC, a devastating magnitude 6-class strike-slip earthquake occurred along an unidentified and unmapped fault in Pidie Jaya, northern Sumatra. We analysed the possible fault using continuous Global Positioning System (GPS) observation available in the region. In our investigation, we searched for the fault source parameters of the north- and south-dipping left-lateral faults and the west- and east-dipping right-lateral faults. We identified that the fault responsible for the earthquake was located offshore, with a southwest-northeast direction. We also computed the Coulomb failure stress and compared the result with the distribution of the aftershocks. In this study, we demonstrated that the result of the geological field survey conducted soon after the mainshock was attributed to the secondary effects of ground shaking and near-surface deformation, and not surface faulting. The newly identified offshore fault proposed by this study calls for further investigation of the corresponding submarine morphological attributes in this particular region.

scholarly journals Three-dimensional deformation time series of glacier motion from multiple-aperture DInSAR observation

2019 ◽  
Vol 93 (12) ◽  
pp. 2651-2660 ◽  
Author(s):  
Sergey Samsonov
Keyword(s):  
Time Series ◽  
Surface Deformation ◽  
Ground Deformation ◽  
Three Dimensional ◽  
Data Sets ◽  
Ice Flow ◽  
The North ◽  
Glacier Ice ◽  
Deformation Component ◽  
3D Deformation

AbstractThe previously presented Multidimensional Small Baseline Subset (MSBAS-2D) technique computes two-dimensional (2D), east and vertical, ground deformation time series from two or more ascending and descending Differential Interferometric Synthetic Aperture Radar (DInSAR) data sets by assuming that the contribution of the north deformation component is negligible. DInSAR data sets can be acquired with different temporal and spatial resolutions, viewing geometries and wavelengths. The MSBAS-2D technique has previously been used for mapping deformation due to mining, urban development, carbon sequestration, permafrost aggradation and pingo growth, and volcanic activities. In the case of glacier ice flow, the north deformation component is often too large to be negligible. Historically, the surface-parallel flow (SPF) constraint was used to compute the static three-dimensional (3D) velocity field at various glaciers. A novel MSBAS-3D technique has been developed for computing 3D deformation time series where the SPF constraint is utilized. This technique is used for mapping 3D deformation at the Barnes Ice Cap, Baffin Island, Nunavut, Canada, during January–March 2015, and the MSBAS-2D and MSBAS-3D solutions are compared. The MSBAS-3D technique can be used for studying glacier ice flow at other glaciers and other surface deformation processes with large north deformation component, such as landslides. The software implementation of MSBAS-3D technique can be downloaded from http://insar.ca/.

Slip Distribution of the 2020 Mw6.9 Samos Earthquake Using a Bayesian Approach

2021 ◽  
Author(s):  
Figen Eskikoy ◽  
Semih Ergintav ◽  
Uğur Dogan ◽  
Seda Özarpacı ◽  
Alpay Özdemir ◽  
...  
Keyword(s):  
Surface Deformation ◽  
Fault Plane ◽  
Source Parameters ◽  
Geodetic Data ◽  
Slip Distribution ◽  
Double Difference ◽  
Gps Data ◽  
Analysis Tool ◽  
Finite Fault ◽  
Continuous Gps

<p>On 2020 October 30, an M<sub>w</sub>6.9 earthquake struck offshore Samos Island. Severe structural damages were observed in Greek Islands and city of Izmir (Turkey). 114 people lost their lives and more than a thousand people were injured in Turkey. The earthquake triggered local tsunami. Significant seismic activity occurred in this region following the earthquake and ~1800 aftershocks (M>1) were recorded by KOERI within the first three days. In this study, we analyze the slip distribution and aftershocks of the 2020 earthquake.</p><p>For the aftershock relocations, the continuous waveforms were collected from NOA, Disaster and Emergency Management Authority of Turkey (AFAD) and KOERI networks. The database   was created based on merged catalogs from AFAD and KOERI. For estimating optimized aftershock location distribution, the P and S phases of the aftershocks are picked manually and relocated with double difference algorithm. In addition, source mechanisms of aftershocks M>4 are obtained from regional body and surface waveforms.</p><p>The surface deformation of the earthquake was obtained from both descending and ascending orbits of the Sentinel-1 A/B and ALOS2 satellites. Since the rupture zone is beneath the Gulf of Kusadası, earthquake related deformation in the interferograms can only be observed on the northern part of the Samos Island. We processed all possible pairs chose the image pairs with the lowest noise level.</p><p>In this study, we used 25 continuous GPS stations which are compiled from TUSAGA-Aktif in Turkey and NOANET in Greece. In addition to continuous GPS data, on 2020 November 1, GPS survey was initiated and the earthquake deformation was measured on 10 GNSS campaign sites (TUTGA), along onshore of Turkey.</p><p>The aim of this study is to estimate the spatial and temporal rupture evolution of the earthquake from geodetic data jointly with near field displacement waveforms. To do so, we use the Bayesian Earthquake Analysis Tool (BEAT).</p><p>As a first step of the study, rectangular source parameters were estimated by using GPS data. In order to estimate the slip distribution, we used both ascending and descending tracks of Sentinel-1 data, ALOS2 and GPS displacements. In our preliminary geodetic data based finite fault model, we used the results of focal mechanism and GPS data inversion solutions for the initial fault plane parameters. The slip distribution results indicate that earthquake rupture is ~35 km long and the maximum slip is ~2 m normal slip along a north dipping fault plane. This EW trending, ~45° north dipping normal faulting system consistent with this tectonic regime in the region. This seismically active area is part of a N-S extensional regime and controlled primarily by normal fault systems.</p><p><strong>Acknowledgements</strong></p><p>This work is supported by the Turkish Directorate of Strategy and Budget under the TAM Project number 2007K12-873.</p>

3D Visual Technology of Geo-Deformation Disasters Induced by Mining Subsidence Based on ArcGIS Engine

2012 ◽  
Vol 500 ◽  
pp. 428-436 ◽  
Author(s):  
Ke Ming Yang ◽  
Jun Ting Ma ◽  
Bo Pang ◽  
Yi Bin Wang ◽  
Ran Wang ◽  
...  
Keyword(s):  
Coal Mining ◽  
Surface Deformation ◽  
Ground Deformation ◽  
Ground Surface ◽  
Horizontal Displacement ◽  
Mining Subsidence ◽  
Underground Coal Mining ◽  
Vertical Movements ◽  
Probability Integral Method ◽  
Arcgis Engine

Mining subsidence often produces significant horizontal and vertical movements at the ground surface, the surface deformation induced by underground coal mining can be predicted by probability integral method, and the surface geo-deformation disasters can be visualized based on GIS components. A three dimensional (3D) visualizing system of surface geo-deformation information is designed and developed with ArcGIS Engine and C# in the study. According to the surface deformation-predicted data induced by underground coal mining in Guobei Coalmine of Huaibei mine field, the extents and degrees of ground deformation disasters are visualized in 3D views for surface vertical subsidence, slope, curvature, horizontal displacement and horizontal strain based on the GIS-developed application platform.

scholarly journals Copernicus Sentinel-1 MT-InSAR, GNSS and Seismic Monitoring of Deformation Patterns and Trends at the Methana Volcano, Greece

2020 ◽  
Vol 10 (18) ◽  
pp. 6445 ◽  
Author(s):  
Theodoros Gatsios ◽  
Francesca Cigna ◽  
Deodato Tapete ◽  
Vassilis Sakkas ◽  
Kyriaki Pavlou ◽  
...  
Keyword(s):  
Land Surface ◽  
Surface Deformation ◽  
Ground Deformation ◽  
Local Population ◽  
Seasonal Fluctuation ◽  
Satellite System ◽  
Deformation Monitoring ◽  
Persistent Scatterers ◽  
Geothermal Activity ◽  
Global Navigation Satellite

The Methana volcano in Greece belongs to the western part of the Hellenic Volcanic Arc, where the African and Eurasian tectonic plates converge at a rate of approximately 3 cm/year. While volcanic hazard in Methana is considered low, the neotectonic basin constituting the Saronic Gulf area is seismically active and there is evidence of local geothermal activity. Monitoring is therefore crucial to characterize any activity at the volcano that could impact the local population. This study aims to detect surface deformation in the whole Methana peninsula based on a long stack of 99 Sentinel-1 C-band Synthetic Aperture Radar (SAR) images in interferometric wide swath mode acquired in March 2015–August 2019. A Multi-Temporal Interferometric SAR (MT-InSAR) processing approach is exploited using the Interferometric Point Target Analysis (IPTA) method, involving the extraction of a network of targets including both Persistent Scatterers (PS) and Distributed Scatterers (DS) to augment the monitoring capability across the varied land cover of the peninsula. Satellite geodetic data from 2006–2019 Global Positioning System (GPS) benchmark surveying are used to calibrate and validate the MT-InSAR results. Deformation monitoring records from permanent Global Navigation Satellite System (GNSS) stations, two of which were installed within the peninsula in 2004 (METH) and 2019 (MTNA), are also exploited for interpretation of the regional deformation scenario. Geological, topographic, and 2006–2019 seismological data enable better understanding of the ground deformation observed. Line-of-sight displacement velocities of the over 4700 PS and 6200 DS within the peninsula are from −18.1 to +7.5 mm/year. The MT-InSAR data suggest a complex displacement pattern across the volcano edifice, including local-scale land surface processes. In Methana town, ground stability is found on volcanoclasts and limestone for the majority of the urban area footprint while some deformation is observed in the suburban zones. At the Mavri Petra andesitic dome, time series of the exceptionally dense PS/DS network across blocks of agglomerate and cinder reveal seasonal fluctuation (5 mm amplitude) overlapping the long-term stable trend. Given the steepness of the slopes along the eastern flank of the volcano, displacement patterns may indicate mass movements. The GNSS, seismological and MT-InSAR analyses lead to a first account of deformation processes and their temporal evolution over the last years for Methana, thus providing initial information to feed into the volcano baseline hazard assessment and monitoring system.

scholarly journals THE USE OF ARRAY PROCESSORS FOR NUMERICAL MODELLING OF TIDAL ESTUARY DYNAMICS

1980 ◽  
Vol 1 (17) ◽  
pp. 142
Author(s):  
D. Prandle ◽  
E.R. Funke ◽  
N.L. Crookshank ◽  
R. Renner
Keyword(s):  
Numerical Model ◽  
Numerical Modelling ◽  
Hybrid Model ◽  
Numerical Models ◽  
Computer Time ◽  
Scale Model ◽  
Hybrid Modelling ◽  
Tidal Propagation ◽  
Physical Scale ◽  
Array Processors

The use of array processors for the numerical modelling of estuarine systems is discussed here in the context of "hybrid modelling", however, it is shown that array processors may be used to advantage in independent numerical simulations. Hybrid modelling of tidal estuaries was first introduced by fiolz (1977) and later by Funke and Crookshank (1978). In a hybrid model, tidal propagation in an estuary is simulated by dynamically linking an hydraulic (or physical) scale model of part of the estuary to a numerical model of the remaining part in a manner such that a free interchange of flow occurs at the interface(s). Typically, the elevation of the water surface at the boundary of the scale model is measured and transmitted to the numerical model. In return, the flow computed at the boundary of the numerical model is fed directly into the scale model. This approach enables the extent of the scale model to be limited to the area of immediate interest (or to that area where flow conditions are such that they can be most accurately simulated by a scale model). In addition, since the region simulated by the numerical model can be extended almost indefinitely, the problems of spurious reflections from downstream boundaries can be eliminated. In normal use, numerical models are evaluated on the basis of computing requirements, cost and accuracy. The computer time required to simulate one tide cycle is, in itself, seldom of interest except in so far as it affects the above criteria. However in hybrid modelling this parameter is often paramount since concurrent operation of the numerical and scale models requires that the former must keep pace with the latter. The earlier hybrid model of the St. Lawrence (Funke and Crookshank, 1978) involved a one-dimensional numerical model of the upstream regions of the river. However, future applications are likely to involve extensive two-dimensional numerical simulation.

Sign in / Sign up

Export Citation Format

Share Document