Difference between revisions of "JSG T.25"

From Icctwiki
Jump to: navigation, search
m (Pendl moved page JSG0.12 to JSG T.25)
(Introduction)
Line 8: Line 8:
 
===Introduction===
 
===Introduction===
  
Efficient numerical methods and HPC (high performance computing) facilities provide new opportunities in many applications in geodesy. The goal of the JSG is to apply numerical methods and/or HPC techniques mostly for gravity field modelling and nonlinear filtering of various geodetic data. The discretization numerical methods like the finite element method (FEM), finite volume method (FVM) and boundary element method (BEM) or the meshless methods like the method of fundamental solutions (MFS) or singular boundary method (SOR) can be efficiently used to solve the geodetic boundary value problems and nonlinear diffusion filtering, or to process e.g. the GOCE observations. Their parallel implementations and large-scale parallel computations on clusters with distributed memory using the MPI (Message Passing Interface) standards allows to solve such problems in spatial domains while obtaining high-resolution numerical solutions.  
+
The seismic tomography is primarily used to provide images of the Earth’s inner structure based on the analysis of seismic waves due to earthquakes and (controlled) explosions. This technique involves several different methods for processing P-, S- and surface waves on the principle of solving inverse problems for finding locations of reflection and refraction of wave pathways in order to create topographic models. In this way, 3D models of P- and S-wave seismic velocity anomalies are obtained which can be interpreted as structural, thermal or compositional variations inside the Earth. Focusing on the Earth’s density structure, the conversion between seismic velocities and mass densities are adopted to construct regional or global seismic density models of the crust and the mantle. Two major limiting aspects restrict possibilities of recovering Earth’s density structure realistically. The first one is practical. Since active seismic experiments are relatively expensive, large parts of the world are not yet covered sufficiently by seismic surveys, most remarkably most of world’s oceans as well as remote parts of Antarctica, Greenland, Africa and South America. The other aspect is of a theoretical nature. The determination of mass density from seismic data could be ambiguous while affected by many uncertainties, meaning that the relationship between seismic velocities and mass densities is not unique. Actually, the density structure inside the Earth is controlled by many factors such a thermal state or mineral composition.  
  
Our JSG is also open for researchers dealing with the classical approaches of gravity field modelling (e.g. the spherical or ellipsoidal harmonics) that are using high performance computing to speed up their processing of enormous amount of input data. This includes large-scale parallel computations on massively parallel architectures as well as heterogeneous parallel computations using graphics processing units (GPUs).
+
Gravity data has been used to interpolate the information about the Earth’s density structure (or density interfaces) where seismic data coverage is uneven or sparse. The National Geospatial-Intelligence Agency in conjunction with its partners from around the world has begun to develop a new global gravitational model, EGM2020, which should be released to public in 2020. EGM2020 should significantly improve the accuracy (as well as the actual resolution) of the global Earth’s gravity field. This will be achieved by incorporating new data sources and procedures. Updated satellite gravity information from the GOCE and GRACE missions will better support the lower harmonics, globally. Multiple new acquisitions (terrestrial, airborne and shipborne) of gravimetric data over specific regions, will provide improved global coverage and resolution over the land as well as for coastal and some oceanic areas. Ongoing accumulation of satellite altimetry data will contribute to refinement and accuracy improvement of the marine gravity field, most notably in polar and near-coastal regions. A significant improvement is also anticipated over large remote regions in Africa, South America, Greenland and Antarctica. EGM2020 will provide opportunities to improve the current knowledge about the Earth’s inner structure and processes particularly in regions with a low seismic data coverage. Gravimetric interpretation of the Earth’s inner density structure is, however, a non-unique problem because infinity many density configurations could be attributed just to the one gravity field solution. Moreover, the gravity inversion is (in a broader mathematical context) an ill-posed problem.  
  
Applications of the aforementioned numerical methods for gravity field modelling involve a detailed discretization of the real Earth’s surface considering its topography. It naturally leads to the oblique derivative problem that needs to be treated. In case of FEM or FVM, unstructured meshes above the topography will be constructed. The meshless methods like MFS or SBM that are based on the point-masses modelling can be applied for processing the gravity gradients observed by the GOCE satellite mission. To reach precise and high-resolution solutions, an elimination of far zones’ contributions is practically inevitable. This can be performed using the fast multipole method or iterative procedures. In both cases such an elimination process improves conditioning of the system matrix and a numerical stability of the problem.  
+
To overcome partially theoretical deficiencies and practical restrictions of both, seismic and gravimetric methods for the recovery of the Earth’s inner density structure, techniques for a combined or constrained inversions of gravity and seismic data are optimally applied, while incorporating additional geophysical, geological and geodynamic constraints. Many such methods already exist or could be developed and further improved within the framework of scientific activities of members of this (multidisciplinary) study group over the next four years. This is achievable, given their expertise in the field of geodesy, geophysics, mathematics and to some extent also geology.  
The aim of the JSG is also to investigate and develop nonlinear filtering methods that allow adaptive smoothing, which effectively reduces the noise while preserves main structures in data. The proposed approach is based on a numerical solution of partial differential equations using a surface finite volume method. It leads to a semi-implicit numerical scheme of the nonlinear diffusion equation on a closed surface where the diffusivity coefficients depend on a combination of the edge detector and a mean curvature of the filtered function. This will avoid undesirable smoothing of local extremes.
+
We expect that our research activities will substantially contribute to the current knowledge of the lithospheric structure, focusing on continental regions of Africa and South America and other continents where seismic data are sparse. Our ongoing research already involves Antarctica and central part of Eurasia. Moreover, a special attention will be given to study the lithospheric structure beneath the Indian Ocean, which is probably the most complex, but the least understood. Despite the lithosphere is the most heterogeneous layer inside the Earth, large lateral structural irregularities are still present even deeper within the mantle below the lithosphere-asthenosphere boundary that are mainly attributed to the mantle convection pattern. The combined gravity and seismic data will be exploited in order to improve existing global or continental-scale mantle density models. A further improvement of the knowledge on the Earth’s inner structure is important, among many other subjects, also for a better understanding of the response of the lithosphere to the mantle convection. This involves numerous study topic, including but not limited to the compensation stage of the crust/lithosphere, the lithospheric strength, mechanisms behind the oceanic subduction, the relation between the mantle convection pattern and the global tectonic configuration (and its spatio-temporal variations), the glacial isostatic adjustment, volcanic processes, or geo-hazard. The members of this study group will address some of these aspects within the following overall objectives.
  
 
===Objectives===
 
===Objectives===

Revision as of 14:09, 9 June 2020

JSG T.25: Advanced computational methods for recovery of high-resolution gravity field models

Chairs: Robert Čunderlík (Slovakia)
Affiliation: Comm. 2 and GGOS

Introduction

The seismic tomography is primarily used to provide images of the Earth’s inner structure based on the analysis of seismic waves due to earthquakes and (controlled) explosions. This technique involves several different methods for processing P-, S- and surface waves on the principle of solving inverse problems for finding locations of reflection and refraction of wave pathways in order to create topographic models. In this way, 3D models of P- and S-wave seismic velocity anomalies are obtained which can be interpreted as structural, thermal or compositional variations inside the Earth. Focusing on the Earth’s density structure, the conversion between seismic velocities and mass densities are adopted to construct regional or global seismic density models of the crust and the mantle. Two major limiting aspects restrict possibilities of recovering Earth’s density structure realistically. The first one is practical. Since active seismic experiments are relatively expensive, large parts of the world are not yet covered sufficiently by seismic surveys, most remarkably most of world’s oceans as well as remote parts of Antarctica, Greenland, Africa and South America. The other aspect is of a theoretical nature. The determination of mass density from seismic data could be ambiguous while affected by many uncertainties, meaning that the relationship between seismic velocities and mass densities is not unique. Actually, the density structure inside the Earth is controlled by many factors such a thermal state or mineral composition.

Gravity data has been used to interpolate the information about the Earth’s density structure (or density interfaces) where seismic data coverage is uneven or sparse. The National Geospatial-Intelligence Agency in conjunction with its partners from around the world has begun to develop a new global gravitational model, EGM2020, which should be released to public in 2020. EGM2020 should significantly improve the accuracy (as well as the actual resolution) of the global Earth’s gravity field. This will be achieved by incorporating new data sources and procedures. Updated satellite gravity information from the GOCE and GRACE missions will better support the lower harmonics, globally. Multiple new acquisitions (terrestrial, airborne and shipborne) of gravimetric data over specific regions, will provide improved global coverage and resolution over the land as well as for coastal and some oceanic areas. Ongoing accumulation of satellite altimetry data will contribute to refinement and accuracy improvement of the marine gravity field, most notably in polar and near-coastal regions. A significant improvement is also anticipated over large remote regions in Africa, South America, Greenland and Antarctica. EGM2020 will provide opportunities to improve the current knowledge about the Earth’s inner structure and processes particularly in regions with a low seismic data coverage. Gravimetric interpretation of the Earth’s inner density structure is, however, a non-unique problem because infinity many density configurations could be attributed just to the one gravity field solution. Moreover, the gravity inversion is (in a broader mathematical context) an ill-posed problem.

To overcome partially theoretical deficiencies and practical restrictions of both, seismic and gravimetric methods for the recovery of the Earth’s inner density structure, techniques for a combined or constrained inversions of gravity and seismic data are optimally applied, while incorporating additional geophysical, geological and geodynamic constraints. Many such methods already exist or could be developed and further improved within the framework of scientific activities of members of this (multidisciplinary) study group over the next four years. This is achievable, given their expertise in the field of geodesy, geophysics, mathematics and to some extent also geology. We expect that our research activities will substantially contribute to the current knowledge of the lithospheric structure, focusing on continental regions of Africa and South America and other continents where seismic data are sparse. Our ongoing research already involves Antarctica and central part of Eurasia. Moreover, a special attention will be given to study the lithospheric structure beneath the Indian Ocean, which is probably the most complex, but the least understood. Despite the lithosphere is the most heterogeneous layer inside the Earth, large lateral structural irregularities are still present even deeper within the mantle below the lithosphere-asthenosphere boundary that are mainly attributed to the mantle convection pattern. The combined gravity and seismic data will be exploited in order to improve existing global or continental-scale mantle density models. A further improvement of the knowledge on the Earth’s inner structure is important, among many other subjects, also for a better understanding of the response of the lithosphere to the mantle convection. This involves numerous study topic, including but not limited to the compensation stage of the crust/lithosphere, the lithospheric strength, mechanisms behind the oceanic subduction, the relation between the mantle convection pattern and the global tectonic configuration (and its spatio-temporal variations), the glacial isostatic adjustment, volcanic processes, or geo-hazard. The members of this study group will address some of these aspects within the following overall objectives.

Objectives

The main objectives of the study group are as follows:

  • to develop algorithms for detailed discretization of the real Earth’s surface including the possibility of adaptive refinement procedures,
  • to create unstructured meshes above the topography for the FVM or FEM approach,
  • to develop the FVM, BEM or FEM numerical models for solving the geodetic BVPs that will treat the oblique derivative problem,
  • to develop numerical models based on MFS or SBM for processing the GOCE observations,
  • to develop parallel implementations of algorithms using the standard MPI procedures,
  • to perform large-scale parallel computations on clusters with distributed memory,
  • to investigate and develop methods for nonlinear diffusion filtering of data on the Earth’s surface where the diffusivity coefficients depend on a combination of the edge detector and a mean curvature of the filtered function,
  • to derive the semi-implicit numerical schemes for the nonlinear diffusion equation on closed surfaces using the surface FVM,
  • and to apply the developed nonlinear filtering methods to real geodetic data.

Program of Activities

  • Active participation at major geodetic workshops and conferences.
  • Organization of group working meetings at main international symposia.
  • Organization of conference sessions.

Members

Róbert Čunderlík (Slovakia), chair
Karol Mikula (Slovakia), vice-chair

Jan Martin Brockmann (Germany)
Walyeldeen Godah (Poland)
Petr Holota (Czech Republic)
Michal Kollár (Slovakia)
Marek Macák (Slovakia)
Zuzana Minarechová (Slovakia)
Otakar Nesvadba (Czech Republic)
Wolf-Dieter Schuh (Germany)