A5019260523- Seismic Imaging and Inversion Techniques
- Seismic Waves and Analysis
- earthquake and tectonic studies
- Seismology and Earthquake Studies
- Geophysical and Geoelectrical Methods
- Electromagnetic Simulation and Numerical Methods
- Geophysical Methods and Applications
- Advanced Numerical Methods in Computational Mathematics
- Electromagnetic Scattering and Analysis
- Earthquake Detection and Analysis
- Geological Modeling and Analysis
- Hydraulic Fracturing and Reservoir Analysis
- Seismic Performance and Analysis
- Drilling and Well Engineering
- Numerical methods for differential equations
- High-pressure geophysics and materials
- Reservoir Engineering and Simulation Methods
- Geophysics and Sensor Technology
- Scientific Computing and Data Management
- Underwater Acoustics Research
- Distributed and Parallel Computing Systems
- Landslides and related hazards
- Model Reduction and Neural Networks
- Soil Moisture and Remote Sensing
- Geological and Geochemical Analysis
Universitat Politècnica de Catalunya
2015-2024
Barcelona Supercomputing Center
2015-2024
CNS Research (United States)
2016
Ludwig-Maximilians-Universität München
2006-2010
Institut de Ciències del Mar
2010
Vietnam Academy of Science and Technology
2010
Lawrence Berkeley National Laboratory
2010
Consejo Superior de Investigaciones Científicas
2009
University of Stuttgart
2008
University of Trento
2008
We present a new numerical method to solve the heterogeneous anelastic, seismic wave equations with arbitrary high order accuracy in space and time on 3-D unstructured tetrahedral meshes. Using velocity–stress formulation provides linear hyperbolic system of source terms that is completed by additional for anelastic functions including strain history material. These result from rheological model generalized Maxwell body permit incorporation realistic attenuation properties viscoelastic...
We present a nodal finite-element method that can be used to compute in parallel highly accurate solutions for 3-D controlled-source electromagnetic forward-modelling problems anisotropic media. Secondary coupled-potential formulation of Maxwell's equations allows avoid the singularities introduced by sources, while completely unstructured tetrahedral meshes and mesh refinement support an representation geological bathymetric complexity improve solution accuracy. Different complex iterative...
Accurate and efficient numerical methods to simulate dynamic earthquake rupture wave propagation in complex media fault geometries are needed address fundamental questions dynamics, integrate seismic geodetic data into emerging approaches for source inversion, generate realistic physics‐based scenarios hazard assessment. Modeling of spontaneous by a high‐order discontinuous Galerkin (DG) method combined with an arbitrarily derivatives (ADER) time integration was introduced two dimensions de...
Finite-difference methods for modeling seismic waves are known to be inaccurate when including a realistic topography, due the large dispersion errors that appear in modelled surface and scattering introduced by staircase approximation topography. As consequence, alternatives finite-difference have been proposed circumvent these issues. We present new numerical scheme 3D elastic wave propagation presence of strong This is based upon staggered grid Lebedev type, or fully (FSG). It uses...
We have developed a new numerical method to solve the heterogeneous poroelastic wave equations in bounded three-dimensional domains. This is discontinuous Galerkin that achieves arbitrary high-order accuracy on unstructured tetrahedral meshes for low-frequency range and inviscid case. By using Biot’s Darcy’s dynamic laws, we built scheme can successfully model propagation fluid-saturated porous media when anisotropy of pore structure allowed. Zero-inflow fluxes are used as absorbing boundary...
We present a new numerical method to solve the heterogeneous elastic anisotropic wave equation with arbitrary high-order accuracy in space and time on unstructured tetrahedral meshes. Using most general Hooke's tensor we derive velocity-stress formulation leading linear hyperbolic system which accounts for variation of material properties depending direction. This approach allows accurate modelling even crystalline symmetry class, triclinic anisotropy, as no interpolation particular mesh...
We introduce the application of an arbitrary high‐order derivative (ADER) discontinuous Galerkin (DG) method to simulate earthquake rupture dynamics. The ADER‐DG uses triangles as computational cells which simplifies process discretization very complex surfaces and volumes by using external automated tools. Discontinuous methods are well suited for solving dynamic problems in velocity‐stress formulation variables naturally at interface between two elements. Therefore, fault has be honored...
We present a new numerical method to solve the heterogeneous anelastic seismic wave equations with arbitrary high order of accuracy in space and time on unstructured triangular tetrahedral meshes two three dimensions, respectively. Using velocity-stress formulation provides linear hyperbolic system source terms that is completed by additional for functions including strain history material. These result from rheological model generalized Maxwell body permit incorporation realistic...
Adjoint methods are a key ingredient of gradient-based full-waveform inversion schemes. While being conceptually elegant, they face the challenge massive memory requirements caused by opposite time directions forward and adjoint simulations necessity to access both wavefields simultaneously for computation sensitivity kernel. To overcome this bottleneck, we have developed lossy compression techniques that significantly reduce with only small computational overhead. Our approach is tailored...
The EU Center of Excellence for Exascale in Solid Earth (ChEESE) develops exascale transition capabilities the domain Earth, an area geophysics rich computational challenges embracing different approaches to (capability, capacity, and urgent computing). first implementation phase project (ChEESE-1P; 2018–2022) addressed scientific technical seismology, tsunami science, volcanology, magnetohydrodynamics, order understand phenomena, anticipate impact natural disasters, contribute risk...
We present an elaborate preconditioning scheme for Krylov subspace methods which has been developed to improve the performance and reduce execution time of parallel node-based finite-element solvers three-dimensional electromagnetic numerical modelling in exploration geophysics.This new preconditioner is based on algebraic multigrid that uses different basic relaxation methods, such as Jacobi, symmetric successive over-relaxation Gauss-Seidel, smoothers wave-front algorithm create groups,...
SUMMARY We present a parallel and high-order Nédélec finite element solution for the marine controlled-source electromagnetic (CSEM) forward problem in 3-D media with isotropic conductivity. Our Python code is implemented on unstructured tetrahedral meshes, which support multiple-scale structures bathymetry general CSEM modelling applications. Based primary/secondary field approach, we solve diffusive form of Maxwell’s equations low-frequency domain. investigate accuracy performance...
We present a routine for 3D magnetotelluric (MT) modeling based upon high-order edge finite element method (HEFEM), tailored and unstructured tetrahedral meshes, high-performance computing (HPC). This implementation extends the PETGEM modeller capabilities, initially developed active-source electromagnetic methods in frequency-domain. assess accuracy, robustness, performance of code using set reference models by MT community well-known reported workshops. The scale geological properties...
We present a quantitative accuracy analysis of the Discontinuous Galerkin Finite-Element method for simulation seismic wave propagation on tetrahedral meshes. Several parameters are responsible results, such as chosen approximation order, spatial discretization, that is, number elements per wavelength, and distance waves due to numerical dispersion dissipation. As error norm we choose time–frequency representation envelope phase misfit seismograms assess resulting since this provides time...
Abstract During the nineteenth and twentieth centuries observational seismologists recorded primarily earthquake-induced translational wave field, while rotational motion still remains poorly observed investigated. We aim to further understand ground its relation with a special emphasis on near few wavelengths away from hypocenter, where damage related might need be considered. A broad picture of available values amplitudes their variability is obtained by gathering most published data...
Southeast Spain experiences relatively low seismicity rates, characterized by slow seismic deformation. However, historical records highlight the significant impact of moderate to large earthquakes on local communities, such as 1518 Vera (Almería) earthquake (Mw 6.4) and 2011 Lorca 5.2). Thus, events pose a considerable risk region, in spite their infrequent occurrence. Given lack comprehensive data this kind events, study contributes towards physics-based hazard model for (SE)...
Urgent Computing (UC) refers to the use of High-Performance (HPC) and Data Analytics (HPDA) Artificial Intelligence (AI) modules during or immediately following emergencies. It typically integrates complex end-to-end workflows with scalable computing resources, where multiple model realizations are necessary account for input uncertainties, all under strict time-to-solution constraints. Enabling urgent HPC in unpredictable events such as earthquakes can significantly enhance resilience...
Earthquakes are among the most frequent yet unpredictable natural hazards, posing substantial risk to human safety and infrastructure globally, particularly, when large-magnitude earthquakes occur. This highlights urgent need develop innovative alternative methodologies for rapidly assessing intensity of ground shaking following an earthquake.This study explores application Machine Learning Estimator Ground Shaking Maps (MLESmap) methodology in New Zealand, a region characterized...
Rotational motions in homogeneous anisotropic elastic media are studied under the assumption of plane wave propagation. The main goal is to investigate influences anisotropy behavior rotational wavefield. focus on P-waves that theoretically do not generate motion isotropic media. By using Kelvin–Christoffel equation, expressions obtained body waves as a function propagation direction and coefficients modulus matrix. As result, amplitudes rotation rates their radiation patterns quantified it...