Difference between pages "IC SG2" and "IC SG3"

From Icctwiki
(Difference between pages)
Jump to: navigation, search
(Introduction)
 
 
Line 1: Line 1:
<big>'''JSG 0.2: Gravity field modelling in support of world height system realization'''</big>
+
<big>'''JSG 0.12: Advanced computational methods for recovery of high-resolution gravity field models'''</big>
  
Chair:''P. Novák (Czech Republic)''<br>
+
Chairs: ''Robert Čunderlík (Slovakia)''<br>
Affiliation:''Comm. 2, 1 and GGOS''
+
Affiliation: ''Comm. 2 and GGOS''
  
 
__TOC__
 
__TOC__
 +
 
===Introduction===
 
===Introduction===
  
Description of the Earth’s gravity field still remains a major research topic in geodesy. The main goal is to provide reliable global models covering all spatially-temporal frequencies of its scalar parameterization through the gravity potential. Detailed and accurate gravity field models are required for proper positioning and orientation of geodetic sensors (data geo-referencing). Geometric properties of the gravity field are then studied including those of its equipotential surfaces and their respective surface normals, since they play a fundamental role in definition and realization of geodetic reference systems. Gravity field models will be applied for definition and realization of a vertical reference system (currently under construction) that will support studies of the Earth system.
+
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.
This study group is an entity of the Inter-Commission Committee on Theory. It is affiliated to Commissions 1 (Reference Frames) and 2 (Gravity Field); its close co-operation with GGOS Theme 1 “Unified Global Height System” is anticipated. It aims at bringing together scientists concerned namely with theoretical aspects in the areas of interest specified below.
+
 
 +
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).
 +
 
 +
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.  
 +
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.
  
 
===Objectives===
 
===Objectives===
  
* Considering different types and large amounts of gravity-related data available today, large variety of gravity field models and the ongoing IAG project of realizing a world height system (WHS), this study group shall focuses on theoretical aspects related to the following (non-exhaustive to WHS) list of problems:
+
The main objectives of the study group are as follows:
* To study available gravity field models in terms of their available resolution, accuracy and stability for the pur-pose of WHS realization.
+
* to develop algorithms for detailed discretization of the real Earth’s surface including the possibility of adaptive refinement procedures,
* To define a role of a conventional model of the Earth’s gravity field (EGM) to be used for WHS realization in-cluding its scale parameters.
+
* to create unstructured meshes above the topography for the FVM or FEM approach,
* To study relations between an adopted conventional EGM and parameters of a geocentric reference ellipsoid of revolution approximating a time invariant equipoten-tial surface of the adopted EGM aligned to reduced observables of mean sea level.
+
* to develop the FVM, BEM or FEM numerical models for solving the geodetic BVPs that will treat the oblique derivative problem,
* To study theoretical aspects of various methods proposed for WHS definition and realization including investiga-tions on tidal system effects.
+
* to develop numerical models based on MFS or SBM for processing the GOCE observations,
To investigate combination of heterogeneous gravity field observables by using spatial inversion, spherical radial functions, collocation, wavelets, etc. and by taking into account their sampling geometry, spectral and stochastic properties.
+
* to develop parallel implementations of algorithms using the standard MPI procedures,
* To investigate methods of gravity field modelling based on combination of global gravitational models, ground and airborne gravity, GNSS/levelling height differences, altimetry data, deflections of the vertical, etc.
+
* to perform large-scale parallel computations on clusters with distributed memory,
* To study stable, accurate and efficient methods for con-tinuation of gravity field parameters including space-borne observables of type GRACE and GOCE.
+
* 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 advance theory and methods for solving various initial and boundary value problems (I/BVP) in geodesy.
+
* to derive the semi-implicit numerical schemes for the nonlinear diffusion equation on closed surfaces using the surface FVM,
* To study methods for gravity potential estimation based on its measured directional derivatives (gravity, gravity gradients) by exploiting advantages of simultaneous con-tinuation and inversion of observations.
+
* and to apply the developed nonlinear filtering methods to real geodetic data.
* To investigate requirements for gravity data (stochastic properties, spatially-temporal sampling, spectral content etc.) in terms of their specific geodetic applications.
 
  
 
===Program of Activities===
 
===Program of Activities===
  
Active participation at major geodetic conferences and meetings.
+
* Active participation at major geodetic workshops and conferences.  
Organizing a session at the Hotine-Marussi Symposium 2013.
+
* Organization of group working meetings at main international symposia.
Co-operation with affiliated IAG Commissions and GGOS.
+
* Organization of conference sessions.
Electronic exchange of ideas and thoughts through a SG web page.
 
Monitoring activities of SG members and external individuals related to SG.
 
Compiling bibliography in the area of SG interest.
 
  
 
===Members===
 
===Members===
  
'' '''Pavel Novák (Czech Republic), chair'''<br />Hussein Abd-Elmotaal (Egypt)<br />Robert Čunderlík (Slovakia)<br />Heiner Denker (Germany)<br />Will Featherstone (Australia)<br />René Forsberg (Denmark)<br />Bernhard Heck (Germany)<br />Jianliang Huang (Canada)<br />
+
'' '''Róbert Čunderlík (Slovakia), chair <br /> Karol Mikula (Slovakia), vice-chair''' <br /> Jan Martin Brockmann (Germany) <br /> Walyeldeen Godah (Poland) <br /> Petr Holota (Czech Republic) <br /> Michal Kollár (Slovakia) <br /> Marek Macák (Slovakia) <br />  
Christopher Jekeli (USA)<br />Dan Roman (USA)<br />Fernando Sansò (Italy)<br />Michael G Sideris (Canada)<br />Lars Sjöberg (Sweden)<br />
+
Zuzana Minarechová (Slovakia) <br /> Otakar Nesvadba (Czech Republic) <br /> Wolf-Dieter Schuh (Germany) <br />''
Robert Tenzer (New Zealand)<br />Yan-Ming Wang (USA)<br />''
 

Latest revision as of 12:07, 24 April 2016

JSG 0.12: Advanced computational methods for recovery of high-resolution gravity field models

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

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.

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).

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. 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.

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)